Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nonlinear waves”

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 703 records · Page 39

Multifidelity validation of digital surrogates using variable-density turbulent mixing models

In this study, ensembles of experimental data are presented and utilized to compare and validate two models used in the simulation of variable-density (Atwood = 0.22), compressible turbulent mixing. Though models of this kind (Reynolds averaged NavierStokes and large-eddy simulations) have been validated extensively with more canonical flows in previous studies, here the present approach offers novelty in the complexity of the geometry, the ensemble-based validation, and the uniformity of the computational framework on which the models are tested. Moreover, all experimental and computational tasks were completed by the authors which has led to a tightly coupled experimental configuration with its “digital twin.” The experimental divergent-shock-tube facility and its data acquisition methods are described and replicated in simulation space. A 2D Euler model which neglects the turbulent mixing at the interface is optimized to experimental data using a Gaussian process. This model then serves as the basis for both the 2D RANS and 3D LES studies that make comparisons to the mixing-layer data from the experiment. A relatively simple RANS model is shown to produce good agreement with experimental data only at late flow development times. The LES ensembles generally show good agreement with experimental data but display sensitivity to the characterization of initial conditions. Resolution-dependent behavior is also observed for certain higher-order statistics of interest. Overall, the LES model successfully captures the effects of divergent geometry, compressibility, and combined nonlinear instabilities inherent to the problem. The successful prediction of mixing width and its growth rate highlight the existence of three distinct regimes in the development of the instability, each with similarities to previously studied instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The saturation mechanism of thermal instability

The literature on thermal instability (TI) reveals that even for a simple homogeneous plasma, the nonlinear outcome ranges from a gentle reconfiguration of the initial state to an explosive one, depending on whether the condensations that form evolve in an isobaric or nonisobaric manner. After summarizing the recent developments on the linear and nonlinear theory of TI, here we derive several general identities from the evolution equation for entropy that reveal the mechanism by which TI saturates; whenever the boundary of the instability region (the Balbus contour) is crossed, a dynamical change is triggered that causes the comoving time derivative of the pressure to change the sign. This event implies that the gas pressure force reverses direction, slowing the continued growth of condensation. For isobaric evolution, this “pressure reversal” occurs nearly simultaneously for every fluid element in condensation and a steady state is quickly reached. For nonisobaric evolution, the condensation is no longer in mechanical equilibrium and the contracting gas rebounds with greater force during the expansion phase that accompanies the gas reaching the equilibrium curve. The cloud then pulsates because the return to mechanical equilibrium becomes wave mediated. We show that both the contraction rebound event and subsequent pulsation behavior follow analytically from an analysis of the new identities. Our analysis also leads to the identification of an isochoric TI zone and makes it clear that unless this zone intersects the equilibrium curve, isochoric modes can only become unstable if the plasma is in a state of thermal non-equilibrium.

79 ASTRONOMY AND ASTROPHYSICS↗

Development of five-moment two-fluid modeling for Z-pinch physics

The Z-pinch m = 0 instability as well as its stabilization by radially sheared axial flow is studied using the nonlinear ideal five-moment twofluid (5M2F) model with an extension of that model to include Braginskii heat and momentum transport. Using the ideal 5M2F model, linear growth rate results are compared with prior work using MHD and Hall MHD. At small normalized wavenumber, 1 < k a < 4 , where a is the effective pinch radius, 5M2F results agree with Hall MHD within 20% in scenarios without radially sheared axial flow. With the sheared flow and focusing on ka = 10/3, agreement with Hall MHD is excellent. In the limit of small ion inertial length, results also match with MHD. A comparison with PIC modeling of shear-free m = 0 stability focuses on a plasma scenario based on recent experimental results. In a scan of mode wavenumber, ideal 5M2F results are qualitatively similar to PIC: the growth rate rises to a peak at a moderate wavenumber and declines at a large wavenumber in contrast to MHD results, which show the saturation of the growth rate with the increasing wavenumber rather than a decline. The peak normalized 5M2F growth rate is &#x3B3; &#x3C4; A = 1.5 , where sA is the Alfven transit time across the pinch. The peak occurs at normalized wavenumber ka = 10. For comparison, PIC results have a peak growth of &#x3B3; &#x3C4; A = 0.77 at ka = 5. Including Braginskiibased closure of the 5M2F model does not qualitatively change the ideal results in this particular case. Nonlinear saturation is studied using the 5M2F model with the dissipative Braginskii-based closure in cases with pinch-edge sheared-flow speed equal to half the Alfven speed. Nonlinear mixing due to the sheared flow yields a quasi-steady state after modest losses of pinch ion inventory and pinch thermal energy, approximately 30% and 10%, respectively. 5M2F modeling captures the essential physics of m = 0 instability and offers a computationally tractable route to high-fidelity modeling of 3D Z-pinch behavior, including m = 1 instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Angle-resolved nonresonant two-photon single ionization of argon using 9.3-eV photons produced via high-order harmonic generation

