Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reduce order modeling”

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,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Quantum simulation of quantum phase transitions using the convex geometry of reduced density matrices

Transitions of many-particle quantum systems between distinct phases at absolute-zero temperature, known as quantum phase transitions, require an exacting treatment of particle correlations. Here in this work, we present a general quantum-computing approach to quantum phase transitions that exploits the geometric structure of reduced density matrices. While typical approaches to quantum phase transitions examine discontinuities in the order parameters, the origin of phase transitions—their order parameters and symmetry breaking—can be understood geometrically in terms of the set of two-particle reduced density matrices (2-RDMs). The convex set of 2-RDMs provides a comprehensive map of the quantum system including its distinct phases as well as the transitions connecting these phases. Because 2-RDMs can potentially be computed on quantum computers at nonexponential cost, even when the quantum system is strongly correlated, they are ideally suited for a quantum-computing approach to quantum phase transitions. We compute the convex set of 2-RDMs for a Lipkin-Meshkov-Glick spin model on IBM superconducting-qubit quantum processors. Even though computations are limited to few-particle models due to device noise, comparisons with a classically solvable 1000-particle model reveal that the finite-particle quantum solutions capture the key features of the phase transitions including the strong correlation and the symmetry breaking.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Oxidation of 3D-printed SiC in air and steam environments

Here, the high-temperature oxidation of additively manufactured and chemically vapor infiltrated (3D-printed SiC) has been compared to chemical vapor deposited (CVD) SiC. 100-h isothermal exposures were conducted at 1425° and 1300°C at 1 atm under both dry air and steam environments. A SiC reaction tube was utilized to reduce silica volatility. After steam oxidation at 1425° and 1300°C, on the 3D-printed SiC surface, which was intrinsically rougher than the CVD surface, scales were 70%–90% thicker at the convex regions compared to concave/flat regions. In the convex regions, large cracks perpendicular to the oxidizing interface were observed. After dry air oxidation, scale thicknesses were comparable between 3D-printed SiC and CVD SiC, regardless of geometry. Finite element modeling, conducted to elucidate the relationship between SiC geometry and ß- to α-cristobalite transformation stress, determined cristobalite transformation tensile stresses to be on the order of 103 MPa during cool down, assuming a 6 vol% reduction. Compared to flat SiC substrates, tensile transformation stresses were elevated at concave regions and relaxed at convex regions. Combined with specimen mass gain (accounting for the rougher surface) of 3D-printed SiC being 15%–32% higher for 3D-printed SiC after 1300°C and 1425°C steam oxidation, the work presented concludes that the increased oxidation of 3D-printed SiC is primarily caused by tensile hoop stresses driven by oxidation volume expansion. Lastly, the efficacy of the 3D-printing method is demonstrated through the production of tristructural isotropic imbedded 3D-printed SiC fuel forms.

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

Accelerating the solution of the S N equations with highly anisotropic scattering using the Fokker-Planck approximation

The discrete ordinates method can model forward-peaked transport problems accurately. However, convergence of discrete ordinates solution can become arbitrarily slow upon use of standard iterative procedures like source iteration. Standard zeroth and first moment-based acceleration methods like nonlinear diffusion acceleration and diffusion synthetic acceleration are ineffective in accelerating such problems because these methods do not correct higher order Legendre-moments of angular flux. In this paper, we explore the idea of using the Fokker-Planck approximation as a preconditioner to accelerate the convergence of solution of forward-peaked transport problems. Finally, we observe that such preconditioning can successfully reduce the iteration count by up to four orders of magnitude, and solver runtime by up to three orders of magnitude when compared to diffusion synthetic acceleration.

42 ENGINEERING↗
DOE OSTI2022

Comparison of spherical harmonics method and discrete ordinates method for radiative transfer in a turbulent jet flame

