Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical 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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2012

NUMERICAL SIMULATIONS OF THE MAGNETIC RAYLEIGH-TAYLOR INSTABILITY IN THE KIPPENHAHN-SCHLUeTER PROMINENCE MODEL. I. FORMATION OF UPFLOWS

The launch of the Hinode satellite led to the discovery of rising plumes, dark in chromospheric lines, that propagate from large ({approx}10 Mm) bubbles that form at the base of quiescent prominences. The plumes move through a height of approximately 10 Mm while developing highly turbulent profiles. The magnetic Rayleigh-Taylor instability was hypothesized to be the mechanism that drives these flows. In this study, using three-dimensional (3D) MHD simulations, we investigate the nonlinear stability of the Kippenhahn-Schlueter prominence model for the interchange mode of the magnetic Rayleigh-Taylor instability. The model simulates the rise of a buoyant tube inside the quiescent prominence model, where the interchange of magnetic field lines becomes possible at the boundary between the buoyant tube and the prominence. Hillier et al. presented the initial results of this study, where upflows of constant velocity (maximum found 6 km s{sup -1}) and a maximum plume width Almost-Equal-To 1.5 Mm which propagate through a height of approximately 6 Mm were found. Nonlinear interaction between plumes was found to be important for determining the plume dynamics. In this paper, using the results of ideal MHD simulations, we determine how the initial parameters for the model and buoyant tube affect the evolution of instability. We find that the 3D mode of the magnetic Rayleigh-Taylor instability grows, creating upflows aligned with the magnetic field of constant velocity (maximum found 7.3 km s{sup -1}). The width of the upflows is dependent on the initial conditions, with a range of 0.5-4 Mm which propagate through heights of 3-6 Mm. These results are in general agreement with the observations of the rising plumes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

NUMERICAL SIMULATIONS OF SUPERNOVA DUST DESTRUCTION. II. METAL-ENRICHED EJECTA KNOTS

Following our previous work, we investigate through hydrodynamic simulations the destruction of newly formed dust grains by sputtering in the reverse shocks of supernova remnants. Using an idealized setup of a planar shock impacting a dense, spherical clump, we implant a population of Lagrangian particles into the clump to represent a distribution of dust grains in size and composition. We vary the relative velocity between the reverse shock and ejecta clump to explore the effects of shock heating and cloud compression. Because supernova ejecta will be metal-enriched, we consider gas metallicities from Z/Z{sub Sun} = 1 to 100 and their influence on the cooling properties of the cloud and the thermal sputtering rates of embedded dust grains. We post-process the simulation output to calculate grain sputtering for a variety of species and size distributions. In the metallicity regime considered in this paper, the balance between increased radiative cooling and increased grain erosion depends on the impact velocity of the reverse shock. For slow shocks (v{sub shock} {<=} 3000 km s{sup -1}), the amount of dust destruction is comparable across metallicities or in some cases is decreased with increased metallicity. For higher shock velocities (v{sub shock} {>=} 5000 km s{sup -1}), an increase in metallicity from Z/Z{sub Sun} = 10 to 100 can lead to an additional 24% destruction of the initial dust mass. While the total dust destruction varies widely across grain species and simulation parameters, our most extreme cases result in complete destruction for some grain species and only 44% dust mass survival for the most robust species. These survival rates are important in understanding how early supernovae contribute to the observed dust masses in high-redshift galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

ANOMALOUS TRANSPORT OF HIGH-ENERGY COSMIC RAYS IN GALACTIC SUPERBUBBLES. I. NUMERICAL SIMULATIONS

