Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear analysis”

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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

New Solids in As-O-Mo, As(P)-O-Mo(W) and As(P)-O-Nb(W) Systems That Exhibit Nonlinear Optical Properties

Interactions between well-mixed fine powders of As 2 O 3 , P 2 O 5 , MoO 3 , WO 3 and Nb 2 O 5 at different stoichiometry in quartz ampoules under vacuum at ~1000 °C in the presence of metallic molybdenum (or niobium), over several weeks, led to shiny dichroic crystalline materials being formed in cooler parts of the reaction vessel. An addition of small quantities of metals-Mo or Nb-was made with the aim of partially reducing their highly oxidized Mo(VI), W(VI) or Nb(V) species to corresponding Mo(V), W(V) and Nb(IV) centers, in order to form mixed valence solids. Sublimed crystals of four new compounds were investigated using a variety of techniques, with prime emphasis on the X-ray analysis, followed by spectroscopy (diffusion reflectance, IR, Raman and EPR), second harmonic generation (SHG), thermal analysis under N 2 and air atmosphere, and single crystals electrical conductivity studies. The results evidenced the formation of new complex solids of previously unknown compositions and structures. Three out of four compounds crystallized in non-centrosymmetric space groups and represent layered 2D polymeric puckered structures that being stacked on each other form 3D lattices. All new solids exhibit strong second-harmonic-generation (SHG effect; based on YAG 1064 nm tests with detection of 532 nm photons), and a rare photosalient effect when crystals physically move in the laser beam. Single crystals’ electrical conductivity of the four new synthesized compounds was measured, and the results showed their semiconductor behavior. Values of band gaps of these new solids were determined using diffusion reflectance spectroscopy in the visible region. Aspects of new solids’ practical usefulness are discussed.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Density Measurements of Molten LiF–BeF 2 and LiF–BeF 2 –LaF 3 Salt Mixtures by Neutron Radiography

The densities of eutectic (LiF) 2 –BeF 2 and mixtures of this salt (FLiBe) with LaF 3 were measured by dilatometry and by neutron attenuation from 673 K to 1,073 K. Because LaF 3 has a limited solubility in FLiBe, it was necessary to determine the amount of LaF 3 in solution before the density could be determined. The FLiBe density determination was favorably benchmarked against the literature data. A simple comparison was not available for the LaF 3 –FLiBe mixtures, so extrapolation of published data was necessary based on analysis using the Molten Salt Thermal Properties Database-Thermochemistry, or MSTDB-TC, developed by the US Department of Energy. Solubilities for LaF 3 in FLiBe ranged from 1 to 4 mol % over 673 to 1,073 K. The salt system was heated and cooled over 24 h to evaluate potential changes in composition and hysteresis during the measurement. Changes in the meniscus were observed, and these were included in the correction for density determinations. Salt surface tension may have led to supersaturation of LaF 3 in the salt because the solubility curve was nonlinear with respect to the inverse temperature, as would be expected for an ideal system. Surface tension measurements are currently underway to test this hypothesis.

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

Local and Global Sensitivity Analysis of a Reactive Transport Model Simulating Floodplain Redox Cycling

Reactive transport models (RTMs) are essential tools that simulate the coupling of advective, diffusive, and reactive processes in the subsurface, but their complexity makes them difficult to understand, develop and improve without accompanying statistical analyses. Although global sensitivity analysis (SA) can address these issues, the computational cost associated with most global SA techniques limits their use with RTMs. In this study, we apply distance-based generalized sensitivity analysis (DGSA), a novel and computationally efficient method of global SA, to a floodplain-scale RTM and compare DGSA results to those from local SA. Our test case focuses on the impact of 17 uncertain environmental parameters on spatially and temporally variable redox conditions within a floodplain aquifer. The input parameters considered include flow and diffusion rates, geochemical reaction rates, and the spatial distribution of sediment facies. Sensitivity was evaluated for three distinct components of the model response, encompassing both multidimensional and categorical output. Parameter rankings differ between local SA and DGSA, due to nonlinear effects of individual parameters and interaction effects between parameters. DGSA results show that fluid residence time, which is controlled by aquifer permeability, generally exerts a stronger control on redox conditions than do geochemical reaction rates. Sensitivity indices also demonstrate that sulfate reduction is key for establishing and maintaining reducing conditions throughout the aquifer. Furthermore, these results provide insights into the key drivers of heterogeneous redox processes within floodplain aquifers, as well as the main sources of uncertainty when modeling complex subsurface systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Generation of shear flows induced by AE / EPM in LHD plasma

