Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reduced order systems”

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 289 records · Page 16

Deep learning of dynamically responsive chemical Hamiltonians with semiempirical quantum mechanics

Conventional machine-learning (ML) models in computational chemistry learn to directly predict molecular properties using quantum chemistry only for reference data. While these heuristic ML methods show quantum-level accuracy with speeds several orders of magnitude faster than traditional quantum chemistry methods, they suffer from poor extensibility and transferability; i.e., their accuracy degrades on large or new chemical systems. Incorporating quantum chemistry frameworks into the ML models directly solves this problem. Here we take the structure of semiempirical quantum mechanics (SEQM) methods to construct dynamically responsive Hamiltonians. SEQM methods use empirical parameters fitted to experimental properties to construct reduced-order Hamiltonians, facilitating much faster calculations than ab initio methods but with compromised accuracy. By replacing these static parameters with machine-learned dynamic values inferred from the local environment, we greatly improve the accuracy of the SEQM methods. Trained on molecular energies and atomic forces, these dynamically generated Hamiltonian parameters show a strong correlation with atomic hybridization and bonding. Trained with only about 60,000 small organic molecular conformers, the resulting model retains interpretability, extensibility, and transferability when testing on much larger chemical systems and predicting various molecular properties. Overall, this work demonstrates the virtues of incorporating physics-based descriptions with ML to develop models that are simultaneously accurate, transferable, and interpretable.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Optimal, Reliable Building-Integrated Energy Storage (Cooperative Research and Development Final Report)

The team will advance the commercial readiness of behind-the-meter (BTM) energy storage (ES) systems by employing health-conscious controls that guarantee lifetime and optimize the ES system's value stream when integrated with onsite renewable energy generation. Specifically, the team will develop ES controls that increase the net present value (NPV) of photovoltaics (PV) by 50% in markets where net-metering policies are being replaced by variable electricity pricing structures. The team will also reduce the risk of achieving a 10-year ES warranty lifetime by at least one order of magnitude. The successful two-year project will develop the controls and system enabling Eaton to commercialize the technology by 2021. The developments achieved through this project may enable wide scale adoption of stationary energy storage benefitting the public.

25 ENERGY STORAGE↗

MOSCATO Solver Development and Integration Plan

During FY21, we conducted ongoing development work for the MOSCATO (Molten Salt Chemistry and Transport) solver. The code development work primarily consisted of transitioning capabilities from the original version of the solver, which was written in OpenFOAM, into Nek5000. In doing so, a fast, highly parallelizable solver was created that is capable of complex chemistry and corrosion simulations for engineering-scale molten salt systems. The Nek5000 version of MOSCATO is now fully featured and capable of higher-fidelity simulations than were previously possible. Demonstration cases including a thermal convection loop have been simulated to test these new capabilities. Although capable of large-scale simulations, MOSCATO is not well-suited to parametric studies of complete reactor geometries. These types of simulations are instead better handled by reduced-order modeling codes such as ORNL’s Mole code. Reduced-order simulation tools like Mole, however, are dependent on high fidelity correlations to account for complex, coupled three-dimensional phenomena that they do not directly simulate. Tools such as MOSCATO must therefore be used to create these correlations, as suitable empirical relationships are not available for most molten salt systems. Toward that end, we used the Nek-derived version of MOSCATO to create new mass transfer correlations for three relevant cases including tubular, tube bank, and subchannel geometries. These new correlations are more accurate than any existing ones and can be readily integrated into any reduced-order modeling tools that are targeting full-scale MSR simulations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Broad-spectrum CRISPR-mediated inhibition of SARS-CoV-2 variants and endemic coronaviruses in vitro

A major challenge in coronavirus vaccination and treatment is to counteract rapid viral evolution and mutations. Here we demonstrate that CRISPR-Cas13d offers a broad-spectrum antiviral (BSA) to inhibit many SARS-CoV-2 variants and diverse human coronavirus strains with >99% reduction of the viral titer. We show that Cas13d-mediated coronavirus inhibition is dependent on the crRNA cellular spatial colocalization with Cas13d and target viral RNA. Cas13d can significantly enhance the therapeutic effects of diverse small molecule drugs against coronaviruses for prophylaxis or treatment purposes, and the best combination reduced viral titer by over four orders of magnitude. Using lipid nanoparticle-mediated RNA delivery, we demonstrate that the Cas13d system can effectively treat infection from multiple variants of coronavirus, including Omicron SARS-CoV-2, in human primary airway epithelium air-liquid interface (ALI) cultures. Our study establishes CRISPR-Cas13 as a BSA which is highly complementary to existing vaccination and antiviral treatment strategies.