We present a simple continuous-time random-walk model for the transport of energetic particles accelerated by a collection of supernova explosions in a galactic superbubble, developed to simulate and highlight signatures of anomalous transport on the particles' evolution and their spectra in a multi-shock context. We assume standard diffusive shock acceleration (DSA) theory for each shock encounter. The superbubble (an OB stars association) is idealized as a heterogeneous region of particle sources and sinks bounded by a random surface. The model is based on two coupled stochastic differential equations and is applied for protons and alpha particles. Using characteristic values for a typical bubble, our simulations suggest that acceleration and transport in the bubble may be sub-diffusive. In addition, a spectral break in the particles' evolution and spectra is evident located at Almost-Equal-To 10{sup 15} eV for protons and Almost-Equal-To 3 Multiplication-Sign 10{sup 15} eV for alphas. Our simulations are consistent with a bubble's mean magnetic field strength of Almost-Equal-To 1 {mu}G and a shock separation distance {approx}0.1 Multiplication-Sign the characteristic radius of the bubble. The simulations imply that the diffusion coefficient (for the elementary shock acceleration process) is {approx}< 10{sup 27} cm{sup 2} s{sup -1} at 1 GeV/c. While the sub-diffusive transport is readily attributed to the stochastic nature of the acceleration time according to DSA theory, the spectral break appears to be an artifact of transport in a finite medium. These simulations point to a new and intriguing phenomenon associated with the statistical nature of collective acceleration of high-energy cosmic rays in galactic superbubbles.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

A RECIPE TO PROBE ALTERNATIVE THEORIES OF GRAVITATION VIA N-BODY NUMERICAL SIMULATIONS. I. SPIRAL GALAXIES

A way to probe alternative theories of gravitation is to study if they could account for the structures of the universe. We therefore modified the well-known Gadget-2 code to probe alternative theories of gravitation through galactic dynamics. As an application, we simulate the evolution of spiral galaxies to probe alternative theories of gravitation whose weak field limits have a Yukawa-like gravitational potential. These simulations show that galactic dynamics can be used to constrain the parameters associated with alternative theories of gravitation. It is worth stressing that the recipe given in this study can be applied to any other alternative theory of gravitation in which the superposition principle is valid.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

NUMERICAL SIMULATION OF SOLAR MICROFLARES IN A CANOPY-TYPE MAGNETIC CONFIGURATION

Microflares are small activities in the solar low atmosphere; some are in the low corona while others are in the chromosphere. Observations show that some of the microflares are triggered by magnetic reconnection between the emerging flux and a pre-existing background magnetic field. We perform 2.5-dimensional, compressible, resistive magnetohydrodynamic simulations of the magnetic reconnection with gravity considered. The background magnetic field is a canopy-type configuration that is rooted at the boundary of the solar supergranule. By changing the bottom boundary conditions in the simulation, a new magnetic flux emerges at the center of the supergranule and reconnects with the canopy-type magnetic field. We successfully simulate the coronal and chromospheric microflares whose current sheets are located at the corona and the chromosphere, respectively. The microflare with a coronal origin has a larger size and a higher temperature enhancement than the microflare with a chromospheric origin. In the microflares with coronal origins, we also found a hot jet ({approx}1.8 Multiplication-Sign 10{sup 6} K), which is probably related to the observational extreme ultraviolet or soft X-ray jets, and a cold jet ({approx}10{sup 4} K), which is similar to the observational H{alpha}/Ca surges. However, there is only a H{alpha}/Ca bright point in the microflares that have chromospheric origins. The study of parameter dependence shows that the size and strength of the emerging magnetic flux are the key parameters that determine the height of the reconnection location, and they further determine the different observational features of the microflares.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

NUMERICAL SIMULATIONS OF THE MAGNETIC RAYLEIGH-TAYLOR INSTABILITY IN THE KIPPENHAHN-SCHLUeTER PROMINENCE MODEL. II. RECONNECTION-TRIGGERED DOWNFLOWS

The launch of the Hinode satellite has allowed high-resolution observations of supersonic bright downflows in quiescent prominences, known as prominence knots. We present observations in the Ca II H spectral line using the Solar Optical Telescope on board the Hinode satellite of a descending plasma knot of size {approx}900 km. The knot initially undergoes ballistic motion before undergoing impulsive accelerations at the same time as experiencing increases in intensity. We also present a subset of our three-dimensional magnetohydrodynamic simulations, performed to investigate the nonlinear stability of the Kippenhahn-Shlueter prominence model to the magnetic Rayleigh-Taylor instability in which interchange reconnection occurs. The interchange reconnection in the model breaks the force balance along the field lines which initiates the downflows. The downflows propagate with a downward fluid velocity of {approx}15 km s{sup -1} and a characteristic size of {approx}700 km. We conclude that the observed plasma blob and the simulated downflow are driven by the breaking of the force balance along the magnetic field as a result of a change in magnetic topology caused by reconnection of the magnetic field.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

