next up previous contents index
Next: 10.2.4 Results Up: DIMEFEM: High-level Portable Previous: 10.2.2 Operations and Elements

10.2.3 Navier-Stokes Solver


We have implemented a sophisticated incompressible flow solver using DIME and DIMEFEM. The algorithm is described more completely in [Bristeau:87a]. The evolution equation for an incompressible Newtonian fluid of viscosity n is

We use a three-stage operator-split scheme whereby for each time step of length dt, the equation is integrated

The parameter is .

Each of these three implicit steps involves the solution of either a Stokes problem: 

or the nonlinear problem:\

where is a parameter inversely proportional to the time step. We solve the Navier-Stokes equation,  and consequently also these subsidiary problems, with given velocity at the boundary (Dirichlet boundary conditions).

Guy Robinson
Wed Mar 1 10:19:35 EST 1995