Engineering Papers⌕ Search

Engineering topics

Xu, Zhijie

Publications and source records attributed to Xu, Zhijie.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Advanced CO 2 Capture Solvent Systems for Dynamic Power Generation: Quarterly Research Performance Progress Report, QR4 (Q4FY24)

We developed an integrated Computational Fluid Dynamics (CFD) model to simulate the multi-physics coupled cooling process of mixed gas by cold water within a Direct Contact Cooler (DCC) equipped with a rotating packing bed (RPB). The model captures the interactions between fluid dynamics, heat transfer, mass transport, and phase transitions, while accounting for key operational variables such as RPB rotational speed and the mass flow rates of both liquid and gas. The CFD model has been validated using experimental data, specifically by comparing predicted outflow gas and liquid temperatures to measured results. Our findings demonstrate the significant effects of RPB rotational speed and mass flow rates on cooling performance, providing valuable insights for optimizing DCC efficiency in industrial applications.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

pnnl/EZBattery

Cell performance optimization is important for improving the system efficiency of a redox flow battery. To gain better insights into key controlling factors Of system efficiency, this work first proposed a theoretical model for a unit cell by extending a two-dimensional analytic model to a full cell. The model is then used for cell performance optimization after validating it with experimental and numerical modeling data.

Bao, Jie↗
DOE OSTI2024

Evaluating large scale aqueous organic redox flow battery performance with a hybrid numerical and machine learning framework

Aqueous organic redox flow battery (AORFB) is a promising cost-competitive technology for large-scale energy storage. Among existing work, the dihydroxyphenazine (DHP)-based AORFB has demonstrated high energy density and low-capacity degradation in 10 cm$^2$ cells during lab tests. However, its commercial-scale performance in more complex environments remains unknown, posing a barrier to commercialization. To address this gap, this work presents a comprehensive performance evaluation of a 780 cm$^2$ DHP-based AORFB by combining a physics-based numerical model, machine learning (ML)-based surrogate models, and ML-derived sensitivity quantification. Specifically, we first select 12 key battery parameters that include 10 physicochemical and 2 operation quantities, then select 6 performance metrics that include energy efficiency (EE), discharging capacity, charging energy, and power losses due to concentration, activation, and ohmic over-potentials. With such selection, 12800 combinations of the 12 parameters are subsequently generated using the Latin Hypercube Sampling method. Such combinations, together with 38 pre-defined State of Charge, are then integrated to a validated AORFB model developed in COMSOL to compute the performance metrics. With both input parameters and performance metrics, 60 deep neural network (DNN) surrogate models are then trained to approximate the relationship between the 10 physicochemical quantities and 6 performance metrics at each flow rate and current density. Sensitivity scores are then calculated based on the DNN models. Two additional sensitivity analysis tools, i.e., MARS, and SHAP, are also used to cross-validate the sensitivity scores from the DNN. The results demonstrate that 1) the standard potential ranks first in controlling EE and charging energy, 2) the membrane conductivity is most critical for power loss and EE, and 3) specific area and reaction rate control activation power loss.

25 ENERGY STORAGE↗
DOE OSTI2024

On the statistical theory of self-gravitating collisionless dark matter flow: High order kinematic and dynamic relations

Dark matter, if it exists, accounts for five times as much as ordinary baryonic matter. To better understand the self-gravitating collisionless dark matter flow on different scales, a statistical theory involving kinematic and dynamic relations must be developed for different types of flow, e.g., incompressible, constant divergence, and irrotational flow. This is mathematically challenging because of the intrinsic complexity of dark matter flow and the lack of a self-closed description of flow velocity. Here, this paper extends our previous work on second-order statistics Xu to kinematic relations of any order for any type of flow. Dynamic relations were also developed to relate statistical measures of different orders. The results were validated by N-body simulations. On large scales, we found that (i) third-order velocity correlations can be related to density correlation or pairwise velocity; (ii) the pth-order velocity correlations follow ∝ a (p+2)/2 for odd p and ∝ a p/2 for even p, where a is the scale factor; (iii) the overdensity δ is proportional to density correlation on the same scale, $\langle$δ$\rangle$∝$\langle$δδ'$\rangle$; (iv) velocity dispersion on a given scale r is proportional to the overdensity on the same scale. On small scales, (i) a self-closed velocity evolution is developed by decomposing the velocity into motion in haloes and motion of haloes; (ii) the evolution of vorticity and enstrophy are derived from the evolution of velocity; (iii) dynamic relations are derived to relate second- and third-order correlations; (iv) while the first moment of pairwise velocity follows $\langle$Δu L $\rangle$=-Har (H is the Hubble parameter), the third moment follows $\langle$(Δu L ) 3 $\rangle$ ∝ ε u ar that can be directly compared with simulations and observations, where ε u ≈ 10 -7 m 2 /s 3 is the constant rate for energy cascade; (v) the pth order velocity correlations follow ∝ a (3p-5)/4 for odd p and ∝ a 3p/4 for even p. Finally, the combined kinematic and dynamic relations lead to exponential and one-fourth power-law velocity correlations on large and small scales, respectively.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Important Role of Ion Flux Regulated by Separators in Lithium Metal Batteries

