Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrodynamics simulations”

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 667 records · Page 37

Detection of the Compressed Primary Stellar Wind in eta Carinae*

A series of three Hubble Space Telescope Space Telescope Imaging Spectrograph (HST/STIS) spectroscopic mappings, spaced approximately one year apart, reveal three partial arcs in [Fe II] and [Ni II] emissions moving outward from η Carinae. We identify these arcs with the shell-like structures, seen in the 3D hydrodynamical simulations, formed by compression of the primary wind by the secondary wind during periastron passages.

spectroscopic mappings↗

The 0.5-2.22 micrometer Scattered Light Spectrum of the Disk around TW Hya: Detection of a Partially Filled Disk Gap at 80 AU*

We present a 0.5-2.2 micrometer scattered light spectrum of the circumstellar disk around TW Hya from a combination of spatially resolved Hubble Space Telescope STIS spectroscopy and NICMOS coronagraphic images of the disk. We investigate the morphology of the disk at distances greater than 40 AU over this wide range of wavelengths, and identify the presence of a depression in surface brightness at approximately 80 AU that could be caused by a gap in the disk. Additionally, we quantify the surface brightness, azimuthal symmetry, and spectral character of the disk as a function of radius. Our analysis shows that the scattering efficiency of the dust is largely neutral to blue over the observed wavelengths. We model the disk as a steady α-disk with an ad hoc gap structure. The thermal properties of the disk are selfconsistently calculated using a three-dimensional radiative transfer code that uses ray tracing to model the heating of the disk interior and scattered light images. We find a good fit to the data over a wide range of distances from the star if we use a model disk with a partially filled gap of 30% depth at 80 AU and with a self-similar truncation knee at 100 AU. The origin of the gap is unclear, but it could arise from a transition in the nature of the disk's dust composition or the presence of a planetary companion. Based on scalings to previous hydrodynamic simulations of gap-opening criteria for embedded proto-planets, we estimate that a planetary companion forming the gap could have a mass between 6 and 28 solar mass.

radiative transfer↗

New Evidence for Efficient Collisionless Heating of Electrons at the Reverse Shock of a Young Supernova Remnant

Although collisionless shocks are ubiquitous in astrophysics, certain key aspects of them are not well understood. In particular, the process known as collisionless electron heating, whereby electrons are rapidly energized at the shock front, is one of the main open issues in shock physics. Here, we present the first clear evidence for efficient collisionless electron heating at the reverse shock of Tycho's supernova remnant (SNR), revealed by Fe K diagnostics using high-quality X-ray data obtained by the Suzaku satellite. We detect K beta (3p yields 1s) fluorescence emission from low-ionization Fe ejecta excited by energetic thermal electrons at the reverse shock front, which peaks at a smaller radius than Fe K alpha (2p yields 1s) emission dominated by a relatively highly ionized component. Comparisons with our hydrodynamical simulations imply instantaneous electron heating to a temperature 1000 times higher than expected from Coulomb collisions alone. The unique environment of the reverse shock, which is propagating with a high Mach number into rarefied ejecta with a low magnetic field strength, puts strong constraints on the physical mechanism responsible for this heating and favors a cross-shock potential created by charge deflection at the shock front. Our sensitive observation also reveals that the reverse shock radius of this SNR is about 10% smaller than the previous measurement using the Fe K alpha morphology from the Chandra observations. Since strong Fe K beta fluorescence is expected only from low-ionization plasma where Fe ions still have many 3p electrons, this feature is key to diagnosing the plasma state and distribution of the immediate postshock ejecta in a young SNR.

Fe K diagnostics↗

Atomic Physics of Shocked Plasma in Winds of Massive Stars

High resolution diffraction grating spectra of X-ray emission from massive stars obtained with Chandra and XMM-Newton have revolutionized our understanding of their powerful, radiation-driven winds. Emission line shapes and line ratios provide diagnostics on a number of key wind parameters. Modeling of resolved emission line velocity profiles allows us to derive independent constraints on stellar mass-loss rates, leading to downward revisions of a factor of a few from previous measurements. Line ratios in He-like ions strongly constrain the spatial distribution of Xray emitting plasma, confirming the expectations of radiation hydrodynamic simulations that X-ray emission begins moderately close to the stellar surface and extends throughout the wind. Some outstanding questions remain, including the possibility of large optical depths in resonance lines, which is hinted at by differences in line shapes of resonance and intercombination lines from the same ion. Resonance scattering leads to nontrivial radiative transfer effects, and modeling it allows us to place constraints on shock size, density, and velocity structure

