Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nonlinear 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,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

The Importance of an Ensemble Approach for Modeling Aerosol‐Convection Interactions

Aerosol-convection interactions modulate cloud microphysics, thermodynamics, and updraft intensity, contributing to climate-scale aerosol-radiative forcing. However, quantifying aerosol indirect effects in mixed-phase deep convection remains challenging due to uncertainties in parameterized physics and initial conditions driving nonlinear evolution of convective processes. This study investigates the convective updraft sensitivity to random initial temperature perturbations using an idealized ensemble modeling framework informed by in situ thermodynamic and aerosol observations from the DOE TRACER field campaign. We analyze the impact of small-scale initial perturbations on updraft velocity and supersaturation, and determine the ensemble size required to minimize stochastic internal variability. Results show that minor thermodynamic perturbations can produce updraft variability comparable to aerosol-induced changes reported in prior work. An ensemble of 10 members sufficiently reduces variability, enabling robust investigation of aerosol-related updraft invigoration and informing the design of future ensemble-based aerosol-convection interaction studies in environments with significant spatiotemporal mesoscale thermodynamic and aerosol heterogeneity.

aerosol-convection interactions↗
DOE OSTI2024

Unraveling the magnetic ground state in the alkali-metal lanthanide oxide Na 2 ⁢Pr⁡O 3

Here, a comprehensive set of muon spin spectroscopy and neutron scattering measurements supported by ab initio and model Hamiltonian simulations have been used to investigate the magnetic ground state of Na 2 ⁢PrO 3 . μSR reveals a Néel antiferromagnetic order below T N ~ 4.9K, with a small static magnetic moment m static ≤ 0.22 μ B /Pr collinearly aligned along the c axis. Inelastic neutron measurements reveal the full spectrum of crystal field excitations and confirm that the Pr 4+ ground-state wave function deviates significantly from the Γ 7 limit that is relevant to the Kitaev model. Single- and two-magnon excitations are observed in the ordered state below T N =4.6K and are well described by nonlinear spin wave theory from the Néel state using a magnetic Hamiltonian with Heisenberg exchange J=1 meV and symmetric anisotropic exchange Γ/J=0.1, corresponding to an XY model. Intense two magnon excitations are accounted for by g-factor anisotropy g z /g ± = 1.29. A fluctuating moment δm 2 = 0.57 ⁢(22) ⁢μ$^2_B$/Pr extracted from the energy and momentum integrated inelastic neutron signal is reduced from expectations for a local J = 1/2 moment with average g factor g avg ≈ 1.1. Together, the results demonstrate that the small moment in Na 2 ⁢PrO 3 arises from crystal field and covalency effects and the material does not exhibit significant quantum fluctuations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A comparative study of two numerical approaches for solving Kim–Kim–Suzuki phase-field models

Among the standard multi-phase multi-component phase-field (PF) methods, the Kim–Kim–Suzuki (KKS) method has the advantage of decoupling interfacial energy from bulk energy and solving concentration as the conserved variable. There are two approaches to numerically solving a KKS method: the global solution approach (GSA) solves all variables in a global system simultaneously, and the local solution approach (LSA) solves phase concentrations locally using a Newton solver. This work compares the performance of LSA and GSA for solving four KKS models of increasing complexity with the finite element method using the MOOSE framework. The solution accuracy, degrees of freedom (DOFs), memory usage, and computational efficiency are compared. We find that GSA and LSA generate similar solutions, with a maximum difference of only 0.34%. For each model, LSA has a lower number of DOFs, utilizes less memory, and less wall time. Additionally, the savings of memory and wall time in LSA increase with increasing mesh density of the same model and are more pronounced in models with higher dimensionality and more nodes. However, GSA is easier to implement in existing codes and can better solve highly nonlinear systems by utilizing sophisticated solvers.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Effect of net direct current on the properties of radio frequency sheaths: simulation and cross-code comparison