NONRESONANT INTERACTION OF CHARGED ENERGETIC PARTICLES WITH LOW-FREQUENCY NONCOMPRESSIVE TURBULENCE: NUMERICAL SIMULATION

A new method for simulating the three-dimensional dynamics of charged energetic particles in very broadband noncompressive magnetic turbulence is introduced. All scales within the primary inertial range of the turbulence observed in the solar wind near 1 AU are now included for the independent computations of both the particle dynamics and the turbulent magnetic field lines (MFLs). While previous theories of resonant particle pitch-angle (PA) scattering and transport in interplanetary magnetic fields had favored interpreting the observed depletions in the electron PA distributions (PADs) around 90 Degree-Sign PA as evidence of poor scattering at low PA cosines, the computed particle dynamics reveal a very different reality. The MFL directions now vary on many scales, and the PADs are depleted around 90 Degree-Sign PA due to nonresonant filtering of the particles that propagate at too large an angle to the local magnetic field. Rather than being too weak, the scattering through 90 Degree-Sign PA is actually so strong that the particles (electrons and protons/ions) are reflected and trapped in the turbulent magnetic fields. While the low-frequency nonresonant turbulence produces ubiquitous magnetic traps that only let through particles with the most field-aligned velocities, higher-frequency near-gyroscale turbulence, when present, enhances particle transport by allowing the particles to navigate between magnetic traps. Finally, visualizing both particle trajectories and MFLs in the very same turbulence reveals a powerful tool for understanding the effects of the turbulent fields on the particle dynamics and cross-field transport. Some cross-field-line scattering, strongly amplified by MFL dispersal, results in a strong cross-field scattering of the particles. From this visualization, it also appears that near-gyroscale turbulence, previously known as gyroresonant turbulence, does not resonantly interact with the particles. The interaction between particles and fields at or near the gyroscale, though potentially strong, does not actually involve the periodic driving of a true resonance.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

DIRECT NUMERICAL SIMULATION OF RADIATION PRESSURE-DRIVEN TURBULENCE AND WINDS IN STAR CLUSTERS AND GALACTIC DISKS

The pressure exerted by the radiation of young stars may be an important feedback mechanism that drives turbulence and winds in forming star clusters and the disks of starburst galaxies. However, there is great uncertainty in how efficiently radiation couples to matter in these high optical depth environments. In particular, it is unclear what levels of turbulence the radiation can produce, and whether the infrared radiation trapped by the dust opacity can give rise to heavily mass-loaded winds. In this paper, we report a series of two-dimensional flux-limited diffusion radiation-hydrodynamics calculations performed with the code ORION in which we drive strong radiation fluxes through columns of dusty matter confined by gravity in order to answer these questions. We consider both systems where the radiation flux is sub-Eddington throughout the gas column, and those where it is super-Eddington at the midplane but sub-Eddington in the atmosphere. In the latter, we find that the radiation-matter interaction gives rise to radiation-driven Rayleigh-Taylor instability, which drives supersonic turbulence at a level sufficient to fully explain the turbulence seen in Galactic protocluster gas clouds, and to make a non-trivial contribution to the turbulence observed in starburst galaxy disks. However, the instability also produces a channel structure in which the radiation-matter interaction is reduced compared to time-steady analytic models because the radiation field is not fully trapped. For astrophysical parameters relevant to forming star clusters and starburst galaxies, we find that this effect reduces the net momentum deposition rate in the dusty gas by a factor of {approx}2-6 compared to simple analytic estimates, and that in steady state the Eddington ratio reaches unity and there are no strong winds. We provide an approximation formula, appropriate for implementation in analytic models and non-radiative simulations, for the force exerted by the infrared radiation field in this regime.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

COMPOSITION STRUCTURE OF INTERPLANETARY CORONAL MASS EJECTIONS FROM MULTISPACECRAFT OBSERVATIONS, MODELING, AND COMPARISON WITH NUMERICAL SIMULATIONS

