Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reduced Order Models”

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

An adaptive scalable fully implicit algorithm based on stabilized finite element for reduced visco-resistive MHD

The magnetohydrodynamics (MHD) equations are continuum models used in the study of a wide range of plasma physics systems, including the evolution of complex plasma dynamics in tokamak disruptions. However, efficient numerical solution methods for MHD are extremely challenging due to disparate time and length scales, strong hyperbolic phenomena, and nonlinearity. Additionally, therefore the development of scalable, implicit MHD algorithms and high-resolution adaptive mesh refinement strategies is of considerable importance. In this work, we develop a high-order stabilized finite-element algorithm for the reduced visco-resistive MHD equations based on the MFEM finite element library (mfem.org). The scheme is fully implicit, solved with the Jacobian-free Newton-Krylov (JFNK) method with a physics-based preconditioning strategy. Our preconditioning strategy is a generalization of the physics-based preconditioning methods in Chacón et al. (2002) to adaptive, stabilized finite elements. Algebraic multigrid methods are used to invert sub-block operators to achieve scalability. A parallel adaptive mesh refinement scheme with dynamic load-balancing is implemented to efficiently resolve the multi-scale spatial features of the system. Our implementation uses the MFEM framework, which provides arbitrary-order polynomials and flexible adaptive conforming and non-conforming meshes capabilities. Results demonstrate the accuracy, efficiency, and scalability of the implicit scheme in the presence of large scale disparity. The potential of the AMR approach is demonstrated on an island coalescence problem in the high Lundquist-number regime (≥ 10 7 ) with the successful resolution of plasmoid instabilities and thin current sheets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1977

Vortex interactions and decay in aircraft wakes

The dynamic interactions of aircraft wake vortices are investigated using both inviscid and viscous models. The phenomenon of vortex merging resulting in the rapid aging of a vortex wake is examined in detail. It is shown that the redistribution of vorticity from convection and diffusion during merging is a mechanism effective in reducing the hazard of a wake. Inviscid computations show that the merging phenomenon may be sensitive to small changes in spanwise load distribution and that the fuselage vortex shed from the wing-fuselage junction can play a significant role in promoting merging of wing-tip and flap vortices. Vortex-wake-merging computations using a second-order closure model of turbulent transport indicate that a low-hazard wake occurs when the generating aircraft flap and wing-tip vortices of the same strength and sign. This optimum is achieved when the flap vortex is located outboard approximately 40% of the distance to the tip vortex.

Bilanin, A. J.↗
DOE OSTI2025

Efforts to stabilize composite localization elements

Localization finite elements seek to provide a robust framework for modeling ductile failure. They utilize the same constitutive model as the bulk material through the introduction of a length scale in a specialized deformation gradient that regularizes displacement discontinuity. Similar to many other elements, localization elements exhibit locking and associated pressure oscillations under incompressible plastic flow, which is a critical issue when attempting to model pressure-driven damage evolution. These issues can be drastically improved through what are essentially reduced integration techniques for the Jacobian and pressure, but there seem to be pressure-related instabilities that persist and are specific to localization elements. This memo summarizes recent efforts to mitigate and understand this problem, mostly for the 12-node composite wedge localization element in particular. At this point, it remains unclear whether the pressure fields within any localization element can be sufficiently stabilized in order to properly model failures that include softening or fracture.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Scalable and compact magnetocaloric heat pump technology

