Search NASA⌕ Search

SEARCH · Search NASA

Results for “STRATIFICATION”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Insights into the broadband emission of the TeV blazar Mrk 501 during the first X-ray polarization measurements

Aims.We present the first multiwavelength study of Mrk 501 that contains simultaneous very-high-energy (VHE)γ-ray observations and X-ray polarization measurements from the Imaging X-ray Polarimetry Explorer (IXPE). Methods.We used radio-to-VHE data from a multiwavelength campaign carried out between March 1, 2022, and July 19, 2022 (MJD 59639 to MJD 59779). The observations were performed by MAGIC,Fermi-LAT,NuSTAR,Swift(XRT and UVOT), and several other instruments that cover the optical and radio bands to complement the IXPE pointings. We characterized the dynamics of the broadband emission around the X-ray polarization measurements through its multiband fractional variability and correlations, and compared changes observed in the polarization degree to changes seen in the broadband emission using a multi-zone leptonic scenario. Results.During the IXPE pointings, the VHE state is close to the average behavior, with a 0.2–1 TeV flux of 20%–50% of the emission of the Crab Nebula. Additionally, it shows low variability and a hint of correlation between VHEγ-rays and X-rays. Despite the average VHE activity, an extreme X-ray behavior is measured for the first two IXPE pointings, taken in March 2022 (MJD 59646 to 59648 and MJD 59665 to 59667), with a synchrotron peak frequency > 1 keV. For the third IXPE pointing, in July 2022 (MJD 59769 to 59772), the synchrotron peak shifts toward lower energies and the optical/X-ray polarization degrees drop. All three IXPE epochs show an atypically low Compton dominance in theγ-rays. The X-ray polarization is systematically higher than at lower energies, suggesting an energy stratification of the jet. While during the IXPE epochs the polarization angles in the X-ray, optical, and radio bands align well, we find a clear discrepancy in the optical and radio polarization angles in the middle of the campaign. Such results further support the hypothesis of an energy-stratified jet. We modeled broadband spectra taken simultaneous to the IXPE pointings, assuming a compact zone that dominates in the X-rays and the VHE band, and an extended zone stretching farther downstream in the jet that dominates the emission at lower energies.NuSTARdata allow us to precisely constrain the synchrotron peak and therefore the underlying electron distribution. The change between the different states observed in the three IXPE pointings can be explained by a change in the magnetization and/or the emission region size, which directly connects the shift in the synchrotron peak to lower energies with the drop in the polarization degree.

Astronomy & Astrophysics↗
DOE OSTI2024

Performance of high-order Godunov-type methods in simulations of astrophysical low Mach number flows

High-order Godunov methods for gas dynamics have become a standard tool for simulating different classes of astrophysical flows. Their accuracy is mostly determined by the spatial interpolant used to reconstruct the pair of Riemann states at cell interfaces and by the Riemann solver that computes the interface fluxes. In most Godunov-type methods, these two steps can be treated independently, so that many different schemes can in principle be built from the same numerical framework. Because astrophysical simulations often test out the limits of what is feasible with the computational resources available, it is essential to find the scheme that produces the numerical solution with the desired accuracy at the lowest computational cost. However, establishing the best combination of numerical options in a Godunov-type method to be used for simulating a complex hydrodynamic problem is a nontrivial task. In fact, formally more accurate schemes do not always outperform simpler and more diffusive methods, especially if sharp gradients are present in the flow. For this work, we used our fully compressible Seven-League Hydro (SLH) code to test the accuracy of six reconstruction methods and three approximate Riemann solvers on two- and three-dimensional (2D and 3D) problems involving subsonic flows only. We considered Mach numbers in the range from 10 −3 to 10 −1 , which are characteristic of many stellar and geophysical flows. In particular, we considered a well-posed, 2D, Kelvin–Helmholtz instability problem and a 3D turbulent convection zone that excites internal gravity waves in an overlying stable layer. Although the different combinations of numerical methods converge to the same solution with increasing grid resolution for most of the quantities analyzed here, we find that (i) there is a spread of almost four orders of magnitude in computational cost per fixed accuracy between the methods tested in this study, with the most performant method being a combination of a low-dissipation Riemann solver and a sextic reconstruction scheme; (ii) the low-dissipation solver always outperforms conventional Riemann solvers on a fixed grid when the reconstruction scheme is kept the same; (iii) in simulations of turbulent flows, increasing the order of spatial reconstruction reduces the characteristic dissipation length scale achieved on a given grid even if the overall scheme is only second order accurate; (iv) reconstruction methods based on slope-limiting techniques tend to generate artificial, high-frequency acoustic waves during the evolution of the flow; and (v) unlimited reconstruction methods introduce oscillations in the thermal stratification near the convective boundary, where the entropy gradient is steep.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

