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

Degradation of performance in ICF implosions due to Rayleigh–Taylor instabilities: A Hamiltonian perspective

The Rayleigh–Taylor instability (RTI) is an ubiquitous phenomenon that occurs in inertial-confinement-fusion (ICF) implosions and is recognized as an important limiting factor of ICF performance. To analytically understand the RTI dynamics and its impact on ICF capsule implosions, we develop a first-principle variational theory that describes an imploding spherical shell undergoing RTI. The model is based on a thin-shell approximation and includes the dynamical coupling between the imploding spherical shell and an adiabatically compressed fluid within its interior. Using a quasilinear analysis, we study the degradation trends of key ICF performance metrics (e.g., stagnation pressure, residual kinetic energy, and areal density) as functions of initial RTI parameters (e.g., the initial amplitude and Legendre mode), as well as the 1D implosion characteristics (e.g., the convergence ratio). We compare analytical results from the theory against nonlinear results obtained by numerically integrating the governing equations of this reduced model. Our findings emphasize the need to incorporate polar flows in the calculation of residual kinetic energy and demonstrate that higher convergence ratios in ICF implosions lead to significantly greater degradation of key performance metrics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Second harmonic generation in topological insulators under quantizing magnetic fields

We theoretically investigate the second harmonic generation (SHG) of topological insulator surface states in a perpendicular magnetic field. Our theory is based on the microscopic expression of the second-order magneto-optical conductivity developed from the density matrix formalism, taking into account hexagonal warping effects on the surface states' band structure. Using numerically exact Landau-level energies and wave functions including hexagonal warping, we calculate the spectrum of SHG conductivities under normal incidence for different values of magnetic field and chemical potential. The imaginary parts of the SHG conductivities show prominent resonant peaks corresponding to one-photon and two-photon inter-Landau-level transitions. Treating the hexagonal warping term perturbatively, these transitions are clarified analytically within a perturbation theory from which approximate selection rules for the allowable optical transitions for SHG are determined. Furthermore, our results show extremely high SHG susceptibility that is easily tunable by magnetic field and doping level for topological surface states in the far-infrared regime, exceeding that of many conventional nonlinear materials. This work highlights the key role of hexagonal warping effects in generating second-order optical responses and provides insights into the nonlinear magneto-optical properties of the topological insulators.

Landau levels↗
DOE OSTI2024

High-resolution particle-in-cell simulations of two-dimensional Bernstein–Greene–Kruskal modes

We present two-dimensional (2D) particle-in-cell (PIC) simulations of 2D Bernstein–Greene–Kruskal modes, which are exact nonlinear steady-state solutions of the Vlasov–Poisson equations, on a 2D plane perpendicular to a background magnetic field, with a cylindrically symmetric electric potential localized on the plane. PIC simulations are initialized using analytic electron distributions and electric potentials from the theory. We confirm the validity of such solutions using high-resolutions up to a 20482 grid. We show that the solutions are dynamically stable for a stronger background magnetic field, while keeping other parameters of the model fixed, but become unstable when the field strength is weaker than a certain value. When a mode becomes unstable, we observe that the instability begins with the excitation of azimuthal electrostatic waves that ends with a spiral pattern.

Physics↗
DOE OSTI2025

Wide-Bandgap Semiconductor Amplifiers for Fusion Plasma Heating and Control