We present an analysis of the ionic composition of iron for two interplanetary coronal mass ejections (ICMEs) observed on 2007 May 21-23 by the ACE and STEREO spacecraft in the context of the magnetic structure of the ejecta flux rope, sheath region, and surrounding solar wind flow. This analysis is made possible due to recent advances in multispacecraft data interpolation, reconstruction, and visualization as well as results from recent modeling of ionic charge states in MHD simulations of magnetic breakout and flux cancellation coronal mass ejection (CME) initiation. We use these advances to interpret specific features of the ICME plasma composition resulting from the magnetic topology and evolution of the CME. We find that, in both the data and our MHD simulations, the flux ropes centers are relatively cool, while charge state enhancements surround and trail the flux ropes. The magnetic orientations of the ICMEs are suggestive of magnetic breakout-like reconnection during the eruption process, which could explain the spatial location of the observed iron enhancements just outside the traditional flux rope magnetic signatures and between the two ICMEs. Detailed comparisons between the simulations and data were more complicated, but a sharp increase in high iron charge states in the ACE and STEREO-A data during the second flux rope corresponds well to similar features in the flux cancellation results. We discuss the prospects of this integrated in situ data analysis and modeling approach to advancing our understanding of the unified CME-to-ICME evolution.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

NUMERICAL SIMULATION OF THREE-DIMENSIONAL ASYMMETRIC RECONNECTION AND APPLICATION TO A PHYSICAL MECHANISM OF PENUMBRAL MICROJETS

Three-dimensional (3D) component reconnection, where reconnecting field lines are not perfectly anti-parallel, is studied with a 3D magnetohydrodynamic simulation. In particular, we consider the asymmetry of the field strength of the reconnecting field lines. As the asymmetry increases, the generated reconnection jet tends to be parallel to stronger field lines. This is because weaker field lines have higher gas pressure in the initial equilibrium, and hence the gas pressure gradient along the reconnected field lines is generated, which accelerates the field-aligned plasma flow. This mechanism may explain penumbral microjets and other types of jets that are parallel to magnetic field lines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

NUMERICAL SIMULATIONS OF SPICULE ACCELERATION

Observations in the H{alpha} line of hydrogen and the H and K lines of singly ionized calcium on the solar limb reveal the existence of structures with jet-like behavior, usually designated as spicules. The driving mechanism for such structures remains poorly understood. Sterling et al. shed some light on the problem mimicking reconnection events in the chromosphere with a one-dimensional code by injecting energy with different spatial and temporal distributions and tracing the thermodynamic evolution of the upper chromospheric plasma. They found three different classes of jets resulting from these injections. We follow their approach but improve the physical description by including non-LTE cooling in strong spectral lines and non-equilibrium hydrogen ionization. Increased cooling and conversion of injected energy into hydrogen ionization energy instead of thermal energy both lead to weaker jets and smaller final extent of the spicules compared with Sterling et al. In our simulations we find different behavior depending on the timescale for hydrogen ionization/recombination. Radiation-driven ionization fronts also form.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

NUMERICAL SIMULATIONS OF Z-PINCH EXPERIMENTS TO CREATE SUPERSONIC DIFFERENTIALLY ROTATING PLASMA FLOWS

The physics of accretion disks is of fundamental importance for understanding of a wide variety of astrophysical sources that includes protostars, X-ray binaries, and active galactic nuclei. The interplay between hydrodynamic flows and magnetic fields and the potential for turbulence-producing instabilities is a topic of active research that would benefit from the support of dedicated experimental studies. Such efforts are in their infancy, but in an effort to push the enterprise forward we propose an experimental configuration which employs a modified cylindrical wire array Z-pinch to produce a rotating plasma flow relevant to accretion disks. We present three-dimensional resistive magnetohydrodynamic simulations which show how this approach can be implemented. In the simulations, a rotating plasma cylinder or ring is formed, with typical rotation velocity {approx}30 km s{sup -1}, Mach number {approx}4, and Reynolds number in excess of 10{sup 7}. The plasma is also differentially rotating. Implementation of different external magnetic field configurations is discussed. It is found that a modest uniform vertical field of 1 T can affect the dynamics of the system and could be used to study magnetic field entrainment and amplification through differential rotation. A dipolar field potentially relevant to the study of accretion columns is also considered.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

NUMERICAL SIMULATIONS OF THE POSSIBLE ORIGIN OF THE TWO SUB-PARSEC SCALE AND COUNTERROTATING STELLAR DISKS AROUND SgrA*

