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

Observation of Nonaxisymmetric Standard Magnetorotational Instability Induced by a Free-Shear Layer

The standard magnetorotational instability (SMRI) with a magnetic field component parallel to the rotation axis is widely believed to be responsible for the fast accretion in astronomical disks. In conventional base flows with a Keplerian profile or an ideal Couette profile, most studies focus on axisymmetric SMRI, since excitation of nonaxisymmetric SMRI in such flows requires a magnetic Reynolds number (Rm) more than an order of magnitude larger. Here, we report that, in a magnetized Taylor-Couette flow, nonaxisymmetric SMRI with an azimuthal mode number m = 1 can be triggered by a free-shear layer in the base flow at Rm ≳1, the same threshold as for axisymmetric SMRI. Global linear analysis reveals that the free-shear layer reduces the required Rm, possibly by introducing an extremum in the vorticity of the base flow. Nonlinear simulations validate the results from linear analysis and confirm that a novel instability recently discovered experimentally [Wang et al., Nat. Commun. 13, 4679 (2022)] is the nonaxisymmetric m =1 SMRI. Further, our finding has astronomical implications as free-shear layers are ubiquitous in celestial systems, such as the disk-star boundary layer, the solar tachocline, and the edge of planet-opened gaps in protoplanetary disks.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Entanglement and seniority

In quantum mechanics entanglement is the most striking phenomenon which has no counterpart in classical systems. Although, different approaches have already been developed to study correlations in the case of indistinguishable particles, the exploration of the so-called mode entanglement is still in its initial stage in nuclear physics. Study of mode-entanglement in the seniority model, derivation of analytic formulas for the one-body reduced density matrix of states with seniority $ν$ = 0, 1, 2, and $ν$ = 3, and also determination of the particle number dependence of the one-body reduced density matrix for arbitrary seniority. In addition, comparison of the predictions of the seniority model with the results of the shell model to gain insight into the structure and correlations of the ground and lowest yrast states. In the seniority model the analytic results are calculated using the quasi-spin formalism. The numerical shell model calculations are carried out in the standard shell model framework using an inert core. The applied realistic effective interactions are derived earlier using the $G$-matrix formalism. The density matrix renormalization group method is also applied in order to directly calculate the mode entropies. In the seniority model simple analytical expressions are given for the mode entropies. The peculiar behavior of the half-filled shells and the seniority-zero states are revealed. Numerical results are presented for the lightest stable calcium isotopes and for 94 Ru nucleus. For 94 Ru, the seniority model accounts for the 0$g$ 9/2 mode entropies, but seniority mixing is important for certain yrast states. Furthermore, interaction induced quantum fluctuations decrease the occupation of the 0$f$ 5/2 , 1$p$ 3/2 and 1$p$ 1/2 shells, and amount in non-negligible mode entropies on these shells, too, clearly outside the scope of the simple (0$g$ 9/2 ) 4 seniority model. The 0$f$ 7/2 shell-based seniority model is more accurate for the Ca isotopes, but seniority mixing is substantial for some 44 Ca yrast states, too. Mode and one-body entanglement entropies are useful tools to investigate the structure of quantum correlations in nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Electron–Phonon and Spin–Lattice Coupling in Atomically Thin Layers of MnBi 2 Te 4

MnBi 2 Te 4 represents a new class of magnetic topological insulators in which novel quantum phases emerge at temperatures higher than those found in magnetically doped thin films. Here in this paper, we investigate how couplings between electron, spin, and lattice are manifested in the phonon spectra of few-septuple-layer thick MnBi 2 Te 4 . After categorizing phonon modes by their symmetries, we study the systematic changes in frequency, line width, and line shape of a spectrally isolated A1g mode. The electron–phonon coupling increases in thinner flakes as manifested in a broader phonon line width, which is likely due to changes of the electron density of states. In 4- and 5-septuple thick samples, the onset of magnetic order below the Néel temperature is concurrent with a transition to an insulating state. We observe signatures of a reduced electron–phonon scattering across this transition as reflected in the reduced Fano parameter. Finally, spin–lattice coupling is measured and modeled from temperature-dependent phonon frequency.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Variational multiscale reinforcement learning for discovering reduced order closure models of nonlinear spatiotemporal transport systems

