Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MHD 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.

At least 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Fast Particle-Wave Interactions and Alfvén Eigenmodes in JET Tokamak Plasmas

This document serves as the closeout report of DOE Grant Award No. DE-FG02-99ER54563 with project period 1 April 2015 through 31 March 2020. The project comprised the international collaboration between MIT and EU scientists on the JET facility to improve our understanding of the physics of energetic particle-wave interactions by measuring the damping rates of stable Alfvén Eigenmodes (AEs) and unstable energetic particle driven modes. More specifically, this project involved the continued participation of MIT, the Culham Centre for Fusion Energy, the Swiss Plasma Center, and theorists from various European laboratories and from UC Irvine. Past contributions from the University of São Paulo are gratefully acknowledged. The Alfvén Eigenmode Active Diagnostic (AEAD) on JET discharges was successfully upgraded in this period and has extracted physically useful information, in primarily deuterium plasma discharges. It is expected that such experiments would be continued during the DT campaign in CY 2021 to assess the damping rates of similar modes in the presence of alpha particles, a product of fusing burning plasma. Only JET would carry out such experiments in the near term in the world. It is important to note that this collaboration was continued under DOE Grant Award No. DE-SC0014264 from 1 April 2020 through the present. In this grant period, the AEAD was upgraded with a set of individual amplifiers for each of a set of six antennas in two toroidally opposite locations. These new amplifiers allowed targeted selection of the antennas’ toroidal spectrum for AEs of interest with toroidal mode number | n | ≤ 20. The resonant detection and measurement of the damping rates of AEs was obtained from magnetic probes that could be compared with theory and simulations. This is a key topic of investigation for ITER and all other next-step fusion experiments, where such modes will interact with energetic particles produced by fusion reactions (alpha particles), Neutral Beam Injection (NBI) and Ion Cyclotron Resonance Heating (ICRH). The majority of the observations during this period were of Toroidal Alfvén Eigenmodes (TAEs) as these are most commonly observed on JET; a number of measurements were made during dedicated TAE experiments. However, a new set of lower frequency band filters were procured with the goal of studying Geodesic Acoustic Modes (GAMs), Beta (Acoustic) AEs (BAE/BAAEs), and Reverse Shear AEs (RSAEs), also predicted by theory. Commissioning and optimizations were successfully completed in this period, and the diagnostic has been in operation during the more recent JET campaigns. Initial results have been obtained in dedicated TAE experiments and successfully compared to drift-kinetic theory. A wide range of theoretical studies have been undertaken in support of the upcoming JET campaigns experiments. To supplement the ongoing use of ideal MHD codes, such as MISHKA during the studies of TAEs, gyrokinetic simulations of low frequency AEs were performed in collaboration with UC Irvine. The Gyrokinetic Toroidal Code (GTC) was used to determine the structure, frequencies and stability of AEs in JET plasmas. Thus, a solid scientific foundation has been laid for future DT campaigns in JET in the 2021 operation period when the damping rates of relevant modes in the presence of alpha particles could be assessed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1986

Numerical simulation of flare energy build-up and release via Joule dissipation

A new numerical MHD model is developed to study the evolution of an active region due to photospheric converging motion, which leads to magnetic-energy buildup in the form of electric current. Because this new MHD model has incorporated finite conductivity, the energy conversion occurs from magnetic mode to thermal mode through Joule dissipation. In order to test the causality relationship between the occurrence of flare and photospheric motion, a multiple-pole configuration with neutral point is used. Using these results it is found that in addition to the converging motion, the initial magnetic-field configuration and the redistribution of the magnetic flux at photospheric level enhance the possibility for the development of a flare.

Wu, S. T.↗
DOE OSTI2022

Kinetic Simulations of Instabilities and Particle Acceleration in Cylindrical Magnetized Relativistic Jets

Relativistic magnetized jets, such as those from AGN, GRBs, and XRBs, are susceptible to current- and pressure-driven MHD instabilities that can lead to particle acceleration and nonthermal radiation. Here, we investigate the development of these instabilities through 3D kinetic simulations of cylindrically symmetric equilibria involving toroidal magnetic fields with electron–positron pair plasma. Generalizing recent treatments by Alves et al. and Davelaar et al., we consider a range of initial structures in which the force due to toroidal magnetic field is balanced by a combination of forces due to axial magnetic field and gas pressure. We argue that the particle energy limit identified by Alves et al. is due to the finite duration of the fast magnetic dissipation phase. We find a rather minor role of electric fields parallel to the local magnetic fields in particle acceleration. In all investigated cases, a kink mode arises in the central core region with a growth timescale consistent with the predictions of linearized MHD models. In the case of a gas-pressure-balanced (Z-pinch) profile, we identify a weak local pinch mode well outside the jet core. We argue that pressure-driven modes are important for relativistic jets, in regions where sufficient gas pressure is produced by other dissipation mechanisms.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Yield and PDV Measurements of the Godiva-IV Critical Assembly (Preliminary Report)

