Search NASA⌕ Search

SEARCH · Search NASA

Results for “VOF”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

Spray Combustion Modeling with VOF and Finite-Rate Chemistry

A spray atomization and combustion model is developed based on the volume-of-fluid (VOF) transport equation with finite-rate chemistry model. The gas-liquid interface mass, momentum and energy conservation laws are modeled by continuum surface force mechanisms. A new solution method is developed such that the present VOF model can be applied for all-speed range flows. The objectives of the present study are: (1) to develop and verify the fractional volume-of-fluid (VOF) cell partitioning approach into a predictor-corrector algorithm to deal with multiphase (gas-liquid) free surface flow problems; (2) to implement the developed unified algorithm in a general purpose computational fluid dynamics (CFD) code, Finite Difference Navier-Stokes (FDNS), with droplet dynamics and finite-rate chemistry models; and (3) to demonstrate the effectiveness of the present approach by simulating benchmark problems of jet breakup/spray atomization and combustion. Modeling multiphase fluid flows poses a significant challenge because a required boundary must be applied to a transient, irregular surface that is discontinuous, and the flow regimes considered can range from incompressible to highspeed compressible flows. The flow-process modeling is further complicated by surface tension, interfacial heat and mass transfer, spray formation and turbulence, and their interactions. The major contribution of the present method is to combine the novel feature of the Volume of Fluid (VOF) method and the Eulerian/Lagrangian method into a unified algorithm for efficient noniterative, time-accurate calculations of multiphase free surface flows valid at all speeds. The proposed method reformulated the VOF equation to strongly couple two distinct phases (liquid and gas), and tracks droplets on a Lagrangian frame when spray model is required, using a unified predictor-corrector technique to account for the non-linear linkages through the convective contributions of VOF. The discontinuities within the sharp interface will be modeled as a volume force to avoid stiffness. Formations of droplets, tracking of droplet dynamics and modeling of the droplet breakup/evaporation, are handled through the same unified predictor-corrector procedure. Thus the new algorithm is non-iterative and is flexible for general geometries with arbitrarily complex topology in free surfaces. The FDNS finite-difference Navier-Stokes code is employed as the baseline of the current development. Benchmark test cases of shear coaxial LOX/H2 liquid jet with atomization/combustion and impinging jet test cases are investigated in the present work. Preliminary data comparisons show good qualitative agreement between data and the present analysis. It is indicative from these results that the present method has great potential to become a general engineering design analysis and diagnostics tool for problems involving spray combustion.

Chen, Yen-Sen↗
NASA NTRS2015

Scale Model Acoustic Test Validation of IOP-SS Water Prediction using Loci-STREAM-VoF

The Scale Model Acoustic Test (SMAT) is a 5% scale test of the Space Launch System (SLS), which is currently being designed at Marshall Space Flight Center (MSFC). SMAT consists of a 5% scale representation of the ignition overpressure sound-suppression system (IOP-SS) that is being tested to quantify the water flow and induced air entrainment in and around the mobile launcher exhaust hole. This data will be compared with computational fluid dynamics (CFD) simulations using the newly developed Loci-STREAM Volume of Fluid (VoF) methods. Compressible and incompressible VoF methods have been formulated, and are currently being used to simulate the water flow of SMAT IOP-SS. The test data will be used to qualitatively and quantitatively assess and validate the VoF methods.

Nielsen, Tanner↗
DOE OSTI2024

Capturing thin structures in VOF simulations with two-plane reconstruction

A novel interface reconstruction strategy for volume of fluid (VOF) methods is introduced that represents the liquid-gas interface as two planes that co-exist within a single computational cell. In comparison to the piecewise linear interface calculation (PLIC), this new algorithm greatly improves the accuracy of the reconstruction, in particular when dealing with thin structures such as films. The placement of the two planes requires the solution of a non-linear optimization problem in six dimensions, which has the potential to be overly expensive. Further, an efficient solution to this optimization problem is presented here that exploits two key ideas: an algorithm for extracting multiple plane orientations from transported surface data, and an efficient and mass-conserving distance-finding algorithm that accounts for two planes with arbitrary orientation. Additionally, a simple and robust strategy is presented to accurately represent the surface tension forces produced at the interface of subgrid-thickness films. The performance of this new VOF reconstruction is demonstrated on several test cases that illustrate the capability to handle arbitrarily thin films.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1995