In order to understand, predict and control ion cyclotron range of frequency (ICRF) interactions with tokamak scrape-off layer plasmas, computational tools which can model radio frequency (RF) sheaths are needed. In particular, models for the effective surface impedance and DC rectified sheath potentials may be coupled with full wave RF simulation codes to predict self-consistent wave fields near surfaces and the resulting power dissipation and plasma-material interactions from ion sputtering. In this study, previous work assuming zero net DC current flow through the sheath is generalized to allow the surface to collect net positive or negative current, as is often observed in experiments. The waveforms, DC potential and RF admittance are investigated by means of analytical theory, nonlinear fluid and particle-in-cell (PIC) codes. Cross-code comparisons provide detailed model verification and elucidate the roles of ion and electron kinetics. When the sheath draws negative (positive) DC current, the voltage rectification is reduced (increased) compared with the zero-current case, and the magnitude of both the real and imaginary parts of the admittance are increased (reduced). In conclusion, a previous four-input parametrization of the sheath rectification and admittance properties is generalized to include a fifth parameter describing the DC sheath current.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Artificial Intelligence/Machine Learning Technologies for Advanced Reactors (Workshop Summary Report)

A workshop on artificial intelligence and machine learning (AI/ML) for advanced reactors (AR) was held October 5-6, 2021. The workshop was to be attended in-person at ANL but COVID restrictions forced the workshop to go virtual. The objectives of the workshop were to identify the most promising AI/ML opportunities for improving advanced reactor design, optimizing plant performance, and enhancing economic competitiveness and to develop an understanding of the scientific, engineering and licensing challenges facing their application. The workshop planning committee included GAIN, EPRI and NEI and members of three national laboratories (ANL, INL, and ORNL). The workshop was attended by more than 200 individuals representing academic and scientific institutions and the nuclear power industry. The definition put forth for an AI/ML system was one that perceives its environment and takes actions that maximize its chance of achieving its goals. In this report AI/ML refers to next generation algorithms that include deep learning, statistical analysis and data analytics and associated scientific computing and their potential application to the design, licensing, operation and maintenance of ARs. These methods typically incorporate models built from process data and may also include data generated by simulations that represent the behavior of a system. The workshop was organized in response to the growing interest in application of AI/ML for improving the economic competitiveness of nuclear energy. Increasingly more resources are being allocated to investigating the benefits of AI/ML methods. The DOE created the Artificial Intelligence & Technology Office to promote their development. And within the Office of Nuclear Energy, resources have been allocated to explore and understand the potential benefits of AI/ML. Additionally, the national laboratories are strategically positioned with DOE computing facilities such as Summit, Perlmutter, Aurora and Frontier that support large-scale simulations, hybrid HPC models with AI surrogates, and the exploration of new types of generative models emerging from multi-model data streams and sources. The workshop was organized with members of the AR community to understand the effort and to identify the level of interest and progress in this emerging technology. The workshop discussions focused on identifying opportunities for AI/ML across diverse areas of the nuclear industry and identifying current scientific and engineering challenges for advanced reactors that might be addressed through transformational uses of AI/ML. Discussion panels focused on four high-interest technical domains for advanced reactors: design, maintenance and operations, energy storage, and materials. The results of those discussions are summarized in this report. This includes opportunities that were identified for exploiting AI techniques and methods to improve the efficacy and efficiency of reactor analysis and to improve the operation and optimization of advanced reactors. Advanced reactor developers expressed an interest in learning more about AI/ML methods and their application. This included understanding whether ML methods can provide an advantage over existing nonlinear data regression methods for collapsing high-fidelity simulation results into faster running models. A consensus emerged that AR advances planned for the next decade will benefit from the use of AI/ML tools. The need exists to understand and model complex systems across length scales and modalities. AI/ML is a tool for discovery that can yield a set of engineering principles for use by nuclear engineers, licensing bodies, and operators to solve problems in plant design, safety analyses, autonomous operation, and predictive maintenance. While AI/ML represents a new set of tools, an awareness by the nuclear community of the full potential is still in the early stages so there is a need to increase awareness. It appears that the wide-spread adoption of AI/ML tools for ARs would be facilitated by future educational workshops that describe foundational methods and capabilities and describe successful applications.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Mshock HED Council Request FY24 [Slides]

