Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nonlinear”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Precision Tests of the Nonlinear Mode Coupling of Anisotropic Flow via High-Energy Collisions of Isobars

Valuable information on dynamics of expanding fluids can be inferred from the response of such systems to perturbations in their initial geometry. We apply this technique in high-energy 96 Ru+ 96 Ru and 96 Zr+ 96 Zr collisions to scrutinize the expansion dynamics of the quark-gluon plasma, where the initial geometry perturbations are sourced by the differences in deformations and radial profiles between 96 Ru and 96 Zr, and the collective response is captured by the change in anisotropic flow V n between the two collision systems. Using a transport model, we analyze how the nonlinear coupling between lower-order flow harmonics V 2 and V 3 to the higher-order flow harmonics V 4 and V 5 , expected to scale as V 4 NL = χ 4 V 2 2 and V 5NL = χ 5 V 2 V 3 , gets modified as one moves from 96 Ru+ 96 Ru to 96 Zr+ 96 Zr systems. We find that these scaling relations are valid to high precision: variations of order 20% in V 4NL and V 5NL due to differences in quadrupole deformation, octupole deformation, and nuclear skin modify χ 4 and χ 5 by about 1–2%. Percent-level deviations are however larger than the expected experimental uncertainties and could be measured. Therefore, collisions of isobars with different nuclear structures are a unique tool to isolate subtle nonlinear effects in the expansion of the quark-gluon plasma that would be otherwise impossible to access in a single collision system.

Physics↗
DOE OSTI2022

The optimal beam-loading in two-bunch nonlinear plasma wakefield accelerators

Abstract Due to the highly nonlinear nature of the beam-loading, it is currently not possible to analytically determine the beam parameters needed in a two-bunch plasma wakefield accelerator for maintaining a low energy spread. Therefore in this paper, by using the Broyden–Fletcher–Goldfarb–Shanno algorithm for the parameter scanning with the code QuickPIC and the polynomial regression together with k -fold cross-validation method, we obtain two fitting formulas for calculating the parameters of tri-Gaussian electron beams when minimizing the energy spread based on the beam-loading effect in a nonlinear plasma wakefield accelerator. One formula allows the optimization of the normalized charge per unit length of a trailing beam to achieve the minimal energy spread, i.e. the optimal beam-loading. The other one directly gives the transformer ratio when the trailing beam achieves the optimal beam-loading. A simple scaling law for charges of drive beams and trailing beams is obtained from the fitting formula, which indicates that the optimal beam-loading is always achieved for a given charge ratio of the two beams when the length and separation of two beams and the plasma density are fixed. The formulas can also help obtain the optimal plasma densities for the maximum accelerated charge and the maximum acceleration efficiency under the optimal beam-loading respectively. These two fitting formulas will significantly enhance the efficiency for designing and optimizing a two-bunch plasma wakefield acceleration stage.

Physics↗
DOE OSTI2023

A solvable model of a nonlinear extension of quantum mechanics

We introduce a particular nonlinear generalization of quantum mechanics which has the property that it is exactly solvable in terms of the eigenvalues and eigenfunctions of the Hamiltonian of the usual linear quantum mechanics problem. Here, we hope that this simple example will elucidate some of the issues of interpreting nonlinear generalization of quantum mechanics that have been put forth to resolve questions about quantum measurement theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Joint velocity and density reconstruction of the Universe with nonlinear differentiable forward modeling

Reconstructing the initial conditions of the Universe from late-time observations has the potential to optimally extract cosmological information. Due to the high dimensionality of the parameter space, a differentiable forward model is needed for convergence, and recent advances have made it possible to perform reconstruction with nonlinear models based on galaxy (or halo) positions. In addition to positions, future surveys will provide measurements of galaxies' peculiar velocities through the kinematic Sunyaev-Zel'dovich effect (kSZ), type Ia supernovae, the fundamental plane relation, and the Tully-Fisher relation. Here we develop the formalism for including halo velocities, in addition to halo positions, to enhance the reconstruction of the initial conditions. We show that using velocity information can significantly improve the reconstruction accuracy compared to using only the halo density field. Here we study this improvement as a function of shot noise, velocity measurement noise, and angle to the line of sight. We also show how halo velocity data can be used to improve the reconstruction of the final nonlinear matter overdensity and velocity fields. We have built our pipeline into the differentiable Particle-Mesh FlowPM package, paving the way to perform field-level cosmological inference with joint velocity and density reconstruction. This is especially useful given the increased ability to measure peculiar velocities in the near future.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Nonlinear two-fluid modeling of plasma response to RMPs for the ELM control in the ITER baseline

