Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MAGNETIC COMPRESSION”

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 1,567 records · Page 87
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Unraveling sources of emission heterogeneity in Silicon Vacancy color centers with cryo-cathodoluminescence microscopy

Diamond color centers have proven to be versatile quantum emitters and exquisite sensors of stress, temperature, electric and magnetic fields, and biochemical processes. Among color centers, the silicon-vacancy (SiV - ) defect exhibits high brightness, minimal phonon coupling, narrow optical linewidths, and high degrees of photon indistinguishability. Yet the creation of reliable and scalable SiV - -based color centers has been hampered by heterogeneous emission, theorized to originate from surface imperfections, crystal lattice strain, defect symmetry, or other lattice impurities. Here, we advance high-resolution cryo-electron microscopy combined with cathodolumines cence spectroscopy and 4D scanning transmission electron microscopy (STEM) to elucidate the structural sources of heterogeneity in SiV - emission from nanodiamond with sub-nanometer-scale resolution. Our diamond nanoparticles are grown directly on TEM membranes from molecular-level seedings, representing the natural formation conditions of color centers in diamond. We show that individual subcrystallites within a single nanodiamond exhibit distinct zero-phonon line (ZPL) energies and differences in brightness that can vary by 0.1 meV in energy and over 70% in brightness. These changes are correlated with the atomic-scale lattice structure. We find that ZPL blue shifts result from tensile strain, while ZPL red shifts are due to compressive strain. We also find that distinct crystallites host distinct densities of SiV - emitters and that grain boundaries impact SiV - emission significantly. Finally, we interrogate nanodiamonds as small as 40 nm in diameter and show that these diamonds exhibit no spatial change to their ZPL energy. Our work provides a foundation for atomic-scale structure-emission correlation, e.g., of single atomic defects in a range of quantum and two-dimensional materials.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Pressure-enhanced interplay between lattice, spin, and charge in the mixed perovskite La 2 FeMnO 6

Spin crossover plays a central role in the structural instability, net magnetic moment modification, metallization, and even in superconductivity in corresponding materials. Most reports on the pressure-induced spin crossover with a large volume collapse have so far focused on compounds with a single transition metal. Here we report a comprehensive high-pressure investigation of a mixed Fe-Mn perovskite La 2 FeMnO 6 . Under pressure, the strong coupling between Fe and Mn leads to a combined valence/spin transition: Fe 3+ (S=5/2) → Fe 2 + (S = 0) and Mn 3+ (S = 2) → Mn 4+ (S = 3/2), with an isostructural phase transition. The spin transitions of both Fe and Mn are offset by ~20 GPa of the onset pressure, and the lattice collapse occurs in between. Interestingly, Fe 3+ ion shows an abnormal behavior when it reaches a lower valence state (Fe 2+ ) accompanied by a +0.5 eV energy shift in the Fe K -absorption edge at 15 GPa. This process is associated with the charge-spin-orbital state transition from high spin Fe 3+ to low spin Fe 2+ , caused by the significantly enhanced t 2 g -e g crystal field splitting in the compressed lattice under high pressure. Density functional theory calculations confirm the energy preference of the high-pressure state with charge redistribution accompanied by spin state transition of Fe ions. Moreover, La 2 FeMnO 6 maintains semiconductor behaviors even when the pressure reached 144.5 GPa as evidenced by the electrical transport measurements, despite the huge resistivity decreasing seven orders of magnitude compared with that at ambient pressure. The investigation carried out here demonstrates high flexibility of double perovskites and their good potentials for optimizing the functionality of these materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Salvage Stereotactic Body Radiotherapy (SBRT) Following In-Field Failure of Initial SBRT for Spinal Metastases

