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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Liquid and liquid–gas flows at all speeds

All speed flows and in particular low Mach number flow algorithms are addressed for the numerical approximation of the Kapila et al. [1] multiphase flow model. This model is valid for fluid mixtures evolving in mechanical equilibrium but out of temperature equilibrium and is efficient for material interfaces computation separating miscible and non-miscible fluids. In this context, the interface is considered as a numerically diffused zone, captured as well as all present waves (shocks, expansion waves). The same flow model can be used to solve cavitating and boiling flows [2]. Many applications occurring with liquid–gas interfaces and cavitating flows involve a very wide range of Mach number, from 10{sup −3} to supersonic (and even hypersonic) conditions with respect to the mixture sound speed. It is thus important to address numerical methods free of restrictions regarding the Mach number. To do this, a preconditioned Riemann solver is built and embedded into the Godunov explicit scheme. It is shown that this method converges to exact solutions but needs too small time steps to be efficient. An implicit version is then derived, first in one dimension and second in the frame of 2D unstructured meshes. Two-phase flow preconditioning is then addressed in the frame of the Saurel et al. [3] algorithm. Modifications of the preconditioned Riemann solver are needed and detailed. Convergence of both single phase and two-phase numerical solutions are demonstrated with the help of single phase and two-phase steady nozzle flow solutions. Last, the method is illustrated by the computation of real cavitating flows in Venturi nozzles. Vapour pocket size and instability frequencies are reproduced by the model and method without using any adjustable parameter.

42 ENGINEERING↗
DOE OSTI2022

Experimental measurements of gas-particle flows in large-scale strippers

Measurements of gas–solid particle flow in a large-scale stripper unit are reported. The experiments target industrial (pilot) scale measurements of gas-particle flow that can be used to validate numerical methods for modeling multiphase flows, such as computational fluid dynamics coupled to the discrete element method (CFD-DEM). Specifically, experiments were performed with Geldart Group B (533 µm) glass beads in a stripper unit. In general, previous CFD-DEM validation datasets were performed in benchtop experimental systems and limited to O(10 5 ) particles. Here, in this work, the number of particles in the isolated section of the stripper is O(10 7 ) for all tests. Measurements are reported for a ~2 m tall section of the ~1 m diameter stripper. Radial mass flux profiles measured at five axial positions and three axial pressure drop measurements are reported for six operating conditions.

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

Simulation of Boiling Flow in a Helical Coil Using a Homogeneous Equilibrium Model with the Nek-2P CFD Code

Recent developments in the field of nuclear energy have lead to an increased interest in the design of helically coiled steam generators. In particular, their compact size, and high heat transfer rates make them ideal for small modular reactors. This has lead to a renewed effort to investigate the characteristics of fluid flow and heat transfer in and around helical coils both numerically and experimentally. The focus of this work is on the development of a consistent modeling, formulation, and simulation framework for boiling flows in helical coils using the highly scalable Nek-2P CFD code. As a first step, two phase flow is accounted for using a homogeneous equilibrium model. Results for both single and multiphase flows are compared to results obtained from Star-CCM+. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2018

Performance Assessment of a Generic Nuclear Waste Repository in Shale - 18392