Numerical modeling, combining the toroidal ideal MHD code GPEC and the nonlinear two-fluid MHD code TM1, was used for comprehensive studies of the plasma response to resonant magnetic perturbations (RMPs) with toroidal mode number n = 1–5 for controlling edge-localized modes (ELMs) in ITER for the standard operation scenario (15 MA Q = 10). Several issues related to RMP ELM control are investigated, including the optimization of the RMP coils configuration, the evaluation of the magnitude of density pump-out and the q95 windows of ELM suppression. GPEC calculates the magnetic response, which consistently includes the very important edge kink/peeling response to static magnetic perturbations. Furthermore, GPEC two-dimensional scans of the relative coil current phasing among the three rows of internal coils, at fixed coil current amplitude, reveal the optimal phasing for the RMP coil configuration with n = 1–5, respectively. The poloidal half wavelength of resonant mode at the edge of plasma calculated by GPEC indicates that the midplane row coils have the best resonant coupling with the plasma for n = 2, while the upper and lower row coils have the best resonant coupling with the plasma for n = 3. Based on the plasma kinetic equilibrium and the GPEC calculations of the magnetic response, TM1 was used to simulate the conditions for RMP field penetration in the ITER pedestal. TM1 shows magnetic island formation at the foot of ITER pedestal with RMP coil current threshold ranging from 4 kAt to 8 kAt with n = 2 to 4. These magnetic islands at the pedestal-foot lead to density pump-out, the magnitude of which scales as ${I}_{\text{RMP}}^{0.5}$ and ranges from 5% to 20% at the pedestal-top when scanning the coil current from 4 to 60 kAt. The density pump-out is found to be weaker for higher n RMP. The nonlinear TM1 simulations also show field penetration at the pedestal-top, where the threshold of RMP coil current depends on the q 95 . The alignment of the magnetic island and the location of the pedestal-top decreases the height and width of the pedestal to suppress ELMs. Simulations by two-dimensional scans of RMP coil current and q 95 reveal the accessible q 95 windows of ELM suppression for both n = 3 and 4 RMPs. The predicted q 95 windows of ELM suppression are very similar to the ones in currently operating tokamaks and the required RMP coil current for ELM suppression is less than 40–50 kAt, which is well within the designed capability for ITER. In addition, the simulations indicate that wide q 95 windows of ELM suppression may be accessible in ITER by operating with dominant n = 4 (or n = 5) RMPs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Global linear and nonlinear gyrokinetic simulations of tearing modes

Abstract To better understand multi-scale interactions between global tearing modes and microturbulence in the Madison Symmetric Torus (MST) reversed-field pinch (RFP), the global gyrokinetic code Gene is modified to describe global tearing mode instability via a shifted Maxwellian (SM) distribution consistent with experimental equilibria. The implementation of the SM is tested and benchmarked by comparisons with different codes and models. Good agreement is obtained in code-code and code-theory comparisons. Linear stability of tearing modes of a non-reversed MST discharge is studied. A collisionality scan is performed to the lowest order unstable modes ( n = 5, n = 6) and shown to behave consistently with theoretical scaling. The nonlinear evolution is simulated, and saturation is found to arise from mode coupling and transfer of energy from the most unstable tearing mode to small-scale stable modes mediated by the m = 2 tearing mode. The work described herein lays the foundation for nonlinear simulation and analysis of the interaction of tearing modes and gyroradius-scale instabilities in RFP plasmas.

Physics↗
DOE OSTI2024

Divertor-safe nonlinear burn control based on a SOLPS parameterized core-edge model for ITER