Purpose: We report our experience in salvaging spinal metastases initially irradiated with stereotactic body radiation therapy (SBRT), who subsequently progressed with imaging-confirmed local tumor progression, and were re-irradiated with a salvage second SBRT course to the same level. Methods and Materials: From a prospective database, 56 metastatic spinal segments in 40 patients were identified as having been irradiated with a salvage second SBRT course to the same level. In addition, 24 of 56 (42.9%) segments had initially been irradiated with conventional external beam radiation therapy before the first course of SBRT. Local control (LC) was defined as no progression on magnetic resonance imaging at the treated segment, and calculated according to the competing risk model. Overall survival (OS) was evaluated for each patient treated by use of the Kaplan-Meier method. Results: The median salvage second SBRT total dose and number of fractions was 30 Gy in 4 fractions (range, 20-35 Gy in 2-5 fractions), and for the first course of SBRT was 24 Gy in 2 fractions (range, 20-35 Gy in 1-5 fractions). The median follow-up time after salvage second SBRT was 6.8 months (range, 0.9-39 months), the median OS was 10.0 months, and the 1-year OS rate was 48%. A longer time interval between the first and second SBRT courses predicted for better OS (P=.02). The crude LC was 77% (43/56), the 1-year LC rate was 81%, and the median time to local failure was 3.0 months (range, 2.7-16.7 months). Of the 13 local failures, 85% (11/13) and 46% (6/13) showed progression within the epidural space and paraspinal soft tissues, respectively. Absence of baseline paraspinal disease predicted for better LC (P<.01). No radiation-induced vertebral compression fractures or cases of myelopathy were observed. Conclusion: A second course of spine SBRT, most often with 30 Gy in 4 fractions, for spinal metastases that failed initial SBRT is a feasible and efficacious salvage treatment option.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2026

Development of a broadband hard x-ray radiography platform for pulsed-power experiments

In this article, we develop and demonstrate a broadband hard x-ray radiography platform at the Zebra Pulsed Power Laboratory that integrates point-projection radiography, bremsstrahlung measurements, and hard x-ray pinhole imaging, designed to diagnose current-driven, cylindrically compressed matter. Initial laser-pulsed-power coupled experiments revealed that intense background radiation generated during 1 MA Zebra current shots overwhelmed laser-produced hard x-rays, obscuring radiographic images. Using combined spectral and spatial diagnostics, we identify energetic electrons accelerated by return currents as the dominant source of background hard x-rays, with electron energies inferred to be 3–4 MeV based on Monte Carlo simulations, and demonstrate mitigation through modifications to the radiation shielding and return-current configuration. The diagnostic platform was validated using a wire-pinch hard x-ray source, allowing radiographs of static 1-mm-diameter aluminum wires to be obtained while simultaneously measuring x-ray source spectra and spatial emission distributions within a single shot. Measured wire transmission profiles were quantitatively reconstructed using radiation transport simulations that incorporate an experimentally inferred two-temperature exponential x-ray spectrum from bremsstrahlung signal analysis and spatially distributed emission sources identified by pinhole imaging. Agreement between measured and simulated transmission profiles demonstrates the validity of the radiographic and x-ray source characterization approach, establishing this diagnostic platform as a promising tool for diagnosing magnetically driven, high-density plasmas relevant to warm dense matter and inertial fusion energy research.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2011

Magnetic Field Strength Fluctuations in the Heliosheath: Voyager 1 Observations during 2009

We analyze the "microscale fluctuations" of the magnetic field strength Beta on a scale of several hours observed by Voyager I (VI) in the heliosheath during 2009. The microscale fluctuations of Beta range from coherent to stochastic structures. The amplitude of microscale fluctuations of Beta during 1 day is measured by the standard deviation (SD) of 48 s averages of B. The distribution of the daily values of SD is lognormal. SD(t) from day of year (DOY) I to 331, 2009, is very intermittent. SD(t) has a 1/f or "pink noise" spectrum on scales from I to 100 days, and it has a broad multi fractal spectrum f(alpha) with 0.57 much < alpha much < 1.39. The time series of increments SD(t + r) - SD(t) has a pink noise spectrum with alpha' = 0.88 +/- 0.14 on scales from 1 to 100 days. The increments have a Tsallis (q-Gaussian) distribution on scales from 1 to 165 days, with an average q = 1.75 +/- 0.12. The skewness S and kurtosis K have Gaussian and lognormal distributions, respectively. The largest spikes in K(t) and Set) are often associated with a change in Beta across a data gap and with identifiable physical structures. The "turbulence" observed by VI during 2009 was weakly compressible on average but still very intermittent, highly variable, and highly compressible at times. The turbulence observed just behind the termination shock by Voyager 2 was twice as strong. These observations place strong constraints on any model of "turbulence" in the heliosheath.