This paper discusses power electronics developed under the ARPA-E GAMOW program to support nuclear fusion power production. The goal of this project was to develop and assess the potential for wide-bandgap (WBG) semiconductor devices in power electronics to enable high-efficiency and high-voltage solid-state systems for fusion plasma generation, heating, and control. The power electronics use an architecture in which multiple high-power boards can be combined to produce megawatt-level power, where using multiple boards provides high reliability. Two main areas of power electronics boards are developed in this project for fusion plasma heating and control applications: (1) pulse generation and control and (2) radiofrequency generation. The first area is for boards capable of driving high-voltage millisecond pulses at high duty cycles. The envisioned application of these pulses is in plasma control of magnetohydrodynamic instabilities, plasma position, and edge-localized modes. Pulse-width modulation allows for the implementation of a wide variety of linear and nonlinear control systems. The boards developed for this project could actuate control coils based on digital input signals and can be parallelized to provide megawatts of output power. The design of the pulse generator is a low-side load switch. A load switch was designed and constructed that utilized 2-kV-rated field-effect transistor (FET)-based cascodes developed by Qorvo under this project to perform initial testing of these cascodes. The second area is being implemented using class E amplifiers with WBG devices and a reactance steering network to handle inductive or capacitive plasma loads. Applications include ion cyclotron resonance heating (ICRH) and high-harmonic fast-wave (HHFW) heating. A class E reactance steering network is demonstrated in modeling and experiment with a resistive-inductive load that models an inductively-coupled plasma. Power combining of boards with class E reactance steering networks is also simulated and demonstrated experimentally, to enable scaling up to high power. Modeling of high-power-density cooling and remaining useful life is conducted to enable reliable, effectively cooled high-power electronics for fusion applications.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Protected fermionic zero modes in periodic gauge fields

It is well known that macroscopically normalizable zero-energy wave functions of spin-$\frac{1}{2}$ particles in a two-dimensional inhomogeneous magnetic field are spin-polarized and exactly calculable with degeneracy equaling the number of flux quanta linking the whole system. Here, extending this argument to massless Dirac fermions subjected to magnetic fields that have zero net flux but are doubly periodic in real space, we show that there exist only two Bloch-normalizable zero-energy eigenstates, one for each spin flavor. This result is immediately relevant to graphene multilayer systems subjected to doubly periodic strain fields, which at low energies enter the Hamiltonian as periodic pseudogauge vector potentials. Furthermore, we explore various related settings including nonlinearly dispersing band structure models and systems with singly periodic magnetic fields.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Second harmonic Bessel-Gauss beam shaping with elliptic axicon aberrations

The second harmonic (SH) of an axicon generated Bessel-Gauss beam is created through the nonlinear interaction of photons with a crystal, resulting in the energy doubling of the output photons. In this work, we show experimentally that in addition to frequency doubling, the SH of Bessel-Gauss beams under asymmetric aberrations from an elliptic axicon exhibit intriguing beam formation. Particularly, the central region of the SH beam profile is composed of two central spots of various geometries surrounded by nested ellipses; one of which is the configuration of two central gamma dots with similar radius knotted by nested ellipses for a zeroth-order Bessel-Gauss pump. These SH beams consistently maintain their spatial profile throughout propagation, reminiscing the behavior of screw dislocations in wave patterns. Our numerical simulations produce beam dynamics consistent with that of experiments and further implicate the remarkable interweaving of bright spots with dark vortices. This is especially noticeable when the beams dynamically oscillate along the optical axis, resulting in the genesis of spatially polarized beams with a knotted framework. The insights gained from our study establish a novel paradigm for exploring interactions of Bessel-like beams with vortex dynamics. This, in turn, has the potential to spark innovations in optical applications, fostering new methodologies to harness and manipulate complex light structures. Our experimental findings could spur new methods of generating logical states of light and new opportunities for material processing control. Published by the American Physical Society 2025

Wang, Tianhong (ORCID:0000000169104201)↗
DOE OSTI2022

Estimating epistemic uncertainty in soil parameters for nonlinear site response analyses: Introducing the Latin Hypercube Sampling technique