The generation of shear flows (SFs) by Alfven Eigenmodes (AEs) and energetic particle modes (EPMs) have important effects on the operation of future nuclear fusion reactors, because SFs regulate the saturation of the AEs/EPMs, the transport of EPs and thermal plasma, as well as the formation of transport barriers among other consequences. The aim of this study is the analysis of SFs generation during the saturation phase of AEs and EPMs in LHD plasma. Experiments performed in the 23rd and 24th LHD experimental campaigns are dedicated to explore the destabilization of AEs/EPMs in discharges with different heating patterns, thermal plasma and magnetic field configurations. In particular, the shots 176490 and 179697 show the destabilization of MHD bursts and energetic-ion-driven resistive interchange modes (EIC), respectively. Charge exchange spectroscopy measurements in both discharges indicate that the generation of SFs by AE/EPM is uncorrelated with the perturbation induced by the neutral beam injector (NBI). Nonlinear simulations performed using the gyro-fluid code FAR3d show the generation of zonal structures, especially SFs, induced during the saturation phase of Toroidal Alfven Eigenmodes (TAEs) triggered in the MHD burst as well as by the 1/1 EIC in the bursting phase. The simulations indicate that SFs are caused by the radial electric fields powered by energy transfers from the unstable AE/EPM towards the thermal plasma. The strongest SFs are measured during the EIC bursting phase once the 1/1 EPM overlaps with nearby resonances at the plasma periphery. Likewise, the largest SFs during the MHD burst are observed once TAEs radially overlap in the inner-middle plasma region.

AE↗
DOE OSTI2021

Probing Intra-Halo Light with Galaxy Stacking in CIBER Images

We study the stellar halos of 0.2 ≲ $z$ ≲ 0.5 galaxies with stellar masses spanning $M$* ~ 10 10.5 to 10 12 $M$ ⊙ (approximately $L$* galaxies at this redshift) using imaging data from the Cosmic Infrared Background Experiment (CIBER). A previous CIBER fluctuation analysis suggested that intra-halo light (IHL) contributes a significant portion of the near-infrared extragalactic background light (EBL), the integrated emission from all sources throughout cosmic history. In this work, we carry out a stacking analysis with a sample of ~30,000 Sloan Digital Sky Survey (SDSS) photometric galaxies from CIBER images in two near-infrared bands (1.1 and 1.8 $μ$m) to directly probe the IHL associated with these galaxies. We stack galaxies in five sub-samples split by brightness and detect an extended galaxy profile beyond the instrument point-spread function (PSF) derived by stacking stars. We jointly fit a model for the inherent galaxy light profile plus large-scale one- and two-halo clustering to measure the extended galaxy IHL. We detect nonlinear one-halo clustering in the 1.8 $μ$m band at a level consistent with numerical simulations. Furthermore, by extrapolating the fraction of extended galaxy light, we measure to all galaxy mass scales, we find ~30%/15% of the total galaxy light budget from galaxies is at radius $r$ > 10/20 kpc, respectively. These results are new at near-infrared wavelengths at the $L$* mass scale and suggest that the IHL emission and one-halo clustering could have appreciable contributions to the amplitude of large-scale EBL background fluctuations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

TEM turbulence in simulation and experiment in quasisymmetric and degraded-symmetry configurations of the HSX stellarator

The Helically Symmetric Experiment (HSX) has demonstrated reduced neoclassical transport and flow damping with quasisymmetry, but the difference of turbulent transport between quasisymmetric and degraded-symmetry configurations has not been definitively assessed. Profile analysis presented here shows that anomalous transport at the mid-radius is larger with degraded quasisymmetry when temperature and density profiles are matched. An initial comparison (${T_\mathrm{i}/T_\mathrm{e}} = 1$) between the quasisymmetric and degraded symmetry configurations finds that trapped-electron mode (TEM) turbulence can support larger heat fluxes in the degraded symmetry configuration despite smaller linear growth rates. However, further simulations at the experimental temperature ratio ($T_\mathrm{i}/T_\mathrm{e} = 0.2$) do not reproduce this effect. Experimental measurements of heat transport and density fluctuation amplitude are compared to nonlinear gyrokinetic simulations at the experimental temperature ratio. The heat flux and density fluctuation amplitude in simulation reproduce a substantial dependence on the density gradient, and the simulated heat flux matches measurements within experimental uncertainties, indicating that TEM turbulence is responsible for the anomalous transport at the steep density gradient in HSX.

