Search NASA⌕ Search

SEARCH · Search NASA

Results for “density distribution”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Evaluating Transfer and Pumping of Slurries from Pulsed Jet Mixed Vessels

Due to gravity, solids in slurries will settle if density differences between the solids and liquid are positive (i.e., particle has a negative buoyant force) unless rheological properties and flow conditions are adequate to overcome the gravitational effects. The rate of settling depends on the force balance of the particle, which includes the surface forces associated with fluid rheology. Given the same fluid and solid properties, particles that are larger and denser tend to settle faster. When pumping slurry into a vessel at concentrations precluding hindered settling with insufficient mixing, particle and density distributions can result in preferential settling, creating stratification in the solids concentration within the vessel. For vessels with transfer line inlets located in the lower portion of the tank, the stratified solids concentration may be detrimental to the transfer system performance. Elevated concentrations of solids in the slurry entrained at the inlet to the transfer line can result in the effective viscosity or slurry bulk density exceeding the design limits of the pump. These conditions could result in plugging of the transfer line or onset of cavitation of the pumps because of excessive pressure drop. These conditions can be exacerbated with periodic inlet conditions existing at the transfer line inlet. Periodic conditions can result when vessel mixing is intermittent such as with pulsed jet mixers (PJM). The transfer line inlet conditions are impacted by the periodic nature of the PJM operations with respect to suspension of solids and their transport to the inlet of the transfer line. A scaling approach is presented, and corresponding test requirements are developed for assessing the prevention of plugging the pipeline. Line plugging mechanisms are addressed that do not include plugging due to steady-state high-density slurry entering the transfer line and reducing the net positive suction head available (NPSHA) at the pump inlet to below that required for pump operation. Items considered include the transition to reduced relative flow velocities, such that the critical pipe velocity for solids deposition, Ucd, is not maintained, and segregation of heavy solids during the transport. The recommended requirements to prevent plugging include: • Limits for viscosity and density for entrained slurry to prevent the pressure drop in the pipeline from exceeding pump capacity. • Limits for viscosity and density for entrained slurry to prevent the net positive suction head available (NPSHA) from falling below the net positive suction head required (NPSHR) for operating the pump. • Transfer line velocity and flow rate requirements to maintain solids in suspension, while avoiding line plugging that results from deposition of solids within the transfer line. This paper describes the development of the scaling and testing requirements to verify that proposed approaches for transfer and pump out are appropriately developed

periodic flow, mixing, pipeline plugging, scaling,↗
DOE OSTI2023

A numerical solver for investigating the space charge effect on the electric field in liquid argon time projection chambers

This paper reports the development of a numerical solver aimed to simulate the interaction between the space charge (i.e. ions) distribution and the electric field in liquid argon time projection chamber (LArTPC) detectors. The ion transport equation is solved by a time-accurate, cell-centered finite volume method and the electric potential equation by a continuous finite element method. The electric potential equation updates the electric field which provides the drift velocity to the ion transport equation. The ion transport equation updates the space charge density distribution which appears as the source term in the electric potential equation. The interaction between the space charge distribution and the electric field is numerically simulated within each physical time step. The convective velocity in the ion transport equation can include the background flow velocity in addition to the electric drift velocity. The numerical solver has been parallelized using the Message Passing Interface (MPI) library. Numerical tests show and verify the capability and accuracy of the current numerical solver. It is planned that the developed numerical solver, together with a Computational Fluid Dynamics (CFD) package which provides the flow velocity field, can be used to investigate the space charge effect on the electric field in large-scale particle detectors.

42 ENGINEERING↗
DOE OSTI2021

Understanding Solid Electrolyte Interphase Nucleation and Growth on Lithium Metal Surfaces

Experiments and theory are needed to decode the exact structure and distribution of components of a passivation layer formed at the anode surface of Li metal batteries, known as the Solid Electrolyte Interphase (SEI). Due to the inherent dynamic behavior as well as the lithium reactivity, the SEI structure and its growth mechanisms are still unclear. This study uses molecular simulation and computational chemistry tools to investigate the initial nucleation and growth dynamics of LiOH and Li2O that provide us with thermodynamics and structural information about the nucleating clusters of each species. Following the most favorable pathways for the addition of each of the components to a given nascent SEI cluster reveals their preferential nucleation mechanisms and illustrates different degrees of crystallinity and electron density distribution that are useful to understand ionic transport through SEI blocks.