A performance assessment (PA) is an important component of a comprehensive safety analysis for a nuclear waste repository. In a PA, probabilistic simulations of the total repository system are performed, and results are evaluated against performance metrics. Uncertainty and sensitivity analyses help prioritize additional research and model development. The United States Department of Energy (DOE) has been developing a state of the art PA simulation software tool-kit, Geologic Disposal Safety Assessment (GDSA) Framework, that couples increasingly higher fidelity models of subsystem processes into total system PA simulations. Over the past several years, PAs of generic repositories in three geologic media (salt, shale, and crystalline rock) have demonstrated ongoing developments in capability. The current PA of a nuclear waste repository in a generic shale formation showcases GDSA Framework, including capabilities in domain discretization (Cubit), multi-physics simulations (PFLOTRAN), uncertainty and sensitivity analysis (Dakota), and visualization (Paraview). The generic shale reference case considers the disposal of 22,000 metric tons heavy metal of commercial spent nuclear fuel (SNF) in a generic shale formation. PA simulations account for the thermal load and radionuclide inventory of the waste form, components of the engineered barrier system including waste package and buffer, and components of the natural barrier system including the host rock shale and underlying and overlying stratigraphic units. Two repository layouts are considered, one for emplacement of waste packages containing 12 pressurized water reactor (PWR) assemblies, and one for 4-PWR waste package emplacement. Model domains are half-symmetry and contain between 7 and 22 million grid cells. Grid refinement captures the detail of individual waste packages, emplacement drifts, access drifts, and shafts. Simulations are run in a high-performance computing (HPC) environment on 512 or 2048 processes. The governing equations describing coupled heat and fluid flow and reactive transport are solved with PFLOTRAN, an open-source, massively parallel multiphase flow and reactive transport code. Additional simulated processes include waste package degradation; waste form dissolution; radioactive decay and ingrowth in aqueous, solid, and adsorbed phases; and a simple biosphere based on a pumping well. Simulations are run to 106 y; dose and radionuclide concentrations are observed within aquifers at a point approximately 5 km downgradient of the repository. Dakota is used to sample likely ranges of input parameters including waste form and waste package degradation rates and properties of engineered and natural materials to quantify uncertainty in predicted concentrations and sensitivity to input parameters. Given the assumptions of the reference case, results of 12-PWR and 4-PWR simulations are very similar. Dose at the pumping well 5-km downgradient of the repository does not exceed 10-9 Sv/y. I-129 is the largest contributor to dose. The dose due to Cl-36, the next largest contributor, is orders of magnitude smaller. At the pumping well, the concentration of I-129 5-km downgradient does not exceed 10-14 mol/L. Rank correlation coefficients indicate that the maximum concentration of I-129 within the aquifer is sensitive to shale (repository host rock) porosity and aquifer permeability. These results affirm that HPC-capable codes can be used to simulate important multi-physics couplings directly in a total system performance assessment of a geologic repository, and can be used in prioritization of future research and model development. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2014

A comparative study of lock-in-amplifiers and improved duffing chaotic oscillators for the detection and processing of weak signals

The measurement of multiphase flow parameters is of great importance in a wide range of industries. In the measurement of multiphase, the signals from the sensors are extremely weak and often buried in strong background noise. It is thus desirable to develop effective signal processing techniques that can detect the weak signal from the sensor outputs. In this paper, two methods, i.e., lock-in-amplifier (LIA) and improved Duffing chaotic oscillator are compared to detect and process the weak signal. For sinusoidal signal buried in noise, the correlation detection with sinusoidal reference signal is simulated by using LIA. The improved Duffing chaotic oscillator method, which based on the Wigner transformation, can restore the signal waveform and detect the frequency. Two methods are combined to detect and extract the weak signal. Simulation results show the effectiveness and accuracy of the proposed improved method. The comparative analysis shows that the improved Duffing chaotic oscillator method can restrain noise strongly since it is sensitive to initial conditions.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

A Scalable Multigrid Reduction Framework for Multiphase Poromechanics of Heterogeneous Media

Simulation of multiphase poromechanics involves solving a multiphysics problem in which multiphase flow and transport are tightly coupled with the porous medium deformation. To capture this dynamic interplay, fully implicit methods, also known as monolithic approaches, are usually preferred. The main bottleneck of a monolithic approach is that it requires solution of large linear systems that result from the discretization and linearization of the governing balance equations. Because such systems are nonsymmetric, indefinite, and highly ill-conditioned, preconditioning is critical for fast convergence. Recently, most efforts in designing efficient preconditioners for multiphase poromechanics have been dominated by physics-based strategies. Current state-of-the-art “black-box” solvers such as algebraic multigrid (AMG) are ineffective because they cannot effectively capture the strong coupling between the mechanics and the flow subproblems, as well as the coupling inherent in the multiphase flow and transport process. In this work, we develop an algebraic framework based on multigrid reduction (MGR) that is suited for tightly coupled systems of PDEs. Using this framework, the decoupling between the equations is done algebraically through defining appropriate interpolation and restriction operators. One can then employ existing solvers for each of the decoupled blocks or design a new solver based on knowledge of the physics. We demonstrate the applicability of our framework when used as a “black-box” solver for multiphase poromechanics. Here, we show that the framework is flexible to accommodate a wide range of scenarios, as well as efficient and scalable for large problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

Smoothed Particle Hydrodynamics pore-scale simulations of unstable immiscible flow in porous media