Burlaga, L. F.↗
DOE OSTI2025

Pressure suppresses the density wave order in kagome metal LuNb 6 ⁢Sn 6

The density waves that develop in kagome metals ScV 6 ⁢Sn 6 and LuNb 6 ⁢Sn 6 at low temperature appear to arise from underfilled atomic columns within a V-Sn or Nb-Sn scaffolding. Compressing this network with applied pressure in ScV 6 ⁢Sn 6 suppressed the structural transition temperature by constraining atomic rattling and inhibiting the shifts that define the structural modulation. We predicted that the density wave transition in LuNb 6⁢ Sn 6 at 68 K would be suppressed by pressure as well. In this Letter, we examine the pressure dependence of the density wave transition by measuring resistance vs temperature up to 2.26 GPa. We found the transition temperature is smoothly depressed and disappears around 1.9 GPa. In conclusion, this result not only addresses our prediction, but strengthens the rattling chains origin of structural instabilities in the HfFe 6⁢ Ge 6 -type kagome metals.

Charge density waves↗
DOE OSTI2019

Tensor-decomposition techniques for ab initio nuclear structure calculations: From chiral nuclear potentials to ground-state energies

The computational resources needed to generate the ab initio solution of the nuclear many-body problem for increasing mass number and/or accuracy necessitates innovative developments to improve upon (i) the storage of many-body operators and (ii) the scaling of many-body methods used to evaluate nuclear observables. The storing and efficient handling of many-body operators with high particle ranks is currently one of the major bottlenecks limiting the applicability range of ab initio studies with respect to mass number and accuracy. Recently, the application of tensor decomposition techniques to many-body tensors has proven highly beneficial to reduce the computational cost of ab initio calculations in quantum chemistry and solid-state physics. The impact of applying state-of-the-art tensor factorization techniques to modern nuclear Hamiltonians derived from chiral effective field theory is investigated. Subsequently, the error induced by the tensor decomposition of the input Hamiltonian on ground-state energies of closed-shell nuclei calculated via second-order many-body perturbation theory is benchmarked. The first proof-of-principles application of tensor-decomposition techniques to the nuclear Hamiltonian is performed. Two different tensor formats are investigated by systematically benchmarking the approximation error on matrix elements stored in various bases of interest. The analysis is achieved while including normal-ordered three-nucleon interactions that are nowadays used as input to the most advanced ab initio calculations in medium-mass nuclei. With the aid of the factorized Hamiltonian, the second-order perturbative correction to ground-state energies is decomposed and the scaling properties of the underlying tensor network are discussed. The employed tensor formats are found to lead to an efficient data compression of two-body matrix elements of the nuclear Hamiltonian. In particular, the sophisticated tensor hypercontraction scheme yields low tensor ranks with respect to both harmonic-oscillator and Hartree-Fock single-particle bases. It is found that the tensor rank depends on the two-body total angular momentum J for which one performs the decomposition, which is itself directly related to the sparsity of the corresponding tensor. Furthermore, including normal-ordered two-body contributions originating from three-body interactions does not compromise the efficient data compression. Ultimately, the use of factorized matrix elements authorizes controlled approximations of the exact second-order ground-state energy corrections. In particular, a small enough error is obtained from low-rank factorizations in He4,O16, and Ca40. It is presently demonstrated that tensor-decomposition techniques can be efficiently applied to systematically approximate the nuclear many-body Hamiltonian in terms of lower-rank tensors. Beyond the input Hamiltonian, tensor-decomposition techniques can be envisioned to scale down the cost of state-of-the-art nonperturbative many-body methods in order to extend ab initio studies to (i) higher precisions, (ii) larger masses, and (iii) nuclei of doubly open-shell character.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