The stability of stratified spatially periodic shear flows at low Péclet number

This work addresses the question of the stability of stratified, spatially periodic shear flows at low Péclet number but high Reynolds number. This little-studied limit is motivated by astrophysical systems, where the Prandtl number is often very small. Furthermore, it can be studied using a reduced set of “low-Péclet-number equations” proposed by Lignières [“The small-Péclet-number approximation in stellar radiative zones,” Astron. Astrophys. 348, 933–939 (1999)]. Through a linear stability analysis, we first determine the conditions for instability to infinitesimal perturbations. We formally extend Squire’s theorem to the low-Péclet-number equations, which shows that the first unstable mode is always two-dimensional. We then perform an energy stability analysis of the low-Péclet-number equations and prove that for a given value of the Reynolds number, above a critical strength of the stratification, any smooth periodic shear flow is stable to perturbations of arbitrary amplitude. In that parameter regime, the flow can only be laminar and turbulent mixing does not take place. Finding that the conditions for linear and energy stability are different, we thus identify a region in parameter space where finite-amplitude instabilities could exist. Using direct numerical simulations, we indeed find that the system is subject to such finite-amplitude instabilities. We determine numerically how far into the linearly stable region of parameter space turbulence can be sustained.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

A whole-path importance-sampling scheme for Feynman path integral calculations of absolute partition functions and free energies

Using Feynman path integrals, a molecular partition function can be written as a double integral with the inner integral involving all closed paths centered at a given molecular configuration, and the outer integral involving all possible molecular configurations. In previous work employing Monte Carlo methods to evaluate such partition functions, we presented schemes for importance sampling and stratification in the molecular configurations that constitute the path centroids, but we relied on free-particle paths for sampling the path integrals. At low temperatures, the path sampling is expensive because the paths can travel far from the centroid configuration. In this paper, we now present a scheme for importance sampling of whole Feynman paths based on harmonic information from an instantaneous normal mode calculation at the centroid configuration, which we refer to as harmonically guided whole-path importance sampling (WPIS). We obtain paths conforming to our chosen importance function by rejection sampling from a distribution of free-particle paths. Sample calculations on CH 4 demonstrate that at a temperature of 200 K, about 99.9% of the free-particle paths can be rejected without integration, and at 300 K, about 98% can be rejected. Finally, we also show that it is typically possible to reduce the overhead associated with the WPIS scheme by sampling the paths using a significantly lower-order path discretization than that which is needed to converge the partition function.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2016

Effect of simple solutes on the long range dipolar correlations in liquid water

Intermolecular correlations in liquid water at ambient conditions have generally been characterized through short range density fluctuations described through the atomic pair distribution functions. Recent numerical and experimental results have suggested that such a description of order or structure in liquid water is incomplete and there exist considerably longer ranged orientational correlations in water that can be studied through dipolar correlations. In this study, using large scale classical, atomistic molecular dynamics simulations using TIP4P-Ew and TIP3P models of water, we show that salts such as sodium chloride (NaCl), potassium chloride (KCl), caesium chloride (CsCl), and magnesium chloride (MgCl{sub 2}) have a long range effect on the dipolar correlations, which cannot be explained by the notion of structure making and breaking by dissolved ions. Observed effects are explained through orientational stratification of water molecules around ions and their long range coupling to the global hydrogen bond network by virtue of the sum rule for water. The observations for single hydrophilic solutes are contrasted with the same for a single methane (CH{sub 4}) molecule. We observe that even a single small hydrophobe can result in enhancement of long range orientational correlations in liquid water, contrary to the case of dissolved ions, which have been observed to have a reducing effect. The observations from this study are discussed in the context of hydrophobic effect.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

Analysis of dynamics, stability, and flow fields' structure of an accelerated hydrodynamic discontinuity with interfacial mass flux by a general matrix method