Abstract Polyolefin separators are the most common separators used in rechargeable lithium (Li)‐ion batteries. However, the influence of different polyolefin separators on the performance of Li metal batteries (LMBs) has not been well studied. By performing particle injection simulations on the reconstructed three‐dimensional pores of different polyethylene separators, it is revealed that the pore structure of the separator has a significant impact on the ion flux distribution, the Li deposition behavior, and consequently, the cycle life of LMBs. It is also discovered that the homogeneity factor of Li‐ion toward Li metal electrode is positively correlated to the longevity and reproducibility of LMBs. This work not only emphasizes the importance of the pore structure of polyolefin separators but also provides an economic and effective method to screen favorable separators for LMBs.

25 ENERGY STORAGE↗
DOE OSTI2023

Experimental measurement of the effective contact angle for solvent/packing interactions in a structured packed column for CO 2 capture

The contact angle is a critical factor for determining the effective mass transfer area for carbon dioxide (CO 2 ) capture via the chemical absorption process in a packed column, and thus the overall capture efficiency of the packed column. Many widely used commercial packings involve microscale features (perforation, corrugation, etc.) that may also affect the wetting behavior. This study proposes a systematic method of using a modified Wilhelmy plate to measure the effective contact angle to characterize the solvent and featured packing interaction. In lieu of computational efforts relying on assumptions guided by semi-informed correlation, the proposed method directly measures the effective contact angle as a function of the solvent and packing thermophysical and hydrodynamic properties. The characterization of the effective contact angle is then integrated in the computational fluid dynamics modeling to reduce uncertainty in the prediction of effective mas transfer area. Experiments were conducted for stainless steel coupons using water and aqueous sodium hydroxide (NaOH) solvent for verification of the proposed experimental protocol. The surface tension of the aqueous NaOH solvent was altered using surfactant and antifoam. The effective contact angle increases with the increased value of the surface tension for flat stainless steel sheets. On the other hand, effective contact angles do not vary in Mellapak coupons for aqueous monoethanolamine (MEA) and NaOH solvent. In this case, surface textures and sheet design play a dominant role in the surface tension of solvents. Furthermore, CO 2 loading has a significant effect on the contact angle for Mellapak coupons. As expected, the contact angle decreases with the increasing temperature of CO 2 capture solvents (MEA, EEMPA). The effective contact angle measurement using the Wilhelmy plate method can provide more accurate and efficient solvents for characterizing the solvent-packing surface interactions. Subsequently, the present method enhances the accuracy of the prediction of the effective mass transfer area in carbon capture by solvent absorption.

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

Dark matter halo mass functions and density profiles from mass and energy cascade

