Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MHD”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Nonlinear MHD modeling of n = 1 RMP -induced pedestal transport and mode coupling effects on ELM suppression in KSTAR

We report fully suppressing edge-localized modes (ELMs), e.g., with resonant magnetic perturbations (RMPs), is essential to reach and sustain high-performance steady-state H-mode plasmas because large ELMs can significantly reduce the lifetime of divertor components in future tokamak reactors. RMP-driven ELM suppression in KSTAR has been modeled by coupling the neoclassical transport code PENTRC to the nonlinear 3D MHD code JOREK. We have found that the radial transport from the combined effects of the kink-peeling, tearing response, and neoclassical toroidal viscosity can explain the pedestal degradation observed in experiments. In addition, it has been found that the RMP response can increase the inter-ELM heat flux on the lower outer divertor by redistributing the heat transport between the divertor plates. In addition to the degraded pedestal, ELM suppression is also attributable to the RMP-induced mode interactions. While the linear stability of peeling-ballooning mode (PBMs) improves owing to the degraded pedestal, the PBM and RMP interaction increases the spectral transfer between edge harmonics, preventing catastrophic growth and the crash of unstable modes. Here, it turns out that the magnetic islands near the pedestal top can play a vital role in mediating the mode interactions.

mode coupling↗
DOE OSTI2013

Effects of Magnetic Field on the Turbulent Wake of a Cylinder in MHD Channel Flow

Results from a free-surface MHD flow experiment are presented detailing the modi cation of vortices in the wake of a circular cylinder with its axis parallel to the applied magnetic fi eld. Experiments were performed with a Reynolds number near Re ~ 104 as the interaction parameter, N = |j x B| / |ρ (υ • ∇), was increased through unity. By concurrently sampling the downstream fluid velocity at sixteen cross-stream locations in the wake, it was possible to extract an ensemble of azimuthal velocity profi les as a function of radius for vortices shed by the cylinder at varying strengths of magnetic field. Results indicate a signi cant change in vortex radius and rotation as N is increased. The lack of deviations from the vortex velocity pro file at high magnetic fi elds suggests the absence of small-scale turbulent features. By sampling the wake at three locations downstream in subsequent experiments, the decay of the vortices was examined and the effective viscosity was found to decrease as N-049±0.4. This reduction in effective viscosity is due to the modi cation of the small-scale eddies by the magnetic fi eld. The slope of the energy spectrum was observed to change from a k-1.8 power-law at low N to a k-3.5 power-law for N > 1. Together, these results suggest the flow smoothly transitioned to a quasi-two-dimensional state in the range 0 < N < 1.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Combustion Synthesis of Boride-Based Electrode Materials for Magnetohydrodynamic (MHD) Direct Power Extraction

The goal of the present project was to develop an advanced, low-cost manufacturing technique for fabrication of boride-based ultrahigh-temperature ceramics (UHTCs) that possess all the required properties to function as sustainable electrodes in magnetohydrodynamic (MHD) direct power extraction applications. Specifically, the project investigated use of mechanical activation-assisted self-propagating high-temperature synthesis (MASHS) followed by pressureless sintering for the fabrication of ceramic materials based on ZrB 2 and HfB 2 from inexpensive raw materials ZrO 2 , HfO 2 , and B 2 O 3 , with Mg as a reactant and NaCl or MgO as an inert diluent.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2020

A Virtual Element Method for the Full MHD system

In this manuscript we present a novel discretization for the incompressible MHD system. Our approach follows the framework of the Virtual Element Method and offers two main advantages. The method can be implemented in unstructured meshes making it highly versatile and capable of handling a wide array of problems involving interfaces, free-boundaries or adaptive refinements on the mesh. The second advantage involves the divergence of the magnetic field, our approach guarantees that it remains solenoidal. We include a theoretical proof of the condition on the magnetic field as well as energy estimates and a well-posedness study. The latter sheds light as to the stability properties of the method.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2014

MHD Modeling of Coronal Loops: the Transition Region Throat

Context. The expansion of coronal loops in the transition region may considerably influence the diagnostics of the plasma emission measure. The cross-sectional area of the loops is expected to depend on the temperature and pressure, and might be sensitive to the heating rate. Aims. The approach here is to study the area response to slow changes in the coronal heating rate, and check the current interpretation in terms of steady heating models. Methods. We study the area response with a time-dependent 2D magnetohydrodynamic (MHD) loop model, including the description of the expanding magnetic field, coronal heating and losses by thermal conduction, and radiation from optically thin plasma. We run a simulation for a loop 50 Mm long and quasi-statically heated to about 4 millikelvin. Results. We find that the area can change substantially with the quasi-steady heating rate, e.g., by approx. 40% at 0.5 millikelvin as the loop temperature varies between 1 millikelvin and 4 millikelvin, and, therefore, affects the interpretation of the differential emission measure vs. temperature (DEM(T)) curves.