We present a high-resolution simulation of an idealized model to explain the origin of the two young, counterrotating, sub-parsec scale stellar disks around the supermassive black hole SgrA* at the center of the Milky Way. In our model, the collision of a single molecular cloud with a circumnuclear gas disk (similar to the one observed presently) leads to multiple streams of gas flowing toward the black hole and creating accretion disks with angular momentum depending on the ratio of cloud and circumnuclear disk material. The infalling gas creates two inclined, counterrotating sub-parsec scale accretion disks around the supermassive black hole with the first disk forming roughly 1 Myr earlier, allowing it to fragment into stars and get dispersed before the second counterrotating disk forms. Fragmentation of the second disk would lead to the two inclined, counterrotating stellar disks which are observed at the Galactic center. A similar event might be happening again right now at the Milky Way Galactic center. Our model predicts that the collision event generates spiral-like filaments of gas, feeding the Galactic center prior to disk formation with a geometry and inflow pattern that is in agreement with the structure of the so-called mini spiral that has been detected in the Galactic center.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

NUMERICAL SIMULATIONS OF CORONAL MASS EJECTION ON 2011 MARCH 7: ONE-TEMPERATURE AND TWO-TEMPERATURE MODEL COMPARISON

During Carrington rotation (CR) 2107, a fast coronal mass ejection (CME; >2000 km s{sup -1}) occurred in active region NOAA 11164. This event is also associated with a solar energetic particle event. In this study, we present simulations of this CME with one-temperature (1T) and two-temperature (2T: coupled thermodynamics of the electron and proton populations) models. Both the 1T and 2T models start from the chromosphere with heat conduction and radiative cooling. The background solar wind is driven by Alfven-wave pressure and heated by Alfven-wave dissipation in which we have incorporated the balanced turbulence at the top of the closed field lines. The magnetic field of the inner boundary is set up using a synoptic map from Solar Dynamics Observatory/Helioseismic and Magnetic Imager. The Titov-Demoulin flux-rope model is used to initiate the CME event. We compare the propagation of fast CMEs and the thermodynamics of CME-driven shocks in both the 1T and 2T CME simulations. Also, the synthesized white light images are compared with the Solar and Heliospheric Observatory/Large Angle and Spectrometric Coronagraph observations. Because there is no distinction between electron and proton temperatures, heat conduction in the 1T model creates an unphysical temperature precursor in front of the CME-driven shock and makes the shock parameters (e.g., shock Mach number, compression ratio) incorrect. Our results demonstrate the importance of the electron heat conduction in conjunction with proton shock heating in order to produce the physically correct CME structures and CME-driven shocks.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

STELLAR DYNAMOS AND CYCLES FROM NUMERICAL SIMULATIONS OF CONVECTION

We present a series of kinematic axisymmetric mean-field αΩ dynamo models applicable to solar-type stars, for 20 distinct combinations of rotation rates and luminosities. The internal differential rotation and kinetic helicity profiles required to calculate source terms in these dynamo models are extracted from a corresponding series of global three-dimensional hydrodynamical simulations of solar/stellar convection, so that the resulting dynamo models end up involving only one free parameter, namely, the turbulent magnetic diffusivity in the convecting layers. Even though the αΩ dynamo solutions exhibit a broad range of morphologies, and sometimes even double cycles, these models manage to reproduce relatively well the observationally inferred relationship between cycle period and rotation rate. On the other hand, they fail in capturing the observed increase of magnetic activity levels with rotation rate. This failure is due to our use of a simple algebraic α-quenching formula as the sole amplitude-limiting nonlinearity. This suggests that α-quenching is not the primary mechanism setting the amplitude of stellar magnetic cycles, with magnetic reaction on large-scale flows emerging as the more likely candidate. This inference is coherent with analyses of various recent global magnetohydrodynamical simulations of solar/stellar convection.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Two-dimensional numerical simulations of supercritical accretion flows revisited

