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.
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.
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.
Methods and apparatus for additive manufacturing utilizing multifunctional composite materials, and articles made therefrom
A method of depositing a multiphase material. The method includes providing a Continuous Multifunctional Composite (CMC) phase containing at least one continuous element in a polymeric matrix, passing the CMC phase through a feeding system containing a cutting system, producing a predetermined length of the CMC phase, providing a flow a molten polymer such that the molten polymer and the CMC phase are merged into a continuous co-extrusion nozzle so as to produce a co-extruded multiphase material, and depositing the co-extruded multiphase material onto a surface. An apparatus for depositing a multiphase material. The apparatus contains a co-extrusion nozzle, a means to introduce a CMC phase and a molten polymer into the co-extrusion nozzle, such that the molten polymer and the CMC phase are co-extruded and deposited on a surface. An article containing a CMC phase containing continuous elements embedded in a polymer resin forming a multiphase structure.
Methods and apparatus for additive manufacturing utilizing multifunctional composite materials, and articles made therefrom
A method of depositing a multiphase material. The method includes providing a Continuous Multifunctional Composite (CMC) phase containing at least one continuous element in a polymeric matrix, passing the CMC phase through a feeding system containing a cutting system, producing a predetermined length of the CMC phase, providing a flow a molten polymer such that the molten polymer and the CMC phase are merged into a continuous co-extrusion nozzle so as to produce a co-extruded multiphase material, and depositing the co-extruded multiphase material onto a surface. An apparatus for depositing a multiphase material. The apparatus contains a co-extrusion nozzle, a means to introduce a CMC phase and a molten polymer into the co-extrusion nozzle, such that the molten polymer and the CMC phase are co-extruded and deposited on a surface. An article containing a CMC phase containing continuous elements embedded in a polymer resin forming a multiphase structure.
Methods and apparatus for additive manufacturing utilizing multifunctional composite materials, and articles made therefrom
A method of depositing a multiphase material. The method includes providing a Continuous Multifunctional Composite (CMC) phase containing at least one continuous element in a polymeric matrix, passing the CMC phase through a feeding system containing a cutting system, producing a predetermined length of the CMC phase, providing a flow a molten polymer such that the molten polymer and the CMC phase are merged into a continuous co-extrusion nozzle so as to produce a co-extruded multiphase material, and depositing the co-extruded multiphase material onto a surface. An apparatus for depositing a multiphase material. The apparatus contains a co-extrusion nozzle, a means to introduce a CMC phase and a molten polymer into the co-extrusion nozzle, such that the molten polymer and the CMC phase are co-extruded and deposited on a surface. An article containing a CMC phase containing continuous elements embedded in a polymer resin forming a multiphase structure.
Realistic fuel spray modeling for gasoline direct injection engine applications
Fuel spray modeling plays a critical role during modern gasoline direct injection (GDI) engine development due to fuel injection’s dominant impact on engine performance and emissions as well as the complex physical processes involved. In engineering three-dimensional (3D) computational fluid dynamics (CFD) simulations, the liquid-phase fuel atomization, evaporation, and mixing are usually modeled with the discrete droplet model (DDM) adopting a Lagrangian approach for multiphase CFD simulations. To this end, general practices heavily depend on the reduced order characterization of the injector nozzle flow. However, such simplified injector modeling may lead to insufficient representations of the complex spray dynamics. To tackle this problem, this study proposes a novel workflow to numerically evaluate GDI sub-cooled and flash-boiling sprays under engine-relevant conditions using a side-mounted GDI injector together with real gasoline fuel properties. The workflow introduces a one-way coupling (OWC) method leveraging high-fidelity nozzle flow simulations to provide realistic boundary conditions to the Lagrangian injector model. The proposed workflow was first verified in a constant volume chamber (CVC) environment and then implemented in a practical GDI engine setup to study spray morphology, fuel-air mixing, and wall-wetting propensity. In addition, detailed comparison was performed between the OWC method and the conventional rate of injection (ROI) routine. In conclusion, quantitative analysis of spray characteristics was conducted to highlight possible source of discrepancies of the conventional ROI method.
A Kinetic Model of the Long-Term Corrosion of Glass-Ceramic Materials
Multiphase waste forms show promise for increased waste loading and for the ability to dispose of contaminated solid and particulate waste through direct densification. However, achieving predictive capability for long-term durability of multiphase waste forms, and thus assessing their possible deployment, requires expanding the current, limited knowledge base. Here, we describe the development of a corrosion model of a two-phase waste form consisting of crystals of known volume fraction embedded in a glass matrix. This model accounts for the dissolution of both the crystalline and glass phases as well as the hydration of the glass phase through an ion exchange reaction. Because of the large difference in solubility between the two phases, the reactive surface of the crystalline phase is a function of the extent of dissolution of the glass phase in this model. Model parameterization was performed using corrosion data, such as from single-pass flow-through tests, for the individual phases. The parameterized corrosion model was evaluated against static dissolution test data for a glass-ceramic multiphase waste form. This evaluation demonstrated the model’s ability to reproduce the time-dependent release of key tracers of glass and crystalline phase dissolution. Hence, the development of a kinetic model provides a pathway for long-term durability predictions and thus the use of multiphase waste forms in nuclear cleanup missions.
Utilization of event shape in search of the chiral magnetic effect in heavy-ion collisions
The search for the chiral magnetic effect (CME) has been a subject of great interest in the field of high-energy heavy-ion collision physics, and various observables have been proposed to probe the CME. Experimental observables are often contaminated with background contributions arising from collective motions (specifically elliptic flow) of the collision system. We present a method study of event-shape engineering (ESE) that projects the CME-sensitive γ 112 correlator and its variations (γ 132 and γ 123 ) to a class of events with minimal flow. We discuss the realization of the zero-flow mode, the sensitivity on the CME signal, and the corresponding statistical significance for Au + Au, Ru + Ru, and Zr + Zr collisions at √s NN = 200 GeV with a multiphase transport model, as well as a new event generator, event-by-event anomalous-viscous fluid dynamics.
Special Issue: To Commemorate the Legacy of Malcolm J. Andrews—Scientist, Educator, and Leader
This special issue celebrates the scientific legacy of Dr. Malcolm J. Andrews, former Editor of the ASME Journal of Fluids Engineering (2010–2018). Malcolm devoted his career to the development and applications of computational, experimental, and analytical tools to investigate and solve diverse, complex fluid dynamics, and heat transfer problems. Malcolm’s most notable research accomplishments include: (i) contributions to innovative experiments, simulations, and theory to understand buoyancydriven turbulent mixing due to Rayleigh–Taylor instability; (ii) development of the multiphase particle-in-cell method for understanding fragmentation in spray jets, coal-water slurry, and other dense particulate flows; (iii) mixing in nano- and biomaterials in applications ranging from protein glues to carbon nanotube composites; (iv) determining dispersion in proteins to improve therapeutic drug delivery via micro-encapsulation; (v) innovative methods for storing information using nanopinned dipole magnetic vortices in superconducting materials; and (vi) modeling and design modifications for shell-and-tube heat exchangers, as well as heat transfer in rotating systems such as turbine blades. The results of his research have been published in a long series of trend-setting journal and conference papers, book chapters, and patents.
Influence of Clay Wettability Alteration on Relative Permeability
Abstract Understanding the wettability of porous materials is important to model fluid flow in the subsurface. One of the critical factors that influences wetting in real reservoirs is the composition of geologic materials. The wetting properties for clay minerals can have a particularly strong impact on flow and transport. In this work, we analyze the chemical composition of a Mt. Simon sandstone core to resolve the microscopic structure of clay regions and assess how alterations to the local wetting properties influence multiphase transport based on core flooding experiments and relative permeability simulations. We show that whichever fluid has greater affinity toward clay minerals will tend to accumulate within these high surface area regions, leading to dramatic shifts in the relative permeability. This work establishes the essential importance of the mineral composition and associated wetting properties in the modeling of flow and transport in reservoir‐scale systems.
Accelerating Multiphase Simulations With Denoising Diffusion Model Driven Initializations
This study introduces a hybrid fluid simulation approach that integrates generative diffusion models with physics‐based simulations, aiming at reducing the computational costs of flow simulations while still honoring all the physical properties of interest. Pore‐scale simulations enhance our understanding of applications such as assessing hydrogen and storage efficiency in underground reservoirs. Nevertheless, they are computationally expensive and the presence of non‐unique solutions can require multiple simulations within a single geometry. To overcome the computational cost hurdle, we propose a method that couples generative diffusion models and physics‐based simulations. While training the data‐driven model, we simultaneously generate initial conditions and perform physics‐based simulations using these. This integrated approach enables us to receive real‐time feedback on a single compute node equipped with both CPUs and GPUs. By efficiently managing these processes within a single compute node, we can continuously monitor performance and halt training once the model meets the specified criteria. To test our model, we generate realizations in a real Berea sandstone fracture which shows that our technique is up to 4.4 times faster than commonly used flow simulation initializations.
Velocity structure functions in multiphase turbulence: interpreting kinematics of Hα filaments in cool-core clusters
ABSTRACT The central regions of cool-core galaxy clusters harbour multiphase gas, with gas temperatures ranging from $10$ to $10^7\, \mathrm{K}$. Feedback from active galactic nuclei jets prevents the gas from undergoing a catastrophic cooling flow. However, the exact mechanism of this feedback energy input is unknown, mainly due to the lack of velocity measurements of the hot-phase gas. However, recent observations have measured the velocity structure functions (VSFs) of the cooler molecular (${\sim} 10\, \mathrm{K}$) and Hα filaments (${\sim} 10^4\, \mathrm{K}$) and used them to indirectly estimate the motions of the hot phase. In the first part of this study, we conduct high-resolution (3843–15363 resolution elements) simulations of homogeneous isotropic subsonic turbulence, without radiative cooling. We analyse the second-order velocity structure functions (VSF2) in these simulations and study the effects of varying spatial resolution, the introduction of magnetic fields, and the effect of projection along the line of sight (LOS) on it. In the second part of the study, we analyse high-resolution (7683 resolution elements) idealized simulations of multiphase turbulence in the intracluster medium from the companion study Mohapatra et al. We compare the VSF2 for both the hot ($T\sim 10^7\, \mathrm{K}$) and cold ($T\sim 10^4\, \mathrm{K}$) phases and find that their amplitude depends on the density contrast between the phases. They have similar scaling with separation, but introducing magnetic fields steepens the VSF2 of only the cold phase. We also find that projection along the LOS steepens the VSF2 for the hot phase and mostly flattens it for the cold phase.
Uncertainty Quantification for Multiphase Computational Fluid Dynamics Closure Relations with a Physics-Informed Bayesian Approach
Multiphase Computational Fluid Dynamics (MCFD) based on the two-fluid model is considered a promising tool to model complex two-phase flow systems. MCFD simulation can predict local flow features without resolving interfacial information. As a result, the MCFD solver relies on closure relations to describe the interaction between the two phases. Those empirical or semi-mechanistic closure relations constitute a major source of uncertainty for MCFD predictions. In this paper, we leverage a physics-informed uncertainty quantification (UQ) approach to inversely quantify the closure relations’ model form uncertainty in a physically consistent manner. This proposed approach considers the model form uncertainty terms as stochastic fields that are additive to the closure relation outputs. Combining dimensionality reduction and Gaussian processes, the posterior distribution of the stochastic fields can be effectively quantified within the Bayesian framework with the support of experimental measurements. As this UQ approach is fully integrated into the MCFD solving process, the physical constraints of the system can be naturally preserved in the UQ results. Here, in a case study of adiabatic bubbly flow, we demonstrate that this UQ approach can quantify the model form uncertainty of the MCFD interfacial force closure relations, thus effectively improving the simulation results with relatively sparse data support.
Multiphase Processing of the Water-Soluble and Insoluble Phases of Biomass Burning Organic Aerosol
Biomass burning is one of the most significant sources of organic aerosol in the atmosphere. Biomass burning organic aerosol (BBOA) has been observed to undergo liquid– liquid phase separation (LLPS) to give core–shell morphology with the hydrophobic phase encapsulating the hydrophilic phase, potentially impacting the evolution of light-absorbing components, i.e., brown carbon (BrC), through multiphase processes. Here, we demonstrate how multiphase processing differs between the watersoluble (i.e., hydrophilic) and insoluble (i.e., hydrophobic) phases of BBOA in terms of reactive uptake of ozone in a coated-wall flow tube. Effects of relative humidity (RH) and ultraviolet (UV) irradiation were investigated. Experimental timeseries were used to inform simulations using multilayer kinetic modeling. Among non-irradiated thin films, the uptake coefficient was greatest for the water-soluble phase at 75% RH (3 × 10 –5 , corresponding to a diffusion coefficient of BrC, D BrC , of 3 × 10 –9 cm 2 s –1 ) and least for the same phase at 0% RH (1 × 10 –5 , corresponding to D BrC of 1 × 10 –10 cm 2 s –1 ). The uptake coefficient for the water-insoluble phase fell between these two (about 1.5 × 10 –5 ), regardless of RH, and the corresponding D BrC increased only slightly (8 × 10 –10 cm 2 s –1 at 0% RH to 9 × 10 –10 cm 2 s –1 at 75% RH). The uptake coefficients of both phases at 0% RH decreased significantly after UV irradiation, consistent with a transition from viscous liquid to solid and supported by qualitative microscopy observations. Modeling multiphase ozone oxidation of primary BrC components in the atmosphere demonstrated, first, that LLPS may extend the lifetime of water-soluble BBOA encapsulated by water-insoluble species by a factor of 1.5 at moderate to high RH and, also, that UV irradiation may extend the lifetime of both phases by more than a factor of 2.5.
Advanced Computational Modeling of High-Level Waste Vitrification at the Hanford Site
The U.S. Department of Energy (DOE) has selected vitrification for stabilizing legacy tank waste at the Hanford site, where radioactive waste from plutonium production was historically stored in underground tanks. This waste will be separated into low-activity waste (LAW) and high-level waste (HLW) fractions and processed at the Waste Treatment and Immobilization Plant (WTP). At WTP, glass melters are used for the vitrification of radioactive tank waste, transforming it into a stable borosilicate glass form for safe long-term storage. The melter vessel is constructed from highly durable and heat-resistant materials, where the vitrification process occurs. The main regions that are modeled are the melt pool, plenum, cold cap, riser/discharge chamber, and surrounding structure with insulation layers. Forced convection induced by air bubblers at the base of the melter ensure uniform temperature distribution and provide heat to the cold cap layer. The cold cap is a region of reacting batch feed that floats on top of the molten glass and is where the batch-to-glass reactions occur. Joule heating provided by electrodes mounted along the vertical walls of the melter and immersed directly in the glass, generates the necessary heat for the net endothermic conversion processes that occur in the cold cap. The high temperatures, radioactivity, and opaque nature of the glass prevent direct observation inside the melters. Therefore, computational models are essential for providing insight into factors that affect melter throughput. Thermocouples in the plenum provide operators with plenum temperature measurements. Operational adjustments include bubbling rate, voltage supplied to the electrodes, feed adjustments, and glass removal rate. Different computational fluid dynamics (CFD) models have been developed, each serving a specific purpose. There are CFD models of different scale melters, as well as models that capture the two-phase flow interfaces of rising bubbles in the molten glass or models with a simplified molten glass region so that the surrounding structure and plenum can be feasibly incorporated. Pilot-scale melter models have been developed to serve as validation of the methods employed in the simulation of the full-scale WTP melters. Models incorporating resolved bubbling are used to develop momentum source terms to implement into a single phase, multi-region, steady-state flow model that is being validated by measured process parameters such as glass production rate, voltage, input power, plenum temperatures, etc. The resolved bubbling model uses the multiphase volume of fluid approach to model the system with a high-resolution interface capturing scheme to maintain sharp interfaces between the molten glass and the air phase. The suite of CFD models is continually being improved to incorporate more realistic physics and achieve faster turnaround time. For example, an incremental controller is implemented to automatically adjust electrode voltage within the simulation to a molten glass set point temperature of 1150°C. Newer models feature improved meshes to ensure conformal meshes between regions and eliminate unnecessary mesh refinement in areas that are not of interest (such as boundary layers in offgas ports). Instead of explicitly modeling the structural, refractory, and insulation layers of the melter, a thermal resistance approach is used with published correlations used for boundary conditions. The development of robust and efficient CFD models will be instrumental in enabling the WTP to successfully fulfill its mission of safely stabilizing legacy nuclear waste.
Flow-plane decorrelations in heavy-ion collisions with multiple-plane cumulants
The azimuthal correlations between local flow planes at different (pseudo)rapidities (η) may reveal important details of the initial nuclear matter density distributions in heavy-ion collisions. Extensive experimental measurements of a factorization ratio (r 2 ) and its derivative (F 2 ) have shown evidence of the longitudinal flow-plane decorrelation. However, nonflow effects also affect this observable and prevent a quantitative understanding of the phenomenon. In this paper, to distinguish decorrelation and nonflow effects, we propose a new cumulant observable, T 2 , which largely suppresses nonflow. Here, the technique sensitivity to different initial-state scenarios and nonflow effects are tested with a simple Monte Carlo model, and in the end, the method is applied to events simulated by a multiphase transport model (AMPT) for Au+Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. We also emphasize that a distinct decorrelation signal requires not only the right sign of an observable but also its proper dependence on the η window of the reference flow plane to be consistent with the pertinent decorrelation picture.
Fractionation of Lignin for Selective Shape Memory Effects at Elevated Temperatures
We report a facile approach to control the shape memory effects and thermomechanical characteristics of a lignin-based multiphase polymer. Solvent fractionation of a syringylpropane-rich technical organosolv lignin resulted in selective lignin structures having excellent thermal stability coupled with high stiffness and melt-flow resistance. The fractionated lignins were reacted with rubber in melt-phase to form partially networked elastomer enabling selective programmability of the material shape either at 70 °C, a temperature that is high enough for rubbery matrix materials, or at an extremely high temperature, 150 °C. Utilizing appropriate functionalities in fractionated lignins, tunable shape fixity with high strain and stress recovery, particularly high-stress tolerance were maintained. Detailed studies of lignin structures and chemistries were correlated to molecular rigidity, morphology, and stress relaxation, as well as shape memory effects of the materials. The fractionation of lignin enabled enrichment of specific lignin properties for efficient shape memory effects that broaden the materials’ application window. Electron microscopy, melt-rheology, dynamic mechanical analysis and ultra-small angle neutron scattering were conducted to establish morphology of acrylonitrile butadiene rubber (NBR)-lignin elastomers from solvent fractionated lignins.