Abstract A central challenge in the computational modeling and simulation of a multitude of science applications is to achieve robust and accurate closures for their coarse-grained representations due to underlying highly nonlinear multiscale interactions. These closure models are common in many nonlinear spatiotemporal systems to account for losses due to reduced order representations, including many transport phenomena in fluids. Previous data-driven closure modeling efforts have mostly focused on supervised learning approaches using high fidelity simulation data. On the other hand, reinforcement learning (RL) is a powerful yet relatively uncharted method in spatiotemporally extended systems. In this study, we put forth a modular dynamic closure modeling and discovery framework to stabilize the Galerkin projection based reduced order models that may arise in many nonlinear spatiotemporal dynamical systems with quadratic nonlinearity. However, a key element in creating a robust RL agent is to introduce a feasible reward function, which can be constituted of any difference metrics between the RL model and high fidelity simulation data. First, we introduce a multi-modal RL to discover mode-dependant closure policies that utilize the high fidelity data in rewarding our RL agent. We then formulate a variational multiscale RL (VMRL) approach to discover closure models without requiring access to the high fidelity data in designing the reward function. Specifically, our chief innovation is to leverage variational multiscale formalism to quantify the difference between modal interactions in Galerkin systems. Our results in simulating the viscous Burgers equation indicate that the proposed VMRL method leads to robust and accurate closure parameterizations, and it may potentially be used to discover scale-aware closure models for complex dynamical systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Effect of parallel flow on resonant layer responses in high beta plasmas

Abstract Resonant layers in a tokamak respond to non-axisymmetric magnetic perturbations by amplifying the mode amplitude and balancing the plasma rotation through magnetic reconnection and force balance, respectively. This resonant response can be characterized by local layer parameters and especially by a single quantity in the linear regime, the so-called inner-layer Δ. The computation of Δ under two-fluid drift-MHD formalism has been progressed by reducing the order of the system in the phase space, where the shielding current is approximated as being only carried by electrons, a posteriori . In this study, we relax the approximation and compute Δ accounted for by the parallel flow associated with the ion shielding current. The posteriori is numerically verified in great agreement with the original SLAYER developed in a previous paper (J.-K. Park 2022 Phys. Plasmas 29 072506). Extending the resonant layer response theory to high β plasmas, our research findings answer two important questions: how the parallel flow influences the resonant layer response and why the parallel flow effect appears in high β plasmas. The complicated plasma compression in high β regime allows the parallel flow response to give rise to the ion shielding current, which not only shifts the zero-crossing condition of the ExB flow but also enhances the field penetration threshold. Technically, the Riccati matrix transformation method is adapted to handle the numerical stiffness due to the increased order of the system. The high fidelity of this numerical method makes use of further extension of the model to higher-order systems to take other physical phenomena into account. This work is envisaged to predict the resonant layer response under high β fusion reactor conditions.

Lee, Yeongsun (ORCID:000000034474416X)↗
DOE OSTI2021

A comparative study of Higgs boson production from vector-boson fusion