Abstract For ITER operations, the range of desirable burning-plasma regimes with high fusion power output will be restricted by various operational constraints. These constraints include the saturation of ITER’s various heating and fueling actuators such as the neutral beam injectors, the ion and electron cyclotron heating systems, the gas puffing system, and the deuterium–tritium pellet injectors. In addition to these actuator constraints, the H-mode power threshold, divertor detachment, and the heat load on the divertor targets may apply limitations to ITER’s operational space. In this work, Plasma Operation Contour (POPCON) plots that map the aforementioned constraints to the temperature-density space are used to investigate which constraints are most limiting towards accessing regimes with high fusion power output. The presented POPCON plots are based on a control-oriented core-edge model that couples the nonlinear density and energy response models for the core-plasma region with SOLPS4.3 parameterizations for conditions in the edge-plasma regions (scrape-off-layer and divertor). Using this control-oriented core-edge model, a nonlinear burn controller, which aims to regulate the plasma temperature and density in the core-plasma region, is constructed in this work. This controller is augmented with an online optimization scheme that governs the control references such that the plasma can be guided towards regimes with high fusion powers while protecting the divertor targets from dangerously high heat loads. A closed-loop simulation study illustrates the capability of this burn control scheme.

Physics↗
DOE OSTI2024

Nonlinear modeling of ELM mitigation with RMP on HL-2A

Abstract Nonlinear modeling of mitigation of the edge localized mode (ELM) with resonant magnetic perturbation (RMP) is performed for the HL-2A tokamak, utilizing the three-dimensional (3D) magnetohydrodynamic code JOREK. Based on the 3D equilibrium established after application of the n = 1 ( n is the toroidal mode number) RMP at 4.9 kAt coil current with odd parity, ELM mitigation is successfully simulated consistent with the experimental result. Nonlinear simulations show strong mode coupling among toroidal Fourier harmonics, allowing redistribution of the magnetic energy such that the most unstable toroidal mode saturates at a lower level. This magnetic energy cascade offers an explanation of the RMP-induced ELM mitigation achieved in HL-2A. Detailed examination of the simulation results shows persistent resonant field screening even during the ELM mitigation phase. Finite plasma resistivity however does enable partial penetration of the resonant field thus modifying the edge magnetic topology and characteristics of the edge transport. Plasma radial profiles undergo pronounced changes around the pedestal region, when the magnetic energy of the most unstable toroidal mode reaches the maximum value. Systematic scans of the applied RMP coil current with the JOREK simulations find a threshold value of around 4.5 kAt required for achieving the ELM mitigation on HL-2A.

Physics↗
DOE OSTI2024

Prediction of performance and turbulence in ITER burning plasmas via nonlinear gyrokinetic profile prediction

Burning plasma performance, transport, and the effect of hydrogen isotope (H, D, D-T fuel mix) on confinement has been predicted for ITER baseline scenario (IBS) conditions using nonlinear gyrokinetic profile predictions. Accelerated by surrogate modeling (Rodriguez-Fernandez et al 2022 Nucl. Fusion 62 076036), high fidelity, nonlinear gyrokinetic simulations performed with the CGYRO code (Candy et al 2016 J. Comput. Phys. 324 73), were used to predict profiles of T i , T e , and n e while including the effects of alpha heating, auxiliary power (NBI + ECH), collisional energy exchange, and radiation losses inside of $r/a$ = 0.9. Predicted profiles and resulting energy confinement are found to produce fusion power and gain that are approximately consistent with mission goals ($P_\textrm{fusion} = 500$ MW at Q = 10) for the baseline scenario and exhibit energy confinement that is within 1σ of the H-mode energy confinement scaling. The power of the surrogate modeling technique is demonstrated through the prediction of alternative ITER scenarios with reduced computational cost. These scenarios include conditions with maximized fusion gain and an investigation of potential resonant magnetic perturbation (RMP) effects on performance with a minimal number of gyrokinetic profile iterations required (3–6). These predictions highlight the stiff ITG nature of the core turbulence predicted in the ITER baseline and demonstrate that $Q \gt$ 17 conditions may be accessible by reducing auxiliary input power while operating in IBS conditions. Prediction of full kinetic profiles allowed for the projection of hydrogen isotope effects around ITER baseline conditions. The gyrokinetic fuel ion species was varied from H, D, and 50/50 D-T and kinetic profiles were predicted. Results indicate that a weak or negligible isotope effect will be observed to arise from core turbulence in IBS conditions. The resulting energy confinement, turbulence, and density peaking, and the implications for ITER operations will be discussed.

