SEARCH · Search NASA

Results for “Image coding”

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 1,711 records · Page 95

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Scalable Hybrid Learning Techniques for Scientific Data Compression

Data compression is becoming critical for storing scientific data because many scientific applications need to store large amounts of data and post process this data for scientific discovery. Unlike image and video compression algorithms that limit errors to primary data (PD), scientists require compression techniques that accurately preserve derived quantities of interest (QoIs). Here, this article presents a physics-informed compression technique implemented as an end-to-end, scalable, GPU-based pipeline for data compression that addresses this requirement. Our hybrid compression technique combines machine learning techniques and standard compression methods. Specifically, we combine an autoencoder, an error-bounded lossy compressor to provide guarantees on raw data error, and a constraint satisfaction post-processing step to preserve the QoIs within a minimal error (generally less than floating point error). The effectiveness of the data compression pipeline is demonstrated by compressing nuclear fusion simulation data generated by a large-scale fusion code, XGC, which produces hundreds of terabytes of data in a single day. Our approach works within the ADIOS framework and results in compression by a factor of more than 150 while requiring only a few percent of the computational resources necessary for generating the data, making the overall approach highly effective for practical scenarios.

ITER↗

DOE OSTI2013

Short-Term Monitoring Results for Advanced New Construction Test House -- Roseville, California

A builder (K. Hovnanian Homes), design consultant, and trades collaborated to identify a systems integrated measures package for a 2,253-ft2 slab-on-grade ranch house to achieve a modeled energy savings of 60% with respect to the Building America House Simulation Protocols (Hendron, R. and Engebrecht, C. 'Building America House Simulation Protocols.' Golden, CO: National Renewable EnergyLaboratory, 2010) while minimizing construction costs and without requiring changes to the drawing that would impact local code or zoning approval. The key building improvements were applying R-10 insulation to the slab edge, increasing exterior wall cavity insulation from R-13 to R-15, and increasing attic insulation from R-30 to R-38. Also, the air handling unit was relocated from the attic toconditioned space, and ductwork was relocated along the attic floor with an insulated bulkhead built above it. Short-term testing results showed that duct air leakage was low due to short duct runs and the placement of ductwork in conditioned space. However, during commissioning, the lack of access for servicing the ductwork and dampers in the bulkhead area prevented retroactive balancing ofindividual branches, resulting in significant differences between specified and measured airflow values for some duct runs. Thermal imaging results performed on the house when operating in both heating and cooling modes validated historic stratification issues of ceiling supply registers with high supply air temperatures. Long-term monitoring results will be detailed in a future report.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2020

Challenges in simulating ground interacting nuclear explosions

This paper summarizes recent above-ground nuclear explosion simulations as part of a broader effort to better characterize conditions within a fireball that may influence the chemical evolution of bomb materials and other materials entrained from the local explosion environment. A critical component of this work is validation against historic footage of atmospheric testing, requiring that we understand how the frequency-dependent sensitivity of the utilized film footage influences data captured in such images. We focus first on the early physics of a nuclear explosion in the atmosphere before discussing some of the technical challenges we seek to capture in late-time models that include more complex emplacement conditions and subsurface features. We discuss required physics packages (compressible hydrodynamics, radiation transport, as well as necessary ancillary tables such as equations of state (EOS) and opacities). Additionally, we note reasonable “shortcuts” one may make and their limitations, e.g., using ideal gas EOS, replacing spectrally resolved radiation with spectrally averaged radiation, and exchanging deterministic transport with diffusion. We then discuss an approach to achieving an equilibrated initial stress state for problems where buoyancy and subsurface lithostatic stress are important. Our methodology is presented in the context of LLNL’s ALE3D multiphysics code but may readily be implemented in other codes. In this paper, we start with a description of the challenges of NUDET simulations, followed by a presentation of the simulated intensity (flux) as it would appear on an analysis of the Dixie test. We then progressively introduce additional complexity in subsequent sections (near-surface burst and gravity initialization) before concluding.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

