Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Multiphase”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Implementing contact angle boundary conditions for second-order Phase-Field models of wall-bounded multiphase flows

In the present work, a general formulation is proposed to implement the contact angle boundary conditions for the second-order Phase-Field models, which is applicable to N-phase (N ≥ 2) moving contact line problems. To remedy the issue of mass change due to the contact angle boundary condition, a source term or Lagrange multiplier is added to the original second-order Phase-Field models, which is determined by the consistent and conservative volume distribution algorithm so that the summation of the order parameters and the consistency of reduction are not influenced. To physically couple the proposed formulation to the hydrodynamics, especially for large-density-ratio problems, the consistent formulation is employed. The reduction-consistent conservative Allen-Cahn models are chosen as examples to illustrate the application of the proposed formulation. The numerical scheme that preserves the consistency and conservation of the proposed formulation is employed to demonstrate its effectiveness. Results produced by the proposed formulation are in good agreement with the exact and/or asymptotic solutions. The proposed method captures complex dynamics of moving contact line problems having large density ratios.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

MPLBM-UT: Multiphase LBM library for permeable media analysis

MPLBM-UT is a specialized lattice-Boltzmann library that makes running single- and two-phase flow simulations in porous media accessible to everyone. We provide a suite of tools to pre-process computational domains for simulation, to set up custom boundary conditions, to run simulations, to post-process simulation outputs, and to visualize simulation results and data. All of these tools are easily accessible to users through the mplbm_utils Python package included in and automatically installed with MPLBM-UT. The high-performance, highly parallel library Palabos is used as the solver backend. MPLBM-UT is easily deployed in a variety of systems, from laptops to supercomputer clusters. MPLBM-UT also features multiple examples and benchmark templates that allow for fast prototyping of different porous media problems. We also provide an interface for reading in different file types and downloading domains from the Digital Rocks Portal to perform simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Evolution of the age-included nearest pair distribution in disperse multiphase flows

The age of the nearest particle pair is introduced as the difference between the current time and the most recent time when the nearest particle pair was formed. The evolution equation for the age-included nearest pair distribution function is derived. With the assumption of random destruction of the nearest particle pairs, the evolution equation predicts the exponential probability distribution of the ages of the nearest particle pairs. Particle-resolved numerical simulations with moving particles are performed to verify this prediction. The equation is then used to derive the evolution equation for the particle–fluid–particle (PFP) stress, which is known to be related to hyperbolicity of the two-fluid equations. It is found that the relaxation time of the age probability distribution is also the relaxation time for the PFP stress. Guided by the closure terms in the PFP stress evolution equation, we study kinematics of the nearest particle pairs in the particle-resolved simulations for flows caused by sedimentation of the particles with initially isotropic and homogeneous particle distributions. At the steady states, the particle Reynolds numbers are around 20. Anisotropy and inhomogeneity of particle distributions are seen to develop in these flows. The mean distances to the nearest particles and evolution of the distribution of the Voronoi cell volumes are studied. We also found the PFP stress is closely related to the changes in these inter-particle scale quantities.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Optimal surface-tension isotropy in the Rothman-Keller color-gradient lattice Boltzmann method for multiphase flow

