Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Multiphase flows”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Spatial characterization of heterogeneous nanopore surfaces from XCT scans of Niobrara shale

The multiphase fluid flow through natural and geo-architected nanopore structures is fundamentally controlled by the geometry of such nanopore structures, pore wettability and interfacial tension. Furthermore, multiphase fluid flow through the nanopore structures is compromised when the pore surfaces are heterogeneous, having surface regions of different composition, polarity, and wetting characteristics. In this regard, compositional analysis of the pore channel surface has been a big challenge in the past, since 3D analysis is required. Currently, due to advances in nano X-ray CT, with a voxel resolution of 20 nm, it is possible to identify nanopore network structures and quantify their surface composition in 3D. In this paper, the development of procedures for detailed analysis to quantify the composition of the pore channel surface regions, including the corresponding area, is described using nanoCT images. Image analysis procedures, including preprocessing of the 3D tomographic images, correction for the partial volume effect (PVE), and the use of a new algorithm for determination of pore channel surface composition and area are reported on. Videos of the nanopore structure have been provided as supplementary information to further describe the heterogeneity of the nanopore surfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Fluid characterization using acoustics

A measurement system and a method for determining steam quality (i.e. vapor mass fraction) measurements of multiphase fluid flowing through pipes are described. An acoustic sensor device consists of an acoustic transmitter and an acoustic receiver that are designed to be attached to a pipe. The acoustic transmitter and the acoustic receiver are exposed to an interior space of the pipe through openings in a wall of the pipe. Acoustic waves generated by the transmitter and captured by the receiver traverse the multiphase fluid flowing into the pipe. Swept-frequency acoustic interferometry (SFAI) technique is used to measure ultrasonic acoustic properties of a fluid. Machine-learning techniques based on principal component analysis, support vector machine regression and support vector machine classification are used for determining steam quality.

47 OTHER INSTRUMENTATION↗
DOE OSTI2016

A new efficient formulation of the HLLEM Riemann solver for general conservative and non-conservative hyperbolic systems

In this paper a new, simple and universal formulation of the HLLEM Riemann solver (RS) is proposed that works for general conservative and non-conservative systems of hyperbolic equations. For non-conservative PDE, a path-conservative formulation of the HLLEM RS is presented for the first time in this paper. The HLLEM Riemann solver is built on top of a novel and very robust path-conservative HLL method. It thus naturally inherits the positivity properties and the entropy enforcement of the underlying HLL scheme. However, with just the slight additional cost of evaluating eigenvectors and eigenvalues of intermediate characteristic fields, we can represent linearly degenerate intermediate waves with a minimum of smearing. For conservative systems, our paper provides the easiest and most seamless path for taking a pre-existing HLL RS and quickly and effortlessly converting it to a RS that provides improved results, comparable with those of an HLLC, HLLD, Osher or Roe-type RS. This is done with minimal additional computational complexity, making our variant of the HLLEM RS also a very fast RS that can accurately represent linearly degenerate discontinuities. Our present HLLEM RS also transparently extends these advantages to non-conservative systems. For shallow water-type systems, the resulting method is proven to be well-balanced. Several test problems are presented for shallow water-type equations and two-phase flow models, as well as for gas dynamics with real equation of state, magnetohydrodynamics (MHD & RMHD), and nonlinear elasticity. Since our new formulation accommodates multiple intermediate waves and has a broader applicability than the original HLLEM method, it could alternatively be called the HLLI Riemann solver, where the “I” stands for the intermediate characteristic fields that can be accounted for. -- Highlights: •New simple and general path-conservative formulation of the HLLEM Riemann solver. •Application to general conservative and non-conservative hyperbolic systems. •Inclusion of sub-structure and resolution of intermediate characteristic fields. •Well-balanced for single- and two-layer shallow water equations and multi-phase flows. •Euler equations with real equation of state, MHD equations, nonlinear elasticity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Efficient implicit LES method for the simulation of turbulent cavitating flows

We present a numerical method for efficient large-eddy simulation of compressible liquid flows with cavitation based on an implicit subgrid-scale model. Phase change and subgrid-scale interface structures are modeled by a homogeneous mixture model that assumes local thermodynamic equilibrium. Unlike previous approaches, emphasis is placed on operating on a small stencil (at most four cells). The truncation error of the discretization is designed to function as a physically consistent subgrid-scale model for turbulence. We formulate a sensor functional that detects shock waves or pseudo-phase boundaries within the homogeneous mixture model for localizing numerical dissipation. In smooth regions of the flow field, a formally non-dissipative central discretization scheme is used in combination with a regularization term to model the effect of unresolved subgrid scales. The new method is validated by computing standard single- and two-phase test-cases. Comparison of results for a turbulent cavitating mixing layer obtained with the new method demonstrates its suitability for the target applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