HSX↗
DOE OSTI2024

Measurement of stability diagrams in the IOTA ring at Fermilab

Nonlinear focusing elements can enhance the stability of particle beams in high-energy colliders through Landau Damping, by means of the tune spread which is introduced. Here we discuss an experiment at Fermilab's Integrable Optics Test Accelerator (IOTA) which investigates the influence of nonlinear focusing elements, such as octupoles, on the beam’s transverse stability. In this experiment, we employ an anti-damper, an active transverse feedback system, as a controlled mechanism to induce coherent beam instability. By utilizing the anti-damper we can examine the impact of a nonlinear focusing element on the beam's transverse stability. The stability diagram, a tool used to determine the system's stability, is measured using a recently demonstrated method at the LHC. The experiment at IOTA adds insight towards this stability diagram measurement method by supplying a reduced machine impedance to investigate the machine impedance’s effect on the stability diagram, as well as by aiming to map out the full stability diagram by using a large phase range of the anti-damper. From this experiment in IOTA, we present the first results of stability diagram analysis with varying octupole currents.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Impact of equilibrium radial electric field on energy loss process after pedestal collapse

An impact of the equilibrium radial electric field on energy loss processes after pedestal collapse is numerically investigated using the BOUT++ framework. Using linear stability analysis, the resistive ballooning mode is shown to be stabilized by the radial shear of the equilibrium radial electric field. On the other hand, the energy loss level after the pedestal collapse increases if the equilibrium radial electric field is taken into account. The spatio-temporal and phase diagram analyses show that the equilibrium radial electric field partially cancels the fluctuation-driven toroidally axisymmetric radial electric field and weakens the E × B shearing rate after pedestal collapse, weakening the turbulence suppression by vortex shearing. Finally, the equilibrium radial electric field therefore increases turbulence intensity in nonlinear cyclic oscillations among pressure gradient, E × B shearing rate, and turbulence intensity, which gives rise to subsequent bursts of turbulent transport and increases the energy loss level.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Landau-Ginzburg theory of charge density wave formation accompanying lattice and electronic long-range ordering

Here we propose an analytical Landau-Ginzburg (LG) theory of the charge density waves coupled with lattice and electronic long-range order parameters. Examples of long-range order include the electronic wave function of superconducting Cooper pairs, structural distortions, electric polarization, and magnetization. We formulate the LG free energy density as a power expansion with respect to the charge density and other long-range order parameters as well as their spatial gradients and biquadratic coupling terms. We introduced a biquadratic coupling between the charge density gradient and long-range order parameters as well as nonlinear higher gradients of the long-range order parameters. The biquadratic gradient coupling is critical to the appearance of different spatially modulated phases in charge-ordered ferroics and high-temperature superconductors. We derived the thermodynamic conditions for the stability of the spatially modulated phases, which are the intertwined spatial waves of charge density and lattice/electronic long-range order. The analytical expressions for the energies of different phases, corresponding order parameters, charge density waves amplitudes, and modulation periods obtained in this paper can be employed to guide the comprehensive physical explanation, deconvolution, and Bayesian analysis of experimental data on quantum materials ranging from charge-ordered ferroics to high-temperature superconductors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Synchrophasor Measurements-based Events Detection Using Deep Learning

Deep learning algorithms have been developed for phasor measurement units (PMUs) analysis aiming at providing grid operators to observe and react to significant real-time changes in the grid associated with multiple factors (e.g., power generation and load variations, different type of faults, and equipment mailfunction), or for offline post-event system diagnostics. In this study, a Long Short-Term Memory (LSTM)-based deep neural network (DNN) is adopted and evaluated to identify the most appropriate model configurations for event detection and longer-term anomalous pattern extraction. The proposed DNN model shows the potential on long-term predictions with the ability to capture nonlinear and nonstationary mixture complex patterns in PMU datasets. Real-world PMU in the WECC system were used for model development and validation.

deep learning↗
DOE OSTI2021

Bandwidth extension and conversion efficiency improvements beyond phase matching limitations using cavity-enhanced OPCPA