Here, in this study, we systematically compared the accuracy and computational cost of two popular solution methods for the radiative transfer equation (RTE): the spherical harmonics method (P N ) and the discrete ordinates method (DOM). We first investigated convergence characteristics of different orders of P N and DOM in a series of 1D homogeneous configurations with varying optical thicknesses. Both solvers perform better for optically thicker cases. The accuracy of P N methods increases with its order, , but the gain in accuracy reduces with the increase in , i.e., improvement of P 7 over P 5 is less than that of P 3 over P 1 . This decreasing trend becomes more prominent as the optical thickness decreases. On the other hand, DOM’s accuracy increases almost linearly with the increase in the number of ordinates (or polar angles in this study) in all cases. While comparing the directional profile of radiative intensity, both solvers perform better when the radiative intensity is more isotropic. These solvers were then connected with a full spectrum k-distribution (FSK) spectral model and used to perform radiation-coupled simulations of a turbulent jet flame in an axi-symmetric cylindrical domain. Results obtained from P 1 to P 7 approximations for P N , and 2 x 4, 4 x 4, 4 x 8, 8 x 8 finite angles for DOM are compared with that from an optically thin model, and a reference solution from line-by-line (LBL) photon Monte Carlo (PMC) method. The choice of radiation solver shows a noticeable impact on the temperature distribution of the flame. The P N solvers lead to slightly higher radiant fractions and the DOM solvers lead to slightly lower radiant fractions than the PMC benchmark solution. Finally, the computational costs of each of these solvers are also reported and an intermittent evaluation / time blending scheme to improve the computational efficiency of radiation solvers in radiation-coupled simulations are also demonstrated.

42 ENGINEERING↗
DOE OSTI2020

Effect of ammonia addition on suppressing soot formation in methane co-flow diffusion flames

Due to issues surrounding carbon dioxide emissions from carbon-containing fuels, there is growing interest in ammonia (NH 3 ) as an alternative combustion fuel. One attractive method of burning NH 3 is to co-fire it with hydrocarbons, such as natural gas, and in this case soot formation is possible. To begin understanding the influence of NH 3 on soot formation when co-fired with hydrocarbons, soot volume fractions and mole fractions of gas-phase species were computationally and experimentally interrogated for CH 4 flames with up to 40% NH 3 by volumetric fuel fraction. Mole fractions of gas-phase species, including C 2 H 2 and C 6 H 6 , were measured with on-line electron impact mass spectrometry, and soot volume fractions were obtained via color-ratio pyrometry. The simulations employed a detailed chemical mechanism developed for capturing nitrogen interactions with hydrocarbons during combustion. The results are compared to findings in N 2 -CH 4 flames, in order to separate thermal and dilution effects from the chemical influence of NH 3 on soot formation. Experimentally, C 2 H 2 concentrations were found to decrease slightly for the NH 3 -CH 4 flames relative to N 2 -CH 4 flames, and a stronger suppression of C 6 H 6 was found for NH 3 relative to N 2 additions. The measured results show a strong suppression of soot with the addition of NH 3 , with soot concentrations reduced by over a factor of 10 with addition of up to 20% or more NH 3 by mole fraction. The model satisfactorily captured relative differences in maximum centerline C 2 H 2 , C 6 H 6 , and soot concentrations with addition of N 2 , but was unable to match measured differences in NH 3 -CH 4 flames. Furthermore, these results highlight the need for an improved understanding of fuel-nitrogen interactions with higher hydrocarbons to enable accurate models for predicting particulate emissions from NH 3 /hydrocarbon combustion.

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

Experimental Quantification of Spin–Phonon Coupling in Molecular Qubits Using Inelastic Neutron Scattering