How does ion temperature gradient turbulence depend on magnetic geometry? Insights from data and machine learning

Magnetic geometry has a significant effect on the level of turbulent transport in fusion plasmas. Here, we model and analyse this dependence using multiple machine learning methods and a dataset of >200 000 nonlinear gyrokinetic simulations of ion-temperature-gradient turbulence in diverse non-axisymmetric geometries. The dataset is generated using a large collection of both optimised and randomly generated stellarator equilibria. At fixed gradients and other input parameters, the turbulent heat flux varies between geometries by several orders of magnitude. Trends are apparent among the configurations with particularly high or particularly low heat flux. Regression and classification techniques from machine learning are then applied to extract patterns in the dataset. Due to a symmetry of the gyrokinetic equation, the heat flux and regressions thereof should be invariant to translations of the raw features in the parallel coordinate, similar to translation invariance in computer vision applications. Multiple regression models including convolutional neural networks (CNNs) and decision trees can achieve reasonable predictive power for the heat flux in held-out test configurations, with highest accuracy for the CNNs. Using Spearman correlation, sequential feature selection and Shapley values to measure feature importance, it is consistently found that the most important geometric lever on the heat flux is the flux surface compression in regions of bad curvature. The second most important geometric feature relates to the magnitude of geodesic curvature. These two features align remarkably with surrogates that have been proposed based on theory, while the methods here allow a natural extension to more features for increased accuracy. The dataset, released with this publication, may also be used to test other proposed surrogates, and we find that many previously published proxies do correlate well with both the heat flux and stability boundary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Effect of Strongly Magnetized Electrons and Ions on Heat Flow and Symmetry of Inertial Fusion Implosions

This Letter presents the first observation on how a strong, 500 kG, externally applied B field increases the mode-two asymmetry in shock-heated inertial fusion implosions. Using a direct-drive implosion with polar illumination and imposed field, we observed that magnetization produces a significant increase in the implosion oblateness (a 2.5 × larger P2 amplitude in x-ray self-emission images) compared with reference experiments with identical drive but with no field applied. Here, the implosions produce strongly magnetized electrons (ω e τ e >> 1) and ions (ω i τ i > 1) that, as shown using simulations, restrict the cross field heat flow necessary for lateral distribution of the laser and shock heating from the implosion pole to the waist, causing the enhanced mode-two shape.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2012

Magnetic Field Strength Fluctuations in the Heliosheath: Voyager 1 Observations During 2009

We analyze the microscale fluctuations of the magnetic field strength B on a scale of several hours observed by Voyager1 (V1) in the heliosheath during 2009. The microscale fluctuations of B range from coherent to stochastic structures. The amplitude of microscale fluctuations of B during 1 day is measured by the standard deviation (SD) of 48 s averages of B. The distribution of the daily values of SD is lognormal. SD(t) from day of year (DOY) 1 to 331, 2009, is very intermittent. SD(t) has a 1/f or "pink noise" spectrum on scales from 1 to 100 days, and it has a broad multifractal spectrum f(alpha) with 0.57 less than or equal to alpha less than or equal to 1.39. The time series of increments SD(t + tau) -- SD(t) has a pink noise spectrum with alpha(1) = 0.88 plus or minus 0.14 on scales from 1 to 100 days. The increments have a Tsallis (q-Gaussian) distribution on scales from 1 to 165 days, with an average q = 1.75 plus or minus 0.12. The skewness S and kurtosis K have Gaussian and lognormal distributions, respectively. The largest spikes in K(t) and S(t) are often associated with a change in B across a data gap and with identifiable physical structures. The "turbulence" observed by V1 during 2009 was weakly compressible on average but still very intermittent, highly variable, and highly compressible at times. The turbulence observed just behind the termination shock by Voyager 2 was twice as strong. These observations place strong constraints on any model of turbulence in the heliosheath.