The Rothman-Keller color-gradient (CG) lattice Boltzmann method is a popular method to simulate two-phase flow because of its ability to deal with fluids with large viscosity contrasts and a wide range of interfacial tensions. Here, two fluids are labeled red and blue, and the gradient in the color difference is used to compute the effect of interfacial tension. It is well known that finite-difference errors in the color-gradient calculation lead to anisotropy of interfacial tension and errors such as spurious currents. Here, we investigate the accuracy of the CG calculation for interfaces between fluids with several radii of curvature and find that the standard CG calculations lead to significant inaccuracy. Specifically, we observe significant anisotropy of the color gradient of order 7% for high curvature of an interface such as when a pinchout occurs. We derive a second order accurate color gradient and find that the diagonal nearest neighbors can be weighted differently than in the usual color-gradient calculation such that anisotropy is minimized to a fraction of a percent. The optimal weights that minimize anisotropy for the smallest radius of curvature interface are found to be w = (0.298, 0.284, 0.275) for diagonal nearest neighbors for the cases of the interface smoothing parameter β = (0.5, 0.7, 0.99), somewhat higher than the w = 0.25 value derived by Leclaire et al. [Leclaire, Reggio, and Trepanier, Computers and Fluids 48, 98 (2011)] based on obtaining isotropic errors to second order. We find that use of these optimal w values yields over a factor of 10 decrease in anisotropy and over a factor of 30 decrease in mean anisotropy relative to using the standard w = 1 value. And we find a factor of about 2 decrease in the anisotropic error and up to factor 15 decrease in mean anisotropic error relative to the choice of w = 0.25 for small radius of curvature interfaces. The improved CG calculations will allow the method to be more reliably applied to studies of phenomenology and pore scale processes such as viscous and capillary fingering, and droplet formation where surface-tension isotropy of narrow fingers and small droplets plays a crucial role in correctly capturing phenomenology. We present an example illustrating how different phenomena can be captured using the improved color-gradient method. Namely, we present simulations of a wetting fluid invading a fluid filled pipe where the viscosity ratio of fluids is unity in which droplets form at the transition to fingering using the improved CG calculations that are not captured using the standard CG calculations. We present an explanation of why this is so which relates to anisotropy of the surface tension, which inhibits the pinchouts needed to form droplets.

58 GEOSCIENCES↗
DOE OSTI2016

Uncertainty quantification tools for multiphase gas-solid flow simulations using MFIX

Computational fluid dynamics (CFD) has been widely studied and used in the scientific community and in the industry. Various models were proposed to solve problems in different areas. However, all models deviate from reality. Uncertainty quantification (UQ) process evaluates the overall uncertainties associated with the prediction of quantities of interest. In particular it studies the propagation of input uncertainties to the outputs of the models so that confidence intervals can be provided for the simulation results. In the present work, a non-intrusive quadrature-based uncertainty quantification (QBUQ) approach is proposed. The probability distribution function (PDF) of the system response can be then reconstructed using extended quadrature method of moments (EQMOM) and extended conditional quadrature method of moments (ECQMOM). The report first explains the theory of QBUQ approach, including methods to generate samples for problems with single or multiple uncertain input parameters, low order statistics, and required number of samples. Then methods for univariate PDF reconstruction (EQMOM) and multivariate PDF reconstruction (ECQMOM) are explained. The implementation of QBUQ approach into the open-source CFD code MFIX is discussed next. At last, QBUQ approach is demonstrated in several applications. The method is first applied to two examples: a developing flow in a channel with uncertain viscosity, and an oblique shock problem with uncertain upstream Mach number. The error in the prediction of the moment response is studied as a function of the number of samples, and the accuracy of the moments required to reconstruct the PDF of the system response is discussed. The QBUQ approach is then demonstrated by considering a bubbling fluidized bed as example application. The mean particle size is assumed to be the uncertain input parameter. The system is simulated with a standard two-fluid model with kinetic theory closures for the particulate phase implemented into MFIX. The effect of uncertainty on the disperse-phase volume fraction, on the phase velocities and on the pressure drop inside the fluidized bed are examined, and the reconstructed PDFs are provided for the three quantities studied. Then the approach is applied to a bubbling fluidized bed with two uncertain parameters, particle-particle and particle-wall restitution coefficients. Contour plots of the mean and standard deviation of solid volume fraction, solid phase velocities and gas pressure are provided. The PDFs of the response are reconstructed using EQMOM with appropriate kernel density functions. The simulation results are compared to experimental data provided by the 2013 NETL small-scale challenge problem. Lastly, the proposed procedure is demonstrated by considering a riser of a circulating fluidized bed as an example application. The mean particle size is considered to be the uncertain input parameter. Contour plots of the mean and standard deviation of solid volume fraction, solid phase velocities, and granular temperature are provided. Mean values and confidence intervals of the quantities of interest are compared to the experiment results. The univariate and bivariate PDF reconstructions of the system response are performed using EQMOM and ECQMOM.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Multiphase Flow Calculations in CartaBlanca

This presentation describes CartaBlanca, a multi-material and multi-physics code and its application in traditional fluid dynamics and advanced multi-material interactions.

42 ENGINEERING↗
DOE OSTI2017