Magnetocaloric heat pumping (MCHP) promises to be more efficient than traditional vapor compression while also eliminating the deleterious effects of gaseous refrigerants. While MCHP devices have shown the temperature spans and efficiencies needed for different heating and cooling applications, they struggle to become commercially viable due to their large size and mass, and resultant high cost. This paper evaluates a baseline MCHP device and explores methods to boost its system power density (SPD). The key components of the baseline system are the gadolinium packed-particle bed active magnetic regenerator (AMR) and a magnetic source composed of permanent magnets and high permeability magnetic steel. To enhance the SPD, the paper evaluates maximizing the AMR volume, opting for first-order magnetocaloric materials, optimizing the magnet and AMR geometry, and reducing the size of magnets and magnetic steel parts. At larger thermal powers, increasing the AMR diameter and the number of magnetic poles were evaluated. Using finite element models, solid models, and estimates of magnetocaloric material performance, thermal powers ranging from 37 W to 44 kW at a nominal 10 K temperature span were projected, and SPD was estimated to improve from 6 W/kg to 81 W/kg. Neglecting end effects, an upper limit of 114 W/g is estimated. Compared to SPD of off-the-shelf compressors with similar environment temperatures, MCHP power density using gadolinium is competitive up to roughly 200 W of cooling power. This is extended to 1 kW when using LaFeSi alloys and up to 3 kW in the limiting case. In conclusion, these results indicate that the performance and mass of MCHP can match that of compressors, which is a critical step toward cost-competitive magnetocaloric technology.

42 ENGINEERING↗
DOE OSTI2026

Alchemy: A Model-Based Approach for 2D to 3D Autonomous Nuclear System Design

Engineering design of nuclear power plant (NPP) piping and equipment systems frequently bypasses crucial 2D system planning, instead moving straight to 3D modeling. This often leads to designs that exceed building envelope constraints, forcing expensive and time-consuming redesigns. When 2D modeling is employed, it typically involves labor-intensive manual workflows that convert 2D drawings into 3D models, resulting in inefficiencies and errors across design iterations. These workflows further suffer from poor software interoperability and dependence on proprietary software ecosystems, thereby contributing to schedule delays and cost overruns. This paper presents Alchemy, an autonomous framework that transforms 2D system definitions into Industry Foundation Classes (IFC)-compliant 3D building information models (BIMs) for expediting nuclear facility design at the conceptual preliminary phase. Using a model-based approach, the framework treats the 2D system diagram as the central reference model employed to automatically generate all subsequent outputs, ensuring consistency between the system definition and the resulting physical design. A web-based interface enables engineers to define hierarchical system topologies including associated equipment, geometric properties, and connectivity requirements. A two-phase equipment layout optimization algorithm automatically computes collision-free spatial configurations within predefined building envelopes. An artificial intelligence (AI)-assisted pipe routing module then generates orthogonal, collision-free routing paths, allowing the user to select either an A* search-based method or an Ant Colony Optimization (ACO)-based method. All outputs are authored natively in IFC format, relying on open-source technologies and standardized formats in order to ensure extensibility and eliminate proprietary software dependencies. The proposed framework is validated on two representative pressurized-water reactor (PWR)-based case studies, for which it autonomously generates IFC-compliant 3D models in minutes, drastically reducing workflows that typically require hours of manual effort. The generated model demonstrates topologically correct equipment placement, physically plausible spatial relationships, and collision-free pipe routing consistent with known PWR loop configurations. This work represents a foundational step toward digital engineering for nuclear facility preliminary design, with future ongoing development targeting design code compliance and expanded system complexity.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Near and long-term perspectives on strategies to decarbonize China's heavy-duty trucks through 2050

China needs to drastically reduce carbon dioxide (CO 2 ) emissions from heavy-duty trucks (HDTs), a key emitter in the growing transport sector, in order to address energy security concerns and meet its climate targets. We address existing research gaps by modeling feasibility, applicability, and energy and emissions impacts of multiple decarbonization strategies at different points in time. China still relies heavily on coal power, so impacts of new HDT technologies depend on the timing of their introduction relative to progress toward non-fossil power. We use a bottom-up model to simulate HDT energy consumption and CO 2 emissions through 2050. Results show that beginning to deploy battery electric and fuel-cell HDTs as early as 2020 and 2035, respectively, could achieve significant and the largest CO 2 emissions reduction by 2050 with a decarbonized power sector. However, viable near-term strategies-improving efficiency and logistics, switching to liquefied natural gas-could halve HDTs' current diesel consumption and CO 2 emissions by 2050. Our results underscore the need for a mix of near- and long-term policy and technology options to decarbonize China's HDTs.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Rational surfaces, flows and radial structure in the TJ-II stellarator

