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 649 records · Page 36

Prediction of residual stresses in additively manufactured parts using lumped capacitance and classical lamination theory

Several industries are interested in Laser Powder Bed Fusion (L-PBF) Additively Manufactured (AM) metal parts because their designs can be made arbitrarily complex while retaining bulk-type material properties. However, the residual stresses (RS) and distortions caused by the heat gradients inherent to L-PBF processes are detrimental to the structural integrity of the parts and must be taken into consideration during the part design cycle. Predicting the state of stresses in as-built 3D printed parts is a difficult problem that is typically approached with the use of transient thermomechanical Finite Element Models (FEMs). However, the nonlinearities associated with AM processes are difficult to capture in these FEMs without increasing the computational cost of the simulation, limiting their ability to be incorporated into practical design cycles. This work presents a novel analytical framework that combines lumped capacitance nonlinear heat transfer with time dependent classical lamination theory to efficiently and accurately predict RS in as-built L-PBF parts without the need of FEMs. The simulation was compared to Neutron Diffraction (ND) residual strain measurements taken at Oak Ridge National Laboratories (ORNL) as well as Synchrotron X-ray Diffraction (XRD) strain data published by the National Institute of Standards and Technology (NIST). The simulation predictions and the experimental data showed excellent agreement for the in-plane strain directions, and general agreement for the out of plane strain component, highlighting an area where further development can be implemented.

42 ENGINEERING↗

Analysis of degradation in residential battery energy storage systems for rate-based use-cases

This article examines the impact of residential battery energy storage (BES) systems' operational modes on the life (i.e. usable energy capacity) of the battery under several climatic conditions and battery chemistries. The sharp increase in residential BES installations has been a result of decreasing costs of batteries, increase in rate structures motivated applications such as solar self-consumption and time-of-use energy management, and customers purchasing these systems for backup power. While these different modes of operations provide a combination of increased bill savings, reliability of supply, and energy sustainability to the customer, their operational characteristics vary significantly between use-cases. Functional life of BES systems is known to be strongly dependent on their operating conditions. In this article, we analyze the operation of residential BES systems under different rate-based use-cases, for different battery chemistries and cell designs, and under different environmental conditions. This is conducted by simulating the control of BES operations using rate-based cycling algorithms and analyzing the prognosis of multiple battery lifetime models that consider complex nonlinear dependencies of operational stress factors such as state-of-charge, depth-of-discharge, and temperature on degradation. Significant variations in battery life are observed owing to the differences in characteristics of the uses-cases coupled with environmental conditions and battery chemistries. Such differences lead to the conclusion that choice of residential BES technologies and chemistries should account for their intended use-cases.

25 ENERGY STORAGE↗

Operator inference for non-intrusive model reduction with quadratic manifolds

This paper proposes a novel approach for learning a data-driven quadratic manifold from high-dimensional data, then employing this quadratic manifold to derive efficient physics-based reduced-order models. The key ingredient of the approach is a polynomial mapping between high-dimensional states and a low-dimensional embedding. This mapping consists of two parts: a representation in a linear subspace (computed in this work using the proper orthogonal decomposition) and a quadratic component. The approach can be viewed as a form of data-driven closure modeling, since the quadratic component introduces directions into the approximation that lie in the orthogonal complement of the linear subspace, but without introducing any additional degrees of freedom to the low-dimensional representation. Combining the quadratic manifold approximation with the operator inference method for projection-based model reduction leads to a scalable non-intrusive approach for learning reduced-order models of dynamical systems. Further, applying the new approach to transport-dominated systems of partial differential equations illustrates the gains in efficiency that can be achieved over approximation in a linear subspace.

42 ENGINEERING↗

Mathematical Programming Models for Shale Oil & Gas Development: A Review and Perspective

