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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Numerical Analysis of Water Harvesting from Saturated Vapor with the Electrospray Water Droplet Injection in The Solar-Driven Desalination Applications

Solar-driven seawater desalination systems provide freshwater through environmentally friendly and carbon-neutral processes. Air humidification and dehumidification desalination (HDD) systems extract water directly from moist air. This paper focuses on a new solar-driven seawater desalination system that shares some common traits with HDD systems but in which the air is eliminated from the process. Seawater is vaporized by solar thermal radiation in high-performance solar panels. Water vapor flows to a new electro-condensation chamber and passes through a series of electrosprays, which inject small nuclei of freshwater droplets into the vapor clouds. Here, due to dielectrophoresis and electrohydrodynamic flows, the vapor molecules are captured by the charged droplets, and vapor is condensed at the droplet surface. The water vapor condensation was modeled by using the multiphase flow numerical model. The CFD model, which was also experimentally validated, used SprayFoam solver modules available in the Open FOAM® CFD open-source software. In preliminary results, 802 grams/(hr-m2) (grams per hour per unit area) of water was harvested directly from moist air at 24C and 90% relative humidity. The simulation results indicated that freshwater productivity could increase up to 4,000 grams/(hr-m2) for fully saturated air conditions. This represents about 5 times the water productivity rate of current conventional HDD systems today. In addition, energy consumption could be reduced by 80 percent compared to conventional HDD systems.

Electrospray, electro-condensation, solar, droplet↗
DOE OSTI2026

Reduced-Order CFD Modeling to Support Waste Loading Optimization in Hanford WTP Vitrification

The U.S. DOE Hanford Site stores over 56 million gallons of radioactive liquid tank waste that must be treated and immobilized for long-term disposal The Waste Treatment and Immobilization Plant (WTP) will vitrify this waste by feeding it into Joule-heated melters, where it is incorporated into a stable borosilicate glass Computational fluid dynamics (CFD) simulations of glass melters can provide insight into the maximum achievable waste loading under varying melter operating conditions Fully resolved VOF multiphase simulations were used as the reference model to capture bubble-driven convection in the melter, including bubble formation, rise behavior, and induced glass melt circulation Effective bubble column diameter and rise velocity were extracted from the resolved simulations, compared with empirical correlations, and refit across relevant viscosity and gas flow rate conditions Explicit gas–liquid interface tracking was replaced with a single-phase momentum source term model, enabling faster steady-state CFD simulations while preserving the dominant hydrodynamic effects of bubbling New empirical correlations were developed for effective bubble column diameter and bubble rise velocity by fitting resolved simulation data across expected melter viscosity and gas flow rate ranges, providing improved inputs for the momentum source term model compared with existing literature correlations The momentum source term model reduced fluid-domain mesh size by 89% and achieved an 8.4× computational speedup relative to resolved bubbling simulations The validated momentum source term approach enables prediction of process-relevant heat transfer behavior in the integrated melter model, including heat transfer from the molten glass to the cold cap, plenum, refractory walls, and surrounding structural regions under varying melter operating conditions

12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗
DOE OSTI2022

Fundamental Effects of Al and Ta on Microstructure and Phase Transformations in the Al–Cr–Mo–Ta–Ti Refractory Complex Concentrated Alloy System

The effect of aluminum and tantalum concentrations on a refractory metal complex concentrated alloy is reported, particularly with respect to their effect on microstructure and phase composition of the alloy in cast and annealed form. Alloys with an equiatomic composition, (AlCrMoTaTi), an aluminum-lean composition (Al 0.75 CrMoTaTi), and a tantalum-lean composition (AlCrMoTa 0.75 Ti) are produced via arc melting. The alloys exhibit multiphase structures, confirmed by X-ray diffraction, microstructural characterization, and thermal analysis. The minor off-equiatomic adjustments of aluminum and tantalum in this alloy system did not drastically alter the prevalence of the Cr–Ta-based Laves phase. Correlations between thermodynamic predictions and observed phase transformations via thermal analysis are improved upon refinement of calculations removing impractical intermediate phases. As a result, experimental findings provide information for the refinement of thermodynamic modeling and deliver additional insight into the optimization of alloy compositions within this five-component system.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Shape Persistent, Highly Conductive Ionogels from Ionic Liquids Reinforced with Cellulose Nanocrystal Network

