Search NASA⌕ Search

SEARCH · Search NASA

Results for “3D 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,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

SYNTHETIC OBSERVATIONS OF MAGNETIC FIELDS IN PROTOSTELLAR CORES

The role of magnetic fields in the early stages of star formation is not well constrained. In order to discriminate between different star formation models, we analyze 3D magnetohydrodynamic simulations of low-mass cores and explore the correlation between magnetic field orientation and outflow orientation over time. We produce synthetic observations of dust polarization at resolutions comparable to millimeter-wave dust polarization maps observed by the Combined Array for Research in Millimeter-wave Astronomy and compare these with 2D visualizations of projected magnetic field and column density. Cumulative distribution functions of the projected angle between the magnetic field and outflow show different degrees of alignment in simulations with differing mass-to-flux ratios. The distribution function for the less magnetized core agrees with observations finding random alignment between outflow and field orientations, while the more magnetized core exhibits stronger alignment. We find that fractional polarization increases when the system is viewed such that the magnetic field is close to the plane of the sky, and the values of fractional polarization are consistent with observational measurements. The simulation outflow, which reflects the underlying angular momentum of the accreted gas, changes direction significantly over over the first ∼0.1 Myr of evolution. This movement could lead to the observed random alignment between outflows and the magnetic fields in protostellar cores.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Observational Signatures of Transverse Magnetohydrodynamic Waves and Associated Dynamic Instabilities in Coronal Flux Tubes

Magnetohydrodynamic (MHD) waves permeate the solar atmosphere and constitute potential coronal heating agents. Yet, the waves detected so far may be but a small subset of the true existing wave power. Detection is limited by instrumental constraints but also by wave processes that localize the wave power in undetectable spatial scales. In this study, we conduct 3D MHD simulations and forward modeling of standing transverse MHD waves in coronal loops with uniform and non-uniform temperature variation in the perpendicular cross-section. The observed signatures are largely dominated by the combination of the Kelvin–Helmholtz instability (KHI), resonant absorption, and phase mixing. In the presence of a cross-loop temperature gradient, we find that emission lines sensitive to the loop core catch different signatures compared to those that are more sensitive to the loop boundary and the surrounding corona, leading to an out-of-phase intensity and Doppler velocity modulation produced by KHI mixing. In all of the considered models, common signatures include an intensity and loop width modulation at half the kink period, a fine strand-like structure, a characteristic arrow-shaped structure in the Doppler maps, and overall line broadening in time but particularly at the loop edges. For our model, most of these features can be captured with a spatial resolution of 0.″33 and a spectral resolution of 25 km s{sup −1}, although we do obtain severe over-estimation of the line width. Resonant absorption leads to a significant decrease of the observed kinetic energy from Doppler motions over time, which is not recovered by a corresponding increase in the line width from phase mixing and KHI motions. We estimate this hidden wave energy to be a factor of 5–10 of the observed value.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Vortex and Sink Flows in Eruptive Flares as a Model for Coronal Implosions

Eruptive flares are sudden releases of magnetic energy that involve many phenomena, several of which can be explained by the standard 2D flare model and its realizations in 3D. We analyze a 3D magnetohydrodynamics simulation, in the framework of this model, that naturally explains the contraction of coronal loops in the proximity of the flare sites, as well as the inflow toward the region above the cusp-shaped loops. We find that two vorticity arcs located along the flanks of the erupting magnetic flux rope are generated as soon as the eruption begins. The magnetic arcades above the flux rope legs are then subjected to expansion, rotation, or contraction depending on which part of the vortex flow advects them. In addition to the vortices, an inward-directed magnetic pressure gradient exists in the current sheet below the magnetic flux rope. It results in the formation of a sink that is maintained by reconnection. We conclude that coronal loop apparent implosions observed during eruptive flares are the result of hydromagnetic effects related to the generation of vortex and sink flows when a flux rope moves in a magnetized environment.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Observational Signatures of a Kink-unstable Coronal Flux Rope Using Hinode /EIS

