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 1,261 records · Page 70
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2015

A NUMERICAL SIMULATION OF COSMIC-RAY MODULATION NEAR THE HELIOPAUSE

Based on a hybrid galactic cosmic-ray transport model, which incorporated MHD global heliospheric data into Parker’s cosmic-ray transport equation, we studied the behavior of the transport of galactic cosmic rays and the corresponding gradients in their flux near the heliopause (HP). We found that, (1) by increasing the ratio of the parallel diffusion coefficient to the perpendicular diffusion coefficient in the interstellar magnetic field of the outer heliosheath, the simulated radial flux near the HP increases as well. As the ratio multiplying factor reached 10{sup 10}, the radial flux experienced a sudden jump near the HP, similar to what Voyager 1 observed in 2012. (2) The effect of changing the diffusion coefficients’ ratio on the radial flux variation depends on the energy of the cosmic rays, the lower the energy, the more pronounced the effect is. (3) The magnitude of the diffusion coefficients also affect the radial flux near the HP, the modulation beyond the HP varies by adjusting the magnitude multiplying factor.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

NUMERICAL SIMULATIONS OF SUNSPOT DECAY: ON THE PENUMBRA–EVERSHED FLOW–MOAT FLOW CONNECTION

We present a series of high-resolution sunspot simulations that cover a timespan of up to 100 hr. The simulation domain extends about 18 Mm in depth beneath the photosphere and 98 Mm horizontally. We use open boundary conditions that do not maintain the initial field structure against decay driven by convective motions. We consider two setups: a sunspot simulation with penumbra, and a “naked-spot” simulation in which we removed the penumbra after 20 hr through a change in the magnetic top boundary condition. While the sunspot has an Evershed outflow of 3–4 km s{sup −1}, the naked spot is surrounded by an inflow of 1–2 km s{sup −1} in close proximity. However, both spots are surrounded by an outflow on larger scales with a few 100 m s{sup −1} flow speed in the photosphere. While the sunspot has an almost constant magnetic flux content for the simulated timespan of three to four days, the naked spot decays steadily at a rate of 10{sup 21} Mx day{sup −1}. A region with reduced downflow filling factor, which is more extended for the sunspot, surrounds both spots. The absence of downflows perturbs the upflow/downflow mass flux balance and leads to a large-scale radially overturning flow system; the photospheric component of this flow is the observable moat flow. The reduction of the downflow filling factor also inhibits the submergence of magnetic field in the proximity of the spots, which stabilizes them against decay. While this effect is present for both spots, it is more pronounced for the sunspot and explains the almost stationary magnetic flux content.

79 ASTRONOMY AND ASTROPHYSICS↗