This study quantifies the effects of epistemic uncertainty in soil parameters on nonlinear (NL) site response analysis (SRA) results, validated against the data recorded at a well-instrumented geotechnical downhole array located in Japan. To this end, a one-dimensional soil column model of the Service Hall Array (SHA) near the Kashiwazaki-Kariwa Nuclear Power Plant (KKNPP) is developed using the finite element (FE) program LS-DYNA. The dynamic stress–strain relationship is characterized by a modified two-stage hyperbolic (MTH) NL backbone curve formulation capable of capturing soil behavior at both small- and large-shear strains. The model is then validated against the ground motion recordings to capture the model bias. The uncertainties associated with the shear-wave velocity profile (a small-strain soil property) and soil shear strength (a large-strain soil property) are incorporated in NL SRA to quantify their separate and joint randomization effects on the results. This study proposes using the Latin Hypercube Sampling (LHS) method as an efficient alternative to commonly used methods, such as Standard Monte Carlo (SMC), to account for uncertainty propagation in such reliability analysis. Both low-intensity and design-level records from the recordings at the SHA are applied to study the contribution of the small- and large-strain NL dynamic soil properties. Results from 46,200 NL FE analyses (23,100 per input motion) are presented. Measured and predicted site response, using recorded ground motions at this downhole array, is compared to assess the significance of soil parameter uncertainty on the observed ground motion dispersions. It is demonstrated that increasing the number of soft realizations and implementing higher level earthquake intensity lead to higher ground motion dispersion. Unlike past studies in randomization of Vs profiles with the SMC method, the LHS method is shown to have no significant effect on the predicted median surface response spectra and amplification factors (AFs) for this case study.

Engineering↗
DOE OSTI2020

Pattern formation in a coupled membrane-bulk reaction-diffusion model for intracellular polarization and oscillations

Reaction-diffusion systems have been widely used to study spatio-temporal phenomena in cell biology, such as cell polarization. Coupled bulk-surface models naturally include compartmentalization of cytosolic and membrane-bound polarity molecules. Here we study the distribution of the polarity protein Cdc42 in a mass-conserved membrane-bulk model, and explore the effects of diffusion and spatial dimensionality on spatio-temporal pattern formation. We first analyze a one-dimensional (1-D) model for Cdc42 oscillations in fission yeast, consisting of two diffusion equations in the bulk domain coupled to nonlinear ODEs for binding kinetics at each end of the cell. In 1-D, our analysis reveals the existence of symmetric and asymmetric steady states, as well as anti-phase relaxation oscillations typical of slow-fast systems. We then extend our analysis to a two-dimensional (2-D) model with circular bulk geometry, for which species can either diffuse inside the cell or become bound to the membrane and undergo a nonlinear reaction-diffusion process. We also consider a nonlocal system of PDEs approximating the dynamics of the 2-D membrane-bulk model in the limit of fast bulk diffusion. In all three model variants we find that mass conservation selects perturbations of spatial modes that simply redistribute mass. In 1-D, only anti-phase oscillations between the two ends of the cell can occur, and in-phase oscillations are excluded. In higher dimensions, no radially symmetric oscillations are observed. Instead, the only instabilities are symmetry-breaking, either corresponding to stationary Turing instabilities, leading to the formation of stationary patterns, or to oscillatory Turing instabilities, leading to traveling and standing waves. Codimension-two Bogdanov—Takens bifurcations occur when the two distinct instabilities coincide, causing traveling waves to slow down and to eventually become stationary patterns. Our work clarifies the effect of geometry and dimensionality on behaviors observed in ma.ss-conserved cell polarity models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Atmospheric clearing of solid-particle debris using femtosecond filaments

Through nonlinear self-focusing, femtosecond pulses can propagate several kilometers beyond diffraction limits, forming an ionization channel in air known as a laser filaments. It has been demonstrated that in the wake of the filament, aerosols can be effectively cleared to improve the transmission of subsequent laser pulses or secondary light sources, pertinent to applications in atmospheric sensing. However, the current understanding of aerosol clearing is founded on interactions with droplets to simulate fogs and clouds and thus does not extend to solid particles or atmospheric debris. Using optical trapping, we isolate both graphite and silica microparticles and directly measure the subsequent displacement caused by the filament using time-resolved shadowgraphy. The shock wave from the filament is demonstrated to propel particles away from the filament, directly contributing to atmospheric debris clearing. Particles exposed to the laser light in either the intense filament core or the surrounding energy reservoir are axially displaced along the beam path. Furthermore, it is found that the optomechanical properties of the particle largely influence the axial displacement induced by laser exposure through mechanisms such as radiation pressure, mass ejection from ablation or optical damage, and particle deagglomeration.

42 ENGINEERING↗
DOE OSTI2021

Visualization of Fast Ion Phase-Space Flow Driven by Alfvén Instabilities