Numerical modeling of spray combustion with an advanced VOF method

This paper summarizes the technical development and validation of a multiphase computational fluid dynamics (CFD) numerical method using the volume-of-fluid (VOF) model and a Lagrangian tracking model which can be employed to analyze general multiphase flow problems with free surface mechanism. The gas-liquid interface mass, momentum and energy conservation relationships are modeled by continuum surface mechanisms. A new solution method is developed such that the present VOF model can be applied for all-speed flow regimes. The objectives of the present study are to develop and verify the fractional volume-of-fluid cell partitioning approach into a predictor-corrector algorithm and to demonstrate the effectiveness of the present approach by simulating benchmark problems including laminar impinging jets, shear coaxial jet atomization and shear coaxial spray combustion flows.

Chen, Yen-Sen↗
NASA NTRS1994

Atomization simulations using an Eulerian-VOF-Lagrangian method

This paper summarizes the technical development and validation of a multiphase computational fluid dynamics (CFD) numerical method using the volume-of-fluid (VOF) model and a Lagrangian tracking model which can be employed to analyze general multiphase flow problems with free surface mechanism. The gas-liquid interface mass, momentum and energy conservations are modeled by continuum surface mechanisms. A new solution method is developed such that the present VOF model can be applied for all-speed flow regimes. The objectives of the present study are to develop and verify the fractional volume-of-fluid cell partitioning approach into a predictor-corrector algorithm and to demonstrate the effectiveness of the present innovative approach by simulating benchmark problems including the coaxial jet atomization.

Chen, Yen-Sen↗
NASA NTRSDate not supplied

Launching Liquid Slugs Through Tube Networks: Investigating Multiphase Liquid Slug Acceleration Behaviors in Spacecraft Propulsion Systems using 3D Multiphase VOF CFD Simulations and High-Speed Flow-Visualization Validation Experiments

Propellant management in spacecraft propulsion systems is complex. In many current and future crewed and robotic spacecraft propulsion systems, a number of unintended or unavoidable scenarios can occasionally lead to liquid “slugs” of propellant accumulating in unwanted locations and thereafter being rapidly accelerated through complex and arbitrary tubing networks otherwise filled with vapor or gas. To develop better understanding of the multiphase fluid dynamics associated with these “slug launch” scenarios, NASA Marshal Space Flight Center’s (MSFC) ER42 branch – the propulsion fluid dynamics research branch – conducted 3D multiphase transient Volume of Fluid (VOF) CFD simulations using an in-house tool called LOCI-Stream VOF on semi-canonical tube networks with straight sections, bends, and bend complexes. In addition, high-speed flow visualization validation experiments using propellants were conducted on similar tube network geometries at the NASA White Sands Test Facility (WSTF). These simulations and flow-visualization experiments found these scenarios lead to a variety of unique multiphase fluid dynamic effects involving the leading and trailing edges of the slug in both straight sections and curved bends as well as the progressive erosion, acceleration, and eventual “punch-through” of the driving gas through the liquid slug. Understanding these multiphase fluid dynamic effects has led to unique and often counter-intuitive observations which have been of critical relevance to important engineering concerns associated with these “slug launch” scenarios.

Multiphase Fluid Dynamics↗
DOE OSTI2020

Materials Data on VOF by Materials Project

