Search NASA⌕ Search

SEARCH · Search NASA

Results for “nonlinear propagation”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Oblique instability of quasi-parallel whistler waves in the presence of cold and warm electron populations

Whistler waves propagating nearly parallel to the ambient magnetic field experience a nonlinear instability due to transverse currents when the background plasma has a population of sufficiently low energy electrons. Intriguingly, this nonlinear process may generate oblique electrostatic waves, including whistlers near the resonance cone with properties resembling oblique chorus waves in the Earth’s magnetosphere. Focusing on the generation of oblique whistlers, earlier analysis of the instability is extended here to the case where low-energy background plasma consists of both a “cold” population with energy of a few eV and a “warm” electron component with energy of the order of 100 eV. This is motivated by spacecraft observations in the Earth’s magnetosphere where oblique chorus waves were shown to interact resonantly with the warm electrons. The main new results are: 1) the instability producing oblique electrostatic waves is sensitive to the shape of the electron distribution at low energies. In the whistler range of frequencies, two distinct peaks in the growth rate are typically present for the model considered: a peak associated with the warm electron population at relatively low wavenumbers and a peak associated with the cold electron population at relatively high wavenumbers; 2) overall, the instability producing oblique whistler waves near the resonance cone persists (with a reduced growth rate) even in the cases where the temperature of the cold population is relatively high, including cases where cold population is absent and only the warm population is included; 3) particle-in-cell simulations show that the instability leads to heating of the background plasma and formation of characteristic plateau and beam features in the parallel electron distribution function in the range of energies resonant with the instability. The plateau/beam features have been previously detected in spacecraft observations of oblique chorus waves. However, they have been attributed to external sources and have been proposed to be the mechanism generating oblique chorus. In the present scenario, the causality link is reversed and the instability generating oblique whistler waves is shown to be a possible mechanism for formation of the plateau and beam features.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Measuring the growth of Alfvén wave parametric decay instability using counter-propagating waves: Theory and simulations

The parametric decay instability (PDI) of Alfvén waves—where a pump Alfvén wave decays into a backward-propagating child Alfvén wave and a forward ion acoustic wave—is a fundamental nonlinear wave-wave interaction and holds significant implications for space and laboratory plasmas. However, to date there has been no direct experimental measurement of PDI. Here, we propose a novel and experimentally viable scheme to quantify the growth of Alfvén wave PDI on a linear device using a large pump Alfvén wave and a small counter-propagating seed Alfvén wave, with the seed-wave frequency tuned to match the backward Alfvén wave generated by standard PDI. Using hybrid simulations, we show that energy transfer from the pump to the seed reduces the latter's spatial damping. By comparing seed-wave amplitudes with and without the pump wave, this damping reduction can be used as a direct and reliable proxy for PDI growth. The method is validated in our simulations across a range of plasma and wave parameters and agrees well with theoretical predictions. Notably, the scheme exhibits no threshold for PDI excitation and is, in principle, readily implementable under current laboratory conditions. Finally, this scheme is a critical step toward solving the challenge of experimentally accessing Alfvén wave PDI and provides an elegant method that may be used to validate fundamental theories of parametric instabilities in controlled laboratory settings.

heliosphere↗
DOE OSTI2022

Enhancing high-fidelity nonlinear solver with reduced order model

Abstract We propose the use of reduced order modeling (ROM) to reduce the computational cost and improve the convergence rate of nonlinear solvers of full order models (FOM) for solving partial differential equations. In this study, a novel ROM-assisted approach is developed to improve the computational efficiency of FOM nonlinear solvers by using ROM’s prediction as an initial guess. We hypothesize that the nonlinear solver will take fewer steps to the converged solutions with an initial guess that is closer to the real solutions. To evaluate our approach, four physical problems with varying degrees of nonlinearity in flow and mechanics have been tested: Richards’ equation of water flow in heterogeneous porous media, a contact problem in a hyperelastic material, two-phase flow in layered porous media, and fracture propagation in a homogeneous material. Overall, our approach maintains the FOM’s accuracy while speeding up nonlinear solver by 18–73% (through suitable ROM-assisted FOMs). More importantly, the proximity of ROM’s prediction to the solution space leads to the improved convergence of FOMs that would have otherwise diverged with default initial guesses. We demonstrate that the ROM’s accuracy can impact the computational efficiency with more accurate ROM solutions, resulting in a better cost reduction. We also illustrate that this approach could be used in many FOM discretizations (e.g., finite volume, finite element, or a combination of those). Since our ROMs are data-driven and non-intrusive, the proposed procedure can easily lend itself to any nonlinear physics-based problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