Physics↗

Late-Time Evolution of Composite Supernova Remnants: Deep Chandra Observations and Hydrodynamical Modeling of a Crushed Pulsar Wind Nebula in SNR G327.1-1.1

In an effort to better understand the evolution of composite supernova remnants (SNRs) and the eventual fate of relativistic particles injected by their pulsars, we present a multifaceted investigation of the interaction between a pulsar wind nebula (PWN) and its host SNR G327.1-1.1. Our 350 ks Chandra X-ray observations of SNR G327.1-1.1 reveal a highly complex morphology; a cometary structure resembling a bow shock, prong-like features extending into large arcs in the SNR interior, and thermal emission from the SNR shell. Spectral analysis of the non-thermal emission offers clues about the origin of the PWN structures, while enhanced abundances in the PWN region provide evidence for mixing of supernova ejecta with PWN material. The overall morphology and spectral properties of the SNR suggest that the PWN has undergone an asymmetric interaction with the SNR reverse shock(RS) that can occur as a result of a density gradient in the ambient medium and or a moving pulsar that displaces the PWN from the center of the remnant. We present hydrodynamical simulations of G327.1-1.1 that show that its morphology and evolution can be described by a approx. 17,000 yr old composite SNR that expanded into a density gradient with an orientation perpendicular to the pulsar's motion. We also show that the RSPWN interaction scenario can reproduce the broadband spectrum of the PWN from radio to gamma-ray wavelengths. The analysis and modeling presented in this work have important implications for our general understanding of the structure and evolution of composite SNRs.

Evolution↗

Parsec-Scale Obscuring Accretion Disk with Large-Scale Magnetic Field in AGNs

A magnetic field dragged from the galactic disk, along with inflowing gas, can provide vertical support to the geometrically and optically thick pc (parsec) -scale torus in AGNs (Active Galactic Nuclei). Using the Soloviev solution initially developed for Tokamaks, we derive an analytical model for a rotating torus that is supported and confined by a magnetic field. We further perform three-dimensional magneto-hydrodynamic simulations of X-ray irradiated, pc-scale, magnetized tori. We follow the time evolution and compare models that adopt initial conditions derived from our analytic model with simulations in which the initial magnetic flux is entirely contained within the gas torus. Numerical simulations demonstrate that the initial conditions based on the analytic solution produce a longer-lived torus that produces obscuration that is generally consistent with observed constraints.

galaxies↗

A Deep X-Ray Look at Abell 2142-Viscosity Constraints from Kelvin-Helmholtz Eddies, a Displaced Cool Peak that Makes a Warm Core, and A Possible Plasma Depletion Layer

cold fronts in detail. We find that the southern cold front exhibits well-developed Kelvin–Helmholtz (KH) eddies seen in the sky plane. Comparing their wavelength and amplitude with those in hydrodynamic simulations of cold fronts in viscous gas, and estimating the gas tangential velocity from centripetal acceleration, we constrain the effective viscosity to be at most 1/5 of Spitzer isotropic viscosity, but consistent with full Braginskii anisotropic viscosity for magnetized plasma. While the northwestern front does not show obvious eddies, its shape and the structure of its brightness profile suggest KH eddies seen in projection. The southern cold front continues in a spiral to the center of the cluster, ending with another cold front only 12 kpc from the gas density peak. The cool peak itself is displaced ∼30 kpc from the brightest cluster galaxy (BCG) (the biggest such offset among centrally peaked clusters), while the X-ray emission on a larger scale is still centered on the BCG, indicating that the BCG is at the center of the gravitational potential and the cool gas is sloshing in it. The specific entropy index of the gas in the peak (K≈49 keV sq.cm) makes A2142 a rare “warm core”; apparently the large displacement of the cool peak by sloshing is the reason. Finally, we find a subtle narrow, straight channel with a 10% drop in X-ray brightness, aligned with the southern cold front—possibly a plasma depletion layer in projection.