VOF is zeta iron carbide-derived structured and crystallizes in the monoclinic P2_1/c space group. The structure is three-dimensional. there are two inequivalent V3+ sites. In the first V3+ site, V3+ is bonded to three O2- and three F1- atoms to form a mixture of edge and corner-sharing VO3F3 octahedra. The corner-sharing octahedra tilt angles range from 46–57°. There is one shorter (1.94 Å) and two longer (2.00 Å) V–O bond length. There are a spread of V–F bond distances ranging from 2.07–2.17 Å. In the second V3+ site, V3+ is bonded to three O2- and three F1- atoms to form a mixture of edge and corner-sharing VO3F3 octahedra. The corner-sharing octahedra tilt angles range from 43–55°. There is one shorter (1.93 Å) and two longer (2.02 Å) V–O bond length. There are a spread of V–F bond distances ranging from 2.04–2.13 Å. There are two inequivalent O2- sites. In the first O2- site, O2- is bonded in a distorted trigonal planar geometry to three V3+ atoms. In the second O2- site, O2- is bonded in a 3-coordinate geometry to three V3+ atoms. There are two inequivalent F1- sites. In the first F1- site, F1- is bonded in a distorted trigonal planar geometry to three V3+ atoms. In the second F1- site, F1- is bonded in a distorted trigonal planar geometry to three V3+ atoms.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Materials Data on VOF by Materials Project

VOF is beta Vanadium nitride-derived structured and crystallizes in the orthorhombic Ama2 space group. The structure is three-dimensional. there are three inequivalent V3+ sites. In the first V3+ site, V3+ is bonded to three O2- and three F1- atoms to form a mixture of edge and corner-sharing VO3F3 octahedra. The corner-sharing octahedra tilt angles range from 48–52°. There is one shorter (1.93 Å) and two longer (1.99 Å) V–O bond length. There are two shorter (2.10 Å) and one longer (2.14 Å) V–F bond lengths. In the second V3+ site, V3+ is bonded to three O2- and three F1- atoms to form a mixture of edge and corner-sharing VO3F3 octahedra. The corner-sharing octahedra tilt angles range from 47–54°. There is one shorter (1.92 Å) and two longer (1.98 Å) V–O bond length. There are one shorter (2.13 Å) and two longer (2.14 Å) V–F bond lengths. In the third V3+ site, V3+ is bonded to three O2- and three F1- atoms to form a mixture of edge and corner-sharing VO3F3 octahedra. The corner-sharing octahedra tilt angles range from 47–54°. There are a spread of V–O bond distances ranging from 1.92–2.00 Å. There are a spread of V–F bond distances ranging from 2.10–2.15 Å. There are three inequivalent O2- sites. In the first O2- site, O2- is bonded in a distorted trigonal planar geometry to three V3+ atoms. In the second O2- site, O2- is bonded in a distorted trigonal planar geometry to three V3+ atoms. In the third O2- site, O2- is bonded in a distorted trigonal planar geometry to three V3+ atoms. There are three inequivalent F1- sites. In the first F1- site, F1- is bonded in a 3-coordinate geometry to three V3+ atoms. In the second F1- site, F1- is bonded in a 3-coordinate geometry to three V3+ atoms. In the third F1- site, F1- is bonded in a distorted trigonal planar geometry to three V3+ atoms.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Materials Data on VOF by Materials Project

