SEARCH · Search NASA

Results for “column generation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 235 records · Page 13

All sources NASA NTRS DOE OSTI

NASA NTRS2014

Demonstration of AIRS Total Ozone Products to Operations to Enhance User Readiness

Cyclogenesis is a key forecast challenge at operational forecasting centers such as WPC and OPC, so these centers have a particular interest in unique products that can identify key storm features. In some cases, explosively developing extratropical cyclones can produce hurricane force, non-convective winds along the East Coast and north Atlantic as well as the Pacific Ocean, with the potential to cause significant damage to life and property. Therefore, anticipating cyclogenesis for these types of storms is crucial for furthering the NOAA goal of a "Weather Ready Nation". Over the last few years, multispectral imagery (i.e. RGB) products have gained popularity among forecasters. The GOES-R satellite champion at WPC/OPC has regularly evaluated the Air Mass RGB products from GOES Sounder, MODIS, and SEVIRI to aid in forecasting cyclogenesis as part of ongoing collaborations with SPoRT within the framework of the GOES-R Proving Ground. WPC/OPC has used these products to identify regions of stratospheric air associated with tropopause folds that can lead to cyclogenesis and hurricane force winds. RGB products combine multiple channels or channel differences into multi-color imagery in which different colors represent a particular cloud or air mass type. Initial interaction and feedback from forecasters evaluating the legacy Air Mass RGBs revealed some uncertainty regarding what physical processes the qualitative RGB products represent and color interpretation. To enhance forecaster confidence and interpretation of the Air Mass RGB, NASA SPoRT has transitioned a total column ozone product from AIRS retrievals to the WPC/OPC. The use of legacy AIRS demonstrates future JPSS capabilities possible with CrIS or OMPS. Since stratospheric air can be identified by anomalous potential vorticity and warm, dry, ozone-rich air, hyperspectral infrared sounder ozone products can be used in conjunction with the Air Mass RGB for identifying the role of stratospheric air in explosive cyclogenesis and hurricane force wind events. Currently, forecasters at WPC/OPC are evaluating the Air Mass RGB imagery in conjunction with the AIRS total column ozone to aid forecasting cyclogenesis and high wind forecasts. One of the limitations of the total ozone product is that it is difficult for forecasters to determine whether elevated ozone concentrations are related to stratospheric air or climatologically high values of ozone in certain regions. To address this limitation, SPoRT created an AIRS ozone anomaly product which calculates the percent of normal ozone based on a global stratospheric ozone mean climatology. With the knowledge that ozone values 125 percent of normal and greater typically represent stratospheric air; the anomaly product can be used with the total column ozone product to confirm regions of stratospheric air on the Air Mass RGB. This presentation describes the generation of these products along with forecaster feedback concerning the use of the AIRS ozone products in conjunction with the Air Mass RGB product for the unique forecast challenges WPC/OPC face. Additionally examples of CrIS ozone and anomaly products will be shown to further demonstrate the utility and capability of JPSS in forecasting unique events.

Berndt, Emily↗

NASA NTRS2015

Development and Application of Hyperspectral Infrared Ozone Retrieval Products for Operational Meteorology

Berndt, Emily↗

NASA NTRS2019

The Broad Absorption Line Tidal Disruption Event iPTF15af: Optical and Ultraviolet Evolution

We present multiwavelength observations of the tidal disruption event (TDE) iPTF15af, discovered by the intermediate Palomar Transient Factory survey at redshift z = 0.07897. The optical and ultraviolet (UV) light curves of the transient show a slow decay over 5 months, in agreement with previous optically discovered TDEs. It also has a comparable blackbody peak luminosity of L(sub peak) approx. = 1.5 x 10(exp 44) erg s(exp -1). The inferred temperature from the optical and UV data shows a value of (3–5) × 10(exp 4) K. The transient is not detected in X-rays up to L(sub X) < 3 x 10(exp 42) erg s(exp -1) within the first 5 months after discovery. The optical spectra exhibit two distinct broad emission lines in the He ii region, and at later times also Hα emission. Additionally, emission from [N iii] and [O iii] is detected, likely produced by the Bowen fluorescence effect. UV spectra reveal broad emission and absorption lines associated with high-ionization states of N v, C iv, Si iv, and possibly P v. These features, analogous to those of broad absorption line quasars (BAL QSOs), require an absorber with column densities N(sub H) > 10(exp 23) cm(exp -2). This optically thick gas would also explain the nondetection in soft X-rays. The profile of the absorption lines with the highest column density material at the largest velocity is opposite that of BAL QSOs. We suggest that radiation pressure generated by the TDE flare at early times could have provided the initial acceleration mechanism for this gas. Spectral UV line monitoring of future TDEs could test this proposal.

