Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiphase fluid 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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Assessment of model parameters in MFiX particle-in-cell approach

The limitations in numerical treatment of solids-phase in conventional methods like Discrete Element Model and Two-Fluid Model have facilitated the development of alternative techniques such as Particle-In-Cell (PIC). However, a number of parameters are involved in PIC due to its empiricism. In this work, global sensitivity analysis of PIC model parameters is performed under three distinct operating regimes common in chemical engineering applications, viz. settling bed, bubbling fluidized bed and circulating fluidized bed. Simulations were performed using the PIC method in Multiphase Flow with Interphase eXchanges (MFiX) developed by National Energy Technology Laboratory (NETL). A non-intrusive uncertainty quantification (UQ) based approach is applied using Nodeworks to first construct an adequate surrogate model and then identify the most influential parameters in each case. This knowledge will aid in developing an effective design of experiments and determine optimal parameters through techniques such as deterministic or statistical calibration.

01 COAL, LIGNITE, AND PEAT↗
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 OSTI2024

Machine Learning-Enhanced Multiphase CFD for Carbon Capture Modeling Run Data

Repository for the data generated as part of the 2023-2024 ALCC project "Machine Learning-Enhanced Multiphase CFD for Carbon Capture Modeling." The data was generated with MFIX-Exa's CFD-DEM model. The problem of interest is gravity driven, particle-laden, gas-solid flow in a triply-periodic domain of length 2048 particle diameters with an aspect ratio of 4. The mean particle concentration ranges from 1% to 40% and the Archimedes number ranges from 18 to 90. The particle-to-fluid density ratio, particle-particle restitution and friction coefficients and domain aspect ratio are held constant at values of 1000, 0.9, 0.25 and 4, respectively. This research used resources of the National Energy Research Scientific Computing Center, a DOE Office of Science User Facility supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231 using NERSC award ALCC-ERCAP0025948.

AMReX↗
DOE OSTI2021

A Physics-Constrained Deep Learning Model for Simulating Multiphase Flow in 3D Heterogeneous Porous Media

Physics-based simulators for multiphase flow in porous media emulate nonlinear processes with coupled physics, and usually require extensive computational resources for software development, maintenance and simulation execution. As a result, a huge demand exists for fast modeling of coupled processes in a wide range of subsurface applications including geological sequestration, hydrocarbon recovery and geothermal energy extraction. In this work, an efficient physics-constrained deep learning model is developed for solving multiphase flow in 3-Dimensional (3D) heterogeneous porous media. The model fully leverages the spatial topology predictive capability of convolutional neural networks, specifically U-Net with successive contracting and expansive steps, and is coupled with an efficient continuity-based smoother to predict flow responses that need spatial continuity. Furthermore, the transient regions are penalized to steer the training process such that the model can accurately capture flow in these regions. The model takes inputs including properties of porous media, fluid properties and well controls, and predicts the temporal-spatial evolution of the state variables (pressure and saturation). While maintaining the continuity of fluid flow, the 3D spatial domain is decomposed into 2D images for reducing training cost, and the decomposition results in an increased number of training data samples and better training efficiency. Additionally, a surrogate model is separately constructed as a postprocessor to calculate well flow rate based on the predictions of state variables from the deep learning model. We use the example of CO 2 injection into saline aquifers, and apply the physics-constrained deep learning model that is trained from physics-based simulation data and emulates the physics process. The model performs prediction with a speedup of ~ 1400 times compared to physics-based simulations, and the average temporal errors of predicted pressure and saturation plumes are 0.27% and 0.099% respectively. Furthermore, water production rate is efficiently predicted by a surrogate model for well flow rate, with a mean error less than 5%. Therefore, with its unique scheme to cope with the fidelity in fluid flow in porous media, the physics-constrained deep learning model can become an efficient predictive model for computationally demanding inverse problems or other coupled processes.

58 GEOSCIENCES↗
DOE OSTI2021

Flow reduction of hydrocarbon liquid in silica nanochannel: Insight from many-body dissipative particle dynamics simulations

