Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reduced order modeing”

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 235 records · Page 13

Simulations of ion cyclotron anisotropy instabilities in the terrestrial magnetosheath

Enhanced transverse magnetic fluctuations observed below the proton cyclotron frequency in the terrestrial magnetosheath have been identified as due to the proton cyclotron and helium cyclotron instabilities driven by the T-perpendicular greater than T-parallel condition of the sheath ions. One-dimensional hybrid computer simulations are used here to examine the nonlinear properties of these two growing modes at relatively weak fluctuation energies and for wave vectors parallel to the background magnetic field. Second-order theory predicts fluctuating magnetic field energies at saturation of the proton cyclotron anisotropy instability in semiquantitative agreement with the simulation results. Introduction of the helium component enhances the wave-particle exchange rate for proton anisotropy reduction by that instability, thereby reducing the saturation energy of that mode. The simulations demonstrate that wave-particle interactions by the proton cyclotron and helium cyclotron instabilities lead to the anticorrelation observed by Anderson and Fuselier (1993).

Gary, S. P.↗

Mitigating Sodium Ordering for Enhanced Solid Solution Behavior in Layered NaNiO 2 Cathodes

Abstract The O‐type layered nickel oxides suffer from undesired cooperative Jahn–Teller distortion stemming from Ni 3+ ions and undergo multiple biphasic structural transformations during the insertion/extraction of large Na + ions, posing a significant challenge to stabilize the structural integrity. We present here a systematic investigation of the impact of substituting 5 % divalent (Mg 2+ ) or trivalent (Al 3+ or Co 3+ ) ions for Ni 3+ to alleviate Na + ion ordering and perturb the Jahn–Teller effect to enhance structural stability. We gauge a fundamental understanding of the Mg−O and Na−O or Mg−O−Na bonding interactions, noting that the ionicity of the Mg−O bond deshields the electronic cloud of oxygen from Na + ions. Furthermore, calculations of the Van Vleck distortion modes reveal a relaxation of NiO 6 octahedra from Jahn–Teller distortion and a reduced electron density at the interlayer with Mg 2+ substitution. Long‐range ( operando X‐ray diffraction) and short‐range (magic angle spinning nuclear magnetic resonance) structural analyses provide insights into reduced ordering, allowing a stable continuous solid solution. Overall, Mg‐substitution results in a high‐capacity retention of ~96 % even after 100 cycles, showcasing the potential of this strategy for overcoming the structural instabilities and enhancing the performance of sodium‐ion batteries.

Sada, Krishnakanth↗

Simulations of Multi-Mode Combustion Regimes Realizable In a Gasoline Direct Injection Engine

Lean and dilute gasoline compression ignition (GCI) operation in spark ignition (SI) engines are an attractive strategy to attain high fuel efficiency and low NOx levels. However, this combustion mode is often limited to low-load engine conditions due to the challenges associated with autoignition controllability. In order to overcome this constrain, multi-mode (MM) operating strategies, consisting of advanced compression ignition (ACI) at low load and conventional SI at high load, have been proposed. In this three-dimensional computational fluid dynamics study, the concept of multi-mode combustion using two RON98 gasoline fuel blends (Co-Optima Alkylate and E30) in a gasoline direct injection (GDI) engine were explored. To this end, a new reduced mechanism for simulating the kinetics of E30 fuel blend is introduced in this study. To cover the varying engine load demands for multi-mode engines, primary combustion dynamics observed in ACI and SI combustion modes was characterized and validated against experimental measurements. In order to implement part-load conditions, a strategy of mode transition between SI and ACI combustion (i.e., mixed-mode combustion) was then explored numerically by creating a virtual test condition. The results obtained from the mixed-mode simulations highlight an important feature that deflagrative flame propagation regime coexists with ignition-assisted end-gas autoignition. This study also identifies a role of turbulent flow property adjacent to premixed flame front in characterizing the mixed-mode combustion. The employed hybrid combustion model was verified to perform simulations aiming at suitable range of multi-mode engine operations.

Kim, Sayop↗

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.↗

Time Resolved Irradiance of an Integrating Sphere Illuminated by a Mode-Locked Optical Parametric Oscillator