The Role of Stratigraphy and Loading History in Generating Complex Compaction Bands in Idealized Field-Scale Settings

The Buckskin Gulch locality in Utah is a landmark example of compaction localization. The outcrop of this locality involves distinct stratigraphic heterogeneity and was exposed to complex loading history. It features multiple sets of deformation bands with different kinematics and orientation. Similar formations were seen in the Valley of Fire, Nevada, and the Orange quarry, France, among other localities. The formation of such complex structures, their propagation mechanisms, and frequency is affected by numerous local and ambient factors whose impacts are not yet fully understood. The simulation of the above-mentioned localities is not feasible because of the limited amount of available information. This work, instead, investigates from a geomechanics standpoint how the interplay among material nonlinearity, outcrop stratigraphy, and loading history interconnects with specific spatiotemporal patterns of compaction band propagation. Our study shows that the system stratigraphy can be responsible for the emergence of coexisting compaction bands with different inclination and kinematics. Specifically, we show that stiffness contrasts induce nonlocal stress changes which may favor the initiation of secondary structures with different compaction localization characteristics. Furthermore, systems of inclined compaction bands induced by burial increase display secondary, noncontemporaneous sets of vertical compaction bands under the effects of postburial shortening. Our results indicate that stages of intermediate burial decrease prior to tectonic shortening can promote the formation of such complex systems. Despite the simplifications involved in our analyses, these findings show how geomechanics computations complement field observations and could provide a mechanics-based validation of site-specific reconstruction hypothesis.

58 GEOSCIENCES↗
DOE OSTI2022

Design of a supervisory control system for autonomous operation of advanced reactors

Advanced reactors to be deployed in the coming decades will face deregulated energy markets, and may adopt flexible operation to boost profitability. To aid in the transition from baseload to flexible operation paradigm, autonomous operation is sought. This work focuses on the control aspect of autonomous operation. Specifically, a hierarchical control system is designed to support constraint enforcement during routine operational transients. Within the system, data-driven modeling, physics-based state observation, and classical control algorithms are integrated to provide an adaptable and robust solution. A 320 MW Fluoride-cooled High-temperature Pebble-bed Reactor is the design basis for demonstrating the proposed control system. The hierarchical control system consists of a supervisory layer and low-level layer. The supervisory layer receives requests to change the system's operating conditions (e.g., the current reactor power to meet a load -follow), and accepts or rejects them based on constraints that have been assigned. Constraints are issued to keep the plant within an optimal operating region. The low-level layer interfaces with the actuators of the system to fulfill requested changes, while maintaining tracking and regulation duties. Further, to accept requests at the supervisory layer, the Reference Governor algorithm was adopted. To model the dynamics of the reactor, a system identification algorithm, Dynamic Mode Decomposition, was utilized. To estimate the evolution of process variables that cannot be directly measured (e.g., the propagation of delayed neutron precursors), the Unscented Kalman Filter, incorporating a nonlinear model of nuclear dynamics, was adopted. The composition of these algorithms led to a numerical demonstration of constraint enforcement during a 40% power drop transient (at a rate of 5 %/min). Uncontrolled secondary-side temperatures were successfully constrained. Adaptability of the proposed system was demonstrated by modifying the constraint values, and enforcing them during the transient. Robustness was also demonstrated by enforcing constraints under noisy environments.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

A resistive MHD model and simulation on plasma flow evolution in the presence of resonant magnetic perturbation in a tokamak

Nonaxisymmetric magnetic fields such as the intrinsic error field and the externally applied resonant magnetic perturbation (RMP) in a tokamak are known to influence the plasma momentum transport and flow evolution through plasma response, which itself strongly depends on the plasma flow as well. The nonlinear interaction between plasma response and flow has been previously modeled in the conventional error field theory with the “no-slip” condition, which has been recently extended to allow the “free-slip” condition. In this work, we further target this specific process and numerically simulate the nonlinear plasma response and flow evolution in the presence of a single-helicity RMP in a circular-shaped model tokamak configuration, based on the full resistive MHD model in the initial-value code NIMROD. Time evolution of the parallel (to k) flow or “slip frequency” profile and its asymptotic steady state obtained from the NIMROD simulations are compared with both conventional and extended nonlinear response theories. Here, k is the wave vector of the propagating island. Good agreement with the extended theory with free-slip condition has been achieved for the parallel flow profile evolution in response to RMP in all resistive regimes, whereas the difference from the conventional theory with the no-slip condition tends to diminish as the plasma resistivity approaches zero.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Synthetically non-Hermitian nonlinear wave-like behavior in a topological mechanical metamaterial