We develop a general matrix method to analyze from a far field the dynamics of an accelerated interface between incompressible ideal fluids of different densities with interfacial mass flux and with negligible density variations and stratification. We rigorously solve the linearized boundary value problem for the dynamics conserving mass, momentum, and energy in the bulk and at the interface. We find a new hydrodynamic instability that develops only when the acceleration magnitude exceeds a threshold. This critical threshold value depends on the magnitudes of the steady velocities of the fluids, the ratio of their densities, and the wavelength of the initial perturbation. The flow has potential velocity fields in the fluid bulk and is shear-free at the interface. The interface stability is set by the interplay of inertia and gravity. For weak acceleration, inertial effects dominate, and the flow fields experience stable oscillations. For strong acceleration, gravity effects dominate, and the dynamics is unstable. For strong accelerations, this new hydrodynamic instability grows faster than accelerated Landau-Darrieus and Rayleigh-Taylor instabilities. For given values of the fluids' densities and their steady bulk velocities, and for a given magnitude of acceleration, we find the critical and maximum values of the initial perturbation wavelength at which this new instability can be stabilized and at which its growth is the fastest. The quantitative, qualitative, and formal properties of the accelerated conservative dynamics depart from those of accelerated Landau-Darrieus and Rayleigh-Taylor dynamics. New diagnostic benchmarks are identified for experiments and simulations of unstable interfaces.

Flow dynamics↗
DOE OSTI2019

Laser-ion acceleration using mixed compositions: Tailoring the target for each species

Particle-in-cell simulations of laser- ion acceleration demonstrate marked discrepancies in the acceleration experienced by the different ion species in complex target compositions, especially when the target becomes relativistically transparent to the pulse during irradiation. Beginning with proton contaminants in a carbon target, we show how the higher charge-to-mass ratio of the protons results in species stratification and late-time suppression of the carbon acceleration. The target normal sheath acceleration (TNSA) primarily experienced by the protons can be exploited to mitigate this tamping by using a shaped rear surface of the target, leaving the break-out afterburner-driven carbons to accelerate close to the laser axis and then experience less tamping during a late-time TNSA phase. Here, we then explore preferentially accelerating the lighter species in a mixed composition target, particularly focusing on deuteron beam applications. We examine three different target compositions with the same areal electron density, CD 2 , CH 2 , and 7 LiD, and propose an alternative shaping of the rear surface of 7 LiD to increase the number of high-energy deuterons in the beam.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Adjustments to the law of the wall above an Amazon forest explained by a spectral link

Modification to the law of the wall represented by a dimensionless correction function ϕ RSL (z/h) is derived using atmospheric turbulence measurements collected at two sites in the Amazon in near-neutral stratification, where z is the distance from the forest floor and h is the mean canopy height. The sites are the Amazon Tall Tower Observatory for z/h∈ [1,2.3] and the Green Ocean Amazon (GoAmazon) site for z/h∈ [1,1.4]. Here, a link between the vertical velocity spectrum E ww (k) (k is the longitudinal wavenumber) and ϕ RSL is then established using a co-spectral budget (CSB) model interpreted by the moving-equilibrium hypothesis. The key finding is that ϕ RSL is determined by the ratio of two turbulent viscosities and is given as ν t,BL /ν t,RSL, where ν t,RSL = (1/A)∫$^{∞}_{0}$ τ(k)E ww (k)dk, ν t,BL = k v (z−d)u * , τ(k) is a scale-dependent decorrelation time scale between velocity components, A = C R /(1−C I ) = 4.5 is predicted from the Rotta constant C R = 1.8, and the isotropization of production constant C I = 3/5 given by rapid distortion theory, k v is the von Kármán constant, u * is the friction velocity at the canopy top, and d is the zero-plane displacement. Because the transfer of energy across scales is conserved in E ww (k) and is determined by the turbulent kinetic energy dissipation rate (ε), the CSB model also predicts that ϕ RSL scales with L BL /L d , where L BL is the length scale of attached eddies to z = d, and L d = u$_{*}^{3}$/ε is a macro-scale dissipation length.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Compressible Rayleigh–Taylor instability subject to isochoric initial background state

The effect of compressibility on the single-mode Rayleigh–Taylor instability is examined using two (2D) and three-dimensional (3D) direct numerical simulations. To isolate compressibility from background stratification effects, this work employs a constant density profile on each side of the interface. The numerical simulations are performed at various Reynolds numbers using the gas kinetic method for static Mach numbers up to M = 0.4. The most important finding is that compressibility acting in isolation enhances the instability and perturbations grows faster with increasing Mach number, unlike previous results with background isothermal state, which show suppression of the instability at higher static Mach numbers. In addition, compressibility is also shown to increase the bubble-spike asymmetry. While the instability grows faster for the 3D case, the findings are qualitatively similar in 2D and 3D. Furthermore, the dynamical reasons underlying the effect of compressibility are elucidated by examining the evolution of vorticity and turbulent kinetic energy transport equations.

