Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Core logs”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Facies Analysis of the Prairie Du Chien Group in the Illinois Basin and Analogous Rocks in Missouri and Kentucky

Funded in 2023 by the U.S. Department of Energy’s Phase II Carbon Storage Assurance Facility Enterprise (CarbonSAFE) initiative, a Heidelberg Materials cement plant in Mitchell, Indiana, is currently being evaluated as a potential Carbon Capture and Storage (CCS) subsurface injection site. The Heidelberg CCS project targets the middle to upper Prairie du Chien Group (Early Ordovician) in southwestern Indiana. Assessment of reservoir feasibility requires collection of field data, seismic surveys, well-log correlation, geologic modeling, characterization well drilling, well testing, and reservoir simulation. However, the proposed Heidelberg CCS site is in a data-limited region, lacking both outcrop analogs and deep wells penetrating the target interval, which makes geologic modelling difficult prior to drilling a characterization well. To directly address this problem, the present study was undertaken to understand the sedimentologic composition and stratigraphic architecture of the Prairie du Chien Group from analogous outcrops and cores in the Illinois Basin and adjacent regions.

Ali, Shah Bilawal [Univ. of Illinois at Urbana-Cha↗
DOE OSTI2015

CALIBRATION OF THE MIXING-LENGTH THEORY FOR CONVECTIVE WHITE DWARF ENVELOPES

A calibration of the mixing-length parameter in the local mixing-length theory (MLT) is presented for the lower part of the convection zone in pure-hydrogen-atmosphere white dwarfs. The parameterization is performed from a comparison of three-dimensional (3D) CO5BOLD simulations with a grid of one-dimensional (1D) envelopes with a varying mixing-length parameter. In many instances, the 3D simulations are restricted to the upper part of the convection zone. The hydrodynamical calculations suggest, in those cases, that the entropy of the upflows does not change significantly from the bottom of the convection zone to regions immediately below the photosphere. We rely on this asymptotic entropy value, characteristic of the deep and adiabatically stratified layers, to calibrate 1D envelopes. The calibration encompasses the convective hydrogen-line (DA) white dwarfs in the effective temperature range 6000 ≤ T {sub eff} (K) ≤15, 000 and the surface gravity range 7.0 ≤ log g ≤ 9.0. It is established that the local MLT is unable to reproduce simultaneously the thermodynamical, flux, and dynamical properties of the 3D simulations. We therefore propose three different parameterizations for these quantities. The resulting calibration can be applied to structure and envelope calculations, in particular for pulsation, chemical diffusion, and convective mixing studies. On the other hand, convection has no effect on the white dwarf cooling rates until there is a convective coupling with the degenerate core below T {sub eff} ∼ 5000 K. In this regime, the 1D structures are insensitive to the MLT parameterization and converge to the mean 3D results, hence they remain fully appropriate for age determinations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

A Chandra study of the image power spectra of 41 cool core and non-cool core galaxy clusters

In this work we propose a new diagnostic to segregate cool core (CC) clusters from non-CC (NCC) clusters by studying the two-dimensional power spectra of the X-ray images observed with the Chandra X-ray observatory. Our sample contains 41 members (z=0.01--0.54) which are selected from the Chandra archive when a high photon count, an adequate angular resolution, a relatively complete detector coverage, and coincident CC–NCC classifications derived with three traditional diagnostics are simultaneously guaranteed. We find that in the log–log space the derived image power spectra can be well represented by a constant model component at large wavenumbers, while at small wavenumbers a power excess beyond the constant component appears in all clusters, with a clear tendency that the excess is stronger in CC clusters. By introducing a new CC diagnostic parameter, i.e., the power excess index (PEI), we classify the clusters in our sample and compare the results with those obtained with three traditional CC diagnostics. We find that the results agree with each other very well. By calculating the PEI values of the simulated clusters, we find that the new diagnostic works well at redshifts up to 0.5 for intermediately sized and massive clusters with a typical Chandra or XMM-Newton pointing observation. The new CC diagnostic has several advantages over its counterparts, e.g., it is free of the effects of the commonly seen centroid shift of the X-ray halo caused by merger event, and the corresponding calculation is straightforward, almost irrelevant to the complicated spectral analysis.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2000

Physiologic Responses Produced by Active and Passive Personal Cooling Vests

Personal thermoregulatory systems which provide chest cooling are used in the industrial and aerospace environments to alleviate thermal stress. However, little information is available regarding the physiologic and circulatory changes produced by routine operation of these systems. The objectives of this study were to document and compare the subjects' response to three cooling vests in their recommended configurations. The Life Enhancement Tech (LET) lightweight active cooling vest with cap, the MicroClimate Systems Change of Phase garment (MCS), and the Steele Vest were each used to cool the chest regions of 12 male and 8 female Healthy subjects (21 to 69 yr.) in this study. The subjects, seated in an upright position at normal room temperature (approx. 22 C), were tested for 60 min. with one of the cooling garments. The LET active garment had an initial coolant fluid inlet temperature of 60 F, and was ramped down to 50 F. Oral, right and left ear canal temperatures were logged manually every 5 min. Arm, leg, chest and rectal temperatures; heart rate; and respiration were recorded continuously on a U.F.I., Inc. Biolog ambulatory monitor. For men, all three vests had similar, significant cooling effects. Decreases in the average rectal temperature, oral temperature, and ear canal temperatures were approximately 0.2 C, 0.2 C and 0.1 C, respectively. In contrast to the men, the female subjects wearing the MCS and Steel vests had similar cooling responses in which the core temperature remained elevated and oral and ear canal temperatures did not drop. The LET active garment cooled most of the female subjects in this study; rectal, oral and ear temperature decreased about 0.2 C, 0.3 C and 0.3 C, respectively. These results show that the garment configurations tested do not elicit a similar thermal response in all subjects. A gender difference is evident. The LET active garment configuration was most effective in decreasing temperatures of the female subjects; the MCS vest was least effective. For male subjects, the three vests appear to be more nearly equivalent. The active garment system under study included a cooling cap, which may account for some of the difference in response.

Ku, Yu-Tsuan E.↗
NASA NTRS1972

Chemical fractionation in iron meteorites and its interpretation.

Published analyses of trace and minor elements in iron meteorites have been compiled and the distributions interpreted with the chemical groups defined by Wasson. When each element is plotted against Ni on log scales, groups are often clearly resolved with all the members of a group falling within the limits of sampling and analytical error on a straight line. Two fractionations have occurred, a primary event which established the bulk composition of each group and a secondary event which fractionated the elements within each group. Group I appears to have escaped the secondary fractionation. An examination of possible fractionation mechanisms suggests that the secondary process took place during solidification of iron cores in the parent bodies. The elements Ir, Os, Pt, Ru and Rh would be enriched in the early Ni-poor solid whilst As, Au, Co, Mo, P, Pd and Sb would concentrate in later solid and produce the observed positive correlations with Ni.-

Scott, E. R. D.↗
NASA NTRS2004

Discovery of Extremely Embedded X-ray Sources in the R Coronae Australis Star Forming Core

We detected three extremely embedded X-ray sources in the R Corona Australis (R CrA) star forming core, IRS 7 region. Two weak X-ray sources are associated with the VLA centimeter radio sources 10E & W, whereas the third brightest source detected in the two XMM-Newton observations on March 2003 has no counterpart at any wavelengths. The large K-band upper-limit (19.4m) measured with the University of Hawaii 88-inch Telescope and strong absorption derived in X-rays (N(sub H) approx. 2.8 x 10(exp 23)/sq cm equivalent to A(sub v) approx. 180 m) indicate that the source is younger than typical Class I protostars, i.e. a Class 0 protostar or an intermittent phase between Class 0 and Class I protostars. The X-ray luminosity was less than one thirtieth (log L(sub x) less than or approx. equals 29.3 ergs/s) in the former Chandra observation in October 2000, which suggests that the X-ray activity, probably generated by magnetic activity, is triggered by an intermittent mass accretion episode such as FU Ori type outbursts. Because the source was detected at high significance in the XMM-Newton observations (approx. 2,000 cnts), X-ray properties of such young protostars can be well investigated for the first time. The light curves were constant in the 1st observation and increased linearly by a factor of two during 30 ksec in the 2nd observation. Both spectra showed iron K lines originated in hot thin-thermal plasma and fluorescence by cold gas. They can be reproduced by an absorbed thin-thermal plasma model with a Gaussian component at 6.4 keV (kT approx. 3-4 keV, L(sub x) approx. 7-20 x 10(exp 30) ergs/s). The rising timescale of the light curves in the 2nd observation was too slow for magnetically generated X-ray flares, whereas large equivalent width of the fluorescence iron K line in the 1st observation (approx. 810 eV) requires strong partial covering of the X-ray source. These results suggest that a confined hot (perhaps accretion) spot on the protostellar core was behind the star in the 1st observation and just appeared in the 2nd observation due to the core rotation with period of greater than or approx. 22.8 days, which is much slower than the break-up velocity previously assumed for young protostars. This means that the source had quiescent X-ray activity during the observations with an order of magnitude stronger level than Class I (older) protostars. We also consider whether the X-ray source associated with 10E could be shock heated plasma by a collision of a jet emanating from the tentative Class 0 object.

Hamaguchi, Ken-Ji↗
DOE OSTI2014

Tracing the outflow of a z = 0.334 FeLoBAL: New constraints from low-ionization absorbers in FBQS J1151+3822

We show for the first time that FBQS J1151+3822 is an iron low-ionization broad absorption line quasar (FeLoBAL QSO), the second-brightest and second-closest known example of this class. He I* and Fe II together act as an effective analytical tool, allowing us to obtain useful kinematic constraints from photoionization models of the outflow without needing to assume any particular acceleration model. The main outflow's log ionization parameter is –1.5, the log hydrogen density (cm{sup –3}) 5.5-8, the log hydrogen column density (cm{sup –2}) 21.7-21.9, the absorption radius 7.2-127 pc, and the kinetic luminosity is 0.16%-4.5% of the bolometric luminosity. We obtain line-of-sight covering fractions of ∼0.25 for strong Fe II, ∼0.5 for He I*, and ∼0.6 for Mg II. Narrower and shallower absorption lines from weaker Fe II and Mn II with an outflow velocity of ∼3400 km s{sup –1} have appeared between 2005 and 2011, suggesting that dense cores may have condensed inside the main outflow. Consideration of the literature might suggest that the FBQS J1151+3822 outflow is a member of a rare and distinct subclass of FeLoBALs with high densities and correspondingly small absorption radii. We find, however, that such outflows are not necessarily a distinct subclass, and that their apparent rarity could be a symptom of selection bias in studies using density-sensitive lines.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Physical property characteristics of gas hydrate-bearing reservoir and associated seal sediments collected during NGHP-02 in the Krishna-Godavari Basin, in the offshore of India

India’s National Gas Hydrate Program Expedition 02 (NGHP-02) was conducted to better understand geologic controls on gas hydrate occurrence and morphology, targeting coarse-grained sediments along the lower continental slope offshore eastern India. This study combines seismic, logging-while-drilling data, and a petroleum system approach to provide a regional geologic and lithologic context for: 1) gas hydrate morphology and distribution, and 2) effects of fine-grained sediments (clays and other grains smaller than 63µm) on gas production in NGHP-02 Area B in the Krishna-Godavari Basin. Area B seismic data show a buried anticline/syncline structure with strong reflectors, R1 and R2, that delineate two of the five lithologic units: Unit I (shallowest), II (R1), III, IV (R2) and V (below the bas of gas hydrate stability). Throughout Area B, gas hydrate morphology depends on its placement within these units. Specifically, core- and grain-scale measurements indicate fines content exerts a primary control on the gas hydrate distribution and morphology. Units I, II and III are generally fine-grained. On the anticline crest, these units host primarily grain-displacing gas hydrate veins linked to pore-occupying gas hydrate in thin, localized, coarser-grained deposits. Diatoms in Unit III increase porosity with depth, reaching ~70% where it contacts Unit IV, the gas hydrate reservoir associated with reflector R2. The Unit III lithology and porosity allow fluid and dissolved-phase methane to escape Unit IV and form gas hydrate in the fine-grained overburden. Within Unit IV, fine-grained layers are interbedded with coarser-grained gas hydrate reservoir sands, and the fines content even in the sands is high enough to impact hydraulic and mechanical properties during gas production. Fluid motion during gas production can mobilize fines, which can then clog pore throats, limiting production rates. In conclusion, pore-water freshening during gas hydrate dissociation can increase fines mobilization, particularly given the smectite identified in the fine-grained interbeds.