Blagorodnova, N.↗

NASA NTRS2011

3D Orbit Visualization for Earth-Observing Missions

This software visualizes orbit paths for the Orbiting Carbon Observatory (OCO), but was designed to be general and applicable to any Earth-observing mission. The software uses the Google Earth user interface to provide a visual mechanism to explore spacecraft orbit paths, ground footprint locations, and local cloud cover conditions. In addition, a drill-down capability allows for users to point and click on a particular observation frame to pop up ancillary information such as data product filenames and directory paths, latitude, longitude, time stamp, column-average dry air mole fraction of carbon dioxide, and solar zenith angle. This software can be integrated with the ground data system for any Earth-observing mission to automatically generate daily orbit path data products in Google Earth KML format. These KML data products can be directly loaded into the Google Earth application for interactive 3D visualization of the orbit paths for each mission day. Each time the application runs, the daily orbit paths are encapsulated in a KML file for each mission day since the last time the application ran. Alternatively, the daily KML for a specified mission day may be generated. The application automatically extracts the spacecraft position and ground footprint geometry as a function of time from a daily Level 1B data product created and archived by the mission s ground data system software. In addition, ancillary data, such as the column-averaged dry air mole fraction of carbon dioxide and solar zenith angle, are automatically extracted from a Level 2 mission data product. Zoom, pan, and rotate capability are provided through the standard Google Earth interface. Cloud cover is indicated with an image layer from the MODIS (Moderate Resolution Imaging Spectroradiometer) aboard the Aqua satellite, which is automatically retrieved from JPL s OnEarth Web service.

Jacob, Joseph C.↗

DOE OSTI2022

Real-space density kernel method for Kohn–Sham density functional theory calculations at high temperature

Kohn–Sham density functional theory calculations using conventional diagonalization based methods become increasingly expensive as temperature increases due to the need to compute increasing numbers of partially occupied states. In this work, we present a density matrix based method for Kohn–Sham calculations at high temperatures that eliminates the need for diagonalization entirely, thus reducing the cost of such calculations significantly. Specifically, we develop real-space expressions for the electron density, electronic free energy, Hellmann–Feynman forces, and Hellmann–Feynman stress tensor in terms of an orthonormal auxiliary orbital basis and its density kernel transform, the density kernel being the matrix representation of the density operator in the auxiliary basis. Using Chebyshev filtering to generate the auxiliary basis, we next develop an approach akin to Clenshaw–Curtis spectral quadrature to calculate the individual columns of the density kernel based on the Fermi operator expansion in Chebyshev polynomials and employ a similar approach to evaluate band structure and entropic energy components. We implement the proposed formulation in the SPARC electronic structure code, using which we show systematic convergence of the aforementioned quantities to exact diagonalization results, and obtain significant speedups relative to conventional diagonalization based methods. Finally, we employ the new method to compute the self-diffusion coefficient and viscosity of aluminum at 116 045 K from Kohn–Sham quantum molecular dynamics, where we find agreement with previous more approximate orbital-free density functional methods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1985

Electrostatic ion-cyclotron waves in a nonuniform magnetic field

The properties of electrostatic ion-cyclotron waves excited in a single-ended cesium Q machine with a nonuniform magnetic field are described. The electrostatic ion-cyclotron waves are generated in the usual manner by drawing an electron current to a small exciter disk immersed in the plasma column. The parallel and perpendicular (to B) wavelengths and phase velocities are determined by mapping out two-dimensional wave phase contours. The wave frequency f depends on the location of the exciter disk in the nonuniform magnetic field, and propagating waves are only observed in the region where f is approximately greater than fci, where fci is the local ion-cyclotron frequency. The parallel phase velocity is in the direction of the electron drift. From measurements of the plasma properties along the axis, it is inferred that the electron drift velocity is not uniform along the entire current channel. The evidence suggests that the waves begin being excited at that axial position where the critical drift velocity is first exceeded, consistent with a current-driven excitation mechanism.