The signatures of energy release and energy transport for a kink-unstable coronal flux rope are investigated via forward modeling. Synthetic intensity and Doppler maps are generated from a 3D numerical simulation. The CHIANTI database is used to compute intensities for three Hinode /EIS emission lines that cover the thermal range of the loop. The intensities and Doppler velocities at simulation-resolution are spatially degraded to the Hinode /EIS pixel size (1″), convolved using a Gaussian point-spread function (3″), and exposed for a characteristic time of 50 s. The synthetic images generated for rasters (moving slit) and sit-and-stare (stationary slit) are analyzed to find the signatures of the twisted flux and the associated instability. We find that there are several qualities of a kink-unstable coronal flux rope that can be detected observationally using Hinode /EIS, namely the growth of the loop radius, the increase in intensity toward the radial edge of the loop, and the Doppler velocity following an internal twisted magnetic field line. However, EIS cannot resolve the small, transient features present in the simulation, such as sites of small-scale reconnection (e.g., nanoflares).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Synchrotron Intensity Gradients as Tracers of Interstellar Magnetic Fields

On the basis of the modern understanding of MHD turbulence, we propose a new way of using synchrotron radiation: using synchrotron intensity gradients (SIGs) for tracing astrophysical magnetic fields. We successfully test the new technique using synthetic data obtained with 3D MHD simulations and provide the demonstration of the practical utility of the technique by comparing the directions of magnetic fields that are obtained with PLANCK synchrotron intensity data to the directions obtained with PLANCK synchrotron polarization data. We demonstrate that the SIGs can reliably trace magnetic fields in the presence of noise and can provide detailed maps of magnetic field directions. We also show that the SIGs are relatively robust for tracing magnetic fields while the low spatial frequencies of the synchrotron image are removed. This makes the SIGs applicable to the tracing of magnetic fields using interferometric data with single-dish measurement absent. We discuss the synergy of using the SIGs together with synchrotron polarization in order to find the actual direction of the magnetic fields and quantify the effects of Faraday rotation as well as with other ways of studying astrophysical magnetic fields. We test our method in the presence of noise and the resolution effects. We stress the complementary nature of the studies using the SIG technique and those employing the recently introduced velocity gradient techniques that trace magnetic fields using spectroscopic data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Numerical Study on Outflows in Seyfert Galaxies I: Narrow Line Region Outflows in NGC 4151

The origin of narrow line region (NLR) outflows remains unknown. In this paper, we explore the scenario in which these outflows are circumnuclear clouds driven by energetic accretion disk winds. We choose the well-studied nearby Seyfert galaxy NGC 4151 as an example. By performing 3D hydrodynamical simulations, we are able to reproduce the radial distributions of velocity, mass outflow rate, and kinetic luminosity of NLR outflows in the inner 100 pc deduced from spatial resolved spectroscopic observations. The demanded kinetic luminosity of disk winds is about two orders of magnitude higher than that inferred from the NLR outflows, but is close to the ultrafast outflows (UFO) detected in the X-ray spectrum and a few times lower than the bolometric luminosity of the Seyfert. Our simulations imply that the scenario is viable for NGC 4151. The existence of the underlying disk winds can be confirmed by their impacts on higher density ISM, e.g., shock excitation signs, and the pressure in NLR.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Scaling Relations of Starburst-driven Galactic Winds

Using synthetic absorption lines generated from 3D hydrodynamical simulations, we explore how the velocity of a starburst-driven galactic wind correlates with the star formation rate (SFR) and SFR density. We find strong correlations for neutral and low ionized gas, but no correlation for highly ionized gas. The correlations for neutral and low ionized gas only hold for SFRs below a critical limit set by the mass loading of the starburst, above which point the scaling relations flatten abruptly. Below this point the scaling relations depend on the temperature regime being probed by the absorption line, not on the mass loading. The exact scaling relation depends on whether the maximum or mean velocity of the absorption line is used. We find that the outflow velocity of neutral gas can be up to five times lower than the average velocity of ionized gas, with the velocity difference increasing for higher ionization states. Furthermore, the velocity difference depends on both the SFR and mass loading of the starburst. Thus, absorption lines of neutral or low ionized gas cannot easily be used as a proxy for the outflow velocity of the hot gas.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Magnetohydrodynamical Effects on Nuclear Deflagration Fronts in Type Ia Supernovae