42 ENGINEERING↗
DOE OSTI2024

Hierarchical Ion Interactions in the Direct Air Capture of CO 2 at Air/Aqueous Interfaces

The direct air capture (DAC) of CO 2 using aqueous solvents is plagued by slow kinetics and interfacial barriers that limit effectiveness in combating climate change. Functionalizing air/aqueous surfaces with charged amphiphiles shows promise in accelerating DAC; however, insight into these interfaces and how they evolve in time remains poorly understood. Specifically, competitive ion interactions between DAC reagents and reaction products feedback onto the interfacial structure, thereby modulating interfacial chemical composition and overall function. In this work, we probe the role of glycine amino acid anions (Gly - ), an effective CO 2 capture reagent, that promotes the organization of cationic oligomers at air/aqueous interfaces. These surfaces are probed with vibrational sum frequency generation spectroscopy and molecular dynamics simulations. Our findings demonstrate that the competition for surface sites between Gly - and captured carbonaceous anions (HCO 3 - , CO 3 2 -, carbamates) drives changes in surface hydration, which in turn tunes oligomer ordering. This phenomenon is related to a hierarchical ordering of anions at the surface that are electrostatically attracted to the surface and their ability to compete for interfacial water. These results point to new ways to tune interfaces for DAC via stratification of ions based on relative surface propensities and specific ion effects.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Climate-driven oscillation of phosphorus and iron limitation in the North Pacific Subtropical Gyre

The supply of nutrients is a fundamental regulator of ocean productivity and carbon sequestration. Nutrient sources, sinks, residence times, and elemental ratios vary over broad scales, including those resulting from climate-driven changes in upper water column stratification, advection, and the deposition of atmospheric dust. These changes can alter the proximate elemental control of ecosystem productivity with cascading ecological effects and impacts on carbon sequestration. Here, we report multidecadal observations revealing that the ecosystem in the eastern region of the North Pacific Subtropical Gyre (NPSG) oscillates on subdecadal scales between inorganic phosphorus (P i ) sufficiency and limitation, when P i concentration in surface waters decreases below 50–60 nmol⋅kg −1 . In situ observations and model simulations suggest that sea-level pressure changes over the northwest Pacific may induce basin-scale variations in the atmospheric transport and deposition of Asian dust-associated iron (Fe), causing the eastern portion of the NPSG ecosystem to shift between states of Fe and P i limitation. Our results highlight the critical need to include both atmospheric and ocean circulation variability when modeling the response of open ocean pelagic ecosystems under future climate change scenarios.

Pacific Decadal Oscillation↗
DOE OSTI2022

Microbial functional diversity across biogeochemical provinces in the central Pacific Ocean

Enzymes catalyze key reactions within Earth’s life-sustaining biogeochemical cycles. Here, we use metaproteomics to examine the enzymatic capabilities of the microbial community (0.2 to 3 µm) along a 5,000-km-long, 1-km-deep transect in the central Pacific Ocean. Eighty-five percent of total protein abundance was of bacterial origin, with Archaea contributing 1.6%. Over 2,000 functional KEGG Ontology (KO) groups were identified, yet only 25 KO groups contributed over half of the protein abundance, simultaneously indicating abundant key functions and a long tail of diverse functions. Vertical attenuation of individual proteins displayed stratification of nutrient transport, carbon utilization, and environmental stress. The microbial community also varied along horizontal scales, shaped by environmental features specific to the oligotrophic North Pacific Subtropical Gyre, the oxygen-depleted Eastern Tropical North Pacific, and nutrient-rich equatorial upwelling. Some of the most abundant proteins were associated with nitrification and C1 metabolisms, with observed interactions between these pathways. The oxidoreductases nitrite oxidoreductase (NxrAB), nitrite reductase (NirK), ammonia monooxygenase (AmoABC), manganese oxidase (MnxG), formate dehydrogenase (FdoGH and FDH), and carbon monoxide dehydrogenase (CoxLM) displayed distributions indicative of biogeochemical status such as oxidative or nutritional stress, with the potential to be more sensitive than chemical sensors. Enzymes that mediate transformations of atmospheric gases like CO, CO 2 , NO, methanethiol, and methylamines were most abundant in the upwelling region. Here, we identified hot spots of biochemical transformation in the central Pacific Ocean, highlighted previously understudied metabolic pathways in the environment, and provided rich empirical data for biogeochemical models critical for forecasting ecosystem response to climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