03 NATURAL GAS↗
DOE OSTI2014

The Potosi Reservoir Model 2013c, Property Modeling Update

As part of a larger project co-funded by the United States Department of Energy (US DOE) to evaluate the potential of formations within the Cambro-Ordovician strata above the Mt. Simon as potential targets for carbon sequestration in the Illinois and Michigan Basins, the Illinois Clean Coal Institute (ICCI) requested Schlumberger to evaluate the potential injectivity and carbon dioxide (CO2) plume size of the Cambrian Potosi Formation. The evaluation of this formation was accomplished using wireline data, core data, pressure data, and seismic data from this project as well as two other separately funded projects: the US DOE-funded Illinois Basin–Decatur Project (IBDP) being conducted by the Midwest Geological Sequestration Consortium (MGSC) in Macon County, Illinois, and the Illinois Industrial Carbon Capture and Sequestration (ICCS) project funded through the American Recovery and Reinvestment Act. In 2010, technical performance evaluations on the Cambrian Potosi Formation were performed through reservoir modeling. The data included formation tops from mud logs, well logs from the Verification Well #1 (VW1) and the Injection Well (CCS1), structural and stratigraphic formation from three dimensional (3D) seismic data, and field data from several waste water injection wells for Potosi Formation. The intention was for 2.2 million tons per annum (2 million tonnes per annum [MTPA]) of CO2 to be injected for 20 years. In the Task Error! Reference source not found., the 2010 Potosi heterogeneous model (referred to as the "Potosi Dynamic Model 2010") was re-run using a new injection scenario of 3.5 million tons per annum (3.2 MTPA) for 30 years. The extent of the Potosi Dynamic Model 2010, however, appeared too small for the new injection target. The models size was insufficient to accommodate the evolution of the plume. The new model, Potosi Dynamic Model 2013a, was built by extending the Potosi Dynamic Model 2010 grid to 30 by 30 mi (48 by 48 km), while preserving all property modeling workflows and layering. This model was retained as the base case. In the preceding Task [1], the Potosi reservoir model was updated to take into account the new data from the Verification Well #2 (VW2) which was drilled in 2012. The porosity and permeability modeling was revised to take into account the log data from the new well. Revisions of the 2010 modeling assumptions were also done on relative permeability, capillary pressures, formation water salinity, and the maximum allowable well bottomhole pressure. Dynamic simulations were run using the injection target of 3.5 million tons per annum (3.2 MTPA) for 30 years. This dynamic model was named Potosi Dynamic Model 2013b. In this Task, a new property modeling workflow was applied, where seismic inversion data guided the porosity mapping and geobody extraction. The static reservoir model was fully guided by PorosityCube interpretations and derivations coupled with petrophysical logs from three wells. The two main assumptions are: porosity features in the PorosityCube that correlate with lost circulation zones represent vugular zones, and that these vugular zones are laterally continuous. Extrapolation was done carefully to populate the vugular facies and their corresponding properties outside the seismic footprint up to the boundary of the 30 by 30 mi (48 by 48 km) model. Dynamic simulations were also run using the injection target of 3.5 million tons per annum (3.2 MTPA) for 30 years. This new dynamic model was named Potosi Dynamic Model 2013c. Reservoir simulation with the latest model gives a cumulative injection of 43 million tons (39 MT) in 30 years with a single well, which corresponds to 40% of the injection target. The injection rate is approx. 3.2 MTPA in the first six months as the well is injecting into the surrounding vugs, and declines rapidly to 1.8 million tons per annum (1.6 MTPA) in year 3 once the surrounding vugs are full and the CO2 start to reach the matrix. After, the injection rate declines gradually to 1.2 million tons per annum (1.1 MTPA) in year 18 and stays constant. This implies that a minimum of three (3) wells could be required in the Potosi to reach the injection target. The injectivity evaluated in this Task was higher compared to the preceding Task, since the current facies modeling (guided by the porosity map from the seismic inversion) indicated a higher density of vugs within the vugular zones. 5 As the CO2 follows the paths where vugs interconnection exists, a reasonably large and irregular plume extent was created. After 30 years of injection, the plume extends 13.7 mi (22 km) in E-W and 9.7 mi (16 km) in N-S directions. After injection finishes, the plume continues to migrate laterally, mainly driven by the remaining pressure gradient. After 60 years post-injection, the plume extends 14.2 mi (22.8 km) in E-W and 10 mi (16 km) in N-S directions, and remains constant as the remaining pressure gradient has become very low. Should the targeted cumulative injection of 106 million tons (96 MT) be achieved; a much larger plume extent could be expected. The increase of reservoir pressure at the end of injection is approximately 1,200 psia (8,274 kPa) around the injector and gradually decreases away from the well. The reservoir pressure increase is less than 10 psia (69 kPa) beyond 14 mi (23 km) away from injector. Should the targeted cumulative injection of 106 million tons (96 MT) be achieved; a much larger areal pressure increase could be expected. The reservoir pressure declines rapidly during the first 30 years post injection and the initial reservoir pressure is nearly restored after 100 years post-injection. The present evaluation is mainly associated with uncertainty on the vugs permeability and interconnectivity. The use of porosity mapping from seismic inversion might have reduced the uncertainty on the lateral vugs body distributions. However, major uncertainties on the Potosi vugs permeability remains. Therefore, injection test and pressure interference test among the wells could be considered to evaluate the local vugs permeability, extent, and interconnectivity. Facies modeling within the Potosi has yet to be thoroughly addressed. The carbonates during the time of deposition are believed to be regionally extensive. However, it may be worth delineating the reservoir with other regional wells or modern day analogues to understand the extent of the Potosi. More specifically, the model could incorporate lateral changes or trends if deemed necessary to represent facies transition. Data acquisitions to characterize the fracture pressure gradient, the formation water properties, the relative permeability, and the capillary pressure could also be considered in order to allow a more rigorous evaluation of the Potosi storage performance. A simulation using several injectors could also be considered to determine the required number of wells and appropriate spacing to achieve the injection target while taking into account the pressure interference.