Topological mechanical metamaterials have enabled new ways to control stress and deformation propagation. Exemplified by Maxwell lattices, they have been studied extensively using a linearized formalism. Herein, we study a two-dimensional topological Maxwell lattice by exploring its large deformation quasi-static response using geometric numerical simulations and experiments. We observe spatial nonlinear wave-like phenomena such as harmonic generation, localized domain switching, amplification-enhanced frequency conversion, and solitary waves. We further map our linearized, homogenized system to a non-Hermitian, nonreciprocal, one-dimensional wave equation, revealing an equivalence between the deformation fields of two-dimensional topological Maxwell lattices and nonlinear dynamical phenomena in one-dimensional active systems. Our study opens a regime for topological mechanical metamaterials and expands their application potential in areas including adaptive and smart materials and mechanical logic, wherein concepts from nonlinear dynamics may be used to create intricate, tailored spatial deformation and stress fields greatly transcending conventional elasticity.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Complex Dependence of Calcite Crack Kinetics on Salinity: The Role of DLVO and Hydration Forces

Abstract Subcritical crack growth (SCG) plays an important role in many geological processes such as delayed earth rupture and rock weathering. The complex dependency of SCG on the in‐crack fluid chemistry, however, is still poorly understood. In this study, we utilize the newly developed surface force‐based fracture theory (SFFT) to elucidate the relative contributions of surface forces and solute transport to the crack growth kinetics of calcite in NaCl solutions. Expanding on Barenblatt's cohesive crack model, SFFT introduces an effective stress intensity at the crack tip that encompasses all the relevant intermolecular forces across the crack in addition to the external far‐field stresses. The nonlinear system of equations portraying the crack opening profile, the solute distribution in a propagating crack, and the crack growth velocity are numerically solved via an implicit scheme. After carefully calibrating the model for calcite‐water systems, the SFFT is used to predict the SCG response of calcite at different NaCl concentrations, based on various hypotheses. These predictions are then compared to existing SCG data from the literature. We demonstrate that the experimentally observed variation of SCG rate with NaCl concentration cannot be explained solely by DLVO forces (electrostatic and Van der Waals interactions). This can be remediated by introducing an exponentially decaying hydration force with a nonlinear, nonmonotonic dependence on NaCl concentration. Furthermore, we demonstrate that accounting for both diffusive and advective transport of ions is important in explaining the absence of a stage‐II SCG response for calcite in electrolyte solutions. Plain Language Summary Subcritical crack growth (SCG) refers to the slow propagation of cracks in materials under a stress below the threshold for catastrophic failure. SCG is a key process in many geological events, for example, delayed earth ruptures and rock weathering. New initiatives such as underground CO 2 and H 2 storage in carbonate reservoirs further call for better understanding of SCG in carbonate minerals subjected to varying fluid chemistry. This study examines the SCG of calcite, a key mineral found in carbonate rocks, intergranular cement in sandstones, and filling material in mineral veins and faults, determining their deformation and strength. A mathematical model is developed to describe how the crack opens and propagates, how solutes (like salts) distribute within the crack, and how the crack surfaces interact with each other. We used the model to predict calcite SCG in water at different salt concentrations and compared it with experimental data. Our results revealed that the hydration force is the dominating factor in determining the complex, non‐linear dependency of SCG on salinity. We also found that both the movement of ions by diffusion and by bulk water flow are crucial for explaining the SCG rates, especially when the cracks grow quickly. Key Points Surface Force‐Based Fracture Theory predicts the complex subcritical crack growth patterns of calcite crystals immersed in NaCl solutions Results highlight the dominant role of hydration forces in altering the fracture behavior of calcite compared to VdW and electric double‐layer forces Advective solute transport explains the absence of stages‐II and ‐III subcritical crack growth responses in solid‐liquid systems

DLVO↗
DOE OSTI2025

Quantum Nonlinear Acoustic Hall Effect and Inverse Acoustic Faraday Effect in Dirac Insulators

Here, we propose to realize the quantum nonlinear Hall effect and the inverse Faraday effect through the acoustic wave in a time-reversal invariant but inversion broken Dirac insulator. We focus on the acoustic frequency much lower than the Dirac gap such that the interband transition is suppressed and these effects arise solely from the intrinsic valley-contrasting band topology. The corresponding acoustoelectric conductivity and magnetoacoustic susceptibility are both proportional to the quantized valley Chern number and independent of the quasiparticle lifetime. The linear and nonlinear components of the longitudinal and transverse topological currents can be tuned by adjusting the polarization and propagation directions of the surface acoustic wave. The static magnetization generated by a circularly polarized acoustic wave scales linearly with the acoustic frequency as well as the strain-induced charge density. Our results unveil a quantized nonlinear topological acoustoelectric response of gapped Dirac materials, like hexagonal boron nitride and transition-metal dichalcogenide, paving the way toward room-temperature acoustoelectric devices due to their large band gaps.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Data-Driven Approach for Uncertainty Propagation and Reachability Analysis in Dynamical Systems