Electronic spin superposition states enable nanoscale sensing through their sensitivity to the local environment, yet their sensitivity to vibrational motion also limits their coherence times. In molecular spin systems, chemical tunability and atomicscale resolution are accompanied by a dense, thermally accessible phonon spectrum that introduces efficient spin relaxation pathways. Despite extensive theoretical work, there is little experimental consensus on which vibrational energies dominate spin relaxation or how molecular structure controls spin−phonon coupling (SPC). We present a fully experimental method to quantify SPC coefficients by combining temperature-dependent vibrational spectra from inelastic neutron scattering with spin relaxation rates measured by electron paramagnetic resonance. We apply this framework to two model S = 1/2 systems, copper(II) phthalocyanine (CuPc) and copper(II) octaethylporphyrin (CuOEP). Two distinct relaxation regimes emerge: below 40 K, weakly coupled lattice modes below 50 cm −1 dominate, whereas above 40 K, optical phonons above ∼185 cm −1 become thermally populated and drive relaxation with SPC coefficients nearly 3 orders of magnitude larger. Structural distortions in CuOEP that break planar symmetry soften the crystal lattice and enhance anharmonic scattering but also raise the energy of stretching modes at the molecular core where the spins reside. This redistributes vibrational energy toward the molecular periphery and out of plane, ultimately reducing SPC relative to CuPc and enabling room-temperature spin coherence in CuOEP. Although our method does not provide mode-specific SPC coefficients, it quantifies contributions from distinct spectral regions and establishes a broadly applicable, fully experimental link between crystal structure, lattice dynamics, and spin relaxation.