Adushita, Yasmin↗
DOE OSTI2018

Not all moissanites are created equal: New constraints on moissanite from metamorphic rocks of Bulgaria

Here, terrestrial moissanite (SiC) is widely reported as an ultra-high pressure mineral occurring in kimberlites, diamonds and ultramafic/mafic rocks of mantle origin. However, the conditions of crystallization remain largely unknown. Moreover, dozens of SiC occurrences have been reported from continental crust sources such as granitoids, andesite–dacite volcanic rocks and their breccia, metasomatic and metamorphic rocks, and even limestones. The validity of many of these reports is still debated primarily due to possible contaminations from the widespread use of synthetic SiC abrasives in samples preparation. Indeed, reports of well-documented in-situ occurrences of moissanite in association with co-existing minerals are still scarce. The only condition of moissanite formation that is agreed upon is that extremely reducing media are required (e.g. 4.5–6 log units below the iron-wustite buffer). Here, we report the new occurrence of moissanite that was found in-situ within the garnet–staurolite–mica schists of Topolovgrad metamorphic group of Triassic age in Southern Bulgaria. The 10–300 μm moissanite crystals are situated within 0.1–1.2 mm isolated clusters, filled with amorphous carbon and nanocrystalline graphite. Most of moissanite crystals are 15R (rhombohedral) and 6H (hexagonal) polytypes, and one prismatic crystal, found within them, exhibits unusual concentric polytypical zoning with core (15R), intermediate zone (6H) and rim (3C-cubic). Experimental data show that this type of polytypical zonation is likely due to a decrease in temperature (or/and pressure?) and changes in Si/C ratio. Indeed, amphibolite facies metamorphism (500–580°C – garnet–staurolite zone) followed by a subsequent cooling during the retrograde stage of green schist facies metamorphism (~400–500°C) could have provided a change in temperature. The SiC containing clusters exhibit evidence that they are pre-metamorphic, and we hypothesize that their protolith was a “black shale” material likely rich in carbon, hydrocarbon and terrigenous silica. The latter served as a source of isolated chemically-reduced media, which is required for SiC formation. Other concepts to explain moissanite occurrences in metasedimentary rocks are also discussed. Importantly, our findings show that the formation conditions of moissanite are likely more variable than previously recognized.

