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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Mapping in situ the assembly and dynamics in aqueous supramolecular polymers

Supramolecular polymers, bonded through directional non-covalent interactions, closely mimic dynamic behaviors of biological nanofibers. However, the complexity of assembly pathways makes it highly challenging to unravel the nature of supramolecular dynamics in aqueous environments. Here we introduce a precise combinatorial titration methodology to probe in situ the assembly of peptide amphiphiles (PAs). This approach reveals a binary assembly mechanism governed by equilibrium between spheroidal micelles and β-sheet polymers. Weakening hydrogen bonding shifts the equilibrium towards micelles and decreases the internal structural order of filamentous polymers, promoting supramolecular dynamics. Extending this methodology to two-component copolymerization systems, we find a surprising tendency to form blocky nanostructures with reduced internal phase separation as the mismatch in peptide sequence decreases. Interestingly, while well-mixed copolymers acquire different dynamics, mismatched ones retain the characteristic supramolecular motion of their homopolymer counterparts. These critical insights into supramolecular dynamics offer strategies to tailor the dynamic functions of supramolecular nanomaterials.

Molecular self-assembly↗
DOE OSTI2023

DT-HYDRO

The software solves the time dependent, one dimensional (1D) coupled mass and momentum balance equations governing the elastic flow of water through the penstock, turbine and draft tube in a hydroelectric facility using a high order finite volume based method. The numerical method is based on the Kurganov-Tadmor central method paired with the Monotonic Upstream-centered Scheme for Conservation Laws (MUSCL). This solution method accurately resolves the fast transient behavior of the flow, including water hammer. Additionally, the software estimates the full 3D flow field within the turbine chamber in real time, a feat that is made possible by leveraging pre-computed CFD results by utilizing a reduced order modeling method based on an efficient randomized singular value decomposition (SVD) driven proper orthogonal decomposition (POD) with POD-mode weight regression. The reduced order model of the 3D flow is directly coupled to the 1D elastic flow model so the entire flow field through the penstock and turbine system is resolved quickly and with high fidelity.