Wang, Qian H. S.↗

3D Radiative Hydrodynamics Modeling of Convection of Stars to Probe Their Interiors and Photospheric Properties

The dramatic flow of data from the Kepler and K2 missions opens the opportunity to significantly improve our knowledge of stellar interiors, surface dynamics, and structure. However, interpretation of these observations is a challenging task because it depends on tiny effects that can be studied only with advanced first-principles modeling. We present results of 3D time-dependent radiative hydrodynamic simulations of stellar outer convection zones and atmospheres taking into account chemical composition, radiative transfer, turbulence effects, and a realistic equation of state for main sequence stars. We will discuss properties of convective structure and dynamics, convective overshoot, effects of magnetic fields and rotation, as well as the potential influence of turbulent surface dynamics on high-precision RV measurements.

Photospheric↗

Synthetic Observations of Helioseismic and Magnetic Imager

Synthetic observations of the FeI line (6173A), which corresponds to the operational wavelength of Helioseismic and Magnetic Imager (HMI) on board of the NASA Solar Dynamics Observatory (SDO). Synthetic observations represent 4D maps for four angular distances from the solar disk center (0deg), 30deg, 45deg, and 60deg at 30degrees latitude in range +/- 0.75 Angstrom from the reference line of 6173.341 Angstrom, number of points across the spectral line is 256, spatial resolution of 100km, and temporal resolution of seconds. Synthetic data set obtained using SPINOR code (Frutiger et al. 2000). Synthetic observations computed from hydrodynamic simulations, which cover 25Mm of upper layers of the convection zone and 1Mm of atmosphere, the horizontal size of the domain of 80x80 Mm. Simulation results obtained the “StellarBox” code (Wray et al., 2015, 2018) with imposed rotation rate at 30deg latitude using F-plane approximation.

Irina Kitiashvili↗

Tidal Disruptions of Main-sequence Stars. IV. Relativistic Effects and Dependence on Black Hole Mass