magnetohydrodynamics (MHD)↗
NASA NTRS2014

MHD Modelling of Coronal Loops: Injection of High-Speed Chromospheric Flows

Context. Observations reveal a correspondence between chromospheric type II spicules and bright upward-moving fronts in the corona observed in the extreme-ultraviolet (EUV) band. However, theoretical considerations suggest that these flows are probably not the main source of heating in coronal magnetic loops. Aims. We investigate the propagation of high-speed chromospheric flows into coronal magnetic flux tubes and the possible production of emission in the EUV band. Methods. We simulated the propagation of a dense 104 K chromospheric jet upward along a coronal loop by means of a 2D cylindrical MHD model that includes gravity, radiative losses, thermal conduction, and magnetic induction. The jet propagates in a complete atmosphere including the chromosphere and a tenuous cool (approximately 0.8 MK) corona, linked through a steep transition region. In our reference model, the jet initial speed is 70 km per second, its initial density is 10(exp 11) per cubic centimeter, and the ambient uniform magnetic field is 10 G. We also explored other values of jet speed and density in 1D and different magnetic field values in 2D, as well as the jet propagation in a hotter (approximately 1.5 MK) background loop. Results. While the initial speed of the jet does not allow it to reach the loop apex, a hot shock-front develops ahead of it and travels to the other extreme of the loop. The shock front compresses the coronal plasma and heats it to about 10(exp 6) K. As a result, a bright moving front becomes visible in the 171 Angstrom channel of the SDO/AIA mission. This result generally applies to all the other explored cases, except for the propagation in the hotter loop. Conclusions. For a cool, low-density initial coronal loop, the post-shock plasma ahead of upward chromospheric flows might explain at least part of the observed correspondence between type II spicules and EUV emission excess.

Sun: chromosphere↗
DOE OSTI2017

Global three-dimensional simulation of Earth's dayside reconnection using a two-way coupled magnetohydrodynamics with embedded particle-in-cell model: initial results: 3D MHD-EPIC simulation of magnetosphere

Here, we perform a three-dimensional (3D) global simulation of Earth's magnetosphere with kinetic reconnection physics to study the flux transfer events (FTEs) and dayside magnetic reconnection with the recently developed magnetohydrodynamics with embedded particle-in-cell model (MHD-EPIC). During the one-hour long simulation, the FTEs are generated quasi-periodically near the subsolar point and move toward the poles. We also find the magnetic field signature of FTEs at their early formation stage is similar to a ‘crater FTE’, which is characterized by a magnetic field strength dip at the FTE center. After the FTE core field grows to a significant value, it becomes an FTE with typical flux rope structure. When an FTE moves across the cusp, reconnection between the FTE field lines and the cusp field lines can dissipate the FTE. The kinetic features are also captured by our model. A crescent electron phase space distribution is found near the reconnection site. A similar distribution is found for ions at the location where the Larmor electric field appears. The lower hybrid drift instability (LHDI) along the current sheet direction also arises at the interface of magnetosheath and magnetosphere plasma. Finally, the LHDI electric field is about 8 mV/m and its dominant wavelength relative to the electron gyroradius agrees reasonably with MMS observations.

58 GEOSCIENCES↗
DOE OSTI2024

HL-2A's ELM cycle simulations by integrating BOUT++'s drift MHD and transport code

A new integrating model has been developed to couple tokamak edge multiscale magnetohydrodynamic (MHD) events and transport simulations, such as edge-localized mode (ELM) cycles. As a proof of principle, we first start from a set of three-field two-fluid model equations, which includes the pressure, current, and vorticity. Here, the equations are separated into the slowly evolving part of the axisymmetric component by taking a time average of the axisymmetric component. The time-averaged fluxes, which are quadratic in fluctuating quantities, act as driven terms for the time-averaged axisymmetric quantities that determine the plasma transport, and therefore the large-scale evolution of the plasma profiles. Then the HL-2A's ELM cycles are simulated using the model. Good agreements of ELM size and pedestal recovery time have been achieved for the solutions obtained from the coupled simulation compared with experiment. For one ELM cycle simulation, the coupled code can achieve a speedup of a factor of up to 30 over standalone code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Effects of Injected Current Streams on MHD Equilibrium Reconstruction of Local Helicity Injection Plasmas in a Spherical Tokamak