The conversion efficiency and phase matching bandwidth of ultrafast optical parametric amplification (OPA) are constrained by the dispersion and nonlinear coefficient of the employed crystal as well as pulse shaping effects. In our work we show that an enhancement cavity resonant with the pump seeded at the full repetition rate of the pump laser can automatically reshape the small-signal gain in optical parametric chirped-pulse amplification (OPCPA) to achieve close-to-optimal operation. This new method termed cavity-enhanced OPCPA or C-OPCPA significantly increases both the gain bandwidth and the conversion efficiency, in addition to boosting gain for high-repetition-rate amplification. The goal in C-OPCPA is to arrive at a condition of impedance matching at all temporal coordinates, such that, in the absence of linear losses, all the incident pump power is dissipated in the nonlinear loss element, i.e., converted to signal and idler. The use of a low finesse enhancement cavity resonant with a low average power (<1W) and a high repetition rate (78MHz) pump source is shown to achieve more than 50% conversion efficiency into signal and idler from the coupled pump in an optical parametric process, whereas an equivalent amount of pump power in a single-pass configuration leads to negligible conversion. Additionally, the gain bandwidth is extended by a factor of 3-4 beyond the phase-matching limit. Our empirical observations are corroborated by a numerical analysis of depletion optimizing the single-pass case, which assesses the underlying impedance matching that is responsible for the observed performance improvements.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Anderson acceleration stability in NDA-accelerated k-eigenvalue problems

Anderson acceleration (AA) has been used to improve the stability and convergence rate of multiphysics iterative methods for reactor analysis. Most applications studied assume a tightly converged solution for the different physics problems, and AA is usually applied to state variables like temperature, density, and heat generation rate. In this paper, we study the theoretical performance of AA in NDA-accelerated k-eigenvalue problems. The problems and algorithms studied are simplified from the coupled iteration scheme adopted by MPACT and many other high-fidelity whole-core reactor codes. Compared to previous analyses of AA for these iteration schemes, we study the case with a partially converged neutronics solution and possibly partially converged nonlinear diffusion acceleration (NDA)/coarse mesh finite difference (CMFD) solutions. We observe that the performance of the iteration scheme with AA is very sensitive to the initial guess and is affected by the partially converged CMFD solutions. When the NDA solution is fully converged, using AA cannot achieve the optimal convergence rate in large-sized problems. Conversely, if the NDA solution is partially converged, the iteration scheme with AA can diverge or converge extremely slowly. It is found that the loss of robustness for AA is due to the fact that it is applied to the iterative subspace of state variables rather than the fundamental unknowns of the governing equations. To improve the robustness, the scalar flux should also be considered in the implementation of AA. After considering the residuals of flux, we observe that the stability is regardless of the partial convergence of NDA solutions. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Robust Stabilization of Inverter-Based Resources Using Virtual Resistance-Based Control

This letter proposes a virtual resistance-based nonlinear control to stabilize and robustify the current layer of inverter-based resources, subject to the grid voltage disturbances, and the grid parameter uncertainties. A class of virtual resistances is proposed and analyzed using concepts from dissipative systems theory. Moreover, specific nonlinear virtual resistance-based controllers are derived, with their corresponding performance analytically bounded. Here, the theoretical and simulation results show that the proposed nonlinear virtual resistance-based controllers significantly reduce the L 2 gain of the closed-loop error system to grid voltage variation and parametric uncertainties. Significant improvements of the transient and steady-state current responses are also demonstrated over linear virtual resistance counterparts.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Chapter 4 - Recent Advances in Identification of Differential Equations from Noisy Data: IDENT Review

Differential equations and numerical methods are extensively used to model various real-world phenomena in science and engineering. With modern developments, we aim to find the underlying differential equation from a single observation of time-dependent data. If we assume that the differential equation is a linear combination of various linear and nonlinear differential terms, then the identification problem can be formulated as solving a linear system. The goal then reduces to finding the optimal coefficient vector that best represents the time derivative of the given data. We review some recent works on the identification of differential equations. We find some common themes for the improved accuracy: (i) The formulation of linear system with proper denoising is important, (ii) how to utilize sparsity and model selection to find the correct coefficient support needs careful attention, and (iii) there are ways to improve the coefficient recovery. We present an overview and analysis of recent developments on the topic.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Efficient data-driven regression for reduced-order modeling of spatial pattern formation