The data taken in Run II at the Large Hadron Collider have started to probe Higgs boson production at high transverse momentum. Future data will provide a large sample of events with boosted Higgs boson topologies, allowing for a detailed understanding of electroweak Higgs boson plus two-jet production, and in particular the vector-boson fusion mode (VBF). We perform a detailed comparison of precision calculations for Higgs boson production in this channel, with particular emphasis on large Higgs boson transverse momenta, and on the jet radius dependence of the cross section. We study fixed-order predictions at next-to-leading order and next-to-next-to-leading order QCD, and compare the results to NLO plus parton shower (NLOPS) matched calculations. The impact of the NNLO corrections on the central predictions is mild, with inclusive scale uncertainties of the order of a few percent, which can increase with the imposition of kinematic cuts. We find good agreement between the fixed-order and matched calculations in non-Sudakov regions, and the various NLOPS predictions also agree well in the Sudakov regime. We analyze backgrounds to VBF Higgs boson production stemming from associated production, and from gluon-gluon fusion. At high Higgs boson transverse momenta, the Δy$_{jj}$ and/or m$_{jj}$ cuts typically used to enhance the VBF signal over background lead to a reduced efficiency. We examine this effect as a function of the jet radius and using different definitions of the tagging jets. QCD radiative corrections increase for all Higgs production modes with increasing Higgs boson p$_{T}$, but the proportionately larger increase in the gluon fusion channel results in a decrease of the gluon-gluon fusion background to electroweak Higgs plus two jet production upon requiring exclusive two-jet topologies. We study this effect in detail and contrast in particular a central jet veto with a global jet multiplicity requirement.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Laser beam pointing stabilization using analog position-sensitive diodes

Laser wire scanners have been used for ion beam profile and emittance measurements at the Spallation Neutron Source linear accelerator. Due to propagation distances exceeding 100 m, reliable measurements require laser beam stabilization, previously accomplished with a feedback control loop using a digital camera as a position sensor [Hardin et al., Opt. Express 19, 2874 (2011)]. Here, this paper presents an upgraded pointing stabilization system utilizing analog position-sensitive diodes (PSDs) as the position detector, optimized with an optical diffuser. The new system significantly surpasses the previous system by offering an order-of-magnitude improvement in radiation tolerance and a bandwidth limited only by the burst repetition rate of the laser beam. Operating with a 60-Hz burst mode laser, the new system effectively suppresses laser beam drifts up to 30 Hz, compared to the previous system’s 4 Hz limit. In addition, the amplitude of spectral components below 0.1 Hz was reduced by a factor of more than 100, which is over seven times greater than was achieved by the previous system.

Data acquisition↗
DOE OSTI2020

New Septum Dump Quadrupole Magnet Design

The purpose of the septum dump quadrupole (SDQ) magnet is to expand the beam size on the septum dump of the DARHT Axis 2 downstream transport in the same manner as the S2 magnet expands the beam on the beam stop for mode 1 operations of DARHT Axis 2. This is necessary in order to minimize the heating of the carbon beam dump and suppress desorption of ions. Initially, ions are generated through stimulated desorption and ionization of the monolayers of water vapor stuck to the beam dump surface. These ions propagate upstream at high speeds and act as an electrostatic focusing lens by reducing the beam space charge. As the beam size on the dump decreases, the dump is sufficiently heated so that a surface plasma is formed and much greater quantities of ions are generated through space charge limited emission. The interaction of these ions with the beam causes significant disruption of the beam.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Solid-state laser refrigeration of core-shell polystyrene microspheres

Microlaser designs based on the coupling of whispering gallery modes (WGMs) with the upconversion processes which take place within lanthanide-doped nanoparticles (UCNPs) have been demonstrated and shown to have many valuable qualities, such as high Q factors and low lasing thresholds. One obstacle that these microlaser designs still face is the challenges caused by photothermal heating of the gain medium, which could be solved through the design of a radiation balanced microlaser. In this work, WGM microresonators composed of 5 µm diameter polystyrene spheres are fabricated with a layer of Yb3+-doped NaYF4 UCNPs in order to test if the anti-Stokes cooling properties of the UCNPs can cool the microresonator and its environment under laser irradiation. We find via calibrated mean fluorescence spectroscopy that the UCNPs can cool their local environment by as much as 23 °C and significantly reduce the heating of the aqueous environment surrounding the microresonator, showing promise for inclusion in a design for a radiation balanced microlaser.

whispering gallery modes, optical refrigeration, m↗
DOE OSTI2025

On reduced modelling of the modulational dynamics in magnetohydrodynamics

