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 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
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 OSTI2015

A numerical comparison of precipitating turbulent flows between large‐eddy simulation and one‐dimensional turbulence

This study presents the results of computational fluid dynamics simulations of a multiphase, reacting, turbulent mixing layer in an idealized geometry. The purpose is to compare large‐eddy simulation (LES) to one‐dimensional turbulence (ODT) and examine the trends of the flow under differing mixing conditions. Aqueous streams are mixed together to precipitate polymorphs of calcium carbonate. The polymorphs of calcium carbonate are tracked numerically using population balance equations (PBE). Each PBE contains all of the relevant physical models to track the particle evolution including nucleation, growth, and aggregation. A simple subgrid mixing model that is convenient for use with PBEs was implemented in the LES code. The higher spatial resolution achievable with ODT allowed an investigation on the effect of resolution on the mixing‐model error. The Reynolds number of the flow is varied and is shown to cause a decrease in average particle sizes with higher mixing rates. © 2015 American Institute of Chemical Engineers AIChE J , 61: 3185–3197, 2015

Abboud, Alex W.↗
DOE OSTI2018

Advanced trust region optimization strategies for glass box / black box models

We present an improved trust region filter (TRF) method for optimization of combined glass box/black box systems. Glass box systems refer to models that are easily expressed in an algebraic modeling language, providing cheap and accurate derivative information. By contrast, black box systems may be computationally expensive and derivatives are unavailable. The TRF method, as first introduced in our previous work (Eason and Biegler, AIChE J . 2016; 62:3124–3136), is able to handle hybrid systems containing both glass and black box components, which can frequently arise in chemical engineering, for example, when a multiphase reactor model is included in a flow sheet optimization problem. We discuss several recent modifications in the algorithm such as the sampling region, which maintains the algorithm's global convergence properties without requiring the trust region to shrink to zero in the limit. To benchmark the development of this optimization method, a test set of problems is generated based on modified problems from the CUTEr and COPS sets. The modified algorithm demonstrates improved performance using the test problem set. Finally, the algorithm is implemented within the Pyomo environment and demonstrated on a rigorous process optimization case study for carbon capture. © 2018 American Institute of Chemical Engineers AIChE J , 64: 3934–3943, 2018

Eason, John P.↗
DOE OSTI2019

Experimental investigation on multiscale hydrodynamics in a novel gas–Liquid–Solid three phase jet-Loop reactor

Multiphase flow hydrodynamics in a novel gas–liquid–solid jet-loop reactor (JLR) were experimentally investigated at the macroscales and mesoscales. The chord length distribution was measured by an optical fiber probe and transformed for bubble size distribution through the maximum entropy method. The impacts of key operating conditions (superficial gas and liquid velocity, solid loading) on hydrodynamics at different axial and radial locations were comprehensively investigated. JLR was found to have good solid suspension ability owing to the internal circulation of bubbles and liquid flow. The gas holdup, axial liquid velocity, and bubble velocity increase with gas velocity, while liquid velocity has little influence on them. Compared with the gas–liquid JLRs, solids decrease the gas holdup and liquid circulation, reduces the bubble velocity and delays the flow development due to the enhanced interaction between bubbles and particles (Stokes number >1). This work also provides a benchmark data for computational fluid dynamics (CFD) model validation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2015

Impact of drilling fluid rheology and wellbore pressure on rock cuttings removal performance: numerical investigation

Drilling fluids exhibit non-Newtonian behavior which is generally described in the literature by two rheological models: Power Law and Herschel-Bulkley models. We have used Eulerian- Eulerian multiphase flow model to evaluate the drilling fluid flow patterns and filter cake formation patterns for both models for ultra-deep drilling conditions. Simulated drilling fluid flow patterns clearly show that the Herschel - Bulkley model exhibits fewer vortices around the drilling pipe, indicating better rock cuttings removal performance. In addition to mud rheology, we have analyzed the pressure drop variation with the filter cake formation on the wellbore during ultradeep drilling process. We have shown that filter cake thickness at the wall linearly increases with differential pressure between the well and formation. This implies that the differential pressure needs to be optimized. Very thin filter cake may lead to excessive fluid seepage into the permeable subterranean vertical zone resulting in a loss of fluid circulation and hence in a loss of rock cuttings removal. Yet, very thick filter cake reduces the effective diameter of the hole which may lead to a poor rock cuttings removal performance and stuck pipe. The mathematical model presented here needs to be experimentally validated when filter cake thickness data in deep drilling conditions become available. The model could be useful to provide drilling engineers insight on the filter cake formation, drilling fluid flow pattern, and give guide-line to optimize the filter cake as field experimental data are scarce.

02 PETROLEUM↗
DOE OSTI2018

A Solvent-Free Synthesis of Lignin-Derived Renewable Carbon with Tunable Porosity for Supercapacitor Electrodes

Synthesis of multiphase materials from lignin, a biorefinery coproduct, offers limited success owing to the inherent difficulty in controlling dispersion of these renewable hyperbranched macromolecules in the product or its intermediates. Effective use of the chemically reactive functionalities in lignin, however, enables tuning morphologies of the materials. Here in this study, we bind lignin oligomers with a rubbery macromolecule followed by thermal crosslinking to form a carbon precursor with phase contrasted morphology at submicron scale. The solvent-free mixing is conducted in a high-shear melt mixer. With this, the carbon precursor is further modified with potassium hydroxide for a single-step carbonization to yield activated carbon with tunable pore structure. A typical precursor with 90 % lignin yields porous carbon with 2120 m 2 g -1 surface area and supercapacitor with 215 F g -1 capacitance. Lastly, the results show a simple route towards manufacturing carbon-based energy-storage materials, eliminating the need for conventional template synthesis.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Modeling CO 2 flow in support of a shallow subsurface controlled leakage field test