Abstract Shape‐persistent, conductive ionogels where both mechanical strength and ionic conductivity are enhanced are developed using multiphase materials composed of cellulose nanocrystals and hyperbranched polymeric ionic liquids (PILs) as a mechanically strong supporting network matrix for ionic liquids with an interrupted ion‐conducting pathway. The integration of needlelike nanocrystals and PIL promotes the formation of multiple hydrogen bonding and electrostatic ionic interaction capacitance, resulting in the formation of interconnected networks capable of confining a high amount of ionic liquid (≈95 wt%) without losing its self‐sustained shape. The resulting nanoporous and robust ionogels possess outstanding mechanical strength with a high compressive elastic modulus (≈5.6 MPa), comparable to that of tough, rubbery materials. Surprisingly, these rigid materials preserve the high ionic conductivity of original ionic liquids (≈7.8 mS cm −1 ), which are distributed within and supported by the nanocrystal network‐like rigid frame. On the one hand, such stable materials possess superior ionic conductivities in comparison to traditional solid electrolytes; on the other hand, the high compression resistance and shape‐persistence allow for easy handling in comparison to traditional fluidic electrolytes. The synergistic enhancement in ion transport and solid‐like mechanical properties afforded by these ionogel materials make them intriguing candidates for sustainable electrodeless energy storage and harvesting matrices.

Lee, Hansol↗
DOE OSTI2025

Freestanding BaTiO 3 ‐Au Vertically Aligned Nanocomposite toward Flexible Multi‐Sensing Platform

Abstract Flexible and wearable sensors show enormous potential for personalized healthcare devices by real‐time monitoring of an individual's health. Typically, a single functional material is selected for one sensor to sense a particular physical signal while multiple materials will be selected for multi‐mode sensing. Vertically aligned nanocomposites (VANs) have recently demonstrated various material combinations and novel coupled multifunctionalities that are hard to achieve in any single‐phase material alone, including multiphase multiferroics, magneto‐optic coupling, and strong magnetic and optical anisotropy. Integrating these novel VANs into wearable sensors shows enormous potential in multi‐mode sensing owing to their multifunctional nature. In this work, the transfer of VANs onto polydimethylsiloxane as a novel flexible chemical and pressure sensor is demonstrated. For this demonstration, the classical BaTiO 3 ‐Au VAN with combined plasmonic and piezoelectric properties is used to demonstrate a multi‐sensing mechanism. A thin water‐soluble buffer of Sr 3 Al 2 O 6 serves as a buffer layer for the epitaxial growth and transfer process. The electrical output based on the piezoelectric responses and identifying 4‐mercaptobenzoic acid by surface‐enhanced Raman spectroscopy reveal great potential for free‐standing VANs in a wearable multifunctional sensing platform.