We have conducted a series of high-resolution numerical experiments using the Pair-Wise Force Smoothed Particle Hydrodynamics (PF-SPH) multiphase flow model. First, we derived analytical expressions relating parameters in the PF-SPH model to the surface tension and static contact angle. Next, we used the model to study viscous fingering, capillary fingering, and stable displacement of immiscible fluids in porous media for a wide range of capillary numbers and viscosity ratios. We demonstrated that the steady state saturation profiles and the boundaries of viscous fingering, capillary fingering, and stable displacement regions compare favorably with micromodel laboratory experimental results. For displacing fluid with low viscosity, we observed that the displacement pattern changes from viscous fingering to stable displacement with increasing injection rate. When a high viscosity fluid is injected, transition behavior from capillary fingering to stable displacement occurred as the flow rate was increased. These observation also agree with the results of the micromodel laboratory experiments.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2024

Effect of the volume fraction gradient on the phase interaction force model for disperse two-phase flows

In this work, the effects of the particle volume fraction gradient on fluid-particle interactions are studied. The phase interaction force is decomposed into three terms. For the first term, namely the symmetrized force density, we present theoretical reasoning and numerical evidence to assume that it is independent of the particle volume fraction gradient. The second term is the particle volume fraction gradient times a newly introduced diffusion stress. The third term is the divergence of the particle-fluid-particle (PFP) stress. If this assumption of independence of the particle volume fraction gradient for the first term can be verified, to the first order of the ratio of the mean distance between particles to the macroscopic lengthscale, all three terms can be studied and modeled in flows with uniform particle distributions. Models thus obtained are applicable to statistically inhomogeneous flows, with the second and third terms accounting for statistical inhomogeneity. To verify this assumption, numerical simulations of flows passing fixed arrays of particles are performed. Both uniform and nonuniform particle volume fractions are studied and compared for disperse multiphase flows with the particle Reynolds numbers ranging from 1 to 100, and particle volume fraction ranging from 1% to 26% in statistically steady states. It is found that the symmetrized force (first) term can be well approximated by the drag force obtained from studies of uniform flows. The diffusion stress is positive along the flow direction and negative in the directions perpendicular to the flow. In the case of moving particles, this stress could potentially cause particle aggregation in the flow direction and dispersion in the directions perpendicular to the flow. Finally, the diffusion stress is only important when there is a volume fraction gradient, while the PFP stress can be important in inhomogeneous flows with either nonuniform particle concentrations or nonuniform average relative velocities between the phases.

42 ENGINEERING↗
DOE OSTI2018

Methane Hydrate Formation in Thick Sandstones by Free Gas Flow

We show with a two-dimensional multiphase flow and multicomponent transport model that free gas flow is a viable mechanism to form concentrated methane hydrate in meter-scale, dipping sandstones far above the base of the hydrate stability zone (BHSZ). In this model, gas preferentially flows updip along the top of sandstone due to buoyancy. This drives hydrate formation, increasing the local salinity to the stability limit and developing three-phase (gas, liquid, and hydrate) equilibrium above the BHSZ. With time, the gas and the hydrate solidification front (HSF) advance together updip. Behind the HSF, hydrate continues to form as the elevated salinity diffuses away. High hydrate saturations reduce the sediment permeability significantly. As a result, as the gas and HSF move updip, they are also pushed perpendicularly from the top to the base of the sandstone. The hydrate system ultimately self-seals itself due to reduced permeability across the entire thickness of the sandstone. Gas starts to retreat downdip and accumulates below the BHSZ. With this model high hydrate saturations form far above the BHSZ at high methane supply rates while hydrate is concentrated at the BHSZ or no hydrate forms at low methane supply rates. This study provides further insights into hydrate formation by free gas flow, which can be used to design the best strategies for economic and environmental production of methane from hydrate reservoirs.

58 GEOSCIENCES↗
DOE OSTI2018

Addressing Two-Phase Flow Maldistribution in Microchannel Heat and Mass Exchangers