Abstract Halo abundance and structure play a central role for modeling structure formation and evolution. Without relying on a spherical or ellipsoidal collapse model, we analytically derive the halo mass function and cuspy halo density (inner slope of −4/3) based on the mass and energy cascade theory in dark matter flow. The hierarchical halo structure formation leads to halo or particle random walk with a position-dependent waiting time $$\tau _g$$ τ g . First, the inverse mass cascade from small to large scales leads to the halo random walk in mass space with $$\tau _g\propto m_h^{-\lambda }$$ τ g ∝ m h - λ , where $$m_h$$ m h is the halo mass and $$\lambda$$ λ is a halo geometry parameter with predicted value of 2/3. The corresponding Fokker-Planck solution for halo random walk in mass space gives rise to the halo mass function with a power-law behavior on small scale and exponential decay on large scale. This can be further improved by considering two different $$\lambda$$ λ for haloes below and above a critical mass scale $$m_h^*$$ m h ∗ , i.e. a double- $$\lambda$$ λ halo mass function. Second, a double- $$\gamma$$ γ density profile can be derived based on the particle random walk in 3D space with a position-dependent waiting time $$\tau _g \propto \Phi (r)^{-1} \propto r^{-\gamma }$$ τ g ∝ Φ ( r ) - 1 ∝ r - γ , where $$\Phi$$ Φ is the gravitational potential and r is the particle distance to halo center. Theory predicts $$\gamma =2/3$$ γ = 2 / 3 that leads to a cuspy density profile with an inner slope of −4/3, consistent with the predicted scaling laws from energy cascade. The Press-Schechter mass function and Einasto density profile are just special cases of proposed models. The small scale permanence can be identified due to the scale-independent rate of mass and energy cascade, where density profiles of different halo masses and redshifts converge to the $$-4/3$$ - 4 / 3 scaling law ( $$\rho _h \propto r^{-4/3}$$ ρ h ∝ r - 4 / 3 ) on small scales. Theory predicts the halo number density scales with halo mass as $$\propto m_h^{-1.9}$$ ∝ m h - 1.9 , while the halo mass density scales as $$\propto m_h^{4/9}$$ ∝ m h 4 / 9 . Results were compared against the Illustris simulations. This new perspective provides a theory for nearly universal halo mass functions and density profiles.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

A hybrid analytical and numerical model for cross-over and performance decay in a unit cell vanadium redox flow battery

Developing an accurate and efficient model for cross-over is critical for improving the long-term performance of redox flow batteries (RFBs). Here, this work presents a hybrid analytical and numerical model that combines a two-dimensional analytical solution to the active species, a one-dimensional analytical model for cross-over mechanisms, and a zero-dimensional numerical model for outlet concentrations of reactants. By comparing with experiment over 41 cycles (ca. 144 h), the model reported a mean voltage difference of 0.0089 V, mean time difference of 14 s per cycle (of 3.5 h), maximum relative difference for capacity and energy of 1.34% and 1.63%, respectively. The predicted mean concentrations for V 2+ , V 3+ , and VO$_{2}^{+}$ in membrane are 65%~77% of measured values from the literature. Upon validation, the model reproduced behaviors in electrolyte imbalance similar to those observed in experiments and numerical models, and revealed the control of cross-over, self-discharge, stoichiometry of side reactions, and Coulombic efficiency on electrolyte imbalance. The model also demonstrates excellent computational efficiency for simulating 41 cycles (around 37,000 points) within 3~4 s. The demonstrated efficiency and accuracy for predicting cross-over and its impacts on voltage, capacity&energy decay, membrane concentrations, and electrolyte imbalance makes it a reliable tool for optimizing RFBs’ long-term performance.

25 ENERGY STORAGE↗
DOE OSTI2023

Packed Column Simulations

CO2 Capture CFD (C3FD) dataset contains STAR-CCM+ simulations of a CO2-capturing solvent flowing across packing structures in a column. This dataset contains 50 2D and 50 3D simulations, each of which spans 500 timesteps. The 2D and 3D simulation meshes have 150,073 and 3,129,815 nodes, respectively. Each node contains information about its physical location, liquid volume fraction, momentum, and node type.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2023

Universal scaling laws and density slopes for dark matter haloes