Tsai, Benson Kunhung [School of Materials Engineer↗
DOE OSTI2025

Improving Ionic Conformality Across Polymer Electrolyte|Electrode Interfaces

Maintaining uniform ionic transport at electrode|electrolyte interfaces, i.e., ionic conformality, remains challenging in polymer electrolyte (PE)-based solid-state batteries. Morphological conformality does not necessarily imply ionic conformality. In PEs, which typically consist of a mechanically supporting component and distinct ionically conductive components, the rearrangement or depletion of mobile ion-conductive domains at interfaces can disrupt ionic transport pathways. Such localized ionic depletion contributes to interfacial instability and capacity degradation in high-voltage lithium-metal batteries. Herein, an electrolyte design approach aimed at minimizing interfacial heterogeneities is demonstrated through compositional adjustments, characterized by spatially resolved structural and chemical X-ray techniques and NMR diffusometry to elucidate ion transport dynamics. This approach improves ionic conformality at electrode interfaces, enhancing cycling stability in Li||LiNi 0.8 Co 0.1 Mn 0.1 O 2 (NMC811) coin and pouch cells cycled at high voltages. These results contribute to understanding interfacial behaviors in multiphase PEs and inform strategies for improving stability across solid-state battery interfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Electronic Excitation‐Driven β‐Ga 2 O 3 Metastability Transformation and Self‐Organization Mechanism: β→κ/γ/δ Phases

Here, irradiation-driven multiphase self-organization presents emergent opportunities for the customization of nanoscale engineering properties, dynamically tuning strain-field distributions and interfacial electronic structures. Responding to intense electronic excitation-induced energy deposition, the dominant phase transformations, with varying Gibbs free energy $Δ,G^o_f$ are confirmed as β → κ → γ → δ that are located in specific microregions for Gallium (III) oxide (Ga 2 O 3 ), as follows: (i) Surface-localized interstitial accumulation under compressive stress triggers β → δ via semi-coherent interface formation. (ii) Tensile stress within latent tracks drives vacancy-mediated oxygen layer truncation (4/12 periodicity along ⟨0001⟩), stabilizing coherent 4H (ABCB) κ and 3C (ABC) β (ABC) interfaces through strain-compensated octahedral distortion. (iii) Screw dislocation-mediated lattice relaxation induces β → γ via cation disordering (Ga 3+ occupancy at β-interstitial sites), forming metastable spinel γ with mixed occupancy across 16d/8a Wyckoff sites. Irradiation-driven β-Ga 2 O 3 →κ/γ/δ transitions, as mechanistically revealed via inelastic thermal spike (i-TS) calculations and molecular dynamics simulations, induce defect-mediated nonlinear photoresponse, critical for optoelectronic engineering.

Han, Xinqing [Shandong Univ., Jinan (China)]↗
DOE OSTI2024

Phase Transformation Driven by Oxygen Vacancy Redistribution as the Mechanism of Ferroelectric Hf 0.5 Zr 0.5 O 2 Fatigue

As a promising candidate for nonvolatile memory devices, the hafnia-based ferroelectric system has recently been a hot research topic. Although significant progress has been made over the past decade, the endurance problem is still an obstacle to its final application. In perovskite-based ferroelectrics, such as the well-studied Pb[Zr x Ti 1-x ]O 3 (PZT) family, polarization fatigue has been discussed within the framework of the interaction of charged defects (such as oxygen vacancies) with the moving domains during the switching process, particularly at the electrode-ferroelectric interface. Armed with this background, a hypothesis is set out to test that a similar mechanism can be in play with the hafnia-based ferroelectrics. The conducting perovskite La-Sr-Mn-O is used as the contact electrode to create La 0.67 Sr 0.33 MnO 3 / Hf 0.5 Zr 0.5 O 2 (HZO)/ La 0.67 Sr 0.33 MnO 3 capacitor structures deposited on SrTiO 3 -Si substrates. Nanoscale X-ray diffraction is performed on single capacitors, and a structural phase transition from polar o-phase toward non-polar m-phase is demonstrated during the bipolar switching process. The energy landscape of multiphase HZO has been calculated at varying oxygen vacancy concentrations. Based on both theoretical and experimental results, it is found that a polar to non-polar phase transformation caused by oxygen vacancy redistribution during electric cycling is a likely explanation for fatigue in HZO.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Highly Active and Durable Air Electrodes for Reversible Protonic Ceramic Electrochemical Cells Enabled by an Efficient Bifunctional Catalyst

The commercialization of reversible protonic ceramic electrochemical cells is hindered by the lack of highly active and durable air electrodes exposed to high concentration of steam under operating conditions. We report findings that dramatically enhance the electrocatalytic activity and stability of a conventional (La 0.6 Sr 0.4 ) 0.95 Co 0.2 Fe 0.8 O 3-δ (LSCF) air electrode by a multiphase catalyst coating composed of a conformal Pr 1-x Ba x CoO 3-δ thin film and exsolved BaCoO 3-δ nanoparticles, are reported. At 600 °C, the catalyst coating decreases the polarization resistance of the LSCF air electrode by a factor of 25 (from 1.09 to 0.043 Ω cm 2 ) in air and the degradation rate by two orders of magnitude (from 1.0 × 10 -2 to 1.8 × 10 -4 Ω cm2 h -1 in humidified air with 30 vol% H2O). Further, a single cell with the catalyst-coated LSCF air electrode at 600 °C demonstrates a high peak power density of 1.04 W cm -2 in the fuel cell mode and a high current density of 1.82 A cm -2 at 1.3 V in the electrolysis mode. The significantly enhanced performance of the LSCF air electrode is attributed mainly to the high rate of surface oxygen exchange, fast surface proton diffusion, and the rapid H 2 O and O 2 dissociation on the catalysts.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Manipulating Na/TM Ratio‐Driven Structural Heterogeneity of O3‐NaNi 1/3 Fe 1/3 Mn 1/3 O 2 Cathode for High‐Voltage Sodium‐Ion Batteries

The stability of O3-type NaNi 1/3 Fe 1/3 Mn 1/3 O 2 under high-voltage cycling is dictated by how synthesis encodes lattice strain and redox heterogeneity. Here, in this study, the role of Na:TM stoichiometry is systematically resolved by tuning the NaOH:precursor ratio during solid-state synthesis. The stoichiometric condition (Na:TM = 1.00) yields minimized microstrain, enabling uniform O3–P3 phase evolution and homogeneous multi-metal redox with preserved octahedral symmetry. In contrast, Na-excess compositions inherit disordered intermediates and heterogeneous distortion fields that trigger abrupt multiphase transitions and promote localized charge redistribution. In situ XRD captures the divergence in phase-transition pathways, TXM resolves particle-level redox heterogeneity, and XANES corroborates a stronger and more reversible Fe redox contribution at stoichiometry, shifting to diminished Fe participation and spatially inhomogeneous redox at higher Na content. These results establish Na:TM stoichiometry as a critical synthesis parameter controlling both structural coherence and redox stability. Electrochemically, the stoichiometric composition exhibits smooth voltage profiles with minimal polarization growth and retains nearly 80% of its initial capacity after 100 cycles even at an extended 4.2 V cutoff, whereas Na-excess compositions show significantly reduced initial coulombic efficiency and rapid voltage fade. Precise stoichiometric tuning provides a scalable route to defect-suppressed O3 frameworks, enabling structurally resilient, high-voltage sodium-layered cathodes.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Open source implementation of glued sphere discrete element method and nonspherical biomass fast pyrolysis simulation

Here, in this research, a glued-sphere Discrete Element Method (DEM) method was implemented in the open-source, Computational Fluid Dynamics (CFD) software MFiX. The implementation was verified using a cylinder-wall collision and then validated by simulating the packing and fluidization of non-spherical particles. The validated code was applied to simulate fast pyrolysis of non-spherical biomass particles in a bubbling fluidized bed. The glued sphere occupancy ratio was proposed to quantify the quality of shape resolution using glued sphere. Shape resolution showed significant influence on the packing height in the simulation of particle packing and an occupancy ratio of 80 % was recommended. Its influence is minor in fully fluidized bed but can be 8 times higher in fixed fluidized bed. Three tested heat transfer models predicted similar yields of elongated biomass fast pyrolysis. The solver developed in this research can be used to simulate other multiphase reacting flows involving non-spherical particles.

42 ENGINEERING↗
DOE OSTI2022

The influence of random packed column parameters on the liquid holdup and interfacial area

Abstract Carbon dioxide capture via solvent absorption in packed columns has emerged as a potential technology to mitigate coal‐fired power plant CO 2 emissions. Parameters, including packing types, solvent properties, and operating conditions, could potentially affect the packed column CO 2 capture efficiency. To understand the importance of those parameters and help packed column optimization, a design of experiments (DoEs) method was proposed to generate input parameter matrix. Combined with multiphase computational fluid dynamics (CFD), the random packed column parameter influence on the liquid holdup and interfacial area can be efficiently investigated. Surrogate‐based sensitivity analysis shows that the solvent flow rate and contact angle are key factors dictating liquid holdup and interfacial area. On the other hand, solvent viscosity has a marginal impact on the interfacial area. The sensitivity scores were calculated for each input parameter to guide the selection of dimensionless numbers for the liquid holdup and interfacial area correlation development.

42 ENGINEERING↗
DOE OSTI2021

Risk-based area of review estimation in overpressured reservoirs to support injection well storage facility permit requirements for CO 2 storage projects

This paper by the Energy & Environmental Research Center presents a workflow and modeling approach for delineating a risk-based area of review (AOR) to support a U.S. Environmental Protection Agency (EPA) Class VI permit for a carbon dioxide (CO 2 ) storage project. The approach combines semianalytical solutions for estimating formation fluid leakage through a hypothetical leaky wellbore with the results of numerical reservoir simulations to define the AOR. The modeling utilizes 1) semianalytical solutions from the peer-reviewed literature for formation fluid leakage through abandoned wellbores by Raven (1990) and Avci (1994), 2) a FORTRAN model compiled and described in Cihan et al. (2011, 2012) called ASLMA (Analytical Solution for Leakage in Multilayered Aquifers), and 3) a computational framework for estimating a risk-based AOR first proposed by Oldenburg et al. (2014, 2016). Therefore, the approach builds upon well-established research and underlying hydrogeological principles that have been upheld for nearly three decades. Moreover, the ASLMA model has been broadly applied to an array of storage projects. The work presented herein extends these earlier works using a custom wrapper written in the software environment, R (R Core Team, 2020), which was developed to perform multiple runs of the ASLMA model using given ranges for one or more input parameters. In addition, the current work simulates the pressure buildup within the storage reservoir in response to CO 2 injection using a compositional simulator to better accommodate the temporospatial evolution of pressure buildup within the storage reservoir that is more accurately modeled using a heterogeneous geologic model and a compositional simulator that accounts for the multiphase interactions. The workflow is demonstrated using a case study for a 180,000-metric-ton-per-year storage project located in the PCOR (Plains CO 2 Reduction) Partnership region. For the storage project evaluated here, under the scenario where the leaky wellbore is open to a saline aquifer (thief zone) between the overlying seal (cap rock) and the underground sources of drinking water (USDW), the risk-based AOR essentially collapses to the areal extent of the CO 2 plume in the storage reservoir because the pressure buildup in the storage reservoir beyond the CO 2 plume is insufficient to drive formation fluids up a hypothetical leaky wellbore into the USDW. However, even under the conservative assumption that the leaky wellbore is not open to a thief zone, beyond the areal extent of the CO 2 plume, the incremental leakage is less than 400 m 3 over 20 years, which represents ~0.0001% or less of the total volume of water contained within the USDW rock volume. As discussed in the text, the threshold criterion for defining the risk-based AOR is site-specific and should be informed by the results of the sensitivity analysis and available site characterization data. The approach outlined in this paper is designed to be protective of USDWs and, therefore, comply with the Safe Drinking Water Act requirements and provisions for the U.S. EPA Class VI Underground Injection Control (UIC) Program (Class VI Rule) and North Dakota Administrative Code Chapter 43-05-01.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

