Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiphase flow”

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 1,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
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.↗
DOE OSTI2016

Mineral assemblage transformation of a metakaolin-based waste form after geopolymer encapsulation

Current plans for nuclear waste vitrification at the Hanford Tank Waste Treatment and Immobilization Plant (WTP) lack the capacity to treat all of the low activity waste (LAW) that is not encapsulated in the vitrified product. Fluidized Bed Steam Reforming (FBSR) is one of the supplemental technologies under consideration to fill this gap. The FBSR process results in a granular product mainly composed of feldspathoid mineral phases that encapsulate the LAW and other contaminants of concern (COCs). In order to better understand the characteristics of the FBSR product, characterization testing has been performed on the granular product as well as the granular product encapsulated in a monolithic geopolymer binder. The non-radioactive simulated tank waste samples created for use in this study are the result of a 2008 Department of Energy sponsored Engineering Scale Technology Demonstration (ESTD) in 2008. These samples were created from waste simulant that was chemically shimmed to resemble actual tank waste, and rhenium has been used as a substitute for technetium. Another set of samples was created by the Savannah River Site Bench-Scale Reformer (BSR) using a chemical shim of Savannah River Site Tank 50 waste in order to simulate a blend of 68 Hanford tank wastes. This paper presents results from coal and moisture removal tests along with XRD, SEM, and BET analyses showing that the major mineral components are predominantly sodium aluminosilicate minerals and that the mineral product is highly porous. Results also show that the materials pass the short-term leach tests: the Toxicity Characteristic Leaching Procedure (TCLP) and Product Consistency Test (PCT).

Feldspathoid Minerals↗
DOE OSTI2023

Constraining the Reaction Rate of Criegee Intermediates with Carboxylic Acids during the Multiphase Ozonolysis of Aerosolized Alkenes

Criegee intermediates (CIs) are central to the ozonolysis of unsaturated organic compounds, controlling the formation of small fragmentation reaction products and higher molecular weight oligomeric compounds through competing unimolecular and bimolecular pathways. In particular, the reaction of thermalized CI with carboxylic acids (RCOOH) to produce oligomeric α-acyloxyalkyl hydroperoxides (AAHPs) is of interest to the chemical transformation of organic compounds in the atmosphere. In the gas phase, the CI + RCOOH reaction proceeds at supercollisional rates, consistent with a barrierless reaction mechanism. However, the rate of this reaction in condensed phases remains uncertain. Here, we report the results of aerosol flow tube studies of the ozonolysis of submicron organic aerosol containing mixtures of an internal n-alkene (Z-9-tricosene, tri) and a saturated acid (2-hexyldecanoic acid, HDA). The decay of the reactants and formation of reaction products are monitored using atmospheric pressure chemical ionization mass spectrometry (APCI–MS). A 6-fold decrease in the reactive uptake of ozone is observed with increasing acid concentration, consistent with a previous study using alcohol additives. The decay of HDA with ozone exposure shows that HDA scavenges CIs at rates comparable to competing unimolecular pathways involving CIs. The decay kinetics of HDA are analyzed using a simple kinetic model in order to constrain the ratio of the CI + RCOOH and unimolecular CI rate constants. Empirical fitting to this model, combined with theoretical estimates of the unimolecular loss rate, yields a CI + RCOOH rate constant of 1.85 ± 0.27 × 10 –19 cm 3 molec –1 s –1 , which is six orders of magnitude slower than the expected diffusion limit in a liquid organic matrix. Stochastic kinetic simulations are conducted with different values of the CI + RCOOH rate constant to validate this result and show that the diffusion-limited rate is indeed much too rapid to reproduce the experimental results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

Above zone pressure interpretation for leaky well characterization and its identification from leaky caprock/fault