Fast ion phase-space flow, driven by Alfven eigenmodes (AEs), is measured by an imaging neutral particle analyzer in the DIII-D tokamak. The flow firstly appears near the minimum safety factor at the injection energy of neutral beams, and then moves radially inward and outward by gaining and losing energy, respectively. The flow trajectories in phase space align well with the intersection lines of the constant magnetic moment surfaces and constant E – (ω/n)P ζ surfaces, where E, P ζ are energy and toroidal canonical momentum of ions; ω and n are angular frequencies and toroidal mode numbers of AEs. It is found that the flow is so destructive that the thermalization of fast ions is no longer observed in regions of strong interaction. Here, the measured phase-space flow is consistent with nonlinear hybrid kinetic-magnetohydrodynamics simulation. Calculations of the relatively narrow phase-space islands reveal that fast ions must transition between different flow trajectories to experience large-scale phase-space transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Coupling magnons to an opto-electronic parametric oscillator

Hybrid magnonic systems have emerged as versatile modular components for quantum signal transduction and sensing applications owing to their capability of connecting distinct quantum platforms. To date, the majority of the magnonic systems have been explored in a local, near-field scheme, due to the close proximity required for realizing a strong coupling between magnons and other excitations. This constraint greatly limits the applicability of magnons in developing remotely coupled, distributed quantum network systems. On the contrary, opto-electronic architectures hosting self-sustained oscillations have been a unique platform for long-haul signal transmission and processing. Here, we integrated an opto-electronic oscillator with a magnonic oscillator consisting of a microwave waveguide and a Y 3 Fe 5 O 12 (YIG) sphere, and demonstrated strong and coherent coupling between YIG’s magnon modes and the opto-electronic oscillator’s characteristic photon modes—revealing the hallmark anti-crossing gap in the measured spectrum. In particular, the photon mode is produced on-demand via a nonlinear, parametric process as stipulated by an external seed pump. Both the internal cavity phase and the external pump phase can be precisely tuned to stabilize either degenerate or nondegenerate auto-oscillations. Our result lays out a new, hybrid platform for investigating the long-distance coupling and nonlinearity in coherent magnonic phenomena, which may be found useful in constructing the future “distributed hybrid magnonic systems.”

36 MATERIALS SCIENCE↗
DOE OSTI2022

Proper orthogonal decomposition based reduced-order modeling of flux-Limited gray thermal radiation

Here, in this work, a proper orthogonal decomposition (POD) based reduced-order model (ROM) is developed to solve gray, flux-limited thermal radiation diffusion. We focus on the variable opacity radiation penetration benchmark posed by Olson, Auer, and Hall. The T -3 relationship for opacity in conjunction with high-temperature radiation penetrating an initially cold material produces a strong thermal radiation shock. This class of problems is particularly challenging for standard POD-based reduced-order modeling due to the nonlinearities presented by 1) the T 4 source term, and 2) flux-limited diffusion operator. To address these challenges and develop a cost competitive ROM, we employ a “hyper-reduction” technique through discrete empirical interpolation (DEIM) and allow for adaptive reduced-order projections through principal interval decomposition (PID). Performance of the proposed methodology is quantified by comparing the cost savings and accuracy relative to a full-order computation. Reference solutions and snapshot data are obtained through a full-order calculation performed by the University of Chicago maintained astrophysics code, FLASH. For consistency and potential extensibility, the developed ROM is also implemented in FLASH. We find that in the initialization regime, where the thermal radiation wave is initially created by the warming of the material, this class of problems is highly reducible and suitable for POD-based ROMs. However, the strong convective nature of the wave propagation regime is less reducible and more challenging to create an efficient ROM.

42 ENGINEERING↗
DOE OSTI2026

Fokker-Planck simulations of fast ion ICRF and electron EC heating in a mirror plasma using CQL3D-m