Open field line currents are intrinsic to DC helicity injection plasma startup and pose a challenge for inferring the plasma equilibrium with standard reconstruction analysis. Local helicity injection (LHI) is a type of DC helicity injection which uses small, modular current sources to drive force-free current along helical field lines to produce tokamak plasmas. MHD modeling and magnetic measurements during LHI indicate the injected current streams remain coherent as helical structures on the outboard edge of a core toroidal plasma that is tokamak-like in a toroidally averaged sense. To extract core plasma equilibrium properties, external magnetic diagnostics corrected for contributions from the injected current streams are fitted by a standard Grad-Shafranov equilibrium code. An iterative approach for estimating and subtracting the stream contributions from the diagnostic signals is described and applied to a model equilibrium database to reduce systematic errors introduced by the streams. Convergence is usually attained with 2 to 4 iterations, with derived equilibrium parameters matching the prescribed axisymmetric core values to within estimated experimental uncertainties. Accurate recovery of core parameters occurs when the ratio of the net toroidal windup current from the streams to the core plasma current is less than 0.2, which is typically satisfied in most experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Finite element solution of MHD power-law fluid with slip velocity effect and non-uniform heat source/sink

The laminar steady boundary layer flow of an incompressible MHD Power-law fluid past a continuously moving surface is investigated numerically. The study involves the influence of surface slip and non-uniform heat source/sink on the flow and heat transfer. The governing boundary layer flow equations are transformed into non-dimensional, non-linear coupled ordinary differential equations with the help of suitable similarity transformations. The Galerkin finite element method is implemented to crack the resulting system. The impact of different involved physical parameters is exhibited on the dimensionless velocity profile, temperature distributions and rate of heat transfer in graphical and tabular forms for pseudoplastic and dilatant fluids. The local Nusselt number is found to be the decreasing function of slip parameter, temperature and space-dependent heat sink parameter whereas it increases with increasing values of temperature and space dependent heat source parameter. The problem has important application in attaining the sustainable heat transfer rate for the cooling of fluids, especially in heat exchangers used frequently in chemical industry, in order to increase the trustworthiness of a system, as it removes high heat loads from these systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

A scalable, fully implicit algorithm for the reduced two-field low-β extended MHD model

Here, we demonstrate a scalable fully implicit algorithm for the two-field low-β extended MHD model. This reduced model describes plasma behavior in the presence of strong guide fields, and is of significant practical impact both in nature and in laboratory plasmas. The model displays strong hyperbolic behavior, as manifested by the presence of fast dispersive waves, which make a fully implicit treatment very challenging. In this study, we employ a Jacobian-free Newton–Krylov nonlinear solver, for which we propose a physics-based preconditioner that renders the linearized set of equations suitable for inversion with multigrid methods. As a result, the algorithm is shown to scale both algorithmically (i.e., the iteration count is insensitive to grid refinement and timestep size) and in parallel in a weak-scaling sense, with the wall-clock time scaling weakly with the number of cores for up to 4096 cores. For a 4096 × 4096 mesh, we demonstrate a wall-clock-time speedup of ~6700 with respect to explicit algorithms. The model is validated linearly (against linear theory predictions) and nonlinearly (against fully kinetic simulations), demonstrating excellent agreement.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Gradient-based optimization of 3D MHD equilibria

Using recently developed adjoint methods for computing the shape derivatives of functions that depend on magnetohydrodynamic (MHD) equilibria (Antonsen et al. , J. Plasma Phys. , vol. 85, issue 2, 2019; Paul et al. , J. Plasma Phys. , vol. 86, issue 1, 2020), we present the first example of analytic gradient-based optimization of fixed-boundary stellarator equilibria. We take advantage of gradient information to optimize figures of merit of relevance for stellarator design, including the rotational transform, magnetic well and quasi-symmetry near the axis. With the application of the adjoint method, we reduce the number of equilibrium evaluations by the dimension of the optimization space ( ${\sim }50\text {--}500$ ) in comparison with a finite-difference gradient-based method. We discuss regularization objectives of relevance for fixed-boundary optimization, including a novel method that prevents self-intersection of the plasma boundary. We present several optimized equilibria, including a vacuum field with very low magnetic shear throughout the volume.

Physics↗
DOE OSTI2019

Studying Dawn‐Dusk Asymmetries of Mercury's Magnetotail Using MHD‐EPIC Simulations

Abstract MESSENGER has observed a lot of dawn‐dusk asymmetries in Mercury's magnetotail, such as the asymmetries of the cross‐tail current sheet thickness and the occurrence of flux ropes, dipolarization events, and energetic electron injections. In order to obtain a global pictures of Mercury's magnetotail dynamics and the relationship between these asymmetries, we perform global simulations with the magnetohydrodynamics with embedded particle‐in‐cell (MHD‐EPIC) model, where Mercury's magnetotail region is covered by a PIC code. Our simulations show that the dawnside current sheet is thicker, the plasma density is larger, and the electron pressure is higher than the duskside. Under a strong interplanetary magnetic field driver, the simulated reconnection sites prefer the dawnside. We also found the dipolarization events and the planetward electron jets are moving dawnward while they are moving toward the planet, so that almost all dipolarization events and high‐speed plasma flows concentrate in the dawn sector. The simulation results are consistent with MESSENGER observations.