A modified many-body dissipative particle dynamics (mDPD) model recently developed for realistic mesoscale multiphase flow simulations is rigorously parameterized, calibrated, and applied for elucidating the flow mechanisms of hydrocarbon liquids (i.e., heptane in this work) in amorphous silica cylindrical nanochannels with inner diameters ranging from 4.5 to 22.5 nm. The simulation results suggest the presence of a strong threshold of pressure gradient under which heptane cannot be driven to flow. The threshold for the 4.5 nm diameter pore is 10 to 100 times as high as for the 9–22.5 nm diameter pore, highlighting a remarkable nanoconfinement effect. Fluid viscosity is found to exhibit a shear-thinning phenomenon with intensity to weaken with increasing channel diameter — a phenomenon not observed in nanochannel flow of liquid water and gas in literature. Most remarkably, the radial profiles of average longitudinal flow velocity fitted by the modified Hagen-Poiseuille equation showed a negative slip length (–2.5% to –0.5% relative to the diameter) and a reduction of apparent permeability by 16% to 23%. Furthermore, this finding suggests silica nanochannels tend to deter hydrocarbon flow, a phenomenon that is opposed to the flow enhancement reported in most of the prior nanochannel flow studies in literature.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Diffusion-to-Imbibition Transition in Water Sorption in Nanoporous Media: Theoretical Studies

The ability to predict multiphase fluid transport in nanoporous rocks such as shales is critical for many geoscience applications, for example unconventional hydrocarbon production, geologic carbon sequestration, and nuclear waste disposal. When the pore sizes approach nanoscales, the impact of the molecular interaction forces between fluids and solids becomes increasingly important. These forces can alter macroscopic fluid phase behavior and control transport. Recent experimental studies have shown that capillary condensation and subsequent imbibition of liquid water can occur in hydrophilic nanoporous media even if the vapor phase is at a critical relative humidity (rhcrit) well below vapor saturation. This study presents a theoretical investigation of the processes controlling adsorption, capillary condensation and imbibition in nanoporous media, using the square-gradient classical density functional theory. The proposed theoretical model explicitly includes the relevant interaction forces among fluids and solids in macroscopic porous media. Application of the model to a relative-humidity-controlled water adsorption experiment is presented to demonstrate the impact of water-pore wall attractive forces on multiphase water behavior in a hydrophilic silicon nanoporous medium. Here, the model represents well the measured time-dependent evolution of the water imbibition front inside the nanoporous medium and also explains the diffusion-like water transport regimes observed at rh < rh crit and the imbibition-like flow regimes observed at rh > rh crit . The study furthermore gives an insight on hysteresis phenomenon in adsorption and desorption isotherms.

58 GEOSCIENCES↗
DOE OSTI2023

Solving multiphysics-based inverse problems with learned surrogates and constraints

Abstract Solving multiphysics-based inverse problems for geological carbon storage monitoring can be challenging when multimodal time-lapse data are expensive to collect and costly to simulate numerically. We overcome these challenges by combining computationally cheap learned surrogates with learned constraints. Not only does this combination lead to vastly improved inversions for the important fluid-flow property, permeability, it also provides a natural platform for inverting multimodal data including well measurements and active-source time-lapse seismic data. By adding a learned constraint, we arrive at a computationally feasible inversion approach that remains accurate. This is accomplished by including a trained deep neural network, known as a normalizing flow, which forces the model iterates to remain in-distribution, thereby safeguarding the accuracy of trained Fourier neural operators that act as surrogates for the computationally expensive multiphase flow simulations involving partial differential equation solves. By means of carefully selected experiments, centered around the problem of geological carbon storage, we demonstrate the efficacy of the proposed constrained optimization method on two different data modalities, namely time-lapse well and time-lapse seismic data. While permeability inversions from both these two modalities have their pluses and minuses, their joint inversion benefits from either, yielding valuable superior permeability inversions and CO 2 plume predictions near, and far away, from the monitoring wells.

Yin, Ziyi (ORCID:0000000250248771)↗
DOE OSTI2021

A microfluidic study of transient flow states in permeable media using fluorescent particle image velocimetry