Gurecky, William [Oak Ridge National Laboratory (O↗
DOE OSTI2020

Enhancing Reliability of Classified Shipments Between the Nuclear Security Enterprise

The extent of this project, as related to by the Nuclear Security Enterprise (NSE), is related to the current shipping and logistical issues and optimization opportunities between several different NSE sites that are the basis of the enterprise. This current scope is part of every site’s, described as Design Agency (DA) and Production Agency (PA), normal operation to meet the needs and mission of the National Nuclear Security Administration (NNSA). The NSE has a broad mission scope encompassing several differing DA and PA sites to meet the NNSA’s mission. These sites require frequent interaction and shipments of material and items between one another in order to meet the mission needs of the enterprise. Currently, there is no methodology that integrates schedules, standards, orders, processes, and best practices to optimize transfers between the differing sites. A centralized location and standard shipping process will optimize and improve the systems behavior, reducing lag time and increasing throughput at sites today. Ultimately this project proposes a central logistical platform, tiered escalation process, and standardized shipping forms to reduce taxpayer costs to meet the NSE’s mission more efficiently. This increase in processing logistics and process accuracy of shipments will save the NSE, not only from a monetary value, but it will reduce the frustrations of individuals working these outdated and over processed systems to improve their quality of life.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2025

Pressure Gain, Stability, and Operability of Methane/Syngas Based RDEs Under Steady and Transient Conditions (Final Project Report)

The scope of this work addresses key issues associated with losses associated with the detonation wave and other processes internal to the RDE operation, as well as it develops modeling tools for the evaluation of these losses and exhaust emissions in RDEs. The main challenge in studying RDEs is that RDE performance is highly reliant on the specifics of the design so much so that simple/canonical systems alone cannot provide useful engineering information, but practical RDE designs are sufficiently complex and involve extreme operational environments that detailed access either experimentally (laser diagnostics, for instance) or computationally (direct numerical simulations) are as yet to become practical. To overcome this challenge, we have conducted a combined experimental/simulation/analytical study investigating key phenomena that control the characteristics of operation of RDEs. As a result, the study has developed tools and methods that can be used to evaluate performance and design approaches using reduced-physics models, with the assumptions validated using detailed simulations, and the model prediction tested using experimental observations. The specific objectives of the research were: (1) Develop and demonstrate a low-loss fully axial injection concept, taking advantage of stratification effects to alter the detonation structure and position the wave favorably within the combustor; (2) Obtain stability and operability characteristics of an RDE across operating conditions to aid in the development of operability and performance rules for the operations of other systems; and (3) Develop quantitative metrics for performance gain as well as quantitative description of the loss mechanisms through a combination of diagnostics development, reduced-order modeling, and detailed simulations. The work conducted here has made contribution on design of low-loss inlets that has broad application within the power generation industry for use with pressure gain combustion. The operability and stability of different designs, while focusing on axial air inlet designs, has been analyzed. The effect of nozzle and injection conditions was studied. Models and simulations of exhaust emissions, focusing on NOx emission has been developed and used to investigate how operation of the RDE affect NOx production using Lagrangian analysis of RDE simulations. This work has built on previous programs, with the goal of further understanding operation of RDEs and elevate the readiness of design consideration. In addition, a suite of diagnostic and modeling tools have been developed to obtain quantitative metrics on performance based on measurements, which can be readily transferred to other experimental configurations.

08 HYDROGEN↗
DOE OSTI2021

A holistic approach to the optimization of neutron beam transport at the LANSCE facility

An advanced neutron beam collimation system has been developed as part of an effort to optimize neutron beam transport at the Weapons Neutron Research (WNR) Facility, located within the Los Alamos Neutron Science Center (LANSCE). The goal of this work was to develop and demonstrate techniques to tailor neutron delivery in order to provide maximum available flux on sample with a specific beam profile, while simultaneously reducing unwanted background. A holistic approach was taken, upgrading the facility spallation target, facility metrology infrastructure, and flight path shutter insert to support the implementation of this advanced collimation system. Here, modern instruments and software were employed to conduct facility surveys, characterization of as-built geometry of critical components, 3D ray tracing and neutron transport calculations. Beam images and flux measurements were acquired after installation to evaluate the performance of the collimation system and to demonstrate consistent agreement with MCNP simulations which showed an 87% increase of flux on sample while providing suppression of background neutron impingement on the sample frame by 10 3 .

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Superconductivity and phase diagrams of CaK(Fe 1- x Mn x ) 4 As 4 single crystals

In order to study the effects of Mn substitution on the superconducting and magnetic ground state of CaKFe 4 As 4 (T c = 35 K), members of the CaK (Fe 1–x Mn x ) 4 As 4 series have been synthesized by high-temperature solution growth in single-crystalline form and characterized by elemental analysis, thermodynamic, and transport measurements. These measurements show that the superconducting transition temperature decreases monotonically and is finally suppressed below 1.8 K as x is increased from 0 to 0.036. For x values greater than 0.016, signatures of a magnetic transition can be detected in both thermodynamic and transport measurements in which kinklike features allow for the determination of the transition temperature T* that increases as Mn substitution increases. Here, a temperature-composition (T–x) phase diagram is constructed, revealing a half-dome of superconductivity with the magnetic transition temperature T* appearing near 26 K for x ~ 0.017 and rising slowly up to 33 K for x ~ 0.036. In addition to the creation of the T–x phase diagram for CaK(Fe 1–x Mn x ) 4 As 4 , specific-heat data are used to track the jump in specific heat at T c ; the CaK(Fe 1–x Mn x ) 4 As 4 data do not follow the scaling of ΔC p with $T$$^{3}_{c}$ as many of the other Fe-based superconducting systems do. These data suggest that, as magnetic pair breaking is present, the jump in C p for a given T c is reduced. Elastoresistivity coefficients 2m 66 and m 11 – m 12 as a function of temperature are also measured. 2m 66 and m 11 – m 12 are qualitatively similar to CaK(Fe 1–x Ni x ) 4 As 4 . This may indicate that the magnetic order in Mn-substituted system may be still the same as CaK(Fe 1–x Mn x ) 4 As 4 . Superconductivity of CaK(Fe 1–x Mn x ) 4 As 4 is also studied as a function of magnetic field. A clear change in H' c2 (T)/T c , where H' c2 (T) is dH' c2 (T)/dT, at x ~ 0.015 is observed and probably is related to change of the Fermi surface due to magnetic order. Coherence lengths and the London penetration depths are also calculated based on H c1 and H c2 data. Coherence lengths as the function of x also show the changes near x = 0.015, again consistent with Fermi-surface changes associated with the magnetic ordering seen for higher-x values.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

A Voltage Injection-Based Current Harmonics Suppression Strategy for Six-Phase PMSM With Nonsinusoidal Back EMF

This article presents a voltage injection-based current harmonic minimization method for wide bandgap device-based six-phase permanent magnet synchronous machine (PMSM) drives. The presence of low-order harmonics is common in the back EMF of PMSMs with low slot and high pole numbers. The back EMF harmonics generate low-order harmonic currents that increase torque ripples. Furthermore, zero-sequence current (ZSC) is also observed if all six phases are connected to a common neutral point. ZSC does not affect the torque ripple but reduces the system efficiency. A special technique is required to reduce ZSC for variable speed drives as direct ac component to synchronous reference frame (SRF) dc transformation is not possible for ZSC. In this article, a voltage injection-based current harmonics reduction method is proposed. The injected harmonic voltage coefficients are calculated from the harmonic estimations of back EMF and updated using closed-loop harmonic compensators to improve steady-state performance. For ZSC reduction, third harmonics power controller, and for 6n ± 1 (n = 1, 2,...) order harmonic current reduction, SRF current controllers are presented in this article. Furthermore, the current harmonics minimization method is experimentally verified using a GaN-based three-level active neutral point clamped inverter-driven six-phase PMSM.

42 ENGINEERING↗
DOE OSTI2024

Circular Economy for Photovoltaics in Service of Energy Transition

The challenge of energy transition is immediate and immense; current projections target 75 TW of photovoltaics (PV) capacity by 2050. While any transition to renewable energy technology is preferable to the current fossil-based system, it is ideal to improve the sustainability of PV to minimize negative environmental and social impacts. Circular economy (CE) has been proposed as a method to improve the sustainability of PV, especially for emerging materials like perovskites. CE is a set of actions, principles, and systems which aim to design out waste and keep products and materials in use, to reduce environmental impacts and enable sustainable development. At the most basic level, CE is "reduce, reuse, recycle", the R-actions, in ranked order. CE of a PV technology can be metricized in a variety of ways, such as the Material Circularity Indicator (Smith and Jones, Ellen MacArthur Foundation, 2019) or recycling rates. Unfortunately, standard CE metrics have several shortcomings for measuring renewable energy technologies in the context of deployment for energy transition (Figge 2018, Saidani 2019): 1) Only measure mass flows; 2) De-prioritization of the use phase in favor of mass circularity when scoring; and 3) Tight focus on a single product scale The use phase and energy flows of PV are key to energy transition, and therefore need to be quantified. Additionally, correlating product-scale to system-scale is necessary for quantifying the environmental impacts of energy transition. Life Cycle Assessment (LCA) can address some of these concerns, but also focuses on a single product scale and has trouble capturing the dynamics of system-scale energy transition, such as the interaction of module lifetime with manufacturing demands for energy transition deployment schedules. Therefore, we developed an open-source Python-based system dynamics model to quantify the mass, energy and carbon impacts of CE R-actions for PV technologies in the energy transition; PV in the CE (PV ICE) (Ovaitt & Mirletz 2021). The tool captures supply chains from material extraction through end of life, incorporating 5 circular end of life pathways. PV ICE takes in any evolving bill of materials, module properties and deployment schedule to support researchers and decision makers with data-backed insights. In this work, we quantify and compare proposed CE sustainable PV module designs and lifecycle management strategies, spanning currently commercialized technologies, government and industry technology targets, and several low Technology Readiness Level (TRL) emerging PV technologies, including perovskites. Our analyses capture the projected evolutions of lifetime, efficiency and material circularity of these PV technologies, as well as their material supply chains. Our analyses emphasize the importance of examining a suite of metrics to identify priorities and tradeoffs, and inform design or lifecycle management decisions holistically. Previous analyses have demonstrated the central importance of PV module lifetime to support energy transition while minimizing impacts. High levels of material circularity (>90%) enable minimizing lifecycle wastes, can reduce virgin material demands if paired with improving efficiency, but demonstrate tradeoffs in energy return on investment. In the fervor of new material and technology development, it is important to remember that CE is not the end goal; decarbonization and energy transition are the end goal. CE should be used in service to improve the sustainability of PV, and R-actions evaluated for their usefulness and efficacy to this end.