ASSESSMENT OF MULTIPHASE CFD WITH ZERO CLOSURE MODEL FOR BOILING WATER REACTOR FUEL ASSEMBLIES

Extension of existing boiling closures to the high void fraction conditions of BWRs is one of challenges for M-CFD applications. This paper presents the assessment of an Eulerian M-CFD model with Zero closure model that has been specifically designed for two-phase flow and heat transfer applications in BWR fuel assemblies. The baseline Zero closure aims at including the simplest and most robust representation of the key physical mechanisms, rather than aiming at absolute accuracy of the predictions. The capabilities of M-CFD with Zero were evaluated against the international OECD/NRC BFBT benchmark data through the comparative analyses of local void fraction, sub-channel void fraction error and exit quality. INL high performance computing system has been employed to demonstrate the scalability of the M-CFD application to nuclear reactor fuel assemblies. The M-CFD simulation results show a good agreement with the experimental data for both quantitative and qualitative results, with the highest uncertainty being observed in the prediction of the local void fraction profile at low void fraction condition. This paper presents a recent view of the best practices for M-CFD of BWR fuel assemblies to provide a useful reference to identify and overcome common challenges in the applications of M-CFD tools. We provide recommendations for geometric modeling, mesh generation and solver configurations to optimize the methodology application.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Development of a Computational Framework for Multiphysics Multiphase Species Tracking using NEAMS Tools

This report implements a high-fidelity multiphysics modeling framework using the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program tools to track isotopic species in Molten Salt Reactors (MSRs), with a specific focus on the 91-depletion chain within the Molten Salt Reactor Experiment (MSRE). The model integrates neutronics, thermal-hydraulics, depletion, and thermochemistry to simulate the production, transport, and phase transitions of isotopes under steady-state and transient conditions. The main findings reveal that isotopes such as bromine-91 largely remain in the liquid phase, while others, including krypton-91and yttrium-91, transition to the gas phase, significantly influencing the reactor’s radiological source term. The study also shows that during transients, like a reactivity insertion transient, rapid void formation and the expansion of the liquid-gas interface led to substantial transfers of dissolved isotopes into the gas phase, altering isotope distribution and largely increasing the source term in the off-gas system. Additionally, the research highlights that short-lived isotopes dominate the initial off-gas response during transients, while longer-lived isotopes determine the equilibrium state, underscoring the necessity of dynamic simulations for accurate species tracking and reactor safety analysis. The developed methodology will be applied in the future to the tracking of a larger number of species and introduce other species tracking mechanisms, such as deposition and plating.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2017

Physical mechanisms for multiphase flow associated with hydrate formation

Many Arctic hydrate reservoirs such as those of the Prudhoe Bay and Kuparuk River area on the Alaska North Slope (ANS) are believed originally to be natural gas accumulations converted to hydrate accumulations after being placed in the gas hydrate stability zone (GHSZ) in response to ancient climate cooling. In this paper, the implications of a previously described mechanistic model for the transport of gaseous and aqueous phases are studied using a transient 1-D transport model during the conversion of a gas reservoir to a hydrate reservoir. The mechanistic model predicts/explains the vertical profile of hydrate saturation in “converted free gas” hydrate reservoirs. The initial gas phase saturation with depth is estimated from the profile of capillary entry pressure, which is estimated from grain size distributions measured in cores. The gas accumulation is assumed to be disconnected from its original source so that methane transport occurs only within it. As the base of the GHSZ descends through the sediment, hydrate forms within the GHSZ. The net volume reduction associated with hydrate formation creates a “sink” which drives flow of gaseous and aqueous phases to the hydrate formation zone. Mechanisms by which this fluid movement could have occurred are analyzed. Flow driven by saturation gradients plays a key role in creating reservoirs of large hydrate saturations, as observed in Mount Elbert stratigraphic test well in the Milne Point Unit of Alaska North Slope (ANS). Viscous-dominated pressure-driven flow of gaseous and aqueous phases cannot explain large hydrate saturations originated from large-saturation gas accumulations. The mode of hydrate formation for a wide range of rate of hydrate formation, the rate of descent of the base of GHSZ, and host sediment characteristics are analyzed and characterized based on dimensionless groups. The proposed transport model is also consistent with field data from hydrate-bearing sand units in Mount Elbert well. Finally, results show that not only the petrophysical properties of the host sediment but also the rate of hydrate formation and the rate of temperature cooling at the surface contribute greatly to the final hydrate saturation profiles.