A high dynamic range, scinillator based detector, the MHD-240, was used to measure the gamma and neutron flux of several prompt critical bursts of the Godiva-IV assembly. The measurements provided the alpha and FWHM for several different burst temperatures. Calibrations of the MHD240 gamma and neutron sensitivity along with MCNP based simulation of the Godiva-IV assembly will allow us to also determine the total fission yield of the burst integrated over approximately 5 ms. Work is ongoing to check the calibrations of the MHD-240, and improve the MCNP simulations. This is a preliminary report, and the fission yields reported here are not final and are likely inaccurate. The measured alpha and FWHM are not impacted by the calibration issues and are considered final.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2008

Magnetotail Current Sheet Thinning and Magnetic Reconnection Dynamics in Global Modeling of Substorms

Magnetotail current sheet thinning and magnetic reconnection are key elements of magnetospheric substorms. We utilized the global MHD model BATS-R-US with Adaptive Mesh Refinement developed at the University of Michigan to investigate the formation and dynamic evolution of the magnetotail thin current sheet. The BATSRUS adaptive grid structure allows resolving magnetotail regions with increased current density up to ion kinetic scales. We investigated dynamics of magnetotail current sheet thinning in response to southwards IMF turning. Gradual slow current sheet thinning during the early growth phase become exponentially fast during the last few minutes prior to nightside reconnection onset. The later stage of current sheet thinning is accompanied by earthward flows and rapid suppression of normal magnetic field component $B-z$. Current sheet thinning set the stage for near-earth magnetic reconnection. In collisionless magnetospheric plasma, the primary mechanism controlling the dissipation in the vicinity of the reconnection site is non-gyrotropic effects with spatial scales comparable with the particle Larmor radius. One of the major challenges in global MHD modeling of the magnetotail magnetic reconnection is to reproduce fast reconnection rates typically observed in smallscale kinetic simulations. Bursts of fast reconnection cause fast magnetic field reconfiguration typical for magnetospheric substorms. To incorporate nongyritropic effects in diffusion regions we developed an algorithm to search for magnetotail reconnection sites, specifically where the magnetic field components perpendicular to the local current direction approaches zero and form an X-type configuration. Spatial scales of the diffusion region and magnitude of the reconnection electric field are calculated self-consistently using MHD plasma and field parameters in the vicinity of the reconnection site. The location of the reconnection sites and spatial scales of the diffusion region are updated during the simulations. Such an approach allows quantifying the interaction between large-scale global magnetospheric dynamics and microphysical processes in diffusion regions localized near reconnection sites. To clarify the role of smallscale non-MHD effects in diffusion region on the global magnetospheric dynamic and to test different models of dissipation we perform simulations with steady southward IMF driving of the magnetosphere.

Kuznetsova, M. M.↗
DOE OSTI2024

Interpretative 3D MHD modelling of deuterium SPI into a JET H-mode plasma

Abstract The pre-thermal quench (pre-TQ) dynamics of a pure deuterium ( D 2 ) shattered pellet injection (SPI) into a 3 MA / 7 MJ JET H-mode plasma is studied via 3D non-linear MHD modelling with the JOREK code. The interpretative modelling captures the overall evolution of the measured density and radiated power. The simulations also identify the importance of the drifts of ablation plasmoids towards the tokamak low field side (LFS) and the impurities in the background plasma in fragment penetration, assimilation, radiative cooling and MHD activity in D 2 SPI experiments. It is found that plasmoid drifts lead to an about 70% reduction of the central line-integrated density (compared to a simulation without drifts) in the JET D 2 SPI discharge considered. Impurities that pre-exist before the SPI as well as those from possible impurity influxes related to the SPI are shown to dominate the radiation in the considered discharge. With inputs from JOREK simulations, modelling with the Lagrangian particle-based pellet code PELOTON reproduces the deviation of the SPI fragments in the direction of the major radius as observed by the fast camera. This confirms the role of rocket effects and plasmoid drifts in the considered discharge and reinforces the validity of the JOREK modelling. The limited core density rise due to plasmoid drifts and the strong radiative cooling and MHD activity with impurities (depending on their species and concentration) could limit the effectiveness of LFS D 2 SPI in runaway electron avoidance and are worth considering in the design of the ITER disruption mitigation system.