Chen, Yuxi↗
DOE OSTI2012

Time-dependent MHD simulations of the solar wind outflow using interplanetary scintillation observations

Numerical modeling of the heliosphere is a critical component of space weather forecasting. The accuracy of heliospheric models can be improved by using realistic boundary conditions and confirming the results with in situ spacecraft measurements. To accurately reproduce the solar wind (SW) plasma flow near Earth, we need realistic, time-dependent boundary conditions at a fixed distance from the Sun. We may prepare such boundary conditions using SW speed and density determined from interplanetary scintillation (IPS) observations, magnetic field derived from photospheric magnetograms, and temperature estimated from its correlation with SW speed. In conclusion, we present here the time-dependent MHD simulation results obtained by using the 2011 IPS data from the Solar-Terrestrial Environment Laboratory as time-varying inner boundary conditions and compare the simulated data at Earth with OMNI data (spacecraft-interspersed, near-Earth solar wind data).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Real time MHD mode control using ECCD in KSTAR: Plan and requirements

For a high-performance, advanced tokamak mode in KSTAR, we have been developing a real-time control system of MHD modes such as sawtooth and Neo-classical Tearing Mode (NTM) by ECH/ECCD. The active feedback control loop will be also added to the mirror position and the real-time detection of the mode position. In this year, for the stabilization of NTM that is crucial to plasma performance we have implemented open-loop ECH antenna control system in KSTAR Plasma Control System (PCS) for ECH mirror movement during a single plasma discharge. KSTAR 170 GHz ECH launcher which was designed and fabricated by collaboration with PPPL and POSTECH has a final mirror of a poloidally and toroidally steerable mirror. The poloidal steering motion is only controlled in the real-time NTM control system and its maximum steering speed is 10 degree/sec by DC motor. However, the latency of the mirror control system and the return period of ECH antenna mirror angle are not fast because the existing launcher mirror control system is based on PLC which is connected to the KSTAR machine network through serial to LAN converter. In this paper, we present the design of real time NTM control system, ECH requirements, and the upgrade plan.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

MHD simulation of RF current drive in MST

Auxiliary heating and current drive using RF waves such as the electron Bernstein wave (EBW) promises to advance the performance of the reversed field pinch (RFP). In previous computational work [1], a hypothetical edge-localized current drive is shown to suppress the tearing activity which governs the macroscopic transport properties of the RFP. The ideal conditions for tearing stabilization include a reduced toroidal induction, and precise width and radial position of the Gaussian-shaped external current drive. In support of the EBW experiment on the Madison Symmetric Torus, an integrated modeling scheme now incorporates ray tracing and Fokker-Plank predictions of auxiliary current into single fluid MHD. Simulations at low Lundquist number (S ∼ 10{sup 4}) generally agree with the previous work; significantly more burdensome simulations at MST-like Lundquist number (S ∼ 3×10{sup 6}) show unexpected results. The effect on nonlinearly saturated current profile by a particular RF-driven external force decreases in magnitude and widens considerably as the Lundquist number increases toward experimental values. Simulations reproduce the periodic current profile relaxation events observed in experiment (sawteeth) in the absence of current profile control. Reduction of the tearing mode amplitudes is still observable; however, reduction is limited to periods between the large bursts of magnetic activity at each sawtooth. The sawtoothing pattern persists with up to 10 MW of externally applied RF power. Periods with prolonged low tearing amplitude are predicted with a combination of external current drive and a reduced toroidal loop voltage, consistent with previous conclusions. Finally, the resistivity profile is observed to have a strong effect on the optimal externally driven current profile for mode stabilization.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Quasi 3D ECE imaging system for study of MHD instabilities in KSTAR

A second electron cyclotron emission imaging (ECEI) system has been installed on the KSTAR tokamak, toroidally separated by 1/16th of the torus from the first ECEI system. For the first time, the dynamical evolutions of MHD instabilities from the plasma core to the edge have been visualized in quasi-3D for a wide range of the KSTAR operation (B{sub 0} = 1.7∼3.5 T). This flexible diagnostic capability has been realized by substantial improvements in large-aperture quasi-optical microwave components including the development of broad-band polarization rotators for imaging of the fundamental ordinary ECE as well as the usual 2nd harmonic extraordinary ECE.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