60 APPLIED LIFE SCIENCES↗

A fractional calculus framework for open quantum dynamics: From Liouville to Lindblad to memory kernels

Open quantum systems exhibit dynamics ranging from unitary evolution to irreversible dissipation. While the Gorini–Kossakowski–Sudarshan–Lindblad equation uniquely characterizes Markovian completely positive and trace-preserving (CPTP) evolution, many physical platforms display non-Markovian features such as algebraic relaxation and coherence backflow. Fractional calculus provides a natural way to model such long-memory behavior through power-law temporal kernels introduced by fractional time derivatives. Here, we develop a unified framework that embeds fractional master equations within the broader hierarchy of open-system formalisms. The fractional equation forms a structured subclass of memory-kernel models, reduces to the Lindblad form at unit order, and, through Bochner–Phillips subordination, admits a CPTP representation as an average over Lindblad semigroups. Its resolvent structure further connects fractional dynamics to established non-Markovian approaches, including Nakajima–Zwanzig kernels and hierarchical equations of motion, providing a compact surrogate for long-memory effects. This formulation positions fractional calculus as a rigorous and practical language for modeling non-Markovian quantum dynamics in chemical physics and physical chemistry, providing a CPTP-preserving, computationally efficient surrogate for structured condensed-phase environments where long-time memory and dissipation play a central role.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Preconditioning for Hyper-reduction in Reduced Order Models

Many projection-based reduced order models (pROM) that utilize the governing equation and data to accelerate physical simulations can be applied to nonlinear dynamical systems. To avoid full order model (FOM) scale update for each time step, hyperreduction techniques are developed to sample high dimensional nonlinear terms. Our study aims to investigate if preconditioning the least-squares problem used for the nonlinear approximation can improve the robustness of the condition number of the problem while achieving high accuracy. In our study, we use the row-normalization matrix motivated by the Christoffel function as the preconditioner and solve the corresponding weighted least-squares problem. Numerical results for Lagrangian hydrodynamics examples are analyzed to explore how the preconditioner works compared to existing hyper-reduction techniques.

97 MATHEMATICS AND COMPUTING↗

Towards Compact Neural Networks via End-to-End Training: A Bayesian Tensor Approach with Automatic Rank Determination

Post-training model compression can reduce the inference costs of deep neural networks, but uncompressed training still consumes enormous hardware resources and energy. To enable low-energy training on edge devices, it is highly desirable to directly train a compact neural network from scratch with a low memory cost. Low-rank tensor decomposition is an effective approach to reduce the memory and computing costs of large neural networks. However, directly training low-rank tensorized neural networks is a very challenging task because it is hard to determine a proper tensor rank a priori, and the tensor rank controls both model complexity and accuracy. Here, this paper presents a novel end-to-end framework for low-rank tensorized training. We first develop a Bayesian model that supports various low-rank tensor formats (e.g., CANDECOMP/PARAFAC, Tucker, tensor-train, and tensor-train matrix) and reduces neural network parameters with automatic rank determination during training. Then we develop a customized Bayesian solver to train large-scale tensorized neural networks. Our training methods shows orders-of-magnitude parameter reduction and little accuracy loss (or even better accuracy) in the experiments. On a very large deep learning recommendation system with over 4.2 ×10 9 model parameters, our method can reduce the parameter number to 1.6 ×10 5 automatically in the training process (i.e., by 2.6 ×10 4 times) while achieving almost the same accuracy. Code is available at https://github.com/colehawkins/bayesian-tensor-rank-determination.

compact neural networks↗

Pressurized Gas-Driven Elemental Redistribution Enables Ultrastable PtNi Catalysts for Heavy-Duty Vehicles

Gas-driven element redistribution, characterized by the preferential enrichment of one element at the surface relative to the bulk, is frequently observed in multicomponent alloys. Using L1 0 -ordered PtNi as a model system, we reveal that gas pressure plays a critical role in governing adsorption-driven surface composition during annealing in reducing gases: low pressure favors Pt surface segregation, while high pressure facilitates Ni surface enrichment. In this study, we developed a high-pressure nitriding (HPN) strategy that modulates the surface structure and composition of PtNi catalysts. The resulting HPN-PtNi exhibits enhanced performance and durability in membrane electrode assemblies for heavy-duty fuel cell applications, maintaining a high current density of 1.19 A cm −2 at 0.7 V after 90,000 voltage cycles. Through a combination of experimental and theoretical analyses, we reveal that the HPN process forms additional stabilizing Ni−N bonds and induces elemental redistribution with Ni surface enrichment and a Ni-deficient Pt subsurface. These modifications alter the atomic coordination environment of the ordered PtNi phase. This work presents a generalizable strategy to design robust and high-performing Pt-based catalysts by controlling gas-pressure-driven elemental redistribution and dopant incorporation.