Smalls scale challenges suggest some missing pieces in our current understanding of dark matter. In this work, a cascade theory for dark matter is proposed to provide extra insights, similar to the cascade phenomenon in hydrodynamic turbulence. The kinetic energy is cascaded in dark matter from small to large scales involves a constant rate $ε_u$ (≈–4.6×10 –7 m 2 / s 3 ). Confirmed by N-body simulations, the energy cascade leads to a two-thirds law for kinetic energy $v^2_r$ on scale $\textit{r}$ such that $v^2_r∝(ε_ur)^{2/3}$. Equivalently, a four-thirds law can be established for mean halo density $ρ_s$ enclosed in the scale radius $r_s$ such that $ρ_s ∝ε^{2/3}_uG^{–1}r^{–4/3}_s$, which was confirmed by galaxy rotation curves. Critical properties of dark matter might be obtained by identifying key constants on relevant scales. First, the largest halo scale $r_l$ can be determined by $–u^3_0/ε_u$, where u0 is the velocity dispersion. Second, the smallest scale rη is dependent on the nature of dark matter. For collisionless dark matter, $r_η∝(–Gℏ/ε_u)^{1/3}$ ≈ 10 –13 m , where ℏ is the Planck constant. An uncertainty principle for momentum and acceleration fluctuations is also postulated. For self-interacting dark matter, $r_η∝ε^2_uG^{–3}(σ/m)^3$, where σ/m is the cross-section of interaction. On halo scale, the energy cascade leads to an asymptotic density slope γ = –4/3 for fully virialized haloes with a vanishing radial flow, which might explain the nearly universal halo density. Based on the continuity equation, halo density is analytically shown to be closely dependent on the radial flow and mass accretion, such that simulated haloes can have different limiting slopes. A modified Einasto density profile is proposed accordingly.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Understanding of Low‐Porosity Sulfur Electrode for High‐Energy Lithium–Sulfur Batteries

Abstract The lithium–sulfur (Li–S) battery is a promising technology for large‐scale energy storage and vehicle electrification due to its high theoretical energy density and low cost. Reducing the sulfur cathode porosity has been identified recently as a viable strategy for improving the cell practical energy density and minimizing pore‐filling electrolytes to extend cell life at lean electrolyte conditions. Direct use of a low‐porosity cathode for Li–S battery results in poor electrode wetting, nonuniform electrode reactions, and thus early cell failure. To understand and mitigate the barriers associated with the use of low‐porosity electrodes, multiscale modeling is performed to predict electrode wetting, electrolyte diffusion, and their impacts on sulfur reactions in Li–S cells by explicitly considering the electrode wettability impacts and electrode morphologies. The study elucidates the critical impact of low tortuosity and large channel pore design for promoting electrode wetting and species diffusion. It is suggested that the secondary particle size should be comparable with the electrode thickness to effectively promote electrolyte wettability and sulfur reactivity. This study provides new insights into the low‐porosity electrode material and designs and is expected to accelerate the development of practical high‐energy Li–S batteries.

25 ENERGY STORAGE↗
DOE OSTI2023

Microstructure-process relationships in monolithic U-10Mo fuel foil single-pass rolling: A parametric simulation study

In this work, a previously validated coupling of Kinetic Monte Carlo (KMC) Potts Model and finite element method (FEM) simulations was implemented to investigate the effects of microstructural features in as-cast and homogenized monolithic U-10Mo foils on the emergent microstructure after rolling and reheating. Parameters that could potentially affect recrystallization behavior of the rolled U-10Mo foils were considered: grain size distribution, uranium carbide (UC) size distribution, UC volume fraction, spatial distribution of UC, and rolling reduction magnitude. Grain structure and the magnitude of rolling reduction have the strongest influence on recrystallization kinetics and the fabricated grain size distribution. The UC distribution had only a weak effect on the recrystallization kinetics and final microstructures. While particle-stimulated nucleation (PSN) occurred in simulation more frequently as grain size increased, its incidence did not appear to considerably affect the recrystallization kinetics or grain size distribution.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Characterization of Electrochemical Behavior for Aqueous Organic Redox Flow Batteries

Use of aqueous redox flow batteries with organic redox-active materials holds great promise for large-scale and sustainable energy storage. The development of low-cost, highly efficient aqueous redox flow batteries lies in a comprehensive understanding of the electrochemical behaviors of redox-active compounds. In this work, an alkaline redox battery with organic dihydroxyphenazine sulfonate (DHPS) anolyte and ferro-/ferricyanide (Fe(CN) 6 ) catholyte is investigated as a typical example of aqueous redox flow batteries using organic redox-active materials. The electrochemical kinetics of DHPS and Fe(CN) 6 are separately characterized using the symmetrical cell design. The resistance components are calculated directly from the experimental measurement. The key kinetic parameters are extracted and compared for DHPS and Fe(CN) 6 electrolytes. The extracted parameters are validated with symmetrical and full flow cell simulations at different operating conditions. Key parameters and internal loss are also compared with all-vanadium redox flow batteries, representing current state of the art. In addition, our extracted key parameters from a symmetrical flow cell are compared with the measured key parameters by cyclic voltammetry, a widely deployed electroanalytical technique. The cell performance prediction of DHPS anolyte on a 780 cm 2 interdigitated cell is made and found the power density is peaked at 475 mW cm -2 at our measurement condition.