In this paper, we propose a data-driven approach for uncertainty propagation and reachability analysis in a dynamical system. The proposed approach relies on the linear lifting of a nonlinear system using linear Perron-Frobenius (P-F) and Koopman operators. Here, the uncertainty can be characterized in terms of the moments of a probability density function. We demonstrate how the P-F and Koopman operators are used for propagating the moments. Time-series data is used for the finite-dimensional approximation of the linear operators, thereby enabling data-driven approach for moment propagation. Simulation results are presented to demonstrate the effectiveness of the proposed method.

42 ENGINEERING↗
DOE OSTI2021

Single-shot measurement of the nonlinear refractive index of air at 92 µm with a picosecond terawatt CO 2 laser

We developed a simple, accurate single-shot method to determine the nonlinear refractive index of air by measuring the evolution of the spatial shape of a laser beam propagating through the atmosphere. A distinctive feature of this new method, which relies on a modified Fresnel propagation model for data analysis, is the use of a hard aperture for producing a welldefined, high-quality beam from a comparatively nonuniform quasi-flat-top beam which is typical for highpeak-power lasers. The nonlinear refractive index of air for a very short (2 ps) long-wave infrared (LWIR) laser pulse was measured for the first time, yielding n 2 = 3.0 × 10 -23 m 2 /W at 9.2 µm. This result is 40% lower than a corresponding measurement with longer (200 ps) LWIR pulses at a similar wavelength.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Detection of Forced Change Within Combined Climate Fields Using Explainable Neural Networks

Abstract Assessing forced climate change requires the extraction of the forced signal from the background of climate noise. Traditionally, tools for extracting forced climate change signals have focused on one atmospheric variable at a time, however, using multiple variables can reduce noise and allow for easier detection of the forced response. Following previous work, we train artificial neural networks to predict the year of single‐ and multi‐variable maps from forced climate model simulations. To perform this task, the neural networks learn patterns that allow them to discriminate between maps from different years—that is, the neural networks learn the patterns of the forced signal amidst the shroud of internal variability and climate model disagreement. When presented with combined input fields (multiple seasons, variables, or both), the neural networks are able to detect the signal of forced change earlier than when given single fields alone by utilizing complex, nonlinear relationships between multiple variables and seasons. We use layer‐wise relevance propagation, a neural network explainability tool, to identify the multivariate patterns learned by the neural networks that serve as reliable indicators of the forced response. These “indicator patterns” vary in time and between climate models, providing a template for investigating inter‐model differences in the time evolution of the forced response. This work demonstrates how neural networks and their explainability tools can be harnessed to identify patterns of the forced signal within combined fields.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Strong-field QED limitations on TeV-class plasma wakefield accelerators

We demonstrate that quantum and classical radiation effects can become non-negligible for TeV-class beams propagating through plasma channels typical of staged plasma accelerators. Although the quantum nonlinearity parameter χ e remains small under currently envisioned experimental conditions, the cumulative influence of radiation over long acceleration distances can lead to significant modifications to the beam’s energy spread, emittance, and polarization. Our analytic models, validated by particle-in-cell simulations, highlight that for standard Gaussian beams, the orbit-induced energy spread dominates over quantum stochastic effects but can be mitigated by tailoring the beam profile, for example, through ring-shaped transverse distributions. In regimes where the radiation reaction approaches the accelerating force, the emittance may be cooled, forming distinctive ring-shaped phase-space structures. Finally, we analyze the influence of the radiation effect on spin transport inside the wakefield. These findings underscore the importance of considering both classical and quantum radiation dynamics in the design and optimization of future high-intensity plasma accelerators.

accelerator↗
DOE OSTI2022

Gaussian processes meet NeuralODEs: a Bayesian framework for learning the dynamics of partially observed systems from scarce and noisy data

We present a machine learning framework (GP-NODE) for Bayesian model discovery from partial, noisy and irregular observations of nonlinear dynamical systems. The proposed method takes advantage of differentiable programming to propagate gradient information through ordinary differential equation solvers and perform Bayesian inference with respect to unknown model parameters using Hamiltonian Monte Carlo sampling and Gaussian Process priors over the observed system states. This allows us to exploit temporal correlations in the observed data, and efficiently infer posterior distributions over plausible models with quantified uncertainty. The use of the Finnish Horseshoe as a sparsity-promoting prior for free model parameters also enables the discovery of parsimonious representations for the latent dynamics. A series of numerical studies is presented to demonstrate the effectiveness of the proposed GP-NODE method including predator–prey systems, systems biology and a 50-dimensional human motion dynamical system. This article is part of the theme issue ‘Data-driven prediction in dynamical systems’.