DOE OSTI2014

Pre-ignition confinement and deflagration violence in LX-10 and PBX 9501

In thermal explosions of the nitramine octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX)-based explosives LX-10 and PBX-9501, the pre-ignition spatial and temporal heating profile defines the ignition location. The ignition location then determines the extent of inertial confinement and the violence of the resulting deflagration. In this work, we present results of experiments in which ∼23 g cylinders of LX-10 and PBX 9501 in thin-walled aluminum confinement vessels were subjected to identical heating profiles but which presented starkly different energy release signatures. Post-explosion LX-10 containment vessels were completely fragmented, while the PBX 9501 vessels were merely ruptured. Flash x-ray radiography images show that the initiation location for the LX-10 is a few mm farther from the end caps of the vessel relative to the initiation location of PBX 9501. This difference increases deflagration confinement for LX-10 at the time of ignition and extends the pressurization time during which the deflagration front propagates in the explosive. The variation in the initiation location, in turn, is determined by the thermal boundary conditions, which differ for these two explosives because of the larger coefficient of thermal expansion and greater thermal stability of the Viton binder in LX-10 relative to the estane and bis(2,2-dinitropropyl) acetal/formal binder of the PBX 9501. The thermal profile and initiation location were modeled for LX-10 using the hydrodynamics and structures code ALE3D; results indicate temperatures in the vicinity of the ignition location in excess of 274 °C near the time of ignition. The conductive burn rates for these two explosives, as determined by flash x-ray radiography, are comparable in the range 0.1–0.2 mm/μs, somewhat faster than rates observed by strand burner experiments for explosives in the temperature range 150–180 °C and pressures up to 100 MPa. The thinnest-wall aluminum containment vessels presented here rupture at lower pressures, in the range 10 MPa, suggesting that moderately higher temperatures and pressures are present near the deflagration front. For these explosives, however the most important property for determining deflagration violence is the degree of inertial confinement.

42 ENGINEERING↗

DOE OSTI2025

Speeding-up fuzzing through directional seeds