Here, in this paper, we provide a comprehensive review of mathematical programming models for shale oil & gas development, and we offer a perspective on outstanding research opportunities. We distinguish contributions in five major topic areas, namely: (1) development planning, (2) water management, (3) production optimization, (4) supplies, gathering & processing, and (5) life cycle analysis & sustainability. We highlight how various types of mathematical programming models (i.e., linear programs, nonlinear programs, mixed-integer linear programs, mixed-integer nonlinear programs) have been proposed primarily by the Process Systems Engineering community to address the respective decision-making problems, and we highlight instances of successful deployment in industry. Finally, based on a critical assessment of the existing body of work, we identify opportunities for future research across the major topic areas.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Adjoint DSMC for nonlinear spatially-homogeneous Boltzmann equation with a general collision model

We derive an adjoint method for the Direct Simulation Monte Carlo (DSMC) method for the spatially homogeneous Boltzmann equation with a general collision law. This generalizes our previous results in Caflisch et al., which was restricted to the case of Maxwell molecules, for which the collision rate is constant. The main difficulty in generalizing the previous results is that a rejection sampling step is required in the DSMC algorithm in order to handle the variable collision rate. We find a new term corresponding to the so-called score function in the adjoint equation and a new adjoint Jacobian matrix capturing the dependence of the collision parameter on the velocities. The new formula works for a much more general class of collision models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

PreMevE Update: Forecasting Ultra-Relativistic Electrons Inside Earth's Outer Radiation Belt

Energetic electrons inside Earth's Van Allen belts pose a major radiation threat to spaceborne electronics that often play vital roles in modern society. Ultra-relativistic electrons with energies greater than or equal to two megaelectron-volt (MeV) are of particular interest, and thus forecasting these ≥2 MeV electrons has a significant meaning to all space sectors. Here, we update the latest development of the predictive model for MeV electrons in the outer radiation belt. The new version, called PREdictive MEV Electron (PreMevE)-2E, forecasts ultra-relativistic electron flux distributions across the outer belt, with no need for in situ measurements of the trapped MeV electron population except at the geosynchronous orbit (GEO). Model inputs include precipitating electrons observed in low-Earthorbits by NOAA satellites, upstream solar wind speeds and densities from solar wind monitors, as well as ultra-relativistic electrons measured by one Los Alamos GEO satellite. We evaluated 32 supervised machine learning models that fall into four different classes of linear and neural network architectures, and successfully tested ensemble forecasting by using groups of top-performing models. All models are individually trained, validated, and tested by in situ electron data from NASA's Van Allen Probes mission. It is shown that the final ensemble model outperforms individual models at most L-shells, and this PreMevE-2E model can provide 25-h (~1-day) and 50-h (~2-day) forecasts with high mean performance efficiency and correlation values. Our results also suggest that this new model is dominated by nonlinear components at L-shells <~4 for ultra-relativistic electrons, different from the dominance of linear components for 1 MeV electrons as previously discovered.

79 ASTRONOMY AND ASTROPHYSICS↗

fluxfinder: An R Package for Reproducible Calculation and Initial Processing of Greenhouse Gas Fluxes From Static Chamber Measurements

Fluxes of greenhouse gases are a critical component of the earth's natural climate, but anthropogenic emissions have created an imbalance and resulted in global climate change. Quantifying the emission of these gases is vital to our understanding of their sources and sinks, both natural and anthropogenic. The static chamber method, in which a system of interest is enclosed, and gas concentrations are measured over time, is widely used to estimate fluxes of greenhouse gases. With the development of instruments such as infrared gas analyzers (IRGAs) supporting high-frequency concentration data, there is a growing need for open-source workflows to calculate fluxes. Here we present fluxfinder, an R package designed to support reproducible calculations and processing of greenhouse gas fluxes measured with the static chamber method. The package includes raw data file parsing from widely used IRGAs, metadata matching, unit conversion, flux estimations, and initial quality assurance/quality control (QA/QC). Diagnostic graphical plots provide a transparent way to differentiate between measurement issues and nonlinear behavior. The package is also designed to be easily integrated with the gasfluxes package for further fitting of nonlinear concentration-time models, allowing alternative or additional flux QA/QC. The fluxfinder package offers a flexible workflow that is easily adaptable to promote open and reproducible greenhouse gas flux estimations.