ICF and HEDP experiments experience complicated loadings including successive shocks, but modeling and analytic work has mostly engaged the simpler shock or shock-reshock case. This is because the co propagating case is difficult to achieve with conventional (non-HED) drivers. Successive shocks will be challenging to model in BHR or modal model interpretations: 1) Ex: Consider the case where long wavelength modes re-invert coherently, but short wavelength nonlinear modes are spun up turbulently. 2) Necessitates more advanced diagnostics and analysis than simple “mixing layers.” High resolution spectral information (concomitant with Ω-EP campaigns) and higher-order moment analyses will be necessary to stress and validate the reduced BHR-type models.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

Zonal Flow and Self-regulating Mechanism in a Hydrodynamic Disk

This study addresses key aspects of momentum transport in hydrodynamic disks, which is critical for understanding zonal flow generation and turbulence in compressible hydrodynamic disks. We find that nonlinear momentum/density transport leads to the formation of zonal flows from the Rossby wave instability in disks. We analytically derive the generation and location of zonal flows and describe a modified Taylor identity applicable to compressible disk flows. We further present a self-regulation model, revealing a dynamic interplay between zonal flow and fluctuations driven by Rossby wave instability that regulates the nonlinear saturation state. This theoretical framework contributes insights into the dynamics of disks such as protoplanetary disks, shedding light on the intricate processes governing momentum/density transport and the emergence of zonal flows in the saturation of protoplanetary disks.

Astrophysical fluid dynamics↗
DOE OSTI2025

The Observed and Projected Changes of Global Monsoons: Current Status and Future Perspectives

The global monsoon system, encompassing the Asian-Australian, African, and American monsoons, sustains two-thirds of the world’s population by regulating water resources and agriculture. Monsoon anomalies pose severe risks, including floods and droughts. Recent research associated with the implementation of the Global Monsoons Model Intercomparison Project under the umbrella of CMIP6 has advanced our understanding of its historical variability and driving mechanisms. Observational data reveal a 20th-century shift: increased rainfall pre-1950s, followed by aridification and partial recovery post-1980s, driven by both internal variability (e.g., Atlantic Multidecadal Oscillation) and external forcings (greenhouse gases, aerosols), while ENSO drives interannual variability through ocean-atmosphere interactions. Future projections under greenhouse forcing suggest long-term monsoon intensification, though regional disparities and model uncertainties persist. Models indicate robust trends but struggle to quantify extremes, where thermodynamic effects (warming-induced moisture rise) uniformly boost heavy rainfall, while dynamical shifts (circulation changes) create spatial heterogeneity. Volcanic eruptions and proposed solar radiation modification (SRM) further complicate predictions: tropical eruptions suppress monsoons, whereas high-latitude events alter cross-equatorial flows, highlighting unresolved feedbacks. The emergent constraint approach is booming in terms of correcting future projections and reducing uncertainty with respect to the global monsoons. Critical challenges remain. Model biases and sparse 20th-century observational data hinder accurate attribution. The interplay between natural variability and anthropogenic forcings, along with nonlinear extreme precipitation risks under warming, demands deeper mechanistic insights. Additionally, SRM’s regional impacts and hemispheric monsoon interactions require systematic evaluation. Addressing these gaps necessitates enhanced observational networks, refined climate models, and interdisciplinary efforts to disentangle multiscale drivers, ultimately improving resilience strategies for monsoon-dependent regions.

climate extreme events↗
DOE OSTI2021

Operation of Natural Gas Pipeline Networks With Storage Under Transient Flow Conditions