15R-SiC↗
DOE OSTI2018

A Communication-Avoiding 3D LU Factorization Algorithm for Sparse Matrices

We propose a new algorithm to improve the strong scalability of right-looking sparse LU factorization on distributed memory systems. Our 3D sparse LU algorithm uses a three-dimensional MPI process grid, aggressively exploits elimination tree parallelism and trades off increased memory for reduced per-process communication. We also analyze the asymptotic improvements for planar graphs (e.g., from 2D grid or mesh domains) and certain non-planar graphs (specifically for 3D grids and meshes). For planar graphs with n vertices, our algorithm reduces communication volume asymptotically in n by a factor of O{log n} and latency by a factor of O{log n}. For non-planar cases, our algorithm can reduce the per-process communication volume by 3× and latency by O{n^1/3} times. In all cases, the memory needed to achieve these gains is a constant factor. We implemented our algorithm by extending the 2D data structure used in superLU. Our new 3D code achieves speedups up to 27× for planar graphs and up to 3.3× for non-planar graphs over the baseline 2D superLU when run on 24,000 cores of a Cray XC30.

Sao, Piyush↗
NASA NTRS1983

A note on compressibility and energy cascade in turbulent molecular clouds

Observed velocity-size correlations are reexamined in the light of an improved theory of turbulent energy cascade that is developed. It is shown that observed velocity-size correlations cannot be compared with the Kolmogorov law, which is based on incompressible turbulent flow. The fact that the log v-log(l/rho) scaling law (v the turbulent velocity, l the associated region size, and rho the fluid density) predicted for compressible energy cascade is always steeper than that observed in molecular clouds indicates the injection rather than the dissipation of mechanical energy at smaller scales of motion. It is also shown that the concept of strict energy cascade may not be generally applicable in the interstellar medium. The agreement between theory and observation turns out to be best for small cool clouds and cloud cores, suggesting that, for these regions at least, the dominant process in establishing the observed v-l-rho correlation is a turbulent energy cascade.