Cartier, S. L.↗

NASA NTRS1990

Flow fields of low pressure vent exhausts

The flow field produced by low pressure gas vents are described based on experimental data obtained from tests in a large vacuum chamber. The gas density, pressure, and flux at any location in the flow field are calculated based on the vent plume description and the knowledge of the flow rate and velocity of the venting gas. The same parameters and the column densities along a specified line of sight traversing the plume are also obtained and shown by a computer generated graphical representation. The fields obtained with a radically scanning Pitot probe within the exhausting gas are described by a power of the cosine function, the mass rate, and the distance from the exit port. The field measurements were made for gas at pressures ranging from 2 to 50 torr venting from pipe fittings with diameters to 3/16 to 1-1/2 inches I.D. (4.76 to 38.1 mm). The N2 mass flow rates ranged from 2E-4 to 3.7E-1 g/s.

Scialdone, John J.↗

DOE OSTI2024

Task-specific sensor optical designs

A method and system architecture for designing a compressive sensing matrix for machine learning includes receiving an image associated with a classification task and; generating a sensing matrix. The sensing matrix includes an array of nonzero elements of the image. A prism array of prism elements is in communication with the sensing matrix. A row of values corresponding with an input angle of the prism array is mapped to a respective column corresponding with a detector. Then the detector detects light refracted at an output angle dictated by the physical shape of the prism element. A physical model of the detector is fabricated and generates a compressed representation of the image. A machine learning classification algorithm is applied to the compressed representation of the image and generates an optimized non-invertible final determination of the image.

Birch, Gabriel Carlisle↗

NASA NTRS1996

Development of Collaborative Research Initiatives to Advance the Aerospace Sciences-via the Communications, Electronics, Information Systems Focus Group

The primary goal of the Adaptive Vision Laboratory Research project was to develop advanced computer vision systems for automatic target recognition. The approach used in this effort combined several machine learning paradigms including evolutionary learning algorithms, neural networks, and adaptive clustering techniques to develop the E-MOR.PH system. This system is capable of generating pattern recognition systems to solve a wide variety of complex recognition tasks. A series of simulation experiments were conducted using E-MORPH to solve problems in OCR, military target recognition, industrial inspection, and medical image analysis. The bulk of the funds provided through this grant were used to purchase computer hardware and software to support these computationally intensive simulations. The payoff from this effort is the reduced need for human involvement in the design and implementation of recognition systems. We have shown that the techniques used in E-MORPH are generic and readily transition to other problem domains. Specifically, E-MORPH is multi-phase evolutionary leaming system that evolves cooperative sets of features detectors and combines their response using an adaptive classifier to form a complete pattern recognition system. The system can operate on binary or grayscale images. In our most recent experiments, we used multi-resolution images that are formed by applying a Gabor wavelet transform to a set of grayscale input images. To begin the leaming process, candidate chips are extracted from the multi-resolution images to form a training set and a test set. A population of detector sets is randomly initialized to start the evolutionary process. Using a combination of evolutionary programming and genetic algorithms, the feature detectors are enhanced to solve a recognition problem. The design of E-MORPH and recognition results for a complex problem in medical image analysis are described at the end of this report. The specific task involves the identification of vertebrae in x-ray images of human spinal columns. This problem is extremely challenging because the individual vertebra exhibit variation in shape, scale, orientation, and contrast. E-MORPH generated several accurate recognition systems to solve this task. This dual use of this ATR technology clearly demonstrates the flexibility and power of our approach.

Knasel, T. Michael↗

NASA NTRS1981

Evaluation of SLAR and thematic mapper MSS data for forest cover mapping using computer-aided analysis techniques

A set of training statistics for the 30 meter resolution simulated thematic mapper MSS data was generated based on land use/land cover classes. In addition to this supervised data set, a nonsupervised multicluster block of training statistics is being defined in order to compare the classification results and evaluate the effect of the different training selection methods on classification performance. Two test data sets, defined using a stratified sampling procedure incorporating a grid system with dimensions of 50 lines by 50 columns, and another set based on an analyst supervised set of test fields were used to evaluate the classifications of the TMS data. The supervised training data set generated training statistics, and a per point Gaussian maximum likelihood classification of the 1979 TMS data was obtained. The August 1980 MSS data was radiometrically adjusted. The SAR data was redigitized and the SAR imagery was qualitatively analyzed.