This article introduces a report of the effects of magnetic fields on non-distributed nuclear burning fronts as a possible solution to a fundamental problem for the thermonuclear explosion of a Chandrasekhar mass (${M}_{\mathrm{Ch}}$) white dwarf (WD), the currently favored scenario for the majority of Type Ia SNe. All existing 3D hydrodynamical simulations predict strong global mixing of the burning products due to Rayleigh–Taylor (RT) instabilities, which contradicts observations. As a first step toward studying the flame physics, we present a set of computational magnet-hydrodynamic models in rectangular flux tubes, resembling a small inner region of a WD. We consider initial magnetic fields up to ${10}^{12}\,{\rm{G}}$ of various orientations. We find an increasing suppression of RT instabilities starting at about ${10}^{9}\,{\rm{G}}$. The front speed tends to decrease with increasing magnitude up to about ${10}^{11}\,{\rm{G}}$. For even higher fields new small-scale, finger-like structures develop, which increase the burning speed by a factor of 3 to 4 above the field-free RT-dominated regime. We imply that the new instability may provide sufficiently accelerated energy production during the distributed burning regime to go over the Chapman–Jougey limit and trigger a detonation. Furthermore, we discuss the possible origins of high magnetic fields during the final stage of the progenitor evolution or the explosion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Identifying Anticyclonic Vortex Features Produced by the Rossby Wave Instability in Protoplanetary Disks

Several nearby protoplanetary disks have been observed to display large-scale crescents in the (sub)millimeter dust continuum emission. One interpretation is that these structures correspond to anticyclonic vortices generated by the Rossby wave instability within the gaseous disk. Such vortices have local gas overdensities and are expected to concentrate dust particles with a Stokes number around unity. This process might catalyze the formation of planetesimals. Whereas recent observations showed that dust crescents are indeed regions where millimeter-size particles have abnormally high concentration relative to the gas and smaller grains, no observations have yet shown that the gas within the crescent region counterrotates with respect to the protoplanetary disk. Here we investigate the detectability of anticyclonic features through measurement of the line-of-sight component of the gas velocity obtained with ALMA. We carry out 2D hydrodynamic simulations and 3D radiative transfer calculations of a protoplanetary disk characterized by a vortex created by the tidal interaction with a massive planet. As a case study, the disk parameters are chosen to mimic the IRS 48 system, which has the most prominent crescent observed to date. We generate synthetic ALMA observations of both the dust continuum and 12 CO emission around the frequency of 345 GHz. We find that the anticyclonic features of the vortex are weak but can be detected if both the source and the observational setup are properly chosen. We provide a recipe for maximizing the probability of detecting such vortex features and present an analysis procedure to infer their kinematic properties.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Magnetically Induced Current Piston for Generating Extreme-ultraviolet Fronts in the Solar Corona

Single-pulse, globally propagating coronal fronts, called Extreme-ultraviolet (EUV) waves, were first confirmed in 1995 by the Extreme-ultraviolet Imaging Telescope and every observed EUV wave since has been associated with a coronal mass ejection (CME). The physical mechanism underlying these waves has been debated for two decades with wave or pseudo-wave theories being advocated. We introduce a hybrid model where EUV waves are compressional fronts driven by a reverse electric current layer induced by the time-dependent CME core current. The reverse current layer flows in a direction opposite to the CME core current and is an eddy current layer necessary to maintain magnetic flux conservation above the layer. Repelled by the core current, the reverse current layer accelerates upward so it acts as a piston that drives a compressional perturbation in the coronal regions above. Given a sufficiently fast piston speed, the compressional perturbation becomes a shock that separates from the piston when the piston slows down. Since the model relates the motion of the EUV front to CME properties, the model provides a bound for the core current of an erupting CME. The model is supported and motivated by detailed results from both laboratory experiments and ideal 3D magnetohydrodynamic simulations. Overlaps and differences with other models and spacecraft observations are reported on.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Pressure–Strain Interaction as the Energy Dissipation Estimate in Collisionless Plasma