VOF is beta Vanadium nitride-derived structured and crystallizes in the monoclinic P2/m space group. The structure is three-dimensional. there are five inequivalent V3+ sites. In the first V3+ site, V3+ is bonded to three O2- and three F1- atoms to form VO3F3 octahedra that share corners with eight VO2F4 octahedra and edges with two equivalent VO3F3 octahedra. The corner-sharing octahedra tilt angles range from 47–52°. There is one shorter (1.89 Å) and two longer (1.98 Å) V–O bond length. There are two shorter (2.13 Å) and one longer (2.17 Å) V–F bond lengths. In the second V3+ site, V3+ is bonded to two O2- and four F1- atoms to form VO2F4 octahedra that share corners with eight VO4F2 octahedra and edges with two equivalent VO2F4 octahedra. The corner-sharing octahedra tilt angles range from 49–50°. Both V–O bond lengths are 1.93 Å. All V–F bond lengths are 2.07 Å. In the third V3+ site, V3+ is bonded to two equivalent O2- and four equivalent F1- atoms to form VO2F4 octahedra that share corners with eight equivalent VO3F3 octahedra and edges with two equivalent VO2F4 octahedra. The corner-sharing octahedra tilt angles range from 49–52°. Both V–O bond lengths are 1.94 Å. All V–F bond lengths are 2.07 Å. In the fourth V3+ site, V3+ is bonded to four O2- and two F1- atoms to form VO4F2 octahedra that share corners with eight VO3F3 octahedra and edges with two equivalent VO4F2 octahedra. The corner-sharing octahedra tilt angles range from 47–50°. There are two shorter (2.01 Å) and two longer (2.02 Å) V–O bond lengths. There are one shorter (2.10 Å) and one longer (2.14 Å) V–F bond lengths. In the fifth V3+ site, V3+ is bonded to four equivalent O2- and two equivalent F1- atoms to form VO4F2 octahedra that share corners with eight equivalent VO2F4 octahedra and edges with two equivalent VO4F2 octahedra. The corner-sharing octahedra tilt angles range from 49–50°. All V–O bond lengths are 2.01 Å. Both V–F bond lengths are 2.13 Å. There are four inequivalent O2- sites. In the first O2- site, O2- is bonded in a distorted trigonal planar geometry to three V3+ atoms. In the second O2- site, O2- is bonded in a distorted trigonal planar geometry to three V3+ atoms. In the third O2- site, O2- is bonded in a distorted trigonal planar geometry to three V3+ atoms. In the fourth O2- site, O2- is bonded in a distorted trigonal planar geometry to three V3+ atoms. There are four inequivalent F1- sites. In the first F1- site, F1- is bonded in a distorted trigonal planar geometry to three V3+ atoms. In the second F1- site, F1- is bonded in a distorted trigonal planar geometry to three V3+ atoms. In the third F1- site, F1- is bonded in a 3-coordinate geometry to three V3+ atoms. In the fourth F1- site, F1- is bonded in a distorted trigonal planar geometry to three V3+ atoms.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Materials Data on VOF by Materials Project

Computed materials data using density functional theory calculations. These calculations determine the electronic structure of bulk materials by solving approximations to the Schrodinger equation. For more information, see https://materialsproject.org/docs/calculations

36 MATERIALS SCIENCE↗
NASA NTRS2003

Vector-Ordering Filter Procedure for Data Reduction

The vector-ordering filter (VOF) technique involves a procedure for sampling a large population of data vectors to select a subset of data vectors that fully characterize the state space of the large population. The VOF technique enables a large reduction of the volume of data that must be handled in the automated monitoring system and method discussed in the two immediately preceding articles. In so doing, the VOF technique enables the development of data-driven mathematical models of a monitored asset from sets of data that would otherwise exceed the memory capacities of conventional engineering computers. Data-driven mathematical models have been shown to offer high fidelity for purposes of control and monitoring of assets. In practice, a collection of asset-operating observations is acquired with the intention that the collection contain observations characteristic of the full dynamic range of operation of the asset. Often, such a collection contains an extremely large number of observations, many of which are redundant. The VOF technique fills the need for a means to extract, from the original collection of observational data, a reduced data matrix that excludes redundant data while maintaining the full statistical character and dynamic range of the original data. The reduced data matrix can then be used as the input data for development of a mathematical model of the monitored asset, or as training data for a neural-network substitute for an explicit mathematical model of the asset. Alternatively, the reduced data matrix can, itself, be used directly as a mathematical model of the monitored asset, as is commonly done in multivariate state-estimation techniques. The original data are collected from the asset over a range of operating states and are put in matrix form. Each column vector in the original data matrix represents the signal values acquired at a particular operational state of the asset. Thus, the number of columns of the original data matrix equals the number of observed states and the number of rows in this matrix equals the number of signals acquired at each observation. In the VOF technique, one extracts the reduced data matrix from the original data matrix through the selection of a representative subset of the column (state) vectors.