Pressure transient analysis has long been used for reservoir characterization. Above-zone (AZ) pressure has been recently investigated for inferring leakage pathway characteristics in leakage events from subsurface injection operations. The recorded pressure in the AZ should be purely related to leakage and therefore it can be safely inverted to deduce leakage characteristics. It is crucial to evaluate fluid leakage through abandoned wells to plan for further measurements of leakage prevention. However, the recorded AZ leakage signal may not be related to leaky well(s). Therefore, identification and spatial investigation of well leakage is required for leakage evaluation. Here in this paper, we propose a pressure interpretation method for early detection of leaky pathways, applying two observation points in the AZ. We distinguish leaky well, fault and caprock based on their corresponding flow regime identification. We show that the pressure difference of the two observation wells can be applied as a proxy for unknown leakage rate, which is crucial for leakage identification as well as characterization. Results show that the estimated location of the leaky well, leakage coefficient, and the leakage rate are in good agreement with the actual values. The estimated leakage coefficient of the leaky well can be used to evaluate well leakage in multiphase systems such as carbon dioxide leakage in deep saline aquifers.

02 PETROLEUM↗
DOE OSTI2015

Liquid jet breakup regimes at supercritical pressures

Previously, a theory has been presented that explains how discrete vapor–liquid interfaces become diminished at certain high-pressure conditions in a manner that leads to well known qualitative trends observed from imaging in a variety of experiments. Rather than surface tension forces, transport processes can dominate over relevant ranges of conditions. In this paper, this framework is now generalized to treat a wide range of fuel-oxidizer combinations in a manner consistent with theories of capillary flows and extended corresponding states theory. Different flow conditions and species-specific molecular properties are shown to produce distinct variations of interfacial structures and local free molecular paths. These variations are shown to occur over the operating ranges in a variety of propulsion and power systems. Despite these variations, the generalized analysis reveals that the envelope of flow conditions at which the transition from classical sprays to diffusion-dominated mixing occurs exhibits a characteristic shape for all liquid–gas combinations. As a result, for alkane-oxidizer mixtures, it explains that these conditions shift to higher pressure flow conditions with increasing carbon number and demonstrates that, instead of widely assumed classical spray atomization, diffusion-dominated mixing may occur under relevant high-pressure conditions in many modern devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRSDate not supplied

Overview of Extractive and Optical PM Measurements from the NASA SE-5 High Pressure Flame Diagnostics Laboratory

Soot emissions measurements of multiphase non-premixed flames are sparse in the literature due to the complexities of sampling in such environments. We present an overview from on-going work at the NASA SE-5 High Pressure Flame Diagnostics laboratory focusing on soot and gaseous emissions at elevated pressures. SE-5 is a gaseous and liquid fueled high pressure flame tube facility capable of operating at pressures up to 880 psia, cooling capacity of 4,000,000 BTU/hr, liquid fuel flow rates up to 2 GPH, and an equivalence ratio range of 0.2 (fuel lean) to 4 (fuel rich). SE-5 provides optical access via four UV-grade fused silica optical windows of the primary reaction zone enabling non-intrusive optical diagnostics of the flame.

Soot↗
DOE OSTI2015

GMC COLLISIONS AS TRIGGERS OF STAR FORMATION. I. PARAMETER SPACE EXPLORATION WITH 2D SIMULATIONS

We utilize magnetohydrodynamic (MHD) simulations to develop a numerical model for giant molecular cloud (GMC)–GMC collisions between nearly magnetically critical clouds. The goal is to determine if, and under what circumstances, cloud collisions can cause pre-existing magnetically subcritical clumps to become supercritical and undergo gravitational collapse. We first develop and implement new photodissociation region based heating and cooling functions that span the atomic to molecular transition, creating a multiphase ISM and allowing modeling of non-equilibrium temperature structures. Then in 2D and with ideal MHD, we explore a wide parameter space of magnetic field strength, magnetic field geometry, collision velocity, and impact parameter and compare isolated versus colliding clouds. We find factors of ∼2–3 increase in mean clump density from typical collisions, with strong dependence on collision velocity and magnetic field strength, but ultimately limited by flux-freezing in 2D geometries. For geometries enabling flow along magnetic field lines, greater degrees of collapse are seen. We discuss observational diagnostics of cloud collisions, focussing on {sup 13}CO(J = 2–1), {sup 13}CO(J = 3–2), and {sup 12}CO(J = 8–7) integrated intensity maps and spectra, which we synthesize from our simulation outputs. We find that the ratio of J = 8–7 to lower-J emission is a powerful diagnostic probe of GMC collisions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Mechanically Robust, Sodium-Ion Conducting Membranes for Nonaqueous Redox Flow Batteries