We present an efficient data-driven regression approach for constructing reduced-order models (ROMs) of reaction-diffusion systems exhibiting pattern formation. The ROMs are learned non-intrusively from available training data of physically accurate numerical simulations. The method can be applied to general nonlinear systems through the use of polynomial model form, while not requiring knowledge of the underlying physical model, governing equations, or numerical solvers. The process of learning ROMs is posed as a low-cost least-squares problem in a reduced-order subspace identified via Proper Orthogonal Decomposition (POD). Numerical experiments on classical pattern-forming systems–including the Schnakenberg and Mimura–Tsujikawa models–demonstrate that higher-order surrogate models significantly improve prediction accuracy while maintaining low computational cost. The proposed method provides a flexible, non-intrusive model reduction framework, well suited for the analysis of complex spatio-temporal pattern formation phenomena.

Data-driven modeling↗
DOE OSTI2022

Investigation of core impurity transport in DIII-D diverted negative triangularity plasmas

Abstract Tokamak operation at negative triangularity has been shown to offer high energy confinement without the typical disadvantages of edge pedestals (Marinoni et al 2021 Nucl. Fusion 61 116010). In this paper, we examine impurity transport in DIII-D diverted negative triangularity experiments. Analysis of charge exchange recombination spectroscopy reveals flat or hollow carbon density profiles in the core, and impurity confinement times consistently shorter than energy confinement times. Bayesian inferences of impurity transport coefficients based on laser blow-off injections and forward modeling via the Aurora package (Sciortino et al 2021 Plasma Phys. Control. Fusion 63 112001) show core cross-field diffusion to be higher in L-mode than in H-mode. Impurity profile shapes remain flat or hollow in all cases. Inferred radial profiles of diffusion and convection are compared to neoclassical, quasilinear gyrofluid, and nonlinear gyrokinetic simulations. Heat transport is observed to be better captured by reduced turbulence models with respect to particle transport. State-of-the-art gyrokinetic modeling compares favorably with measurements across multiple transport channels. Overall, these results suggest that diverted negative triangularity discharges may offer a path to a highly-radiative L-mode scenario with high core performance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

An In-Depth Analysis of the Transformation of Tin Foil Anodes during Electrochemical Cycling in Lithium-Ion Batteries

Tin foils have an impressive lithium-storage capacity more than triple that of graphite anodes, and their adoption could facilitate a drastic improvement in battery energy density. However, implementation of a dense foil electrode architecture represents a significant departure from the standard blade-cast geometry with a distinct electrochemical environment, and this has led to confusion with regards to the first cycle efficiency of the system. In this work, we investigate the unique behavior of a tin active material in a foil architecture to understand its performance as an anode. We find shallow cycling of the foil results in an irreversible formation (< 40%) due to diffusional trapping, but intermediate and complete utilization allows for a remarkably reversible formation reaction (> 90%). This striking nonlinearity stems from an in situ transformation from bulk metal to porous electrode that occurs during formation cycles and defines electrode-level lithium-transport on subsequent cycles. As a result, an alternative cycling procedure for assessing the stability of foils is proposed to account for this chemomechanical effect.

25 ENERGY STORAGE↗
DOE OSTI2025

Multi-scale, Multi-disciplinary, and Multi-agent Explainable AI with Koopman-Undergirded Learning, Prediction, and Analysis (M3EA KULPA) (Project Closeout Report)

The goal of this project was to develop and use domain-aware machine learning formulations, based on the Koopman Operator (KO), for modelling multi-scale, multi-disciplinary (e.g., multi-physics), and/or multi-agent systems. The project developed these formulations for the following cases: • Systems with dynamics at two separate time scales, • Systems with a bi-level hierarchical control structure, • Systems with bi-level hierarchical control and dynamics at two separate time scales (the lower level controls operating at the faster time scale), and • Systems with n separate but interacting agents/disciplines (with/without control, respectively); the controls for each agent could include bi-level hierarchical control and dynamics at two separate time scales as described above. The project then defined a set of dynamical systems consisting of different nonlinear oscillators that could be used to test these different formulations and then subsequently learned the KO models for those systems. With the KO models, we were able to do the following: • Quantify system stability, including both long-term and transient behavior, • Quantify the effects of feedbacks between the different time scales and agents/disciplines in terms of those feedbacks’ effects on system stability, • Replace a standard Proportional-Integral (PI) control in the hierarchical control structure with a KO-based Linear-Quadratic Regular (LQR), a form of optimal control, • Calculate optimal supervisory control policies a) with and without time scale separated dynamics at the lower level control levels and b) with both PI and KO-based LQR lower level control policies, and • Calculate dynamic Nash equilibria for multi-agent systems where each agent makes its own control decisions.

97 MATHEMATICS AND COMPUTING↗