Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parametric analysis”

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 OSTI2026

Topological Control of Triply Periodic Minimal Surfaces for Thermal Design and Advanced Manufacturing: A Gyroid Case Study

Recently, there has been a heightened interest in using triply periodic minimal surfaces (TPMSs) in the design of compact process engineering components. The benefits of high surface area per unit volume, modular form, and inherent periodicity provide a holistic self-supporting network and flow-conducive features. Applications of importance include thermal power management, biomimetic scaffolds and structures, and feasibility of advanced manufacturing. This study presents a novel approach to the manipulation of the characteristic Schwarz-G, or gyroid TPMS, for thermal design in the context of advanced manufacturing. The study presents relationships between design parameters and resulting surface area as a target response using the characteristic equation of a gyroid. Through parametric control, the characteristic equation is manipulated to produce a 20-fold increase in achievable area over a baseline design characteristic of 25.4 mm through controlled combinations of design parameters. A second relationship is presented as a function of the maximum area achieved and manipulated design parameters. Through the analysis, the study presents a framework to identify and maximize the achievable area of TPMSs for advanced manufacturing and thermal management applications.

gyroid↗
DOE OSTI2025

An analysis of parameter compression and Full-Modeling techniques with Velocileptors for DESI 2024 and beyond

In anticipation of forthcoming data releases of current and future spectroscopic surveys, we present the validation tests and analysis of systematic effects within velocileptors modeling pipeline when fitting mock data from the AbacusSummit N-body simulations. We compare the constraints obtained from parameter compression methods to the direct fitting (Full-Modeling) approaches of modeling the galaxy power spectra, and show that the ShapeFit extension to the traditional template method is consistent with the Full-Modeling method within the standard ΛCDM parameter space. We show the dependence on scale cuts when fitting the different redshift bins using the ShapeFit and Full-Modeling methods. We test the ability to jointly fit data from multiple redshift bins as well as joint analysis of the pre-reconstruction power spectrum with the post-reconstruction BAO correlation function signal. We further demonstrate the behavior of the model when opening up the parameter space beyond ΛCDM and also when combining likelihoods with external datasets, namely the Planck CMB priors. Finally, we describe different parametrization options for the galaxy bias, counterterm, and stochastic parameters, and employ the halo model in order to physically motivate suitable priors that are necessary to ensure the stability of the perturbation theory.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Engineering Layer For System Analysis

ELSA offers various utility classes and methods to streamline the definition of regions, materials, and geometries in nuclear simulations. Key features include generating OpenMC regions, managing material properties, and providing convenient abstractions for complex geometrical and physical configurations. Additionally, ELSA supports the creation of submodels, enabling users to build modular and reusable components for their simulations. The codebase also includes robust extrusion and revolution capabilities, facilitating the efficient creation of 3D parametric geometries from 2D profiles through linear and rotational transformations.