Fleck, R. C., Jr.↗
NASA NTRS2001

Determining CDOM Absorption Spectra in Diverse Aquatic Environments Using a Multiple Pathlength, Liquid Core Waveguide System

We evaluated the accuracy, sensitivity and precision of a multiple pathlength, liquid core waveguide (MPLCW) system for measuring colored dissolved organic matter (CDOM) absorption in the UV-visible spectral range (370-700 nm). The MPLCW has four optical paths (2.0, 9.8, 49.3, and 204 cm) coupled to a single Teflon AF sample cell. Water samples were obtained from inland, coastal and ocean waters ranging in salinity from 0 to 36 PSU. Reference solutions for the MPLCW were made having a refractive index of the sample. CDOM absorption coefficients, aCDOM, and the slope of the log-linearized absorption spectra, S, were compared with values obtained using a dual-beam spectrophotometer. Absorption of phenol red secondary standards measured by the MPLCW at 558 nm were highly correlated with spectrophotometer values and showed a linear response across all four pathlengths. Values of aCDOM measured using the MPLCW were virtually identical to spectrophotometer values over a wide range of concentrations. The dynamic range of aCDOM for MPLCW measurements was 0.002 - 231.5 m-1. At low CDOM concentrations spectrophotometric aCDOM were slightly greater than MPLCW values and showed larger fluctuations at longer wavelengths due to limitations in instrument precision. In contrast, MPLCW spectra followed an exponential to 600 nm for all samples.