Hoffer, R. M.↗

NASA NTRS2020

Dust and gas absorption in the High Mass X-ray Binary IGR J16318−4848

Context. With an absorption column density on the order of 10(exp 24)per sq. cm, IGR J16318-4848 is one of the most extreme cases of a highly obscured High Mass X-ray Binary. Besides the overall continuum absorption, the source spectrum exhibits a strong iron and nickel fluorescence line complex at 6.4 keV. Previous empirical modeling of these features and comparison with radiative transfer simulations raised questions about the structure and covering fraction of the absorber and the profile of the fluorescence lines. Aims. We aim at a self-consistent description of the continuum absorption, the absorption edges and the fluorescent lines to constrain properties of the absorbing material, such as ionization structure and geometry. We further investigate the effects of dust absorption on the observed spectra and the possibility of fluorescent emission from dust grains. Methods. We use XMM-Newton and NuSTAR spectra to first constrain empirically the incident continuum and fluorescence lines. Next, we use XSTAR to construct a customized photoionization model where we vary the ionization parameter, column density, and covering fraction. In a third step, we model the absorption and fluorescence in a dusty olivine absorber and employ both, a simple analytical model for the fluorescent line emission and Monte Carlo radiative transfer, spectral shapes and line fluxes that are very close to the data are generated. Results. Our empirical spectral modeling is in agreement with previous works. Our second model, the single gas absorber does not describe the observational data. In particular, irrespective of the ionization state or column density of the absorber, a much higher covering fraction than previously estimated is needed to produce the strong fluorescence lines and the large continuum absorption. A dusty, spherical absorber (modeled as consisting of olivine dust, although the nature of dust cannot be constrained) is well able to produce the observed continuum absorption and edges. Conclusions. A dense, dusty absorber in the direct vicinity of the source consisting of dust offers a consistent description of both the strong continuum absorption and the strong emission features in the X-ray spectrum of IGR J16318􀀀4848. In particular, for low optical depth of individual grains the dust will contribute significantly to the fluorescent emission, which is the case for typical densities and grain size distribution models.

Ralf Ballhausen↗

DOE OSTI2021

Processed Lab Data for Neural Network-Based Shear Stress Level Prediction

Machine learning can be used to predict fault properties such as shear stress, friction, and time to failure using continuous records of fault zone acoustic emissions. The files are extracted features and labels from lab data (experiment p4679). The features are extracted with a non-overlapping window from the original acoustic data. The first column is the time of the window. The second and third columns are the mean and the variance of the acoustic data in this window, respectively. The 4th-11th column is the the power spectrum density ranging from low to high frequency. And the last column is the corresponding label (shear stress level). The name of the file means which driving velocity the sequence is generated from. Data were generated from laboratory friction experiments conducted with a biaxial shear apparatus. Experiments were conducted in the double direct shear configuration in which two fault zones are sheared between three rigid forcing blocks. Our samples consisted of two 5-mm-thick layers of simulated fault gouge with a nominal contact area of 10 by 10 cm^2. Gouge material consisted of soda-lime glass beads with initial particle size between 105 and 149 micrometers. Prior to shearing, we impose a constant fault normal stress of 2 MPa using a servo-controlled load-feedback mechanism and allow the sample to compact. Once the sample has reached a constant layer thickness, the central block is driven down at constant rate of 10 micrometers per second. In tandem, we collect an AE signal continuously at 4 MHz from a piezoceramic sensor embedded in a steel forcing block about 22 mm from the gouge layer The data from this experiment can be used with the deep learning algorithm to train it for future fault property prediction.

15 GEOTHERMAL ENERGY↗

NASA NTRS1989

Flow fields of low pressure vent exhausts

The flow field produced by low pressure gas vents are described based on experimental data obtained from tests in a large vacuum chamber. The gas density, pressure, and flux at any location in the flow field are calculated based on the vent plume description and the knowledge of the flow rate and velocity of the venting gas. The same parameters and the column densities along a specified line of sight traversing the plume are also obtained and shown by a computer-generated graphical representation. The fields obtained with a radially scanning Pitot probe within the exhausting gas are described by a power of the cosine function, the mass rate and the distance from the exit port. The field measurements were made for gas at pressures ranging from 2 to 50 torr venting from pipe fittings with diameters of 3/16 inch to 1-1/2 inches I.D. (4.76 mm to 38.1 mm). The N(2) mass flow rates ranged from 2E-4 to 3.7E-1 g/s.

