Flow curvature effects on the parallel velocity shear driven instability: MHD simulations
Not provided.
SEARCH · Engineering Papers
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.
Not provided.
Not provided.
The solar wind emanates from the hot and tenuous solar corona. Earlier studies using 1.5-dimensional simulations show that Alfven waves generated in the photosphere play an important role in coronal heating through the process of nonlinear mode conversion. In order to understand the physics of coronal heating and solar wind acceleration together, it is important to consider the regions from photosphere to interplanetary space as a single system. We performed 2.5-dimensional, self-consistent magnetohydrodynamic simulations, covering from the photosphere to the interplanetary space for the first time. We carefully set up the grid points with spherical coordinates to treat the Alfven waves in the atmosphere with huge density contrast and successfully simulate the solar wind streaming out from the hot solar corona as a result of the surface convective motion. The footpoint motion excites Alfven waves along an open magnetic flux tube, and these waves traveling upward in the non-uniform medium undergo wave reflection, nonlinear mode conversion from Alfven mode to slow mode, and turbulent cascade. These processes lead to the dissipation of Alfven waves and acceleration of the solar wind. It is found that the shock heating by the dissipation of the slow-mode wave plays a fundamental role in the coronal heating process, whereas the turbulent cascade and shock heating drive the solar wind.
We present a pair of three-dimensional magnetohydrodynamical simulations of intermittent jets from a central active galactic nucleus (AGN) in a galaxy cluster extracted from a high-resolution cosmological simulation. The selected cluster was chosen as an apparently relatively relaxed system, not having undergone a major merger in almost 7 Gyr. Despite this characterization and history, the intracluster medium (ICM) contains quite active 'weather'. We explore the effects of this ICM weather on the morphological evolution of the AGN jets and lobes. The orientation of the jets is different in the two simulations so that they probe different aspects of the ICM structure and dynamics. We find that even for this cluster, which can be characterized as relaxed by an observational standard, the large-scale, bulk ICM motions can significantly distort the jets and lobes. Synthetic X-ray observations of the simulations show that the jets produce complex cavity systems, while synthetic radio observations reveal bending of the jets and lobes similar to wide-angle tail radio sources. The jets are cycled on and off with a 26 Myr period using a 50% duty cycle. This leads to morphological features similar to those in 'double-double' radio galaxies. While the jet and ICM magnetic fields are generally too weak in the simulations to play a major role in the dynamics, Maxwell stresses can still become locally significant.
Alfven waves may be difficult to excite at the photosphere due to low-ionization fraction and suffer near-total reflection at the transition region (TR). Yet they are ubiquitous in the corona and heliosphere. To overcome these difficulties, we show that they may instead be generated high in the chromosphere by conversion from reflecting fast magnetohydrodynamic waves, and that Alfvenic TR reflection is greatly reduced if the fast reflection point is within a few scale heights of the TR. The influence of mode conversion on the phase of the reflected fast wave is also explored. This phase can potentially be misinterpreted as a travel speed perturbation with implications for the practical seismic probing of active regions.
Stellar winds are believed to be the dominant factor in the spin-down of stars over time. However, stellar winds of solar analogs are poorly constrained due to observational challenges. In this paper, we present a grid of magnetohydrodynamic models to study and quantify the values of stellar mass loss and angular momentum loss rates as a function of the stellar rotation period, magnetic dipole component, and coronal base density. We derive simple scaling laws for the loss rates as a function of these parameters, and constrain the possible mass loss rate of stars with thermally driven winds. Despite the success of our scaling law in matching the results of the model, we find a deviation between the 'solar dipole' case and a real case based on solar observations that overestimates the actual solar mass loss rate by a factor of three. This implies that the model for stellar fields might require a further investigation with additional complexity. Mass loss rates in general are largely controlled by the magnetic field strength, with the wind density varying in proportion to the confining magnetic pressure B {sup 2}. We also find that the mass loss rates obtained using our grid models drop much faster with the increase in rotation period than scaling laws derived using observed stellar activity. For main-sequence solar-like stars, our scaling law for angular momentum loss versus poloidal magnetic field strength retrieves the well-known Skumanich decline of angular velocity with time, Ω{sub *}∝t {sup –1/2}, if the large-scale poloidal magnetic field scales with rotation rate as B{sub p}∝Ω{sub ⋆}{sup 2}.
Giant flares on soft gamma-ray repeaters that are thought to take place on magnetars release enormous energy in a short time interval. Their power can be explained by catastrophic instabilities occurring in the magnetic field configuration and the subsequent magnetic reconnection. By analogy with the coronal mass ejection events on the Sun, we develop a theoretical model via an analytic approach for magnetar giant flares. In this model, the rotation and/or displacement of the crust causes the field to twist and deform, leading to flux rope formation in the magnetosphere and energy accumulation in the related configuration. When the energy and helicity stored in the configuration reach a threshold, the system loses its equilibrium, the flux rope is ejected outward in a catastrophic way, and magnetic reconnection helps the catastrophe develop to a plausible eruption. By taking SGR 1806–20 as an example, we calculate the free magnetic energy released in such an eruptive process and find that it is more than 10{sup 47} erg, which is enough to power a giant flare. The released free magnetic energy is converted into radiative energy, kinetic energy, and gravitational energy of the flux rope. We calculated the light curves of the eruptive processes for the giant flares of SGR 1806–20, SGR 0526–66, and SGR 1900+14, and compared them with the observational data. The calculated light curves are in good agreement with the observed light curves of giant flares.
Measurements from the Spectral Irradiance Monitor (SIM) on board the Solar Radiation and Climate Experiment mission indicate that solar spectral irradiance at visible and IR wavelengths varies in counter phase with the solar activity cycle. The sign of these variations is not reproduced by most of the irradiance reconstruction techniques based on variations of surface magnetism employed so far, and it is not yet clear whether SIM calibration procedures need to be improved or if instead new physical mechanisms must be invoked to explain such variations. We employ three-dimensional magnetohydrodynamic simulations of the solar photosphere to investigate the dependence of solar radiance in SIM visible and IR spectral ranges on variations of the filling factor of surface magnetic fields. We find that the contribution of magnetic features to solar radiance is strongly dependent on the location on the disk of the features, which are negative close to disk center and positive toward the limb. If features are homogeneously distributed over a region around the equator (activity belt), then their contribution to irradiance is positive with respect to the contribution of HD snapshots, but decreases with the increase of their magnetic flux for average magnetic flux larger than 50 G in at least two of the visible and IR spectral bands monitored by SIM. Under the assumption that the 50 G snapshots are representative of quiet-Sun regions, we thus find that the Spectral Irradiance can be in counter-phase with the solar magnetic activity cycle.
We present the results of axisymmetric magnetohydrodynamic simulations investigating the launching of jets and outflows from a magnetically diffusive accretion disk. The time evolution of the disk structure is self-consistently taken into account. In contrast to previous works, we have applied spherical coordinates for the numerical grid, implying substantial benefits concerning the numerical resolution and the stability of the simulation. Due to the new setup, we were able to run simulations for more than 150,000 dynamical times on a domain extending 1500 inner disk radii with a resolution of up to 24 cells per disk height in the inner disk. Depending on the disk magnetization, jet launching occurs in two different but complementary regimes—jets driven predominantly by centrifugal or magnetic forces. These regimes differ in the ejection efficiency concerning mass, energy, and angular momentum. We show that it is the actual disk magnetization and not so much the initial magnetization that best describes the disk-jet evolution. Considering the actual disk magnetization, we also find that simulations starting with different initial magnetization evolve in a similar—typical—way, due to advection and diffusion, as the magnetic flux in the disk evolves in time. Exploring a new, modified diffusivity model, we confirm the self-similar structure of the global jet-launching disk, obtaining power laws for the radial profiles of the disk physical variables such as density, magnetic field strength, or accretion velocity.
We present the results of a large suite of three-dimensional models of the collapse of magnetic molecular cloud cores using the adaptive mesh refinement code Enzo2.2 in the ideal magnetohydrodynamics approximation. The cloud cores are initially either prolate or oblate, centrally condensed clouds with masses of 1.73 or 2.73 M {sub ☉}, respectively. The radial density profiles are Gaussian, with central densities 20 times higher than boundary densities. A barotropic equation of state is used to represent the transition from low density isothermal phases, to high density optically thick phases. The initial magnetic field strength ranges from 6.3 to 100 μG, corresponding to clouds that are strongly to marginally supercritical, respectively, in terms of the mass to magnetic flux ratio. The magnetic field is initially uniform and aligned with the clouds' rotation axes, with initial ratios of rotational to gravitational energy ranging from 10{sup –4} to 0.1. Two significantly different outcomes for collapse result: (1) formation of single protostars with spiral arms, and (2) fragmentation into multiple protostar systems. The transition between these two outcomes depends primarily on the initial magnetic field strength, with fragmentation occurring for mass to flux ratios greater than about 14 times the critical ratio for prolate clouds. Oblate clouds typically fragment into several times more clumps than prolate clouds. Multiple, rather than binary, system formation is the general rule in either case, suggesting that binary stars are primarily the result of the orbital dissolution of multiple protostar systems.
In this paper we present a numerical study of the time evolution of solar prominences embedded in sheared magnetic arcades. The prominence is represented by a density enhancement in a background-stratified atmosphere and is connected to the photosphere through the magnetic field. By solving the ideal magnetohydrodynamic equations in three dimensions, we study the dynamics for a range of parameters representative of real prominences. Depending on the parameters considered, we find prominences that are suspended above the photosphere, i.e., detached prominences, but also configurations resembling curtain or hedgerow prominences whose material continuously connects to the photosphere. The plasma-β is an important parameter that determines the shape of the structure. In many cases magnetic Rayleigh-Taylor instabilities and oscillatory phenomena develop. Fingers and plumes are generated, affecting the whole prominence body and producing vertical structures in an essentially horizontal magnetic field. However, magnetic shear is able to reduce or even to suppress this instability.
The solar atmosphere is extremely dynamic, and many important phenomena develop on small scales that are unresolved in observations with the Helioseismic and Magnetic Imager (HMI) instrument on the Solar Dynamics Observatory. For correct calibration and interpretation of the observations, it is very important to investigate the effects of small-scale structures and dynamics on the HMI observables, such as Doppler shift, continuum intensity, spectral line depth, and width. We use 3D radiative hydrodynamics simulations of the upper turbulent convective layer and the atmosphere of the Sun, and a spectro-polarimetric radiative transfer code to study observational characteristics of the Fe i 6173 Å line observed by HMI in quiet-Sun regions. We use the modeling results to investigate the sensitivity of the line Doppler shift to plasma velocity, and also sensitivities of the line parameters to plasma temperature and density, and determine effective line formation heights for observations of solar regions located at different distances from the disk center. These estimates are important for the interpretation of helioseismology measurements. In addition, we consider various center-to-limb effects, such as convective blueshift, variations of helioseismic travel-times, and the “concave” Sun effect, and show that the simulations can qualitatively reproduce the observed phenomena, indicating that these effects are related to a complex interaction of the solar dynamics and radiative transfer.
Using a three-dimensional (3D) magnetohydrodynamic model, we simulate the magnetic reconnection in a single current sheet. We assume a finite guide field, a random perturbation on the velocity field, and uniform resistivity. Our model enhances the reconnection rate relative to the classical Sweet–Parker model in the same configuration. The efficiency of magnetic energy conversion is increased by interactions between the multiple tearing layers coexisting in the global current sheet. This interaction, which forms a positive-feedback system, arises from coupling of the inflow and outflow regions in different layers across the current sheet. The coupling accelerates the elementary reconnection events, thereby enhancing the global reconnection rate. The reconnection establishes flux tubes along each tearing layer. Slow-mode shocks gradually form along the outer boundaries of these tubes, further accelerating the magnetic energy conversion. Such a positive-feedback system is absent in two-dimensional simulations, 3D reconnection without a guide field, and reconnection under a single perturbation mode. We refer to our model as the “shock-evoking positive-feedback” model.
Without any doubt, solar flaring loops possess a multithread internal structure that is poorly resolved, and there are no means to observe heating episodes and thermodynamic evolution of the individual threads. These limitations cause fundamental problems in numerical modeling of flaring loops, such as selection of a structure and a number of threads, and an implementation of a proper model of the energy deposition process. A set of one-dimensional (1D) hydrodynamic and two-dimensional (2D) magnetohydrodynamic models of a flaring loop are developed to compare energy redistribution and plasma dynamics in the course of a prototypical solar flare. Basic parameters of the modeled loop are set according to the progenitor M1.8 flare recorded in AR 10126 on 2002 September 20 between 09:21 UT and 09:50 UT. The nonideal 1D models include thermal conduction and radiative losses of the optically thin plasma as energy-loss mechanisms, while the nonideal 2D models take into account viscosity and thermal conduction as energy-loss mechanisms only. The 2D models have a continuous distribution of the parameters of the plasma across the loop and are powered by varying in time and space along and across the loop heating flux. We show that such 2D models are an extreme borderline case of a multithread internal structure of the flaring loop, with a filling factor equal to 1. Nevertheless, these simple models ensure the general correctness of the obtained results and can be adopted as a correct approximation of the real flaring structures.
Explore the source record for details and available documents.
Abstract not provided
One of the primary tools foreseen for heating, current drive (CD) and q-profile control in future fusion reactors such as ITER and a Fusion Nuclear Science Facility is the neutral beam injection (NBI). However, fast ions from NBI may also provide the drive for energetic particle-driven instabilities (e.g. Alfvénic modes (AEs)), which in turn redistribute fast ions in both space and energy, thus hampering the control capabilities and overall efficiency of NB-driven current. Based on experiments on the NSTX tokamak (M. Ono et al 2000 Nucl. Fusion 40 557), the effects of AEs and other low-frequency magneto-hydrodynamic instabilities on NB-CD efficiency are investigated. When looking at the new fast ion transport model, which accounts for particle transport in phase space as required for resonant AE perturbations, is utilized to obtain consistent simulations of NB-CD through the tokamak transport code TRANSP. It is found that instabilities do indeed reduce the NB-driven current density over most of the plasma radius by up to ~50%. Moreover, the details of the current profile evolution are sensitive to the specific model used to mimic the interaction between NB ions and instabilities. Finally, implications for fast ion transport modeling in integrated tokamak simulations are briefly discussed.
We introduce CAFE, a new independent code designed to solve the equations of relativistic ideal magnetohydrodynamics (RMHD) in three dimensions. We present the standard tests for an RMHD code and for the relativistic hydrodynamics regime because we have not reported them before. The tests include the one-dimensional Riemann problems related to blast waves, head-on collisions of streams, and states with transverse velocities, with and without magnetic field, which is aligned or transverse, constant or discontinuous across the initial discontinuity. Among the two-dimensional (2D) and 3D tests without magnetic field, we include the 2D Riemann problem, a one-dimensional shock tube along a diagonal, the high-speed Emery wind tunnel, the Kelvin–Helmholtz (KH) instability, a set of jets, and a 3D spherical blast wave, whereas in the presence of a magnetic field we show the magnetic rotor, the cylindrical explosion, a case of Kelvin–Helmholtz instability, and a 3D magnetic field advection loop. The code uses high-resolution shock-capturing methods, and we present the error analysis for a combination that uses the Harten, Lax, van Leer, and Einfeldt (HLLE) flux formula combined with a linear, piecewise parabolic method and fifth-order weighted essentially nonoscillatory reconstructors. We use the flux-constrained transport and the divergence cleaning methods to control the divergence-free magnetic field constraint.