[petsc-users] How to solve N-S equation with SIMPLE-like prediction-correction method?

[Zhang Qingwen]

Hi, all,

In the field of CFD, the prediction-correction method (e.g., the SIMPLE algorithm) has been widely used to solved the Navier-Stokes equations with finite volume method.

However, I have no idea about how it can be implemented with PETSc.

Let me describe a simplified case of creeping flow (e.g., mantle flow), in which the inertia term in the N-S equ. can be neglected and reduced to the so-called Stokes equ.
This equation only contains a diffusion term and source term of pressure gradient, and there are two coupled equations to solved:

  *   The equation of momentum conservation
  *   The equation of mass conservation

The SIMPLE does not solve these coupled equations simultaneously. Instead, this can be done in a iterative fashion with SEPERATE steps until the solutions of velocity and pressure are accurate enough:

1. set initial u*, v*, P*
2. solve the momentum equation to get the velocity
3. solve the continuity equation to get the pressure correction P'
4. correct the velocity and the pressure with P' to get the new u, v, P
5. update the solutions: u*=u, v*=v, P*=P
6. go to step 2 if the solutions are not accurate enough

(A document is attached to clearify the governing equation and the SIMPLE algorithm, please check the attachment)

There is a cavity flow example using TS that demonstrate how to solve N-S equations simultaneously,
which solve four coupled equations with nonlinear iterations inside SNES and TS, and does not involves separate prediction-correction steps like the SIMPLE algorithm above.

So, any hints about how can the SIMPLE algorithm be implemented in PETSc?
Particularly, how can I solve these equations in separate steps, so I can solve the governing equations in a prediction-correction fashion?

Best,
Qw
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20220301/06e9f97b/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: SIMPLE-like_algorithm.pdf
Type: application/pdf
Size: 227159 bytes
Desc: SIMPLE-like_algorithm.pdf
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20220301/06e9f97b/attachment-0001.pdf>