Miller, Richard L.↗
DOE OSTI2018

Rancor Microworld Simulation Environment for Nuclear Process Control

The Rancor Microworld is a simplified nuclear power plant simulator. It represents the process of an integrated pressurized light water reactor core with a secondary Rankine cycle turbine/generator. The reactor core is governed by simplified core physics. Operators have the ability to raise and lower control rods to change the reactivity and consequently generate more or less electricity. The microworld provides a digital human machine interface for operating and monitoring the plant efficiently to maximize profits within the simulation. The simulation supports multi-unit operations analogous to small modular reactor plants and supports faults on various components within each unit to simulate some common training scenarios such as reactor and turbine trip events. The microworld includes experimental modules to log process parameters and both automated and user control inputs and events, administer situation awareness probes, and interface with eye tracking software.

Ulrich, ThomasA.↗
NASA NTRS2000

Determining CDOM Absorption Spectra in Diverse Coastal Environments Using a Multiple Pathlength, Liquid Core Waveguide System. Measuring the Absorption of CDOM in the Field Using a Multiple Pathlength Liquid Waveguide System

We evaluated the accuracy, sensitivity and precision of a multiple pathlength, liquid core waveguide (MPLCW) system for measuring colored dissolved organic matter (CDOM) absorption in the UV-visible spectral range (370-700 nm). The MPLCW has four optical paths (2.0, 9.8, 49.3, and 204 cm) coupled to a single Teflon AF sample cell. Water samples were obtained from inland, coastal and ocean waters ranging in salinity from 0 to 36 PSU. Reference solutions for the MPLCW were made having a refractive index of the sample. CDOM absorption coefficients, a(sub CDOM), and the slope of the log-linearized absorption spectra, S, were compared with values obtained using a dual-beam spectrophotometer. Absorption of phenol red secondary standards measured by the MPLCW at 558 nm were highly correlated with spectrophotometer values (r > 0.99) and showed a linear response across all four pathlengths. Values of a(sub CDOM) measured using the MPLCW were virtually identical to spectrophotometer values over a wide range of concentrations. The dynamic range of a(sub CDOM) for MPLCW measurements was 0.002 - 231.5/m. At low CDOM concentrations (a(sub 370) < 0.1/m) spectrophotometric a(sub CDOM) were slightly greater than MPLCW values and showed larger fluctuations at longer wavelengths due to limitations in instrument precision. In contrast, MPLCW spectra followed an exponential to 600 nm for all samples. The maximum deviation in replicate MPLCW spectra was less than 0.001 absorbance units. The portability, sampling, and optical characteristics of a MPLCW system provide significant enhancements for routine CDOM absorption measurements in a broad range of natural waters.