Science & Technology - Other Topics↗
DOE OSTI2024

Apparent fracture toughness estimation of additively manufactured alumina from as-printed chevron-notched test specimens

Ceramic vat photopolymerization (VPP) is a digital light processing method used to make additively manufactured green ceramic components that are then sintered. Despite the ever-advancing maturation of this method's green-state process, some of the material properties of sintered VPP-processed ceramics are still not well defined or understood. One example is Mode I fracture toughness, K Ic . In this study, attention was devoted to K Ic measurement using net-shape chevron-notched bend bars that were VPP-processed and the examination of whether valid K Ic measurement could occur by testing them. Further, stable crack propagation was observed in >60% of tested samples, indicated by a smooth nonlinear transition through the measured maximum force prior to final fracture. However, their results cannot yet be responsibly referred to as K Ic because of continuing violations of other prerequisites needed for valid K Ic testing.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Optical bistability in PECVD silicon-rich nitride

We present a study of optical bi-stability in a 3.02 refractive index at 1550nm plasma enhanced chemical vapor deposition (PECVD) silicon-rich nitride (SRN) film, as it pertains to bi-stable switching, memory applications, and thermal sensing applications. In this work we utilize an SRN ring resonator device, which we first characterize at low-power and then compare thermo-optic coefficients, (2.12 ± 0.125) × 10 −4 /°C, obtained from thermal-heating induced resonance shifts to optically induced resonance shifts as well as estimated propagation loss and absorption. We then measure the time response of this nonlinearity demonstrating the relaxation time to be 18.7 us, indicating the mechanism to be thermal in nature. Finally, we demonstrate bi-stable optical switching.

42 ENGINEERING↗
DOE OSTI2020

Bayesian differential programming for robust systems identification under uncertainty

This paper presents a machine learning framework for Bayesian systems identification from noisy, sparse and irregular observations of nonlinear dynamical systems. The proposed method takes advantage of recent developments in differentiable programming to propagate gradient information through ordinary differential equation solvers and perform Bayesian inference with respect to unknown model parameters using Hamiltonian Monte Carlo sampling. This allows an efficient inference of the posterior distributions over plausible models with quantified uncertainty, while the use of sparsity-promoting priors enables the discovery of interpretable and parsimonious representations for the underlying latent dynamics. A series of numerical studies is presented to demonstrate the effectiveness of the proposed methods, including nonlinear oscillators, predator–prey systems and examples from systems biology. Taken together, our findings put forth a flexible and robust workflow for data-driven model discovery under uncertainty. All codes and data accompanying this article are available at https://bit.ly/34FOJMj .

Science & Technology - Other Topics↗
DOE OSTI2023

Inversion-based correction of Double-Torsion (DT) subcritical crack growth tests for crack profile geometry

Because of its simplicity and the ability to produce a stable, slow-propagating crack, the Double-Torsion (DT) method has been used widely for investigating the critical and subcritical propagation of a slow-propagating tensile (mode-I) crack. However, to determine the complex relationship between the crack velocity $\mathcal{v_c}$ vs. the strain energy release rate $\mathscr{G}$ (or the stress intensity factor K) from laboratory measurements, several corrections must be made to account for the impact of sample and crack geometry. Particularly, DT test typically produces a crack with a curved edge profile instead of a straight line, causing the local $\mathcal{v_c}$ and $\mathscr{G}$ vary along the crack front. The experimentally measured $\mathcal{v_c}$ and $\mathscr{G}$ data merely reflect collective, averaged behavior of the crack. This makes inversion for the intrinsic, “true” crack growth kinetics necessary, based upon the knowledge of the crack geometry. Simple and effective correction methods have been proposed and validated for the slow, chemical-reaction-controlled part (Region I) of the $\mathcal{v_c}$-$\mathscr{G}$ curve. However, reliable methods for the highly nonlinear, transport-dominated part (Region II) and its sudden transition to the dynamic propagation part (Region III) are still lacking. Here we propose a method for determining the intrinsic $\mathcal{v_c}$-$\mathscr{G}$ relationship across all three Regions based upon DT test data, using a simple model function and its numerical inversion. The performance of this approach is examined and demonstrated using both synthetic and laboratory data for subcritical crack growth in soda lime glass.

42 ENGINEERING↗