Wilson, Stephanie J.↗

Resolving ECRH deposition broadening due to edge turbulence in DIII-D

Microwave heat pulse propagation experiments have demonstrated a correlation between millimeter-scale turbulence and deposition profile broadening of electron cyclotron (EC) waves on the DIII-D tokamak. In a set of discharges in DIII-D, a variation in edge density fluctuations on the mm-scale is associated with 40%–150% broader deposition profiles, expressed in terms of normalized minor radius, as compared with equilibrium ray tracing. The 1D power profile is determined from transport analysis of the electron temperature response to EC power modulation using perturbative analysis with a square wave power modulation at 20–70 Hz, producing a series of Fourier harmonics that are fit collectively to resolve transport. Fitting an integrated heat flux expressed in the Fourier basis of the modulation to diffusive, convective, and coupled transport terms in a linear model can resolve the broadened EC deposition width from the power perturbation to resolve a broadening in each case. The best fit degree of beam broadening observed scales approximately linearly with the Doppler backscattering measured fluctuation level in the steep gradient region. Quantifying the effect of edge fluctuation broadening on EC current drive power needs of future devices will require 3D full-wave codes that can be validated on the current generation of machines. These DIII-D experiments provide a quantitative measure of fluctuation effects and a dataset to benchmark full-wave simulations that can model and eventually predict nonlinear effects neglected by 1D equilibrium beam and ray tracing.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The role of collective elasticity on activated structural relaxation, yielding, and steady state flow in hard sphere fluids and colloidal suspensions under strong deformation

In this work, we theoretically study the effect of external deformation on activated structural relaxation and aspects of the nonlinear mechanical response of glassy hard sphere fluids in the context of elastically collective nonlinear Langevin equation theory. This microscopic force-based approach describes activated relaxation as a coupled local–nonlocal event involving caging and longer range collective elasticity, with the latter becoming more important and ultimately dominant with increasing packing fraction under equilibrium conditions. The central new question we address is how this physical picture of activated relaxation, and the relative importance of local caging vs collective elasticity physics, depends on external deformation. Theoretical predictions are presented for deformation-induced enhancement of mobility, the onset of relaxation speed up at remarkably low values of stress, strain, or shear rate, apparent power law thinning of the steady state structural relaxation time and viscosity, a non-vanishing activation barrier in the shear thinning regime, an apparent Herschel–Bulkley form of the rate dependence of the steady state shear stress, exponential growth of different measures of a dynamic yield or flow stress with the packing fraction, and reduced fragility and dynamic heterogeneity under deformation. The results are contrasted with experiments and simulations, and qualitative or better agreement is found. An overarching conclusion is that deformation strongly reduces the importance of longer range collective elastic effects relative to the local caging aspect for most, but not all, physical questions, with deformation-dependent fragility and dynamic heterogeneity phenomena being qualitatively sensitive to collective elasticity. Overall, nonlinear rheology is predicted to be a more local problem than quiescent structural relaxation, albeit with deformation-modified activated processes still important.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Consistent coupling algorithms for coupled core-edge simulations of plasma turbulence

Two schemes for coupling gyrokinetic simulations of microturbulence in tokamaks are proposed. The first scheme is based on an additive Schwarz domain decomposition. In this work, we show that, because the goal of turbulence is long-time averages of the dynamics rather than temporal accuracy, the iteration to self-consistency across domains, which is typically required by Schwarz schemes, can be avoided, thereby accelerating the computation. Second, we propose a coupling scheme that relies entirely on the addition of source terms, leaving the boundary conditions arbitrary. The practical motivations for such a scheme are discussed, and forms of the source terms that ensure consistency and stability are derived. The schemes are tested on a nonlinear, one-dimensional model problem, and the first scheme is further tested on the Hasegawa–Wakatani model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Deterministic scattering of relativistic electron beams by off-resonant circularly polarized electromagnetic waves