Miller, Richard L.↗
NASA NTRS2011

A Combined Subaru/VLT/MMT 1-5 Micrometer Study of Planets Orbiting HR 8799: Implications For Atmospheric Properties, Masses and Formation

We present new 1-1.25 micron (z and J band) Subaru/IRCS and 2 micron (K band) VLT/NaCo data for HR 8799 and a rereduction of the 3-5 micron MMT/Clio data first presented by Hinz et al. Our VLT/NaCo data yield a detection of a fourth planet at a projected separation of approximately 15 AU--"HR 8799e ." We also report new, albeit weak detections of HR 8799b at 1.03 micron and 3.3 micron. Empirical comparisons to field brown dwarfs show that at least HR 8799b and HR 8799c, and possibly HR 8799d, have near-to-mid-IR colors/ magnitudes significantly discrepant from the L/T dwarf sequence. Standard cloud deck atmosphere models appropriate for brown dwarfs provide only (marginally) statistically meaningful fits to HR 8799b and c for unphysically small radii. Models with thicker cloud layers not present in brown dwarfs reproduce the planets' spectral energy distributions far more accurately and without the need for resealing the planets' radii. Our preliminary modeling suggests that HR 8799b has log(g) = 4-4.5, T(sub eff) = 900 K. while HR 8799c, d, and (by inference) e have log(g) = 4-4.5, T(sub eff) = 1000-1200 K. Combining results from planet evolution models and new dynamical stability limits implies that the masses of HR 8799b, c, d, and e are 6-7 M(sub j), 7-10 M(sub j), 7-10 M(sub j), and 7-10 M(sub j). "Patchy" cloud prescriptions may provide even better fits to the data and may lower the estimated surface gravities and masses. Finally, contrary to some recent claims, forming the HR 8799 planets by core accretion is still plausible, although such systems are likely rare.