Multiphase flow phenomena in single micro and minichannels have been widely studied. Characteristics of two-phase flow through a large array of microchannels are investigated here. An air–water mixture is used to represent the two phases flowing through a microchannel array representative of those employed in practical applications. Flow distribution of the air and water flow across 52 parallel microchannels of 0.4 mm hydraulic diameter is visually investigated using high-speed photography. Two microchannel configurations are studied and compared, with mixing features incorporated into the second configuration. Slug and annular flow regimes are observed in the channels. Void fractions and interfacial areas are calculated for each channel from these observations. The flow distribution is tracked at various lengths along the microchannel array sheets. Statistical distributions of void fraction and interfacial area along the microchannel array are measured. The design with mixing features yields improved flow distribution. Void fraction and interfacial area change along the length of the second configuration, indicating a change in fluid distribution among the channels. The void fraction and interfacial area results are used to predict the performance of different microchannel array configurations for heat and mass transfer applications. Results from this study can help inform the design of compact thermal-fluid energy systems.

Engineering↗
DOE OSTI2022

Implementing General Framework in MFIX for Radiative Heat Transfer in Gas-Solid Reacting Flows (Final Report)

The current knowledge gap between radiation and multiphase flow communities is the greatest barrier to developing predictive models for radiation, the dominant mode of heat transfer in gas-solid reacting flows. The objectives of this research are (1) to develop and implement a general framework to support the integration of modern radiation models for gas-solid reacting flows; (2) to implement a methodology for developing new multiphase radiation models with accuracy and efficiency commensurate with their importance in energy-related applications, and (3) to demonstrate the accuracy and efficiency of the radiation models under typical gas-solid reacting flow conditions.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2024

Hydrogen, Methane, Brine Flow Behavior, and Saturation in Sandstone Cores During H 2 and CH 4 Injection and Displacement

Large-scale underground hydrogen storage (UHS) is a critical component in the emerging hydrogen economy. Knowledge of multiphase flow behavior involving hydrogen in storage reservoir formations is crucial to characterizing hydrogen transport properties and essential for the deliverability and storage operations of UHS. There are still many gaps in fully understanding hydrogen–methane–brine multiphase phase flow that require further investigation. In this work, H 2 and CH 4 were injected through brine-saturated sandstone cores using a tri-axial core holder system fitted with flow rate meters and pressure transducers, while the effluent gas concentrations were analyzed using an online micro gas chromatograph. Brine displacement, permeability, and gas breakthrough curves were measured. We studied the flow behavior of hydrogen and methane in sandstone cores through testing brine displacement by gas injection and comparing the hydrogen displacement of methane with the methane displacement of hydrogen. We also tested the differences between horizontal and vertical flow in brine displacement. The results showed that brine displacement was more efficient in a core with higher permeability and porosity, resulting in a higher initial gas saturation. A higher gas injection rate brought about faster gas breakthrough measured by pore volume and sharper concentration curves. Hydrogen did not exhibit abnormal flow in the sandstone when the flow was horizontal and downward vertical. Gas overriding was observed in brine displacements when the flow was horizontal, with hydrogen showing this behavior more profoundly compared to methane. Downward vertical gas injection induced higher efficiency brine displacement compared to horizontal displacement and resulted in a higher initial gas saturation in the sandstone cores. These findings address critical knowledge gaps regarding gas flow patterns and displacement behaviors during hydrogen injection and recovery phases in UHS facilities using methane as the cushion gas. The insights from this research offer valuable guidance for optimizing UHS systems, ensuring operational efficiency, and advancing sustainable energy solutions in alignment with decarbonization goals.

08 HYDROGEN↗
DOE OSTI2015

Smoothed dissipative particle dynamics with angular momentum conservation

Smoothed dissipative particle dynamics (SDPD) combines two popular mesoscopic techniques, the smoothed particle hydrodynamics and dissipative particle dynamics (DPD) methods, and can be considered as an improved dissipative particle dynamics approach. Despite several advantages of the SDPD method over the conventional DPD model, the original formulation of SDPD by Español and Revenga (2003) [9], lacks angular momentum conservation, leading to unphysical results for problems where the conservation of angular momentum is essential. To overcome this limitation, we extend the SDPD method by introducing a particle spin variable such that local and global angular momentum conservation is restored. The new SDPD formulation (SDPD+a) is directly derived from the Navier–Stokes equation for fluids with spin, while thermal fluctuations are incorporated similarly to the DPD method. We test the new SDPD method and demonstrate that it properly reproduces fluid transport coefficients. Also, SDPD with angular momentum conservation is validated using two problems: (i) the Taylor–Couette flow with two immiscible fluids and (ii) a tank-treading vesicle in shear flow with a viscosity contrast between inner and outer fluids. For both problems, the new SDPD method leads to simulation predictions in agreement with the corresponding analytical theories, while the original SDPD method fails to capture properly physical characteristics of the systems due to violation of angular momentum conservation. In conclusion, the extended SDPD method with angular momentum conservation provides a new approach to tackle fluid problems such as multiphase flows and vesicle/cell suspensions, where the conservation of angular momentum is essential.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Shock-driven dispersal of a corrugated finite-thickness particle layer