Ferney, Paul [Idaho National Laboratory (INL), Ida↗
DOE OSTI2026

Multiprobe constraints on early and late time dark energy

We perform a multiprobe analysis combining cosmic microwave background (CMB) data from Planck and the Atacama Cosmology Telescope (ACT), ACT CMB lensing, and large-scale structure (LSS) measurements from the Dark Energy Spectroscopic Instrument (DESI), including DESI Legacy Imaging Survey (LS) galaxies and baryon acoustic oscillations (BAOs). We present the first 5×2 pt analysis of ACT DR6 lensing, DESI LS, and Planck Integrated Sachs-Wolfe (ISW). Within Λ⁢CDM, this yields 𝑆 8 =𝜎 8 ⁢(Ω 𝑚 /0.3) 0.5 =0.819 ±0.016, in good agreement with primary CMB inferences and provides a sound-horizon-free Hubble constant constraint of 𝐻 0 =70.0±4.4 km s −1 Mpc −1 . Then, combining with CMB primary and BAO, we reconfirm a CMB–BAO discrepancy in the Ω 𝑚 –$\frac{𝐷𝑣}{𝑟𝑑}$ plane, which is heightened when combining BAO with the 5 ×2 pt data vector. We explore two dark-energy extensions that may reconcile this: an early-time modification, early dark energy (EDE), and late-time dynamical dark energy (DDE) parametrized by 𝑤 0 ⁢𝑤 𝑎 . For CMB primary +BAO +5 ×2 pt, we find a 3.3⁢𝜎 preference for DDE over Λ⁢CDM, while EDE is modestly favored at 2.3⁢𝜎. The models address different shortcomings of Λ⁢CDM: DDE relaxes the neutrino mass bound (𝑀 𝜈 <0.17 eV vs <0.050 eV under Λ⁢CDM), making it compatible with neutrino oscillation measurements, while EDE raises the Hubble constant to 𝐻 0 =70.5±1.2 km s −1 Mpc −1 , easing the discrepancy with SH0ES. However, neither model resolves both issues simultaneously. Our analysis indicates that both DDE and EDE remain viable extensions of Λ⁢CDM within current uncertainties and demonstrates the capacity of combined probes to place increasingly stringent constraints on cosmological parameters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Subtleties in the trainability of quantum machine learning models

A new paradigm for data science has emerged, with quantum data, quantum models, and quantum computational devices. This field, called quantum machine learning (QML), aims to achieve a speedup over traditional machine learning for data analysis. However, its success usually hinges on efficiently training the parameters in quantum neural networks, and the field of QML is still lacking theoretical scaling results for their trainability. Some trainability results have been proven for a closely related field called variational quantum algorithms (VQAs). While both fields involve training a parametrized quantum circuit, there are crucial differences that make the results for one setting not readily applicable to the other. In this work, we bridge the two frameworks and show that gradient scaling results for VQAs can also be applied to study the gradient scaling of QML models. Our results indicate that features deemed detrimental for VQA trainability can also lead to issues such as barren plateaus in QML. Consequently, our work has implications for several QML proposals in the literature. In addition, we provide theoretical and numerical evidence that QML models exhibit further trainability issues not present in VQAs, arising from the use of a training dataset. We refer to these as dataset-induced barren plateaus. These results are most relevant when dealing with classical data, as here the choice of embedding scheme (i.e., the map between classical data and quantum states) can greatly affect the gradient scaling.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Advancing the STS Neutron Moderator Design with an Automated Optimization Workflow and Unstructured Mesh Modeling

With the Second Target Station approaching its final design phase, a detailed neutronics evaluation of its critical components is necessary. Optimizing the dimensions of the two cold-source moderators that are at the heart of this facility presents a multi-objective optimization problem for which an accurate geometric description is crucial. We have applied a fully automated optimization workflow in which a detailed unstructured mesh geometry is automatically generated with Attila4MC, starting from a parametrized CREO geometry followed by preprocessing with SpaceClaim. With this geometry, a MCNP run is performed to calculate the brightness metrics, which are subsequently provided to the optimization algorithm in DAKOTA that provides new parameters and drives the optimization loop until convergence. In this paper, we show the results of the analysis that are used for the final design of the cylindrical and tube moderator. The optimization simulations provide a refinement to and confirmation of the conclusions of the previous design iteration. Additional to the optimization, a sensitivity study is performed to study the effect of minor geometry changes, which is important for the final engineering design. In conclusion, with these studies, we demonstrate that the automated workflow and high-fidelity unstructured mesh modeling are efficient tools for a thorough design evaluation.

DAKOTA↗
DOE OSTI2024

Coupled-Channel Analysis of the χ c 1 ( 3872 ) Line Shape with BESIII Data

We perform a study of the χ c 1 ( 3872 ) line shape using the data samples of e + e − → γ χ c 1 ( 3872 ) , χ c 1 ( 3872 ) → D 0 D ¯ 0 π 0 , and π + π − J / ψ collected with the BESIII detector. The effects of the coupled channels and the off-shell D * 0 are included in the parametrization of the line shape. The line shape mass parameter is obtained to be M X = ( 3871.63 ± 0.1 3 − 0.05 + 0.06 ) MeV . Two poles are found on the first and second Riemann sheets corresponding to the D * 0 D ¯ 0 branch cut. The pole location on the first sheet is much closer to the D * 0 D ¯ 0 threshold than the other, and is determined to be 7.04 ± 0.1 5 − 0.08 + 0.07 MeV above the D 0 D ¯ 0 π 0 threshold with an imaginary part − 0.19 ± 0.0 8 − 0.19 + 0.14 MeV . Published by the American Physical Society 2024

Physics↗
DOE OSTI2025

Equilibrium Core Model for Micro Pebble Bed Reactors Using OpenMC

Estimating the equilibrium state for pebble bed reactors (PBRs) presents complex challenges as it requires simultaneous consideration of changes in the pebbles’ movement as well as their fuel compositions. Whereas traditional approaches use multigroup diffusion codes for neutronics calculations of PBRs’ equilibrium state, the double-heterogeneity of PBRs complicates neutron cross-section generation. Continuous-energy Monte Carlo (MC) methods are better suited for detailed PBR analysis because of their natural handling of double-heterogeneity, but they demand substantially more computational resources. Here, this study introduces a novel method for efficiently estimating the equilibrium state in small and micro PBRs with reduced computational cost. The method is anticipated to accelerate the processes of core design and performing parametric studies for utilizing advanced fuel and structural materials. The HTR-10 reactor design was used for validating the method’s predictions and evaluating its computational efficiency. When compared to reference calculation values from the literature, criticality (k-effective) was predicted to be approximately within the margin of error of the MC transport calculation, average core power density (in megawatts per cubic meter) was predicted within 2.5% relative error, and maximum thermal flux (10 13 n/cm 2 .s −1 ) was predicted within 1.8% relative error. The calculated inventory of fission products and fuel composition in the equilibrium core were within 15% and 16.6%, respectively, when compared to reported values from the literature. The difference is attributed to variance in the considered values of the core temperature, which was found to significantly affect the depletion analyses.

Equilibrium core↗
DOE OSTI2025

Polynomial Chaos Surrogate Construction for Random Fields with Parametric Uncertainty

Engineering and applied science rely on computational experiments to rigorously study physical systems. The mathematical models used to probe these systems are highly complex, and sampling-intensive studies often require prohibitively many simulations for acceptable accuracy. Surrogate models provide a means of circumventing the high computational expense of sampling such complex models. In particular, polynomial chaos expansions (PCEs) have been successfully used for uncertainty quantification studies of deterministic models where the dominant source of uncertainty is parametric. We discuss an extension to conventional PCE surrogate modeling to enable surrogate construction for stochastic computational models that have intrinsic noise in addition to parametric uncertainty. We develop a PCE surrogate on a joint space of intrinsic and parametric uncertainty, enabled by Rosenblatt transformations, which are evaluated via kernel density estimation of the associated conditional cumulative distributions. Furthermore, we extend the construction to random field data via the Karhunen–Loève expansion. We then take advantage of closed-form solutions for computing PCE Sobol indices to perform a global sensitivity analysis of the model which quantifies the intrinsic noise contribution to the overall model output variance. Additionally, the resulting joint PCE is generative in the sense that it allows generating random realizations at any input parameter setting that are statistically approximately equivalent to realizations from the underlying stochastic model. The method is demonstrated on a chemical catalysis example model and a synthetic example controlled by a parameter that enables a switch from unimodal to bimodal response distributions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Magic zeroes and hidden symmetries

Selection rules arising from accidental or broken symmetries may be sufficiently obscure that their agency is hidden, leading to the appearance of “magic zeroes” — quantities that are suppressed without apparent recourse to a symmetry explanation. Magic zeroes and their corresponding hidden symmetries may shed new light on parametric hierarchies in the Standard Model and beyond. We identify the hidden symmetry responsible for a recently-discovered magic zero, the vanishing of the putative leading contribution to the anomalous dipole moments of the muon upon integrating out weak doublet and singlet vector-like fermions. Some of the tools involved — spurion analysis leveraging discrete symmetries of the free theory, field redefinitions, spectator fields, and non-supersymmetric non-renormalization theorems — may prove useful in the hunt for new magic zeroes and their hidden symmetries.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Low-pass spectral analysis of time-resolved serial femtosecond crystallography data

Low-pass spectral analysis (LPSA) is a recently developed dynamics retrieval algorithm showing excellent retrieval properties when applied to model data affected by extreme incompleteness and stochastic weighting. In this work, we apply LPSA to an experimental time-resolved serial femtosecond crystallography (TR-SFX) dataset from the membrane protein bacteriorhodopsin (bR) and analyze its parametric sensitivity. While most dynamical modes are contaminated by nonphysical high-frequency features, we identify two dominant modes, which are little affected by spurious frequencies. The dynamics retrieved using these modes shows an isomerization signal compatible with previous findings. We employ synthetic data with increasing timing uncertainty, increasing incompleteness level, pixel-dependent incompleteness, and photon counting errors to investigate the root cause of the high-frequency contamination of our TR-SFX modes. By testing a range of methods, we show that timing errors comparable to the dynamical periods to be retrieved produce a smearing of dynamical features, hampering dynamics retrieval, but with no introduction of spurious components in the solution, when convergence criteria are met. Using model data, we are able to attribute the high-frequency contamination of low-order dynamical modes to the high levels of noise present in the data. Finally, we propose a method to handle missing observations that produces a substantial dynamics retrieval improvement from synthetic data with a significant static component. Reprocessing of the bR TR-SFX data using the improved method yields dynamical movies with strong isomerization signals compatible with previous findings.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Spin Waves in Dirac Semimetal Ca 0.6⁢ Sr 0.4⁢ MnSb 2 Investigated with Neutrons by the Diffraction Method

The tunability of Dirac semimetals with antiferromagnetic Mn layers is of great interest. The observed sign change of interlayer magnetic coupling between CaMnBi 2 and SrMnBi 2 suggests ionic substitution as a potential tuning mechanism. If so, novel behavior near the compensation point could be expected. Here, to explore this, we study a mixed-cation analog, Ca 0.6⁢ Sr 0.4⁢ MnSb 2 , where Bi is replaced by Sb. Conventional inelastic neutron scattering is impractical due to the small crystal size (m ≈ 0.28 g) available for compositional studies; however, we find that using a neutron diffractometer with a wide-angle area detector we can obtain a good quality spin-wave signal, which is shaped by energy-momentum conservation and retains spectroscopic information even without direct energy analysis. Spin-wave modeling reveals an interlayer coupling quantitatively similar to SrMnSb 2 , indicating it is not directly tuned by ionic size and that the sign change in Bi-based compounds likely arises from the observed change in lattice symmetry. Beyond this key insight, our results present an efficient method for parametric and compositional studies of spin dynamics in small crystals.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Thermodynamic Analysis of Thermo-vacuum Clothes Drying Operation

Clothes drying accounts for a significant amount of energy consumed in residential and commercial sectors. A thermal vacuum clothes drying technology (TVCD) is proposed as an advanced clothes dryer which can significantly reduce the energy requirements by expediting the drying process. In the conventional convective clothes dryer, hot dry air is introduced into the drum which gets in direct contact to dry the clothes. This process is energy inefficient since the significant amount of heat and the water carried out with the exhaust stream are wasted. In contrast to the conventional convective drying technique, the drying mechanism of TVCD is through nucleate boiling at low temperature due to reduced vessel pressure. The process is not only efficient but also reduces the required time for drying. The objective of this paper is to develop a comprehensive thermodynamic model to predict the transient drying process of TVCD. The three-stage system model including the detail analysis of individual components is capable of simulating the variation of the water content in the textile under various operational conditions. The preliminary results show that the drying time of 3 lb textile from 70% to 2.5% in TVCD is approximately four times less than the time required in the conventional clothes dryer. Parametric studies are carried out to understand the effect of operation conditions and component geometry on the system performance. Energy consumption of the system is also analyzed.

Yang, Cheng-Min↗
DOE OSTI2023

Comparative Analysis of Confidence Metrics for Nuclear Criticality Safety

Nuclear criticality safety standards provide guidance on the requirements and recommendations to establish confidence in computerized model results used to support operation with fissionable materials. By design, the guidance is not prescriptive, leaving the analysts free to determine how various sources of uncertainties are to be statistically aggregated. This report compares the analyses and key assumptions behind four notable methodologies documented in the nuclear criticality safety literature: the parametric, nonparametric, Whisper, and TSURFER methodologies. Because of the involved use of statistics entangled with heuristic recipes, the results of these methodologies are often difficult to interpret. Also, they are augmented by additional large administrative margins, eliminating the incentive to understand their differences. With the new resurgent wave of advanced nuclear systems focused on economizing operation—including advanced reactors, fuel cycles, and fuel concepts—there is a strong need to develop a clear understanding of uncertainties and their fusion methodologies to reduce uncertainties in a scientifically defensible manner. This report offers a deep dive into the various assumptions of the four noted methodologies, their adequacy, and their limitations, to provide guidance on developing confidence for the emergent nuclear systems. These systems are expected to be challenged by the scarcity of experimental data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Characterization and Valuation of the Uncertainty of Calibrated Parameters in Microsimulation Decision Models

We evaluated the implications of different approaches to characterize the uncertainty of calibrated parameters of microsimulation decision models (DMs) and quantified the value of such uncertainty in decision making. We calibrated the natural history model of CRC to simulated epidemiological data with different degrees of uncertainty and obtained the joint posterior distribution of the parameters using a Bayesian approach. We conducted a probabilistic sensitivity analysis (PSA) on all the model parameters with different characterizations of the uncertainty of the calibrated parameters. We estimated the value of uncertainty of the various characterizations with a value of information analysis. We conducted all analyses using high-performance computing resources running the Extreme-scale Model Exploration with Swift (EMEWS) framework. The posterior distribution had a high correlation among some parameters. The parameters of the Weibull hazard function for the age of onset of adenomas had the highest posterior correlation of -0.958. When comparing full posterior distributions and the maximum-a-posteriori estimate of the calibrated parameters, there is little difference in the spread of the distribution of the CEA outcomes with a similar expected value of perfect information (EVPI) of $\$$653 and $\$$685, respectively, at a willingness-to-pay (WTP) threshold of $\$$66,000 per quality-adjusted life year (QALY). Ignoring correlation on the calibrated parameters’ posterior distribution produced the broadest distribution of CEA outcomes and the highest EVPI of $\$$809 at the same WTP threshold. Different characterizations of the uncertainty of calibrated parameters affect the expected value of eliminating parametric uncertainty on the CEA. Ignoring inherent correlation among calibrated parameters on a PSA overestimates the value of uncertainty.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

AN ACID BAKING APPROACH TO ENHANCE RARE EARTH ELEMENT RECOVERY FROM BITUMINOUS COAL SOURCES

Rare earth elements (REE) are a group of 17 elements typically classified as light and heavy rare earth elements, which play a crucial role in developing the latest technologies for energy, defense, and medical sectors. Even though REEs have been found in more than 200 minerals, only bastnaesite, monazite and xenotime are commercially exploited for REE extraction. However, the recent exponential increase in REE demand has spurred countries such as the United States into research for the extraction of REEs from secondary sources such as coal, acid mine drainage, and coal ash. Several coal sources (e.g., Fire Clay seam coal) across the United States have been identified to contain elevated concentrations of rare earth elements (>600 ppm), and various researchers have investigated the feasibility of both physical and hydrometallurgical extraction techniques for rare earth concentration and subsequent extraction. However, both physical and direct leaching were concluded to be inefficient for RE beneficiation and extraction due to low recoveries. Alternatively, thermal treatment provides promising means for RE recovery from bituminous coal sources. The positive impact of thermal treatment/calcination was established to be due to the decarbonization and dehydroxylation of the clays, which released entrapped rare earth elements within the dominant minerals and converted them into an acid-soluble form. Nonetheless, the improvement in recovery was limited to the light REEs (LREEs) with an insignificant increase in the heavy REEs (HREE). It was demonstrated that the light and heavy REEs in the material were associated with difficult-to-leach minerals such as monazite, xenotime, and zircon, which were not decomposed by simple calcination due to their high thermal stability. Hence, roasting in the presence of chemicals was necessary to ensure the decomposition of those REE-containing minerals. As such, this study was focused on the acid baking treatment of bituminous coals with an aim to enhance REE recovery, especially HREEs. Based on the presence of REE minerals like monazite and xenotime, three pre-leach treatment methods, i.e., 1) roasting, 2) direct acid baking, and 3) acid baking after roasting were investigated. Roasting tests at 600 ⁰C revealed that the recovery of light REEs (LREEs) was enhanced while the recovery of HREEs remained relatively unaffected. LREE and HREE recovery values of 38.3% and 21.3%, respectively, were achieved using a 50 g/L sulfuric acid solution at 5% solid concentration and a solution temperature of 75 ⁰C for 2 hours. Comparatively, direct acid baking at 250 °C provided substantial increased LREE and HREE recovery values to approximately 49.4% and 53.0%, respectively, using an equivalent acid dosage. Recoveries were maximized to 77.0% and 79.6% for LREE and HREE, respectively, by roasting followed by acid baking. Similar results were obtained from the treatment of a second bituminous coal source. Due to strong correlations between REE and Al recovery values, tests were performed on kaolinite and illite, which were prominent clay minerals within the source coals. These experiments revealed that the REE recovery improvements were likely a result of dehydroxylation of clays and subsequent release and decomposition of REE-bearing minerals such as monazite, xenotime and zircon. Subsequently, a parametric study was conducted to identify the impact of acid baking parameters on rare earth element recovery. The factors investigated using a three-level statistical experimental program were acid baking time, acid solution concentration, baking temperature, and acid solution-to-solids ratio which were found to significantly impact REE and contaminant element (Al, Fe, and Ca) recovery. An increase in baking temperature up to around 250 ⁰C improved the light and heavy REE recovery values by more than 50 absolute percentage points relative to performances achieved when direct leaching. As aforementioned, acid baking was needed to both decompose the clay minerals and liberate the REE minerals, which allowed access for the acid to solubilize the REEs. Acid concentration of the solution used for acid baking was studied as a means of minimizing the amount of acid needed to achieve a target REE recovery. However, thermo-gravimetric and differential scanning calorimetry analysis (TGA-DSC) of sulfuric acid under oxidizing atmosphere revealed that the addition of water decreased the evaporation temperature, which explains the lower REE recovery values obtained when using acid concentrations less than 100%. Using pure sulfuric acid at an acid-to-solid ratio of 0.8:1 resulted in recovery values of around 70% for both LREEs and HREEs. The decomposition reaction time was relatively quick with 65% of the TREEs recovered within the first 10 minutes. Following the identification of optimum operating conditions through the parametric study, a systematic leaching study was carried out to examine the impact of leaching parameters, such as solid-to-liquid (S/L) ratio, temperature, and time, on REE recovery using acid baking conditions of 1:1 acid to solids ratio at 250 °C for 30 minutes. The solid-to-liquid ratio was varied from 1-20% by weight at 25 °C, 50 °C, and 75 °C solution temperatures. The results indicated that reaction time and solution temperature considerably impacted the recovery of heavy and light REEs. Interestingly, LREE recovery reduced from 68% at 5% S/L to 58% at 20%S/L, whereas HREE recovery of 78% remained unaffected. The decrease in the LREE recovery was determined to be due to La and Ce precipitation, likely through isomorphic substitution with calcium in gypsum. The kinetic data indicated that 67% LREE and 77% HREE recovery could be obtained within the first 15 minutes of the reaction, suggesting fast reaction kinetics. Furthermore, raising the solution temperature from 25 °C to 75 °C increased the LREE and HREE recovery from 60% and 32% to 67% and 79%, respectively. The kinetic modeling results demonstrated that the rate-limiting step in the LREE dissolution was diffusion and chemical reaction, whereas the HREE extraction was controlled by only chemical reaction. The leaching study concluded that using 20% S/L at 75 °C for 15 minutes maximized LREE and HREE recovery. The lab-scale precipitation study showed that Fe and Al in solution could be removed at pH 4.5 followed by REE precipitation at pH 6.0 using 6 mol/L NaOH. Finally, the bench-scale data was used to develop a process flowsheet for REE recovery from low-grade bituminous coal sources using acid baking. Finally, based on the proposed flowsheet, a concentrated RE-cake obtained through selective precipitation at pH 6.5 was re-leached using HCl at pH 1.5. The resultant leachate was used to identify the impact of various operating parameters on REE recovery and purity with an aim to maximize REE precipitation efficiency while minimizing the oxalic acid dosage. The operating parameters for this investigation were oxalic acid dosage, iron (III) contamination, solution pH and temperature. The resultant model suggested that oxalic acid dosage and reaction pH are the most significant factors for the REE precipitation efficiency, followed by the interaction of oxalic dosage and Fe concentration. Test results indicated that increasing the oxalic acid concentration from 0g/L to 80g/L improved the REE precipitation efficiency from approximately 4.2% to 95.0%. Furthermore, raising the solution pH from 0.5 to 2.5 considerably enhanced the precipitation efficiency from 0.0% to 98.9%. A solution temperature elevation decreased REE recovery, which indicated an exothermic reaction between REEs and oxalate anions. Finally, a high level of Fe contamination adversely impacted REE precipitation efficiency. The speciation analysis revealed that the dominant iron species in the solution system were Fe-(C₂O₄)₃³⁻, Fe-(C₂O₄)²⁻, and Fe-(C₂O₄)⁺, which consumed the majority of the oxalate anions