gyrokinetics↗
DOE OSTI2021

Behavior of solitary waves of coupled nonlinear Schrödinger equations subjected to complex external periodic potentials with odd-$\mathcal{PT}$ symmetry

In this work, we discuss the response of both moving and trapped solitary wave solutions of a two-component nonlinear Schrödinger system in 1 + 1 dimensions to an odd-$\mathcal{PT}$ external periodic complex potential. The dynamical behavior of perturbed solitary waves is explored by conducting numerical simulations of the nonlinear system and using a collective coordinate variational approximation. We present case examples corresponding to choices of parameter values and initial conditions involved therein. The results of the collective coordinate approximation are compared against numerical simulations where we observe qualitatively good agreement between the two. Unlike the case for a single-component solitary wave in a complex periodic $\mathcal{PT}$-symmetric potential, the collective coordinate equations do not have a small oscillation regime, and initially the height of the two components changes in opposite directions often causing instability. We find that the dynamic stability criteria we have used in the one-component case are a good indicator for the onset of dynamic instabilities in the present setup.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Photon Frequency Conversion in High-Q Superconducting Resonators: Axion Electrodynamics, QED, and Nonlinear Meissner Radiation

High-Q superconducting resonators have been proposed and developed as detectors of light-by-light scattering mediated by the hypothesized axion or virtual electron–positron pairs in quantum electrodynamics: the Euler–Heisenberg (EH) interaction. Photon frequency and mode conversion is central to the scheme for detecting such rare events. Superconducting resonators are nonlinear devices. The Meissner screening currents that confine the electromagnetic fields to the vacuum region of a superconducting RF cavity are nonlinear functions of the EM field at the vacuum–superconducting interface, and as a result can generate source currents and frequency conversion of microwave photons in the cavity. In this report we consider photon frequency and mode conversion in superconducting resonators with high quality factors from Meissner currents in single- and dual-cavity setups proposed for axion and QED searches based on light-by-light scattering. In a single cavity with two pump modes, photon frequency conversion by the Meissner screening current dominates photon generation by the EH interaction for cavities with Q≲1012⁠. The Meissner currents also generate background photons that limit the operation of the resonator for axion detection in three-mode, single-cavity setups. We also consider the leakage of photons from pump modes into the signal mode for both axion- and EH-mediated light-by-light scattering. Photon frequency conversion by the EH interaction can compete with Meissner and leakage radiation in ultra-high-Q cavities that are beyond the current state of the art. Meissner radiation and leakage backgrounds can be suppressed in dual-cavity setups with appropriate choices for pump and spectator modes, as well as the single-cavity setup proposed for heterodyne detection of galactic axion dark matter.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Nonlinear elasticity and short-range mechanical coupling govern the rate and symmetry of mouth opening in Hydra

Hydra has a tubular bilayered epithelial body column with a dome-shaped head on one end and a foot on the other. Hydra lacks a permanent mouth: its head epithelium is sealed. Upon neuronal activation, a mouth opens at the apex of the head which can exceed the body column diameter in seconds, allowing Hydra to ingest prey larger than itself. While the kinematics of mouth opening are well characterized, the underlying mechanism is unknown. We show that Hydra mouth opening is generated by independent local contractions that require tissue-level coordination. We model the head epithelium as an active viscoelastic nonlinear spring network. The model reproduces the size, timescale and symmetry of mouth opening. It shows that radial contractions, travelling inwards from the outer boundary of the head, pull the mouth open. Nonlinear elasticity makes mouth opening larger and faster, contrary to expectations. Here, the model correctly predicts changes in mouth shape in response to external forces. By generating innervated : nerve-free chimera in experiments and simulations, we show that nearest-neighbor mechanical signaling suffices to coordinate mouth opening. Hydra mouth opening shows that in the absence of long-range chemical or neuronal signals, short-range mechanical coupling is sufficient to produce long-range order in tissue deformations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Nuclear deformation as a source of the nonlinearity of the King plot in the Yb + ion