High-speed impact of micron-sized diesel drop trains—Splashing dynamics, secondary droplet formation, and effects of pre-existing film thickness

This work investigates surface impingement of monodispersed diesel drop trains using computational simulations. The three-dimensional simulations include high-speed impact of micrometer-sized fuel drops onto initially dry and wetted stainless steel substrates. The drop size and impact velocity are representative of fuel injection conditions in internal combustion engines (ICEs). The drop trains serve as a simplified representation of fuel spray. To accurately capture the film that naturally develops on initially dry substrates, a dynamic diesel-stainless steel contact angle model was developed using experiments of single drop impact under ICE representative temperature. Simulations are presented for a highly splashing case and the effects of a pre-existing film on splashing dynamics were investigated, including the temporal evolution of splashed mass and film thickness. It was concluded that for thin films, the effect of pre-existing film thickness is less significant than in single drop impingement. Secondary droplet characterization was performed on simulation results, leading to instantaneous and time-averaged distributions of secondary droplet size, velocity magnitude, and trajectory angle. It was found that for each drop impingement, approximately 58% of the splashed mass is from the impinging drop itself, while the remainder of the splashed mass is composed of film liquid. From a high resolution case, details of secondary droplet formation are observed and three distinct phases of secondary droplet formation are identified. The detailed analysis of drop train impingement under engine-relevant conditions serves as a first step toward a robust understanding of fuel impingement and the development of cleaner and more efficient ICEs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Wettability and capillary effects: Dynamics of pinch-off in unconstricted straight capillary tubes

We study the interfacial evolution of immiscible two-phase flow within a capillary tube in the partial wetting regime using direct numerical simulation. We investigate the flow patterns resulting from the displacement of a more viscous fluid by a less viscous one under a wide range of wettability conditions. We find that beyond a wettability dependent critical capillary number, a uniform displacement by a less viscous fluid can transition into a growing finger that eventually breaks up into discrete blobs by a series of pinch-off events for both wetting and non-wetting contact angles. This study validates previous experimental observations of pinch-off for wetting contact angles and extends those to non-wetting contact angles. We find that the blob length increases with the capillary number. We observe that the time between consecutive pinch-off events decreases with the capillary number and is greater for more wetting conditions in the displaced phase. Here, we further show that the blob separation distance as a function of the difference between the inlet velocity and the contact line speed collapses into two monotonically decreasing curves for wetting and non-wetting contact angles. For the phase separation in the form of pinch-off, this work provides a quantitative study of the emerging length and time scales and their dependence on the wettability conditions, capillary effects, and viscous forces.

42 ENGINEERING↗
DOE OSTI2018

T2Well/ECO2N Version 1.0: Multiphase and Non-Isothermal Model for Coupled Wellbore-Reservoir Flow of Carbon Dioxide and Variable Salinity Water

T2Well/ECO2N is a coupled wellbore and reservoir model for simulating the dynamics of CO2 injection and leakage through wellbores. It can be seen as an extension to standard TOUGH/ECO2N V2.0, and can be applied to situations relevant to geologic CO2 storage involving upward flow (e.g., leakage) and downward flow (injection). The new simulator integrates a wellbore-reservoir system by assigning the wellbore and reservoir to two different sub-domains in which flow is controlled by appropriate physical laws. In the reservoir, we model flow using a standard multiphase Darcy flow approach. In the wellbores, we use the Drift-Flux Model and related conservation equations for describing transient two-phase non-isothermal wellbore flow of CO2-water mixtures. The mass and thermal energy balance equations are solved numerically by a finite difference scheme with wellbore heat transmission to the surrounding rock handled either semi-analytically or numerically. The momentum balance equation for the flow in the wellbore is solved numerically with a semi-explicit scheme.

Oldenburg, CurtM↗
DOE OSTI2017

Relaxation drag history of shock accelerated microparticles