We study the dynamics of super-Eddington accretion flows by performing two-dimensional radiation-hydrodynamic simulations. Compared with previous works, in this paper we include the T {sub θφ} component of the viscous stress and consider various values of the viscous parameter α. We find that when T {sub θφ} is included, the rotational speed of the high-latitude flow decreases, while the density increases and decreases at the high and low latitudes, respectively. We calculate the radial profiles of inflow and outflow rates. We find that the inflow rate decreases inward, following a power law form of M-dot {sub in}∝r{sup s}. The value of s depends on the magnitude of α and is within the range of ∼0.4-1.0. Correspondingly, the radial profile of density becomes flatter compared with the case of a constant M-dot (r). We find that the density profile can be described by ρ(r)∝r {sup –p} and the value of p is almost same for a wide range of α ranging from α = 0.1 to 0.005. The inward decrease of inflow accretion rate is very similar to hot accretion flows, which is attributed to the mass loss in outflows. To study the origin of outflow, we analyze the convective stability of the slim disk. We find that depending on the value of α, the flow is marginally stable (when α is small) or unstable (when α is large). This is different from the case of hydrodynamical hot accretion flow, where radiation is dynamically unimportant and the flow is always convectively unstable. We speculate that the reason for the difference is because radiation can stabilize convection. The origin of outflow is thus likely because of the joint function of convection and radiation, but further investigation is required.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Numerical simulation of hot accretion flows. III. Revisiting wind properties using the trajectory approach

Previous MHD simulations have shown that wind must exist in black hole hot accretion flows. In this paper, we continue our study by investigating the detailed properties of wind and the mechanism of wind production. For this aim, we make use of a 3D general relativistic MHD simulation of hot accretion flows around a Schwarzschild black hole. To distinguish real wind from turbulent outflows, we track the trajectories of the virtual Lagrangian particles from simulation data. We find two types of real outflows, i.e., a jet and a wind. The mass flux of wind is very significant, and its radial profile can be described by M-dot {sub wind}≈ M-dot {sub BH}(r/20 r{sub s}), with M-dot {sub BH} being the mass accretion rate at the black hole horizon and r{sub s} being the Schwarzschild radius. The poloidal wind speed almost remains constant once they are produced, but the flux-weighted wind speed roughly follows v{sub p,wind}(r)≈0.25v{sub k}(r), with v{sub k}(r) being the Keplerian speed at radius r. The mass flux of the jet is much lower, but the speed is much higher, v{sub p,jet} ∼ (0.3–0.4)c. Consequently, both the energy and momentum fluxes of the wind are much larger than those of the jet. The wind is produced and accelerated primarily by the combination of centrifugal force and magnetic pressure gradient, while the jet is mainly accelerated by the magnetic pressure gradient. Finally, we find that the wind production efficiency ϵ{sub wind}≡ E-dot {sub wind}/ M-dot {sub BH}c{sup 2}∼1/1000 is in good agreement with the value required from large-scale galaxy simulations with active galactic nucleus feedback.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

NUMERICAL SIMULATIONS OF HELICITY CONDENSATION IN THE SOLAR CORONA

The helicity condensation model has been proposed by Antiochos to explain the observed smoothness of coronal loops and the observed buildup of magnetic shear at filament channels. The basic hypothesis of the model is that magnetic reconnection in the corona causes the magnetic stress injected by photospheric motions to collect only at those special locations where prominences are observed to form. In this work we present the first detailed quantitative MHD simulations of the reconnection evolution proposed by the helicity condensation model. We use the well-known ansatz of modeling the closed corona as an initially uniform field between two horizontal photospheric plates. The system is driven by applying photospheric rotational flows that inject magnetic helicity into the corona. The flows are confined to a finite region on the photosphere so as to mimic the finite flux system of a bipolar active region, for example. The calculations demonstrate that, contrary to common belief, opposite helicity twists do not lead to significant reconnection in such a coronal system, whereas twists with the same sense of helicity do produce substantial reconnection. Furthermore, we find that for a given amount of helicity injected into the corona, the evolution of the magnetic shear is insensitive to whether the pattern of driving photospheric motions is fixed or quasi-random. In all cases, the shear propagates via reconnection to the boundary of the flow region while the total magnetic helicity is conserved, as predicted by the model. We discuss the implications of our results for solar observations and for future, more realistic simulations of the helicity condensation process.

79 ASTRONOMY AND ASTROPHYSICS↗