Bickford, Randall L.↗
NASA NTRSDate not supplied

Modeling Autogenous Pressurization and Draining of a Cryogenic Storage Tank in Normal Gravity

A two-phase CFD model for autogenous pressurization and draining of a cryogenic storage tank is presented using both the Sharp Interface and Volume-Of-Fluid (VOF) approaches for capturing the front and the associated interfacial heat, mass and momentum transfer between the liquid and vapor regions. Both models are validated against data provided by the Cryogenic Propellant Storage and Transfer (CPST) Engineering Development Unit (EDU) experiment. The results of the autogenous pressurization are presented first, focusing on the phase change and turbulence effects on the tank pressure and temperature predictions. Both the Sharp Interface (SI-CFD) and VOF (VOF-CFD) multiphase models predict tank pressure during pressurization within 3% of the measured values. The sensitivity of key physical and numerical parameters of the problem are tested using the Sharp Interface model. Effects of the accommodation coefficient (AC), and the computational grid structure are studied. The second part of this paper is devoted to validating the SI-VOF model, with an enhanced capability of moving the liquid-vapor interface, against the draining segment of the EDU experiment. The VOF-CFD model was also used to simulate tank draining and its results were compared with the results of the Sharp Interface model with the moving interface. Both models predict tank pressure during draining within 3.5% of the measured values. Pressure decrease rate is underpredicted by both models during the first 100 seconds of draining but matches the experimental rate for the rest of the simulation.

Computational Fluid Dynamics↗
NASA NTRSDate not supplied

Modeling Autogenous Pressurization and Draining of a Cryogenic Storage Tank in Normal Gravity

A two-phase CFD model for autogenous pressurization and draining of a cryogenic storage tank is presented using both the Sharp Interface and Volume-Of-Fluid (VOF) approaches for capturing the front and the associated interfacial heat, mass and momentum transfer between the liquid and vapor regions. Both models are validated against data provided by the Cryogenic Propellant Storage and Transfer (CPST) Engineering Development Unit (EDU) experiment1. The results of the autogenous pressurization are presented first, focusing on the phase change and turbulence effects on the tank pressure and temperature predictions. Both the Sharp Interface (SI-CFD) and VOF (VOF-CFD) multiphase models predict tank pressure during pressurization within 3% of the measured values. The sensitivity of key physical and numerical parameters of the problem are tested using the Sharp Interface model. Effects of the accommodation coefficient (AC), and the computational grid structure are studied. The second part of this paper is devoted to validating the SI-VOF model, with an enhanced capability of moving the liquid-vapor interface, against the draining segment of the EDU experiment. The VOF-CFD model was also used to simulate tank draining and its results were compared with the results of the Sharp Interface model with the moving interface. Both models predict tank pressure during draining within 3.5% of the measured values. Pressure decrease rate is underpredicted by both models during the first 100 seconds of draining but matches the experimental rate for the rest of the simulation.

Computational Fluid Dynamics↗
DOE OSTI2025

A Finite Element Method for Compressible and Turbulent Multiphase Flow Instabilities with Heat Transfer