Tunable laser based illumination sources have been adopted by a number of recent hyperspectral instruments for pre-launch spectral radiometric calibration. Some of these sources use mode-locked cavity configurations, which produce a pulse train and create an undesirable fluctuation in the radiometric source at the mode locking frequency. Although typically this frequency is several orders of magnitude higher than the instrument response, care must be taken to ensure the resultant calibration is not biased due to detector or electronic saturation effects. The temporal averaging properties of an integrating sphere can be used to reduce high frequency fluctuations to an acceptable level. In this work, the time resolved output of an integrating sphere coupled with a mode-locked source was measured. A 76 cm diameter sphere with 30 cm output aperture was illuminated using a mode-locked pulse train at 76 MHz, with single pulse durations of 12 ps. The time constant of the sphere was found to be 36 ns, and the resultant temporal averaging of pulses produced a maximum time varying irradiance at the output of 20% of the mean. The form of this time variance is a sawtooth pattern, with a sharp rise during the introduction of each pulse into the sphere followed by a slower decay between pulses. A comparison of instrument calibration data generated with this integrating sphere using both a mode locked source and a continuous wave source is also given. .

McAndrew, Brendan↗

Finite-temperature many-body perturbation theory for anharmonic vibrations: Recursions, algebraic reduction, second-quantized reduction, diagrammatic rules, linked-diagram theorem, finite-temperature self-consistent field, and general-order algorithm

A unified theory is presented for finite-temperature many-body perturbation expansions of the anharmonic vibrational contributions to thermodynamic functions, i.e., the free energy, internal energy, and entropy. The theory is diagrammatically size-consistent at any order, as ensured by the linked-diagram theorem proved in this study, and, thus, applicable to molecular gases and solids on an equal footing. It is also a basis-set-free formalism, just like its underlying Bose–Einstein theory, capable of summing anharmonic effects over an infinite number of states analytically. It is formulated by the Rayleigh–Schrödinger-style recursions, generating sum-over-states formulas for the perturbation series, which unambiguously converges at the finite-temperature vibrational full-configuration-interaction limits. Two strategies are introduced to reduce these sum-over-states formulas into compact sum-over-modes analytical formulas. One is a purely algebraic method that factorizes each many-mode thermal average into a product of one-mode thermal averages, which are then evaluated by the thermal Born–Huang rules. Canonical forms of these rules are proposed, dramatically expediting the reduction process. The other is finite-temperature normal-ordered second quantization, which is fully developed in this study, including a proof of thermal Wick’s theorem and the derivation of a normal-ordered vibrational Hamiltonian at finite temperature. The latter naturally defines a finite-temperature extension of size-extensive vibrational self-consistent field theory. These reduced formulas can be represented graphically as Feynman diagrams with resolvent lines, which include anomalous and renormalization diagrams. Two order-by-order and one general-order algorithms of computing these perturbation corrections are implemented and applied up to the eighth order. The results show no signs of Kohn–Luttinger-type nonconvergence.

74 ATOMIC AND MOLECULAR PHYSICS↗

Projections of H-mode access and edge pedestal in the SPARC tokamak

In order to inform core performance projections and divertor design, the baseline SPARC tokamak plasma discharge is evaluated for its expected H-mode access, pedestal pressure and edge-localized mode (ELM) characteristics. A clear window for H-mode access is predicted for full field DT plasmas, with the available 25 MW of design auxiliary power. Additional alpha heating is likely needed for H-mode sustainment. Pressure pedestal predictions in the developed H-mode are surveyed using the EPED model. The projected SPARC pedestal would be limited dominantly by peeling modes and may achieve pressures in excess of 0.3 MPa at a density of approximately 3 × 10 20 m -3 . High pedestal pressure is partially enabled by strong equilibrium shaping, which has been increased as part of recent design iterations. Edge-localized modes (ELMs) with >1 MJ of energy are projected, and approaches for reducing the ELM size, and thus the peak energy fluence to divertor surfaces, are under consideration. The high pedestal predicted for SPARC provides ample margin to satisfy its high fusion gain (Q) mission, so that even if ELM mitigation techniques result in a 2x reduction of the pedestal pressure, Q > 2 is still predicted.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Spectrally accelerated edge and scrape-off layer gyrokinetic turbulence simulations