High-Fidelity Velocity and Concentration Measurements of Turbulent Buoyant Jets

Accurate models of turbulent buoyant flows are essential for the design of the cooling circuit of nuclear reactors and passive safety systems. However, available models fail to fully capture the physics of turbulent mixing when buoyancy becomes predominant with respect to momentum. Therefore, high-fidelity experiments of well-controlled fundamental flows are needed to develop and validate more accurate models. We analyze experiments of positive and negative turbulent buoyant jets, both in uniform and stratified environments, with the aim of understanding the thermal hydraulics of turbulent mixing with variable density and providing high-fidelity data for the development and validation of turbulence models. Non-intrusive, simultaneous particle image velocimetry and laser-induced fluorescence measurements were carried out to acquire instantaneous velocity and concentration fields on a vertical section parallel to the axis of a jet in the self-similar region. The refractive index matching method was applied to measure high-resolution buoyant jets with up to 8.6% density difference. These data are free of the typical errors that characterize optical measurements of buoyancy-driven flows (e.g. natural and mixed convection) where the refractive index of the fluid is inhomogeneous throughout the measurement domain. Turbulent statistics and entrainment of buoyant jets in uniform and stratified environments are presented. These data are compared with non-buoyant jets in a uniform environment, as a reference to investigate the effects of buoyancy and stratification on turbulent mixing. The results will be used for the assessment of current turbulence models and as a basis for the development of a new model that captures turbulent mixing.

laser induced fluorescence↗
DOE OSTI2026

Mathematical Models and Numerical Methods for High-Fidelity Simulation of Ignition of Reactive Mixtures by Nanosecond Plasma Discharges in Realistic Configurations

We present a newly developed framework for the numerical simulation of ignition of reactive mixtures using single or repeated nanosecond discharge pulses. The framework builds upon the AMReX library, using the existing compressible solver PeleC and low-Mach solver PeleLMeX and allowing for adaptive mesh refinement, complex geometries, and execution on next-generation high-performance computing (HPC) systems. High-fidelity elementary models are adopted for weakly-ionised plasma discharges with significant energy deposition, consistent with nanosecond discharge pulses, and then implemented in the solver. The treatment of non-thermal electrons and charged species, thermodynamics of non-equilbrium species, plasma kinetics, limiting time scales, and boundary conditions for charged species are discussed and addressed for computational efficiency. The framework is demonstrated for three relevant applications: single and multi-pulse discharges in air, single pulse ignition of an ethylene/air mixture, and a three-dimensional plasma discharge in air with temperature stratification. The successful application of the framework demonstrates the feasibility of high-fidelity simulation of ignition of air/hydrocarbon mixtures in three-dimensions with multiple discharge pulses.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Turbulence kinetic energy budget and conditional sampling of momentum, scalar, and intermittency fluxes in thermally stratified wind farms

A description of the turbulence kinetic energy budget is investigated to understand the dynamics and turbulence energy transfer between the atmospheric boundary layer and large wind farms. The turbulence kinetic energy budget terms are a dominant factor in explaining the energy supply process in wind farms. Large eddy simulations are used to generate the thermally stratified wind turbine array boundary layer, where stable and unstable conditions are considered. The wakes represent a sink of mean kinetic energy, and turbulence transport redistributes energy in the wake region. The energy balance between the budget terms is dependent on the physical location in the wind farm. During the unstable scenario, turbulence is driven by the production, transport, and buoyancy terms. Further, the contribution of the turbulence kinetic energy budget terms is reduced with increasing atmospheric stability, especially above the shear layer, where the mechanical production is the dominant term. Pressure diffusion redistributes the turbulent energy from the streamwise direction into the other components of the normal Reynolds stresses. This distinctive feature of the changing thermal stratification is presented through the quadrant analysis. The momentum ejection and sweep frequencies are both comparable and dominant. The scalar flux shows an increase in the contribution of the negative correlation in the quadrants. The quadrant analysis is extended to account for intermittency based on the pointwise Holder exponents. As a result, positive correlations are dominant in the wake and in higher layers of the domain for both considered cases.