Abstract In this work, we report on the results obtained by measuring several turbulent quantities well inside the plasma edge by means of a Langmuir probe during dynamical rotational transform scans in the TJ-II stellarator, while applying a radial electric field to the edge plasma using a biasing probe. By calculating the intermittence parameter from floating potential measurements, we are able to identify a major low order rational surface and hence relate the probe measurements to the local value of the rotational transform. Based on the former, we are able to show that the poloidal plasma velocity (and hence radial electric field) has a significant radial structure that is clearly related to the rotational transform profile and in particular the lowest order rational surfaces in the range studied. The poloidal velocity is also affected by the edge biasing. The particle flux Γ was also found to exhibit a radial pattern, as did the flow shear suppression term ω E × B , but the relation of the former to the low-order rational surfaces was less clear. We surmise that this lack of direct correspondence is due to an unknown term in the turbulence evolution equation: the instability growth rate, γ . We make use of a reduced Magnetohydrodynamic turbulence model to interpret the results. Overall, a picture is obtained in which the plasma self-organizes towards a state with a clear radial pattern of the radial electric field, in line with expectations from some numerical studies describing the spontaneous formation of an ‘ E × B staircase’, consisting of alternating layers with fast and slow radial transport. In this state, the radial profiles of various quantities (density, temperature, pressure) will not be smooth.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Extending the Gutzwiller approximation to intersite interactions

In this work, we develop an extension of the Gutzwiller approximation (GA) formalism that includes the effects of Coulomb interactions of arbitrary range (including density density, exchange, pair hopping, and Coulomb-assisted hopping terms). This formalism reduces to the ordinary GA formalism for the multiband Hubbard models in the presence of only local interactions. This is accomplished by combining the 1 / z expansion—where z is the coordination number, and only the leading-order terms contribute in the limit of infinite dimensions—with a P R † P R - I expansion, where P R is the Gutzwiller projector on the site R . Furthermore, the method is conveniently formulated in terms of a Gutzwiller Lagrange function. We apply our theory to the extended single-band Hubbard model. Similarly to the usual Brinkman-Rice mechanism, we find a Mott transition. A valence skipping transition is observed, where the occupation of the empty and doubly occupied states for the Gutzwiller wave function is enhanced with respect to the uncorrelated Slater determinant wave function.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Melting of spatially modulated phases at domain wall/surface junctions in antiferrodistortive multiferroics

A physical understanding of the nature of spatially modulated phases (SMPs) in rare-earth-doped antiferrodistortive (AFD) multiferroics and how they behave close to surfaces and interfaces is lacking. Here the emergence of the antiferroelectric (AFE), ferroelectric (FE), or ferrielectric (AFE-FE) spatial modulation in the vicinity of the morphotropic phase transition in La x Bi 1 – x Fe O 3 ( x ~ 0.2 ) is explored on the atomic level using high-resolution scanning transmission electron microscopy (HRSTEM). The suppression, or “melting,” of the AFE-type SMP in the vicinity of the AFD twin wall/surface junction is revealed by HRSTEM in La 0.22 Bi 0.78 Fe O 3 films and explained by the hybrid approach combining Landau-Ginzburg-Devonshire (LGD) phenomenology and the semimicroscopic four-sublattice model (FSM). The LGD-FSM approach reduces the problem of AFE (or AFE-FE) SMP emergence and stability to the thermodynamic analysis of the free-energy functional with AFE, FE, and AFD long-range order parameters and two master parameters: the FE-AFE coupling strength between four neighboring A sites and the nonstoichiometry factor, which are proportional to the variations of La concentration in La x Bi 1 – x Fe O 3 films. Finally, we establish that the surface-induced melting of SMPs and the associated broadening of AFE AFD domain walls minimize the film free energy under certain conditions imposed on the master parameters and gradient energy below the critical value. The observed behavior provides insight into the origin of SMPs in AFD multiferroics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Deep Reinforcement Learning based Model-free On-line Dynamic Multi-Microgrid Formation to Enhance Resilience

