Search NASA⌕ Search

SEARCH · Search NASA

Results for “Nonlinear approximation”

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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Solitary waves in the nonlinear Dirac equation in the presence of external driving forces

In this paper, we consider the nonlinear Dirac (NLD) equation in (1 + 1) dimensions with scalar–scalar self interaction g 2 /κ + 1 (Ψ¯Ψ) κ + 1 in the presence of external forces as well as damping of the form f(x) - iμγ 0 Ψ, where both f and Ψ are two-component spinors. We develop an approximate variational approach using collective coordinates (CC) for studying the time dependent response of the solitary waves to these external forces. This approach predicts intrinsic oscillations of the solitary waves, i.e. the amplitude, width and phase all oscillate with the same frequency. The translational motion is also affected, because the soliton position oscillates around a mean trajectory. For κ = 1 we solve explicitly the CC equations of the variational approximation for slow moving solitary waves in a constant external force without damping and find reasonable agreement with solving numerically the CC equations. Finally, we then compare the results of the variational approximation with no damping with numerical simulations of the NLD equation for κ = 1, when the components of the external force are of the form f j = r j exp(–iΚx) and again find agreement if we take into account a certain linear excitation with specific wavenumber that is excited together with the intrinsic oscillations such that the momentum in a transformed NLD equation is conserved.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Uncovering extreme nonlinear dynamics in solids through time-domain field analysis

Time-domain analysis of harmonic fields with sub-cycle resolution is now experimentally viable due to the emergence of sensitive, on-chip techniques for petahertz-scale optical-field sampling. We demonstrate how such a time-domain, field-resolved analysis uncovers the extreme nonlinear electron dynamics responsible for high-harmonic generation within solids. Time-dependent density functional theory was used to simulate harmonic generation from a solid-state band-gap system driven by near- to mid-infrared waveforms. Particular attention was paid to regimes where both intraband and interband emission mechanisms play a critical role in shaping the nonlinear response. We show that a time-domain analysis of the harmonic radiation fields identifies the interplay between intra- and interband dynamical processes underlying the nonlinear light generation. With further analysis, we show that changes to the dominant emission regime can occur after only slight changes to the peak driving intensity and central driving wavelength. Furthermore, time-domain analysis of harmonic fields also reveals, for the first time, the possibility of rapid changes in the dominant emission mechanism within the temporal window of the driving pulse envelope. Finally, we examine the experimental viability of performing time-domain analysis of harmonic fields with sub-cycle resolution using realistic parameters.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Average nonlinear dynamics of particles in gravitational pulses: Effective Hamiltonian, secular acceleration, and gravitational susceptibility

Particles interacting with a prescribed quasimonochromatic gravitational wave (GW) exhibit secular (average) nonlinear dynamics that can be described by Hamilton’s equations. We derive the Hamiltonian of this “ponderomotive” dynamics to the second order in the GW amplitude for a general background metric. For the special case of vacuum GWs, we show that our Hamiltonian is equivalent to that of a free particle in an effective metric, which we calculate explicitly. We also show that already a linear GW pulse displaces a particle from its unperturbed trajectory by a finite distance that is independent of the GW phase and proportional to the integral of the pulse intensity.This effect is independent from the nonlinear memory effects that has been known. We calculate the particle displacement analytically and show that our result is in agreement with numerical simulations. We also show how the Hamiltonian of the nonlinear averaged dynamics naturally leads to the concept of the linear gravitational susceptibility of a particle gas with an arbitrary phase-space distribution. We calculate this susceptibility explicitly to apply it, in a follow-up paper, toward studying self-consistent GWs in inhomogeneous media within the geometrical-optics approximation

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Uncertainty propagation and sensitivity analysis for constrained optimization of nuclear waste vitrification

Abstract The vitrification of high‐level waste (HLW) by heating a mixture of glass‐forming chemicals (GFCs) with the waste can be improved using a constrained optimization problem. This study explores how different uncertainty propagation (UP) methods implemented with the optimization process can affect the glass formulation of nuclear waste glasses. UP is the effort of propagating uncertain inputs through a system to understand and quantify output distributions. Uncertainty intervals are crafted from output distributions to inform the optimization algorithm. UP is often implemented with Monte Carlo (MC) sampling for large nonlinear systems, which can be difficult to implement within a constrained optimization algorithm that requires derivative information. Other UP methods often used for optimization under uncertainty (OUU) can be designed to work within an established constrained optimization framework. Methods of UP are evaluated in this study including iterative sampling approaches, first‐order approximations, and surrogate modeling with machine learning (ML). A method of dimensional reduction based on global sensitivity analysis is introduced to support the UP methods for the large dimensionality of the problem. Analytical UP methods able to achieve similar optimums 10 times faster than the baseline MC approach, and produce 93.9% similar output distributions are reported.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2016