Here, we formulate a nonlinear optimal control problem for intraday operation of a natural gas pipeline network that includes storage reservoirs. The dynamics of compressible gas flow through pipes, compressors, reservoirs, and wells are considered. In particular, a reservoir is modeled as a rigid, hollow container that stores gas under isothermal conditions and uniform density, and a well is modeled as a vertical pipe. For each pipe, flow dynamics are described by a coupled partial differential equation (PDE) system in density and mass flux variables, with momentum dissipation modeled using the Darcy–Wiesbach friction approximation. Compressors are modeled as scaling up the pressure of gas between the inlet and outlet. The governing equations for all network components are spatially discretized and assembled into a nonlinear differential-algebraic equation (DAE) system, which synthesizes above-ground pipeline and subsurface reservoir dynamics into a single reduced-order model. We seek to maximize an objective function that quantifies economic profit and network efficiency subject to the flow equations and inequalities that represent operating limitations. The problem is solved using a primal–dual interior point solver, and the solutions are validated in computational experiments and simulations on several pipeline test networks to demonstrate the effectiveness of the proposed methodology.

03 NATURAL GAS↗
DOE OSTI2026

Energy-Optimal Vehicle Longitudinal Motion Control via Pontryagin’s Minimum Principle and Ultra-Local Model

Longitudinal vehicle motion control is essential for enhancing performance and optimizing a vehicle’s energy usage. However, it remains a challenging task due to the nonlinear and uncertain nature of vehicle dynamics, along with varying driving conditions. This paper presents a novel ultra-local optimal control approach based on Pontryagin’s Minimum Principle (PMP) that circumvents the need for detailed system identification by employing an ultra-local model. The control objective is to minimize the total energy consumption under boundary conditions while ensuring smooth traction force generation. The proposed approach is evaluated using a high-fidelity vehicle model in three representative scenarios: (i) nominal driving, (ii) a change in tire road friction coefficient (TRFC) from 0.5 to 0.65 and road slope from 0% to 5% during the maneuver, with target velocity unchanged, and (iii) a change in target velocity from 20 m/s to 0 m/s during the maneuver, while maintaining nominal TRFC and slope conditions. The simulation results demonstrate that the proposed method delivers robust performance, effectively balancing consumption and tracking accuracy in all tested scenarios.