Scialdone, John J.↗

NASA NTRS1995

Cosmic-Ray Detectors With Interdigitated Electrodes

Detectors measure both positions of incidence and energies of incident charged particles. Stack of detector wafers intercept cosmic ray. Measure positions of incidence to determine cosmic-ray trajectory and charge generated within them (proportional to cosmic-ray energy dissipated within them). Interdigital electrode pattern repeated over many rows and columns on tops of detector wafers in stack. Electrode pattern defines pixels within which points of incidence of incident cosmic rays located.

Cunningham, Thomas J.↗

NASA NTRS2000

Chandra Analysis with "CIAO"

CIAO (Chandra Interactive Analysis of Observations) is the name for the suite of tools and applications developed at the CXC for the analysis of data from the Chandra X-Ray Observatory and other missions. We describe some of the uses of CIAO for the analysis of Chandra (and other) data. There are several innovative aspects of CIAO. The user has the choice of a command line or a graphic-oriented approach. the use of the CXC "Data Model" results in: file-format independence: users have the ability of reading and writing files in FITS, QPOE, and IMH format raw and column filtering and binning can be applied to input file for *all* CIAO tools and applications avoiding the need to generate subsidiary data files on disk all filters are recorded in the "data subspace" of output files which is then recognized in the subsequent analysis. CIAO, in its current release, includes several tools both instrument specific (for ACIS, HRC and grating analysis) and for generic data analysis (data manipulation, source detection, transforms & convolutions, response, and timing) and several applications (Prism, Chips, and Sherpa plus Firstlook, ToolAgent and Filter Window in beta version) .

Elvis, Martin↗

NASA NTRS2004

A New Approach to Using a Cloud-Resolving Model to Study the Interactions between Clouds, Precipitation and Aerosols

Recently Grabowski and Khairoutdinov and Randall have proposed and demonstrated, the feasibility of using 2D CRMs (these did not include detailed land processes and used periodic lateral boundary conditions) as 'super-parameterizations or a multi-scale modeling framework MMF' for cloud processes within atmospheric generation circulation models (GCMs). In the MMF, a fine-resolution 2D CRM takes the place of the single-column parameterization that is used in conventional GCMs. An overview of this approach is given by Randall et al.. Since a CRM can explicitly simulate cloud processes at the natural space and time scales of clouddynamical processes, it's possible to compute statistics of cloud fields, including radiances and radar backscatters, directly simulated/produced by the CRM. A seed fund is available at NASA Goddard to build a MMF based on a 2D GCE model and the Goddard fvGCM. A prototype MMF will be developed at end of 2004. The major objectives are to test and evaluate its performance using the GCE model as the super parameterization and compare its results with other super parameterizations, conventional GCMs and observations. Recently, a detailed spectral-bin microphysical scheme was implemented into the Goddard Cumulus Ensemble (WE) model. The GCE model is a CRM and it has been improved and applied to study the precipitation processes in the past two decades. Atmospheric aerosols are also described using number density size-distribution functions. The model is tested by studying the evolution of deep tropical clouds in the west Pacific warm pool region and summertime convection over a mid-latitude continent with different concentrations of CCN: a low 'clean' concentration and a high 'dirty' concentration. The impact of atmospheric aerosol concentration on cloud and precipitation will be investigated. In addition, the similarities and differences between bulk- and spectral-bin microphysics will be presented,

Tao, Wei-Kuo↗

DOE OSTI2020

LLNL Kimberlina 1.2 NUFT Simulations June 2018 (v2)