Physics↗
DOE OSTI2023

Optical design and synthetic analysis of the electron cyclotron emission imaging diagnostic of HL-2M tokamak

The electron cyclotron emission imaging (ECEI) diagnostic is a powerful tool to study the MHD and turbulent transport in magnetically confined fusion plasmas. In this work, the optical system including the local oscillator (LO) coupling and radio frequency (RF) receiving optics, has been designed and analyzed for the HL-2M ECEI diagnostic. The LO optics can illuminate the antenna array and drive the mixer diode to work efficiently, with more than 36% of the beam intensity for the channels which probe the plasma edge relative to those of channels which probe the plasma core. The RF optics aims at guiding the plasma emission signal to the antenna array. To meet different physical requirements, three types of field of view have been achieved for the plasma imaging, with zoom factors of around 1, 1.5 and 2, respectively. The focal surfaces are almost flat, with a maximum off-mapping (defined as the radial distance between the beam waists of the lower-/uppermost antennas and electron cyclotron emission layer) less than 2.5 cm, which can match the electron cyclotron emission layer quite well and suppress the image distortion in the plasma edge region. Furthermore, the impact of possible vibrations and installation error on the focal plane has been evaluated. In addition, the predefined MHD and turbulence perturbations are well reproduced by combining the optical simulation results and synthetic ECEI modeling, which further verifies the good performance of the RF optics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Improved Boundary Conditions for Coupled Geospace Models: An Application in Modeling Spacecraft Surface Charging Environment

Abstract Spacecraft surface charging in the inner magnetosphere often occurs in the pre‐midnight to the dawn sector when electron fluxes of tens of keV increase. Inner magnetosphere ring current models can be used to simulate/predict the spacecraft surface charging environment, with their outer boundary conditions specified either based on observations or provided by other models, such as MHD models. In the latter approach, using MHD quantities, the flux spectrum at the outer boundary is commonly assumed to follow a Kappa or Maxwellian distribution function. Such a method however often departs greatly from the realistic spectrum at E < tens of keV, a crucial energy range in the surface charging anomaly. In order to achieve a better representation of the surface charging environment, we propose to combine the MHD‐parameterized flux spectrum with an empirical electron flux model of E < 40 keV to set the electron flux boundary condition. Results indicate that as opposed to the case where the MHD‐parameterized flux distribution is solely used at the model boundary, simulations with the new boundary condition yields a more intense surface charging environment. The integrated electron flux between 10 < E < 50 keV, a measure of the severity of the surface charging environment, is significantly enhanced by 1‐2 orders of magnitude, leading to a much better agreement with Van Allen Probes measurements. This study hence demonstrates a reasonable solution to the setting of outer boundary conditions for inner magnetosphere models and is recommended for coupled geospace circulation models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Proton Acceleration in Low- β Magnetic Reconnection with Energetic Particle Feedback

Magnetic reconnection regions in space and astrophysics are known as active particle acceleration sites. There is ample evidence showing that energetic particles can take a substantial amount of converted energy during magnetic reconnection. However, there has been a lack of studies understanding the backreaction of energetic particles at magnetohydrodynamical scales in magnetic reconnection. To address this, we have developed a new computational method to explore the feedback by nonthermal energetic particles. This approach considers the backreaction from these energetic particles by incorporating their pressure into magnetohydrodynamics (MHD) equations. The pressure of the energetic particles is evaluated from their distribution evolved through Parker's transport equation, solved using stochastic differential equations (SDEs), so we coin the name MHD-SDE. Applying this method to low-β magnetic reconnection simulations, we find that reconnection is capable of accelerating a large fraction of energetic particles that contain a substantial amount of energy. When the feedback from these particles is included, their pressure suppresses the compression structures generated by magnetic reconnection, thereby mediating particle energization. Consequently, the feedback from energetic particles results in a steeper power-law energy spectrum. These findings suggest that feedback from nonthermal energetic particles plays a crucial role in magnetic reconnection and particle acceleration.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1986

Numerical simulations of comets - Predictions for Comet Giacobini-Zinner

Simulations of Comet Giacobini-Zinner's interaction with solar wind are described and results are presented. The simulations are carried out via the numerical solution of the ideal MHD equations as an initial value problem in a uniform solar wind. The calculations are performed on a Cartesian mesh centered at the comet. Results reveal that the first significant modifications of the solar wind along the ISEE/ICE trajectory will occur 100,000 km from the solar wind comet axis.