The dissipative mechanism in weakly collisional plasma is a topic that pervades decades of studies without a consensus solution. We compare several energy dissipation estimates based on energy transfer processes in plasma turbulence and provide justification for the pressure–strain interaction as a direct estimate of the energy dissipation rate. The global and scale-by-scale energy balances are examined in 2.5D and 3D kinetic simulations. We show that the global internal energy increase and the temperature enhancement of each species are directly tracked by the pressure–strain interaction. The incompressive part of the pressure–strain interaction dominates over its compressive part in all simulations considered. The scale-by-scale energy balance is quantified by scale filtered Vlasov–Maxwell equations, a kinetic plasma approach, and the lag dependent von Kármán–Howarth equation, an approach based on fluid models. We find that the energy balance is exactly satisfied across all scales, but the lack of a well-defined inertial range influences the distribution of the energy budget among different terms in the inertial range. Therefore, the widespread use of the Yaglom relation in estimating the dissipation rate is questionable in some cases, especially when the scale separation in the system is not clearly defined. In contrast, the pressure–strain interaction balances exactly the dissipation rate at kinetic scales regardless of the scale separation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Nature and Scalings of Density Fluctuations of Compressible Magnetohydrodynamic Turbulence with Applications to the Solar Wind

Abstract The solar wind is a magnetized and turbulent plasma. Its turbulence is often dominated by Alfvénic fluctuations and often deemed as nearly incompressible far away from the Sun, as shown by in situ measurements near 1 au. However, for solar wind closer to the Sun, the plasma β decreases (often lower than unity) while the turbulent Mach number M t increases (can approach unity, e.g., transonic fluctuations). These conditions could produce significantly more compressible effects, characterized by enhanced density fluctuations, as seen by several space missions. In this paper, a series of 3D MHD simulations of turbulence are carried out to understand the properties of compressible turbulence, particularly the generation of density fluctuations. We find that, over a broad range of parameter space in plasma β , cross helicity, and polytropic index, the turbulent density fluctuations scale linearly as a function of M t , with the scaling coefficients showing weak dependence on parameters. Furthermore, through detailed spatiotemporal analysis, we show that the density fluctuations are dominated by low-frequency nonlinear structures, rather than compressible MHD eigenwaves. These results could be important for understanding how compressible turbulence contributes to solar wind heating near the Sun.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Formation and Reconnection of Electron Scale Current Layers in the Turbulent Outflows of a Primary Reconnection Site

We simulate with 3D particle in cell, the spontaneous formation of turbulent outflows in an initially laminar 3D reconnecting current layer. We observe the formation of many secondary current layers and reconnection sites in the outflow. The approach we follow is to study each individual feature within the turbulent outflow. To identify all clusters of current in the outflow we use a clustering technique widely used in unsupervised machine learning: density-based spatial clustering of applications with noise. Once the clusters are identified we measure their size and compute reconnection indicators to establish which are undergoing reconnection. With this analysis we establish that the size of the current clusters reaches all the way from its initial system scale down to subelectron skin depth scale. We observe that the smaller current clusters are more prone to reconnecting and to releasing energy. We then find the process of reconnection of the smaller current cluster to be of the recently observed electron-only type that leaves the ions essentially unaffected.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Afterglows from Binary Neutron Star Postmerger Systems Embedded in Active Galactic Nuclei Disks

Abstract The observability of afterglows from binary neutron star mergers occurring within active galactic nuclei (AGN) disks is investigated. We perform 3D GRMHD simulations of a postmerger system and follow the jet launched from the compact object. We use semianalytic techniques to study the propagation of the blast wave powered by the jet through an AGN disk-like external environment, extending to distances beyond the disk scale height. The synchrotron emission produced by the jet-driven forward shock is calculated to obtain the afterglow emission. The observability of this emission at different frequencies is assessed by comparing it to the quiescent AGN emission. In the scenarios where the afterglow could temporarily outshine the AGN, we find that detection will be more feasible at higher frequencies (≳10 14 Hz) and the electromagnetic counterpart could manifest as a fast variability in the AGN emission, on timescales less than a day.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Evolution of Magnetohydrodynamic Turbulence in the Expanding Solar Wind: Residual Energy and Intermittency