Sodium-based batteries are promising for grid-storage applications because of significantly lower cost compared to lithium-based systems. The advancement of solid-state and redox-flow sodium-ion batteries requires sodium-ion exchange membranes with high conductivity, electrochemical stability, and mechanical robustness. This study demonstrates that membranes based on poly(ethylene oxide) (PEO) can meet these requirements. Membranes plasticized with tetraethylene glycol dimethyl ether (TEGDME) achieve high ionic conductivity. Plasticized PEO membranes containing sodium triflate salt (NaTFS) show about 2 orders of magnitude higher conductivity compared to nonplasticized PEO membranes. Results from vibrational spectroscopy and differential scanning calorimetry describe the coordination chemistry in these multiphase materials and explain the mechanisms behind the increased conductivity. The mechanical properties of the membranes improve by addition of 5 wt % sodium carboxymethyl cellulose (CMC) without compromising the conductivity or electrochemical stability against sodium metal. In conclusion, the optimized membrane is an excellent candidate for low-cost energy storage systems that operate over a wide voltage window near ambient temperature.

25 ENERGY STORAGE↗
DOE OSTI2019

Assessment on DNB performance for high pressure subcooled boiling with modified CASL boiling model

Understanding and predicting the Departure from Nucleation Boiling phenomena in the nuclear power industry has been a grand challenge due to its large impact on the reactor safety and operational economics. Over the past several years, CASL THM team has extensively investigated and demonstrated a state of art Multiphase CFD framework to provide a high-fidelity DNB simulation capability based on the two-fluid Eulerian model in conjunction with a wall heat flux partitioning approach. The overarching goal of this task is to develop and demonstrate a robust CFD-based DNB model which is capable of predicting an entire boiling history up to the DNB condition in which the wall boiling mode is transitioning from the nucleate boiling mode to the film boiling mode. Three different flow channel configurations (pipe flow, 5x5 fuel bundle with and without mixing vane tests) were studied and thorough closure sensitivity studies were performed to evaluate the wall boiling and momentum closure effect on the DNB performance. Predicated DNB values were compared to the measurements from high pressure subcooled flow boiling tests resulting in a consistent under-predictive trend for all three validation campaigns in FY16-18. With extensive closure sensitivity study, a modified CASL boiling closure model (i.e. Gen1-FY19) is introduced with the following key changes: 1) improved wall boiling closure for high pressure application, 2) newly calibrated lift and wall lubrication force models with modified interfacial length scale model, and 3) updated liquid phase condensation model calibrated for the subcooled boiling application.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS2005

The nakhlite meteorites: Augite-rich igneous rocks from Mars