EMAPS: An Efficient Multiscale Approach to Plasma Systems with Non-MHD Scale Effects

Global interactions of energetic ions with magnetoplasmas and neutral gases lie at the core of many space and laboratory plasma phenomena ranging from solar wind entry into and transport within planetary magnetospheres and exospheres to fast-ion driven instabilities in fusion devices to astrophysics-in-lab experiments. The ability of computational models to properly account for physical effects that underlie such interactions, namely ion kinetic, ion cyclotron, Hall, collisional and ionization processes is important for the success and planning of experimental research in plasma physics. Understanding the physics of energetic ions, in particular their nonlinear resonance interactions with Alfvén waves, is central to improving the heating performance of magnetically confined plasmas for future energy generation. Fluid models are not adequate for high-beta plasmas as they cannot fully capture ion kinetic and cyclotron physics (e.g., ion behavior in the presence of magnetic nulls, shock structures, plasma interpenetration, etc.). Recent results from global reconnection simulations show that even in a MHD-like regime there may be significant differences between kinetic and MHD simulations. Therefore, kinetic modeling becomes essential for meeting modern day challenges in plasma physics. The hybrid approximation is an intermediate approximation between the fluid and fully kinetic approximations. It eliminates light waves, removes the electron inertial temporal and spatial scales from the problem and enables full-orbit ion kinetics. As a result, hybrid codes have become effective tools for exploring ion-scale driven phenomena associated with ion beams, shocks, reconnection and turbulence that control the large-scale behavior of laboratory and space magnetoplasmas. A number of numerical issues, however, make three-dimensional (3D) large-scale hybrid simulations of inhomogeneous magnetized plasmas prohibitively expensive or even impossible. To resolve these difficulties we have developed a novel Event-driven Multiscale Asynchronous Parallel Simulation (EMAPS) technology that replaces time stepping with self-adaptive update events. Local calculations are carried out only on an “as needed basis”. EMAPS (i) guarantees accurate and stable processing of physical variables in time accurate simulations, and (ii) eliminates unnecessary computation. Applying EMAPS to the hybrid model has resulted in the development of a unique parallel code, dimension-independent (compile-time-configurable) HYPERS (Hybrid Parallel Event-Resolved Simulator) that scales to hundreds of thousands of parallel processors. HYPERS advances electromagnetic fields and particles asynchronously on time scales determined by local physical laws and mesh properties. To achieve high computational accuracy in complex device geometries, HYPERS employs high-fidelity Cartesian grids with masked conductive cells. The HYPERS model includes multiple ion species, energy and momentum conserving ion-ion collisions, and provides a number of approximations for plasma resistivity and vacuum regions. Both local and periodic boundary conditions are allowed. The HYPERS solver preserves zero divergence of magnetic field. The project has demonstrated HYPERS capabilities on a number of applications of interest to fusion and astrophysical plasma physics applications listed below. 1. Theta-pinch formation of FRCs The formation, spontaneous spin-up, and stability of theta-pinch formed field-reversed configurations have been studied self-consistently in 3D. The end-to-end hybrid simulations reveal poloidal profiles of implosion-driven fast toroidal plasma rotation and demonstrate three discharge regimes as a function of experimental parameters: the decaying stable configuration, the tilt unstable configuration, and the nonlinear evolution of a fast growing tearing mode. 2. FRC collisions with magnetic mirrors Interactions of fast plasma streams and objects with magnetic obstacles (dipoles, mirrors, etc) lie at the core of many space and laboratory plasma phenomena ranging from magnetoshells and solar wind interactions with planetary magnetospheres to compact fusion plasmas. HYPERS simulations are compared with data from the MSX experiment (LANL) that focuses on the physics of magnetized collisionless shocks through the acceleration and subsequent stagnation of FRC plasmoids against a strong magnetic mirrors and flux-conserving boundaries. 3. Exploding magnetoplasmas Results from hybrid simulations of two experiments at the LAPD and Nevada Terawatt Facility are discussed where short-pulse lasers are used to ablate solid targets to produce plasmas that expand across external magnetic fields. The first simulation recreates flutelike density striations observed at the leading edge of a carbon plasma and predicts an early destruction of the magnetic cavity in agreement with experimental evidence. In the second simulation a polyethylene target is ablated into a mixture of protons and carbon ions. A mechanism is demonstrated that allows protons to penetrate the magnetic field in the form of a collimated flow. The results are compared to experimental data and single-fluid MHD simulations. The EMAPS framework has the potential for wide application in many other engineering and scientific fields, such as climate models, biological systems, electronic devices, seismic events, oil reservation simulators that all involve advancing solutions of partial differential equations in time where the rate of activity can be adapted widely over the spatial domain depending on locally space/time phenomena (“events”).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Unexpected Scaling of Peak Acceleration for a Yielding Mass-Spring System Subjected to a Triangular Base Acceleration Pulse