This paper presents the first gyrokinetic (GK) simulations of edge and scrape-off layer (SOL) turbulence accelerated by a velocity-space spectral approach in the full-f GK code GENE-X. Building upon the original grid velocity-space discretization, we derive and implement a new spectral formulation and verify the numerical implementation using the method of manufactured solution. We conduct a series of spectral turbulence simulations focusing on the TCV-X21 reference case (Oliveira et al., 2022 [26]) and compare these results with previously validated grid simulations (Ulbl et al., 2023 [25]). The spectral approach reproduces the outboard midplane (OMP) profiles (density, temperature, and radial electric field), dominated by trapped electron mode (TEM) turbulence, with excellent agreement and significantly lower velocity-space resolution. As a consequence, the spectral approach reduces the computational cost (CPUh) by at least an order of magnitude, of approximately 50 for the TCV-X21 case. This enables high-fidelity GK simulations to be performed within a few days on modern CPU-based supercomputers for medium-sized devices and establishes GENE-X as a powerful tool for studying edge and SOL turbulence, moving towards reactor-relevant devices like ITER.

Gyrokinetic↗

Ferroelastic Domain-Induced Electronic Modulation in Halide Perovskites

Lead halide perovskites have emerged as promising materials for optoelectronic applications due to their exceptional properties. In the all-inorganic CsPbBr 3 perovskites, ferroelastic domains formed during phase transitions enhance bulk transport and emissive efficiency. However, the microscopic mechanisms governing carrier dynamics remain poorly understood. In this study, we employ cathodoluminescence (CL) and micro-Raman spectroscopy to image and investigate the electronic properties of the ferroelastic domain walls in CsPbBr 3 single crystals. CL measurements reveal a reduced emissive yield and a slight redshift in emission at the domain walls. Further, micro-Raman studies provide spatially resolved mapping of vibrational modes, exhibiting second-order phonon modes localized at the domain boundaries. In conclusion, our findings suggest that electron–phonon coupling at the twin domain walls plays a critical role in facilitating efficient charge separation, thereby improving the optoelectronic performance of the CsPbBr 3 perovskites.

Carrier dynamics↗

Nonlinear analysis of rotor-bearing systems using component mode synthesis

The method of component mode synthesis is developed to determine the forced response of nonlinear, multishaft, rotor-bearing systems. The formulation allows for simulation of system response due to blade loss, distributed unbalance, base shock, maneuver loads, and specified fixed frame forces. The motion of each rotating component of the system is described by superposing constraint modes associated with boundary coordinates and constrained precessional modes associated with internal coordinates. The precessional modes are truncated for each component and the reduced component equations are assembled with the nonlinear supports and interconnections to form a set of nonlinear system equations of reduced order. These equations are then numerically integrated to obtain the system response. A computer program, which is presently restricted to single shaft systems, has been written and results are presented for transient system response associated with blade loss dynamics with squeeze film dampers, and with interference rubs.

Nelson, H. D.↗

Measurement of thermoacoustic convection heat transfer phenomenon

In this paper the results of an experimental investigation of thermoacoustic convection (TAC) heat transfer phenomenon in both zero-gravity and gravity environment are presented and compared with pure conduction heat transfer. The numerical solutions of the governing equations obtained by others for TAC heat transfer phenomenon are also discussed. The experimental results show that for rapid heating rate at a boundary, the contribution of TAC heat transfer to a gas could be significantly (one order of magnitude) higher than heat transfer rate from pure conduction. The results also show significantly reduced transient time in heat transfer processes involving thermoacoustic convective heat transfer mode in both space and gravity environment.

Parang, M.↗

Variable-Depth Liner Evaluation Using Two NASA Flow Ducts