A research area emerging in the multiphase flow community is the study of shock-driven multiphase instability (SDMI), a gas–particle analog of the traditional fluid-fluid Richtmyer–Meshkov instability (RMI). In this work, we study the interaction of planar air shocks with corrugated glass particle curtains through the use of numerical simulations with an Eulerian–Lagrangian approach. One objective of this study is to compare the simulated particle curtains to a comparable set of shock tube experiments performed to analyze traditional RMI of a gas curtain. The simulations are set to match the experimental shock Mach numbers and perturbation wavelengths (3.6 and 7.2 mm) while also matching the Atwood number of the experiments to the multiphase Atwood number of the simulations. Varying particle diameters are tested in the simulations to explore the impact of particle diameter on the evolution of the particle curtain. This simulation setup allows for a one-to-one comparison between RMI and SDMI under comparable conditions while also allowing for a separate study into the validity of the use of the multiphase Atwood number to compare the single-phase and multiphase instabilities. In particular, we show that the comparison depends on the diameter of the particles (thus, dependent on the Stokes number of the flow). A second objective of this study is to analyze the effect of the initial particle volume fraction on the evolution of the curtain and the behavior of the instability. This is done through analyzing the effect of the multiphase terms of the vorticity evolution equation on the vorticity deposition in SDMI. Also discussed is the effect of the particle diameter on the multiphase generation terms as well as in the baroclinic vorticity generation term in SDMI as the shock passes over the curtain.

42 ENGINEERING↗
DOE OSTI2021

Tomographic X-ray particle tracking velocimetry

We investigate the feasibility of in-laboratory tomographic X-ray particle tracking velocimetry (TXPTV) and consider creeping flows with nearly density matched flow tracers. Specifically, in these proof-of-concept experiments we examined a Poiseuille flow, flow through porous media and a multiphase flow with a Taylor bubble. For a full 360° computed tomography (CT) scan we show that the specially selected 60 micron tracer particles could be imaged in less than 3 seconds with a signal-to-noise ratio between the tracers and the fluid of 2.5, sufficient to achieve proper volumetric segmentation at each time step. In the pipe flow, continuous Lagrangian particle trajectories were obtained, after which all the standard techniques used for PTV or PIV (taken at visible wave lengths) could also be employed for TXPTV data. And, with TXPTV we can examine flows inaccessible with visible wave lengths due to opaque media or numerous refractive interfaces. In the case of opaque porous media we were able to observe material accumulation and pore clogging, and for flow with Taylor bubble we can trace the particles and hence obtain velocities in the liquid film between the wall and bubble, with thickness of liquid film itself also simultaneously obtained from the volumetric reconstruction after segmentation. While improvements in scan speed are anticipated due to continuing improvements in CT system components, we show that for the flows examined even the presently available CT systems could yield quantitative flow data with the primary limitation being the quality of available flow tracers.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1991

Theoretical and computational dynamics of a compressible flow

An introduction to the theoretical and computational fluid dynamics of a compressible fluid is presented. The general topics addressed include: thermodynamics and physical properties of compressible fluids; 1D flow of an inviscid compressible fluid; shock waves; fundamental equations of the dynamics of a compressible inviscid non-heat-conducting and radiating fluid, method of small perturbations, linearized theory; 2D subsonic steady potential flow; hodograph and rheograph methods, exact solutions of 2D insentropic steady flow equations, 2D steady transonic and hypersonic flows, method of characteristics, linearized theory of 3D potential flow, nonlinear theory of 3D compressibe flow, anisentropic (rotational) flow of inviscid compressible fluid, electromagnetogasdynamics, multiphase flows, flows of a compressible fluid with transport phenomena.

Pai, Shih-I↗