The use of bounding scenarios is a common practice that greatly simplifies the design and qualification of structures. However, this approach implicitly assumes that the quantities of interest increase monotonically with the input to the structure, which is not necessarily true for nonlinear structures. This paper surveys the literature for observations of nonmonotonic behavior of nonlinear systems and finds such observations in both the earthquake engineering and applied mechanics literature. Numerical simulations of a single degree-of-freedom mass-spring system with an elastic–plastic spring subjected to a triangular base acceleration pulse are then presented, and it is shown that the relative acceleration of this system scales nonmonotonically with the input magnitude in some cases. The equation of motion for this system is solved symbolically and an approximate expression for the relative acceleration is developed, which qualitatively agrees with the nonmonotonic behavior seen in the numerical results. The nonmonotonicity is investigated and found to be a result of dynamics excited by the discontinuous derivative of the base acceleration pulse, the magnitude of which scales nonmonotonically with the input magnitude due to the fact that the first yield of the spring occurs earlier as the input magnitude is increased. The relevance of this finding within the context of defining bounding scenarios is discussed, and it is recommended that modeling be used to perform a survey of the full range of possible inputs prior to defining bounding scenarios.

42 ENGINEERING↗
DOE OSTI2024

A reduced-order model for nonlinear radiative transfer problems based on moment equations and POD-Petrov-Galerkin projection of the normalized Boltzmann transport equation

A data-driven projection-based reduced-order model (ROM) for nonlinear thermal radiative transfer (TRT) problems is presented. The TRT ROM is formulated by (i) a hierarchy of low-order quasidiffusion (aka variable Eddington factor) equations for moments of the radiation intensity and (ii) the normalized Boltzmann transport equation (BTE). The multilevel system of moment equations is derived by projection of the BTE onto a sequence of subspaces which represent elements of the phase space of the problem. Exact closure for the moment equations is provided by the Eddington tensor. A Petrov-Galerkin (PG) projection of the normalized BTE is formulated using a proper orthogonal decomposition (POD) basis representing the normalized radiation intensity over the whole phase space and time. The Eddington tensor linearly depends on the solution of the normalized BTE. By linear superposition of the POD basis functions, a low-rank expansion of the Eddington tensor is constructed with coefficients defined by the PG projected normalized BTE. The material energy balance (MEB) equation is coupled with the effective gray low-order equations which exist on the same dimensional scale as the MEB equation. The resulting TRT ROM is structure and asymptotic preserving. A detailed analysis of the ROM is performed on the classical Fleck-Cummings (F-C) TRT multigroup test problem in 2D geometry. Numerical results are presented to demonstrate the ROM's effectiveness in the simulation of radiation wave phenomena. Importantly, the ROM is shown to produce solutions with sufficiently high accuracy while using low-rank approximation of the normalized BTE solution. Essential physical characteristics of supersonic radiation wave are preserved in the ROM solutions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

The data-driven future of high-energy-density physics

High-energy-density physics is the field of physics concerned with studying matter at extremely high temperatures and densities. Such conditions produce highly nonlinear plasmas, in which several phenomena that can normally be treated independently of one another become strongly coupled. The study of these plasmas is important for our understanding of astrophysics, nuclear fusion and fundamental physics—however, the nonlinearities and strong couplings present in these extreme physical systems makes them very difficult to understand theoretically or to optimize experimentally. Here we argue that machine learning models and data-driven methods are in the process of reshaping our exploration of these extreme systems that have hitherto proved far too nonlinear for human researchers. Furthermore, from a fundamental perspective, our understanding can be improved by the way in which machine learning models can rapidly discover complex interactions in large datasets. From a practical point of view, the newest generation of extreme physics facilities can perform experiments multiple times a second (as opposed to approximately daily), thus moving away from human-based control towards automatic control based on real-time interpretation of diagnostic data and updates of the physics model. To make the most of these emerging opportunities, we suggest proposals for the community in terms of research design, training, best practice and support for synthetic diagnostics and data analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Heterotrophic respiration by soil microbes in a changing climate