This paper explores structure formation in two-dimensional magnetohydrodynamic (MHD) turbulence as a modulational instability (MI) of turbulent fluctuations. We focus on the early stages of structure formation and consider simple backgrounds that allow for a tractable model of the MI while retaining the full chain of modulational harmonics. This approach allows us to systematically examine the validity of popular closures such as the quasilinear approximation and other low-order truncations. We find that, although such simple closures can provide quantitatively accurate approximations of the MI growth rates in some regimes, they can fail to capture the modulational dynamics in adjacent regimes even qualitatively, falsely predicting MI when the system is actually stable. We find that this discrepancy is due to the excitation of propagating spectral waves (PSWs) which can ballistically transport energy along the modulational spectrum, unimpeded until dissipative scales, thereby breaking the feedback loops that would otherwise sustain MIs. The PSWs can be self-maintained as global modes with real frequencies and drain energy from the primary structure at a constant rate until the primary structure is depleted. To describe these waves within a reduced model, we propose an approximate spectral closure that captures them and MIs on the same footing. We also find that introducing corrections to ideal MHD, conservative or dissipative, can suppress PSWs and reinstate the accuracy of the quasilinear approximation. In this sense, ideal MHD is a ‘singular’ system that is particularly sensitive to the accuracy of the closure within mean-field models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Near-Monochromatic Tuneable Cryogenic Niobium Electron Field Emitter

We report creating, manipulating, and detecting coherent electrons is at the heart of future quantum microscopy and spectroscopy technologies. Leveraging and specifically altering the quantum features of an electron beam source at low temperatures can enhance its emission properties. Here, we describe electron field emission from a monocrystalline, superconducting niobium nanotip at a temperature of 5.9 K. The emitted electron energy spectrum reveals an ultranarrow distribution down to 16 meV due to tunable resonant tunneling field emission via localized band states at a nanoprotrusion’s apex and a cutoff at the sharp low-temperature Fermi edge. This is an order of magnitude lower than for conventional field emission electron sources. The self-focusing geometry of the tip leads to emission in an angle of 3.7°, a reduced brightness of 3.8 × 10 8 A/(m 2 srV), and a stability of hours at 4.1 nA beam current and 69 meV energy width. This source will decrease the impact of lens aberration and enable new modes in low-energy electron microscopy, electron energy loss spectroscopy, and high-resolution vibrational spectroscopy.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Space Charge Simulations of High Intensity Proton Beams in the AGS Booster

Computer simulation studies have been performed to understand the beam behavior and to explore intensity limitations of proton beams in the AGS Booster at higher beam intensities. During the 100 GeV polarized proton operations of RHIC Run 2024, sPHENIX operated in modes with a crossing angle at collisions in order to mitigate beam-beam effects. Three different running modes were employed: (a) sPHENIX operated with a negative (-2 mrad) crossing angle, and STAR operated with 0 mrad. Both experiments were brought into collisions at the start of the store. (b) sPHENIX was brought into collisions with 0 mrad first. Then STAR was brought into collisions after the beam-beam parameter from sPHENIX reduced to below $10 \times 10^{-3}$. (c) sPHENIX operated with a positive (+1.5 mrad) crossing angle, and STAR operated with 0 mrad. Both experiments were brought into collisions at the start of the store. The collisions with a crossing angle of up to $\pm 2$ mrad, as in running modes (a) and (c), lead to large Piwinski angle in the new sPHENIX detector, which reduces luminosity if other parameters are unchanged. There are two ways to compensate the reduction in luminosity: squeeze $\beta^{*}$ if there is sufficient dynamic aperture, or increase the injected beam intensity. The first part of polarized proton operation during RHIC Run 2024 was dedicated to increasing the intensity. Different configurations were tested with crossing angle and lattice adjustments on RHIC. At the same time, new injector configurations were developed and tested in an effort to push for both higher intensity and better quality of the beam injected into RHIC. When the beam intensity is increased, space charge is a concern particularly in the lower energy stages of acceleration, such as during the injection and the early part of the Booster cycle, which could become a dominant effect in limiting the intensity of the beam that can be delivered to RHIC.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Moment expansion of polarized dust SED: A new path towards capturing the CMB B -modes with LiteBIRD