carbon↗
DOE OSTI2026

Time-dependent-bases with local CUR decomposition method for accelerating turbulent combustion simulations

Here, this study presents a novel reduced-order modeling framework, Time-Dependent Bases with Local CUR decomposition (TDB-L-CUR), designed to efficiently and accurately approximate the species transport equations in reacting flow simulations. The method extends the existing TDB-CUR approach for chemically reacting flows (Jung et al. Comput. Methods Appl. Mech. Engrg. 437 (2025) 117758), which leverages matrix decomposition techniques to form a global-in-space, time-dependent low-dimensional manifold. While TDB-CUR performs well in homogeneous systems, it may be less well-suited to spatially heterogeneous systems such as turbulent flames, where higher-rank approximations are typically required. The proposed TDB-L-CUR framework introduces two methodological extensions to the baseline approach. First, it applies unsupervised clustering to partition the physical domain into distinct regions, enabling spatially localized manifold construction, thereby reducing the rank required for the reduced-order representation. Second, it incorporates a computational singular perturbation (CSP)-based scheme for identifying and penalizing fast species, allowing for spatio-temporally adaptive mitigation of chemical stiffness. The proposed framework is validated on a hierarchy of test cases, including a one-dimensional premixed flame, a two-dimensional nonpremixed ignition case with vortex interaction, and a three-dimensional turbulent premixed flame. TDB-L-CUR significantly improves accuracy over TDB-CUR while further reducing computational cost. The fully on-the-fly formulation of TDB-L-CUR (i.e., requiring no offline training or prior knowledge) makes it a robust and scalable tool for reduced-order modeling of reactive flows.