Experimental measurements of the displacements of shock accelerated microparticles from shortly after shock interaction to the particle relaxation time show time-dependent drag coefficients ( $C_{D}$ ) that are much higher than those predicted by quasi-steady and unsteady drag models. Nylon particles with mean diameter of $4~\unicode[STIX]{x03BC}\text{m}$ , accelerated by one-dimensional normal shocks (Mach number $M_{s}=1.2$ , 1.3 and 1.4), have measured $C_{D}$ values that follow a power-law behaviour. The drag is a function of the time-dependent Knudsen number, $Kn^{\ast }=M_{s}/Re_{p}$ , where the particle Reynolds number ( $Re_{p}$ ) is calculated using the time-dependent slip velocity. Also, some portion of the drag can be attributed to quasi-steady forces, but the total drag cannot be predicted by current unsteady force models that are based on the Basset–Boussinesq–Oseen equation and pressure drag. The largest contribution to the total drag is the unsteady component ( $C_{D,us}$ ) until the particle attains $Kn^{\ast }\approx 0.5{-}1.0$ , then the unsteady contribution decays. The quasi-steady component ( $C_{D,qs}$ ) increases almost linearly with $Kn^{\ast }$ , intersects the $C_{D,us}$ at $Kn^{\ast }\approx 2$ and becomes the primary contributor to the drag towards the end of the relaxation zone as $Re_{p}\rightarrow 0$ . Finally, there are currently no analytical models that are able to predict the nonlinear behaviour of the shock accelerated particles during the relaxation phase of the flow.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Machine learning accelerated turbulence modeling of transient flashing jets

We report modeling the sudden depressurization of superheated liquids through nozzles is a challenge because the pressure drop causes rapid flash boiling of the liquid. The resulting jet usually demonstrates a wide range of structures, including ligaments and droplets, due to both mechanical and thermodynamic effects. As the simulation comprises increasingly numerous phenomena, the computational cost begins to increase. One way to moderate the additional cost is to use machine learning surrogacy for specific elements of the calculation. This study presents a machine learning-assisted computational fluid dynamics approach for simulating the atomization of flashing liquids accounting for distinct stages, from primary atomization to secondary breakup to small droplets using the Σ - Y model coupled with the homogeneous relaxation model. Notably, the models for thermodynamic non-equilibrium (HRM) and Σ - Y are coupled, for the first time, with a deep neural network that simulates the turbulence quantities, which are then used in the prediction of superheated liquid jet atomization. The data-driven component of this method is used for turbulence modeling, avoiding the solution of the two-equation turbulence model typically used for Reynolds-averaged Navier-Stokes simulations for these problems. Both the accuracy and speed of the hybrid approach are evaluated, demonstrating adequate accuracy and at least 25% faster computational fluid dynamics simulations than the traditional approach. This acceleration suggests that perhaps additional components of the calculation could be replaced for even further benefit. Published under an exclusive license by AIP Publishing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Pore occupancy of gas hydrate

Methane hydrate deposits are one of the largest fractions of hydrocarbons in the Earth's crust. They are found mainly in ocean sediments, and the configuration of the deposits—at the largest (kilometer) and smallest (micrometer) scales—determines how and when gaseous methane is released. Here, in this work, using thermodynamic arguments, we show how the confined spaces of the sediment allow methane hydrate to coexist with both water and gas to create a three-phase region. We find that the hydrate pore occupancy changes depending on the depth of the three-phase region, due to changes in methane density of the gas phase. We estimate the thickness of three-phase regions (hydrate, gas, and water) within ocean sediments. We further show how this directly predicts how the presence of hydrate affects the flow properties of gaseous methane and water through porous media.

flows in porous media↗
DOE OSTI2025

Competing mechanisms at vibrated interfaces of density-contrast fluids

Fluid-fluid interfacial instability and subsequent fluid mixing are ubiquitous in nature and engineering. The hydrodynamic instability of fluid interfaces has long centered on the pressure gradient-driven long-wavelength Rayleigh-Taylor instability and the resonance-induced short-wavelength Faraday instability. However, neither instability alone can explain the dynamics when both mechanisms are present. We identify a previously unseen multi-modal instability emerging from their coexistence. When the denser fluid is polydimethylsiloxane, the mixed region at a high density contrast (Atwood number = 0.9) spans a vibration amplitude range approximately twice the gravitational acceleration. Using Floquet stability analysis, we show how vibrations govern transitions between the RT and Faraday instabilities, leading to contention between these instabilities rather than resonant enhancement. Here, the initial transient growth is represented by the exponential modal growth of the most unstable Floquet exponent, along with its accompanying periodic behavior. Direct numerical simulations validate these findings and track interface breakup into the multiscale and nonlinear regimes. Specifically, we show that growing RT modes nonlinearly suppresses Faraday responses even when the initial growth rate of the Faraday instability is 3.63 times that of RT, so a bidirectional competition hinders their sustained coexistence.

Direct numerical simulations↗
DOE OSTI2018

Fluid Flow, Dissolution, and Mixing Phenomena in Argon-Stirred Steel Ladles