Brulaga, L. F.↗
NASA NTRS1996

Trabecular Bone Mechanical Properties and Fractal Dimension

Countermeasures for reducing bone loss and muscle atrophy due to extended exposure to the microgravity environment of space are continuing to be developed and improved. An important component of this effort is finite element modeling of the lower extremity and spinal column. These models will permit analysis and evaluation specific to each individual and thereby provide more efficient and effective exercise protocols. Inflight countermeasures and post-flight rehabilitation can then be customized and targeted on a case-by-case basis. Recent Summer Faculty Fellowship participants have focused upon finite element mesh generation, muscle force estimation, and fractal calculations of trabecular bone microstructure. Methods have been developed for generating the three-dimensional geometry of the femur from serial section magnetic resonance images (MRI). The use of MRI as an imaging modality avoids excessive exposure to radiation associated with X-ray based methods. These images can also detect trabecular bone microstructure and architecture. The goal of the current research is to determine the degree to which the fractal dimension of trabecular architecture can be used to predict the mechanical properties of trabecular bone tissue. The elastic modulus and the ultimate strength (or strain) can then be estimated from non-invasive, non-radiating imaging and incorporated into the finite element models to more accurately represent the bone tissue of each individual of interest. Trabecular bone specimens from the proximal tibia are being studied in this first phase of the work. Detailed protocols and procedures have been developed for carrying test specimens through all of the steps of a multi-faceted test program. The test program begins with MRI and X-ray imaging of the whole bones before excising a smaller workpiece from the proximal tibia region. High resolution MRI scans are then made and the piece further cut into slabs (roughly 1 cm thick). The slabs are X-rayed again and also scanned using dual-energy X-ray absorptiometry (DEXA). Cube specimens are then cut from the slabs and tested mechanically in compression. Correlations between mechanical properties and fractal dimension will then be examined to assess and quantify the predictive capability of the fractal calculations.

Hogan, Harry A.↗
DOE OSTI2022

Data-driven model for divertor plasma detachment prediction

We present a fast and accurate data-driven surrogate model for divertor plasma detachment prediction leveraging the latent feature space concept in machine learning research. Our approach involves constructing and training two neural networks: an autoencoder that finds a proper latent space representation (LSR) of plasma state by compressing the multi-modal diagnostic measurements and a forward model using multi-layer perception (MLP) that projects a set of plasma control parameters to its corresponding LSR. By combining the forward model and the decoder network from autoencoder, this new data-driven surrogate model is able to predict a consistent set of diagnostic measurements based on a few plasma control parameters. In order to ensure that the crucial detachment physics is correctly captured, highly efficient 1D UEDGE model is used to generate training and validation data in this study. The benchmark between the data-driven surrogate model and UEDGE simulations shows that our surrogate model is capable of providing accurate detachment prediction (usually within a few per cent relative error margin) but with at least four orders of magnitude speed-up, indicating that performance-wise, it has the potential to facilitate integrated tokamak design and plasma control. Comparing with the widely used two-point model and/or two-point model formatting, the new data-driven model features additional detachment front prediction and can be easily extended to incorporate richer physics. This study demonstrates that the complicated divertor and scrape-off-layer plasma state has a low-dimensional representation in latent space. Understanding plasma dynamics in latent space and utilising this knowledge could open a new path for plasma control in magnetic fusion energy research.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

On hydromagnetic wave interactions in collisionless, high- β plasmas