Abstract Fuzzing is an automated process for discovering inputs in a program that may trigger unexpected behavior. Today, fuzzing has become a standard practice for the discovery of bugs and security vulnerabilities. However, the main issue with such practices is that the exploration of the input space of programs can often be prohibitively expensive. Therefore, several alternative fuzzing strategies have been introduced during the last few years. Some fuzzing techniques rely on human expertise to provide a plausible set of initial input examples, namely, seeds. However, the process of handcrafting seeds for fuzzing purposes often becomes strenuous for humans as it requires a deeper understanding of the Program-Under-Test (PUT). Also, the use of known inputs to programs often does not trigger vulnerable program behavior or may not reach potentially vulnerable code locations. To address those issues, we propose a seed generation framework that enables Human-In-The-Loop (HITL) directed fuzzing where the human assumes a more active role in the creation of seeds that can penetrate and assess desired locations of the PUT. Our proposed framework uses Symbolic Execution (SE) to generate seeds that exercise paths to target program locations. Moreover, our framework enables the visualization of the explored execution paths in the binary of the PUT for the generated seeds. We evaluated our approach on a set of 12 carefully designed C programs with diverse characteristics that mimic real-world programs. The experimental results show the effectiveness of the proposed approach in improving the performance of standard fuzzing tools such as the American Fuzzy Lop ("Image missing" <#comment/> ). Specifically, our solution can generate seeds that substantially enhance the performance of the fuzzer, achieving speedups ranging from $$1.46\times $$ 1.46 × to $$68.53\times $$ 68.53 × for branch conditions, $$1.39\times $$ 1.39 × to $$254.62\times $$ 254.62 × for branch depths, $$14,879.59\times $$ 14 , 879.59 × to $$30,295.88\times $$ 30 , 295.88 × for branch widths over traditional seeds. Additionally, the speedup increases with the number of target function ranging from $$12,260\times $$ 12 , 260 × to $$22,856.07\times $$ 22 , 856.07 × over traditional seeds while only requiring less than 15 seconds on average for the seed generation step.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Wall heating by subcritical energetic electrons generated by the runaway electron avalanche source *

Abstract Subcritical energetic electrons (SEEs) produced by the runaway electron (RE) avalanche source at energies below the runaway threshold are found to be the primary contributor to surface heating of plasma-facing components (PFCs) during final loss events. This finding is supported by theoretical analysis, computational modeling with the Kinetic Orbit Runaway electrons Code (KORC), and qualitative agreement with DIII-D experimental observations. The avalanche source generates significantly more secondary electrons below the runaway threshold, which thermalize rapidly when well-confined. However, during a final loss event, the RE beam impacts the first wall, and SEEs are deconfined before they can thermalize. Additionally, because the energy deposition length decreases faster than energy, the deposited energy density, and thus the maximum PFC surface temperature change, is larger for SEEs than REs. KORC simulations employ an analytic first wall to model particle deconfinement onto a non-axisymmetric wall composed of individual tiles. PFC surface heating is calculated using a 1D model extended to include an energy-dependent deposition length scale. Simulations of DIII-D qualitatively agree with infrared (IR) imaging only when SEEs from the avalanche source are included. These results demonstrate that SEEs are the dominant contributor to PFC surface heating and indicate that the avalanche source plays a critical role in the PFC damage caused during final loss events. The prominence of SEEs also has important implications for interpreting IR imaging, one of the primary diagnostics for RE-wall interaction diagnosis, despite REs dominating the energy and current density. This result improves predictions of wall damage due to post-disruption REs to estimate material lifetime and design RE mitigation systems for ITER and future reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2013

Short-Term Monitoring Results for Advanced New Construction Test House—Roseville, California

A builder (K. Hovnanian® Homes®), design consultant, and trades collaborated to identify a systems integrated measures package for a 2,253-ft 2 slab-on-grade ranch house to achieve a modeled energy savings of 60% with respect to the Building America House Simulation Protocols, while minimizing construction costs and without requiring changes to the drawing that would impact local code or zoning approval. The key building improvements were applying R-10 insulation to the slab edge, increasing exterior wall cavity insulation from R-13 to R-15, and increasing attic insulation from R-30 to R-38. Also, the air handling unit was relocated from the attic to conditioned space, and ductwork was relocated along the attic floor with an insulated bulkhead built above it. Short-term testing results showed that duct air leakage was low due to short duct runs and the placement of ductwork in conditioned space. However, during commissioning, the lack of access for servicing the ductwork and dampers in the bulkhead area prevented retroactive balancing of individual branches, resulting in significant differences between specified and measured airflow values for some duct runs. Thermal imaging results performed on the house when operating in both heating and cooling modes validated historic stratification issues of ceiling supply registers with high supply air temperatures. Long-term monitoring results will be detailed in a future report.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2017

Indirect measurements of slow modes in the Hermean magnetosphere

The aim of this study is to simulate slow mode structures in the Hermean magnetosphere. We use the single fluid MHD code PLUTO and a multipolar expansion of the Northward displaced Hermean magnetic field, to perform simulations mimicking the solar wind and interplanetary magnetic field conditions during 08/09/2011 MESSERGER orbit, susceptible to cross a slow mode structure [F. Pantellini et al, PSS, 112, 1 (2015) and J. Varela et al, PSS, 125, 80 (2016)]. The simulation shows the presence of a slow mode structure nearby the North pole. The satellite trajectory crosses the slow mode compressional front, although the region with the largest pressure and magnetic field lines bending is located Northernmost. The comparison with the MESSENGER data shows a decrease of the magnetic field module in the trajectory region where the slow mode structure is predicted by the simulation. Further analysis of the plasma density measurements by MESSENGER Fast Imaging Plasma Spectrometer are required to confirm the presence of the slow mode structure.

Filippo, Pantellini↗

DOE OSTI2019

High temperature oxidation and microstructural evolution of cold spray chromium coatings on Zircaloy-4 in steam environments

Oxidation kinetics and microstructural evolution of cold sprayed Cr coatings on Zircaloy-4 at 1130–1310 °C in flowing steam at atmospheric pressure have been studied. The study is aimed at understanding of the response of Cr coated Zr-alloy under a steam environment and high temperatures pertinent to design basis accidents (DBAs) and beyond design basis accidents (BDBAs) in light water reactors (LWRs). Surface morphology, microstructure, and phases of post-oxidation test samples were characterized using Scanning Electron Microscopy (SEM), x-ray diffraction (XRD), and Scanning Transmission Electron Microscopy (STEM). Growth kinetics of the Cr-oxide scale and interdiffusion layers between the Cr coating and the Zr-alloy substrate were quantified from cross-sectional SEM images. Cross-sectional analysis showed that the Cr coatings offered a 50 times reduction in oxidation rate over bare Zircaloy-4 in a 1310 °C steam environment. Oxidation kinetics at 1130 °C followed parabolic law (i.e., n ~ 0.5 in power law kinetics) but at 1230 °C and 1310 °C the value of n was suppressed to below 0.5 possibly due to the volatilization of Cr species at the two highest temperatures. Interdiffusion at coating/substrate interface resulted in formation of a brittle Cr 2 Zr or Zr(Fe, Cr) 2 intermetallic compound layer on the order of micrometer in thickness and scattered Cr-rich precipitates were observed well below the interface within the Zr-alloy substrate after cooling. These experimental results could provide data to LWR system simulation codes for better estimation of coping time in the event of accidents.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Interplay between beam-driven chirping modes and plasma confinement transitions in spherical tokamak ST40

Experiments on the high field spherical tokamak ST40 have led to the recent observation of interplay between beam-driven modes of sweeping frequency (chirping modes) and transitions to the enhanced global confinement regime (H-mode) and back to the low confinement regime (L-mode). The H-modes of plasma confinement are identified from decreased intensity of D α signal and from clear distinctions in the edge gradients of the visible plasma boundary (observed as a sharp plasma edge in camera images). The beam-driven chirping modes, identified as ideal magnetohydrodynamics beta-induced Alfvén acoustic eigenmodes modes, are observed in Mirnov coil signals, interferometry, and soft x-ray diagnostics. A moderate amplitude ‘primer’ chirping mode usually precedes an H–L transition. This is followed by a ‘dominant’ chirping mode with higher amplitude during the L-mode. The L–H transition back to the improved confinement occurs on a longer time scale of tens of ms, consistent with the slowing down time scale of fast beam ions. A dramatic decrease in toroidal plasma rotation is systematically observed associated with chirping modes sweeping down to zero frequency. Resonance maps built for the beam-driven chirping modes with the ASCOT (accelerated simulation of charged particle orbits in torodoial devices) code show that the resonant beam ions have orbits near the trapped-passing boundary. The ASCOT modelling assesses how losses of the resonant fast ions caused by the chirping modes with high enough amplitude modify the torque, potentially affecting the plasma rotation.

magnetic confinement fusion↗

DOE OSTI2018

Sandia Optical Fringe Analysis Slope Tool (SOFAST) v. 2018.11.20

SOFAST is used to characterize the surface slope of reflective mirrors for solar applications. SOFAST uses a large monitor or projection screen to display fringe patterns, and a machine vision camera to image the reflection of these patterns in the subject mirror. From these images, a detailed map of surface normals can be generated and compared to design or fitted mirror shapes. SOFAST uses standard Fringe Reflection (Deflectometry) approaches to measure the mirror surface normals. SOFAST uses an extrinsic analysis of key points on the facet to locate the camera and monitor relative to the facet coordinate system. It then refines this position based on the measured surface slope and integrated shape of the mirror facet. The facet is placed into a reference frame such that key points on the facet match the design facet in orientation and position. This is key to evaluating a facet as suitable for a specific solar application. SOFAST reports the measurements of the facet as detailed surface normal location in a format suitable for ray tracing optical analysis codes. SOFAST also reports summary information as to the facet fitted shape (monomial) and error parameters. Useful plots of the error distribution are also presented.

Yellowhair, Julius↗

DOE OSTI2015

Helium Loop Cooling Channel Hydraulic Characterization

New methods for generating 99 Mo are being explored in an effort to eliminate proliferation issues and provide a domestic supply of 99m Tc for medical imaging. Electron accelerating technology is used by sending an electron beam through a series of 100 Mo targets. During this process a large amount of heat is created, which directly affects the operating temperature set for the system. In order to maintain the required temperature range, helium gas is used to serve as a cooling agent that flows through narrow channels between the target disks. Currently we are tailoring the cooling channel entrance and exits to decrease the pressure drop through the targets. Currently all hardware has be procured and manufactured to conduct flow measurements and visualization via solid particle seeder. Pressure drop will be studied as a function of mass flow and diffuser angle. The results from these experiments will help in determining target cooling geometry and validate CFD code results.

42 ENGINEERING↗

DOE OSTI2015

SOFAST: Sandia Optical Fringe Analysis Slope Tool

SOFAST is used to characterize the surface slope of reflective mirrors for solar applications. SOFAST uses a large monitor or projections screen to display fringe patterns, and a machine vision camera to image the reflection of these patterns in the subject mirror. From these images, a detailed map of surface normals can be generated and compared to design or fitted mirror shapes. SOFAST uses standard Fringe Reflection (Deflectometry) approaches to measure the mirror surface normals. SOFAST uses an extrinsic analysis of key points on the facet to locate the camera and monitor relative to the facet coordinate system. It then refines this position based on the measured surface slope and integrated shape of the mirror facet. The facet is placed into a reference frame such that key points on the facet match the design facet in orientation and position. This is key to evaluating a facet as suitable for a specific solar application. SOFAST reports the measurements of the facet as detailed surface normal location in a format suitable for ray tracing optical analysis codes. SOFAST also reports summary information as to the facet fitted shape (monomial) and error parameters. Useful plots of the error distribution are also presented.

Andraka, Charles E.↗

DOE OSTI2012

A statistically defined anthropomorphic software breast phantom

Purpose: Digital anthropomorphic breast phantoms have emerged in the past decade because of recent advances in 3D breast x-ray imaging techniques. Computer phantoms in the literature have incorporated power-law noise to represent glandular tissue and branching structures to represent linear components such as ducts. When power-law noise is added to those phantoms in one piece, the simulated fibroglandular tissue is distributed randomly throughout the breast, resulting in dense tissue placement that may not be observed in a real breast. The authors describe a method for enhancing an existing digital anthropomorphic breast phantom by adding binarized power-law noise to a limited area of the breast. Methods: Phantoms with (0.5 mm){sup 3} voxel size were generated using software developed by Bakic et al. Between 0% and 40% of adipose compartments in each phantom were replaced with binarized power-law noise ({beta} = 3.0) ranging from 0.1 to 0.6 volumetric glandular fraction. The phantoms were compressed to 7.5 cm thickness, then blurred using a 3 Multiplication-Sign 3 boxcar kernel and up-sampled to (0.1 mm){sup 3} voxel size using trilinear interpolation. Following interpolation, the phantoms were adjusted for volumetric glandular fraction using global thresholding. Monoenergetic phantom projections were created, including quantum noise and simulated detector blur. Texture was quantified in the simulated projections using power-spectrum analysis to estimate the power-law exponent {beta} from 25.6 Multiplication-Sign 25.6 mm{sup 2} regions of interest. Results: Phantoms were generated with total volumetric glandular fraction ranging from 3% to 24%. Values for {beta} (averaged per projection view) were found to be between 2.67 and 3.73. Thus, the range of textures of the simulated breasts covers the textures observed in clinical images. Conclusions: Using these new techniques, digital anthropomorphic breast phantoms can be generated with a variety of glandular fractions and patterns. {beta} values for this new phantom are comparable with published values for breast tissue in x-ray projection modalities. The combination of conspicuous linear structures and binarized power-law noise added to a limited area of the phantom qualitatively improves its realism.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

SU-F-T-516: Effects of Inter-Fraction Organ Displacement/deformation On the Delivered Doses to the Heart, Esophagus, and Lungs in Patients Receiving Thoracic Radiotherapy

Purpose: To quantitatively assess the effects of inter-fraction changes in organ shape and location on the delivered dose distribution to the organs at risk (OAR) in lung cancer patients. Methods: This study analyzes treatment data of 10 patients, who were treated to 60Gy in 30 fractions. In each fraction a cone beam CT (CBCT) was acquired. Each CBCT was registered with the planning CT using deformable registration tools within MIM Software. The daily setup shifts were used to translate the planned dose distribution on the deformed planning CT. The structures of lungs, esophagus and heart were re-delineated by a physician on each CBCT. The doses delivered to each OAR, reflecting changes in the position and shape variations, were recomputed. Resultant daily dose volume histograms (DVHs) for OARs were computed and compared to those from the planning CT. Results: Based on the findings of two patients and 24 CBCTs analyzed so far, higher doses are delivered to the lungs and esophagus compared to the treatment plan. The dose differences per fraction between the delivered doses and those in the treatment plan are: for patient 1, lung mean dose = 5.3±1.3cGy and esophagus mean dose = 3.4±3.5cGy. For patient 2, lung mean dose = 12.0±3.9cGy and esophagus mean dose = 34.2±7.5cGy. Regarding the maximum dose to heart, the results varied (−18.9±22.0cGy for patient1 and 53.0±62.2cGy for patient2). Conclusion: The dosimetric effects of inter-fractional anatomical variations could be estimated using deformable image registration and manual organ segmentation for each CBCT. A considerable dose distribution variation between fractions was observed for the OARs. These changes are currently not taken into account while treating the patients and these may explain cases with severe side effects even when the treatment plan looks satisfactory. These results suggest the need for automated daily dose tracking and accumulation.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2014

Hyperspectral Sensors Final Report CRADA No. TC02173.0

This was a collaborative effort between Lawrence Livermore National Security, LLC as manager and operator of Lawrence Livermore National Laboratory (LLNL) and Science Applications International Corporation (SAIC), National Security Space Operations/SRBU, to develop long-wave infrared (LWIR) hyperspectral imaging (HSI) sensors for airborne and potentially ground and space, platforms. The overarching purpose of the CRADA project was to facilitate the transfer of the necessary technology from LLNL to SAIC thereby allowing the U.S. Government to procure LWIR HSI sensors from this company. The purpose of the initial activities involving BIG GREEN sensors was to meet the customer’s delivery schedule for all of the BIG GREEN sensors. BIG GREEN is a complex and sophisticated instrument, and the requisite technology must be imparted to SAIC to allow them to successfully produce the instrument. It was expected that in the first few months of joint work, sensor models (performance, structural, thermal and optical) would be provided to SAIC. Additionally, existing drawings, assembly procedures, test procedures, test equipment definitions, software descriptions and code, and complete Bills of Material (BOM) would be provided. SAIC would analyze the models, purchase material per the BOM, create its own set of production drawings, models and procedures, and identify areas of concern or where further information transfer is required.

47 OTHER INSTRUMENTATION↗

DOE OSTI2016

ON SHOCKS DRIVEN BY HIGH-MASS PLANETS IN RADIATIVELY INEFFICIENT DISKS. II. THREE-DIMENSIONAL GLOBAL DISK SIMULATIONS

Recent high-resolution, near-infrared images of protoplanetary disks have shown that these disks often present spiral features. Spiral arms are among the structures predicted by models of disk–planet interaction and thus it is tempting to suspect that planetary perturbers are responsible for these signatures. However, such interpretation is not free of problems. The observed spirals have large pitch angles, and in at least one case (HD 100546) it appears effectively unpolarized, implying thermal emission of the order of 1000 K (465 ± 40 K at closer inspection). We have recently shown in two-dimensional models that shock dissipation in the supersonic wake of high-mass planets can lead to significant heating if the disk is sufficiently adiabatic. Here we extend this analysis to three dimensions in thermodynamically evolving disks. We use the Pencil Code in spherical coordinates for our models, with a prescription for thermal cooling based on the optical depth of the local vertical gas column. We use a 5M{sub J} planet, and show that shocks in the region around the planet where the Lindblad resonances occur heat the gas to substantially higher temperatures than the ambient gas. The gas is accelerated vertically away from the midplane to form shock bores, and the gas falling back toward the midplane breaks up into a turbulent surf. This turbulence, although localized, has high α values, reaching 0.05 in the inner Lindblad resonance, and 0.1 in the outer one. We find evidence that the disk regions heated up by the shocks become superadiabatic, generating convection far from the planet’s orbit.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

Bringing Fenton Hill into the Digital Age: Data Conversion in Support of the Geothermal Technologies Office Code Comparison Study Challenge Problems

The Geothermal Technologies Office Code Comparison Study (GTO-CCS) was established by the U.S. Department of Energy to facilitate collaboration among members of the geothermal modeling community and to evaluate and improve upon the ability of existing codes to simulate thermal, hydrological, mechanical, and chemical processes associated with complex enhanced geothermal systems (EGS). The first stage of the project, which has been completed, involved comparing simulations for seven benchmark problems that were primarily designed using well-prescribed, simplified data sets. In the second stage, the participating teams are tackling two challenge problems based on the EGS research conducted in hot dry rock (HDR) at Fenton Hill, near Los Alamos, New Mexico. The Fenton Hill project, conducted by Los Alamos National Laboratory (LANL) from 1970 to 1995, was the world’s first HDR demonstration project. One of the criteria for selecting this experiment as the basis for the challenge problems was the amount and availability of data for generating model inputs. The Fenton Hill HDR system consisted of two reservoirs – an earlier Phase I reservoir tested from 1974 to 1981 and a deeper Phase II reservoir tested from 1980 to 1995. Detailed accounts of both phases of the HDR project have been presented in a number of books and reports, including a recently published summary of the lessons learned and a final report with a chronological description of the Fenton Hill project, prepared by LANL. Project documents and records have been archived and made public through the National Geothermal Data System (NGDS). Some of the data acquired from Phase II are available in electronic format readable on modern computers. These include the microseismic data from some of the important experiments (e.g. the massive hydraulic fracturing test conducted in 1983) and the injection/production wellhead data from the circulation tests conducted between 1992-1995. However, much of the data collected during the project, while publicly available, currently only exist in the form of tables or graphs within scanned documents. Therefore, in support of the GTO-CCS, the data needed for developing simulation inputs are being compiled and converted to platform independent, open readable formats so that all participating teams will have access to the same electronic data set. In some cases this requires conversion using optical character recognition, digitizing existing images, and generating the appropriate metadata from project documents. The GTO-Velo knowledge management framework, developed by Pacific Northwest National Laboratory (PNNL), was used for the benchmark problem stage of the comparison study and will also be used as the data repository for the challenge problem data sets. It is staggering and impractical to convert all published data for the Fenton Hill site, so the focus is on data that supports simulations for the three topical areas defined by the study for the challenge problems: 1) reservoir creation/stimulation, 2) reactive and passive transport, and 3) thermal recovery. Conversion of these data provide value not only to GTO-CCS participants, but also to members of the geothermal community at large who may be interested in revisiting the Fenton Hill experiment in the future.

Code Comparison Study↗