In the current study, the multiphase fluid flow in argon-stirred steel ladles is simulated using an Eulerian–Lagrangian two-phase approach. The momentum source and the turbulent kinetic energy source due to the motion of the bubble are considered for the liquid phase. Argon bubbles are treated as discrete phase particles, and the interfacial forces between the bubbles and the liquid phase; the dependence of the gas density and the bubble diameter on the temperature and the static pressure; and the bubble size distribution are considered. When the fluid flow reaches the quasi-steady state, the ferroalloy melting and mixing phenomena is also modeled. The melting time and the trajectory length of each ferroalloy particle are recorded using a user-defined function (UDF). Local mixing time is predicted in the entire computational domain by checking the mixing criteria in every cell. The effects of gas flow rate, porous plug location, and separation angle of two porous plugs on the fluid flow and the mixing phenomena are investigated. The results show that the flow intensity increases, and the mixing time decreases with the increasing gas flow rate. The optimal porous plug’s radial position with one porous plug is 0.50R for its best mixing condition. When two porous plugs are adopted, the separation angle of 90 deg is recommended to improve the flow field and mixing phenomena.

42 ENGINEERING↗
DOE OSTI2017

A positive and entropy-satisfying finite volume scheme for the Baer–Nunziato model

We present a relaxation scheme for approximating the entropy dissipating weak solutions of the Baer–Nunziato two-phase flow model. This relaxation scheme is straightforwardly obtained as an extension of the relaxation scheme designed in for the isentropic Baer–Nunziato model and consequently inherits its main properties. To our knowledge, this is the only existing scheme for which the approximated phase fractions, phase densities and phase internal energies are proven to remain positive without any restrictive condition other than a classical fully computable CFL condition. For ideal gas and stiffened gas equations of state, real values of the phasic speeds of sound are also proven to be maintained by the numerical scheme. It is also the only scheme for which a discrete entropy inequality is proven, under a CFL condition derived from the natural sub-characteristic condition associated with the relaxation approximation. This last property, which ensures the non-linear stability of the numerical method, is satisfied for any admissible equation of state. We provide a numerical study for the convergence of the approximate solutions towards some exact Riemann solutions. The numerical simulations show that the relaxation scheme compares well with two of the most popular existing schemes available for the Baer–Nunziato model, namely Schwendeman–Wahle–Kapila's Godunov-type scheme and Tokareva–Toro's HLLC scheme . The relaxation scheme also shows a higher precision and a lower computational cost (for comparable accuracy) than a standard numerical scheme used in the nuclear industry, namely Rusanov's scheme. Finally, we assess the good behavior of the scheme when approximating vanishing phase solutions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Evaluation of thermodynamic closure models for partially reacted two-phase mixture of condensed phase explosives

One of the key fundamental issues that is crucial in the continuum modeling of reactive flow phenomena is the thermodynamically consistent description of reaction mixture properties. To define the mixture properties, thermodynamic closure rules that relate the properties of the individual reaction components to the mixture properties are required. In the context of reactive two-phase modeling approaches, various strategies to define the thermodynamic closures have been adopted such as pressure temperature (PT) equilibrium between the individual reaction components, pressure (specific) volume (PV) equilibrium, etc. The choice of closure rules determines the relative distribution of specific volume and energy across the reaction components that comprise the mixture. Therefore, depending on the choice of the closure, the mixture thermodynamic behavior can vary. The present work examines the effect of different closure approaches on the thermodynamic properties of the reaction mixture. The analysis is performed for a condensed phase HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine) based plastic bonded explosive (PBX) 9501 explosive using four different thermodynamic closures, viz., PT equilibrium, PV equilibrium, volume temperature (VT) equilibrium, and pressure (P) equilibrium with reactants on an isentrope. Furthermore, the relative variations in the thermodynamic properties of the mixture are analyzed and compared under both compression and expansion loading regimes. It is shown that out of the four closure models, only PT equilibrium and P equilibrium closures lead to a thermodynamically accurate description of the mixture under both compression and expansion.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1976

Non-equilibrium radiation from viscous chemically reacting two-phase exhaust plumes

A knowledge of the structure of the rocket exhaust plumes is necessary to solve problems involving plume signatures, base heating, plume/surface interactions, etc. An algorithm is presented which treats the viscous flow of multiphase chemically reacting fluids in a two-dimensional or axisymmetric supersonic flow field. The gas-particle flow solution is fully coupled with the chemical kinetics calculated using an implicit scheme to calculate chemical production rates. Viscous effects include chemical species diffusion with the viscosity coefficient calculated using a two-equation turbulent kinetic energy model.

Penny, M. M.↗