Soil microbes strongly influence the soil organic carbon (SOC) pool, which globally stores ~2,000 PgC. Specifically, the balance between microbial heterotrophic respiration (R H ), which degrades SOC, and plant–microbe interactions that stabilize SOC determines whether terrestrial ecosystems are a net source or sink of CO 2 to the atmosphere. Here, in this Review, we evaluate how climate change alters these competing processes. R H is approximately half of total soil respiration, at ~50 PgC yr −1 , with 70% occurring in topsoils. Warming accelerates microbial metabolism, with a 10 °C temperature increase estimated to raise R H by ~50%, an effect that is particularly strong in Arctic soils. Warming also reduces soil moisture, further modulating R H , which responds nonlinearly to soil moisture, being limited by saturation and desiccation and meeting a maximum at intermediate levels. Consequently, R H is highly sensitive to future precipitation changes and drought. However, soil management strategies could enhance SOC stocks and persistence under climate change. Bacterial and fungal inoculants can promote SOC production and stabilization, while deep-rooting plants increase SOC inputs to deeper layers that experience lower R H . Agricultural practices and biochar amendments can also enhance SOC and reduce R H . Expanding field trials across regions, climates and soil types would improve empirical understanding of these responses and support better representation of R H in predictive models, enabling more accurate assessments of climate impacts on SOC storage.