This dataset contains the output 6,000, 3-dimensional reactive multi-phase flow and transport aquifer simulations of brine and CO2 leakage into a protective aquiver in California’s San Joaquin Valley and input data files detailing the geologic mesh, aquifer physical properties and CO2 and brine injection rates. This data set was generated as an ongoing effort with the US DOE National Risk Assessment Partnership (NRAP) to evaluate the effectiveness of monitoring techniques to detect brine and CO2 leakage from legacy wells into underground sources of drinking water overlaying a CO2 storage reservoir. Each simulation contains a unique set of input parameters, generated stochastically. The outputs consist of these upper three geologic layers (from top): the Etchegoin, Macoma-Chanac, Santa Margarita-McLure formations. These simulations span the several distances (1, 3 and 6 km or wells W31-0.2, W31-0.5 and W31-1.0, respectively) from the CO2 injector, initiated from bottom hole pressure and saturation to calculate wellbore leakage from the storage reservoir, with low and high regional groundwater gradients and wellbore leakage into 5 leaky nodes. The dataset includes 1,000 unique simulations for each distance, which each contain a unique aquifer heterogeneity, aquifer and caprock permeability, and two model generations are included with a high permeability (prod07) and hybrid permeability (prod09). The range of permeability distributions is listed in Table 1. Each model generation consists of 3,000 simulations. Included in the dataset are the leakage rates determined from 2D wellbore models which utilize the pressure and CO2 saturation from LBL's reservoir simulations, NUFT mesh files with distributed lithology, NUFT rocktab files which describe the material properties for the geologic layers and the NUFT input files and post-processed output 'ntab' files. Each ntab file contains spatial (rows) and temporal (columns) model output tables for each model cell, the locations (x,y,z) and dimensions for each cells (dx, dy, dz). Table 1. Permeability distribution ranges for prod07 and prod09 model generations Geologic Layer: Permeability Range (log10 m^2) prod07 prod09 Etchegoin -12.92 to -10.92 -13.70 to -11.44 Macoma-Chanac -12.72 to -10.72 -13.50 to -11.24 Santa Margarita-McLure -12.70 to -10.70 -13.48 to -11.22 The input files used to generate the model include which are included in the dataset are: Time series of CO2 leakage input into the model (ex: Q_brn.W31-0.2.sim1000.layers123.tab) Time series of CO2 leakage input into the model (ex: Q_CO2.W31-0.2.sim1000.layers123.tab) Physical properties of the aquifer materials detailing the aquifer porosity, solid density, partitioning coefficients, permeabilities and van-Genuchten parameters detailed in a NUFT rocktab file: (ex: sim1000.usnt.rocktab) Numerical mesh and geologic data assigned to each model cell detailed in a NUFT genmsh format (ex: sim1000.mesh_k16.prod07.trans.genmsh) The primary output parameters are: pH (use absolute value) Change in TDS (mg/kg) Change in Pressure (Pa) Change CO2 gas saturation (fraction range 0.0-1.0) for example, the directory /p/lscratchh/mansoor1/nrap/kimberlina/prod09/mainfiles/sim1000/W31- 0.2 contains: sim1000.W31-0.2.trans.pH.red.ntab sim1000.W31-0.2.no_bg.trans.TDS.red.ntab sim1000.W31-0.2.usnt.P.deltabg.red.ntab sim1000.W31-0.2.usnt.CO2_sat.deltabg.red.ntab Each row in the NTAB files consist of model output per numerical grid cell. Each output file contains 33 columns (variables), including the information of numerical records, geologic location and sizes and the simulated parameter values over time. The first 13 variables are about numerical records and relative geologic information for a simulation grid: 1. index: simulation index 2. i: the ith grid of x-axis 3. j: the ith grid of y-axis 4. k: the ith grid of z-axis 5. element_ref: element reference 6. nuft_ind: nuft index 7. x: grid location in the x axis direction 8. y: grid location in the y axis direction 9. z: grid location in the z axis direction 10. dx: grid length in the x axis direction 11. dy: grid length in the y axis direction 12. dz: grid length in the z axis direction 13. volume: volume of the simulation grid The remainder (14, 15, 16...) variables are the simulated parameter values over time, take Pressure as an example, are: 14. 0.0y: initial pressure per cell. 15. 10.0y: simulated pressure at the end of the 10th year. 16. 20.0y: simulated pressure at the end of the 20th year. ... (repeated for every 10 years until 200 years)... The model extends 10,000 m, 5,000 m and 1,411 m in the x,y and z dimensions, respectively. The mesh consists of 164,832 cells with mesh dimensions of 101 x 51 x 32 (nx, ny, nz), with cell dimensions ranging from 100 m laterally (along x and y-axis) and model layers are as designated in the z-axis: Layer 1: atmosphere (1e-30 m thick) Layer 2: upper caprock (10 m thick) Layers 3-13: Etchegoin (536.23 m thck) Layers 14-27: Macoma-Chanac (679.04 m thick) Layers 28-32: Santa Margarita-McLure (185.94 m thick) The wellbore is placed along node i=51, j=26, and extends vertically along 5 nodes from the top to the bottom of the model. Special instructions when extracting files: Each Gzip archive (ex: prod07.sim1000-sim00099.tar.gz) contains 100 simulations. Gzip archives should be transferred into base directories (ie. In Linux: mkdir prod07; mv prod07.*.tar.gz prod07/.) before extracting, or files will be overwritten. Each sub-simulation tree should have the following file structure pattern (using the linux 'tree' command): |-- prod07 | |-- sim0001 | |-- W31-0.2 | | |-- Q_brn.W31-0.2.sim0001.layers123.tab | | |-- Q_co2.W31-0.2.sim0001.layers123.tab | | |-- sim0001.W31-0.2.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.2.trans.pH.red.ntab | | |-- sim0001.W31-0.2.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.2.usnt.P.deltabg.red.ntab | |-- W31-0.5 | | |-- Q_brn.W31-0.5.sim0001.layers123.tab | | |-- Q_co2.W31-0.5.sim0001.layers123.tab | | |-- sim0001.W31-0.5.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.5.trans.pH.red.ntab | | |-- sim0001.W31-0.5.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.5.usnt.P.deltabg.red.ntab | |-- W31-1.0 | | |-- Q_brn.W31-1.0.sim0001.layers123.tab | | |-- Q_co2.W31-1.0.sim0001.layers123.tab | | |-- sim0001.W31-1.0.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-1.0.trans.pH.red.ntab | | |-- sim0001.W31-1.0.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-1.0.usnt.P.deltabg.red.ntab | |-- sim0001.mesh_k16.prod07.trans.genmsh Disclaimer This document was prepared as an account of work sponsored by an agency of the United States government. Neither the United States government nor Lawrence Livermore National Security, LLC, nor any of their employees makes any warranty, expressed or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States government or Lawrence Livermore National Security, LLC. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States government or Lawrence Livermore National Security, LLC, and shall not be used for advertising or product endorsement purposes. Lawrence Livermore National Laboratory is operated by Lawrence Livermore National Security, LLC, for the U.S. Department of Energy, National Nuclear Security Administration under Contract DE-AC52-07NA27344. This report was reviewed and released as LLNL-MI-753464.

