A spectroscopic electric field vector imaging diagnostic for electron cyclotron heating systems
Not Available
Martin, E. H.↗
SEARCH · Search NASA
Results for “3D Simulation”
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.
Hyperdiffusion of dust particles in a turbulent tokamak plasma
In this work, the effect of plasma turbulence on the trajectories of dust particles is investigated for the first time. The dynamics of dust particles is computed using the ad hoc developed Dust Injection Simulator code, using a 3D turbulent plasma background computed with the TOKAM3X code. As a result, the evolution of the particle trajectories is governed by the ion drag force, and the shape of the trajectory is set by the Stokes number St∝a d /n 0 , with a d the dust radius and n 0 the density at the separatrix. The plasma turbulence is observed to scatter the dust particles, exhibiting a hyperdiffusive regime in all cases. The amplitude of the turbulent spread of the trajectories Δr 2 is shown to depend on the ratio Ku/St, with Ku∝u rms the Kubo number and urms the fluctuation level of the plasma flow. These results are compared with a simple analytical model, predicting Δr 2 ∝(Ku/St) 2 t 3 , or Δr 2 ∝(u rms n 0 /a d ) 2 t 3 . As the dust is heated by the plasma fluxes, thermionic emission sets the dust charge, originally negative, to slightly positive values. This results in a substantial reduction of the ion drag force through the suppression of its Coulomb scattering component. The dust grain inertia is then no longer negligible and drives the transition from a hyperdiffusive regime toward a ballistic one.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Ideal MHD induced temperature flattening in spherical tokamaks
This paper extends the analysis first presented in Jardin et al. [Phys. Rev. Lett. 128, 245001 (2022)] to more thoroughly examine the stability of spherical torus equilibrium to ideal magnetohydrodynamic (MHD) infernal modes and their nonlinear consequences. Here we demonstrate that in a 3D resistive magnetohydrodynamic (MHD) simulation of a NSTX discharge, anomalous transport can occur due to these instabilities. Here we generate a family of equilibrium of differing β and use this to show that these instabilities could explain the experimentally observed flattening of the electron temperature profile at modest β. The modes studied in this paper are found to occur with poloidal mode number m and toroidal mode number n when the ratio m/ n is in the range of 1.2–1.5, when the central safety factor is in this range or slightly lower, and when the central region has very low magnetic shear. Our analysis gives some insight as to why the unstable linear growth rates are oscillatory functions of the toroidal mode number n. We present a simulation of an initially stable configuration that passes through a stability boundary at a critical β as it is heated. We also show that a particular NSTX discharge is unstable to these modes over a timescale of several hundred ms. We conclude that these modes must be taken into account when performing predictive modeling. An appendix shows that similar modes can be found in R / a = 4 tokamaks for certain q-profiles and β values.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Radiatively cooled magnetic reconnection experiments driven by pulsed power
We present evidence for strong radiative cooling in a pulsed-power-driven magnetic reconnection experiment. Two aluminum exploding wire arrays, driven by a 20 MA peak current, 300 ns rise time pulse from the Z machine (Sandia National Laboratories), generate strongly driven plasma flows (MA≈7) with anti-parallel magnetic fields, which form a reconnection layer (SL≈120) at the mid-plane. The net cooling rate far exceeds the Alfvénic transit rate (τcool−1/τA−1≫1), leading to strong cooling of the reconnection layer. We determine the advected magnetic field and flow velocity using inductive probes positioned in the inflow to the layer, and inflow ion density and temperature from analysis of visible emission spectroscopy. A sharp decrease in x-ray emission from the reconnection layer, measured using filtered diodes and time-gated x-ray imaging, provides evidence for strong cooling of the reconnection layer after its initial formation. X-ray images also show localized hotspots, regions of strong x-ray emission, with velocities comparable to the expected outflow velocity from the reconnection layer. These hotspots are consistent with plasmoids observed in 3D radiative resistive magnetohydrodynamic simulations of the experiment. X-ray spectroscopy further indicates that the hotspots have a temperature (170 eV) much higher than the bulk layer (≤75 eV) and inflow temperatures (about 2 eV) and that these hotspots generate the majority of the high-energy (>1 keV) emission.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Modeling the contributions to acoustic nonlinearity from complex dislocation networks using 3D dislocation dynamics
Nonlinear ultrasonic parameters are highly sensitive to microstructural features that affect macroscale material behavior, providing a nondestructive means to characterize their evolution. Although dislocations are known to be a strong source of acoustic nonlinearity, establishing quantitative links between the acoustic nonlinearity parameter (β), measured via Second Harmonic Generation, and dislocation morphology—such as dislocation length and density—remains an open challenge. This work advances the numerical modeling of dislocation–β relationships using 3D dislocation dynamics (DD) simulations in two approaches: a “static” method computing strain and stress fields from dislocation configurations in the absence of external loading, and a “quasi-static” method to estimate β from the curvature of dislocation lines under applied load. First, the static method is combined with finite element analysis to investigate a recent assertion that heterogeneous initial strain fields can induce higher harmonic generation in a linear elastic medium; the present results do not corroborate this outcome. Then, the quasi-static method is applied to multiple-dislocation scenarios through parametric studies, revealing behaviors not predicted by analytical models, such as the competing interactions of edge and screw dislocations and the significant influence of applied stress on β. Finally, the simulations are used to model SHG experimental results and validate the hypothesis that β can decrease during plastic deformation, despite increasing dislocation density. As the DD code used here is open-source, it provides a practical platform for future investigation into microstructure–β relationships important to the interpretation of SHG results.
Materials science↗
REALISTIC MODELING OF LOCAL DYNAMO PROCESSES ON THE SUN
Magnetic fields are usually observed in the quiet Sun as small-scale elements that cover the entire solar surface (the “salt-and-pepper” patterns in line-of-sight magnetograms). By using 3D radiative MHD numerical simulations, we find that these fields result from a local dynamo action in the top layers of the convection zone, where extremely weak “seed” magnetic fields (e.g., from a 10{sup −6} G) can locally grow above the mean equipartition field to a stronger than 2000 G field localized in magnetic structures. Our results reveal that the magnetic flux is predominantly generated in regions of small-scale helical downflows. We find that the local dynamo action takes place mostly in a shallow, about 500 km deep, subsurface layer, from which the generated field is transported into the deeper layers by convective downdrafts. We demonstrate that the observed dominance of vertical magnetic fields at the photosphere and horizontal fields above the photosphere can be explained by small-scale magnetic loops produced by the dynamo. Such small-scale loops play an important role in the structure and dynamics of the solar atmosphere and their detection in observations is critical for understanding the local dynamo action on the Sun.
79 ASTRONOMY AND ASTROPHYSICS↗
Probing plasma wakefields using electron bunches generated from a laser wakefield accelerator
Here, we show experimental results of probing the electric field structure of plasma wakes by using femtosecond relativistic electron bunches generated from a laser wakefield accelerator. Snapshots of laser-driven linear wakes in plasmas with different densities and density gradients are captured. The spatiotemporal evolution of the wake in a plasma density up-ramp is recorded. Two parallel wakes driven by a laser with a main spot and sidelobes are identified in the experiment and reproduced in simulations. The capability of this new method for capturing the electron- and positron-driven wakes is also shown via 3D particle-in-cell simulations.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Topological bifurcation of magnetic islands in NSTX-U
The topological bifurcation of magnetic islands is observed in two M3D-C1 linear resistive MHD calculations of model NSTX-U equilibria during the application of small 3D magnetic perturbation fields. In comparing to a DIII-D ELM suppression case, reconstructed from shot 147170 at time slice 3745 ms, we propose that stronger kink response in the NSTX-U plasma may be the primary causes of the bifurcations. Island evolution and bifurcation in both the DIII-D and the two NSTX-U cases with the increase of perturbation current are analyzed. Further, a hypothesis is proposed that such magnetic island topology bears an intrinsic property to confine charged particles and cause their density inside the islands to rise, which is examined by implementing a random kick algorithm in the magnetic field line integration code TRIP3D-GPU to approximate electron collisional transport. Although full consideration of all forces between electrons is missing as it requires 3D large scale kinetic simulation in toroidal geometry, the positive result of this approximate transport simulation does provide first step confirmation towards fully validating the hypothesis.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Pickup Ion Effect of the Solar Wind Interaction with the Local Interstellar Medium
Pickup ions are created when interstellar neutral atoms resonantly exchange charge with the solar wind (SW) ions, especially in the supersonic part of the wind, where they carry most of the plasma pressure. Here we present numerical simulation results of the 3D heliospheric interface treating pickup ions as a separate proton fluid. To satisfy the fundamental conservation laws, we solve the system of equations describing the flow of the mixture of electrons, thermal protons, and pickup ions. To find the density and pressure of pickup ions behind the termination shock, we employ simple boundary conditions that take into account the \emph{Voyager} observations that showed that the decrease in the kinetic energy of the mixture at the termination shock predominantly contributed to the increase in the pressure of pickup ions. We show that this model adequately describes the flow of the plasma mixture and results in a noticeable decrease in the heliosheath width.
79 ASTRONOMY AND ASTROPHYSICS↗
SALPETER NORMALIZATION OF THE STELLAR INITIAL MASS FUNCTION FOR MASSIVE GALAXIES AT z ∼ 1
The stellar initial mass function (IMF) is a key parameter for studying galaxy evolution. Here we measure the IMF mass normalization for a sample of 68 field galaxies in the redshift range 0.7-0.9 within the Extended Groth Strip. To do this we derive the total (stellar + dark matter) mass-to-light [(M/L)] ratio using axisymmetric dynamical models. Within the region where we have kinematics (about one half-light radius), the models assume (1) that mass follows light, implying negligible differences between the slope of the stellar and total density profiles, (2) constant velocity anisotropy (β{sub z}≡1−σ{sub z}{sup 2}/σ{sub R}{sup 2}=0.2), and (3) that galaxies are seen at the average inclination for random orientations (i.e., i = 60°, where i = 90° represents edge-on). The dynamical models are based on anisotropic Jeans equations, constrained by Hubble Space Telescope/Advanced Camera for Surveys imaging and the central velocity dispersion of the galaxies, extracted from good-quality spectra taken by the DEEP2 survey. The population (M/L) are derived from full-spectrum fitting of the same spectra with a grid of simple stellar population models. Recent dynamical modeling results from the ATLAS{sup 3D} project and numerical simulations of galaxy evolution indicate that the dark matter fraction within the central regions of our galaxies should be small. This suggests that our derived total (M/L) should closely approximate the stellar M/L. Our comparison of the dynamical (M/L) and the population (M/L) then implies that for galaxies with stellar mass M {sub *} ≳ 10{sup 11} M {sub ☉}, the average normalization of the IMF is consistent with a Salpeter slope, with a substantial scatter. This is similar to what is found within a similar mass range for nearby galaxies.
79 ASTRONOMY AND ASTROPHYSICS↗
The role of magnetic field geometry in the evolution of neutron star merger accretion discs
ABSTRACT Neutron star mergers are unique laboratories of accretion, ejection, and r-process nucleosynthesis. We used 3D general relativistic magnetohydrodynamic simulations to study the role of the post-merger magnetic geometry in the evolution of merger remnant discs around stationary Kerr black holes. Our simulations fully capture mass accretion, ejection, and jet production, owing to their exceptionally long duration exceeding 4 s. Poloidal post-merger magnetic field configurations produce jets with energies Ejet ∼ (4–30) × 1050 erg, isotropic equivalent energies Eiso ∼ (4–20) × 1052 erg, opening angles θjet ∼ 6–13°, and durations tj ≲ 1 s. Accompanying the production of jets is the ejection of $f_\mathrm{ej}\sim 30\!-\!40{{\ \rm per\ cent}}$ of the post-merger disc mass, continuing out to times >1 s. We discover that a more natural, purely toroidal post-merger magnetic field geometry generates large-scale poloidal magnetic flux of alternating polarity and striped jets. The first stripe, of $E_\mathrm{jet}\simeq 2\times 10^{48}\, \mathrm{erg}$, Eiso ∼ 1051 erg, θjet ∼ 3.5–5°, and tj ∼ 0.1 s, is followed by ≳4 s of striped jet activity with $f_\mathrm{ej}\simeq 27{{\ \rm per\ cent}}$. The dissipation of such stripes could power the short-duration gamma-ray burst (sGRB) prompt emission. Our simulated jet energies and durations span the range of sGRBs. We find that although the blue kilonova component is initially hidden from view by the red component, it expands faster, outruns the red component, and becomes visible to off-axis observers. In comparison to GW 170817/GRB 170817A, our simulations underpredict the mass of the blue relative to red component by a factor of few. Including the dynamical ejecta and neutrino absorption may reduce this tension.
Christie, I. M.↗
From hydrogen to helium: the spectral evolution of white dwarfs as evidence for convective mixing
ABSTRACT We present a study of the hypothesis that white dwarfs undergo a spectral change from hydrogen- to helium-dominated atmospheres using a volume-limited photometric sample drawn from the Gaia-DR2 catalogue, the Sloan Digital Sky Survey (SDSS), and the Galaxy Evolution Explorer (GALEX). We exploit the strength of the Balmer jump in hydrogen-atmosphere DA white dwarfs to separate them from helium-dominated objects in SDSS colour space. Across the effective temperature range from 20 000 to 9000 K, we find that 22 per cent of white dwarfs will undergo a spectral change, with no spectral evolution being ruled out at 5σ. The most likely explanation is that the increase in He-rich objects is caused by the convective mixing of DA stars with thin hydrogen layers, in which helium is dredged up from deeper layers by a surface hydrogen convection zone. The rate of change in the fraction of He-rich objects as a function of temperature, coupled with a recent grid of 3D radiation-hydrodynamic simulations of convective DA white dwarfs – which include the full overshoot region – lead to a discussion on the distribution of total hydrogen mass in white dwarfs. We find that 60 per cent of white dwarfs must have a hydrogen mass larger than MH/MWD = 10−10, another 25 per cent have masses in the range MH/MWD = 10−14–10−10, and 15 per cent have less hydrogen than MH/MWD = 10−14. These results have implications for white dwarf asteroseismology, stellar evolution through the asymptotic giant branch and accretion of planetesimals on to white dwarfs.
Cunningham, Tim↗
Thermal and kinematic properties of ejecta in SN1987A revealed by XRISM
We present an analysis of high-resolution spectra from the shock-heated plasmas in SN 1987A, based on an observation using the Resolve instrument onboard the X-Ray Imaging and Spectroscopy Mission (XRISM). The 1.7–10 keV Resolve spectra are accurately represented by a single-component, plane-parallel shock plasma model, with a temperature of $2.84_{-0.08}^{+0.09}$ keV and an ionization parameter of $2.64_{-0.45}^{+0.58}$ × $10^{11}\,\,{\rm s\,\, cm}^{-3}$. The Resolve spectra are also well reproduced by the 3D magneto-hydrodynamic simulation presented by Orlando et al. (2020, A&A, 636, A22) suggesting substantial contribution from the ejecta. The metal abundances obtained with Resolve align with the Large Magellanic Cloud value, indicating that the X-rays in 2024 originate from “non-metal-rich” shock-heated ejecta and the reverse shock has not reached the inner metal-rich region of ejecta. Doppler widths of the atomic lines from Si, S, and Fe correspond to velocities of 1500–1700 km s$^{-1}$, where the thermal broadening effects in this non-metal-rich plasma are negligible. Therefore, the line broadening seen in Resolve spectra is determined by the large bulk motion of ejecta. For reference, we determined a $90\%$ upper limit on non-thermal emission from a pulsar wind nebula at $4.3 \times 10^{-13}$ erg cm$^{-2}$ s$^{-1}$ in the 2–10 keV range, aligning with NuSTAR findings by Greco et al. (2022, ApJ, 931, 132). Additionally, we searched for the $^{44}$Sc K line feature and found a $1\sigma$ upper limit of $1.0 \times 10^{-6}$ photons cm$^{-2}$ s$^{-1}$, which translates to an initial $^{44}$Ti mass of approximately $2 \times 10^{-4}\, M_{\odot }$, consistent with previous X-ray to soft gamma-ray observations (Boggs et al. 2015, Science, 348, 670; Grebenev et al. 2012, Nature, 490, 373; Leising 2006, ApJ, 651, 1019).
ISM: supernova remnants↗
3D multiphysics modeling of superconducting cavities with a massively parallel simulation suite
Not Available
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Self-Organized Stationary States of Tokamaks
We demonstrate that in a 3D resistive magnetohydrodynamic (MHD) simulation, for some parameters it is possible to form a stationary state in a tokamak where a saturated interchange mode in the center of the discharge drives a near helical flow pattern that acts to non-linearly sustain the configuration by adjusting the central loop voltage through a dynamo action. This could explain the physical mechanism for maintaining stationary non-sawtoothing “hybrid” discharges, often referred to as “flux-pumping”.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Delocalization Transition of a Disordered Axion Insulator
The axion insulator is a higher-order topological insulator protected by inversion symmetry. We show that, under quenched disorder respecting inversion symmetry on average, the topology of the axion insulator stays robust, and an intermediate metallic phase in which states are delocalized is unavoidable at the transition from an axion insulator to a trivial insulator. We derive this conclusion from general arguments, from classical percolation theory, and from the numerical study of a 3D quantum network model simulating a disordered axion insulator through a layer construction. We find the localization length critical exponent near the delocalization transition to be ν = 1.42 ± 0.12. We further show that this delocalization transition is stable even to weak breaking of the average inversion symmetry, up to a critical strength. Finally, we also quantitatively map our quantum network model to an effective Hamiltonian and we find its low-energy k ∙ p expansion.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Geostatistical Mapping of Salinity Conditioned on Borehole Logs, Montebello Oil Field, California
We present a geostatistics-based stochastic salinity estimation framework for the Montebello Oil Field that capitalizes on available total dissolved solids (TDS) data from groundwater samples as well as electrical resistivity (ER) data from borehole logging. Data from TDS samples (n = 4924) was coded into an indicator framework based on falling below four selected thresholds (500, 1000, 3000, and 10,000 mg/L). Collocated TDS-ER data from the surrounding groundwater basin were then employed to produce a kernel density estimator to establish conditional probabilities for ER data (n = 8 boreholes) falling below the selected TDS thresholds within the Montebello Oil Field area. Directional variograms were estimated from these indicator coded data, and 500 TDS realizations from conditional indicator simulation were generated for the subsurface region above the Montebello Oil Field reservoir. Simulations were summarized as 3D maps of median TDS, most likely salinity class, and probability for exceeding each of the specified TDS thresholds. Results suggested TDS was below 500 mg/L in most of the study area, with a trend toward higher values (500 to 1000 mg/L) to the southwest; consistent with the average regional groundwater flow direction. Discrete localized zones of TDS greater than 1000 mg/L were observed, with one of these zones in the greater than 10,000 mg/L range; however, these areas were not prevalent. The probabilistic approach used here is adaptable and is readily modified to include additional data and types and can be employed in time-lapse salinity modeling through Bayesian updating.
54 ENVIRONMENTAL SCIENCES↗
Comprehensive computer model for magnetron sputtering. I. Gas heating and rarefaction
The complex interaction between several variables in magnetron sputtering discharges is a challenge in developing engineering design tools for industrial applications. For instance, at high pressures, rarefaction and gas heating should no longer be neglected for determining several parameters of the process. In this article, we use a comprehensive 3D reactor-scale simulator that incorporates most phenomena of interest in a self-consistent manner to simulate the transport of sputtered particles over a wide range of pressures and powers. Calculations of aluminum deposition rates and metal vapor densities are in reasonable agreement with experiments over a wide range of pressures and powers. Of the elements investigated (Al, Ti, and Cu), copper showed the greatest rarefaction (30%) due to its higher sputtering yield. Titanium, despite a slightly lower sputtering yield than Al, shows a greater rarefaction than aluminum as more particles are reflected from the target as high energy neutrals. In this case, a more efficient energy transfer process is responsible for the higher rarefaction observed in Ti sputtering when compared to Al. The authors also observed that by sputtering at a higher pressure, the probability of electron impact ionization of sputtered particles is increased and speculate about the role of this process in contrast to penning ionization, which is believed to be the dominant ionization mechanism in magnetron sputtering.
36 MATERIALS SCIENCE↗