Velocity fields in flow in permeable media are of great importance to many subsurface processes such as geologic storage of CO 2 , oil and gas extraction, and geothermal systems. Steady-state flow is characterized by velocity fields that do not change significantly over time. The flow field transitions to a new steady state once it experiences a disturbance such as a change in flow rate or in pressure gradient. This transition is often assumed to be instantaneous, which justifies the expression of constitutive relations as functions of instantaneous phase saturations. This work examines the evolution of velocity fields in a surrogate quasi-2D permeable medium using a microfluidic device, a microscopy system, and a high-speed camera. Tracer particles are injected into the medium along with DI water. The evolution of the velocity field is examined by tracing these particles in the captured images using the standard high-density particle image velocimetry algorithm founded on cross-correlation. The results suggest that the transition between steady states for an incompressible fluid takes a finite and non-negligible amount of time that is independent of the magnitude of the change in pressure gradient. The existence of transient states and the nature of the response during these states are readily interpreted by the principle of least action where flow gradually establishes an optimal configuration such that energy dissipation is 2 minimized. The findings provide evidence against the applicability of the assumption that flowing phases relax instantaneously to their steady states and, hence, against the accuracy of the classical multiphase extension of Darcy's law.

42 ENGINEERING↗
DOE OSTI2021

Acoustic gas volume fraction measurement in a multiphase flowing liquid

Apparatus and methods for the measurement of gas volume fraction of produced oil are described. A first method measures the response of a pipe containing the produced oil excited by a source of vibration in the form of an acoustic frequency chirp containing a linearly varying range of frequencies in the tens of kilohertz range encompassing at least one resonant mode of the pipe. As the gas volume fraction increases, the location of the peak maximum of the measured frequency spectrum responsive to the excitation increases in frequency, and the height of the peak maximum increases, thereby permitting a linear calibration curve to be obtained. A second method measures the response of a pipe containing the produced oil to excitation by a continuous source of vibration having a chosen frequency above those which excite flexural vibrations in the pipe and simultaneously excite acoustic waves in the fluid contained in the pipe, known as the coincidence frequency. Gas present in the fluid will interrupt sound propagation or reverberation, thereby generating fluctuations in the amplitude of the measured vibrations of the pipe. The amplitude fluctuation level provides a measure of the gas volume present inside the pipe. A third method measures the response of a pipe containing the produced oil to excitation by a high-bandwidth, short pulse having a chosen center frequency above the coincidence frequency. Gas present in the fluid will interrupt pulse propagation, thereby generating fluctuations in the amplitude of the measured vibrations of the pipe.

Sinha, Dipen N.↗
DOE OSTI2021

An Investigation of the Effects of Volume Fraction on Drag Coefficient of Non-Spherical Particles Using PR-DNS

Prediction of the drag coefficient is required in gas-particle multiphase flow modeling and simulation. Experimental data and correlations on the fixed-bed system of spherical particles with high volume fractions for various possible arrangements are available in the literature. However, the effect of volume fraction on the drag coefficient of non-spherical particles is not well studied. In solving the momentum equation, the volume fraction plays a vital role in determining the flow resistances. In this paper, we study the impact of volume fraction in the range of 0.069 to 0.65 on the drag coefficient using the computational fluid dynamics (CFD) simulation of air for Reynold number in the range of 10 to 10000 using particle resolved direct numerical solution (PR-DNS). Regular non-spherical particles such as a cube, tetrahedron, and spheroids are used in this study since their single particle’s drag coefficient data are available in the literature for comparison. For this work, the simulations are carried out in the Ansys Fluent using polyhedral mesh, which consumes significantly less computational time and power. The study showed the sphericity and volume fraction have significant impact on the bed pressure drop and average drag coefficient of the particles in the bed especially in high Reynolds number regime. The bed of the spheroid experiences the lowest drag being the most streamlined particle, and the particles with the edges result in a large drag coefficient due to flow separation at the discontinuity. The vector plots verify this behavior where large wake regions are observed behind the tetrahedron particle.

Mahyawansi, Pratik↗
DOE OSTI2024

Development of a multiphase PIC model for slurry flow modeling