Four liners are investigated experimentally via tests in the NASA Langley Grazing Flow Impedance Tube. These include an axially-segmented liner and three liners that use reordering of the chambers. Chamber reordering is shown to have a strong effect on the axial sound pressure level profiles, but a limited effect on the overall attenuation. It is also shown that bent chambers can be used to reduce the liner depth with minimal effects on the attenuation. A numerical study is also conducted to explore the effects of a planar and three higher-order mode sources based on the NASA Langley Curved Duct Test Rig geometry. A four-segment liner is designed using the NASA Langley CDL code with a Python-based optimizer. Five additional liner designs, four with rearrangements of the first liner segments and one with a redistribution of the individual chambers, are evaluated for each of the four sources. The liner configuration affects the sound pressure level profile much more than the attenuation spectra for the planar and first two higher-order mode sources, but has a much larger effect on the SPL profiles and attenuation spectra for the last higher-order mode source. Overall, axially variable-depth liners offer the potential to provide improved fan noise reduction, regardless of whether the axially variable depths are achieved via a distributed array of chambers (depths vary from chamber to chamber) or a group of zones (groups of chambers for which the depth is constant).

Jones, M. G.↗

Combustion noise and combustion instabilities in propulsion systems

This paper is concerned with some aspects of non-linear behavior of unsteady motions in combustion chambers. The emphasis is on conditions under which organized oscillations having discrete frequencies may exist in the presence of random motions. In order to treat the two types of motions together, and particularly to investigate coupling between noise and combustion instabilities, the unsteady field is represented as a synthesis of acoustic modes having time-varying amplitudes. Each of the amplitudes are written as the sum of two parts, one associated with the random field and the remainder representing the organized oscillations. After spatial averaging, the general problem is reduced to solution of a set of second-order ordinary differential equations whose structure depends on the sorts of nonlinear processes accounted for. This formulation accommodates any physical process; in particular, terms are included to represent noise sources, although only limited modeling is discussed. Our results suggest that random sources of noise have only small effects on combustion instabilities and seem not to be a cause of unstable motions. However, the coupling between the two sorts of unsteady motions may be important as an essential process in a proposed scheme for noise control. It is now a familiar observation that many nonlinear deterministic systems are capable of exhibiting apparently random motions called 'chaos.' This is a particularly interesting possibility for systems which also executed non-deterministic random motions. In combustion chambers, a nonlinear deterministic system (acoustical motions) exists in the presence of noise produced by flow separation, turbulent motions, and energy released by combustion processes. The last part of the paper is directed to the matter of discovering whether or not chaotic motions exist in combustion systems. Analysis has not progressed sufficiently far to answer the question. We report here recent results of processing data taken in one combustor to determine the dimensions of any attractors in the motions. No evidence has been found for chaos in the strict sense, but the method seems to be an important means of investigating the nonlinear behavior of combustion systems.

F. E. C. Culick↗

The PI-Mode of Project Management

The PI-Mode is NASA's new approach to project management. It responds to the Agency's new policy to develop scientific missions that deliver the highest quality science for a fixed cost. It also attempts to provide more research opportunities by reducing project development times and increasing the number of launches per year. In order to accomplish this, the Principal Investigator is placed at the helm of the project with full responsibility over all aspects of the mission, including instrument and spacecraft development, as well as mission operations and data analysis. This paper intends to study the PI-Mode to determine the strengths and weaknesses of such a new project management technique. It also presents an analysis of its possible impact on the scientific community and its relations with industry, NASA, and other institutions.

Isaac, Dan↗

Hutchinson Trace Estimation for high-dimensional and high-order Physics-Informed Neural Networks

Physics-Informed Neural Networks (PINNs) have proven effective in solving partial differential equations (PDEs), especially when some data are available by seamlessly blending data and physics. However, extending PINNs to high-dimensional and even high-order PDEs encounters significant challenges due to the computational cost associated with automatic differentiation in the residual loss function calculation. Herein, we address the limitations of PINNs in handling high-dimensional and high-order PDEs by introducing the Hutchinson Trace Estimation (HTE) method. Starting with the second-order high-dimensional PDEs, which are ubiquitous in scientific computing, HTE is applied to transform the calculation of the entire Hessian matrix into a Hessian vector product (HVP). This approach not only alleviates the computational bottleneck via Taylor-mode automatic differentiation but also significantly reduces memory consumption from the Hessian matrix to an HVP’s scalar output. We further showcase HTE’s convergence to the original PINN loss and its unbiased behavior under specific conditions. Comparisons with the Stochastic Dimension Gradient Descent (SDGD) highlight the distinct advantages of HTE, particularly in scenarios with significant variability and variance among dimensions. We further extend the application of HTE to higher-order and higher-dimensional PDEs, specifically addressing the biharmonic equation. By employing tensor-vector products (TVP), HTE efficiently computes the colossal tensor associated with the fourth-order high-dimensional biharmonic equation, saving memory and enabling rapid computation. The effectiveness of HTE is illustrated through experimental setups, demonstrating comparable convergence rates with SDGD under memory and speed constraints. Additionally, HTE proves valuable in accelerating the Gradient-Enhanced PINN (gPINN) version as well as the Biharmonic equation. Overall, HTE opens up a new capability in scientific machine learning for tackling high-order and high-dimensional PDEs.