Here, controlled release of carbon dioxide (CO 2 ) into the soil and atmosphere is performed to test detection and monitoring tools, for which several field laboratories were established by a number of institutions worldwide. Numerical simulations of CO 2 behavior in the shallow subsurface region are other forms of validation and verification of the leakage pathways and destinations. These reports aim to improve monitoring and verification of CO 2 in case of unexpected leakages for public assurance. In this work, we introduce the results of a numerical modeling study conducted to simulate the injection of CO 2 as carried out during a field test in Viamão, southern Brazil, where 20 kg day –1 of CO 2 was pumped for 30 days through a vertical well 3 m below ground in an altered granitic soil. Multiphase flow simulations were performed with the TOUGH2/EOS7CA software for unsaturated porous media, using field data and injection parameters, including sensitivity tests to permeability direction, diffusivity, and boundary conditions. Results with increased horizontal permeabilities are in better agreement with the field observations. In this condition, mass balance calculations indicate approximately 90% of injected CO2 (20 kg day –1 during 30 days) remains in the soil after 180 days from injection start, consistent with the measured flow through the soil–atmosphere interface.

03 NATURAL GAS↗
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↗
DOE OSTI2019

Juxtatumoral perinephric fat analysis in clear cell renal cell carcinoma

Purpose: The purpose of the study was to evaluate the feasibility of using contrast-enhanced computed tomography (CECT)-based texture analysis (CTTA) metrics to differentiate between juxtatumoral perinephric fat (JPF) surrounding low-grade (ISUP 1–2) versus high-grade (ISUP 3–4) clear cell renal cell carcinoma (ccRCC). Methods: In this IRB-approved study, we retrospectively queried the surgical database between June 2009 and April 2016 and identified 83 patients with pathologically confirmed ccRCC (low grade: n = 54, mean age = 61.5 years, 18F/35M; high grade n = 30, mean age = 61.7 years, 8F/22M) who also had pre-operative multiphase CT acquisitions. CT images were transferred to a 3D workstation, and nephrographic phase JPF regions were manually segmented. Using an in-house developed Matlab program, a CTTA panel comprising of texture metrics extracted using six different methods, histogram, 2D- and 3D-Gray-level co-occurrence matrix (GLCM) and Gray-level difference matrix (GLDM), and 2D-Fast Fourier Transform (FFT) analyses, was applied to the segmented images to assess JPF textural heterogeneity in low- versus high-grade ccRCC. Univariate analysis and receiver-operator characteristics (ROC) analysis were used to assess interclass differences in texture metrics and their prediction accuracy, respectively. Results: All methods except GLCM consistently revealed increased heterogeneity in the JPF surrounding high- versus low-grade ccRCC. FFT showed increased complexity index (p < 0.01). Histogram analysis showed increased kurtosis and positive skewness in (p < 0.03), and GLDM analysis showed decreased measure of correlation coefficient (MCC) (p < 0.04). Several of the GLCM metrics showed statistically significant (p < 0.04) textural differences between the two groups, but with no consistent trend. ROC analysis showed that MCC in GLCM analysis had an area under the curve of 0.75. Conclusions: Our study suggests that CTTA of ccRCC shows statistically significant textural differences in JPF surrounding high- versus low-grade ccRCC.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2019

Visibility of Hypovascularized Liver Tumors during Intra-Arterial Therapy Using Split-Bolus Single-Phase Cone Beam CT

Purpose: To validate a split-bolus contrast injection protocol for single-phase CBCT in terms of detectability of hypovascular liver tumors compared to digital subtraction angiography (DSA). Materials and Methods: In this retrospective, single-center study, 20 consecutive patients with in total 77 hypovascularized tumors referred for intra-arterial therapy received a split-bolus single-phase CBCT. Two readers rated the visibility of the target tumors scheduled for embolization in CBCT and DSA compared to the pre-interventional multiphasic CT or MRI used as reference on a 3-point scoring system (1 = optimal, 3 = not visible) and catheter-associated artifacts (1 = none, 3 = extended). SNR, CNR and contrast values were derived from 37 target tumors in CBCT and MRI. Statistical analysis included the kappa test to determine interrater reliability, the Friedman’s test for the inter-modality comparison evaluating tumor visibility in DSA and CBCT as well as for quantitative assessment. Post hoc analysis included the Wilcoxon signed-rank test. p values < 0.05 were considered significant. Results: Ninety percentage of target tumors were rated as visible in CBCT and 37.5% in DSA (p < 0.001). 70.1% of pre-interventionally detected hypovascularized tumors were depicted with CBCT and 31.2% by DSA (p < 0.001). 7.8% of known tumors were outside the FOV. Quantitative assessment showed higher image contrasts in CBCT (1.91 ± 7.01) compared to hepatobiliary-phase MRI (0.29 ± 0.14, p = 0.003) and to portal-venous (p.v.) MRI (0.31 ± 0.13, p < 0.001), but higher CNR for MRI (1.18 ± 0.80; 13.92 ± 15.82; 13.79 ± 6.65). Conclusion: In conclusion, the split-bolus single-phase CBCT detects significantly more hypovascularized liver tumors compared to DSA performed through the proper hepatic artery with high image contrasts. Level of Evidence: Level III, diagnostic study.

62 RADIOLOGY AND NUCLEAR MEDICINE↗