Here we present an experimental study on the photoionization dynamics of nonresonant one-color two-photon single valence ionization of neutral argon atoms. Using 9.3-eV photons produced via high harmonic generation and a three-dimensional momentum imaging spectrometer, we detect the photoelectrons and ions produced from nonresonant two-photon ionization in coincidence. Photoionization from the $3p$ orbital produces a photoelectron scattering wave function with $p$- and $f$-partial-wave components, which interfere and result in a photoelectron angular distribution with peak amplitude perpendicular to the vacuum ultraviolet polarization. The comparison between the present results and the two previous sets of theoretical calculations [Phys. Rev. A 44, 324 (1991), and J. Phys. B: At. Mol. Phys. 16, 2737 (1983)] indicates that electron-electron correlation contributes appreciably to the two-photon ionization dynamics.

74 ATOMIC AND MOLECULAR PHYSICS↗

Resistive wall tearing mode disruptions in DIII-D and ITER tokamaks

Disruptions are a serious problem in tokamaks, in which thermal and magnetic energy confinement is lost. This paper uses data from the DIII-D experiment, theory, and simulations to demonstrate that resistive wall tearing modes (RWTMs) produce the thermal quench (TQ) in a typical locked mode shot. Analysis of the linear RWTM dispersion relation shows the parameter dependence of the growth rate, particularly on the resistive wall time. Linear simulations of the locked mode equilibrium show that it is unstable with a resistive wall and stable with an ideally conducting wall. Nonlinear simulations demonstrate that the RWTM grows to sufficient amplitude to cause a complete thermal quench. The RWTM growth time is proportional to the thermal quench time. Here, the nonlinearly saturated RWTM magnetic perturbation amplitude agrees with experimental measurements. The onset condition is that the q = 2 rational surface is sufficiently close to the resistive wall. Collectively, this identifies the RWTM as the cause of the TQ. In ITER, RWTMs will produce long TQ times compared to present-day experiments. ITER disruptions may be significantly more benign than previously predicted.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Instantaneous difference frequency locking observed during toroidicity-induced Alfvén eigenmode coupling in the DIII-D tokamak

In magnetic confinement fusion, toroidicity-induced Alfvén eigenmodes (TAEs) are well-studied, weakly stable solutions of the linearized ideal magnetohydrodynamic equations. Driven unstable by suprathermal populations of energetic particles, TAE pose a key vulnerability to the confinement of high-energy alpha particles generated by fusion reactions. Hence, it is paramount to understand TAE dynamics if a working reactor is to be realized. In this work, we detect and characterize signatures of nonstationary nonlinear coupling between TAE using a novel, time-resolved bispectral analysis; results are supported by analytic signal of band-passed data. Crucially, a stationary phase relationship between two TAE and a nascent low frequency fluctuation is observed precisely when the triple product of magnetic fluctuation amplitudes is enhanced. Local mode number and frequency spectrum, gleaned from beam-emission spectroscopy, corroborates simultaneous satisfaction of nonlinear matching conditions, and provides a tool to identify theorized pathways of energy transfer, e.g. TAE parametric instability.

bispectral analysis↗

Theory of terahertz pulse transmission through ferroelectric nanomembranes

An analytical model is developed to predict the temporal evolution of the lattice polarization in ferroelectric nanomembranes upon the excitation by a terahertz (THz) electromagnetic pulse of an arbitrary waveform and the concurrent transmission of the THz pulse in both linear and nonlinear regimes. It involves the use of the perturbation method to solve the equation of motion for the lattice polarization in both unclamped and strained ferroelectric nanomembranes within the framework of Landau-Ginzburg-Devonshire theory. The model is applicable to perovskite oxides such as BaTiO 3 and SrTiO 3 , wurtzite Al 1−𝑥 ⁢Sc 𝑥 ⁢N, and trigonal LiNbO 3 . Our analytical model provides a theoretical basis for determining the thermodynamic and kinetic parameters of ferroelectric materials through a THz transmission experiment. The calculation results also suggest an approach to reversing the chirality of a circularly polarized THz pulse by harnessing the resonant polarization-photon coupling in ferroelectrics. This capability of chirality reversal, along with the high tunability from a strain applied along any arbitrarily oriented in-plane axis, provides new opportunities for THz wave modulation without relying on complex metasurface designs.