We present a new finite element framework for modeling compressible, turbulent multiphase flows with heat transfer. For two-fluid systems with a free surface, the Volume of Fluid (VOF) method is implemented without the need for interface reconstruction, while turbulence is resolved using a dynamic Vreman large eddy simulation (LES) model. Unlike most two-phase VOF studies, which neglect heat transfer, the present approach incorporates energy transport equations within the VOF formulation to account for heat exchange, an effect particularly important in turbulent flows. Conjugate heat transfer is often challenging in finite volume methods, which require explicit specification of heat fluxes at the solid–fluid interface, limiting accuracy and predictive capability. By contrast, the finite element formulation does not require heat flux inputs, allowing more accurate and robust simulation of heat transfer between solids and fluids. The method is demonstrated through three representative cases. First, a two-fluid instability with a single-mode perturbation is simulated and validated against analytical growth rates. Second, conjugate heat transfer is examined in a high-temperature flow over a cold metal cylinder, with validation performed both quantitatively—via pressure coefficient comparisons with experimental data—and qualitatively using vector field topology. Finally, compressible spray injection and breakup are modeled, demonstrating the ability of the framework to capture interfacial dynamics and atomization under turbulent, high-speed conditions. In the compressible spray injection and breakup case, the results indicate that the finite element formulation achieved higher predictive accuracy and robustness than the finite-volume method. With the same mesh resolution, the FEM reduced the root mean square error (RMSE) and mean absolute percentage error (MAPE) from 6.96 mm and 26.0% (for the FVM) to 4.85 mm and 12.7%, respectively, demonstrating improved accuracy and robustness in capturing interfacial dynamics and heat transfer. The study also introduced vector field topology to visualize and interpret coherent flow structures and instabilities, offering insights beyond conventional scalar-field analyses.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Non-Diffusive Volume Advection with A High Order Interface Reconstruction Method

We show that non-diffusive volume advection in two-dimensions is achieved with several benchmark problems using a newly developed high-order volume of fluids (VOF) interface reconstruction method. (1) A new VOF interface reconstruction method using circular/corner facets (linear facets are a degenerate case of arcs). We create a circular interface facet in each mixed zone by matching neighbor volume with a hybrid Newton’s-bisection method and the local solution is final. In the general case, the new VOF interface reconstruction has 3rd order accuracy and can be easily made seamless. The new method addresses intrinsic issues with Young’s method such as gaps between interface facets in the case of a curved interface, and inability to define curvature nor identify corners. (2) A non-diffusive volume advection scheme. In an ALE advection step, a well-defined interface can be carried over through a Lagrange step and used to compute volume distribution into a relaxed mesh. Then, an interface reconstruction step is performed to redefine the interface in the relaxed mesh. We must point out that the interface carried over is also a solution of interface reconstruction because all the volume fractions in the relaxed mesh are naturally matched. We provide an interface tracking method compatible with our reconstruction scheme, where it is granted to use the prior info as an initial guess to capture sub-mesh resolution features. As a result, we are able to treat multiple facets inside a single mixed cell and obtain highly accurate, non-diffusive solution for advection problems with rather coarse meshes. We show our solutions for two-dimensional incompressible flows with two materials with a) the X + O diagonal translation; b) the Zalesak rotational test; and c) the single vortex spiral test.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2015

CFD Modeling of the Multipurpose Hydrogen Test Bed (MHTB) Self-Pressurization and Spray Bar Mixing Experiments in Normal Gravity: Effect of the Accommodation Coefficient on the Tank Pressure

A CFD model for simulating the self-pressurization of a large scale liquid hydrogen storage tank is utilized in this paper to model the MHTB self-pressurization experiment. The kinetics-based Schrage equation is used to account for the evaporative and condensi ng interfacial mass flows in this model. The effect of the accommodation coefficient for calculating the interfacial mass transfer rate on the tank pressure during tank selfpressurization is studied. The values of the accommodation coefficient which were considered in this study vary from 1.0e-3 to 1.0e-1 for the explicit VOF model and from 1.0e-4 to 1.0e-3 for the implicit VOF model. The ullage pressure evolutions are compared against experimental data. A CFD model for controlling pressure in cryogenic storage tanks by spraying cold liquid into the ullage is also presented. The Euler-Lagrange approach is utilized for tracking the spray droplets and for modeling the interaction between the droplets and the continuous phase (ullage). The spray model is coupled with the VOF model by performing particle tracking in the ullage, removing particles from the ullage when they reach the interface, and then adding their contributions to the liquid. Droplet-ullage heat and mass transfer are modeled. The flow, temperature, and interfacial mass flux, as well as droplets trajectories, size distribution and temperatures predicted by the model are presented. The ul lage pressure and vapor temperature evolutions are compared with experimental data obtained from the MHTB spray bar mixing experiment. The effect of the accommodation coefficient for calculating the interfacial and droplet mass transfer rates on the tank pressure during mixing of the vapor using spray is studied. The values used for the accommodation coefficient at the interface vary from 1.0e-5 to 1.0e-2. The droplet accommodation coefficient values vary from 2.0e-6 to 1.0e-4.

