# Physical formulation and numerical algorithm for simulating N immiscible incompressible fluids involving general order parameters

## Abstract

We present a family of physical formulations, and a numerical algorithm, based on a class of general order parameters for simulating the motion of a mixture of N (N⩾2) immiscible incompressible fluids with given densities, dynamic viscosities, and pairwise surface tensions. The N-phase formulations stem from a phase field model we developed in a recent work based on the conservations of mass/momentum, and the second law of thermodynamics. The introduction of general order parameters leads to an extremely strongly-coupled system of (N−1) phase field equations. On the other hand, the general form enables one to compute the N-phase mixing energy density coefficients in an explicit fashion in terms of the pairwise surface tensions. We show that the increased complexity in the form of the phase field equations associated with general order parameters in actuality does not cause essential computational difficulties. Our numerical algorithm reformulates the (N−1) strongly-coupled phase field equations for general order parameters into 2(N−1) Helmholtz-type equations that are completely de-coupled from one another. This leads to a computational complexity comparable to that for the simplified phase field equations associated with certain special choice of the order parameters. We demonstrate the capabilities of the method developed herein using severalmore »

- Authors:

- Publication Date:

- OSTI Identifier:
- 22382186

- Resource Type:
- Journal Article

- Resource Relation:
- Journal Name: Journal of Computational Physics; Journal Volume: 283; Other Information: Copyright (c) 2014 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved.; Country of input: International Atomic Energy Agency (IAEA)

- Country of Publication:
- United States

- Language:
- English

- Subject:
- 71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; ALGORITHMS; COMPARATIVE EVALUATIONS; COMPUTERIZED SIMULATION; DENSITY; ENERGY DENSITY; FIELD EQUATIONS; LIQUIDS; MIXTURES; MULTIPHASE FLOW; ORDER PARAMETERS; SURFACE TENSION; THERMODYNAMICS; VISCOSITY

### Citation Formats

```
Dong, S., E-mail: sdong@purdue.edu.
```*Physical formulation and numerical algorithm for simulating N immiscible incompressible fluids involving general order parameters*. United States: N. p., 2015.
Web. doi:10.1016/J.JCP.2014.11.039.

```
Dong, S., E-mail: sdong@purdue.edu.
```*Physical formulation and numerical algorithm for simulating N immiscible incompressible fluids involving general order parameters*. United States. doi:10.1016/J.JCP.2014.11.039.

```
Dong, S., E-mail: sdong@purdue.edu. Sun .
"Physical formulation and numerical algorithm for simulating N immiscible incompressible fluids involving general order parameters". United States.
doi:10.1016/J.JCP.2014.11.039.
```

```
@article{osti_22382186,
```

title = {Physical formulation and numerical algorithm for simulating N immiscible incompressible fluids involving general order parameters},

author = {Dong, S., E-mail: sdong@purdue.edu},

abstractNote = {We present a family of physical formulations, and a numerical algorithm, based on a class of general order parameters for simulating the motion of a mixture of N (N⩾2) immiscible incompressible fluids with given densities, dynamic viscosities, and pairwise surface tensions. The N-phase formulations stem from a phase field model we developed in a recent work based on the conservations of mass/momentum, and the second law of thermodynamics. The introduction of general order parameters leads to an extremely strongly-coupled system of (N−1) phase field equations. On the other hand, the general form enables one to compute the N-phase mixing energy density coefficients in an explicit fashion in terms of the pairwise surface tensions. We show that the increased complexity in the form of the phase field equations associated with general order parameters in actuality does not cause essential computational difficulties. Our numerical algorithm reformulates the (N−1) strongly-coupled phase field equations for general order parameters into 2(N−1) Helmholtz-type equations that are completely de-coupled from one another. This leads to a computational complexity comparable to that for the simplified phase field equations associated with certain special choice of the order parameters. We demonstrate the capabilities of the method developed herein using several test problems involving multiple fluid phases and large contrasts in densities and viscosities among the multitude of fluids. In particular, by comparing simulation results with the Langmuir–de Gennes theory of floating liquid lenses we show that the method using general order parameters produces physically accurate results for multiple fluid phases.},

doi = {10.1016/J.JCP.2014.11.039},

journal = {Journal of Computational Physics},

number = ,

volume = 283,

place = {United States},

year = {Sun Feb 15 00:00:00 EST 2015},

month = {Sun Feb 15 00:00:00 EST 2015}

}