36 MATERIALS SCIENCE↗

Collaborative Research: Particle Simulation of Collisionless Magnetic Reconnection Under Finite Guide Field

Magnetic reconnection is believed to be one of the most fundamental physical processes in laboratory and space plasmas. Fundamental physics of collisionless reconnection under general conditions of finite guide magnetic field, realistic electron-to-ion mass ratio, and in 3D geometries is still not adequately understood. The major goal of this project is to study the physics of collisionless magnetic reconnection in a current sheet under a finite guide field and with a realistic ion-to-electron mass ratio mi/me. Under the support of this DoE grant, we have made significant progress in (1) the investigation of linear and nonlinear current sheet instabilities under a wide range of guide field using the 3-D gyro kinetic electron and fully kinetic ion (GeFi) particle simulation model and (2) the investigation of fully nonlinear reconnection process with a strong guide field with the GeFi simulations. Furthermore, (3) we have extended the study to the realistic geometry of the magnetopause using 3D global hybrid simulations. This study, from local to global scales and for 3-D physics at plasma boundaries, has provides advanced theoretical understanding of the fundamental physics of magnetic reconnection in laboratory and space plasma regimes on multi-scales.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Model predictive current profile control in tokamaks by exploiting spatially moving electron cyclotron current drives

The plasma control systems in next-generation tokamaks like ITER will balance competing control objectives to achieve the desired level of performance in advanced scenarios while preventing magnetohydrodynamic instabilities and disruptions. During normal tokamak operation, the points of incidence of the electromagnetic waves generated by the electron cyclotron heating and current drives (EC H&CDs) are usually fixed in space. However, the points of incidence can be modified in real-time by changing the angles of the mirrors that reflect the EC H&CD waves. Altering the points of incidence, in turn, varies the ability of the plasma control system to regulate a plasma property. For instance, changing the EC H&CD wave incidence location may place the power demands necessary to achieve a particular plasma target within saturation limits. Therefore, using the EC H&CD deposition location, which is related to the EC H&CD mirror angle, as a supplementary controllable variable may facilitate access to a given target scenario. However, active scenario-control algorithms have not been designed so far to fully exploit this capability in real time. In this work, a model predictive controller that can handle actuation locations as control inputs is developed. In particular, the controller is designed to regulate both the auxiliary powers and the EC H&CD deposition locations in a pre-defined optimal sense to achieve the control objective of attaining and sustaining a target current profile. Further, the proposed controller is tested for a DIII-D tokamak scenario in nonlinear simulations using the Control Oriented Transport SIMulator (COTSIM).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Resistive hose modes in tokamak runaway electron beams II

Resistive hose instabilities of runaway electron (RE) beams immersed in resistive background plasmas are examined with analytic and numerical calculations. The RE beam-plasma equilibria considered are characteristic of the situation observed post-thermal quench in a tokamak disruption. An analytic linear dispersion relation is presented for the case of a uniform RE current density profile with a sharp boundary in cylindrical geometry. Initial value linear calculations for a more general profile in toroidal geometry find that reducing aspect ratio increases the resistive hose mode growth rates with fixed safety factor profile. Nonlinear calculations in cylindrical and toroidal geometry find that the resistive hose instability-driven fluctuations relax the gradient of the current density profile. In toroidal geometry, changes to the magnetic topology are observed as a result of the resistive hose activity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Viscous Behavior of Clay‐Rich Rocks and Its Role in Focused Fluid Flow

Abstract Focused fluid flow is common in sedimentary basins worldwide, where flow structures often penetrate through sandy reservoir rocks, and clay‐rich caprocks. To better understand the mechanisms forming such structures, the impacts of the viscoelastic deformation and strongly nonlinear porosity‐dependent permeability of clay‐rich materials are assessed from an experimental and numerical modeling perspective. The experimental methods to measure the poroviscoelastic and transport properties of intact and remolded shale have been developed, and the experimental data is used to constrain the numerical simulations. It is demonstrated that viscoelastic deformation combined with nonlinear porosity‐dependent permeability triggers the development of localized flow channels, often imaged as seismic chimneys. The permeability inside a channel increases by several orders of magnitude compared to the background values. In addition, the propagation time scale and the channel size strongly depend on the material properties of the fluid and the rock. The time‐dependent behavior of the clay‐rich rock may play a key role in the long‐term integrity of the subsurface formations.

