A mass and momentum conserving unsplit semiLagrangian framework for simulating multiphase flows
Abstract
In this work, we present a computational methodology for convection and advection that handles discontinuities with second order accuracy and maintains conservation to machine precision. This method can transport a variety of discontinuous quantities and is used in the context of an incompressible gas–liquid flow to transport the phase interface, momentum, and scalars. The proposed method provides a modification to the threedimensional, unsplit, secondorder semiLagrangian flux method of Owkes & Desjardins (JCP, 2014). The modification adds a refined grid that provides consistent fluxes of mass and momentum defined on a staggered grid and discrete conservation of mass and momentum, even for flows with large density ratios. Additionally, the refined grid doubles the resolution of the interface without significantly increasing the computational cost over previous nonconservative schemes. This is possible due to a novel partitioning of the semiLagrangian fluxes into a small number of simplices. The proposed scheme is tested using canonical verification tests, rising bubbles, and an atomizing liquid jet.
 Authors:
 Mechanical and Industrial Engineering, Montana State University, Bozeman, MT 59717 (United States)
 Sibley School of Mechanical and Aerospace Engineering, Cornell University, Ithaca, NY 14853 (United States)
 Publication Date:
 OSTI Identifier:
 22622258
 Resource Type:
 Journal Article
 Resource Relation:
 Journal Name: Journal of Computational Physics; Journal Volume: 332; Other Information: Copyright (c) 2016 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:
 75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND SUPERFLUIDITY; ACCURACY; ADVECTION; COMPUTERIZED SIMULATION; CONVECTION; DENSITY; INTERFACES; LAGRANGIAN FUNCTION; LIQUID FLOW; LIQUIDS; MODIFICATIONS; MULTIPHASE FLOW; RESOLUTION; THREEDIMENSIONAL CALCULATIONS; VERIFICATION
Citation Formats
Owkes, Mark, Email: mark.owkes@montana.edu, and Desjardins, Olivier. A mass and momentum conserving unsplit semiLagrangian framework for simulating multiphase flows. United States: N. p., 2017.
Web. doi:10.1016/J.JCP.2016.11.046.
Owkes, Mark, Email: mark.owkes@montana.edu, & Desjardins, Olivier. A mass and momentum conserving unsplit semiLagrangian framework for simulating multiphase flows. United States. doi:10.1016/J.JCP.2016.11.046.
Owkes, Mark, Email: mark.owkes@montana.edu, and Desjardins, Olivier. Wed .
"A mass and momentum conserving unsplit semiLagrangian framework for simulating multiphase flows". United States.
doi:10.1016/J.JCP.2016.11.046.
@article{osti_22622258,
title = {A mass and momentum conserving unsplit semiLagrangian framework for simulating multiphase flows},
author = {Owkes, Mark, Email: mark.owkes@montana.edu and Desjardins, Olivier},
abstractNote = {In this work, we present a computational methodology for convection and advection that handles discontinuities with second order accuracy and maintains conservation to machine precision. This method can transport a variety of discontinuous quantities and is used in the context of an incompressible gas–liquid flow to transport the phase interface, momentum, and scalars. The proposed method provides a modification to the threedimensional, unsplit, secondorder semiLagrangian flux method of Owkes & Desjardins (JCP, 2014). The modification adds a refined grid that provides consistent fluxes of mass and momentum defined on a staggered grid and discrete conservation of mass and momentum, even for flows with large density ratios. Additionally, the refined grid doubles the resolution of the interface without significantly increasing the computational cost over previous nonconservative schemes. This is possible due to a novel partitioning of the semiLagrangian fluxes into a small number of simplices. The proposed scheme is tested using canonical verification tests, rising bubbles, and an atomizing liquid jet.},
doi = {10.1016/J.JCP.2016.11.046},
journal = {Journal of Computational Physics},
number = ,
volume = 332,
place = {United States},
year = {Wed Mar 01 00:00:00 EST 2017},
month = {Wed Mar 01 00:00:00 EST 2017}
}

Lagrangian particle model for multiphase flows
A Lagrangian particle model for multiphase multicomponent fluid flow, based on smoothed particle hydrodynamics (SPH), was developed and used to simulate the flow of an emulsion consisting of bubbles of a nonwetting liquid surrounded by a wetting liquid. In SPH simulations, fluids are represented by sets of particles that are used as discretization points to solve the NavierStokes fluid dynamics equations. In the multiphase multicomponent SPH model, a modified van der Waals equation of state is used to close the system of flow equations. The combination of the momentum conservation equation with the van der Waals equation of state resultsmore » 
A spectrally refined interface approach for simulating multiphase flows
This paper presents a novel approach to phaseinterface transport based on pseudospectral subgrid refinement of a level set function. In each flow solver grid cell, a set of quadrature points is introduced on which the value of the level set function is known. This methodology allows to define a polynomial reconstruction of the level set function in each cell. The transport is performed using a semiLagrangian technique, removing all constraints on the time step size. Such an approach provides subcell resolution of the phaseinterface and leads to excellent accuracy in the transport, while a reasonable cost is obtained by precomputingmore » 
A nominally secondorder cellcentered Lagrangian scheme for simulating elastic–plastic flows on twodimensional unstructured grids
In this paper, we describe a cellcentered Lagrangian scheme devoted to the numerical simulation of solid dynamics on twodimensional unstructured grids in planar geometry. This numerical method, utilizes the classical elasticperfectly plastic material model initially proposed by Wilkins [M.L. Wilkins, Calculation of elastic–plastic flow, Meth. Comput. Phys. (1964)]. In this model, the Cauchy stress tensor is decomposed into the sum of its deviatoric part and the thermodynamic pressure which is defined by means of an equation of state. Regarding the deviatoric stress, its time evolution is governed by a classical constitutive law for isotropic material. The plasticity model employs themore »