by Yi Zhang
In last post I talked wrote about a time integration for NS equation. As a special case, time update for potential flow is actually easier, due to fact that
1. the governing equation, i.e. the Laplace equation is linear.
2. the solution dependes only on the boundary condition, in which only the free surface boundary condition is the tricky part.
3. Even though the field equation is in Eularian representation, the free surface can be described using Lagrangian representation.
For free surface potential flow, the time updating of the problem is twofold: update the (free surface) geometry and update the boundary conditioin at the free surface. Without mentioning the boundary condition at other part of the domain boundary, we are solving the problem characterized by
Free surface kinetic boundary condition(FSKBC):
Free surface dynamic boundary condition(FSDBC):
The first boundary condition says that the particle velocity for material particle on the free surface, i.e. the velocity of free surface, equals to the flow velocity defined by velocity potential. The second boundary condition is the Bernoulli equation, with constant on RHS assigned to 0, which could be looked as consumed by . Remembering that
FSDBC can be converted to
The plan of mixed-Eularian-Lagrangian(MEL) method, is to solve at time , then use the BC equations to update the boundary condition. One should remember, the LHS of BC equations use material particle as variable, i.e. in Lagrangian representation. Denote particle on the free surface by , MEL is characterized by two steps:
1. Solve field equation with BCs with variable at , i.e. ;
2. Update the BCs at the free surface using
As mentioned at the end of last post, usually the computing cost is characterized by the first step, while the second step concerns the stability and convergence.