Using a suite of fully relativistic hydrodynamic simulations applied to main-sequence stars with realistic internal density profiles, we examine full and partial tidal disruptions across a wide range of black hole mass (10(exp 5)≤M(BH)/M(ʘ)≤5 x 10(exp 7) and stellar mass (0.3≤M(*)/M(ʘ)≤3) as larger M(BH) leads to stronger relativistic effects. For fixed M( ), as M(BH) increases, the ratio of the maximum pericenter distance yielding full disruptions ( ) to its Newtonian prediction rises rapidly, becoming triple the Newtonian value for M(BH)=5 x 10(exp 7)M(ʘ), while the ratio of the energy width of the stellar debris for full disruptions to the Newtonian prediction decreases steeply, resulting in a factor of 2 correction at M(BH)=5 x 10(exp 7)M(ʘ). We provide approximate formulae that express the relativistic corrections of both R(t) and the energy width relative to their Newtonian approximate estimates. For partial disruptions, we find that the fractional remnant mass for a given ratio of the pericenter to R(t) is higher for larger M(BH). These results have several implications. As M(BH) increases above ~10(exp 7)M(ʘ), the cross section for complete disruptions is suppressed by competition with direct capture. However, the cross-section ratio for partial to complete disruptions depends only weakly on M(BH). The relativistic correction to the debris energy width delays the time of peak mass-return rate and diminishes the magnitude of the peak return rate. For M(BH)≳10(exp 7)M(ʘ), the M(BH)-dependence of the full disruption cross section and the peak mass-return rate and time is influenced more by relativistic effects than by Newtonian dynamics.

Taeho Ryu↗

On the Maximum Stellar Rotation to Form a Black Hole without an Accompanying Luminous Transient

The collapse of a massive star with low angular momentum content is commonly argued to result in the formation of a black hole without an accompanying bright transient. Our goal in this Letter is to understand the flow in and around a newly-formed black hole, involving accretion and rotation, via general relativistic hydrodynamics simulations aimed at studying the conditions under which infalling material can accrete without forming a centrifugally supported structure and, as a result, generate no effective feedback. If the feedback from the black hole is, on the other hand, significant, the collapse would be halted and we suggest that the event is likely to be followed by a bright transient. We find that feedback is only efficient if the specific angular momentum of the infalling material at the innermost stable circular orbit exceeds that of geodesic circular flow at that radius by at least≈20%. We use the results of our simulations to constrain the maximal stellar rotation rates of the disappearing massive progenitors PHL293B-LBV and N6946-BH1, and to provide an estimate of the overall rate of disappearing massive stars. We find that about a few percent of single O-type stars with measured rotational velocities are expected to spin below the critical value before collapse and are thus predicted to vanish without a trace.

Scott C Noble↗

Tidal Disruptions of Main Sequence Stars - I. Observable Quantities and their Dependence on Stellar and Black Hole Mass

This paper introduces a series of papers presenting a quantitative theory for the tidal disruption of main-sequence stars by supermassive black holes. Using fully general relativistic hydrodynamics simulations and MESA-model initial conditions, we explore the pericenter-dependence of tidal disruption properties for eight stellar masses(MM0.1510)and six black hole masses( MM105 105BH7). We present here the results most relevant to observations. The effects of internal stellar structure and relativity decouple for both the disruption cross section and the characteristic energy width of the debris. Moreover, the full disruption cross section is almost independent of Må for Må/Me3. Independent of Må, relativistic effects increase the critical pericenter distance for full disruption events by up to a factor of∼3 relative to the Newtonian prediction. The probability of a direct capture is also independent of Må; at MBH/Me;5×106 this probability is equal to the probability of a complete disruption. The breadth of the debris energy distribution ΔE can differ from the standard estimate by factors of0.35−2, depending on Må and MBH, implying a corresponding change(∝(ΔE)−3/2)in the characteristic mass-return timescale. We provide analytic forms, suitable for use in both event rate estimates and parameter inference, to describe all these trends. For partial disruptions, we find a nearly universal relation between the star’s angular momentum and the fraction of Må remaining. Within the “empty loss-cone” regime, partial disruptions must precede full disruptions. These partial disruptions can drastically affect the rate and appearance of subsequent total disruptions.

Taeho Ryu↗

Tidal Disruptions of Main-sequence Stars. IV. Relativistic Effects and Dependence on Black Hole Mass

Using a suite of fully relativistic hydrodynamic simulations applied to main-sequence stars with realistic internal density profiles, we examine full and partial tidal disruptions across a wide range of black hole mass( MM105 105BH7)and stellar mass (MM0.33) as larger MBH leads to stronger relativistic effects. For fixed Må, as MBH increases, the ratio of the maximum pericenter distance yielding full disruptions (t) to its Newtonian prediction rises rapidly, becoming triple the Newtonian value for= ́MM510BH7, while the ratio of the energy width of the stellar debris for full disruptions to the Newtonian prediction decreases steeply, resulting in a factor of 2 correction at= ́MM510BH7. We provide approximate formulae that express the relativistic corrections of both t and the energy width relative to their Newtonian approximate estimates. For partial disruptions, we find that the fractional remnant mass for a given ratio of the pericenter to t is higher for larger MBH. These results have several implications. As MBH increases above~M107, the cross section for complete disruptions is suppressed by competition with direct capture. However, the cross-section ratio for partial to complete disruptions depends only weakly on MBH. The relativistic correction to the debris energy width delays the time of peak mass-return rate and diminishes the magnitude of the peak return rate. ForMM10BH7, the MBH-dependence of the full disruption cross section and the peak mass-return rate and time is influenced more by relativistic effects than by Newtonian dynamics.

Taeho Ryu↗

Radiography in High Mass X-Ray Binaries: Micro-Structure of the Stellar Wind Through Variability of the Column Density

Context. In high mass X-ray binaries, an accreting compact object orbits a high mass star, which loses mass through a dense and in homogeneous wind. Aims. Using the compact object as an X-ray backlight, the time variability of the absorbing column density in the wind can be exploited in order to shed light on the micro-structure of the wind and obtain unbiased stellar mass-loss rates for high mass stars. Methods. We developed a simplified representation of the stellar wind where all the matter is gathered in spherical “clumps” that are radially advected away from the star. This model enables us to explore the connections between the stochastic properties of the wind and the variability of the column density for a comprehensive set of parameters related to the orbit and to the wind micro-structure, such as the size of the clumps and their individual mass. In particular, we focus on the evolution with the orbital phase of the standard deviation of the column density and of the characteristic duration of enhanced absorption episodes. Using the porosity length, we derive analytical predictions and compare them to the standard deviations and coherence time scales that were obtained. Results. We identified the favorable systems and orbital phases to determine the wind micro-structure. The coherence time scale of the column density is shown to be the self-crossing time of a single clump in front of the compact object. We thus provide a procedure to get accurate measurements of the size and of the mass of the clumps, purely based on the observable time variability of the column density. Conclusions. The coherence time scale grants direct access to the size of the clumps, while their mass can be deduced separately from the amplitude of the variability. We further show how monitoring the variability at superior conjunctions can probe the onset of the clump-forming region above the stellar photosphere. If the high column density variations in some high mass X-ray binaries are due to unaccreted clumps which are passing by the line-of-sight, this would require high mass clumps to reproduce the observed peak-to-peak amplitude and coherence time scales. These clump properties are marginally compatible with the ones derived from radiative-hydrodynamics simulations. Alternatively, the following components could contribute to the variability of the column density: larger orbital scale structures produced by a mechanism that has yet to be identified or a dense environment in the immediate vicinity of the accretor, such as an accretion disk, an outflow, or a spherical shell surrounding the magnetosphere of the accreting neutron star.

stars: mass loss↗

Deciphering the nature of the pulsar wind nebula CTB 87 with XMM-Newton

CTB 87 (G74.9+1.2) is an evolved supernova remnant (SNR) which hosts a peculiar pulsar wind nebula (PWN). The X-ray peak is offset from that observed in radio and lies towards the edge of the radio nebula. The putative pulsar, CXOU J201609.2+371110, was first resolved with Chandra and is surrounded by a compact and a more extended X-ray nebula. Here, we use a deep XMM–Newton observation to examine the morphology and evolutionary stage of the PWN and to search for thermal emission expected from a supernova shell or reverse shock interaction with supernova ejecta. We do not find evidence of thermal X-ray emission from the SNR and place an upper limit on the electron density of 0.05 per cu.cm for a plasma temperature kT ∼ 0.8 keV. The morphology and spectral properties are consistent with a ∼20-kyr-old relic PWN expanding into a stellar wind-blown bubble. We also present the first X-ray spectral index map from the PWN and show that we can reproduce its morphology by means of 2D axisymmetric relativistic hydrodynamical simulations.

B Guest↗

Atacama Cosmology Telescope: Combined kinematic and thermal Sunyaev-Zel’dovich measurements from BOSS CMASS and LOWZ halos

The scattering of cosmic microwave background (CMB) photons off the free-electron gas in galaxies and clusters leaves detectable imprints on high resolution CMB maps: the thermal and kinematic Sunyaev-Zel’dovich effects (tSZ and kSZ respectively). We use combined microwave maps from the Atacama Cosmology Telescope DR5 and Planck in combination with the CMASS (mean redshifthzi¼0.55and host halo masshMviri¼3×1013M⊙) and LOWZ (hzi¼0.31,hMviri¼5×1013M⊙) galaxy catalogs from the Baryon Oscillation Spectroscopic Survey (BOSS DR10 and DR12), to study the gas associated with these galaxy groups. Using individual reconstructed velocities, we perform a stacking analysis and reject the no-kSZ hypothes is at 6.5σ, the highest significance to date. This directly translates into a measurement of the electron number density profile, and thus of the gas density profile. Despite the limited signal to noise, the measurement shows at high significance that the gas density profile is more extended than the dark matter density profile, for any reasonable baryon abundance (formally>90σfor the cosmic baryon abundance). We simultaneously measure the tSZ signal, i.e., the electron thermal pressure profile of the same CMASS objects, and reject theno-tSZ hypothesis at10σ. We combine tSZ and kSZ measurements to estimate the electron temperature to20% precision in several aperture bins, and find it comparable to the virial temperature. In a companion paper, we analyze these measurements to constrain the gas thermodynamics and the properties of feedback inside galaxy groups. We present the corresponding LOWZ measurements in this paper, ruling out a null kSZ (tSZ)signal at 2.9ð13.9Þσ, and leave their interpretation to future work. This paper and the companion paper demonstrate that current CMB experiments can detect and resolve gas profiles in low mass halos and at high redshifts, which are the most sensitive to feedback in galaxy formation and the most difficult to measure any other way. They will be a crucial input to cosmological hydrodynamical simulations, thus improving our understanding of galaxy formation. These precise gas profiles are already sufficient to reduce the main limiting theoretical systematic in galaxy-galaxy lensing: baryonic uncertainties. Future such measurements will thus unleash the statistical power of weak lensing from the Rubin, Euclid and Roman observatories. Our stacking software Thumb Stack is publicly available and directly applicable to future Simons Observatory andCMB-S4 data.

Emmanuel Schaan↗

Modeling the Hydrodynamics, Sediment Transport, and Valley Incision of Outlet-forming Floods from Martian Crater Lakes

Open-basin lakes that have outlets, formed where the stored water breached the basin-confining topography (e.g., crater rims), have been recognized on Mars for decades. However, the mechanics involved in formation of these outlets has not been studied in detail. Here, we study this outlet formation process through a series of numerical experiments accomplished by adding erosion and sediment transport routines to an existing hydrodynamic simulation, ANUGA. These experiments are consistent with outlets being carved in single flood events, at least if sediment was available throughout outlet formation as in a transport-limited regime. Peak discharges from these open-basin lakes, including from Jezero crater, were 10^4–10^6 m^3/s, and the outlet-forming floods lasted for days to weeks, depending on lake size. Moreover, we find that sediment transport by suspension was likely important during these floods – outlets could not have reached their observed sizes if flow was limited to bedload transport alone. This in turn may explain why small open-basin lakes are less commonly recognized on Mars, as small lakes had less inherent capacity to transport sediment in suspension and thus incise confining topography catastrophically. Geometry, topographic setting, and earlier fluvial valley development may have played a significant role in influencing the magnitude of outlet-forming floods and resulting outlet canyon incision.

Mars↗

The Asymmetric Inner Disk of the Herbig Ae Star HD 163296 in the Eyes of VLTI/MATISSE: Evidence for a Vortex?

Context.A complex environment exists in the inner few astronomical units of planet-forming disks. High-angular-resolution observa-tions play a key role in our understanding of the disk structure and the dynamical processes at work.Aims.In this study we aim to characterize the mid-infrared brightness distribution of the inner disk of the young intermediate-massstar HD 163296 from early VLTI/MATISSE observations taken in theL- andN-bands. We put special emphasis on the detection ofpotential disk asymmetries.Methods.We use simple geometric models to fit the interferometric visibilities and closure phases. Our models include a smoothedring, a flat disk with an inner cavity, and a 2D Gaussian. The models can account for disk inclination and for azimuthal asymmetriesas well. We also perform numerical hydrodynamical simulations of the inner edge of the disk.Results.Our modeling reveals a significant brightness asymmetry in theL-band disk emission. The brightness maximum of the asym-metry is located at the NW part of the disk image, nearly at the position angle of the semimajor axis. The surface brightness ratio inthe azimuthal variation is3.5±0.2. Comparing our result on the location of the asymmetry with other interferometric measurements,we confirm that the morphology of ther<0.3au disk region is time-variable. We propose that this asymmetric structure, located in ornear the inner rim of the dusty disk, orbits the star. To find the physical origin of the asymmetry, we tested a hypothesis where a vortexis created by Rossby wave instability, and we find that a unique large-scale vortex may be compatible with our data. The half-lightradius of theL-band-emitting region is0.33±0.01au, the inclination is52◦+5◦−7◦, and the position angle is143◦±3◦. Our models predictthat a non-negligible fraction of theL-band disk emission originates inside the dust sublimation radius forμm-sized grains. Refractorygrains or large (&10μm-sized) grains could be the origin of this emission.N-band observations may also support a lack of smallsilicate grains in the innermost disk (r.0.6au), in agreement with our findings fromL-band data.

J Varga↗