25 ENERGY STORAGE↗
DOE OSTI2022

A three-dimensional pore-scale model for redox flow battery electrode design analysis

A three-dimensional (3-D) pore-scale model has been developed to construct multiscale fibrous electrodes for redox flow batteries (RFB). New designs, such as biporous electrodes modify electrode structures by creating secondary pores on single carbon fiber to reduce internal battery resistance. Existing pore-scale models only resolve electrodes to the single carbon fiber scale and cannot incorporate recent multiscale electrode designs into numerical models. Our new model aims to bridge this gap and provide a tool to rapidly screen new electrode configurations. Two multiscale electrodes, laser-perforated and biporous electrodes, were investigated at varying operation conditions with the proposed 3-D pore-scale models. The laser-perforated electrode exhibits a reduced pressure drop, but follows the same permeability correlation compared to the corresponding pristine electrode. For the biporous electrode, the added specific surface area and faster reaction kinetics from the secondary pores are the most influential factors leading to improved battery efficiency. However, operating the biporous electrode in limiting current density conditions should be avoided due to the decreased mass transfer efficiency and a more significant voltage loss. Finally, we believe that our 3-D pore-scale model can accelerate the flow battery electrode design process and provide new insights into electrode geometry optimizations.

25 ENERGY STORAGE↗
DOE OSTI2022

Hydrodynamics of countercurrent flow in an additive-manufactured column with triply periodic minimal surfaces for carbon dioxide capture

Computational investigation of the countercurrent flows in an additively manufactured (AM) column packed with triply periodic minimal surface (TPMS) is challenging. Multiphase flow simulations using the volume of fluid (VOF) method can explain the column’s hydrodynamic characteristics, including interfacial area, liquid holdup, pressure drop, flow regime, etc. The interfacial area is a key factor dictating the mass transfer via the absorption process, and thereby the overall efficiency of the column. The impacts of a solvent's physical properties, liquid and gas loads, and contact angle on the hydrodynamics are extensively explored. The solvents include traditional solvent (monoethanolamine [MEA]) used in industry as well as a water-lean solvent (EEMPA) as prospective solvents for carbon capture. As expected, interfacial area and liquid holdup increase with increasing liquid loads (q L ). The EEMPA exhibits higher values for the interfacial area and liquid holdup compared to those for the MEA. The dry and wet pressure drops per unit length in the column at different gas loads (F G ) are calculated and further compared with the empirical correlations for traditional packings. At low gas loads, i.e., in the preloading regime, the predicted values of the wet and dry pressure drops match well with the corresponding values calculated from available correlations. The gas load has a marginal impact on the interfacial area at its lower value, but the higher value of F G (>2.18 $\sqrt{Pa}$) causes flooding in the column. The effect of varying contact angles on the interfacial area and liquid holdup in the preloading regime is also studied. Both interfacial area and liquid holdup decrease with the increase in the value of the contact angle. Altogether, the effect of the contact angle on interfacial area is more pronounced than other parameters influencing the interfacial area.

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

Differentiation of static and dynamic interfacial area in the structured packed column

Effective mass-transfer area plays a key role for post-combustion carbon capture. The static area, an inactive part relative to mass transfer, needs to be differentiated in the interfacial area. In the present study, CFD method is employed to investigate static and dynamic interfacial areas in Mellapak 500.Y structured packing. Here, a wide range of physical properties and loading conditions are included to understand their effects on both areas. The result shows that the relationship between interfacial area and interfacial velocity at the local scale can be described by gamma distribution with good accuracy. The influences of various physical properties derived from simulation are compared with existing experiment-based correlations for both static and interfacial areas. The result shows the dominant influence of viscosity on the static area, and the most importance of the influence of contact angle on the interfacial area. This study also highlights the importance of viscosity, which is usually ignored.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