We describe the interaction of parallel-propagating Alfvén waves with ion-acoustic waves and other Alfvén waves, in magnetized, high- β collisionless plasmas. This is accomplished through a combination of analytical theory and numerical fluid simulations of the Chew–Goldberger–Low (CGL) magnetohydrodynamic (MHD) equations closed by Landau-fluid heat fluxes. An asymptotic ordering is employed to simplify the CGL-MHD equations and derive solutions for the deformation of an Alfvén wave that results from its interaction with the pressure anisotropy generated either by an ion-acoustic wave or another, larger-amplitude Alfvén wave. The difference in time scales of acoustic and Alfvénic fluctuations at high- β means that interactions that are local in wavenumber space yield little modification to either mode within the time it takes the acoustic wave to Landau damp away. Instead, order-unity changes in the amplitude of Alfvénic fluctuations can result after interacting with frequency-matched acoustic waves. Additionally, we show that the propagation speed of an Alfvén-wave packet in an otherwise homogeneous background is a function of its self-generated pressure anisotropy. This allows for the eventual interaction of separate co-propagating Alfvén-wave packets of differing amplitudes. The results of the CGL-MHD simulations agree well with these predictions, suggesting that theoretical models relying on the interaction of these modes should be reconsidered in certain astrophysical environments. Applications of these results to weak Alfvénic turbulence and to the interaction between the compressive and Alfvénic cascades in strong, collisionless turbulence are also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2002

Photoevaporation of Clumps in Photodissociation Regions

We present the results of an investigation of the effects of Far Ultraviolet (FUV) radiation (6.0eV < hv < 13.6eV) from hot early type OB stars on clumps in star-forming molecular clouds. Clumps in FUV-illuminated regions (or photodissociation regions or PDRs) undergo external heating and photodissociation as they are exposed to the FUV field, resulting in a loss of cold, molecular lump mass as it is converted to warm atomic gas. The heating, if rapid, creates strong photoevaporative mass flows off the clump surfaces, and drives shocks into the clumps, compressing them to high densities. The clumps lose mass on relatively short timescales. The evolution of an individual clump is found to be sensitive to three dimensionless parameters: Nc0, the ratio of the initial column density of the clump to the column N(0) approx. 10(exp 21) cm(exp -2) of a warm FUV-heated surface region; upsilon, the ratio of the sound speed in the heated surface to that in the cold clump material: and t(FUV)t(c), the ratio of the "turn-on time" t(FUV) of the heating flux on a clump to its initial sound crossing-time t(c). The evolution also depends on whether a confining interclump medium exists, or whether the interclump region has negligible pressure, as is the case for turbulence-generated clumps. In this paper, we use spherical 1-D numerical hydrodynamic models as well as approximate analytical models to study the dependence of clump photoevaporation on the physical parameters of the clump, and to derive the dynamical evolution, mass loss rates and photoevaporative timescales of a clump for a variety of astrophysical situations. Turbulent clumps evolve so that their column densities are equal to a critical value determined by the local FUV field, and typically have short photo evaporation timescales, approx. 10(exp 4-5) years for a 1 M(solar mass) clump in a typical star-forming region (Nc0 = 10, upsilon = 10). Clumps with insufficient magnetic pressure support, and in strong FUV fields may be driven to collapse by the compressional effect of converging shock waves. We also estimate the rocket effect on photoevaporating clumps and find that it is significant only for the smallest clumps, with sizes much less than the extent of the PDR itself. Clumps that are confined by all interclump medium may either get completely photoevaporated, or may preserve a shielded core with a warm, dissociated, protective shell that, absorbs the incident FUV flux. We compare our results with observations of some well studied PDRs: the Orion Bar, M17SW NGC 2023 and the Rosette Nebula. The data are consistent with both interpretations of clump origin. turbulence and pressure confinement, with a slight indication for favouring the turbulent model for clumps over pressure-confined

Gorti, Uma↗
DOE OSTI2018

Role of Precipitates in Recrystallization Mechanisms of Nb-Mo Microalloyed Steel