17 WIND ENERGY↗
DOE OSTI2017

Reconsidering barriers to wind power projects: community engagement, developer transparency and place

In 2016, we undertook a nationally representative wind power perceptions survey of individuals living within 8 km of over 600 projects in the United States, generating 1705 telephone, web, and mail responses. We sought information on a variety of topics, including procedural fairness and its relationship to project attitude, the foci of the present analysis. Here, we present a series of descriptive statistics and regression results, emphasizing those residents who were aware of their local project prior to construction. Sample weighting is employed to account for stratification and non-response. We find that a developer being open and transparent, a community being able to influence the outcome, and having a say in the planning process are all statistically significant predictors of a process perceived as being ‘fair,’ with an open and transparent developer having the largest effect. We also find developer transparency and ability to influence outcomes to have statistically significant relationships to a more positive attitude, with those findings holding when aesthetics, landscape, and wind turbine sound considerations are controlled for. The results indicate that jurisdictions might consider developing procedures, which ensure citizens are consulted and heard, and benchmarks or best practices for developer interaction with communities and citizens.

17 WIND ENERGY↗
DOE OSTI2015

Microlensing constraints on broad absorption and emission line flows in the quasar H1413+117

We present new integral field spectroscopy of the gravitationally lensed broad absorption line (BAL) quasar H1413+117, covering the ultraviolet restframe spectral range. We observe strong microlensing signatures in lensed image D, and we use this microlensing to simultaneously constrain both the broad emission and broad absorption line gas. The wavelength independence of image D magnifications across the broad emission lines (BELs) indicates a lower limit on the broad emission line region (BELR) size equal to the Einstein radius (ER) of the system: ≳11 (〈M〉/M{sub ⊙}){sup 0.5} lt-day for a lens redshift of 1.4 and ≳15 (〈M〉/M{sub ⊙}){sup 0.5} lt-day for z{sub L} = 0.94. Lensing simulations verify that the observed wavelength independence is very unlikely for BELRs with significant velocity stratification at size scales below an ER. We perform spectral decomposition to derive the intrinsic BEL and continuum spectrum, subject to BAL absorption. We reconstruct the intrinsic BAL absorption profile, whose features allow us to constrain outflow kinematics in the context of a disk-wind model. We find a very sharp, blueshifted onset of absorption of 1500 km s{sup −1} in both C iv and N v, which may correspond to an inner edge of a disk-wind’s radial outflow. The lower ionization Si iv and Al iii have higher-velocity absorption onsets, consistent with a decreasing ionization parameter with radius in an accelerating outflow. There is evidence of strong absorption in the BEL component, which indicates a high covering factor for absorption over two orders of magnitude in outflow radius.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Solving the excitation and chemical abundances in shocks: The case of HH 1

We present deep spectroscopic (3600–24700 Å ) X-shooter observations of the bright Herbig–Haro object HH 1, one of the best laboratories to study the chemical and physical modifications caused by protostellar shocks on the natal cloud. We observe atomic fine structure lines, H i and He i recombination lines and H{sub 2} ro-vibrational lines (more than 500 detections in total). Line emission was analyzed by means of Non-local Thermal Equilibiurm codes to derive the electron temperature and density, and for the first time we are able to accurately probe different physical regimes behind a dissociative shock. We find a temperature stratification in the range 4000 K ÷ 80,000 K, and a significant correlation between temperature and ionization energy. Two density regimes are identified for the ionized gas, a more tenuous, spatially broad component (density ∼10{sup 3} cm{sup −3}), and a more compact component (density ⩾ 10{sup 5} cm{sup −3}) likely associated with the hottest gas. A further neutral component is also evidenced, having a temperature ≲10,000 K and a density >10{sup 4} cm{sup −3}. The gas fractional ionization was estimated by solving the ionization equilibrium equations of atoms detected in different ionization stages. We find that neutral and fully ionized regions co-exist inside the shock. Also, indications in favor of at least partially dissociative shock as the main mechanism for molecular excitation are derived. Chemical abundances are estimated for the majority of the detected species. On average, abundances of non-refractory/refractory elements are lower than solar of about 0.15/0.5 dex. This indicates the presence of dust inside the medium, with a depletion factor of iron of ∼40%.

79 ASTRONOMY AND ASTROPHYSICS↗