58 GEOSCIENCES↗
DOE OSTI2017

Micro-PIV measurements of multiphase flow of water and liquid CO 2 in 2-D heterogeneous porous micromodels

We present an experimental study of pore-scale flow dynamics of liquid CO 2 and water in a two-dimensional heterogeneous porous micromodel, inspired by the structure of a reservoir rock, at reservoir-relevant conditions (80 bar, 21°C). The entire process of CO 2 infiltration into a water-saturated micromodel was captured using fluorescence microscopy and the micro-PIV method, which together reveal complex fluid displacement patterns and abrupt changes in velocity. The CO 2 front migrated through the resident water in an intermittent manner, forming dendritic structures, termed fingers, in directions along, normal to, and even opposing the bulk pressure gradient. Such characteristics indicate the dominance of capillary fingering through the micromodel. Velocity burst events, termed Haines jumps, were also captured in the heterogeneous micromodel, during which the local Reynolds number was estimated to be ~21 in the CO 2 phase, exceeding the range of validity of Darcy's law. Furthermore, these drainage events were observed to be cooperative (i.e., across multiple pores simultaneously), with the zone of influence of such events extending beyond tens of pores, confirming, in a quantitative manner, that Haines jumps are nonlocal phenomena. After CO 2 completely breaks through the porous section, shear-induced circulations caused by flowing CO 2 were also observed, in agreement with previous studies using a homogeneous porous micromodel. To our knowledge, this study is the first quantitative measurement that incorporates both reservoir-relevant conditions and rock-inspired heterogeneity, and thus will be useful for pore-scale model development and validation.

42 ENGINEERING↗
DOE OSTI2024

Digital Multiphase Composites via Additive Manufacturing

Abstract Mechanical properties of traditional engineering materials are typically coupled to each other, presenting a challenge to practitioners with multi‐dimensional material property requirements. In this work, continuous, independent control over multiple mechanical properties is demonstrated in composite materials realized using additive manufacturing. For the first time, composites additively manufactured from rigid plastic, soft elastomer, and liquid constituents are experimentally characterized, demonstrating materials which span four orders of magnitude in modulus and two orders of magnitude in toughness. By forming analytical mappings between relative concentrations of constituents at the microscale and resulting macroscale material properties, inverse material design is enabled; the method is showcased by printing artifacts with prescribed toughness and elasticity distributions. The properties of these composites are placed in the context of biological tissues, showing they have promise as mechanically plausible tissue mimics.

Smith, Lawrence T.↗
DOE OSTI2015

Multiphase Nanostructure of a Quinary Metal Oxide Electrocatalyst Reveals a New Direction for OER Electrocatalyst Design

Ce-rich mixed metal oxides comprise a recently discovered class of ­electrocatalysts for the oxygen evolution reaction (OER). In particular, at current densities below 10 mA cm -2 , Ni 0.3 Fe 0.07 Co 0.2 Ce 0.43 O x exhibits ­superior activity compared to the corresponding transition metal oxides, despite the relative inactivity of ceria. To elucidate the enhanced activity and underlying catalytic mechanism, detailed structural characterization of this quinary oxide electrocatalyst is reported in this paper. Transmission electron microscopy imaging of cross-section films as-prepared and after electrochemical testing reveals a stable two-phase nanostructure composed of 3–5 nm diameter crystallites of fluorite CeO 2 intimately mixed with 3–5 nm crystallites of transition metal oxides alloyed in the rock salt NiO structure. Dosing experiments demonstrate that an electron flux greater than ≈1000 e Å -2 s -1 causes the inherently crystalline material to become amorphous. A very low dose rate of 130 e Å -2 s -1 is employed for atomic resolution imaging using inline holography techniques to reveal a nanostructure in which the transition metal oxide nanocrystals form atomically sharp boundaries with the ceria nanocrystals, and these results are corroborated with extensive synchrotron X-ray absorption spectroscopy measurements. Finally, ceria is a well-studied cocatalyst for other heterogeneous and electrochemical reactions, and our discovery introduces biphasic cocatalysis as a design concept for improved OER electrocatalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