ReactionMechanismSimulator.jl: A modern approach to chemical kinetic mechanism simulation and analysis

Abstract We present ReactionMechanismSimulator.jl (RMS), a modern differentiable software for the simulation and analysis of chemical kinetic mechanisms, including multiphase systems. RMS has already been applied to problems in combustion, pyrolysis, polymers, pharmaceuticals, catalysis, and electrocatalysis. RMS is written in Julia, making it easy to develop and allowing it to take advantage of Julia's extensive numerical computing ecosystem. In addition to its extensive library of optimized analytic Jacobians, RMS can generate and use Jacobians computed using automatic differentiation and symbolically generated analytic Jacobians. RMS is demonstrated to be faster than Cantera and Chemkin in several benchmarks. RMS also implements an extensive set of features for analyzing chemical mechanisms, including a library of easy‐to‐call plotting functions, molecular structure resolved flux diagram generation, crash analysis, traditional sensitivity analysis, transitory sensitivity analysis, and an automatic mechanism analysis toolkit. RMS implements efficient adjoint and parallel forward sensitivity analyses. We also demonstrate the ease of adding new features to RMS.

Johnson, Matthew S.↗
DOE OSTI2024

Sequential formulation of all‐way coupled finite strain thermoporomechanics for largely deformable gas hydrate deposits