Curse of dimensionality↗

Mitigating Sodium Ordering for Enhanced Solid Solution Behavior in Layered NaNiO 2 Cathodes

The O-type layered nickel oxides suffer from undesired cooperative Jahn–Teller distortion stemming from Ni 3+ ions and undergo multiple biphasic structural transformations during the insertion/extraction of large Na + ions, posing a significant challenge to stabilize the structural integrity. We present here a systematic investigation of the impact of substituting 5 % divalent (Mg 2+ ) or trivalent (Al 3+ or Co 3+ ) ions for Ni 3+ to alleviate Na + ion ordering and perturb the Jahn–Teller effect to enhance structural stability. Here, we gauge a fundamental understanding of the Mg−O and Na−O or Mg−O−Na bonding interactions, noting that the ionicity of the Mg−O bond deshields the electronic cloud of oxygen from Na + ions. Furthermore, calculations of the Van Vleck distortion modes reveal a relaxation of NiO 6 octahedra from Jahn–Teller distortion and a reduced electron density at the interlayer with Mg 2+ substitution. Long-range (operando X-ray diffraction) and short-range (magic angle spinning nuclear magnetic resonance) structural analyses provide insights into reduced ordering, allowing a stable continuous solid solution. Overall, Mg-substitution results in a high-capacity retention of ~96 % even after 100 cycles, showcasing the potential of this strategy for overcoming the structural instabilities and enhancing the performance of sodium-ion batteries.

25 ENERGY STORAGE↗

The Construction of Reduced Order Models for the HYBRID Repository

Reduced order models (ROMs) present an opportunity to reduce the computational burden of IES analysis using HYBRID models without significantly sacrificing data accuracy. Other applications of ROMs are digital twins, stochastic optimization studies, model exporting, and protection of proprietary information. In this paper the capability of creating ROMs of the models within the HYBRID repository will be demonstrated. Specifically, a linear regression ROM of a gas turbine and a dynamic mode decomposition with control (DMDc) ROM of a lithium-ion battery are constructed.

97 MATHEMATICS AND COMPUTING↗

Flight control synthesis for flexible aircraft using Eigenspace assignment

The use of eigenspace assignment techniques to synthesize flight control systems for flexible aircraft is explored. Eigenspace assignment techniques are used to achieve a specified desired eigenspace, chosen to yield desirable system impulse residue magnitudes for selected system responses. Two of these are investigated. The first directly determines constant measurement feedback gains that will yield a close-loop system eigenspace close to a desired eigenspace. The second technique selects quadratic weighting matrices in a linear quadratic control synthesis that will asymptotically yield the close-loop achievable eigenspace. Finally, the possibility of using either of these techniques with state estimation is explored. Application of the methods to synthesize integrated flight-control and structural-mode-control laws for a large flexible aircraft is demonstrated and results discussed. Eigenspace selection criteria based on design goals are discussed, and for the study case it would appear that a desirable eigenspace can be obtained. In addition, the importance of state-space selection is noted along with problems with reduced-order measurement feedback. Since the full-state control laws may be implemented with dynamic compensation (state estimation), the use of reduced-order measurement feedback is less desirable. This is especially true since no change in the transient response from the pilot's input results if state estimation is used appropriately. The potential is also noted for high actuator bandwidth requirements if the linear quadratic synthesis approach is utilized. Even with the actuator pole location selected, a problem with unmodeled modes is noted due to high bandwidth. Some suggestions for future research include investigating how to choose an eigenspace that will achieve certain desired dynamics and stability robustness, determining how the choice of measurements effects synthesis results, and exploring how the phase relationships between desired eigenvector elements effects the synthesis results.

Davidson, J. B.↗