MFIX-Exa is a recently released multiphase CFD code originally developed for the simulation of particle-laden gas-solid flows. Due to its high-performance computing capabilities, MFIX-Exa is an ideal candidate for scale-up studies of slurry reactors, specifically the coarse-grained particle-in-cell (PIC) model with its statistical treatment of the particle phase. Unfortunately, several physical models that were neglected during original development because they are not relevant for high-density ratio gas-solid flows are important in slurry flows where the particle-to-fluid density ratio is near unity. In this preliminary work we focus on the effective (suspension) viscosity. The models of Brinkman (1952), Krieger and Dougherty (1956), and Cheng and Law (2003) are considered. The impact of the effective viscosity model is studied on horizontal pipe flow. The experimental data of Gillies et al. (2002) is used to assess the pressure drop predictions.

Fullmer, William D.↗
DOE OSTI2020

A proximity‐based image‐processing algorithm for colloid assignment in segmented multiphase flow datasets

Summary Colloidal transport and deposition are of both environmental and engineering importance. Easier access to x‐ray microtomography (XMT) coupled with improved imaging resolution has made XMT a unique and viable tool for visualizing and quantifying these processes. Currently, there is scant information in the literature addressing colloid segmentation and analysis in saturated and unsaturated porous media, in particular related to spatial partitioning of colloids. To support this need, an approach to assign segmented colloidal particles and aggregates to different partitioning classes based on their proximity to different phases is presented here. The method uses different markers for each attachment site (e.g. wetting‐nonwetting phase interfaces). An example XMT dataset from a drainage experiment is used to demonstrate the efficacy of the image processing algorithms. Flow conditions, and fluid and colloid properties, can thus be compared to the behaviour of colloids within the porous medium. This algorithm can help elucidate colloidal deposition mechanisms and the importance of different attachment sites, explore the importance of fluid properties, as well as the arrangement and shape of the colloids.

BRUECK, C. L.↗
DOE OSTI2026

Two-phase flow numerical analysis of electrode geometry for alkaline water electrolyzers

Hydrogen is a promising component of a future energy-secure and efficient economy, but its competitiveness depends on reducing production costs. One strategy is to operate alkaline water electrolyzers at higher current densities to increase output. However, this intensifies performance losses due to gas bubble accumulation, which blocks transport pathways and deactivates electrochemically active surfaces. Enhancing bubble evacuation through electrode design is therefore essential. Previous studies have explored various approaches — such as modifying surface morphology, applying sonication or pressure modulation, and introducing surfactants — but these efforts have addressed a limited range of conditions due to the complexity of two-phase flow and electrode geometries. Experiments have also largely been focused on either cell level improvements, which lack the information necessary to isolate each contributing factor, or on modified geometries that are not relevant to practical cell operation. From a modeling perspective, conventional Eulerian multiphase models do not track the complex gas–liquid interfacial dynamics and often neglect surface tension and contact angle effects, reducing their predictive accuracy. To provide insights on the effects of different electrode geometries on the performance of alklaine water electrolyzers this work employs an immersed boundary volume-of-fluid method to simulate bubble behavior in 3D porous electrodes. Multiple base electrode geometries, typically used in practice, with varying porosity are evaluated under a constant surface gas generation rate. Simulation data is analyzed to quantify electrode gas coverage, bubble size dynamics and other relevant metrics. Results show that porosity strongly influences bubble accumulation on electrode surfaces, with higher porosity reducing gas coverage, and its not strictly dependent on the electrode geometry. However, the electrode’s base geometry significantly affects gas accumulation at the separator gap, independent of porosity. A foam electrode geometry resulted in the lowest gas coverage of all electrodes with a median volumetric gas coverage of 11%, but at the cost of a 70% reduction in active area compared with the largest surface area electrode, while gyroid electrodes showed the best trade-off between gas coverage, particularly at the separator surface, and electrochemically active area. In conclusion, the results highlight the need for holistic electrode design strategies.

Alkaline↗
DOE OSTI2023

Effective permeability of an immiscible fluid in porous media determined from its geometric state