25 ENERGY STORAGE↗
DOE OSTI2022

High-pressure structural systematics of dysprosium metal compressed in a neon pressure medium to 182 GPa

Here we present an experimental and theoretical study of dysprosium metal compressed in the soft pressure transmitting medium Ne up to 182 and 300 GPa, respectively. Angle-dispersive x-ray powder diffraction data from each of the high-pressure polymorphs shows anisotropic compression behavior indicating changes to the electron density distribution throughout its polymorphic landscape. We compare the monoclinic (mC4) and orthorhombic (oF16) structures for the collapsed structure for Dy above 82 GPa and verify that the oF16 structure offers a better fit to our data than the previously reported mC4 structure. Further, we have found that the oF16 structure undergoes similar anisotropic compression of its lattice parameters, with a turning point above 160 GPa; suggesting a potential phase transition at pressures much higher than achieved in this study. Density functional theory calculations show the likely candidate for this new high-pressure phase is the isosymmetric oF8 structure, which is predicted to be lower in energy than the oF16 structure above 275 GPa.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Improving the probability tables of the cross section of near closed-shell nuclei [Slides]

The level density distribution of near closed-shell nuclei is much lower than the typical nucleus, therefore, the cross sections show significant fluctuations, and these fluctuations are not predictable. The current methodology used to describe such behavior and construct the probability table of the cross section is based on the extrapolation of the average resonance widths and average resonance spacings from the resonance region and use these parameters to construct the probability table. Although this is a standard and widely used technique, it does not take into account the existing experimental data, such for total and the elastic cross section. Our goal is to extend the current theory and provide a more general approach to compute the probability distribution function of the cross section combining the existing probability tables and the available experimental data. Results will be presented for 90 Zr.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2019

Extended range X-ray pair distribution functions

Here we describe a dual detector system for high-energy X-ray, simultaneous, small and wide-angle X-ray scattering (SAXS and WAXS), designed to extract extended-range pair distribution functions (ER-PDF) for disordered materials. The hardware and software provides continuous reciprocal space coverage over atomic to nanometer length-scales. Details of the varying resolution, splicing of data and normalization on an absolute scale are outlined. In addition, the combination of SAXS and WAXS theory is considered with a view to enabling a direct Fourier transformation of the structure factor spanning multiple length-scales into real space. Important distinctions between the ER-PDF and the pair distance distribution function (PDDF) representations are demonstrated. It is shown that when the SAXS intensity in the structure factor, S(Q), is similar to the WAXS intensity, the contributions to the ER-PDF are minimal. However, when the SAXS S(Q) intensities are substantially stronger than the WAXS, the ER-PDF can provide important structural information on the local, intermediate and nanometer length-scales. Notably, the ER-PDF method provides direct information on particle sizes and their density distributions, overcoming the limitations of PDDF analysis for densely packed systems.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Simulations of nozzle gas flow and gas-puff Z-pinch implosions on the Weizmann Z-pinch

We present simulations of an oxygen gas puff Z-pinch on a University scale generator at the Weizmann Institute of Science. The work accounts for the detailed geometry of the nozzle, the initial neutral gas density distribution, and the subsequent implosion. The modeling results show significant improvement with data for the current at the time of stagnation in comparison with a previous effort [Rosenzweig et al., Phys. Plasmas 27, 022705 (2020)]. As a first step, we performed simulations of the flow of neutral diatomic oxygen from a plenum through a nozzle within a recessed cathode, across a gap, and into the anode with a recessed grounded honeycomb. These simulations show an agreement with the measured initial gas density profiles within the region not blocked by the recesses and accessible to visible measurements. The computed neutral gas flow profile serves as the initial condition for a radiation magnetohydrodynamic simulation of the implosion using the MACH2-TCRE code. By considering the specific details of the nozzle and chamber geometry, we find agreement with the measured current profile, including the inductive notch. The simulations predict that the plasma undergoes a strong pinch within the hidden anode recess. The simulations also predict the strongest radiation pulse occurs within the anode recess and at the time of the observed inductive notch.