Recrystallization behavior of a Nb-Mo microalloyed steel has been investigated using double-hit compression tests. A suitable reheating temperature and soaking time were established for the complete dissolution of the microalloying precipitates prior to hot deformation. The influence of hot deformation conditions on the flow behavior and the effect of alloying elements, Nb in particular, on austenite recrystallization kinetics are highlighted. The strain-induced precipitation has been found to play an important role in hindering recrystallization kinetics. Post-deformation microstructural analysis indicates that strain-induced grain boundary migration (SIMB) is one of the mechanisms for the formation of recrystallized grains. Retardation of recrystallization has been explained by estimating pinning force of Nb(C, N) precipitates and recrystallization driving force at 1000 °C, which is in good agreement with the experimental observations in the present study as well as other reported data. The study brings out a better understanding on the influence of strain-induced precipitates on recrystallization behavior.

36 MATERIALS SCIENCE↗
NASA NTRS2011

A THEMIS Survey of Flux Ropes and Traveling Compression Regions: Location of the Near-Earth Reconnection Site During Solar Minimum

A statistical study of flux ropes and traveling compression regions (TCRs) during the Time History of Events and Macroscale Interactions during Substorms (THEMIS) second tail season has been performed. A combined total of 135 flux ropes and TCRs in the range GSM X approx -14 to -31 R(sub E) were identified, many of these occurring in series of two or more events separated by a few tens of seconds. Those occurring within 10 min of each other were combined into aggregated reconnection events. For the purposes of this survey, these are most likely the products of reconnect ion occurring simultaneously at multiple, closely spaced x-lines as opposed to statistically independent episodes of reconnection. The 135 flux ropes and TCRs were grouped into 87 reconnection events; of these, 28 were moving tailward and 59 were moving Earthward. The average location of the near-Earth x-line determined from statistical analysis of these reconnection events is (X(sub GSM), Y*(sub GSM)) = (-30R(sub E), 5R(sub E)), where Y* includes a correction for the solar aberration angle. A strong east-west asymmetry is present in the tailward events, with >80% being observed at GSM Y* > O. Our results indicate that the Earthward flows are similarly asymmetric in the midtail region, becoming more symmetric inside - 18 R(sub E). Superposed epoch analyses indicate that the occurrence of reconnection closer to the Earth, i.e., X > -20 R(sub E), is associated with elevated solar wind velocity and enhanced negative interplanetary magnetic field B(sub z). Reconnection events taking place closer to the Earth are also far more effective in producing geomagnetic activity, judged by the AL index, than reconnection initiated beyond X approx -25 R(sub E).

Imber, S. M.↗
NASA NTRS2002

Efficient Low Dissipative High Order Schemes for Multiscale MHD Flows

Accurate numerical simulations of complex multiscale compressible viscous flows, especially high speed turbulence combustion and acoustics, demand high order schemes with adaptive numerical dissipation controls. Standard high resolution shock-capturing methods are too dissipative to capture the small scales and/or long-time wave propagations without extreme grid refinements and small time steps. An integrated approach for the control of numerical dissipation in high order schemes for the compressible Euler and Navier-Stokes equations has been developed and verified by the authors and collaborators. These schemes are suitable for the problems in question. Basically, the scheme consists of sixth-order or higher non-dissipative spatial difference operators as the base scheme. To control the amount of numerical dissipation, multiresolution wavelets are used as sensors to adaptively limit the amount and to aid the selection and/or blending of the appropriate types of numerical dissipation to be used. Magnetohydrodynamics (MHD) waves play a key role in drag reduction in highly maneuverable high speed combat aircraft, in space weather forecasting, and in the understanding of the dynamics of the evolution of our solar system and the main sequence stars. Although there exist a few well-studied second and third-order high-resolution shock-capturing schemes for the MHD in the literature, these schemes are too diffusive and not practical for turbulence/combustion MHD flows. On the other hand, extension of higher than third-order high-resolution schemes to the MHD system of equations is not straightforward. Unlike the hydrodynamic equations, the inviscid MHD system is non-strictly hyperbolic with non-convex fluxes. The wave structures and shock types are different from their hydrodynamic counterparts. Many of the non-traditional hydrodynamic shocks are not fully understood. Consequently, reliable and highly accurate numerical schemes for multiscale MHD equations pose a great challenge to algorithm development. In addition, controlling the numerical error of the divergence free condition of the magnetic fields for high order methods has been a stumbling block. Lower order methods are not practical for the astrophysical problems in question. We propose to extend our hydrodynamics schemes to the MHD equations with several desired properties over commonly used MHD schemes.