cryogenic fluid storage↗
DOE OSTI2022

The role of breakup and coalescence in fine-scale bubble-induced turbulence. II. Kinematics

This second part of our research explores the kinematic aspect of fine-scale bubble-induced turbulence (BIT) to (i) present the effect of bubble breakup and coalescence and (ii) compare it against the universal kinematic fine-scale turbulence characteristics reported in the literature. To this end, we simulate a dilute bubbly system of 0.5% void fraction using two distinct numerical simulations. In the volume-of-fluid (VoF) simulation, bubbles undergo breakup and coalescence. In the immersed boundary method (IBM) simulation, however, they act as rigid spheres. We also perform a simulation of classical homogeneous isotropic turbulence (HIT). The first important outcome of this study is that BIT is radically different from HIT in terms of its kinematic fine-scale characteristics. In the vorticity-dominating regions, BIT exhibits a weak vortex stretching. This weak vortex stretching is due to (a) the intermediate strain-rate eigenvalues skewed weakly to positive and (b) the extensive strain-rate eigenvector aligning perpendicular to the vorticity vector. The BIT has, on average, not only a weak enstrophy production but also a weak strain production in strain-dominating regions. The weak strain production is due to (a) the presence of vortex stretching in highly strained fluid elements and (b) the absolute magnitude of compressive strain-rate eigenvalue being as close to the extensive strain-rate eigenvalue. Thus, none of the kinematic fine-scale HIT characteristics is noted for BIT. The second important conclusion is that bubble breakup and coalescence play little to no influence on the kinematics of fine-scale BIT as VoF and IBM simulations produce similar results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Numerical Study of Liquid Piston Compression Using Large-Eddy Simulation and Volume-of-Fluid Approach

Efforts to increase the efficiency of residential and commercial air conditioners and heat pumps have demonstrated that the compressor accounts for most of the system’s electrical energy usage. Therefore, the efficiency of this component should be improved to reduce its energy usage. The US Department of Energy’s Oak Ridge National Laboratory developed a near-isothermal liquid piston compressor (LPC) that uses propylene glycol (PG) to compress CO2. This report presents numerical studies of the LPC in which the compression chamber fills with injected PG from the bottom inlet. Numerical simulations were performed using the large-eddy simulation (LES) with the wall-adapting local eddy-viscosity (WALE) subgrid-scale model coupled with the multiphase volume-of-fluid (VOF) model to simulate the transient interface between gas and liquid and to capture the heat and mass transfer within the compression chamber. In this effort, effects of boundary conditions applied to the LES-VOF calculations (e.g., no wall, an adiabatic wall, and a wall with a heat flux subscribed) to the overall pressure and temperature of CO2 gas as well as the transient evolution of flow and heat transfer evolution within the compression chamber are investigated and discussed. It was found that the LES calculation with no wall have shown no dynamical flow patterns and the volume-averaged temperature of CO2 increased from 305 to 392.7 K, while LES calculations with a constant wall temperature or a wall heat flux had similar increases of CO2 temperatures. Results of LES simulation using a wall heat flux showed different stages in the compression process and revealed dynamical formation and interaction of CO2 gas layers and circulation flow patterns within the chamber that contribute to the overall heat transfer between the solid wall, gas, and liquid surface in the compressor.

Nguyen, Thien D.↗