Physics↗
DOE OSTI2024

A novel Bayesian approach for decomposing the radio emission of quasars: I. Modelling the radio excess in red quasars

ABSTRACT Studies show that both radio jets from the active galactic nuclei (AGNs) and the star formation (SF) activity in quasar host galaxies contribute to the quasar radio emission; yet their relative contributions across the population remain unclear. Here, we present an improved parametric model that allows us to statistically separate the SF and AGN components in observed quasar radio flux density distributions, and investigate how their relative contributions evolve with AGN bolometric luminosity ($L_\mathrm{bol}$) and redshift (z) using a fully Bayesian method. Based on the newest data from LOw-Frequency ARray Two-metre Sky Survey data release 2, our model gives robust fitting results out to $z\sim 4$, showing a quasar host galaxy SF rate (SFR) evolution that increases with bolometric luminosity and with redshift out to $z\sim 4$. This differs from the global cosmic SFR density, perhaps due to the importance of galaxy mergers. The prevalence of radio AGN emissions increases with quasar luminosity, but has little dependence on redshift. Furthermore, our new methodology and large sample size allow us to subdivide our data set to investigate the role of other parameters. Specifically, in this paper, we explore quasar colour and demonstrate that the radio excess in red quasars is due to an enhancement in AGN-related emission, since the host galaxy SF contribution to the total radio emission is independent of quasar colour. We also find evidence that this radio enhancement occurs mostly in quasars with weak or intermediate radio power.

Yue, B. -H (ORCID:0009000989352929)↗
DOE OSTI2021

Heterogeneous Composition of Oxygen-Evolving Complexes in Crystal Structures of Dark-Adapted Photosystem II

Photosystem II (PSII) is a homodimeric protein complex that catalyzes water oxidation at the oxygen-evolving complex (OEC), a heterocubanoid calcium-tetra-manganese cluster. Here, we analyze the omit electron-density peaks of the OEC’s metal ions in five X-ray free-electron laser PSII structures at a resolution between 2.15 Å and 1.95 Å. The omit peaks can be described by the total number of electrons and approximated by the variance of electron density distribution when the distributions are spherically symmetric. We show that the number of electrons of metal centers are different in the two OECs of PSII dimers, implying that either the oxidation states and/or occupancies of individual metal ions are different in the two monomers. In either case, we find that the two OECs of dark-adapted PSII dimers in crystals are not fully synchronized in the same expected S 1 state. Differences in redox states of the OEC in PSII only partially account for the observation that the electron densities integrate to a smaller number of electrons than expected. Differences between the determined and expected relative electron numbers are much larger than the estimated errors, indicating heterogeneity in the OEC composition. Here, the reported findings suggest that two OECs of dark-adapted PSII crystals are typically not fully synchronized in the S 1 state and likely generate heterogeneity in the crystallographical data of PSII intermediate structures.

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

Direct observation of elemental fluctuation and oxygen octahedral distortion-dependent charge distribution in high entropy oxides

The enhanced compositional flexibility to incorporate multiple-principal cations in high entropy oxides (HEOs) offers the opportunity to expand boundaries for accessible compositions and unconventional properties in oxides. Attractive functionalities have been reported in some bulk HEOs, which are attributed to the long-range compositional homogeneity, lattice distortion, and local chemical bonding characteristics in materials. However, the intricate details of local composition fluctuation, metal-oxygen bond distortion and covalency are difficult to visualize experimentally, especially on the atomic scale. Here, we study the atomic structure-chemical bonding-property correlations in a series of perovskite-HEOs utilizing the recently developed four-dimensional scanning transmission electron microscopy techniques which enables to determine the structure, chemical bonding, electric field, and charge density on the atomic scale. The existence of compositional fluctuations along with significant composition-dependent distortion of metal-oxygen bonds is observed. Consequently, distinct variations of metal-oxygen bonding covalency are shown by the real-space charge-density distribution maps with sub-ångström resolution. The observed atomic features not only provide a realistic picture of the local physico-chemistry of chemically complex HEOs but can also be directly correlated to their distinctive magneto-electronic properties.