We develop a numerically stable sequential formulation of thermoporomechanics for largely deformable gas hydrate deposits, extended from the fixed stress split of infinitesimal transformation. Constitutive equations are based on the total Lagrangian approach for both flow and geomechanics, including dynamic full tensor permeability and thermal conductivity updated from the deformation gradient. For space discretization, we take the cell-centered finite volume and node-based finite element method for flow and geomechanics, respectively. Then, we propose a sequential implicit method for all-way coupled thermoporomechanics, where the nonisothermal multiphase flow problem of gas hydrates is solved implicitly first and then the geomechanics problem is solved implicitly at the next step. During solution of the flow problem, we fix the rate of first Pioal total stress for numerical stability as well as apply porosity correction and entropy correction to account for geomechanical effects. We test numerical examples where flow and geomechanics parameters are based on deep oceanic gas hydrate deposits. When applying depressurization, even though the results between the infinitesimal transformation and finite strain geomechanics are similar in the early stages due to small deformation, we find differences between them in the late times as deformation becomes large. Accordingly, permeability and thermal conductivity tensors become nonisotropic full tensors although they are initially isotropic. Furthermore, we identify numerical stability of the developed sequential method from the test cases that exhibit the highly complex coupled gas hydrate systems with large deformation. Thus, the proposed sequential formulation can be applied in largely deformable gas hydrate systems.