Waleed khan, Muhammad [The University of Texas at ↗
DOE OSTI2021

Reduced-order modeling of advection-dominated systems with recurrent neural networks and convolutional autoencoders

A common strategy for the dimensionality reduction of nonlinear partial differential equations (PDEs) relies on the use of the proper orthogonal decomposition (POD) to identify a reduced subspace and the Galerkin projection for evolving dynamics in this reduced space. However, advection-dominated PDEs are represented poorly by this methodology since the process of truncation discards important interactions between higher-order modes during time evolution. In this study, we demonstrate that encoding using convolutional autoencoders (CAEs) followed by a reduced-space time evolution by recurrent neural networks overcomes this limitation effectively. We demonstrate that a truncated system of only two latent space dimensions can reproduce a sharp advecting shock profile for the viscous Burgers equation with very low viscosities, and a six-dimensional latent space can recreate the evolution of the inviscid shallow water equations. Additionally, the proposed framework is extended to a parametric reduced-order model by directly embedding parametric information into the latent space to detect trends in system evolution. Furthermore, our results show that these advection-dominated systems are more amenable to low-dimensional encoding and time evolution by a CAE and recurrent neural network combination than the POD-Galerkin technique.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Dark Matter Constraints from Small-Scale Cosmic Structure

Small-scale cosmic structure provides a powerful test of the fundamental nature of dark matter (DM). A wide range of DM models impact matter clustering on small scales, including warm, fuzzy, and (self-)interacting DM. In these scenarios, DM physics such as free-streaming, wave interference, and self/Standard Model interactions alter the abundance and internal structure of DM halos. Cosmological and astrophysical probes of nonlinear structure---including dwarf galaxies, strong lensing, the Lyman-$α$ forest, stellar streams, and high-redshift galaxies---are therefore sensitive to these effects. Here, we review DM constraints provided by small-scale structure, focusing on observables that probe scales smaller than $\sim 1~\mathrm{Mpc}$, which define the frontier of current measurements. We summarize how these constraints have been translated to limits on microphysical DM models, and we discuss key modeling uncertainties and observational systematics. Finally, we highlight the growing importance of probe combination and simulation-based inference for this field, and we overview upcoming observational facilities that will sharpen small-scale structure tests of DM physics.

Nadler, Ethan O. [UC, San Diego] (ORCID:0000000211↗
DOE OSTI2020

Improved single-particle phase-space distributions for viscous fluid dynamic models of relativistic heavy ion collisions

Comparison of hydrodynamic calculations with experimental data inevitably requires a model for converting the fluid to particles. For shear viscous fluids, several different conversion models have been proposed, such as the quadratic Grad corrections, the Strickland-Romatschke (SR) ansatz, self-consistent shear corrections from linearized kinetic theory, or corrections from the relaxation time approach. We show, via comparison to nonlinear 2 → 2 kinetic theory evolution, that even for low specific shear viscosities 0.03 ≲ $\textit{η/s}$ ≲ 0.2 , these four models have varying accuracy and still significant errors in reproducing particle phase space densities from hydrodynamic fields. Moreover, we demonstrate that the overall reconstruction error of additive shear viscous $f = f_{eq} + δf$ models dramatically reduces by up to one order of magnitude, if one ensures through exponentiation that $\textit{f}$ is always positive. Further, we also illustrate how even more accurate viscous $\textit{δf}$ models can be constructed by incorporating the first time derivatives of hydrodynamic fields to account for the past evolution of the system. Such time derivatives are readily available in hydrodynamic simulations, though usually not included in the output. Though our comparisons are limited to a one-component massless system undergoing a (0 + 1)-dimensional (0 + 1D) longitudinal boost-invariant expansion, we expect that the improved models will be useful in 2 + 1D and 3 + 1D hydrodynamic studies as well.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Generation of Subauroral Longitudinally Extended Emissions Following Intensifications of the Poleward Boundary of the Substorm Bulge and Streamer Production

Observations showing the development of Subauroral Longitudinally Extended Emigerssions on a global scale are presented. It is demonstrated that they occur as a separation of an east-west arc-like band of luminosity that detaches away from the equatorward edge of the auroral distribution following episodes of auroral streamer production. They persist for time intervals on the order of ~30 min and devolve into patchy disjointed segments before they fade. Emissions are seen in both 557.7 nm OI and 391.4 nm N + 2 1NG lines, and 630.0 nm emissions are often observed equatorward. The 391.4 nm emissions are typically weaker and fade away more quickly than the 557.7 nm emissions, suggesting that an auroral precipitation source is initially present but is rapidly depleted as the forms age. All cases are associated with enhanced Subauroral Polarization Stream (SAPS) flows, and one event shows clear association with large-scale Giant Undulations (GUs) and the formation of spur-like forms that fold equatorward under the main oval. A model is proposed for the production of subauroral longitudinally extended forms in which nonlinear growth of SAPS-induced surface waves on the plasmapause results in a disruption of the boundary. It is suggested that hot plasma of plasma sheet origin becomes entrained in the plasmasphere to produce transient precipitation-associated auroral emissions that may decay into STEVE emissions and that cold dense plasma from the plasmasphere becomes entrained onto open drift paths to feed long-lived drainage plumes. This process may occur quasi periodically during intense substorms and storm-time conditions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Axisymmetric simulations of vertical displacement events in tokamaks: A benchmark of M3D-C 1 , NIMROD, and JOREK

A benchmark exercise for the modeling of vertical displacement events (VDEs) is presented and applied to the 3D nonlinear magnetohydrodynamic codes M3D-C 1 , JOREK, and NIMROD. The simulations are based on a vertically unstable NSTX equilibrium enclosed by an axisymmetric resistive wall with a rectangular cross section. A linear dependence of the linear VDE growth rates on the resistivity of the wall is recovered for sufficiently large wall conductivity and small temperatures in the open field line region. We report the benchmark results show good agreement between the VDE growth rates obtained from linear NIMROD and M3D-C 1 simulations and from the linear phase of axisymmetric nonlinear JOREK, NIMROD, and M3D-C 1 simulations. Axisymmetric nonlinear simulations of a full VDE performed with the three codes are compared, and an excellent agreement is found regarding the plasma location and plasma currents, as well as eddy and halo currents in the wall.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Electromagnetic pinned solitons for space debris detection

Electromagnetic “pinned” solitons in the form of stationary nonlinear waves are studied within the framework of an inertial magneto-hydrodynamic model. These structures, that can arise when a charged source moves in a magnetized plasma, have a velocity that is equal to the source velocity and, hence, appear as “pinned” structures that envelope the source. We investigate the excitation of such solitons in the Low Earth Orbit region due to the passage of charged orbital debris objects. The spatial size of these electromagnetic solitons, typically of the order of a few ion skin depths, can be very large in this region. Such solitons can be detected using a variety of ground- or orbit-based radio sounding techniques and may provide a convenient additional means of tracking small sized orbital debris objects that are difficult to spot optically.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Subspace recursive Fermi-operator expansion strategies for large-scale DFT eigenvalue problems on HPC architectures

Quantum mechanical calculations for material modeling using Kohn–Sham density functional theory (DFT) involve the solution of a nonlinear eigenvalue problem for N smallest eigenvector-eigenvalue pairs, with N proportional to the number of electrons in the material system. Here, these calculations are computationally demanding and have asymptotic cubic scaling complexity with the number of electrons. Large-scale matrix eigenvalue problems arising from the discretization of the Kohn–Sham DFT equations employing a systematically convergent basis traditionally rely on iterative orthogonal projection methods, which are shown to be computationally efficient and scalable on massively parallel computing architectures. However, as the size of the material system increases, these methods are known to incur dominant computational costs through the Rayleigh–Ritz projection step of the discretized Kohn–Sham Hamiltonian matrix and the subsequent subspace diagonalization of the projected matrix. This work explores the potential of polynomial expansion approaches based on recursive Fermi-operator expansion as an alternative to the subspace diagonalization of the projected Hamiltonian matrix to reduce the computational cost. Subsequently, we perform a detailed comparison of various recursive polynomial expansion approaches to the traditional approach of explicit diagonalization on both multi-node central processing unit and graphics processing unit architectures and assess their relative performance in terms of accuracy, computational efficiency, scaling behavior, and energy efficiency.

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

Degenerate parametric down-conversion facilitated by exciton-plasmon polariton states in a nonlinear plasmonic cavity

Abstract We study the effect of degenerate parametric down-conversion (DPDC) in an ensemble of two-level quantum emitters (QEs) coupled via near-field interactions to a single surface plasmon (SP) mode of a nonlinear plasmonic cavity. For this purpose, we develop a quantum driven-dissipative model capturing non-equilibrium dynamics of the system in which incoherently pumped QEs have transition frequency tuned near the second-harmonic response of the SPs. Considering the strong coupling regime, i.e. the SP-QE interaction rate exceeds system dissipation rates, we find a critical SP-QE coupling attributed to the phase transition between normal and lasing steady states. Examining fluctuations above the system’s steady states, we predict new elementary excitations, namely, the exciton-plasmon polaritons formed by the two-SP quanta and single-exciton states of QEs. The contribution of two-SP quanta results in the linear scaling of the SP-QE interaction rate with the number of QEs,  o , as opposed to the  o -scaling known for the Dicke and Tavis–Cummings models. We further examine how SP-QE interaction scaling affects the polariton dispersions and power spectra in the vicinity of the critical coupling. For this purpose, we compare the calculation results assuming a finite ensemble of QEs and the model thermodynamic limit. The calculated power spectra predict an interplay of coherent photon emission by QEs near the second-harmonic frequency and correlated photon-pair emission at the fundamental frequency by the SPs (i.e. the photonic DPDC effect).

77 NANOSCIENCE AND NANOTECHNOLOGY↗