Lohaus, Stefan H. [California Institute of Technol↗
DOE OSTI2022

Reduced models for ETG transport in the tokamak pedestal

This paper reports on the development of reduced models for electron temperature gradient (ETG) driven transport in the pedestal. Model development is enabled by a set of 61 nonlinear gyrokinetic simulations with input parameters taken from pedestals in a broad range of experimental scenarios. The simulation data have been consolidated in a new database for gyrokinetic simulation data, the multiscale gyrokinetic database (MGKDB), facilitating the analysis. The modeling approach may be considered a generalization of the standard quasilinear mixing length procedure. The parameter η, the ratio of the density to temperature gradient scale length, emerges as the key parameter for formulating an effective saturation rule. With a single order-unity fitting coefficient, the model achieves an error of 15%. A similar model for ETG particle flux is also described. We also present simple algebraic expressions for the transport informed by an algorithm for symbolic regression.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Medium amplitude parallel superposition (MAPS) rheology. Part 2: Experimental protocols and data analysis

An experimental protocol is developed to directly measure the new material functions revealed by medium amplitude parallel superposition (MAPS) rheology. This protocol measures the medium amplitude response of a material to a simple shear deformation composed of three sine waves at different frequencies, revealing a rich dataset consisting of up to 19 measurements of the third-order complex modulus at distinct three-frequency coordinates. We discuss how the choice of input frequencies influences the features of the MAPS domain studied by the experiment. A polynomial interpolation method for reducing the bias of measured values from spectral leakage and reducing variance due to noise is discussed, including a derivation of the optimal range of amplitudes for the input signal. This leads to the conclusion that conducting the experiment in a stress-controlled fashion possesses a distinct advantage to the strain-controlled mode. The experimental protocol is demonstrated through measurements of the MAPS response of a model complex fluid: a surfactant solution of wormlike micelles. The resulting dataset is indeed large and feature-rich, while still acquired in a time comparable to similar medium amplitude oscillatory shear (MAOS) experiments. We demonstrate that the data represent measurements of an intrinsic material function by studying its internal consistency, compatibility with low-frequency predictions for Coleman–Noll simple fluids, and agreement with data obtained via MAOS amplitude sweeps. Finally, the data are compared to predictions from the corotational Maxwell model to demonstrate the power of MAPS rheology in determining whether a constitutive model is consistent with a material’s time-dependent response.

Lennon, Kyle R. (ORCID:0000000212515461)↗
DOE OSTI2021

Comprehensive Material Characterization and Simultaneous Model Calibration for Improved Computational Simulation Credibility

Computational simulation is increasingly relied upon for high-consequence engineering decisions, and a foundational element to solid mechanics simulations is a credible material model. Our ultimate vision is to interlace material characterization and model calibration in a real-time feedback loop, where the current model calibration results will drive the experiment to load regimes that add the most useful information to reduce parameter uncertainty. The current work investigated one key step to this Interlaced Characterization and Calibration (ICC) paradigm, using a finite load-path tree to incorporate history/path dependency of nonlinear material models into a network of surrogate models that replace computationally-expensive finite-element analyses. Our reference simulation was an elastoplastic material point subject to biaxial deformation with a Hill anisotropic yield criterion. Training data was generated using either a space-filling or adaptive sampling method, and surrogates were built using either Gaussian process or polynomial chaos expansion methods. Surrogate error was evaluated to be on the order of 10 ⁻5 and 10 ⁻3 percent for the space-filling and adaptive sampling training data, respectively. Direct Bayesian inference was performed with the surrogate network and with the reference material point simulator, and results agreed to within 3 significant figures for the mean parameter values, with a reduction in computational cost over 5 orders of magnitude. These results bought down risk regarding the surrogate network and facilitated a successful FY22-24 full LDRD proposal to research and develop the complete ICC paradigm.

36 MATERIALS SCIENCE↗
DOE OSTI2023

The Development of the Pooled Rideshare Acceptance Model (PRAM)

Due to the advancements in real-time information communication technologies and sharing economies, rideshare services have gained significant momentum by offering dynamic and/or on-demand services. Rideshare service companies evolved from personal rideshare, where riders traveled solo or with known individuals, into pooled rideshare (PR), where riders can travel with one to multiple unknown riders. Similar to other shared economy services, pooled rideshare is beneficial as it efficiently utilizes resources, resulting in reduced energy usage, as well as reduced costs for the riders. However, previous research has demonstrated that riders have concerns about using pooled rideshare, especially regarding personal safety. A U.S. national survey with 5385 participants was used to understand human factor-related barriers and user preferences to develop a novel Pooled Rideshare Acceptance Model (PRAM). This model used a covariance-based structural equation model (CB-SEM) to identify the relationships between willingness to consider PR factors (time/cost, privacy, safety, service experience, and traffic/environment) and optimizing one’s experience of PR factors (vehicle technology/accessibility, convenience, comfort/ease of use, and passenger safety), resulting in the higher-order factor trust service. We examined the factors’ relative contribution to one’s willingness/attitude towards PR and user acceptance of PR. Privacy, safety, trust service, and convenience were statistically significant factors in the model, as were the comfort/ease of use factor and the service experience, traffic/environment, and passenger safety factors. The only two non-significant factors in the model were time/cost and vehicle technology/accessibility; it is only when a rider feels safe that individuals then consider the additional non-significant variables of time, cost, technology, and accessibility. Privacy, safety, and service experience were factors that discouraged the use of PR, whereas the convenience factor greatly encouraged the acceptance of PR. Despite the time/cost factor’s lack of significance, individual items related to time and cost were crucial when viewed within the context of convenience. This highlights that while user perceptions of privacy and safety are paramount to their attitude towards PR, once safety concerns are addressed, and services are deemed convenient, time and cost elements significantly enhance their trust in pooled rideshare services. This study provides a comprehensive understanding of user acceptance of PR services and offers actionable insights for policymakers and rideshare companies to improve their services and increase user adoption.

dynamic rideshare system↗
DOE OSTI2020

Energy demand science for a decarbonized society in the context of the residential sector

To develop a decarbonized society, two contradictory requirements must be met: (1) reducing energy demand and (2) creating flexibility in energy demand in order to respond to fluctuations in renewable electricity generation. To help meet these requirements, conventional energy efficiency studies should be extended to incorporate “energy demand science.” This paper presents a definition of “energy demand science” and then reviews the related history and research questions of energy demand science in the context of the residential sector. It then examines three key areas that must be integrated into the next-generation energy demand science: (1) energy demand measurement with detailed granularity and analysis using cutting-edge technology, (2) energy demand modeling that helps clarify the formation mechanism of energy demand, and (3) identification of the factors that influence people's decision making, which represents typical human-dimension research. Energy demand science consists of technical, human, natural environment, demographic, and land-use dimensions, and their integration is key for the establishment of a decarbonized society.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Sparsity-Independent Lyapunov Exponent in the Sachdev-Ye-Kitaev Model

The saturation of a recently proposed universal bound on the Lyapunov exponent has been conjectured to signal the existence of a gravity dual. This saturation occurs in the low-temperature limit of the dense Sachdev-Ye-Kitaev (SYK) model, N Majorana fermions with q body ( q > 2 ) infinite-range interactions. We calculate certain out-of-time-order correlators (OTOCs) for N ≤ 64 fermions for a highly sparse SYK model and find no significant dependence of the Lyapunov exponent on sparsity up to near the percolation limit where the Hamiltonian breaks up into blocks. This provides strong support to the saturation of the Lyapunov exponent in the low-temperature limit of the sparse SYK. A key ingredient to reaching N = 64 is the development of a novel quantum spin model simulation library that implements highly optimized matrix-free Krylov subspace methods on graphical processing units. This leads to a significantly lower simulation time as well as vastly reduced memory usage over previous approaches, while using modest computational resources. Strong sparsity-driven statistical fluctuations require both the use of a much larger number of disorder realizations with respect to the dense limit and a careful finite size scaling analysis. The saturation of the bound in the sparse SYK points to the existence of a gravity analog that would enlarge substantially the number of field theories with this feature. Published by the American Physical Society 2024

Physics↗
DOE OSTI2022

Efficient Region of Attraction Characterization for Control and Stabilization of Load Tap Changer Dynamics

In this article, we study the monitoring and control of long-term voltage stability considering load tap changer (LTC) dynamics. We show that under generic conditions, the LTC dynamics admit a unique stable equilibrium. For the stable equilibrium, we characterize an explicit inner approximation of its largest region of attraction (ROA). Compared to existing results, the computational complexity of the ROA characterization is drastically reduced. We propose a quadratically constrained linear program formulation for the ROA characterization problem. In addition, we formulate a second-order cone program for online voltage stability monitoring and control exploiting the proposed ROA characterization. Finally, we demonstrate the efficacy of the proposed formulations on the ROA characterization and stability monitoring and control using a standard IEEE test system.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

US Fallout Shelter

Being indoors (sheltering) reduces radiation exposures resulting from nuclear fallout or a power-plant accident. However, most US building types, particularly non-residential buildings, lack quantitative estimates. We provide here a high-level modeling analysis of modern US building protection based on a novel radiation protection building attribute taxonomy, a new building protection model, compilation of prior experimental results, and nationally representative building survey data. This approach provides a consistent, quantitative understanding of US building protection and, when combined with the distribution of people among different buildings, can assist in selecting emergency response strategies. We find that indoor radiation protection varies by orders of magnitude. Most people in non-residential buildings are adequately protected, particularly below ground or in the building center. Residential buildings are less protective and those with lighter weight walls lack adequate, above ground protection. This work also highlights key areas where further investigation will improve the current results. These include the improving the understanding of residential basement protection, additional experimental data on non-residential buildings (particularly schools and industrial buildings), and the frequency of brick veneer and stucco exterior residential buildings.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2024

Projective Integral Updates for High-Dimensional Variational Inference

Variational inference is an approximation framework for Bayesian inference that seeks to improve quantified uncertainty in predictions by optimizing a simplified distribution over parameters to stand in for the full posterior. Capturing model variations that remain consistent with training data enables more robust predictions by reducing parameter sensitivity. This work introduces a fixed-point optimization for variational inference that is applicable when every feasible log density can be expressed as a linear combination of functions from a given basis. In such cases, the optimizer becomes a fixed-point of projective integral updates. When the basis spans univariate quadratics in each parameter, the feasible distributions are Gaussian mean-fields and the projective integral updates yield quasi-Newton variational Bayes (QNVB). Other bases and updates are also possible. Since these updates require high-dimensional integration, this work begins by proposing an efficient quasirandom sequence of quadratures for mean-field distributions. Each iterate of the sequence contains two evaluation points that combine to correctly integrate all univariate quadratic functions and, if the mean-field factors are symmetric, all univariate cubics. More importantly, averaging results over short subsequences achieves periodic exactness on a much larger space of multivariate polynomials of quadratic total degree. The corresponding variational updates require four loss evaluations with standard (not second-order) backpropagation to eliminate error terms from over half of all multivariate quadratic basis functions. Furthermore, this integration technique is motivated by first proposing stochastic blocked mean-field quadratures, which may be useful in other contexts. A PyTorch implementation of QNVB allows for better control over model uncertainty during training than competing methods. Experiments demonstrate superior generalizability for multiple learning problems and architectures.

Gaussian mean-field↗
DOE OSTI2026

Designing the Protocols for Programmable Ammonia Catalysis

Programmable catalysis can provide a more energy-efficient and cost-effective route to enhancing commercial ammonia production, a key process in the advancement of renewable energy technologies and the manufacture of fertilizers and basic chemicals. This work explores the computational discovery of optimal forcing protocols to drive such dynamic catalysis models. By employing matrix-free time-stepper methods, coupled with an optimization approach, that integrates Bayesian optimization with a Bayesian continuation strategy to efficiently discover the periodic steady states of such periodically forced systems, we enable the discovery of complex optimal catalyst strain waveforms, while ensuring robust solver convergence. We demonstrate the flexibility of our approach to discover optimized forcing protocols under varying physical constraints on strain modulation or other catalyst operating parameters. We show that these can have a temporal structure more complex than simple step functions. In order to detect undesirable catalytic loops that may correlate with overall reduced performance, we perform a study using graph-theoretical analysis to investigate the dynamics of catalytic kinetic networks formed.

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

Coupling multi-fidelity xRAGE with machine learning for graded inner shell design optimization in double shell capsules

Bayesian optimization has shown promise for the design optimization of inertial confinement fusion targets. Specifically, in Vazirani et al. [Phys. Plasmas 28 , 122709 (2021)], optimal designs for double shell capsules with graded inner shells were identified using one-dimensional xRAGE simulation yield calculations. While the machine learning models were able to accurately learn and predict one-dimensional simulation target performance, using simulations with higher fidelity would improve design optimization and better match with the expected experimental performance. However, higher fidelity physics modeling, i.e., two-dimensional xRAGE simulations, requires significantly larger computational time/cost, usually at least an order of magnitude, in comparison with one-dimensional simulations. This study presents a multi-fidelity Bayesian optimization, in which the machine learning model leverages low-fidelity (one-dimensional xRAGE) and high-fidelity (two-dimensional xRAGE) simulations to more accurately predict “pre-shot” target performance with respect to the expected experimental performance. By building a multi-fidelity Bayesian optimization framework coupled with xRAGE, the low-fidelity and high-fidelity simulations are able to inform one another, such that we have: (1) improved physics modeling in comparison with using low-fidelity simulations alone, (2) reduced computational time/cost in comparison with using high-fidelity simulations alone, and (3) more confidence in the expected performance of optimized targets during real-world experiments. In the future, we plan to use this robust multi-fidelity Bayesian optimization methodology to expedite the design of graded inner shells further and eventually full capsules as a part of the current double shell campaign at the National Ignition Facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Hanford Tank Waste Matrix Impact on Ion Exchange Performance Using Crystalline Silicotitanate

The removal of radiocesium from Hanford tank waste supernate is a critical step in preparing feed for low-activity waste immobilization. This study evaluated cesium ion exchange performance using crystalline silicotitanate (CST) media in a series of tests designed to evaluate the influence of waste matrix variability on capacity and kinetics. Tank waste supernate subsampled from five Hanford double-shell tanks encompassed a range of sodium, hydroxide, nitrate, and nitrite concentrations in order to assess the impact of feed variability on the performance of the ion exchange system. Both equilibrium and dynamic ion exchange tests were conducted to quantify cesium distribution coefficients and breakthrough behavior under prototypic operating conditions. Results indicated that effective cesium capacity varied by up to a factor of five across the matrices tested, with higher sodium concentrations significantly reducing uptake. Kinetic behavior was similarly matrix-dependent, with solution viscosity contributing to a twofold variation in mass-transfer rates. These results demonstrate the strong dependence of CST ion exchange performance on waste composition and must be incorporated into predictive models for future treatment system design and optimization.

Westesen, Amy M.↗