rare earth elements, acid baking, high-temperature↗
DOE OSTI2025

A Regularized Variance-Reduced Modified Extragradient Method for Stochastic Hierarchical Games

We consider an N -player hierarchical game in which the i th player’s objective comprises of an expectation-valued term, parametrized by rival decisions, and a hierarchical term. Such a framework allows for capturing a broad range of stochastic hierarchical optimization problems, Stackelberg equilibrium problems, and leader-follower games. We develop an iteratively regularized and smoothed variance-reduced modified extragradient framework for iteratively approaching hierarchical equilibria in a stochastic setting. We equip our analysis with rate statements, complexity guarantees, and almost-sure convergence results. We then extend these statements to settings where the lower-level problem is solved inexactly and provide the corresponding rate and complexity statements. Our model framework encompasses many game theoretic equilibrium problems studied in the context of power markets. We present a realistic application to the study of virtual power plants, emphasizing the role of hierarchical decision making and regularization. Preliminary numerics suggest that empirical behavior compares well with theoretical guarantees.

Tikhonov regularization↗
DOE OSTI2021

Internally cooled membrane-based absorber for dehumidification and water heating: Validated model and simulation study

This report presents a comprehensive model for advanced membrane-based absorber components accounting for three separate fluid streams, one of them providing internal cooling, with complex flow patterns. The model was implemented with various liquid desiccants and their properties (including LiCl, CaCl 2 , and [emim][OMS] ionic liquid [IL]). The model has been validated using data from two laboratory prototype absorbers, one with CaCl 2 as the working fluid and one with an IL as the working fluid. Entropy analysis was further carried out to verify the model and understand effect of internal cooling on the entropy variation of three fluid streams. The developed model and codes are expected to enable detailed configuration optimization and provide in-depth understanding of three-fluid heat and mass exchanger (HMX) performance. The paper also describes parametric studies of the HMX components in utilizing the latent heat removed in space cooling to heat water and numerically explores the unique benefit of its application in a semi-open absorption heat pump water heater. The results show that a single three-fluid HMX has the potential to achieve simultaneous dehumidification and water heating efficiently at cost effectiveness, particularly in hot and humid climates.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