42 ENGINEERING↗
DOE OSTI2026

Extrusion‐Spheronization of Energetic Materials

The prevailing method to produce plastic‐bonded explosive (PBX) molding powder, or “prills”, is a complex, multiphase, and bespoke process that was developed by the high explosives (HEs) manufacturing industry several decades ago. This work demonstrates the utility of a simpler, widely‐used mechanical process—extrusion‐spheronization—to produce PBX prills. We begin by detailing precautions taken to enable safe remote operation of extrusion‐spheronization equipment with HE. We then perform a study investigating the effect of lacquer solvent composition on the particle packing, pressed density, and compressive strength properties of a 95 wt.% TATB/5 wt.% polymer binder formulation akin to PBX 9502. It was found that increased composition of low vapor pressure solvents caused prolonged retention of the solvent, resulting in tackier materials that would agglomerate and form larger prills. The larger prills also led to lower poured density and tapped density of HE prills and compressive strength of pressed PBX articles. The samples prepared with a 75% propyl acetate/25% butyl acetate lacquer solvent composition exhibited the highest compressive strength. However, it is believed that the prill packing and compressive strength properties are primarily driven by the prill size rather than the chemical composition of the lacquer itself. Extrusion‐spheronization remains a promising method to reliably and repeatably produce HE prills that is less sensitive to feedstock or process variation than traditional methods.

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

Discovering Reactant Supply Pathways at Electrode/PEM Reaction Interfaces Via a Tailored Interface-Visible Characterization Cell

In situ and micro-scale visualization of electrochemical reactions and multiphase transports on the interface of porous transport electrode (PTE) materials and solid polymer electrolyte (SPE) has been one of the greatest challenges for electrochemical energy conversion devices, such as proton exchange membrane electrolyzer cells (PEMECs), CO 2 reduction electrolyzers, PEM fuel cells, etc. Here, an interface-visible characterization cell (IV-CC) is developed to in situ visualize micro-scaled and rapid electrochemical reactions and transports in PTE/SPE interfaces. Taking the PEMEC of a green hydrogen generator as a study case, the unanticipated local gas blockage, micro water droplets, and their evolution processes are successfully visualized on PTE/PEM interfaces in a practical PEMEC device, indicating the existence of unconventional reactant supply pathways in PEMs. Further comprehensive results reveal that PEM water supplies to reaction interfaces are significantly impacted with current densities. Here these results provide critical insights about the reaction interface optimization and mass transport enhancement in various electrochemical energy conversion devices.

30 DIRECT ENERGY CONVERSION↗