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 Nphase 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 stronglycoupled system of (N−1) phase field equations. On the other hand, the general form enables one to compute the Nphase 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) stronglycoupled phase field equations for general order parameters into 2(N−1) Helmholtztype equations that are completely decoupled 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., Email: 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., Email: 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., Email: sdong@purdue.edu. 2015.
"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., Email: 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 Nphase 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 stronglycoupled system of (N−1) phase field equations. On the other hand, the general form enables one to compute the Nphase 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) stronglycoupled phase field equations for general order parameters into 2(N−1) Helmholtztype equations that are completely decoupled 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 = 2015,
month = 2
}

This paper presents a general numerical formulation for simulating physical dispersion that avoids the convergence and gridorientation problems associated with numerical dispersion. The formulation is applicable to a wide range of problems, including enhanced gasdrives with phase discontinuities. A new ninepoint scheme is also presented, along with general rules for any ninepoint scheme to obey Darcy's law. Converged solutions are obtained for displacements at adverse mobility ratios of 5,000 with no gridorientation effects.

Superconducting states of reduced symmetry: General order parameters and physical implications
It is shown that the order parameter in reducedsymmetry superconductors can be represented in terms of a complete set of basisfunction multiplets, which is analogous to a complete set of crystal harmonics. These complete sets are found for several symmetries that are germane to the study of heavyfermion superconductors. The implications of a general order parameter are discussed for several physical properties: line nodes in the energy gap, the Knight shift, surface pairbreaking, and timereversal symmetry. 
Numerical simulation of directcontact evaporation of a drop rising in a hot, less volatile immiscible liquid of higher density  Possibilities and limits of the SOLAVOF/CSF algorithm
A contribution toward the full numerical simulation of directcontact evaporation of a drop rising in a hot, immiscible and less volatile liquid of higher density is presented. Based on a fixedgrid Eulerian description, the classical SOLAVOF method is largely extended to incorporate, for example, three incompressible fluids and liquidvapor phase change. The thorough validation and assessment process covers several benchmark simulations, some which are presented, documenting the multipurpose value of the new code. The directcontact evaporation simulations reveal severe numerical problems that are closely related to the fixedgrid Euler formulation. As a consequence, the comparison to experiments have to bemore »