Fedder, J. A.↗
NASA NTRS1986

Explosive coalescence of magnetic islands

Simulation results from both the EM collisionless particle code and the MHD particle code reveal an explosive reconnection process associated with nonlinear evolution of the coalescence instability. The explosive coalescence is a self-similar process of magnetic collapse, and ensuing amplitude oscillations in the magnetic and electrostatic energies and temperatures are modeled by an equation of motion for the scale factor in the Sagdeev potential. This phenomenon may explain the rapid energy release of a certain class of solar flares during their impulsive phase.

Tajima, T.↗
NASA NTRS2001

Global MHD Models of the Solar Corona

Global magnetohydrodynamic (MHD) models of the solar corona are computationally intensive, numerically complex simulations that have produced important new results over the past few years. After a brief overview of how these models usually work, I will address three topics: (1) How these models are now routinely used to predict the morphology of the corona and analyze Earth and space-based remote observations of the Sun; (2) The direct application of these models to the analysis of physical processes in the corona and chromosphere and to the interpretation of in situ solar wind observations; and (3) The use of results from global models to validate the approximations used to make detailed studies of physical processes in the corona that are not otherwise possible using the global models themselves.

Suess, S. T.↗
NASA NTRS2001

A Global Magnetohydrodynamic Model of Jovian Magnetosphere