Local manifold↗
DOE OSTI2023

Interface flux recovery framework for constructing partitioned heterogeneous time-integration methods

A common approach for the development of partitioned schemes employing different time integrators on different subdomains is to lag the coupling terms in time. This can lead to accuracy issues, especially in multistage methods. Here, in this article, we present a novel framework for partitioned heterogeneous time-integration methods, which allows the coupling of arbitrary multistage and multistep methods without reducing their order of accuracy. At the core of our approach are accurate estimates of the interface flux obtained from the Schur complement of an auxiliary monolithic system. We use these estimates to construct a polynomial-in-time approximation of the interface flux over the current time coupling window. This approximation provides the interface boundary conditions necessary to decouple the subdomain problems at any point within the coupling window. In so doing our framework enables a flexible choice of time-integrators for the individual subproblems without compromising the time-accuracy at the coupled problem level. This feature is the main distinction between our framework and other approaches. To demonstrate the framework, we construct a family of partitioned heterogeneous time-integration methods, combining multistage and multistep methods, for a simplified tracer transport component of the coupled air-sea system in Earth system models. We report numerical tests evaluating accuracy and flux conservation for different pairs of time-integrators from the explicit Runge-Kutta and Adams-Moulton families.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

A linearity preserving nodal variation limiting algorithm for continuous Galerkin discretization of ideal MHD equations

A stabilized continuous Galerkin (CG) method for magnetohydrodynamics (MHD) is presented herein. Ideal, compressible inviscid MHD equations are discretized in space on unstructured meshes using piecewise linear or bilinear finite element bases to get a semi-discrete scheme. Stabilization is then introduced to the semi-discrete method in a strategy that follows the algebraic flux correction paradigm. This involves adding some artificial diffusion to the high order, semi-discrete method and mass lumping in the time derivative term. The result is a low order method that provides local extremum diminishing properties for hyperbolic systems. The difference between the low order method and the high order method is scaled element-wise using a limiter and added to the low order scheme. The limiter is solution dependent and computed via an iterative linearity preserving nodal variation limiting strategy. The stabilization also involves an optional consistent background high order dissipation that reduces phase errors. The resulting stabilized scheme is a semi-discrete method that can be applied to inviscid shock MHD problems and may be even extended to resistive and viscous MHD problems. To satisfy the divergence free constraint of the MHD equations, we add parabolic divergence cleaning to the system. Various time integration methods can be used to discretize the scheme in time. We demonstrate the robustness of the scheme by solving several shock MHD problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