36 MATERIALS SCIENCE↗
DOE OSTI2023

SCALE 6.2 Lattice Physics Performance Assessment

The US Nuclear Regulatory Commission relies on the lattice physics analysis capabilities of the SCALE code system to perform confirmatory licensing analyses. Either SCALE lattice physics code—TRITON/NEWT or Polaris—can be used to generate cross section data used by the PARCS nodal core simulator for full-core neutronics calculations. This report presents an assessment of the accuracy of SCALE lattice physics codes for preparation of lattice physics data that are used to support simulator codes such as the NRC’s PARCS, for UO 2 -mixed oxide (MOX)/Zr fueled light water reactor (LWR) analyses. Due to the nature of lattice physics calculations, critical reactor experiment benchmarks cannot be modeled in explicit detail in a lattice physics code. However, this limitation does not mean that these measurement data are not usable for lattice physics studies. Therefore, either geometry approximations or axial buckling must be implemented to determine the critical water height. These modeling limitations have led to development of a three-phase assessment strategy. In the first phase, selected critical experiment benchmarks are modeled using the SCALE 3D continuous-energy (CE) Monte Carlo (MC) code KENO, which is the most rigorous neutron transport method available in SCALE, with no approximations in the spatial, angle, or energy treatments. Biases and statistical uncertainties in quantities of interest such as k eff and pin power distributions are determined by comparing CE KENO results to experimental data. This first phase in the accuracy assessment is to demonstrate and establish the use of CE KENO as a reference solution for the second phase. In the second phase, 14 numerical test suites are used to compare SCALE lattice physics calculations with CE KENO as a reference solution. In the third phase, SCALE lattice physics depletion calculations are performed, and the spent fuel isotopic results are compared with available radiochemical assay measurements. Isotopic measurement comparisons provide quantitative assessment of isotopic density distribution predictions with the depletion models in SCALE lattice physics codes. This report documents results for all test suites. The assessment was performed using standard production techniques unless otherwise noted. Both TRITON/NEWT and Polaris exhibited acceptable accuracy for most test cases. For the few test cases in which acceptable accuracy criteria were not met, further code and data development are planned. The computer codes used in this assessment are as follows: SCALE 6.2 rev19189 (pre-release of 6.2.1) was used for CE KENO, TRITON/NEWT and Polaris (PWR only) calculations; SCALE 6.2.2 was used for Polaris calculations to address several updates in support of support boiling water reactor (BWR) lattice geometries. For this work, there is no difference between SCALE 6.2.1 and SCALE 6.2.2 beyond the inclusion of the Polaris BWR analysis capability; The ENDF/B-VII.1 continuous energy and 252 group libraries deployed in SCALE 6.2 were generated with AMPX 6.2 that is distributed with SCALE.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Electronic Structure of the Primary Electron Donor P 700 +• in Photosystem I Studied by Multifrequency HYSCORE Spectroscopy at X- and Q-Band

The primary electron donor P 700 of the photosystem I ( PS I) is a heterodimer consisting of two chlorophyll molecules. A series of electron-transfer events immediately following the initial light excitation leads to a stabilization of the positive charge by its cation radical form, P 700 +• . Here, the electronic structure of P 700 +• and, in particular, its asymmetry with respect to the two chlorophyll monomers is of fundamental interest and is not fully understood up to this date. Here, we apply multifrequency X- (9 GHz) and Q-band (35 GHz) hyperfine sublevel correlation (HYSCORE) spectroscopy to investigate the electron spin density distribution in the cation radical P 700 +• of PS I from a thermophilic cyanobacterium Thermosynechococcus elongatus . Six 14 N and two 1 H distinct nuclei have been resolved in the HYSCORE spectra and parameters of the corresponding nuclear hyperfine and quadrupolar hyperfine interactions were obtained by combining the analysis of HYSCORE spectral features with direct numerical simulations. Based on a close similarity of the nuclear quadrupole tensor parameters, all of the resolved 14 N nuclei were assigned to six out of total eight available pyrrole ring nitrogen atoms (i.e., four in each of the chlorophylls), providing direct evidence of spin density delocalization over the both monomers in the heterodimer. Using the obtained experimental values of the 14 N electron-nuclear hyperfine interaction parameters, the upper limit of the electron spin density asymmetry parameter is estimated as R A/B upper = 7.7 ± 0.5, while a tentative assignment of 14 N observed in the HYSCORE spectra yields R B/A = 3.1 ± 0.5.

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