The goal of this project was to develop a new global magnetohydrodynamic model of the interaction of the Jovian magnetosphere with the solar wind. Observations from 28 orbits of Jupiter by Galileo along with those from previous spacecraft at Jupiter, Pioneer 10 and 11, Voyager I and 2 and Ulysses, have revealed that the Jovian magnetosphere is a vast, complicated system. The Jovian aurora also has been monitored for several years. Like auroral observations at Earth, these measurements provide us with a global picture of magnetospheric dynamics. Despite this wide range of observations, we have limited quantitative understanding of the Jovian magnetosphere and how it interacts with the solar wind. For the past several years we have been working toward a quantitative understanding of the Jovian magnetosphere and its interaction with the solar wind by employing global magnetohydrodynamic simulations to model the magnetosphere. Our model has been an explicit MHD code (previously used to model the Earth's magnetosphere) to study Jupiter's magnetosphere. We continue to obtain important insights with this code, but it suffers from some severe limitations. In particular with this code we are limited to considering the region outside of 15RJ, with cell sizes of about 1.5R(sub J). The problem arises because of the presence of widely separated time scales throughout the magnetosphere. The numerical stability criterion for explicit MHD codes is the CFL limit and is given by C(sub max)(Delta)t/(Delta)x less than 1 where C(sub max) is the maximum group velocity in a given cell, (Delta)x is the grid spacing and (Delta)t is the time step. If the maximum wave velocity is C(sub w) and the flow speed is C(sub f), C(sub max) = C(sub w) + C(sub f). Near Jupiter the Alfven wave speed becomes very large (it approaches the speed of light at one Jovian radius). Operating with this time step makes the calculation essentially intractable. Therefore under this funding we have been designing a new MHD model that will be able to compute solutions in the wide parameter regime of the Jovian magnetosphere.

Walker, Raymond J.↗
NASA NTRS2003

A Plasmoid Thruster for Space Propulsion

There are a number of possible advantages to using accelerated plasmoids for in-space propulsion. A plasmoid is a compact plasma structure with an integral magnetic field. They have been studied extensively in controlled fusion research and are classified according to the relative strength of the poloidal and toroidal magnetic field (B(sub p), and B(sub t), respectively). An object with B(sub p), / B(sub t) much greater than 1 is classified as a Field Reversed Configuration (FRC); if B(sub p) approximately equal to B(sub t), it is called a Spheromak. The plasmoid thruster operates by producing FRC-like plasmoids and subsequently ejecting them from the device at a high velocity. The plasmoid is formed inside of a single-turn conical theta-pinch coil. As this process is inductive, there are no electrodes. Similar experiments have yielded plasmoid velocities of at least 50 km/s, and calculations indicate that velocities in excess of 100 km/s should be possible. This concept should be capable of producing Isp's in the range of 5,000 - 15,000 s with thrust densities on the order of 10(exp 5) N per square meters. The current experiment is designed to produce jet powers in the range of 5 - 10 kW, although the concept should be scalable to several MW's. The plasmoid mass and velocity will be measured with a variety of diagnostics, including internal and external B-dot probes, flux loops, Langmuir probes, high-speed cameras and a laser interferometer. Also of key importance will be measurements of the efficiency and mass utilization. Simulations of the plasmoid thruster using MOQUI, a time-dependent MHD code, will be carried out concurrently with experimental testing.

Koelfgen, Syri J.↗
NASA NTRS2002

EVIDENCE OF POSTERUPTION RECONNECTION ASSOCIATED WITH CORONAL MASS EJECTIONS IN THE SOLAR WIND

Using a coupled 2.5-dimensional, time-dependent MHD model of the solar corona and inner heliosphere, we have simulated the eruption and evolution of a coronal mass ejection containing a flux rope all the way from the Sun to 1 AU. Although idealized, we find that the simulation reproduces many generic features of magnetic clouds. In this paper we report on a new, intriguing aspect of these comparisons. Specifically, the results suggest that jetted outflow, driven by posteruptive reconnection underneath the flux rope, occurs and may remain intact out to 1 AU and beyond. We present an example of a magnetic cloud with precisely these signatures and show that the velocity perturbations are consistent with reconnection outflow. We suggest that other velocity and/or density enhancements observed trailing magnetic clouds may be signatures of such reconnection and, in some cases, may not be associated with prominence material, as has previously been suggested.

Riley, Pete↗
NASA NTRS2010

Broken Ergodicity in MHD Turbulence

Ideal magnetohydrodynamic (MHD) turbulence may be represented by finite Fourier series, where the inherent periodic box serves as a surrogate for a bounded astrophysical plasma. Independent Fourier coefficients form a canonical ensemble described by a Gaussian probability density function containing a Hermitian covariance matrix with positive eigenvalues. The eigenvalues at lowest wave number can be very small, resulting in a large-scale coherent structure: a turbulent dynamo. This is seen in computations and a theoretical explanation in terms of 'broken ergodicity' contains Taylor s theory of force-free states. An important problem for future work is the case of real, i.e., dissipative flows. In real flows, broken ergodicity and coherent structure are still expected to occur in MHD turbulence at the largest scale, as suggested by low resolution simulations. One challenge is to incorporate coherent structure at the largest scale into the theory of turbulent fluctuations at smaller scales.

Shebalin, John V.↗
NASA NTRS2010

Multiscale Simulations of Magnetic Island Coalescence

We describe a new interactive parallel Adaptive Mesh Refinement (AMR) framework written in the Python programming language. This new framework, PyAMR, hides the details of parallel AMR data structures and algorithms (e.g., domain decomposition, grid partition, and inter-process communication), allowing the user to focus on the development of algorithms for advancing the solution of a systems of partial differential equations on a single uniform mesh. We demonstrate the use of PyAMR by simulating the pairwise coalescence of magnetic islands using the resistive Hall MHD equations. Techniques for coupling different physics models on different levels of the AMR grid hierarchy are discussed.

Dorelli, John C.↗
NASA NTRS1994

Particle/fluid simulations of an eruptive flare: Identifying the field-aligned currents responsible for the hard x-rays

While magnetohydrodynamics (MHD) can provide a reasonable description of the overall magnetic reconnection that is believed to drive flares, additional, and often separate processes have to be envoked to in order to explain the electron acceleration that is responsible for many of the observed flare emissions. A new model that incorporates the dynamic coronal current sheets, the reconnection site, and possible electron acceleration processes is developed through the use of two-dimensional particle and modified two-fluid simulations. The specific example of an eruptive flare driven by the coalescence of flux tubes supported by prescribed photospheric current elements is evaluated. It is shown that the electrons and ions have differential trajectories through the coronal current sheet which leads to the development of additonal plasma currents that flow around the surface of the current sheet. These surface currents are explicitly neglected in MHD but they are vital to the flare dynamics because they divert current from the coronal current sheet into the chromosphere, producing an effective resistivity that aids the development of fast reconnection. Because the surface currents are in the plane of the magnetic field, electrons in them experience strong acceleration and can account for the observed hard X-ray emissions. Model predictions are compared with observed time profiles of hard X-ray emissions and Doppler shifts seen in soft X-ray line emissions and are able to account for such features as (1) the asymmetry in the rise and decay time of the hard X-rays, (2) the apparent delay between the largest Doppler shifts and the hard X-ray peak, and (3) the relatively low intensity of the blue-shifted component. The use of particle and fluid simulations is important because it provides different, but complementary treatments of the electron acceleration, the global magnetic morphology, and the flare current system.

Winglee, R. M.↗