Accurate characterization of the polarized dust emission from our Galaxy will be decisive in the quest for the cosmic microwave background (CMB) primordial B-modes. An incomplete modeling of its potentially complex spectral properties could lead to biases in the CMB polarization analyses and to a spurious measurement of the tensor-to-scalar ratio r. It is particularly crucial for future surveys like the LiteBIRD satellite, the goal of which is to constrain the faint primordial signal leftover by inflation with an accuracy on the tensor-to-scalar ratio r of the order of 10 –3 . Variations of the dust properties along and between lines of sight lead to unavoidable distortions of the spectral energy distribution (SED) that cannot be easily anticipated by standard component-separation methods. This issue can be tackled using a moment expansion of the dust SED, an innovative parametrization method imposing minimal assumptions on the sky complexity. In the present paper, we apply this formalism to the B-mode cross-angular power spectra computed from simulated LiteBIRD polarization data at frequencies between 100 and 402 GHz that contain CMB, dust, and instrumental noise. The spatial variation of the dust spectral parameters (spectral index β and temperature T) in our simulations lead to significant biases on r (~21 σ r ) if not properly taken into account. Performing the moment expansion in β, as in previous studies, reduces the bias but does not lead to sufficiently reliable estimates of r. We introduce, for the first time, the expansion of the cross-angular power spectra SED in both β and T, showing that, at the sensitivity of LiteBIRD, the SED complexity due to temperature variations needs to be taken into account in order to prevent analysis biases on r. Thanks to this expansion, and despite the existing correlations between some of the dust moments and the CMB signal responsible for a rise in the error on r, we can measure an unbiased value of the tensor-to-scalar ratio with a dispersion as low as σ r = 8.8 × 10 –4 .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Canonical and noncanonical Hamiltonian operator inference

Here, a method for the nonintrusive and structure-preserving model reduction of canonical and noncanonical Hamiltonian systems is presented. Based on the idea of operator inference, this technique is provably convergent and reduces to a straightforward linear solve given snapshot data and gray-box knowledge of the system Hamiltonian. Examples involving several hyperbolic partial differential equations show that the proposed method yields reduced models which, in addition to being accurate and stable with respect to the addition of basis modes, preserve conserved quantities well outside the range of their training data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Mitigation of mode-one asymmetry in laser-direct-drive inertial confinement fusion implosions

Nonuniformities present in the laser illumination and target in laser-driven inertial confinement fusion experiments lead to an asymmetric compression of the target, resulting in an inefficient conversion of shell kinetic energy to thermal energy of the hot-spot plasma. In this paper, the effects of asymmetric compression of cryogenic deuterium tritium laser-direct-drive implosions are examined using a suite of nuclear and x-ray diagnostics on the OMEGA laser. The neutron-averaged hot-spot velocity (~uhs) and apparent ion temperature (Ti) asymmetry are determined from neutron time-of-flight measurements of the primary deuterium tritium fusion neutron energy spectrum, while the areal density (rhoR) of the compressed fuel surrounding the hot spot is inferred from measurements of the scattered neutron energy spectrum. The low-mode perturbations of the hot-spot shape are characterized from x-ray self-emission images recorded along three quasi-orthogonal lines of sight. Implosions with significant mode-1 laser drive asymmetries show large hot-spot velocities (>100 km/s) in a direction consistent with the hot-spot elongation observed in x-ray images, measured Ti asymmetry, and rhoR asymmetry. Laser drive corrections have been applied through shifting the initial target location in order to mitigate the observed asymmetry. With the asymmetry corrected, a more-symmetric hot spot is observed with reduced ~uhs, Ti asymmetry, rhoR asymmetry, and a 30% increase in the fusion yield.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Mitigation of mode-one asymmetry in laser-direct-drive inertial confinement fusion implosions

