Search NASA⌕ Search

SEARCH · Search NASA

Results for “3D equilibrium”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

3D Equilibrium Reconstruction and MHD Studies in Stellarators (Final Report)

This is the final technical report for the grant Three Dimensional Equilibrium Reconstruction and MHD Studies in Stellarators. The major goals of this project were to: further develop V3FIT code capabilities, continue to support experiments using the V3FIT code, continue theoretical support for the CTH experiment at Auburn, particularly in the area of NIMROD simulations of CTH discharges. Dr. Lopez Ortiz has made significant progress on his analysis of NIMROD simulations of CTH MHD phenomena. In particular, his study has focused on the fixed point dynamics and energy flow amongst toroidal Fourier modes in sawtooth oscillations simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Inputs and Outputs for paper "The DESC stellarator code suite. Part 1. Quick and accurate equilibria computations"

Contains the DESC and VMEC input and output files used in the paper, as well as the plotting scripts used to create the figures in the paper (Current with the arXiv Mar 31 2022 version, https://arxiv.org/abs/2203.17173v1): 3D equilibrium codes are vital for stellarator design and operation, and high-accuracy equilibria are also necessary for stability studies. This paper details comparisons of two 3D equilibrium codes, VMEC, which uses a steepest-descent algorithm to reach a minimum-energy plasma state, and DESC, which minimizes the MHD force error in real space directly. Accuracy as measured by final plasma energy and satisfaction of MHD force balance, as well as other metrics, will be presented for each code, along with the computation time. It is shown that DESC is able to achieve more accurate solutions, especially near-axis. DESC's global Fourier-Zernike basis also yields the solution everywhere in the plasma volume, not just on discrete flux surfaces. Further, DESC can compute the same accuracy solution as VMEC in an order of magnitude less time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Surrogate models for plasma displacement and current in 3D perturbed magnetohydrodynamic equilibria in tokamaks

Abstract A numerical database of over one thousand perturbed three-dimensional (3D) equilibria has been generated, constructed based on the MARS-F (Liu et al 2000 Phys. Plasmas 7 3681) computed plasma response to the externally applied 3D field sources in multiple tokamak devices. Perturbed 3D equilibria with the n = 1–4 ( n is the toroidal mode number) toroidal periodicity are computed. Surrogate models are created for the computed perturbed 3D equilibrium utilizing model order reduction (MOR) techniques. In particular, retaining the first few eigenstates from the singular value decomposition (SVD) of the data is found to produce reasonably accurate MOR-representations for the key perturbed quantities, such as the perturbed parallel plasma current density and the plasma radial displacement. SVD also helps to reveal the core versus edge plasma response to the applied 3D field. For the database covering the conventional aspect ratio devices, about 95% of data can be represented by the truncated SVD-series with inclusion of only the first five eigenstates, achieving a relative error (RE) below 20%. The MOR-data is further utilized to train neural networks (NNs) to enable fast reconstruction of perturbed 3D equilibria, based on the two-dimensional equilibrium input and the 3D source field. The best NN-training is achieved for the MOR-data obtained with a global SVD approach, where the full set of samples used for NN training and testing are stretched and form a large matrix which is then subject to SVD. The fully connected multi-layer perceptron, with one or two hidden layers, can be trained to predict the MOR-data with less than 10% RE. As a key insight, a better strategy is to train separate NNs for the plasma response fields with different toroidal mode numbers. It is also better to apply MOR and to subsequently train NNs separately for conventional and low aspect ratio devices, due to enhanced toroidal coupling of Fourier spectra in the plasma response in the latter case.

3D equilibrium↗
DOE OSTI2022

Light-induced dimension crossover dictated by excitonic correlations

In low-dimensional systems with strong electronic correlations, the application of an ultrashort laser pulse often yields novel phases that are otherwise inaccessible. The central challenge in understanding such phenomena is to determine how dimensionality and many-body correlations together govern the pathway of a non-adiabatic transition. To this end, we examine a layered compound, 1T-TiSe 2 , whose three-dimensional charge-density-wave (3D CDW) state also features exciton condensation due to strong electron-hole interactions. We find that photoexcitation suppresses the equilibrium 3D CDW while creating a nonequilibrium 2D CDW. Remarkably, the dimension reduction does not occur unless bound electron-hole pairs are broken. This relation suggests that excitonic correlations maintain the out-of-plane CDW coherence, settling a long-standing debate over their role in the CDW transition. Our findings demonstrate how optical manipulation of electronic interaction enables one to control the dimensionality of a broken-symmetry order, paving the way for realizing other emergent states in strongly correlated systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Physics Basis, Optimization, and Control for Integrated 3D Edge Long-pulse Tokamak Scenarios

The aim of this project with international collaborations (PI: J.-K. Park of Princeton Plasma Physics Laboratory) is to develop a unified physics basis and predictive capability for the control of edge-localized modes (ELMs) with optimized non-axisymmetric (3D) fields, by leveraging the unique research capabilities of international KSTAR tokamak in Korea. At the University of California, Irvine, we use a fusion simulation code GTC to study the effects of the 3D magnetic fields (including both magnetic islands and stochastic magnetic field lines) on microturbulence. We focus on the effects of 3D equilibrium on turbulent and neoclassical transport that leads to ELM suppression. During the project period, the UCI team of this collaboration has made significant progress in understanding turbulent transport in 3D geometry including both tokamaks with RMP and stellarators. Four UCI postdoctoral researchers, Drs. Javier Nicolau, Gyungjin Choi, Pengfei Liu, and Xishuo Wei were partially supported by this project and collaborating with the entire project team. We have published 4 journal papers and gave several invited talks. This project advance physics understand of 3D effects on microturbulence, specially effects of magnetic islands. It will contribute to the DOE SciDAC ISEP project in the area of GTC development for general 3D geometry capability, rigorous verification and validation, and readiness for the emerging exascale computing. Through collaborations with computational scientists, the project helps to advance high performance computing in the area of GPU acceleration, linear solver, and data management.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRSDate not supplied

The Thickness of Electric Current Sheets

The thickness of electric current sheets is extremely important, especially as it relates to the onset of fast magnetic reconnection. Onset determines how much magnetic free energy can build up in a field before it is explosively released. This has implications for coronal heating, CMEs, and many other phenomena on the Sun and throughout the universe. Significant effort has been devoted to the question of whether equilibrium current sheets in realistic geometries have finite or zero thickness. Using a simple force balance analysis, we show why current sheets without a guide field (2D) and with a guide field that is invariant in the guide field direction (2.5D) cannot be in equilibrium if they have both finite thickness and finite width. We then estimate the conditions under which the tension of a curved line-tied guide field can produce a 3D equilibrium sheet that is finite in all dimensions. Finally, we speculate that some quasi-statically evolving current sheets undergoing slow stressing may reach a critical point at which equilibrium is no longer possible. They then rapidly thin and reconnect.

James A Klimchuk↗
DOE OSTI2021

The magnetic structure of long-wavelength magnetohydrodynamic modes in current-carrying stellarator plasmas

Magnetic field fluctuations are observed in current-carrying stellarator plasmas when the rotational transform is close to a rational value at the edge of the plasma. At low plasma pressure, these fluctuations are associated with perturbed currents parallel to the equilibrium magnetic field lines. A model for these magnetohydrodynamic modes in a low-β, three-dimensional stellarator equilibria has been developed. A set of helical current filaments are constrained to mimic the structure of magnetic field lines on rational surfaces derived from three-dimensional (3D) equilibrium reconstructions. Transformation to straight field line coordinates then allows fitting of the poloidal magnetic sensor data to a single harmonic function, which fixes the modeled toroidal mode structure via the field line flow geometry. The developed procedure accurately captures phase and amplitude variation for m/n = 3/2, 3/1, and 4/1 modes in the 3D equilibria of the compact toroidal hybrid experiment.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Pines’ demon observed as a 3D acoustic plasmon in Sr2RuO4

Abstract The characteristic excitation of a metal is its plasmon, which is a quantized collective oscillation of its electron density. In 1956, David Pines predicted that a distinct type of plasmon, dubbed a ‘demon’, could exist in three-dimensional (3D) metals containing more than one species of charge carrier 1 . Consisting of out-of-phase movement of electrons in different bands, demons are acoustic, electrically neutral and do not couple to light, so have never been detected in an equilibrium, 3D metal. Nevertheless, demons are believed to be critical for diverse phenomena including phase transitions in mixed-valence semimetals 2 , optical properties of metal nanoparticles 3 , soundarons in Weyl semimetals 4 and high-temperature superconductivity in, for example, metal hydrides 3,5–7 . Here, we present evidence for a demon in Sr 2 RuO 4 from momentum-resolved electron energy-loss spectroscopy. Formed of electrons in the β and γ bands, the demon is gapless with critical momentum q c = 0.08 reciprocal lattice units and room-temperature velocity v = (1.065 ± 0.12) × 10 5 m s −1 that undergoes a 31% renormalization upon cooling to 30 K because of coupling to the particle–hole continuum. The momentum dependence of the intensity of the demon confirms its neutral character. Our study confirms a 67-year old prediction and indicates that demons may be a pervasive feature of multiband metals.

Science & Technology - Other Topics↗
DOE OSTI2020

Global gyrokinetic particle simulations of microturbulence in W7-X and LHD stellarators

Global gyrokinetic particle simulations of electrostatic ion temperature gradient (ITG) instability show that the most unstable eigenmode is localized to some magnetic fieldlines or discrete locations on the poloidal plane in the Wendelstein 7-X (W7-X) stellarator due to its mirror-like magnetic fields, which vary strongly in the toroidal direction and induce coupling of more toroidal harmonics (n) to form the linear eigenmode than in the Large Helical Device (LHD) stellarator. Nonlinear electrostatic simulation results show that self-generated zonal flows are the dominant saturation mechanism for the ITG instabilities in both the LHD and W7-X. Additionally, radial widths of the fluctuation intensity in both the LHD and W7-X are significantly broadened from the linear phase to the nonlinear phase due to turbulence spreading. Finally, nonlinear spectra in the W7-X are dominated by low-n harmonics, which can be generated both by nonlinear toroidal coupling of high-n harmonics and by linear toroidal coupling with large amplitude zonal flows due to the 3D equilibrium magnetic fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Pressure-driven stochastization of flux surfaces in stellarator equilibria: a review of the experimental observations and their analysis

Stellarator vacuum fields are designed to have a dense set of nested flux surfaces, with only small islands and stochastic regions in the interior of the confinement region. As the plasma pressure is increased, pressure driven currents appear, and they can significantly modify the vacuum field. The question arises whether pressure driven currents can cause significant flux surface breakage. Both the Wendelstein 7-AS and large helical device stellarators had dedicated experimental campaigns to study the β limit in those devices. In this paper, we review the evidence that the flux surfaces in a significant fraction of the plasma volume were stochastized by the pressure-driven currents in the highest β experiments. That phenomenon appears to have been the dominant source of pressure-driven flux surface loss in these experiments.

stellarator equilibrium↗
DOE OSTI2023

Modeling of plasma beta effects on the island divertor transport in the standard configuration of W7-X

The influence of plasma beta effects on the edge plasma transport in the Wendelstein 7-X standard configuration is studied systematically by using EMC3-EIRENE combined with a 3D equilibrium code named HINT. The magnetic topology changes induced by plasma beta effects are significantly reflected in plasma transport behaviors and heat flux patterns on divertor targets. After validating the modeling strategy by comparisons with experimental data, the extended simulations for high performance plasmas show that the threshold separatrix density for accessing the power detachment is reduced in higher beta plasmas. Compared with the vacuum field case, the impurity radiation distributions with finite beta effects are modified in the magnetic island region. The divertor heat flux is distributed more evenly along the toroidal direction on the strike line at the vertical target. The strike line on the horizontal target moves towards the pumping gap with an increase in the plasma beta. In addition, the different pressure profiles with the same central beta also result in a modified heat flux pattern on the divertor targets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

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

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

Physics↗
DOE OSTI2020

STELLOPT

The STELLOPT code is designed to optimize 3D MHD equilibria to a set of target physics parameters encompassing stellarator design and 3D equilibrium reconstructions.

Lazerson, Samuel↗
NASA NTRS1992

Bending response of thin-walled laminated composite cylinders

Closed-form linear response solutions are presented for the bending response of thin-walled laminated composite cylinders. Displacements and intralaminar stresses are computed using the Donnell theory, and interlaminar stresses are computed using 3D equilibrium equations of elasticity. Laminates with stacking sequencies and an aluminum cylinder are used to illustrate the displacement response, intralaminar, and interlaminar stresses for two length-to-radius ratios of 1 and 5. Data obtained indicate that if only fiber breakage causes ultimate failure, material failure will occur on the compression side of the cylinders, at the ends, due to fiber compression. The interlaminar normal stress is an order of magnitude smaller than the interlaminar shear stresses and it is not likely to initiate failure.

Fuchs, H. P.↗
DOE OSTI2020

Edge ICRF simulations in 3D geometry: From MHD equilibrium to coupling determination

We present in this work a consistent numerical scheme that allows the computation of 3D magnetic fields a nd 3D density profiles and their usage in ion cyclotron range of frequencies (ICRF) coupling simulations. We first utilize the PARVMEC code to compute the 3D free-boundary plasma equilibrium in the ideal magnetohydrodynamic (MHD) approximation. Since the PARVMEC solution is only defined within the last closed flux surface (LCFS), the magnetic field domain is extended to the scrape-off layer (SOL) via the BMW code, which computes a divergence-free magnetic field solution arising from the external conductors’ vacuum field a nd t he PARVMEC flux surface currents. This magnetic reconstruction is then used in the EMC3-EIRENE transport code in order to compute 3D density profiles. In the last step, the RAPLICASOL code is utilized to compute the ICRF antenna S-matrices resulting from the 3D density profiles. We exemplify this scheme for the ASDEX Upgrade tokamak. A new implementation of a curved model for the ASDEX Upgrade ICRF 2-strap antenna in RAPLICASOL allows simulations in realistic geometry, without any coordinate transformations.

Lopez, Graham↗
DOE OSTI2020

Equilibrium shapes of facetted 3D metal nanoclusters intercalated near the surface of layered materials

Experimental studies indicate that 3D crystalline metal nanoclusters (NCs) intercalated under the surface of graphite have flat-topped equilibrated shapes. In this work, we characterize the shapes of these facetted NCs sandwiched between a blanketing graphene layer and the underlying graphite substrate. Specifically, we focus on the cases of fcc Cu and hcp Fe NCs. Furthermore, the analysis involves numerical minimization of the system energy for a specified NC volume and NC height, the latter corresponding to the separation between parallel top and bottom facets. Our numerical analysis quantifies how the distance of the side facet planes from center of the nanocluster varies linearly with a natural characteristic linear dimension of the nanocluster. Calculated shapes of fcc Cu and hcp Fe NCs are consistent with the hexagonal footprints observed in scanning tunneling microscopy studies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

FES 2023 Theory & Simulation Performance Target

The Report for the 2023 Theory & Simulation Performance Target (TSPT) is presented herein. The 1st, 2nd and 3rd Quarter Reports are given in Sec. I, Sec. II and Sec. III. These are unchanged from those previously submitted, except that in this final document the equations, figures and citations have been automatically renumbered. The 4th Quarter Report is given in Sec. IV. The following material presents verification calculations of several widely used three-dimensional (3D) magnetohydrodynamic (MHD) equilibrium codes in strongly shaped configurations; for both fixed and free-boundary equilibria; for the so-called vacuum case and for the case of high plasma pressure; and, for magnetic fields that have nested magnetic surfaces and for fields that have a magnetic island. The configurations studied includes the recently discovered “precise-QA” case [Phys. Rev. Lett., 128, 035001 (2002)] and other highly optimized configurations. Because of the strong (extreme) shaping, the latter of these pose challenges to the 3D MHD codes. A new method is introduced for constructing coordinates in arbitrarily shaped toroidal domains, and this can improve the robustness of 3D equilibrium calculations. Beyond the equilibrium calculation, this report demonstrates the capability of performing non-linear extended-MHD modelling of optimized stellarator equilibria. This capability allows the assessment not only of linear stability thresholds but also of the nonlinear evolution of unstable equilibria to new configurations that might be far from the original configuration, with significantly different transport properties and magnetic topology. They also enable the self-consistent determination of pressure profiles in non-integrable magnetic topology through the use of highly anisotropic thermal-transport equations. The calculations presented here represent the first demonstration of extended-MHD calculations in strongly-shaped, optimized quasi-symmetric configurations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Application of machine learning and artificial intelligence to extend EFIT equilibrium reconstruction

Recent progress in the application of machine learning (ML)/artificial intelligence (AI) algorithms to improve the Equilibrium Fitting (EFIT) code equilibrium reconstruction for fusion data analysis applications is presented. A device-independent portable core equilibrium solver capable of computing or reconstructing equilibrium for different tokamaks has been created to facilitate adaptation of ML/AI algorithms. A large EFIT database comprising of DIII-D magnetic, motional Stark effect, and kinetic reconstruction data has been generated for developments of EFIT model-order-reduction (MOR) surrogate models to reconstruct approximate equilibrium solutions. Furthermore, a neural-network MOR surrogate model has been successfully trained and tested using the magnetically reconstructed datasets with encouraging results. Other progress includes developments of a Gaussian process Bayesian framework that can adapt its many hyperparameters to improve processing of experimental input data and a 3D perturbed equilibrium database from toroidal full magnetohydrodynamic linear response modeling using the Magnetohydrodynamic Resistive Spectrum - Feedback (MARS-F) code for developments of 3D-MOR surrogate models.

Gaussian process↗