A mechanism is presented whereby relativistic electron beams localized in phase space are deterministically scattered by coherent circularly polarized electromagnetic waves without stochastic processes. It is shown via an exact single-particle analysis that the condition for maximal scattering is an off-resonant condition, contrary to previous kinetic analyses that predict maximal diffusion or interaction at exact resonance or its harmonics. The mechanism, verified by single-particle simulations, enables a fast, nonlinear redistribution of the beam particles. A possible application of this mechanism to runaway electron suppression is presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Machine learning methods for probabilistic locked-mode predictors in tokamak plasmas

A rotating tokamak plasma can interact resonantly with the external helical magnetic perturbations, also known as error fields. This can lead to locking and then to disruptions. We leverage machine learning (ML) methods to predict the locking events. We use a coupled third-order nonlinear ordinary differential equation model to represent the interaction of the magnetic perturbation and the plasma rotation with the error field. This model is sufficient to describe qualitatively the locking and unlocking bifurcations. Here, we explore using ML algorithms with the simulation data and experimental data, focusing on the methods that can be used with sparse datasets. These methods lead to the possibility of the avoidance of locking in real-time operations. We describe the operational space in terms of two control parameters: the magnitude of the error field and the rotation frequency associated with the momentum source that maintains the plasma rotation. The outcomes are quan- tified by order parameters that completely characterize the state, whether locked or unlocked. We use unsupervised ML methods to classify locked/unlocked states and note the usefulness of a certain normalization of the order parameters. Three supervised ML classifiers are used in suite to estimate the probability of locking in the region of control parameter space with hysteresis, i.e., the set of control parameters for which both locked and unlocked states can exist. The results show that a neural network gives the best estimate of the locking probability. An analogy of the present locking model with the van der Waals equation of state is also provided.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fluid and hybrid simulations of the ionization instabilities in Hall thruster

Low-frequency axial oscillations in the range of 5–50 kHz stand out as a pervasive feature observed in many types of Hall thrusters. While it is widely recognized that the ionization effects play the central role in this mode, as manifested via the large-scale oscillations of neutral and plasma density, the exact mechanism(s) of the instabilities remain unclear. To gain further insight into the physics of the breathing mode and evaluate the role of kinetic effects, a one-dimensional time-dependent full nonlinear low-frequency model describing neutral atoms, ions, and electrons is developed in full fluid formulation and compared to the hybrid model in which the ions and neutrals are kinetic. Both models are quasi-neutral and share the same electron fluid equations that include the electron diffusion, mobility across the magnetic field, and the electron energy evolution. The ionization models are also similar in both approaches. Further, the predictions of fluid and hybrid simulations are compared for different test cases. Two main regimes are identified in both models: one with pure low-frequency behavior and the other one, where the low-frequency oscillations coexist with high-frequency oscillations in the range of 100–200 kHz, with the characteristic time scale of the ion channel fly-by time, 100–200 kHz. The other test case demonstrates the effect of a finite temperature of injected neutral atoms, which has a substantial suppression effect on the oscillation amplitude.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Implicit full-F simulations of neoclassical ion transport

The development of implicit time integration capabilities for axisymmetric full-F continuum simulations of ion neoclassical transport is reported. The approach involves the implicit treatment of the gyrokinetic Vlasov equation coupled to the nonlinear Fokker–Planck collision model in the long-wavelength limit approximation. To facilitate implicit simulations, advanced preconditioning of individual physics operators is developed, and a global multi-physics preconditioner is constructed by adopting an operator splitting methodology. The algorithm is implemented in the finite-volume code COGENT and is applied to study neoclassical transport properties for both the main ion species and the lithium impurity species in the closed-field-line region of the LTX- β tokamak. The implicit COGENT simulations elucidate the role of non-local transport effects, while demonstrating substantial speedup over the corresponding explicit approach.