Nonuniformities present in the laser illumination and target in laser-driven inertial confinement fusion experiments lead to an asymmetric compression of the target, resulting in an inefficient conversion of shell kinetic energy to thermal energy of the hot-spot plasma. In this paper, the effects of asymmetric compression of cryogenic deuterium tritium laser-direct-drive implosions are examined using a suite of nuclear and x-ray diagnostics on the OMEGA laser. The neutron-averaged hot-spot velocity (u → hs ) and apparent ion temperature (T i ) asymmetry are determined from neutron time-of-flight measurements of the primary deuterium tritium fusion neutron energy spectrum, while the areal density (ρR) of the compressed fuel surrounding the hot spot is inferred from measurements of the scattered neutron energy spectrum. Here, the low-mode perturbations of the hot-spot shape are characterized from x-ray self-emission images recorded along three quasi-orthogonal lines of sight. Implosions with significant mode-1 laser drive asymmetries show large hot-spot velocities (>100 km/s) in a direction consistent with the hot-spot elongation observed in x-ray images, measured T i asymmetry, and ρR asymmetry. Laser drive corrections have been applied through shifting the initial target location in order to mitigate the observed asymmetry. With the asymmetry corrected, a more-symmetric hot spot is observed with reduced u → hs , T i asymmetry, ρR asymmetry, and a 30% increase in the fusion yield.

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

Optimization-Based Model Reduction Scheme for Renewable Energy Power Plants Using Standardized Testing Scenarios

This paper presents an optimization-based model reduction scheme for renewable energy (RE) power plants consisting of inverter-based resources (IBRs) operating in grid-following (GFL) or grid-forming (GFM) modes. More importantly, the datasets feeding the optimization-based model reduction scheme are generated and re-used through the standardized grid-interactive testing scenarios. Particularly, the proposed scheme makes use of the power plant point of common coupling (PCC) measurements of various quantities specified by standardized tests (e.g., voltage and frequency ride through) as per IEEE 2800, to estimate the parameters of the reduced-order model such that its dynamic performance aligns with the original detailed power plant model. The proposed model reduction approach does not require the parameters of individual IBRs and using standardized test data as input to the formulated optimization problem simplifies the reduced-order modelling scheme. Extensive case studies following standardized test scenarios verified the remarkable accuracy of the proposed approach.

Yallamilli, Ram S. [Purdue University]↗
DOE OSTI2026

Nonlinear magnetohydrodynamic modeling of ideal ballooning modes in high- β Wendelstein 7-X plasmas

We present nonlinear magnetohydrodynamic (MHD) simulations of high- β Wendelstein 7-X plasmas using the stellarator extension of the M3D- C 1 code, building on the recent work that shows benign saturation of ideal ballooning modes above the designed β limit in the standard configuration [Zhou et al., Phys. Rev. Lett. 133, 135102 (2024)]. First, we examine the results' sensitivity to the parallel thermal conductivity. It is found that while an increased parallel conductivity reduces the linear growth rate, the saturated pressure profile is barely affected. Second, we consider the dependence on the profile shape. It is shown that an equilibrium with a peaked pressure profile and lower β is subject to more significant change than a broad profile with higher β and a larger growth rate, suggesting that benign saturation, or nonlinear stability, is not guaranteed and not dictated by linear growth. Third, we study the influence of the magnetic configuration, with the equilibrium rotational transform varied by adjusting the planar coil current. With similar growth rates, similar magnitudes of profile change are found regardless of the presence of a low-order resonance, which implies that the saturation mechanism is not specific to a resonant or non-resonant mode. These results indicate that MHD stability should still be treated seriously in stellarator operation and design, for which nonlinear modeling using tools like M3D- C 1 can play an instrumental role.

Zhou, Yao [Shanghai Jiao Tong University (China)] ↗