58 GEOSCIENCES↗

Effects of multiply charged ions on microturbulence-driven electron transport in partially magnetized plasmas

Herein, nonlinear interaction between kinetic instabilities in partially magnetized plasmas in the presence of multiply-charged ion streams is investigated using kinetic simulations. It was observed by Hara and Tsikata [Phys. Rev. E 102, 023202 (2020)] that the axial ion-ion two-stream instability (IITSI) due to singly and doubly charged ion streams, coupled with the azimuthal electron cyclotron drift instability (ECDI), enhances cross-field electron transport. In the present study, it is observed that the addition of triply charged ions (as a third ion species) contributes to damping of the excited modes, leading to a reduction in the cross-field electron transport. The net instability-driven electron transport is shown to be a function not only of the azimuthal modes, such as the ECDI, but of the multiple ion species that dictate the development of additional plasma waves. It is found that trapping of the higher ion charge states within the plasma waves results in a broadening of the ion velocity distribution functions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Perspectives and progress on wurtzite ferroelectrics: Synthesis, characterization, theory, and device applications

Wurtzite ferroelectrics are an emerging material class that expands the functionality and application space of wide bandgap semiconductors. Promising physical properties of binary wurtzite semiconductors include a large, reorientable spontaneous polarization, direct band gaps that span from the infrared to ultraviolet, large thermal conductivities and acoustic wave velocities, high mobility electron and hole channels, and low optical losses. The ability to reverse the polarization in ternary wurtzite semiconductors at room temperature enables memory and analog type functionality and quasi-phase matching in optical devices and boosts the ecosystem of wurtzite semiconductors, provided the appropriate combination of properties can be achieved for any given application. In this article, advances in the design, synthesis, and characterization of wurtzite ferroelectric materials and devices are discussed. Highlights include: the direct and quantitative observation of polarization reversal of ~135 μC/cm 2 charge in Al 1-x B x N via electron microscopy, Al 1-x B x N ferroelectric domain patterns poled down to 400 nm in width via scanning probe microscopy, and full polarization retention after over 1000 h of 200 °C baking and a 2× enhancement relative to ZnO in the nonlinear optical response of Zn 1-x Mg x O. In conclusion, the main tradeoffs, challenges, and opportunities in thin film deposition, heterostructure design and characterization, and device fabrication are overviewed.

30 DIRECT ENERGY CONVERSION↗

Theory for anomalous terahertz emission in striped cuprate superconductors

Recent experiments in the doped cuprates La 2−𝑥 ⁢Ba 𝑥⁢ CuO 4 have revealed the emission of anomalous terahertz radiation after impulsive optical excitation. Here, we theoretically investigate the nonlinear electrodynamics of such striped superconductors and explore the origin of the observed radiation. We argue that photoexcitation is converted into a photocurrent by a second-order optical nonlinearity, which is activated by the breaking of inversion symmetry in certain stripe configurations. We point out the importance of including umklapp photocurrents modulated at the stripe periodicity itself, which impulsively drive surface Josephson plasmons and lead to a resonant structure of outgoing radiation, consistent with the experiments. We speculate on the utility of the proposed mechanism in the context of generating tunable terahertz radiation.

36 MATERIALS SCIENCE↗

Ion versus Electron Heating in Compressively Driven Astrophysical Gyrokinetic Turbulence

The partition of irreversible heating between ions and electrons in compressively driven (but subsonic) collisionless turbulence is investigated by means of nonlinear hybrid gyrokinetic simulations. We derive a prescription for the ion-to-electron heating ratio Q i /Q e as a function of the compressive-to-Alfvénic driving power ratio P compr /P AW , of the ratio of ion thermal pressure to magnetic pressure β i , and of the ratio of ion-to-electron background temperatures T i /T e . It is shown that Q i /Q e is an increasing function of P compr /P AW . When the compressive driving is sufficiently large, Q i /Q e approaches ≃P compr /P AW . This indicates that, in turbulence with large compressive fluctuations, the partition of heating is decided at the injection scales, rather than at kinetic scales. Analysis of phase-space spectra shows that the energy transfer from inertial-range compressive fluctuations to sub-Larmor-scale kinetic Alfvén waves is absent for both low and high β i , meaning that the compressive driving is directly connected to the ion-entropy fluctuations, which are converted into ion thermal energy. This result suggests that preferential electron heating is a very special case requiring low βi and no, or weak, compressive driving. Our heating prescription has wide-ranging applications, including to the solar wind and to hot accretion disks such as M87 and Sgr A*.