Multi-microgrid formation (MMGF) is a promising solution for enhancing power system resilience. This paper proposes a new deep reinforcement learning (RL) based model-free on-line dynamic MMGF scheme. Additionally, the dynamic MMGF problem is formulated as a Markov decision process, and a complete deep RL framework is specially designed for the topologytransformable micro-grids. In order to reduce the large action space caused by flexible switch operations, a topology transformation method is proposed and an action-decoupling Q-value is applied. Then, a convolutional neural network (CNN) based multi-buffer double deep Q-network (CM-DDQN) is developed to further improve the learning ability of the original DQN method. The proposed deep RL method provides real-time computing to support the on-line dynamic MMGF scheme, and the scheme handles a long-term resilience enhancement problem using an adaptive on-line MMGF to defend changeable conditions. The effectiveness of the proposed method is validated using a 7-bus system and the IEEE 123-bus system. The results show strong learning ability, timely response for varying system conditions and convincing resilience enhancement.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Auxiliary Space Preconditioning of Finite Element Equations Using a Nonconforming Interior Penalty Reformulation and Static Condensation

We modify the well-known interior penalty finite element discretization method so that it allows for element-by-element assembly. This is possible due to the introduction of additional unknowns associated with the interfaces between neighboring elements. The resulting bilinear form, and a Schur complement (reduced) version of it, are utilized in a number of auxiliary space preconditioners for the original conforming finite element discretization problem. Furthermore, these preconditioners are analyzed on the fine scale, and their performance is illustrated on model second order scalar elliptic problems discretized with high order elements.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Collisional Excitation and Non-LTE Modeling of Interstellar Chiral Propylene Oxide

Abstract The first set of theoretical rotational cross sections for propylene oxide (CH 3 CHCH 2 O) colliding with cold He atoms has been obtained at the full quantum level using a high-accuracy potential energy surface. By scaling the collision reduced mass, rotational rate coefficients for collisions with para-H 2 are deduced in the temperature range 5–30 K. These collisional coefficients are combined with radiative data in a non-LTE radiative transfer model in order to reproduce observations of propylene oxide made toward the Sagittarius B2(N) molecular cloud with the Green Bank and Parkes radio telescopes. The three detected absorption lines are found to probe the cold (∼10 K) and translucent ( n H ∼ 2000 cm −3 ) gas in the outer edges of the extended Sgr B2(N) envelope. The derived column density for propylene oxide is N tot ∼ 3 × 10 12 cm −2 , corresponding to a fractional abundance relative to total hydrogen of ∼2.5 × 10 −11 . The present results are expected to help our understanding of the chemistry of propylene oxide, including a potential enantiomeric excess, in the cold interstellar medium.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1982

On the modelling and simulation of the dynamics and control of large flexible orbiting systems

This paper attempts to review the steps involved in the development of mathematical models that can be used to simulate the in-orbit dynamics of large flexible systems. The use of graph theoretic techniques can often be used to reduce the computational effort involved for calculating the eigenvalues of large ordered systems. Computer generated graphical techniques may provide additional insight into the understanding of elastic modal shape functions of complex systems. Finally the numerical techniques commonly used to develop shape and attitude control laws will be briefly reviewed.

Bainum, P. M.↗
NASA NTRS1982

Pole placement and order reduction in two-time-scale control systems through Riccati iteration