Sjoegreen, Bjoern↗
DOE OSTI2024

Overview of ASDEX upgrade results in view of ITER and DEMO

Experiments on ASDEX Upgrade (AUG) in 2021 and 2022 have addressed a number of critical issues for ITER and EU DEMO. A major objective of the AUG programme is to shed light on the underlying physics of confinement, stability, and plasma exhaust in order to allow reliable extrapolation of results obtained on present day machines to these reactor-grade devices. Concerning pedestal physics, the mitigation of edge localised modes (ELMs) using resonant magnetic perturbations (RMPs) was found to be consistent with a reduction of the linear peeling-ballooning stability threshold due to the helical deformation of the plasma. Conversely, ELM suppression by RMPs is ascribed to an increased pedestal transport that keeps the plasma away from this boundary. Candidates for this increased transport are locally enhanced turbulence and a locked magnetic island in the pedestal. The enhanced D-alpha (EDA) and quasi-continuous exhaust (QCE) regimes have been established as promising ELM-free scenarios. Here, the pressure gradient at the foot of the H-mode pedestal is reduced by a quasi-coherent mode, consistent with violation of the high-n ballooning mode stability limit there. This is suggestive that the EDA and QCE regimes have a common underlying physics origin. In the area of transport physics, full radius models for both L- and H-modes have been developed. These models predict energy confinement in AUG better than the commonly used global scaling laws, representing a large step towards the goal of predictive capability. A new momentum transport analysis framework has been developed that provides access to the intrinsic torque in the plasma core. In the field of exhaust, the X-Point Radiator (XPR), a cold and dense plasma region on closed flux surfaces close to the X-point, was described by an analytical model that provides an understanding of its formation as well as its stability, i.e., the conditions under which it transitions into a deleterious MARFE with the potential to result in a disruptive termination. With the XPR close to the divertor target, a new detached divertor concept, the compact radiative divertor, was developed. Here, the exhaust power is radiated before reaching the target, allowing close proximity of the X-point to the target. No limitations by the shallow field line angle due to the large flux expansion were observed, and sufficient compression of neutral density was demonstrated. With respect to the pumping of non-recycling impurities, the divertor enrichment was found to mainly depend on the ionisation energy of the impurity under consideration. In the area of MHD physics, analysis of the hot plasma core motion in sawtooth crashes showed good agreement with nonlinear 2-fluid simulations. This indicates that the fast reconnection observed in these events is adequately described including the pressure gradient and the electron inertia in the parallel Ohm’s law. Concerning disruption physics, a shattered pellet injection system was installed in collaboration with the ITER International Organisation. Thanks to the ability to vary the shard size distribution independently of the injection velocity, as well as its impurity admixture, it was possible to tailor the current quench rate, which is an important requirement for future large devices such as ITER. Progress was also made modelling the force reduction of VDEs induced by massive gas injection on AUG. The H-mode density limit was characterised in terms of safe operational space with a newly developed active feedback control method that allowed the stability boundary to be probed several times within a single discharge without inducing a disruptive termination. Regarding integrated operation scenarios, the role of density peaking in the confinement of the ITER baseline scenario (high plasma current) was clarified. The usual energy confinement scaling ITER98(p,y) does not capture this effect, but the more recent H20 scaling does, highlighting again the importance of developing adequate physics based models. Advanced tokamak scenarios, aiming at large non-inductive current fraction due to non-standard profiles of the safety factor in combination with high normalised plasma pressure were studied with a focus on their access conditions. A method to guide the approach of the targeted safety factor profiles was developed, and the conditions for achieving good confinement were clarified. Based on this, two types of advanced scenarios (‘hybrid’ and ‘elevated’

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