The CQL3D-m continuum bounce-average Fokker-Planck code is adapted for magnetic mirror plasmas [1] and is now routinely used in no-free-parameter classical integrated modeling of mirror devices [2, 3]. In the present effort, we report on two RF methods of plasma heating in mirror machine. The fast ions (FI) are heated by Fast waves at 2nd-4th harmonic, where FIs originate from neutral beam injection at 45 degrees to the magnetic field. The scenario shows an efficient ion heating near the FI bouncing point. The electrons are heated by X-mode launched from the high magnetic field side towards the resonance. Different from the tokamak applications, CQL3D-m provides an evolving self-consistent ambipolar parallel electric field, which determines the shape of the loss cone and hence an accurate confinement time of both ions and electrons. Also, it includes a description of ion and electron sources and sinks (related to charge exchange and impact ionization) which are updated at every time step. CQL3D-m utilizes a fully nonlinear Coulomb collision operator that is important for the significantly non-Maxwellian ion distributions typically established in mirror plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Don't cross the streams: caustics from fuzzy dark matter

Here we study how tidal streams from globular clusters may be used to constrain the mass of ultra-light dark matter particles, called `fuzzy' dark matter (FDM). A general feature of FDM models is the presence of ubiquitous density fluctuations in bound, virialized dark matter structures, on the scale of the de Broglie wavelength, arising from wave interference in the evolving dark matter distribution. These time-varying fluctuations can disturb the motions of stars, leading to potentially observable signatures in cold thin tidal streams in our own Galaxy. The study of this effect has been hindered by the difficulty in simulating the FDM wavefunction in Milky Way-sized systems. We present a simple method to evolve realistic wavefunctions in nearly static potentials, that should provide an accurate estimate of this granulation effect. We quantify the impact of FDM perturbations on tidal streams, and show that initially, while stream perturbations are small in amplitude, their power spectra exhibit a sharp cutoff corresponding to the de Broglie wavelength of the FDM potential fluctuations. Eventually, when stream perturbations become nonlinear, fold caustics generically arise that lead to density fluctuations with universal behavior. This erases the signature of the de Broglie wavelength in the stream density power spectrum, but we show that it will still be possible to determine the FDM mass in this regime, by considering the fluctuations in quantities like angular momenta or actions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Develop and verify soil/structure interaction for pile/foundation interaction

Phase II of the Offshore Code Comparison Collaboration, Continued, with Correlation and unCertainty (OC6) project was used to verify the implementation of a new soil-structure interaction (SSI) model for use within offshore wind turbine modeling software. The REDWIN Macro-element model implemented and verified in this study enables a computationally efficient way to model the linear and nonlinear SSI problem, including hysteretic damping, of a monopile structure. The modeling approach was integrated into several modeling tools and a series of increasingly complex simulations was conducted using the IEA 10MW reference turbine mounted on a monopile support structure to verify the coupling between the tools and the REDWIN Macro-element SSI model. This campaign includes only numerical verification between various software and modeling approaches so no experimental measurements are available. The load cases (LC) considered include: LC1 – static response of the tower and substructure LC2 – frequency and mode-shape analysis of the tower and substructure LC3 – response of the tower and substructure due to wind-only loading LC4 – response of the tower and substructure due to wave-only loading LC5 – response of the tower and substructure due to wind and wave loading. Detailed properties of the modeled system are found in the following reference, “Bergua, Roger, Amy Robertson, Jason Jonkman, and Andy Platt. 2021. "Specification Document for OC6 Phase II: Verification of an Advanced Soil-Structure Interaction Model for Offshore Wind Turbines.” Golden, CO: National Renewable Energy Laboratory. NREL/TP-5000-79938. https://www.nlr.gov/docs/fy21osti/79938.pdf. Details on the results from the OC6 Phase II project can be found in the following reference, “Bergua R, Robertson A, Jonkman J, et al. OC6 Phase II: Integration and verification of a new soil–structure interaction model for offshore wind design.” Wind Energy. 2022;25(5):793-810. doi:10.1002/we.2698

17 WIND ENERGY↗
DOE OSTI2025

Scaling up the transcorrelated density matrix renormalization group

Explicitly correlated methods, such as the transcorrelated method which shifts a Jastrow or Gutzwiller correlator from the wave function to the Hamiltonian, are designed for high-accuracy calculations of electronic structures, but their application to larger systems has been hampered by the computational cost. We develop improved techniques for the transcorrelated density-matrix renormalization group (DMRG), in which the ground state of the transcorrelated Hamiltonian is represented as a matrix product state (MPS), and demonstrate large-scale calculations of the ground-state energy of the two-dimensional Fermi-Hubbard model. Our developments stem from three technical inventions: (i) constructing matrix product operators (MPOs) of transcorrelated Hamiltonians with low bond dimension and high sparsity, (ii) exploiting the entanglement structure of the ground states to increase the accuracy of the MPS representation, and (iii) optimizing the nonlinear parameter of the Gutzwiller correlator to mitigate the nonvariational nature of the transcorrelated method. Here, we examine systems of size up to 12×12 lattice sites, four times larger than previous transcorrelated DMRG studies, and demonstrate that transcorrelated DMRG yields significant improvements over standard nontranscorrelated DMRG for equivalent computational effort. Transcorrelated DMRG reduces the error of the ground-state energy by 2.4×–14×, with the smallest improvement seen for a small system at half filling and the largest improvement in a dilute closed-shell system.

Density matrix renormalization group↗
DOE OSTI2023

Transient gap generation in BaFe2As2 driven by coherent lattice vibrations

Abstract Iron-based superconductors provide a rich platform to investigate the interplay between unconventional superconductivity, nematicity, and magnetism. The electronic structure and the magnetic properties of iron-based superconductors are highly sensitive to the pnictogen height. Coherent excitation of the A1g phonon by femtosecond laser directly modulates the pnictogen height, which has been used to control the physical properties of iron-based superconductors. Previous studies show that the driven A1g phonon resulted in a transient increase of the pnictogen height in BaFe2As2, favoring an enhanced Fe magnetic moment. However, there are no direct observations on either the enhanced Fe magnetic moments or the enhanced spin-density wave (SDW) gap. Here, we use time-resolved broadband terahertz spectroscopy to investigate the dynamics of BaFe2As2 in the A1g phonon-driven state. Below the SDW transition temperature, we observe a transient gap generation at early-time delays. A similar transient feature is observed in the normal state up to room temperature.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

The Impact of Anomalous Resistivity in Vacuum Contaminant Plasmas on the Electrothermal Instability

This manuscript presents an assessment of the electrothermal instability (ETI) in the presence of anomalous resistivity (AR) in vacuum contaminant plasmas (VCP) when applied to a magnetized liner inertial fusion (MagLIF)-like load.Pulsed-power driven dielectrically coated metallic liners, like in MagLIF, experience the current-driven electrothermal instability which occurs when a material’s resistivity changes with temperature and is subject to ohmic heating. Large scale pulsed-power facilities that use magnetically insulated transmission lines (MITL) have been shown to generate low-density plasma which enters the target chamber and coalesces around the load. The low-density high-temperature vacuum contaminant plasmas (VCP) can parasitically divert current from the load through causing a short in the anode-cathode gap inside the target chamber. Resistive magnetohydrodynamic (MHD) simulations of these VCP experience unphysical runaway ohmic heating due to under predicting the resistivity by using a purely collisional resistivity model.AR provides a physics-based way to address this runaway heating through increasing the resistivity in a proportional way with the drift speed. In this work, 1D simulations probe the effect that AR in VCP has on the magnetic diffusion rate, and 2D simulations show how this effect manifests in the nonlinear striation form of the ETI for a MagLIF-like load. Beryllium and aluminum dielectrically coated liners are used for the 1D and 2D simulations in this work. The1D simulations show that a VCP causes a delay in the current delivery to the load by upwards of 8 ns at 60 ns into atypical current pulse for an 1 × 10 7 A scale pulsed-power accelerator. The 2D simulations show the delay observed in the 1D simulations is reduced substantially (4 ns) by 2D turbulence that disrupts the VCP layer early in time. Here, the 2Dsimulated ETI growth varies across AR models, more so for the beryllium liner than the aluminum liner because the beryllium liner shows an enhanced rate of penetration for the magnetic diffusion wave in comparison to the aluminum liner. The 2D simulations show the bulk dielectric thickness varies across AR models with the Davidson AR model beingthe largest and the Buneman AR model being the smallest, and in connection with the thickness the Rayleigh-Taylor bubble-spike distances varies correspondingly.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