Hello,

I am working on a CFD code in which the Navier Stokes equations are solved for a one phase fluidum. I want to extend this to a two phase model by solving the Navier Stokes equations for each phase seperately (which are coupled by a drag force term). So I get two velocity fields. What should I do with the pressure? Two pressure fields or just one? In case of two pressure fields can I use the SIMPLE method (by Patankar) for each phase? Hopefully somebody has experience with this. Thanks in advance.

Marcel

Except for surface tension effects (Laplace pressure inside a droplet with surface tension and small radius) the pressure fields should be the same.

Hi there,

You may look at the following review papers first. (1) Annu. Rev. Fluid Mech. 31: 567¨C603; (2) J. Comput. Phys. 169: 463¨C502; (3) J. Comput. Phys. 169: 708¨C759.

What you are thinking is probably called the ¡°sharp interface method¡±, and you may find the following papers and the references therein useful. J. Comput. Phys. 148: 366-396. J. Comput. Phys. 171: 822¨C842. J. Comput. Phys. 160: 151-178.

We too have introduced some concepts in dealing with the singularities of discontinuous pressure and normal derivatives of velocity. (J. Comput. Phys. 187: 255¨C273).

Recently, we have a paper in press particularly on the projection scheme with the pressure jump. If you are interested, I can send you the preprint.

Hope this helps.

Tony

Tony, can you send me the preprint? Thanks