3D theory of microscopic instabilities driven by space-charge forces

Microscopic, or short-wavelength, instabilities are known for a drastic reduction of the beam quality and strong amplification of the noise in a beam. Space charge and coherent synchrotron radiation are known to be the leading causes of such instabilities. In this paper, we present a rigorous 3D theory of such instabilities driven by the space-charge forces. We define the condition when our theory is applicable to an arbitrary accelerator system with 3D coupling. Finally, we derive a linear integral equation describing such instability and identify conditions when it can be reduced to an ordinary second-order differential equation.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Neutrinoless double-β decay: Combining quantum Monte Carlo and the nuclear shell model with the generalized contact formalism

Neutrinoless double beta decay searches can determine the Majorana nature of neutrinos, the absolute neutrino mass, and provide invaluable insights on the matter dominance of the universe. However, the uncertainty in the nuclear matrix elements that govern the decay limits the physics reach of these experiments. Here, we devise a novel framework based on the generalized contact formalism that combines the nuclear shell model and quantum Monte Carlo methods and compute the neutrinoless double-beta decay of nuclei used in the most advanced experiments, including 76 Ge, 130 Te, and 136 Xe. Our results cover all relevant terms, including the leading-order short-range operator recognized recently. We validate our method in light nuclei by comparing against accurate variational Monte Carlo results. On heavy systems we obtain reduced nuclear matrix elements compared with previous calculations due to additional correlations captured by quantum Monte Carlo and introduced within the generalized contact formalism, suggesting longer decay half-lives than previously considered. On the other hand, we find an enhancement of the nuclear matrix elements due to the new short-range operator.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Reducing the complexity of finite-temperature auxiliary-field quantum Monte Carlo

The auxiliary-field quantum Monte Carlo (AFMC) method is a powerful and widely used technique for ground-state and finite-temperature simulations of quantum many-body systems. Here we introduce several algorithmic improvements for finite-temperature AFMC calculations of dilute fermionic systems that reduce the computational complexity of most parts of the algorithm. This is principally achieved by reducing the number of single-particle states that contribute at each configuration of the auxiliary fields to a number that is of the order of the number of fermions. Our methods are applicable for both the canonical and grand-canonical ensembles. We demonstrate the reduced computational complexity of the methods for the homogeneous unitary Fermi gas.

67.85.Lm↗
DOE OSTI2026

Reaching quantum critical point by adding nonmagnetic disorder in single crystals of superconductor (Ca 𝑥 ⁢Sr 1−𝑥 ) 3 ⁢Rh 4 ⁢Sn 13

The Remeika series superconductor, (Ca 𝑥 ⁢Sr 1−𝑥 ) 3 ⁢Rh 4 ⁢Sn 13 , shows a rare nonmagnetic quantum critical point (QCP) associated with the continuous charge-density wave (CDW) and structural transition under the “dome” of superconductivity achieved by tuning composition and applying pressure. Here, we use a nonmagnetic pointlike disorder induced by 2.5 MeV electron irradiation to suppress the CDW and drive the system to and even beyond the QCP. This conclusion is based on a clear evolution of temperature-dependent resistivity, 𝜌⁡(𝑇), from the Fermi liquid to the non-Fermi liquid regime with increasing amount of disorder. Starting on the CDW side, below the suggested QCP concentration of 𝑥 𝑐 =0.9, added disorder resulted in a progressively larger linear term and a reduced quadratic term in 𝜌⁡(𝑇). Nearly perfect 𝑇−linear dependence is observed at the dose at which long-range CDW order is suppressed to 𝑇=0, consistent with the expectations. We refine the QCP location in this system and place it in the interval between 𝑥=0.75 and 0.85. Our results strongly support the concept that the disorder can tune the system to the quantum critical regime and even beyond. It follows from the argument by Imry and Ma that any ordered phase is unstable toward quenched disorder. Introduced in a controlled way, this disorder becomes a nonthermal tuning parameter likely applicable to a variety of different systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