Dorf, Mikhail [Lawrence Livermore National Laborat↗

Non-disruptive error field measurement in DIII-D low safety factor plasmas and projection to ITER

Previous experiments in DIII-D (Paz-Soldan et al 2022 Nucl. Fusion 62 126007) introduced a method to identify intrinsic error fields (EFs) in tokamaks with minimal disruption risk by promptly healing driven magnetic islands during the conventional 'compass scan'. This paper presents recent experimental and numerical advancements in extending this approach to low q 95 plasmas, and projects its applicability to ITER. Non-disruptive EF measurement is achieved at q 95 = 4.5 and 3.9 without any initial EF correction (EFC) by reducing the time between the occurrence of the locked mode (LM) and control action to 10 ms and increasing the density 50%–100%. However, 50% correction of the intrinsic EF is required to achieve island healing at q 95 = 3.2 with 10 ms delay for the control action. Nonlinear two-fluid modeling with the TM1 code reproduces the DIII-D experimental observations, indicating that promptly turning off the 3D coil current reduces both magnetic island width and electromagnetic force, while raising the density increases plasma viscosity, facilitating magnetic island healing. The simulations show that for scenarios with q 95 = 3.2, lowering the control action time to 5 ms will lead to island healing without EFC. TM1 simulations are extended to future ITER scenarios with 5 MA and 7.5 MA plasma currents, predicting the dependence of required density rise on action time and EF amplitude. These simulations indicate that, benefiting from the much longer resistive time, island healing can be successfully achieved in ITER when taking control action 100–500 ms after a LM occurrence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Validation of 3D MHD simulations of Ne/D 2 mixed shattered pellet injection in JET

Nonlinear 3D magnetohydrodynamic (MHD) modeling of shattered pellet injection (SPI) in JET tokamak plasmas is performed with the JOREK code. The study focuses on the validation of simulation result with respect to experimental observations, addressing in particular figures of merit for the efficiency of the SPI, as a technique to mitigate thermal loads by injecting radiative impurities like neon during the thermal quench (TQ) phase of disruptions. A set of JET pulses with neon/deuterium atomic mixture ratio in the shattered pellet varying from 10% to 100% provides the experimental data to be compared with. Simulations using different models for the ablation of solid fragments and the radiation of neon impurities are considered. Synthetic diagnostics are employed for a direct quantitative comparison with key JET diagnostic systems providing in particular radiation, electron temperature and density, and magnetic measurements. 3D radiation structures from MHD simulations are analyzed to estimate the toroidal asymmetry of radiation, which is difficult to measure in JET due to the presence of only two toroidally displaced multi-channel bolometry systems. JOREK simulations show that, before the TQ, the bulk of the radiation is concentrated in a toroidal region enclosing the radiating fragments that is missed by both multi-channel bolometers, suggesting a possible underestimation of the total radiation by such diagnostic systems in the pre-TQ phase. On the other hand, JOREK predicts that during the TQ the radiation becomes more toroidally symmetric, with positive implications for the heat load to plasma facing components.

JOREK↗

Dark Energy Survey Year 6 results: Galaxy-galaxy lensing

Here, we present galaxy-galaxy lensing measurements from the full six years of data from the Dark Energy Survey (DES Y6), covering 4031 deg 2 and used in the DES Y6 3×2 pt cosmological analysis. We use the Maglim++ lens sample, containing ∼9 million galaxies divided into six redshift bins, and the Metadetection source catalog, including ∼140 million galaxies divided into four redshift bins. The mean tangential shear signal achieves a total signal-to-noise ratio (S/N) of 173, corresponding to a 17% improvement over DES Y3. After applying the scale cuts used in the cosmological analysis, with 𝑅 min =6 Mpc/ℎ (4 Mpc/ℎ) for the linear (nonlinear) galaxy-bias model, the S/N is reduced to 75 (90). A comprehensive suite of validation tests demonstrates that the measurement is robust against observational and astrophysical systematics at the statistical precision required for the DES Y6 analysis. Although not used in the main cosmological analysis, we extract high signal-to-noise geometric shear-ratio measurements from the galaxy-galaxy lensing signal on small angular scales. These measurements provide an internal consistency check on the photometric redshift distributions and shear calibration used in the 3 ×2 pt analysis.

cosmological parameters↗