A transformation of variables taken from singular perturbations may be applied to two-time-scale linear systems in state space form to reduce the system to block-diagonal form with slow and fast modes decoupled. The transformation is easily computed by applying the new Riccati iteration. The iteration yields a solution to the nonsymmetric algebraic Riccati equation obtained by partitioning the original system matrix A. The numerical procedure is initiated with the trivial iterate L(0) = 0, and is globally convergent to the desired unique time scale decoupling solution. After transformation, the decoupled system may be used in controller design to achieve exact closed loop pole placement in the slow subsystem without altering the poles of the fast subsystem. The decoupled form may also be used to reduce system order by wetting a small parameter to zero. Provided the fast subsystem is stable, the order reduction can be expected to yield a good approximation to the original system. These methods are demonstrated using the 16th order linear model of a turbofan engine.

Anderson, L. R.↗
NASA NTRS1983

Reflectance Measurements, Part 2

Improving the productivity of spectroreflectometer readings taken to monitor induced contamination of lenses, mirrors, and coatings in the space environment are described. More sophisticated modelling techniques were applied to the laboratory data collected last summer in order to produce a more accurate model of the nonlinearities and noise sources of the instrumentation. This would ultimately allow for correction of known anomalies by the digital computer, yielding an anticipated order to magnitude improvement in the accuracy of the reflectance measurements, while reducing the total measuring time. The construction of and programming for a modernized version of the electronics used in Beckman DK-2 spectroreflectometers suffered because of a schedule lag in fabrication of a modified form of the carrier frequency pickoff. Separate development of the electronics and the computer programming was partially accomplished.

Brown, R. A.↗
NASA NTRS1987

Fixed order compensator design based on frequency shaped cost functionals

This paper treats the design of fixed order compensators using frequency shaped cost functionals. The formulation is done in an output feedback setting which exploits an observer canonical form to represent the compensator dynamics. The major advantages of the design approach are that the order of the compensator is fixed by the design process, and it is not necessary to realize the frequency shaping dynamics as an integral part of the compensator design. This avoids the two-step approach of full state feedback design followed by the design of an observer, which in the case of a full order observer requires a compensator of order n + nf where n is the order of the plant and nf is the order of the frequency shaping realization. The formulation also precludes the use of direct feedback of the plant output, thus reducing the effect of sensor noise and improving the robustness to high frequency unmodeled dynamics. A structural vibration model describing the fast dynamics of a lightweight flexible arm is used as an example to illustrate the design procedure.

Calise, Anthony J.↗
NASA NTRS1991

Analysis of opposed-jet hydrogen-air counter flow diffusion flame

An opposed-jet counterflow diffusion-flame configuration is considered for the analysis of a nitrogen-diluted hydrogen-air diffusion flame. A boundary-layer similarity solution is employed in order to reduce the governing equations to a set of equations in one independent variable. The equation set is written in the time-dependent form and solved by the finite-volume time-marching technique. This model uses detailed chemistry and accounts for the variations of Prandtl number and Lewis number as well as the effect of thermal diffusion on the flame. It is noted that a one-step model can predict several features of the flame, while the detailed-chemistry model can be used for fine-tuning the results. The present results indicate that thermal diffusion has negligible effect on the characteristics of the flame.

Ho, Y. H.↗
NASA NTRS2002

Computational Study of Near-limit Propagation of Detonation in Hydrogen-air Mixtures

A computational investigation of the near-limit propagation of detonation in lean and rich hydrogen-air mixtures is presented. The calculations were carried out over an equivalence ratio range of 0.4 to 5.0, pressures ranging from 0.2 bar to 1.0 bar and ambient initial temperature. The computations involved solution of the one-dimensional Euler equations with detailed finite-rate chemistry. The numerical method is based on a second-order spatially accurate total-variation-diminishing (TVD) scheme, and a point implicit, first-order-accurate, time marching algorithm. The hydrogen-air combustion was modeled with a 9-species, 19-step reaction mechanism. A multi-level, dynamically adaptive grid was utilized in order to resolve the structure of the detonation. The results of the computations indicate that when hydrogen concentrations are reduced below certain levels, the detonation wave switches from a high-frequency, low amplitude oscillation mode to a low frequency mode exhibiting large fluctuations in the detonation wave speed; that is, a 'galloping' propagation mode is established.

Yungster, S.↗