79 ASTRONOMY AND ASTROPHYSICS↗

Floquet stability analyses of stratified oscillating boundary layers on adiabatic slopes

The presence of a no-slip, impermeable, adiabatic, sloped boundary in an otherwise quiescent, stably stratified, Boussinesq flow generates baroclinic vorticity within a diffusive boundary layer. Such conditions are typical of the oscillating boundary layers on adiabatic abyssal slopes, sloped lake bathymetry and sloped coastal bathymetry in the absence of high-wavenumber internal waves, mean flows, far-field turbulence on larger scales, and resonant tidal–bathymetric interaction. We investigate the linear stability of the oscillating flow within non-dimensional parameter space typical of the $M_2$ tide and hydraulically smooth, middle-latitude abyssal slopes through Floquet linear stability analysis. The flow dynamics depends on three non-dimensional variables: the Reynolds number for Stokes’ second problem, the Prandtl number, and a frequency ratio that accounts for the resonance conditions ($C$, criticality) of the buoyant restoring force and the tidal forcing. The Floquet analysis results suggest that oscillating laminar boundary layers on adiabatic abyssal slopes are increasingly unstable as Reynolds number, criticality parameter and/or spanwise disturbance wavenumber are increased. We also show that the two-dimensional Floquet linear instability necessarily generates three-dimensional baroclinic vorticity, which suggests that the evolution of the gravitational instabilities may be nonlinear as $t→∞$.

42 ENGINEERING↗

Cooling and instabilities in colliding flows

ABSTRACT Collisional self-interactions occurring in protostellar jets give rise to strong shocks, the structure of which can be affected by radiative cooling within the flow. To study such colliding flows, we use the AstroBEAR AMR code to conduct hydrodynamic simulations in both one and three dimensions with a power-law cooling function. The characteristic length and time-scales for cooling are temperature dependent and thus may vary as shocked gas cools. When the cooling length decreases sufficiently and rapidly, the system becomes unstable to the radiative shock instability, which produces oscillations in the position of the shock front; these oscillations can be seen in both the one- and three-dimensional cases. Our simulations show no evidence of the density clumping characteristic of a thermal instability, even when the cooling function meets the expected criteria. In the three-dimensional case, the nonlinear thin shell instability (NTSI) is found to dominate when the cooling length is sufficiently small. When the flows are subjected to the radiative shock instability, oscillations in the size of the cooling region allow NTSI to occur at larger cooling lengths, though larger cooling lengths delay the onset of NTSI by increasing the oscillation period.

79 ASTRONOMY AND ASTROPHYSICS↗

The influence of Al2O3 on the structural properties of MgSiO3 akimotoite

Abstract Akimotoite, a MgSiO3 polymorph present in the lower transition zone within ultramafic portions of subducting slabs and potentially also in the ambient mantle, will partition some amount of Al, raising the question of how this will affect its crystal structure and properties. In this study, a series of samples along the MgSiO3-Al2O3 (akimotoite-corundum) solid solution have been investigated by means of single-crystal X-ray diffraction to examine their crystal chemistry. Results show a strong nonlinear behavior of the a- and c-axes as a function of Al content, which arises from fundamentally different accommodation mechanisms in the akimotoite and corundum structures. Furthermore, two Al2O3-bearing akimotoite samples were investigated at high pressure to determine the different compression mechanisms associated with Al substitution. Al2O3-bearing akimotoite becomes more compressible at least up to 20 mol% Al2O3, due likely to an increase in compressibility as the Al cation is incorporated into the SiO6 octahedron. This observation is in strong contrast to the stiffer corundum end-member having a KT = 250 GPa, which is larger than that of the akimotoite end-member [KT = 205(1) GPa]. These findings have implications for mineral physics models of elastic properties, which have in the past assumed linear mixing behavior between the MgSiO3 akimotoite and Al2O3 corundum end-members to calculate sound wave velocities for Al-bearing akimotoite at high pressure and temperature.

Geochemistry & Geophysics↗