Study of the Cosmic Rays and Interstellar Medium in Local H I Clouds Using Fermi -LAT Gamma-Ray Observations

An accurate estimate of the interstellar gas density distribution is crucial to understanding the interstellar medium (ISM) and Galactic cosmic rays (CRs). To comprehend the ISM and CRs in a local environment, a study of the diffuse $γ$-ray emission in a midlatitude region of the third quadrant was performed. Specifically, the $γ$-ray data in the 0.1–25.6 GeV energy range of the Fermi Large Area Telescope (LAT) and other interstellar gas tracers such as the HI4PI survey data and the Planck dust thermal emission model were used, and the northern and southern regions were analyzed separately. The variation of the dust emission D em with the total neutral gas column density N H was studied in high dust temperature areas, and the N H /D em ratio was calibrated using $γ$-ray data under the assumption of a uniform CR intensity in the studied regions. The measured integrated $γ$-ray emissivities above 100 MeV are (1.58 ± 0.04) x 10 -26 photons s -1 sr -1 H-atom-1 and (1.59 ± 0.02) x 10-26 photons s-1 sr-1 H-atom -1 in the northern and southern regions, respectively, supporting the existence of a uniform CR intensity in the vicinity of the solar system. While most of the gas can be interpreted to be H I with a spin temperature of T S = 125 K or higher, an area dominated by optically thick H I with T S ~ 40 K was identified.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Redshift inference from the combination of galaxy colours and clustering in a hierarchical Bayesian model – Application to realistic N -body simulations

ABSTRACT Photometric galaxy surveys constitute a powerful cosmological probe but rely on the accurate characterization of their redshift distributions using only broad-band imaging, and can be very sensitive to incomplete or biased priors used for redshift calibration. A hierarchical Bayesian model has recently been developed to estimate those from the robust combination of prior information, photometry of single galaxies, and the information contained in the galaxy clustering against a well-characterized tracer population. In this work, we extend the method so that it can be applied to real data, developing some necessary new extensions to it, especially in the treatment of galaxy clustering information, and we test it on realistic simulations. After marginalizing over the mapping between the clustering estimator and the actual density distribution of the sample galaxies, and using prior information from a small patch of the survey, we find the incorporation of clustering information with photo-z’s tightens the redshift posteriors and overcomes biases in the prior that mimic those happening in spectroscopic samples. The method presented here uses all the information at hand to reduce prior biases and incompleteness. Even in cases where we artificially bias the spectroscopic sample to induce a shift in mean redshift of $\Delta \bar{z} \approx 0.05,$ the final biases in the posterior are $\Delta \bar{z} \lesssim 0.003.$ This robustness to flaws in the redshift prior or training samples would constitute a milestone for the control of redshift systematic uncertainties in future weak lensing analyses.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

SPATIOTEMPORAL TRACKING OF WIDE AREA POWER OUTAGE FROM NIGHT-TIME LIGHT IMAGERY

Monitoring progress of power restoration following extreme events is essential for situational awareness about spatio-temporal distribution of populations without power and to help with response efforts. Because of the proprietary nature of restoration data, and the difficulty in obtaining power outage data from utility companies in near real-time (e.g., during hurricane Maria (2017)), this project used satellite derived nighttime lights data from Visible Infrared Imaging Radiometer Suite (VIIRS) Day/Night Band (DNB) to (i) monitor wide area power outage and (ii) estimate impacted customers over time to assist with restoration. A discussion of the methodology and its implementation during hurricanes Maria (2017) and Eta (2020) is presented in this paper. Future work will focus on calibrating estimated customers based on light intensity and density distribution, and generation of restoration profiles for emergency response.

Kar, Bandana↗
DOE OSTI2023

Simulation tools, first results and experimental status of the MURAVES experiment