Currie, Thayne↗
DOE OSTI2019

A communication-avoiding 3D algorithm for sparse LU factorization on heterogeneous systems

We propose a new algorithm to improve the strong scalability of right-looking sparse LU factorization on distributed memory systems. Our 3D algorithm for sparse LU uses a three-dimensional MPI process grid, exploits elimination tree parallelism, and trades off increased memory for reduced per-process communication. We also analyze the asymptotic improvements for planar graphs (e.g., those arising from 2D grid or mesh discretizations) and certain non-planar graphs (specifically for 3D grids and meshes). For a planar graph with $n$ vertices, our algorithm reduces communication volume asymptotically in $n$ by a factor of $\mathscr{O}$ $\Big(\sqrt{log \ n}\Big)$ and latency by a factor of $\mathscr{O}$ $(log \ n)$. For nonplanar cases, our algorithm can reduce the per-process communication volume by 3× and latency by $\mathscr{O}$ $\Big(n^\frac13\Big)$ times. In all cases, the memory needed to achieve these gains is a constant factor. We implemented our algorithm by extending the 2D data structure used in SuperLU_DIST. Our new 3D code achieves empirical speedups up to 27× for planar graphs and up to 3.3× for non-planar graphs over the baseline 2D SuperLU_DIST when run on 24,000 cores of a Cray XC30. We extend the 3D algorithm for heterogeneous architectures by adding the Highly Asynchronous Lazy Offload (Halo) algorithm for co-processor offload. On 4096 nodes of a Cray XK7 with 32,768 CPU cores and 4096 Nvidia K20x GPUs, the 3D algorithm achieves empirical speedups up to 24× for planar graphs and 3.5× for non-planar graphs over the baseline 2D SuperLU_DIST with co-processor acceleration.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1993

The distance to the Large Magellanic Cloud - Constraints from Cepheids in Large Magellanic Cloud star clusters

We have used recent observational data for the Cepheids in the rich, young LMC clusters NGC 1866 and NGC 2031 to constrain the cluster distances with the mass-equivalency (ME) method. The basis of this approach is to fix the cluster distances by requiring the Cepheid evolutionary and pulsational masses to be equal. Using evolutionary models incorporating a mild amount of core and envelope overshooting along with recent pulsational models, we derive distance moduli of 18.51 +/- 0.21 and 18.32 +/- 0.20 for NGC 1866 and NGC 2031, respectively. The quoted errors are dominated by the uncertainties in the heavy element abundances of the clusters (assumed to be 0.3 dex for both clusters), with a smaller contribution due to the apparently intrinsic spread in the masses of the Cepheids in each cluster. For the ME method, we find that Delta(m-M)(0)/Delta sub Z(1) = 0.69, where Z(1) = log (Z/0.016). This result implies that the cluster distances can be determined to better than +/- 5 percent if the cluster abundances can be measured to better than about +/- 0.15 dex.

Bertelli, Gianpaolo↗