The seven nakhlite meteorites are augite-rich igneous rocks that formed in flows or shallow intrusions of basaltic magma on Mars. They consist of euhedral to subhedral crystals of augite and olivine (to 1 cm long) in fine-grained mesostases. The augite crystals have homogeneous cores of Mg' = 63% and rims that are normally zoned to iron enrichment. The core-rim zoning is cut by iron-enriched zones along fractures and is replaced locally by ferroan low-Ca pyroxene. The core compositions of the olivines vary inversely with the steepness of their rim zoning - sharp rim zoning goes with the most magnesian cores (Mg' = 42%), homogeneous olivines are the most ferroan. The olivine and augite crystals contain multiphase inclusions representing trapped magma. Among the olivine and augite crystals is mesostasis, composed principally of plagioclase and/or glass, with euhedra of titanomagnetite and many minor minerals. Olivine and mesostasis glass are partially replaced by veinlets and patches of iddingsite, a mixture of smectite clays, iron oxy-hydroxides and carbonate minerals. In the mesostasis are rare patches of a salt alteration assemblage: halite, siderite, and anhydrite/ gypsum. The nakhlites are little shocked, but have been affected chemically and biologically by their residence on Earth. Differences among the chemical compositions of the nakhlites can be ascribed mostly to different proportions of augite, olivine, and mesostasis. Compared to common basalts, they are rich in Ca, strongly depleted in Al, and enriched in magmaphile (incompatible) elements, including the LREE. Nakhlites contain little pre-terrestrial organic matter. Oxygen isotope ratios are not terrestrial, and are different in anhydrous silicates and in iddingsite. The alteration assemblages all have heavy oxygen and heavy carbon, while D/H values are extreme and scattered. Igneous sulfur had a solar-system isotopic ratio, but in most minerals was altered to higher and lower values. High precision analyses show mass-independent fractionations of S isotopes. Nitrogen and noble gases are complex and represent three components: two mantle sources (Chas-E and Chas-S), and fractionated Martian atmosphere. The nakhlites are igneous cumulate rocks, formed from basaltic magma at approx.1.3 Ga, containing excess crystals over what would form from pure magma. After accumulation of their augite and olivine crystals, they were affected (to various degrees) by crystallization of the magma, element diffusion among minerals and magma, chemical reactions among minerals and magma, magma movement among the crystals, and post-igneous chemical equilibration. The extent of these modifications varies, from least to greatest, in the order: MIL03346, NWA817, Y000593, Nakhla = Governador Valadares, Lafayette, and NWA998. Chemical, isotopic, and chronologic data confirm that the nakhlites formed on Mars, most likely in thick lava flows or shallow intrusions. Their crystallization ages, referenced to crater count chronologies for Mars, suggest that the nakhlites formed on the large volcanic constructs of Tharsis, Elysium, or Syrtis Major. The nakhlites were suffused with liquid water, probably at approx.620 ma. This water dissolved olivine and mesostasis glass, and deposited iddingsite and salt minerals in their places. The nakhlites were ejected from Mars at approx.10.75Ma by an asteroid impact and fell to Earth within the last 10,000 years. Although the nakhlites are enriched in incompatible elements, their source mantle was strongly depleted. This depletion event was ancient, as the nakhlites source mantle was fractionated while short-lived radionuclides (e.g., t(sub 1/2 = 9 my) were still active. This differentiation event may have been core formation coupled with a magma ocean, as is inferred for the moon.

Treiman, Allan H.↗
DOE OSTI2022

Equation of state and strength of diamond in high-pressure ramp loading

Diamond is used extensively as a component in high-energy-density experiments, but existing equation of state (EOS) models do not capture its observed response to dynamic loading. In particular, in contrast with first-principles theoretical models of equilibrium multiphase EOS, no solid-solid phase changes have been detected, and no general-purpose EOS models match the measured ambient isotherm. We have performed density functional theory (DFT) calculations of the diamond phase to ~10 TPa, well beyond its predicted range of thermodynamic stability, and used these results as the basis of a Mie-Grüneisen EOS. We also performed DFT calculations of the elastic moduli, and calibrated an algebraic elasticity model for use in simulations. We then estimated the flow stress of diamond by comparison with the stress-density relation measured experimentally in ramp-loading experiments. Furthermore, the resulting constitutive model allows us to place a constraint on the Taylor-Quinney factor (the fraction of plastic work converted to heat) from the observation that diamond does not melt on ramp compression.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1989

High-temperature ordered intermetallic alloys III; Proceedings of the Third Symposium, Boston, MA, Nov. 29-Dec. 1, 1988

The present conference on high-temperature ordered intermetallic alloys discusses alloy theory and phase stability, defects and microstructures, mechanical behavior, alloy design and material processing, multiphase and composite materials, nickel aluminides, titanium aluminides, and other ordered intermetallics. Attention is given to compositional effects on processing and properties of nickel aluminides, dynamic compaction and hot-isostatic-pressing of nickel aluminides, improvement of elevated-temperature material properties in Ni3Al, the effect of microalloying B on the high-temperature mechanical properties of Ti3Al, the effects of structure on creep of Ti-53.4 mol pct Al intermetallics, the deformation microstructure in Ni3Si polycrystals strained over the range of temperature of flow stress anomaly, and the microstructure and mechanical properties of dual phase alloys consisting of the intermetallic phases.

Liu, C. T.↗
DOE OSTI2014

Plasmas in Multiphase Media: Bubble Enhanced Discharges in Liquids and Plasma/Liquid Phase Boundaries

In this research project, the interaction of atmospheric pressure plasmas with multi-phase media was computationally investigated. Multi-phase media includes liquids, particles, complex materials and porous surfaces. Although this investigation addressed fundamental plasma transport and chemical processes, the outcomes directly and beneficially affected applications including biotechnology, medicine and environmental remediation (e.g., water purification). During this project, we made advances in our understanding of the interaction of atmospheric pressure plasmas in the form of dielectric barrier discharges and plasma jets with organic materials and liquids. We also made advances in our ability to use computer modeling to represent these complex processes. We determined the method that atmospheric pressure plasmas flow along solid and liquid surfaces, and through endoscopic like tubes, deliver optical and high energy ion activation energy to organic and liquid surfaces, and produce reactivity in thin liquid layers, as might cover a wound. We determined the mechanisms whereby plasmas can deliver activation energy to the inside of liquids by sustaining plasmas in bubbles. These findings are important to the advancement of new technology areas such as plasma medicine

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Development of a coupled thermo-hydro-mechanical model in discontinuous media for carbon sequestration

Geomechanical alteration of porous media is generally ignored for most shallow subsurface applications, whereas CO 2 injection, migration, and trapping in deep saline aquifers will be controlled by coupled multifluid flow, energy transfer, and geomechanical processes. The accurate assessment of the risks associated with potential leakage of injected CO 2 and the design of effective injection systems requires that we represent these coupled processes within numerical simulators. The objectives of this study were to develop a coupled thermal-hydro-mechanical model into a single software, and to examine the coupling of thermal, hydrological, and geomechanical processes for simulation of CO 2 injection into the subsurface for carbon sequestration. A numerical model is developed to couple nonisothermal multiphase hydrological and geomechanical processes for prediction of multiple interconnected processes for carbon sequestration in deep saline aquifers. The geomechanics model was based on Rigid Body-Spring Model (RBSM), one of the discrete methods to model discontinuous rock system. Poisson’s effect that was often ignored by RBSM was considered in the model. The simulation of large-scale and long-term coupled processes in carbon capture and storage projects requires large memory and computational performance. Global Array Toolkit was used to build the model to permit the high performance simulations of the coupled processes. The model was used to simulate a case study with several scenarios to demonstrate the impacts of considering coupled processes and Poisson’s effect for the prediction of CO 2 sequestration.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Predicting CO 2 Plume Migration in Heterogeneous Formations Using Conditional Deep Convolutional Generative Adversarial Network

Numerical simulation of flow and transport in heterogeneous formations has long been studied, especially for uncertainty quantification and risk assessment. The high computational cost associated with running large-scale numerical simulations in a Monte Carlo sense has motivated the development of surrogate models, which aim to capture the important input-output relations of physics-based models but require only a fraction of the cost of full model runs. In this work, we formulate a conditional deep convolutional generative adversarial network (cDC-GAN) surrogate model to learn the dynamic functional mappings in multiphase models. The cDC-GAN belongs to a class of semisupervised learning methods that can be used to learn the data generation processes. Like the original GAN, a main strength of the cDC-GAN is that it includes a self-training scheme for improving the quality of generative modeling in a game theoretic framework, without requiring extensive statistical knowledge and assumptions on input data distributions. In particular, our cDC-GAN model is designed to learn cross-domain mappings between high-dimensional input (e.g., permeability) and output (e.g., phase saturations) pairs, with the ability to incorporate conditioning information (e.g., prediction time). As a use case, we demonstrate the performance of cDC-GAN for predicting the migration of carbon dioxide (CO 2 ) plume in heterogeneous carbon storage reservoirs, which has both numerical and practical significance because of the safe storage requirements now mandated in many countries. Results show that cDC-GAN achieves high accuracy in predicting the spatial and temporal evolution patterns of the injected CO 2 plume, as compared to the original results obtained using a compositional reservoir simulator. The performance of cDC-GAN models, trained using the same number of training samples, stays relatively robust when the level of spatial heterogeneity is increased. Our cDC-GAN is pattern based and is not limited by the underlying physics. Thus, it provides a general framework for developing surrogate models, and for conducting uncertainty analyses for a wide range of physics-based models used in both groundwater and subsurface energy exploration applications.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1986

Shock waves and shock tubes; Proceedings of the Fifteenth International Symposium, Berkeley, CA, July 28-August 2, 1985

A detailed survey is presented of shock tube experiments, theoretical developments, and applications being carried out worldwide. The discussions explore shock tube physics and the related chemical, physical and biological science and technology. Extensive attention is devoted to shock wave phenomena in dusty gases and other multiphase and heterogeneous systems, including chemically reactive mixtures. Consideration is given to techniques for measuring, visualizing and theoretically modeling flowfield, shock wave and rarefaction wave characteristics. Numerical modeling is explored in terms of the application of computational fluid dynamics techniques to describing flowfields in shock tubes. Shock interactions and propagation, in both solids, fluids, gases and mixed media are investigated, along with the behavior of shocks in condensed matter. Finally, chemical reactions that are initiated as the result of passage of a shock wave are discussed, together with methods of controlling the evolution of laminar separated flows at concave corners on advanced reentry vehicles.

Bershader, D.↗
DOE OSTI2022

Two- and three-particle nonflow contributions to the chiral magnetic effect measurement by spectator and participant planes in relativistic heavy ion collisions

Correlation measurements with respect to the spectator and participant planes in relativistic heavy ion collisions were proposed to extract the chiral magnetic effect (CME) from background dominated azimuthal correlators. Here, this paper investigates the effects of two- and three-particle nonflow correlations on the extracted CME signal fraction, f CME . It is found, guided by a multiphase transport (ampt) model and the heavy ion jet interaction generator (hijing) together with experimental data, that the nonflow effects amount to approximately (4 ± 5)% and (–5 ± 3)% without and with pseudorapidity gaps, respectively, in 20–50% centrality Au + Au collisions at √ s NN = 200 GeV.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

BHR equations re-derived with immiscible particle effects

Compressible and variable density turbulent flows with dispersed phase effects are found in many applications ranging from combustion to cloud formation. These types of flows are among the most challenging to simulate. While the exact equations governing a system of particles and fluid are known, computational resources limit the scale and detail that can be simulated in this type of problem. Therefore, a common method is to simulate averaged versions of the flow equations, which still capture salient physics and is relatively less computationally expensive. Besnard et al. [1] developed such a model for variable density miscible turbulence, where ensemble-averaging was applied to the flow equations to yield a set of filtered equations. Besnard et al. further derived transport equations for the Reynolds stresses, the turbulent mass flux, and the density-specific volume covariance, to help close the filtered momentum and continuity equations. We re-derive the exact BHR closure equations which include integral terms owing to immiscible effects. Physical interpretations of the additional terms are proposed along with simple models. The goal of this work is to extend the BHR model to allow for the simulation of turbulent flows where an immiscible dispersed phase is non-trivially coupled with the carrier phase.

20 FOSSIL-FUELED POWER PLANTS↗