The MUon RAdiography of VESuvius (MURAVES) project aims at the study of Mt. Vesuvius, an active and hazardous volcano near Naples, Italy, with the use of muons freely and abundantly produced by cosmic rays. In particular, the MURAVES experiment intends to perform muographic imaging of the internal structure of the summit of Mt. Vesuvius. The challenging measurement of the rock density distribution in its summit by muography, in conjunction with data from other geophysical techniques, can help model possible eruption dynamics. The MURAVES apparatus consists of an array of three independent and identical muon trackers, with a total sensitive area of 3 square meters. In each tracker, a sequence of 4 XY tracking planes made of plastic scintillators is complemented by a 60 cm thick lead wall inserted between the two downstream planes to improve rejection of background from low energy muons. The apparatus is currently acquiring data. This paper presents preliminary results from the analysis of the first data samples acquired with trackers pointing towards Mt. Vesuvius, including the first relative measurement of the density projection of two flanks of the volcano at three different altitudes; we also present the workflow of the simulation chain of the MURAVES experiment and its ongoing developments.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

A small core in Vesta inferred from Dawn’s observations

Vesta’s large-scale interior structure had previously been constrained primarily using the gravity and shape data from the Dawn mission. However, these data alone still allow a wide range of possibilities for the differentiation state of the body. The moment of inertia is arguably the most diagnostic parameter related to the radial density distribution of a planetary body, making it crucial for assessing the body’s state of internal differentiation. Determining the moment of inertia requires additional measurements of the amplitudes of small rotational motions, such as precession and nutation. Here we report an updated estimate of the moment of inertia of Vesta inferred from Dawn’s Doppler tracking via the Deep Space Network and onboard imaging data. The recovered value for Vesta’s normalized polar moment of inertia is $\overline{C}$/MR 2 = 0.4208 ± 0.0047 (where M is the mass of Vesta and R is the reference radius), which is only 6.6% lower than the homogeneous value of 0.4505. This value, combined with the gravity field and global shape, suggests that Vesta’s interior has limited density stratification beneath its howardite–eucrite–diogenite-dominated crust. We propose two possible origin scenarios that are consistent with the observed constraints. In the first scenario, Vesta’s interior did not undergo full differentiation due to late accretion. In the second scenario, Vesta originated as an impact remnant of a larger differentiated body re-accreted with non-chondritic bulk composition produced from a catastrophic impact. Vesta did not experience complete differentiation in either scenario, suggesting that its current state reflects a complex interplay between its accretion timing, thermal evolution, redistribution of 26 Al bearing melt and/or impact processes.

CNEOS 2014-01-08 bolide↗
DOE OSTI2023

The evolution of HCO + in molecular clouds using a novel chemical post-processing algorithm

Modelling the chemistry of molecular clouds is critical to accurately simulating their evolution. To reduce computational cost, 3D simulations generally restrict their chemistry to species with strong heating and cooling effects. Time-dependent information about the evolution of other species is therefore often neglected. We address this gap by post-processing tracer particles in the SILCC-Zoom molecular cloud simulations. Using a chemical network of 39 species and 301 reactions (including freeze-out of CO and H 2 O) and a novel algorithm to reconstruct a density grid from sparse tracer particle data, we produce time-dependent density distributions for various species. We focus upon the evolution of HCO + , which is a critical formation reactant of CO but is not typically modelled on the fly. We find that ∼ 90 per cent of the HCO + content of the cold molecular gas forms in situ around n HCO + ∼ 10 3 –10 4 cm −3 , over a time-scale of approximately 1 Myr. The remaining ∼ 10 per cent forms at high extinction sites, with minimal turbulent mixing out into the less dense gas. We further show that the dominant HCO + formation pathway is dependent on the visual extinction, with the reaction H 3 + + CO contributing 90 per cent of the total HCO + production above A V, 3D = 3. We produce the very first maps of the HCO + column density, N(HCO + ), and show that it reaches values as high as 10 15 cm −2 . We find that 50 per cent of the HCO + mass is located within AV ∼ 10–30 in a density range of 10 3.5 –10 4.5 cm −3 . Our maps of N(HCO + ) are shown to be in good agreement with recent observations of the W49A star-forming region.

79 ASTRONOMY AND ASTROPHYSICS↗