30 DIRECT ENERGY CONVERSION↗

Accelerating structural dynamics simulations with localised phenomena through matrix compression and projection‐based model order reduction

In this work, a novel approach is introduced for accelerating the solution of structural dynamics problems in the presence of localised phenomena, such as cracks. For this category of problems, conventional projection-based Model Order Reduction (MOR) methods are either limited with respect to the range of system configurations that can be represented or require frequent solutions of the Full Order Model (FOM) to update the low-dimensional spaces, in which solutions are represented. In the proposed approach, low-dimensional spaces, constructed for the healthy structure, are enriched with appropriately selected columns of the flexibility matrix of the system. It can be shown that these spaces contain the solution to the original problem for the static case, while their dimension is much smaller. In order to allow their online construction for arbitrary localised features, the full flexibility matrix of the system should be available. To this end, a hierarchical representation is used for the matrices involved, allowing to compute the flexibility matrix efficiently and with reduced memory requirements. The resulting method offers significant speedups, without sacrificing the flexibility and accuracy of the full order model. The performance and limitations of the approach are studied through a series of examples in structural dynamics.

fracture mechanics↗

Quantifying Charge Carrier Localization in PBTTT Using Thermoelectric and Spectroscopic Techniques

Chemically doped poly[2,5-bis(3-alkylthiophen-2-yl)thieno[3,2- b ]thiophene] (PBTTT) shows promise for many organic electronic applications, but rationalizing its charge transport properties is challenging because conjugated polymers are inhomogeneous, with convoluted optical and solid-state transport properties. Herein, we use the semilocalized transport (SLoT) model to quantify how the charge transport properties of PBTTT change as a function of iron(III) chloride (FeCl 3 ) doping level. We use the SLoT model to calculate fundamental transport parameters, including the carrier density needed for metal-like electrical conductivities and the position of the Fermi energy level with respect to the transport edge. We then contextualize these parameters with other polymer-dopant systems and previous PBTTT reports. Additionally, we use grazing incidence wide-angle X-ray scattering and spectroscopic ellipsometry techniques to better characterize inhomogeneity in PBTTT. Our analyses indicate that PBTTT obtains high electrical conductivities due to its quickly rising reduced Fermi energy level, and this rise is afforded by its locally high carrier densities in highly ordered microdomains. Ultimately, this report sets a benchmark for comparing transport properties across polymer-dopant-processing systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A rotor-based multileaf collimator for beam shaping

We introduce a new style of multileaf collimator which employs rotors with angularly dependent radius to control the masking aperture: a rotor-based multileaf collimator (RMLC). Using a padlock-inspired mechanism, a single motor can set dozens of rotors, i.e. leaves, independently. This is especially important for an ultra-high vacuum (UHV) compatible MLC, since this reduces the number of actuators and vacuum feedthroughs required by more than an order of magnitude. This new RMLC will complement previous work employing a UHV compatible MLC with an emittance exchange beamline to create arbitrarily shaped beams on demand. A feed-forward control system which abstracts away the complexity of the RMLC operation, and is adaptable to real beamline conditions, is discussed and demonstrated in simulation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Toward real-time optimization through model reduction and model discrepancy sensitivities

Optimization problems arise in a range of scenarios, from optimal control to model parameter estimation. In many applications, such as the development of digital twins, it is essential to solve these optimization problems within wall-clock-time limitations. However, this is often unattainable for complex systems, such as those modeled by nonlinear partial differential equations. One strategy for mitigating this issue is to construct a reduced-order model (ROM) that enables more rapid optimization. In particular, the use of nonintrusive ROMs—those that do not require access to the full-order model at evaluation time—is popular because they facilitate the computation of optimization solutions within the wall-clock time requirements. However, the optimization solution will be unreliable if the iterates move outside the ROM training data. This article proposes the use of hyper-differential sensitivity analysis with respect to model discrepancy (HDSA-MD) as a computationally efficient tool to augment ROM-constrained optimization and improve its reliability. The proposed approach consists of two phases: (i) an offline phase where several full-order model evaluations are computed to train the ROM, and (ii) an online phase where a ROM-constrained optimization problem is solved, a limited number of full-order model evaluations are computed, and HDSA-MD is used to enhance the optimization solution. Numerical results are demonstrated for two examples, atmospheric contaminant control and wildfire ignition location estimation, in which a ROM is trained offline using inaccurate atmospheric data. In conclusion, the HDSA-MD update yields a significant improvement in the ROM-constrained optimization solution using only one full-order model evaluation online with corrected atmospheric data.