We conduct 3D magnetohydrodynamic simulations of decaying turbulence in the context of the solar wind. To account for the spherical expansion of the solar wind, we implement the expanding box model. The initial turbulence comprises uncorrelated counterpropagating Alfvén waves and exhibits an isotropic power spectrum. Our findings reveal the consistent generation of negative residual energy whenever nonlinear interactions are present, independent of the normalized cross helicity σc and compressibility. The spherical expansion facilitates this process. The resulting residual energy is primarily distributed in the perpendicular direction, with S 2 (b) − S 2 (u) ∝ l ⊥ or equivalently –E r ∝ k$^{–2}_{⟂}$. Here S 2 (b) and S 2 (u) are second-order structure functions of magnetic field and velocity respectively. In most runs, S 2 (b) develops a scaling relation S 2 (b) ∝ $l$$^{1/2}_{⊥}$ (E b ∝ k$^{–3/2}_{⟂}$). In contrast, S 2 (u) is consistently shallower than S 2 (b), which aligns with in situ observations of the solar wind. We observe that the higher-order statistics of the turbulence, which act as a proxy for intermittency, depend on the initial σc and are strongly affected by the expansion effect. Generally, the intermittency is more pronounced when the expansion effect is present. Finally, we find that in our simulations, although the negative residual energy and intermittency grow simultaneously as the turbulence evolves, the causal relation between them seems to be weak, possibly because they are generated on different scales.

Interplanetary turbulence↗
DOE OSTI2016

MODELING OBSERVED DECAY-LESS OSCILLATIONS AS RESONANTLY ENHANCED KELVIN–HELMHOLTZ VORTICES FROM TRANSVERSE MHD WAVES AND THEIR SEISMOLOGICAL APPLICATION

In the highly structured solar corona, resonant absorption is an unavoidable mechanism of energy transfer from global transverse MHD waves to local azimuthal Alfvén waves. Due to its localized nature, direct detection of this mechanism is extremely difficult. Yet, it is the leading theory explaining the observed fast damping of the global transverse waves. However, at odds with this theoretical prediction are recent observations that indicate that in the low-amplitude regime such transverse MHD waves can also appear decay-less, a still unsolved phenomenon. Recent numerical work has shown that Kelvin–Helmholtz instabilities (KHI) often accompany transverse MHD waves. In this work, we combine 3D MHD simulations and forward modeling to show that for currently achieved spatial resolution and observed small amplitudes, an apparent decay-less oscillation is obtained. This effect results from the combination of periodic brightenings produced by the KHI and the coherent motion of the KHI vortices amplified by resonant absorption. Such an effect is especially clear in emission lines forming at temperatures that capture the boundary dynamics rather than the core, and reflects the low damping character of the local azimuthal Alfvén waves resonantly coupled to the kink mode. Due to phase mixing, the detected period can vary depending on the emission line, with those sensitive to the boundary having shorter periods than those sensitive to the loop core. This allows us to estimate the density contrast at the boundary.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Kinematic Evidence for an Embedded Planet in the IM Lupi Disk

Abstract We test the hypothesis that an embedded giant planet in the IM Lupi protostellar disk can produce velocity kinks seen in CO line observations as well as the spiral arms seen in scattered light and continuum emission. We inject planets into 3D hydrodynamics simulations of IM Lupi, generating synthetic observations using Monte Carlo radiative transfer. We find that an embedded planet of 2–3 M Jup can reproduce non-Keplerian velocity perturbations, or “kinks”, in the 12 CO J = 2–1 channel maps. Such a planet can also explain the spiral arms seen in 1.25 mm dust continuum emission and 1.6 μ m scattered-light images. We show that the wake of the planet can be traced in the observed peak velocity map, which appears to closely follow the morphology expected from our simulations and from analytic models of planet–disk interaction.

79 ASTRONOMY AND ASTROPHYSICS↗