We perform atomic relativistic many-body calculations of the field isotope shifts and calculations of corresponding nuclear parameters for all stable even-even isotopes of the Yb + ion. We demonstrate that if we take nuclear parameters of the Yb isotopes from a range of the state of the art nuclear models, which all predict strong quadrupole nuclear deformation, and then calculate nonlinearity of the King plot caused by the difference in the deformation in different isotopes, the result is consistent with the nonlinearity observed in the experiment. In conclusion, the changes of nuclear rms radius between isotopes extracted from experiment are consistent with those obtained in the nuclear calculations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Experimental observation of polarization-resolved nonlinear Thomson scattering of elliptically polarized light

We report experimental results from a study of nonlinear Thomson scattering of elliptically polarized light. Polarization-resolved radiation patterns of the scattered light are measured as a function of the elliptical polarization state of the incident laser light. The relativistic electron trajectory in intense elliptically polarized fields leads to the formation of unique radiated polarization states, which are observed by our measurements and predicted by a theoretical model. The polarization of Thomson scattered light depends strongly on the intensity of the incident light due to nonlinearity. The results are relevant to high-field electrodynamics and to research and development of light sources with novel capabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Simple model for the nonlinear radiation field of a free electron laser

It is shown that Jacobi elliptic function solutions of the nonlinear Duffing equation model the radiation field in a high-gain free electron laser through early saturation. After initial start-up, the field can be expressed equivalently with a hyperbolic secant. The model is derived for arbitrary detuning from resonance, which enables study of the spectral properties in the early nonlinear regime.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Nonlinear time-domain spectroscopy near a band inversion

We develop a theory for the nonlinear time-domain response of a semimetal driven by an ultrafast optical pulse. At quadratic order in the driving field we find that the response near a band inversion transition contains a coherent oscillating component proportional to field intensity with a frequency that can be tuned over a wide range (from terahertz to the near IR) selected by the chemical potential. Finally, we illustrate the effect by calculating the induced current as a function of time for a two-band system with time-reversal symmetry but broken parity symmetry to study semimetals where large Berry curvature near a band inversion transition promotes this nonlinear response.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Exact results for nonlinear Drude weights in the spin-1/2 XXZ chain

Nonlinear Drude weight (NLDW) is a generalization of the linear Drude weight, which characterizes the nonlinear transport in quantum many-body systems. We investigate these weights for the spin-1/2 XXZ chain in the critical regime. The effects of the Dzyaloshinskii-Moriya interaction and an external magnetic field are also studied. Solving the Bethe equations numerically, we obtain these weights for very large system sizes and identify parameter regimes where the weights diverge in the thermodynamic limit. These divergences appear in all the orders studied in this Letter and can be regarded as a generic feature of the NLDWs. We study the origin of these divergences and reveal that they result from nonanalytic finite-size corrections to the ground-state energy. As a result, we compute closed-form expressions for several weights in the thermodynamic limit and find excellent agreement with the numerical results.

1-dimensional systems↗
DOE OSTI2022

High-pressure control of optical nonlinearity in the polar Weyl semimetal TaAs

The transition metal monopnictide family of Weyl semimetals recently has been shown to exhibit anomalously strong second-order optical nonlinearity, which is theoretically attributed to a highly asymmetric polarization distribution induced by their polar structure. We experimentally test this hypothesis by measuring optical second harmonic generation (SHG) from TaAs across a pressure-tuned polar-to-nonpolar structural phase transition. Despite the high-pressure structure remaining noncentrosymmetric, the SHG yield is reduced by more than 60% by 20 GPa as compared to the ambient pressure value. By examining the pressure dependence of distinct groups of SHG susceptibility tensor elements, we find that the yield is primarily controlled by a single element that governs the response along the polar axis. Furthermore, our results confirm a connection between the polar axis and the giant optical nonlinearity of Weyl semimetals and demonstrate pressure as a means to tune this effect in situ.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