PDE-constrained optimization↗

The Localized Active Space Method with Unitary Selective Coupled Cluster

Here, we introduce a hybrid quantum-classical algorithm, the localized active space unitary selective coupled cluster singles and doubles (LAS-USCCSD) method. Derived from the localized active space unitary coupled cluster (LAS-UCCSD) method, LAS-USCCSD first performs a classical LASSCF calculation, then selectively identifies the most important parameters (cluster amplitudes used to build the multireference UCC ansatz) for restoring interfragment interaction energy using this reduced set of parameters with the variational quantum eigensolver method. We benchmark LAS-USCCSD against LAS-UCCSD by calculating the total energies of (H 2 ) 2 , (H 2 ) 4 , and trans-butadiene, and the magnetic coupling constant for a bimetallic compound [Cr 2 (OH) 3 (NH 3 ) 6 ] 3+ . For these systems, we find that LAS-USCCSD reduces the number of required parameters and thus the circuit depth by at least 1 order of magnitude, an aspect which is important for the practical implementation of multireference hybrid quantum-classical algorithms like LAS-UCCSD on near-term quantum computers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Néel-type antiferromagnetic order and magnetic field–temperature phase diagram in the spin-1/2 rare-earth honeycomb compound YbCl 3

Most of the searches for Kitaev materials deal with 4 d / 5 d magnets with spin-orbit-coupled J = 1 / 2 local moments such as iridates and α - RuCl 3 . Here we propose the monoclinic YbCl 3 with a Yb 3 + honeycomb lattice for the exploration of Kitaev physics. We perform thermodynamic, a c susceptibility, angle-dependent magnetic torque, and neutron diffraction measurements on YbCl 3 single crystal. We find that the Yb 3 + ion exhibits a Kramers doublet ground state that gives rise to an effective spin J eff = 1 / 2 local moment. Additionally, the compound exhibits short-range magnetic order below 1.20 K, followed by a long-range Néel-type antiferromagnetic order at 0.60 K, below which the ordered Yb 3 + spins lie in the a c plane with an angle of 16(11) ° away from the a axis. These orders can be suppressed by in-plane and out-of-plane magnetic fields at around 6 and 10 T, respectively. Moreover, the Néel temperature varies nonmonotonically under the out-of-plane magnetic fields, suggesting a reduced spin dimensionality. Finally, together with the strong in-plane magnetic anisotropy and the reduced order moment 0.8(1) μ B at 0.25 K, all indicate that YbCl 3 could be a two-dimensional spin system to proximate the Kitaev physics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Effects of dimensionality on the electronic structure of Ruddlesden-Popper chromates Sr n+1 ⁢Cr n ⁢O 3⁢n+1

Transition-metal oxides host a wide variety of electronic phenomena that can be significantly influenced by the effective dimensionality of the system under consideration. These include charge, spin, and orbital orderings, as well as unconventional superconductivity. In this context, the Ruddlesden-Popper chromates Sr n+1 Cr n O 3n+1 emerge as a particularly intriguing series of materials. Formally, the chromium atom displays a rather special 4+ oxidation state throughout the entire series. However, the effective dimensionality changes from quasi-2D to 3D as n increases from 1 to ∞. As a result, the insulating antiferromagnetic behavior observed for n = 1, 2, 3 transforms into itinerant antiferromagnetism with reduced transition temperature for the n = ∞ end member of the series, i.e., the perovskite SrCrO 3 . Further, distinct orbital orderings with exotic singlet states have been predicted for these systems. However, the lack of single-crystal bulk or thin-film samples has made experimental progress difficult. Here we demonstrate the synthesis of thin films of the perovskite SrCrO 3 and the associated layered chromates via oxide molecular beam epitaxy for n = 1 to n = 5. Our electrical transport measurements reveal a gradual evolution from a strongly insulating state in Sr 2 CrO 4 to a metallic state in the end member SrCrO 3 . X-ray absorption spectroscopy measurements demonstrate a varying hybridization strength of the Cr 4+ valence electrons across the series, helping to explain the trend in conduction. Finally, density functional theory calculations further confirm the observed transport trend and identify additional distortions present in the system.