aquifer↗

DOE OSTI2020

Dust and gas absorption in the high mass X-ray binary IGR J16318-4848

With an absorption column density on the order of 10 24 cm -2 , IGR J16318-4848 is one of the most extreme cases of a highly obscured high mass X-ray binary. In addition to the overall continuum absorption, the source spectrum exhibits a strong iron and nickel fluorescence line complex at 6.4 keV. Previous empirical modeling of these features and comparison with radiative transfer simulations raised questions about the structure and covering fraction of the absorber and the profile of the fluorescence lines. We aim at a self-consistent description of the continuum absorption, the absorption edges, and the fluorescence lines to constrain the properties of the absorbing material, such as ionization structure and geometry. We further investigate the effects of dust absorption on the observed spectra and the possibility of fluorescence emission from dust grains. We used XMM-Newton and NuSTAR spectra to first empirically constrain the incident continuum and fluorescence lines. Next we used XSTAR to construct a customized photoionization model where we vary the ionization parameter, column density, and covering fraction. In the third step we modeled the absorption and fluorescence in a dusty olivine absorber and employed both a simple analytical model for the fluorescence line emission and a Monte Carlo simulation of radiative transfer that generates line fluxes, which are very close to the observational data. Our empirical spectral modeling is in agreement with previous works. Our second model, the single gas absorber does not describe the observational data. In particular, irrespective of the ionization state or column density of the absorber, a much higher covering fraction than previously estimated is needed to produce the strong fluorescence lines and the large continuum absorption. A dusty, spherical absorber (modeled as consisting of olivine dust, although the nature of dust cannot be constrained) is able to produce the observed continuum absorption and edges. A dense, dusty absorber in the direct vicinity of the source consisting of dust offers a consistent description of both the strong continuum absorption and the strong emission features in the X-ray spectrum of IGR J16318-4848. In particular, for low optical depth of individual grains, which is the case for typical volume densities and grain size distribution models, the dust will contribute significantly to the fluorescence emission.

79 ASTRONOMY AND ASTROPHYSICS↗