Jansson, Janet K. [Pacific Northwest National Labo↗
DOE OSTI2017

Kinetic fluid moments closure for a magnetized plasma with collisions

A novel method aimed at a kinetic moments closure for a magnetized plasma with collisions is proposed. The velocity distribution function for each species is expanded in 8 Gaussian Radial Basis Functions (GRBFs) which are essentially shifted Maxwellians at eight representative 3D-velocity points of drift. The vector of 8 fluid moments (for particle density, 3 particle fluxes, total energy density, and 3 energy fluxes) has an 8x8 analytic linear matrix relation to the vector of 8 GRBF density weighs in 3D-real space. The 8 fluid moments with sources for each species are advanced in time while the 8 GRBF weighs are determined from the 8x8 inverse matrix. The two closure moments (for the stress tensor and the energy weighted stress tensor) are linearly determined from the GRBF weighs. Most importantly the velocity moments of the nonlinear Coulomb Fokker-Planck collision operator [Rosenbluth et al, Phys. Rev 107, 1957] are evaluated from the GRBF weights. Generalization from 8 to 12, 16, 20 .., in an energy weighted moment hierarchy is straightforward. The electric field follows from a generalized vorticity (quasi-neutral charge conservation) equation. A strong drift ordering approximation can be applied to eliminate any spuriously unstable high frequency cyclotron motions. A novel weak drift ordering two-time step scheme avoids the vorticity equation by following ion cyclotron motion in time to get the electric field with ion gyroaveraging. Inclusion of low-beta magnetic perturbations is straightforward.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Nonlinear bioheat transfer models and multi-objective numerical optimization of the cryosurgery operations

Numerical computation of the three dimensional problem of the freezing interface propagation during the cryosurgery coupled with the multi-objective optimization methods is used in order to improve the efficiency and safety of the cryosurgery operations performing. Prostate cancer treatment and cutaneous cryosurgery are considered. The heat transfer in soft tissue during the thermal exposure to low temperature is described by the Pennes bioheat model and is coupled with an enthalpy method for blurred phase change computations. The finite volume method combined with the control volume approximation of the heat fluxes is applied for the cryosurgery numerical modeling on the tumor tissue of a quite arbitrary shape. The flux relaxation approach is used for the stability improvement of the explicit finite difference schemes. The method of the additional heating elements mounting is studied as an approach to control the cellular necrosis front propagation. Whereas the undestucted tumor tissue and destucted healthy tissue volumes are considered as objective functions, the locations of additional heating elements in cutaneous cryosurgery and cryotips in prostate cancer cryotreatment are considered as objective variables in multi-objective problem. The quasi-gradient method is proposed for the searching of the Pareto front segments as the multi-objective optimization problem solutions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Density functional theory investigation of the LiIn 1- x Ga x Se 2 solid solution

Here, the electronic structure and optical properties of the LiIn 1-x Ga x Se 2 (x=0, 0.25, 0.5, 0.75, 1) solid solution were studied by density functional theory (DFT) with pure functionals. The exchange-correlation is treated within the local density approximation (LDA) and generalized-gradient approximation (GGA). The electronic structures for each respective compound are discussed in detail. Calculations reveal that gallium incorporation can be used to tune the optical-electrical properties of the solid solution and correlates with the lattice parameter. The band gap trend of the LiIn 1-x Ga x Se 2 system follows a nonlinear behavior between the LiInSe 2 and LiGaSe 2 ternary boundaries. The bowing parameter is estimated to be on the order of 0.1- 0.3 eV at the point. Low-temperature optical absorption revealed a 30% change in the temperature dependence of the band gap for the intermediate compound LiIn 0.6 Ga 0.4 Se 2 compared to ternary boundaries and suggests the heat capacity to be another control element through strain.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Elastic Scattering of Ultrashort Laser Pulses on Noble Gas Atoms

The work deals with a theoretical investigation of elastic scattering of ultrashort electromagnetic pulses (USP) on the noble gas atoms in the far ultra-violet and soft X-ray regions. The investigation is performed within the framework of a dipole approximation and the first-order quantum-mechanics perturbation theory. Major attention is given to the study of the dependence of the total probability of scattering within the pulse time on the USP duration at different carrying frequency values. It is shown that at certain values of the problem parameters this dependence is nonlinear. The calculations are performed for the corrected Gaussian pulses and wavelet pulses without the carrying frequency.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Enhancement of the sensitivity of the light-induced drift effect to interatomic interaction potentials using a mixture of two buffer gases

Based on five different nonempirical (calculated ab initio) interaction potentials for pairs of colliding Li – Ne particles and three different interaction potentials for pairs of colliding Li – Ar particles, we have theoretically investigated the spectral features of the light-induced drift (LID) rate of Li atoms in the buffer Ne gas in a mixture of buffer Ne + Ar gases. The calculations of LID of Li atoms in the buffer Ne gas for two interaction potentials predict anomalous LID and, as a result, strong sensitivity of the spectral shape of the LID lines to the differences in these potentials. For three other potentials (out of the five in question), the shape of the LID line of Li atoms in Ne is insensitive to the shape of the potential, since calculations with these potentials predict a normal LID effect. In this case, as it turned out, by adding a small fraction (approximately 10 %) of Ar to Ne, one can go from normal LID to anomalous LID and thereby radically increase the sensitivity of the LID line shape of Li atoms to the difference in these interaction potentials. The results obtained enable high-precision testing of the interatomic interaction potentials in experiments on anomalous LID. (nonlinear optical phenomena)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Growth mechanism and optical properties of Ti thin films deposited onto fluorine-doped tin oxide glass substrate

In this work, a detailed study of the influence of the thickness on the morphological and optical properties of titanium (Ti) thin films deposited onto rough fluorine-doped tin oxide glass by d.c. magnetron sputtering is carried out. The films were characterized by several methods for composition, crystallinity, morphology, and optical properties. Regardless of the deposition time, all the studied Ti films of 400, 1500, 2000, and 2500 nm in thickness were single crystalline in the α-Ti phase and also very similar to each other with respect to composition. Using the atomic force microscopy (AFM) technique, the authors analyzed the roughness evolution of the Ti films characteristics as a function of the film thickness. By applying the dynamic scaling theory to the AFM images, a steady growth roughness exponent α = 0.72 ± 0.02 and a dynamic growth roughness exponent β = 0.22 ± 0.02 were determined. The value of α and β are consistent with nonlinear growth model incorporating random deposition with surface diffusion. Finally, measuring the reflection spectra of the samples by a spectrophotometer in the spectral range of 300–1100 nm allowed us to investigate the optical properties. The authors observed the increments of the reflection of Ti films with thickness, which by employing the effective medium approximation theory showed an increase in thickness followed by an increase in the volume fraction of metal.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Generative Vulnerability Assessment for Cyber-Physical Systems

Cyber-physical systems (CPS) are highly susceptible to malicious attacks due to their complex dynamics and interconnectivity. A comprehensive understanding of their vulnerabilities is essential for designing effective resilience measures. This paper presents a data-driven attack generative system for evaluating the vulnerability of CPS. The proposed approach formulates the vulnerability assessment problem as determining the feasibility of a specific attack set based on two boundary functions that represent the effectiveness and stealthiness of attacks. The attack generative model is trained using a custom loss function, with two universal approximators designed to learn the effectiveness and stealthiness functions simultaneously. Theoretical results for successful generation and asymptotic convergence of the resulting training algorithm are given. As a result, the proposed approach is evaluated via numerical simulation of an IEEE 14-bus system and gas pipeline systems, demonstrating its viability in learning how to attack nonlinear CPS and identify potential vulnerabilities.

Computer systems organization↗