Based on the phenomenological extension of Darcy's law, two-fluid flow is dependent on a relative permeability function of saturation only that is process or path dependent with an underlying dependency on pore structure and wettability. For applications (fuel cells to underground CO 2 storage), it is imperative to determine the effective phase permeability relationships where the traditional approach is based on the inverse modeling of time-consuming experiments. The underlying reason is that the fundamental upscaling step from pore to Darcy scale, which links the pore structure of the porous medium to the continuum hydraulic conductivities, is not solved. Herein, we develop an artificial neural network (ANN) that relies on fundamental geometrical relationships to determine the mechanical energy dissipation during creeping immiscible two-fluid flow. Finally, the developed ANN is based on a prescribed set of state variables based on physical insights that predict the effective permeability of 4500 unseen pore-scale geometrical states with R 2 = 0.98.

42 ENGINEERING↗
DOE OSTI2022

Experimental Investigation of the Effect of Atomizing Electrospray Nozzles on the Cooling and Reduced Humidification of Air

Separating sensible and latent cooling have shown the potential to reduce energy consumption in A/C systems for building HVAC applications. Several technologies exist based on vapor-compression refrigeration, enthalpy wheels, chemical adsorption and absorption materials, and mechanical cooling. However, their thermodynamics limits and high energy consumption hinder their deployment in hot and humid climatic regions. This paper focuses on separating sensible and latent cooling by using a new thermodynamics process and in-kind (i.e., non-vapor compression-based) technology. Highly electrically charged water droplets were sprayed in the airflow. These droplets attracted water vapor molecules to their surfaces and promoted condensation. The phenomenon was the result of simultaneous dielectrophoresis and electro-diffusion interactions. Studies in the literature have shown that using multiple capillary electrodes reduced air moisture by up to 5% when using nanometer-size droplets in the spray. Unfortunately, these studies were limited to low airflow rates, and the objective of this paper was to investigate how to scale up this in-kind approach to airflows typical of buildings. In the present paper, droplets of micrometer size were utilized to control the humidity for a 5-cfm flow rate. While this airflow was still low for building applications, it was 100 times fold the airflows in the literature studies. The air was tested at 20°C and at 50 and 80% relative humidity. A two-fluid atomizing nozzle produced the droplets in the spray, and high DC electric potential, up to 25 kV, was used to charge the fine droplets electrically. The air atomizing nozzle with high voltage potential resembled an evaporative cooler process. However, a measurable reduction of the absolute humidity of up to 2% was observed compared to the case of the nozzle with no high voltage potential. The entire device had one small nozzle selected from off-the-shelf components and had less than a 9 cm2 footprint area.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
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 OSTI2025

A sticky situation: The heterogeneity and implications of magma mush unlocking

Magma, partially molten rock below Earth's surface, can be made up of both fluid (melt) and particles (crystals). The relative fraction of the fluid and particles matters for the magma's mobility and ability to reach Earth's surface. Here we conduct physical laboratory experiments to examine the mechanisms through which a magma becomes mobile and the speed at which pure melt and infiltrate a mixture of fluids and particles. The results have implications for understanding magma mobility and the mechanisms that move magma to Earth's surface. Poster presented at the 2025 International Association of Volcanology and Chemistry of Earth's Interior Meeting (IAVCEI).

multiphase flows↗
DOE OSTI2023

Multiphase Modeling in a Parallel Plate Fluidized Bed Receiver for Concentrating Solar Power

A novel high temperature particle solar receiver is developed by using a light trapping planar cavity configuration. As particles fall through the cavity, the concentrated solar radiation warms the boundaries of the receiver and in turn heats the particles. Particles flow through the system, forming a packed bed at the lower end, leaving the system from the bottom at a constant flow rate. Air is introduced to the system as the fluidizing medium to improve particle heat transfer and mixing. A laboratory scale cavity receiver is built and a near IR quartz lamp is used to provide flux to the vertical wall of the heat exchanger. The system is modeled using a continuum two-fluid method. The computational model matches the experimental system size and the particle size distribution is assumed monodisperse. A conduction model that accounts for the effects of solid concentration is implemented, and the heat flux boundary condition matches the experimental setup. Radiative heat transfer is estimated using a widely used correlation during the post-processing step to determine an overall heat transfer coefficient. The model is validated against testing data and achieves less than 30% discrepancy and a heat transfer coefficient greater than 1000 W/m2K.

CSP↗