Publication Cover
Numerical Heat Transfer, Part B: Fundamentals
An International Journal of Computation and Methodology
Volume 72, 2017 - Issue 4
344
Views
6
CrossRef citations to date
0
Altmetric
Original Articles

Modeling multiphase flow: Spray breakup using volume of fluids in a dynamics LES FEM method

, &
Pages 285-299 | Received 10 Aug 2017, Accepted 26 Oct 2017, Published online: 17 Nov 2017
 

ABSTRACT

Currently, all commercial software for engine modeling investigates the dispersed droplet phase of the injection process. Understanding the effect of geometry of the injector nozzle, initial jet conditions, fluid properties in the liquid film, breakup, resulting droplet sizes, and distribution are of primary importance to improve fuel efficiency and lower gas emissions. We have developed an innovative computational method and models to make this atomization process more predictive: a multiscale, multiphase fluid simulation, using a volume-of-fluid method implemented in a large eddy simulation algorithm found in the new KIVA-hpFE, a finite element method flow solver for all flow regimes.

Nomenclature

Cp=

specific heat capacity at constant P (J/kg · K)

c=

sound speed (m/s)

E=

total internal energy (J/kg)

e=

any element

fs=

surface tension force

=

unit normal of the interface surface

P=

pressure (Pa)

T=

temperature (K)

tij=

grid-scale (resolved scale) shear stress

=

intermediate velocity

Ve=

volume of computational element

Greek symbols=
ρ=

density (kg/m3)

β=

artificial compressibility

=

mass fraction

σ=

surface tension coefficient

=

Dirac delta function

κ=

surface curvature

μ=

fluid viscosity (

=

turbulent eddy viscosity

ν=

kinematic viscosity

τij=

subgrid-scale stress tensor

ϕ=

volume of fraction

ϕE=

elemental value for volume of fraction

Nomenclature

Cp=

specific heat capacity at constant P (J/kg · K)

c=

sound speed (m/s)

E=

total internal energy (J/kg)

e=

any element

fs=

surface tension force

=

unit normal of the interface surface

P=

pressure (Pa)

T=

temperature (K)

tij=

grid-scale (resolved scale) shear stress

=

intermediate velocity

Ve=

volume of computational element

Greek symbols=
ρ=

density (kg/m3)

β=

artificial compressibility

=

mass fraction

σ=

surface tension coefficient

=

Dirac delta function

κ=

surface curvature

μ=

fluid viscosity (

=

turbulent eddy viscosity

ν=

kinematic viscosity

τij=

subgrid-scale stress tensor

ϕ=

volume of fraction

ϕE=

elemental value for volume of fraction

Additional information

Funding

The DOE’s Office of Energy Efficiency and Renewable Energy (EERE) Advanced Combustion Program (Gurpreet Singh and Leo Breton) is supporting this effort. Los Alamos National Laboratory, an affirmative action/equal opportunity employer, is operated by the Los Alamos National Security, LLC for the National Nuclear Security Administration of the U.S. Department of Energy (DOE) under contract DE-AC52-06NA25396. Los Alamos National Laboratory strongly supports academic freedom and a researcher’s right to publish; as an institution, however, the Laboratory does not endorse the viewpoint of a publication or guarantee its technical correctness.

Reprints and Corporate Permissions

Please note: Selecting permissions does not provide access to the full text of the article, please see our help page How do I view content?

To request a reprint or corporate permissions for this article, please click on the relevant link below:

Academic Permissions

Please note: Selecting permissions does not provide access to the full text of the article, please see our help page How do I view content?

Obtain permissions instantly via Rightslink by clicking on the button below:

If you are unable to obtain permissions via Rightslink, please complete and submit this Permissions form. For more information, please visit our Permissions help page.