36 MATERIALS SCIENCE↗

Advanced Models and Fast Simulation Algorithms for MARS Plants

With the successful commercialization of modular multilevel converters (MMC) and rapid increase in photovoltaic (PV) penetration and energy storage systems (ESS), one of the research areas being explored dynamically is integration of PV systems and ESS systems to MMCs. On this front, an integrated development of PV plants, ESS, and alternating- and direct-current (ac/dc) systems through Multi-Port Autonomous Reconfigurable Solar (MARS) plant is proposed and being researched upon. However, simulation of MARS systems on traditional existing software takes very long time (order of days) due to the presence of a large number of states and non-linear devices. Therefore, development of fast high fidelity models is of paramount importance to test the hardware design, control design and planning of power system expansions. In this context, two research objectives are addressed in this paper: firstly, an ultra fast single- or multi- CPU simulation algorithm to simulate the MARS system based on state-space models, hybrid discretization algorithm with a relaxation technique that reduces the imposed computational burden and a multi-rate method is presented. Secondly, a reduced order model of MARS is also presented. Both the developed techniques are validated with respect to reference PSCAD/EMTDC model.

Marthi, Phani Ratna Vanamali↗

Exploration of multifidelity UQ sampling strategies for computer network applications.

Network modeling is a powerful tool to enable rapid analysis of complex systems that can be challenging to study directly using physical testing. Here, two approaches are considered: emulation and simulation. The former runs real software on virtualized hardware, while the latter mimics the behavior of network components and their interactions in software. Although emulation provides an accurate representation of physical networks, this approach alone cannot guarantee the characterization of the system under realistic operative conditions. Operative conditions for physical networks are often characterized by intrinsic variability (payload size, packet latency, etc.) or a lack of precise knowledge regarding the network configuration (bandwidth, delays, etc.); therefore uncertainty quantification (UQ) strategies should be also employed. UQ strategies require multiple evaluations of the system with a number of evaluation instances that roughly increases with the problem dimensionality, i.e., the number of uncertain parameters. It follows that a typical UQ workflow for network modeling based on emulation can easily become unattainable due to its prohibitive computational cost. In this paper, a multifidelity sampling approach is discussed and applied to network modeling problems. The main idea is to optimally fuse information coming from simulations, which are a low-fidelity version of the emulation problem of interest, in order to decrease the estimator variance. By reducing the estimator variance in a sampling approach it is usually possible to obtain more reliable statistics and therefore a more reliable system characterization. Several network problems of increasing difficulty are presented. For each of them, the performance of the multifidelity estimator is compared with respect to the single fidelity counterpart, namely, Monte Carlo sampling. For all the test problems studied in this work, the multifidelity estimator demonstrated an increased efficiency with respect to MC.

97 MATHEMATICS AND COMPUTING↗

Hollow fiber mid-ID spectrometer with UV laser ablation sampling for fine spatial reoslution of isotope ratios in solids

We describe a system that combines isotope ratio analysis via mid-infrared (Mid-IR) laser absorption spectroscopy with fine spatial resolution sampling using a UV pulsed laser. The UV laser ablates a pit in a solid on the order of 10 microns in diameter. The sample-derived particulates resulting from laser ablation pass through a micro-combustor, and the resulting gas is analyzed using Mid-IR laser absorption spectroscopy in a capillary absorption spectrometer (CAS). The CAS uses a hollow fiber optic waveguide with a reflective inner coating and a small internal volume on the order of 1ml. The hollow fiber both guides the laser light from source to detector and contains the gas sample at reduced pressure. Near unity overlap between the laser beam and sample enables sensitive analysis with ultra-small sample size. A prototype system has been demonstrated to enable stable carbon isotopic analysis (d13C) with 1 per mil precision using <1 picomole of carbon and is currently being used to study nutrient exchange in soil/root/microbial rhizosphere studies. The smaller sample size of this system is enabling fine spatial resolution analysis (on the order of 10 microns), which is roughly an order of magnitude smaller than was possible with an isotope ratio mass spectrometer (IRMS). In addition to rhizosphere studies, the system can provide a useful tool for fine scale isotope analysis with applications in biological, forensic, and environmental science.

capillary absorption spectroscopy, laser ablation,↗