Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “filled”

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 1,315 records · Page 73

Space-tiled colloidal crystals from DNA-forced shape-complementary polyhedra pairing

Generating space-filling arrangements of most discrete polyhedra nanostructures of the same shape is not possible. However, if the appropriate individual building blocks are selected (e.g., cubes), or multiple shapes of the appropriate dimensions are matched (e.g., octahedra and tetrahedra) and their pairing interactions are subsequently forced, space-filled architectures may be possible. With flexible molecular ligands (polyethylene glycol–modified DNA), the shape of a polyhedral nanoparticle can be deliberately altered and used to realize geometries that favor space tessellation. In this work, 10 new colloidal crystals were synthesized from DNA-modified nanocrystal building blocks that differed in shapes and sizes, designed to form space-filling architectures with micron-scale dimensions. The insights and capabilities provided by this new strategy substantially expand the scope of colloidal crystals possible and provide an expanded tool kit for researchers interested in designing metamaterials.

Science & Technology - Other Topics↗

Rapid clay precipitation in explosion-induced fractures

Fractures within the earth control rock strength and fluid flow, but their dynamic nature is not well understood. As part of a series of underground chemical explosions in granite in Nevada, we collected and analyzed microfracture density data sets prior to, and following, individual explosions. Our work shows an ~4-fold increase in both open and filled microfractures following the explosions. Based on the timing of core retrieval, filling of some new fractures occurs in as little as 6 wk after fracture opening under shallow (<100 m) crustal conditions. These results suggest that near-surface fractures may fill quite rapidly, potentially changing permeability on time scales relevant to oil, gas, and geothermal energy production; carbon sequestration; seismic cycles; and radionuclide migration from nuclear waste storage and underground nuclear explosions.

58 GEOSCIENCES↗

Ordered and disordered skyrmion states on a square substrate

We examine the ordering of skyrmions interacting with a square substrate created from a modulation of anisotropy using atomistic simulations. We consider fillings of ƒ = 1.5, 2.0, 2.5, 3.0, and 3.5 skyrmions per potential minimum as a function of magnetic field and sample size. For a filling of ƒ = 2.0, we find various dimer orderings, such as tilted dimer states, as well as antiferromagnetic ordering that is similar to the colloidal dimer ordering seen on square substrates. Furthermore, the ability of the skyrmions to change shape or annihilate produces additional states that do not occur in the colloidal systems. For certain parameters at ƒ = 2.0, half of the skyrmions can annihilate to form a square lattice, or a superlattice of trimers and monomers containing skyrmions of different sizes can form. At lower fields, ordered stretched skyrmion states can appear, and for zero field, there can be ordered stripe states. For ƒ = 3.0, we find ferromagnetic ordered trimers, tilted lattices, columnar lattices, and stretched phases. For fillings of ƒ = 1.5 and 2.5, we find bipartite lattices, different ordered and disordered states, and several extended disordered regions produced by frustration effects.

36 MATERIALS SCIENCE↗

Formation of wide-blocky calcite veins by extreme growth competition

Liassic limestones on the coast of Somerset in the UK contain dense arrays of calcite microveins with a common, but poorly understood microstructure, characterized by laterally wide crystals that form bridges across the vein. We investigated the mechanisms of formation and evolution of these ‘wide-blocky’ vein microstructures using a combination of high-resolution analytical methods, including virtual petrography, optical cathodoluminescence and scanning electron microscopy techniques (e.g. energy-dispersive X-ray spectrometry, back-scattered electron imaging, cathodoluminescence and electron back-scattered diffraction), laboratory experiments and multiphase field modelling. Our results indicate that the studied veins formed in open, fluid-filled fractures, each in a single opening and sealing episode. As shown by the optical and electron back-scattered diffraction images, the vein crystals grew epitaxially on grains of the wall rock and we hypothesize that their growth rates differed depending on whether the crystals were on a wall rock grain substrate that fractured intergranularly (slow growth rates) or transgranularly (rapid growth rates). Our multiphase field models support this hypothesis, showing that wide, blocky crystals only form where there are significant differences in the growth rate and are dependent on the type of seed grain. These results provide strong evidence for extreme growth competition, a process that we propose controls vein-filling in many micritic carbonate reservoirs, as well as demonstrate that the characteristics of the fracture wall can affect the filling processes in syntaxial veins. Supplementary material: The description and images of the studied thin sections are available at https://doi.org/10.6084/m9.figshare.c.5172371 . High-resolution optical microscopy mosaics (under plane-polarized- and crossed polarized light) of the thin section collection in PetroScan file format are available on request from the authors.

Geology↗

Fracture, Dissolution, and Cementation Events in Ordovician Carbonate Reservoirs, Tarim Basin, NW China

Ordovician carbonate rocks of the Yijianfang Formation in the Tabei Uplift, Tarim Basin, contain deeply buried (>6000 m), highly productive oil and gas reservoirs associated with large cavities (>10 m). Previous workers inferred that large cavities are paleocaves (paleokarst) formed near the surface and subsequently buried. Alternately, caves may have formed by dissolution at depth along faults. Using 227 samples from 16 cores, we document textures and cement compositions bearing on cavity histories with petrographic, high-resolution scanning electron microscopy (SEM), isotopic, and fluid inclusion microthermometric observations. Results show that dissolution occurred at depth and was caused by (1) acidic fluids derived from Middle-Late Silurian and/or Devonian-Permian hydrocarbon generation and maturation, (2) high-temperature fluids, of which some were associated with Late Permian igneous activity, and (3) Mg-rich fluids that accompanied Jurassic-Cretaceous deformation and the formation of partially open fractures and stylobreccias (fault breccias). The relative paragenetic sequence of the structure-related diagenesis suggests seven stages of fracturing, dissolution, and cementation. Mottle fabrics in the Yijianfang Formation contain argillaceous carbonate-rich silt and are bioturbation features formed within the marine environment. Those mottled fabrics differ from clearly karstic features in the overlying Lianglitage Formation, which formed by near-surface dissolution and subsequent infilling of cavities by allochthonous sediment. Mottle fabrics are crosscut by compacted fractures filled with phreatic-vadose marine cements and followed by subsequent generations of cement-filled fractures and vugs indicating that some fractures and vugs became cement filled prior to later dissolution events. Calcite cements in fractures and vugs show progressively depleted values of δ18O documenting cement precipitation within the shallow (~220 m), intermediate (~625 m), and deep (~2000 m) diagenetic environments. Deep (mesogenetic) dissolution associated with fractures is therefore the principal source of the high porosity-permeability in the reservoir, consistent with other pieces of evidence for cavities localized near faults.

58 GEOSCIENCES↗

Compositional simulation of hydrocarbon recovery from oil shale reservoirs with diverse initial saturations of fluid phases by various thermal processes

We have studied the hydrocarbon production from oil shale reservoirs filled with diverse initial saturations of fluid phases by implementing numerical simulations of various thermal in-situ upgrading processes. We use our in-house fully functional, fully implicit, and non-isothermal simulator, which describes the in-situ upgrading processes and hydrocarbon recovery by multiphase-multicomponent systems. We have conducted two sets of simulation cases—five-spot well pattern problems and Shell In-situ Conversion Process (ICP) problems. In the five-spot well pattern problems, we have analyzed the effects of initial fluid phase that fills the single-phase reservoir and thermal processes by four cases—electrical heating in the single-phase-aqueous reservoir, electrical heating in the single-phase-gaseous reservoir, hot water injection in the single-phase-aqueous reservoir, and hot CO 2 -injection in the single-phase-gaseous reservoir. In the ICP problems, we have analyzed the effects of initial saturations of fluid phases that fill two-phase-aqueous-and-gaseous reservoir by three cases—initial aqueous phase saturations of 0.16, 0.44, and 0.72. Through the simulation cases, system response and production behavior including temperature profile, kerogen fraction profile, evolution of effective porosity and absolute permeability, phase production, and product selectivity are analyzed. In the five-spot well pattern problems, it is found that the hot water injection in the aqueous phase reservoir shows the highest total hydrocarbon production, but also shows the highest water-oil-mass-ratio. Productions of phases and components show very different behavior in the cases of electrical heating in the aqueous phase reservoir and the gaseous phase reservoir. In the ICP problems, it is found that the speed of kerogen decomposition is almost identical in the cases, but the production behavior of phases and components is very different. It is found that more liquid organic phase has been produced in the case with the higher initial saturation of aqueous phase by the less production of gaseous phase.

04 OIL SHALES AND TAR SANDS↗

An effective-medium model for P-wave velocities of saturated, unconsolidated saline permafrost

To better understand the relationship between P-wave velocities and ice content in saturated, unconsolidated saline permafrost, we constructed an effective-medium model based upon ultrasonic P-wave data that were obtained from earlier laboratory studies. The model uses a two-end-member mixing approach in which an ice-filled, fully frozen end member and a water-filled, fully unfrozen end member are mixed together to form the effective medium of partially frozen sediments. This mixing approach has two key advantages: (1) It does not require parameter tuning of the mixing ratios, and (2) it inherently assumes mixed pore-scale distributions of ice that consist of frame-strengthening (i.e., cementing and/or loadbearing) ice and pore-filling ice. The model-predicted P-wave velocities agree well with our laboratory data, demonstrating the effectiveness of the model for quantitatively inferring ice content from P-wave velocities. The modeling workflow is simple and is largely free of calibration parameters -attributes that ease its application in interpreting field data sets.

54 ENVIRONMENTAL SCIENCES↗

Development of the Los Alamos National Laboratory Cryogenic Pressure Loader

Targets for inertial fusion research and ignition at OMEGA, the National Ignition Facility, LMJ, and future facilities rely on beta-radiation-driven layering of spherical cryogenic DT ice layers contained within plastic or metal shells. Plastic shells will be permeation filled at room temperature then cooled to cryogenic temperatures before removal of the overpressure. The cryogenic pressure loader (CPL) was recently developed at Los Alamos National Laboratory as a testbed for studying the filling and layering of plastic target shells with DT. A technical description of the CPL is provided. The CPL consists of a cryostat, which contains a high-pressure permeation cell, and has optical access for investigating beta layering. The cryostat is housed within a tritium glovebox that contains manifolds for supplying high-pressure DT. The CPL shares some design elements with the cryogenic target handling system at the OMEGA facility to allow testing of tritium issues related to that system. The CPL has the capability to fill plastic targets by permeation to pressures up to 100 MPa and to cool them to 15 K. The CPL will accommodate a range of targets and may be modified for future experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Washington Geothermal Play Fairway Analysis Heat, Permeability, and Fracture Model Data

This submission contains raster and vector data for the entire state of Washington, with specific emphasis on the three geothermal play fairway sites: Mount St. Helens seismic zone (MSHSZ), Wind River valley (WRV), and Mount Baker (MB). Data are provided for 3 major geothermal models: heat, permeability, and fluid-filled fractures, and an additional infrastructure model. Both of the permeability and fluid-filled-fracture models are produced at 200 m and at 2 km depths; the heat model is only produced at the 200 m depth. Values are provided for both model favorability and model confidence. A combined model at 200m and 2 km depths is provided for favorability, confidence, and exploration risk. Raster data are provided in GeoTiff format and have a statewide coverage. Cell size is 104.355 ft; file type is unsigned 8-bit integer (0-255); 0 represents no favorability or confidence; 255 represents maximum favorability or confidence. The NAD83(HARN)/Washington South (ftUS) projection is used (EPSG:2927). Vector data are provided in shapefile or comma-delimited text file formats. Geographic coordinates, where provided, are in WGS84. A readme file accompanies each folder and provides an overview and description of the enclosed data. The heat model combines 5 intermediate raster layers (which are included in the download package): temperature gradient wells, young volcanic vents, hot springs, young intrusive volcanic rocks, and geothermometry. The permeability model combines 8 intermediate raster layers: density of mapped faults, 2D dilation tendency of mapped faults, 2D slip tendency of mapped faults, seismicity, 3D dilation tendency, 3D slip tendency, 3D maximum coulomb shear stress, and 3D slip gradients. The fluid-filled fracture model combines up to 4 intermediate rasters: resistivity from magneto-telluric 3D inversions, seismicity, Vp/Vs anomalies from passive seismic tomography, and Vs anomalies from ambient-noise tomography. A statewide infrastructure model is also provided that formalizes land-use constraints and restrictions relevant for geothermal prospecting and development. This model combines 10 intermediate rasters: areas off limits to drilling, existing or proposed geothermal leases, DNR-owned land, land-use restrictions along the Columbia River Gorge, areas inundated by water, availability of potential process water, proximity to existing roads, proximity to transmission lines, distance from urban areas, and snow-related elevation restrictions. Supporting vector data for the development of each raster layer is provided. For details on the areas of interest and modeling process please see the 'WA_State_Play_Fairway_Phase_2_Technical_Report' in the download package.

15 GEOTHERMAL ENERGY↗

Data From: "Warming and snow loss increase reliance on old groundwater in a Colorado River headwater"

This repository contains the data and code associated with the paper titled "Warming and snow loss increase reliance on old groundwater in a Colorado River headwater," published in Nature Geoscience, 2026. This study seeks to answer how various ages of groundwater interact with mountainous streamflow in mountainous headwaters such as the East River. It includes various model-data processing scripts, primarily for ParFlow-CLM analysis of simulated water years 2015-2021, and two numerical warming experiments (+2.5 and +4.0 degrees C), including run scripts, forcing scripts, and post-processing, as well as comparison to observation datasets, detailed below. This data requires the use of R (.r, .rmd), Python (.py), Jupyter Notebook or Jupyter Lab (.ipynb), ParFLOW-CLM, EcoSLIM. Further information on the use of all file formats mentioned below (e.g. .tff. .nc) are provided within the associated scripts and directory where the files are located. Contents & Usage ASO/: ​​Contains the bash and python scripts used to convert airborne snow observatory (ASO) data (ASO, 2023) in various data formats (georeferenced tiff file, NetCDF, UTM, and to latitude/longitude) then regrided to the ParFlow equivalent grid. Output data are in regrid_regll_data.zip and subsequently visualized and analyzed in plot_and_compare.py for Supplementary Figures A14 and A15. The wksht_ASO_comparison.xlsx spreadsheet is used to calculate the data for Supplementary Figure A16. EcoSLIM/: Contains the scripts and input files to run the EcoSLIM particle tracking simulations (/run_scripts) and the post-processing python script (/plot_scripts/eco_agedist_plots.ipynb). Jasechko et al./: Contains the jupyter notebook (Extract_Elevation.ipynb) to determine the outlet elevations of the 260 watersheds used in Jasechko et al. (2016), and the corresponding table, Table_S1_Watersheds_alt.csv. Used to create Supplementary Information Figure A2. PLM_Wells/: Contains the QA/QC-ed groundwater level time series of the PLM-1 and PLM-6 Monitoring Wells from Faybishenko et al. (2023), reformatted to water years used for Supplementary Figures A19 and and A20. ParFlow/: Contains the input files and run scripts to run ParFlow-CLM (/run_scripts), the python and tool command language (Tcl) scripts to create and distribute the ParFlow forcing simulation files (/forcing), and various scripts and intermediary files to analyze the model outputs (/post_process). SQUIRE/: Contains the processing scripts and intermediary files for the Surface QUantitatIve pRecipitation Estimation (SQUIRE) data (Grover, 2023) used to generate Supplementary Figure A18. USGS_Streamflow/: Contains the raw and gap-filled United States Geological Survey streamflow data (U.S. Geological Survey, 2026) used at the Almont station (site number 09112500). Gap-filling is performed in the R script with data from the Taylor station (site number 09110000). (/USGS_09112500_EAST_RIVER_AT_ALMONT_GAP_FILLED/code_almont_streamflow_gap_fill.Rmd). discharge/: Contains the gap-filled discharge data at the Watershed Function SFA East River pumphouse site (Newcomer et al., 2022) used to generate Supplementary Figure A13 and to compute hourly Nash-Sutcliffe model efficiency coefficients (NSE) in Table A4. snotel_and_flux_tower/: Contains the snow telemetry data (U.S. Department of Agriculture, 2024) from the Butte (site ID 380) and Schofield (site ID 737) stations, reformatted by water year, accessed with the snotelr R package. Used to create Supplementary Figure A17. Also contains the flux tower observational data (FluxTower_Pumphouse_ESS-DIVE.ET_only.h.txt) from Ryken et al. (2022) and sap flux transpiration data (MaxB_Transpiration_5Sites.daily_sums.h.txt) from Ryken (2021), used to create Supplementary Figures A22 and A23, respectively. Raw EcoSLIM model outputs are in excess of 24TB, and are stored on National Energy Research Scientific Computing Center (NERSC) and publicly available via the external link provided in the paper.

atmospheric warming↗

Development of a MODIS-Derived Surface Albedo Data Set: An Improved Model Input for Processing the NSRDB

A significant source of bias in the transposition of global horizontal irradiance to plane-of-array (POA) irradiance arises from inaccurate estimations of surface albedo. The current physics-based model used to produce the National Solar Radiation Database (NSRDB) relies on model estimations of surface albedo from a reanalysis climatalogy produced at relatively coarse spatial resolution compared to that of the NSRDB. As an input to spectral decomposition and transposition models, more accurate surface albedo data from remotely sensed imagery at finer spatial resolutions would improve accuracy in the final product. The National Renewable Energy Laboratory (NREL) developed an improved white-sky (bi-hemispherical reflectance) broadband (0.3-5.0 ..mu..m) surface albedo data set for processing the NSRDB from two existing data sets: a gap-filled albedo product and a daily snow cover product. The Moderate Resolution Imaging Spectroradiometer (MODIS) sensors onboard the Terra and Aqua satellites have provided high-quality measurements of surface albedo at 30 arc-second spatial resolution and 8-day temporal resolution since 2001. The high spatial and temporal resolutions and the temporal coverage of the MODIS sensor will allow for improved modeling of POA irradiance in the NSRDB. However, cloud and snow cover interfere with MODIS observations of ground surface albedo, and thus they require post-processing. The MODIS production team applied a gap-filling methodology to interpolate observations obscured by clouds or ephemeral snow. This approach filled pixels with ephemeral snow cover because the 8-day temporal resolution is too coarse to accurately capture the variability of snow cover and its impact on albedo estimates. However, for this project, accurate representation of daily snow cover change is important in producing the NSRDB. Therefore, NREL also used the Integrated Multisensor Snow and Ice Mapping System data set, which provides daily snow cover observations of the Northern Hemisphere for the temporal extent of the NSRDB (1998-2015). We provide a review of validation studies conducted on these two products and describe the methodology developed by NREL to remap the data products to the NSRDB grid and integrate them into a seamless daily data set.

14 SOLAR ENERGY↗

Interim Report on Mixing During the Casting of LEU-10Mo Plates in the Triple Plate Molds

LEU-10%Mo castings are commonly produced by down blending unalloyed HEU with a DU-12.7%Mo master-alloy. This work uses process modeling to provide insight into the mixing of the unalloyed uranium and U-Mo master alloy during melting and mold filling of a triple plate casting. Two different sets of situations are considered: (1) mixing during mold filling from a compositionally stratified crucible and (2) convective mixing of a compositionally stratified crucible during mold heating. The mold filling simulations are performed on the original Y-12 triple plate mold and the horizontal triple plate mold.

36 MATERIALS SCIENCE↗

LDRD Final Report: Advanced Hohlraum Concepts

Indirect drive inertial confinement fusion (ICF) experiments to date have mostly used cylindrical, laser-heated, gas-filled hohlraums to produce the radiation drive needed to symmetrically implode DT-filled fusion capsules. These hohlraums have generally been unable to produce a symmetric radiation drive through the end of the desired drive pulse, and are plagued with complications due to laser-plasma interactions (LPI) that have made it difficult to predict their performance. In this project we developed several alternate hohlraum concepts. These new hohlraums utilize different hohlraum geometries, radiation shields, and foam materials in an attempt to improve performance relative to cylindrical hohlraums. Each alternate design was optimized using radiation hydrodynamic (RH) design codes to implode a reference DT capsule with a high-density carbon (HDC) ablator. The laser power and energy required to produce the desired time-dependent radiation drive, and the resulting time-dependent radiation symmetry for each new concept were compared to the results for a reference cylindrical hohlraum. Since several of the new designs needed extra laser entrance holes (LEHs), techniques to keep small LEHs open longer, including high-Z foam liners and low-Z wires at the LEH axis, were investigated numerically. Supporting experiments and target fabrication efforts were also done as part of this project. On the Janus laser facility plastic tubes open at one end (halfraums) and filled with SiO 2 or Ta 2 O 5 foam were heated with a single 2w laser. Laser propagation and backscatter were measured. Generally the measured propagation was slower than calculated, and the measured laser backscatter was less than calculated. A comparable, scaled up experiment was designed for the NIF facility and four targets were built. Since low density gold foam was identified as a desirable material for lining the LEH and the hohlraum wall, a technique was developed to produce 550 mg/cc gold foam, and a sample of this material was successfully manufactured.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Initial Computational Study of a New Multi-Hole Hohlraum (the "Midraum")

Existing cylindrical hohlraums with two oppositely positioned laser entrance holes (LEHs) have multiple constraints. Their goal is to produce radiation sources distributed over the sky, as visible from the spherical implosion capsule, with most of the deposition near the zeroes of the fourth Legendre polynomial in cosine of the polar angle. This requires some of the laser light to propagate across the hohlraum to positions near the hohlraum symmetry plane. The ratio of case spherical radius to capsule spherical radius should exceed 3 so that the light doesn’t pass through over-dense ablator plasma. Radiation transport can smooth higher radiation modes. For capsules that demand long pulse lengths, hohlraum walls can blow in and change the position where light is absorbed. This changes the radiation symmetry in a time dependent fashion. This affects both P2 and P4. This wall motion can be reduced by introducing fill gas into the hohlraum. The gas provides back pressure and tamps the wall motion. Adding the fill gas comes at some cost. It leads to increased absorption of laser light along the path. The fill gas adds heat capacity to the system, ultimately requiring more laser energy to meet the radiation flux goals, both in total and particularly in the amount of radiation coming from the vicinity of the capsule waist. Given the existing beam pointing at NIF energy from the outer beams must be transferred into the inner beams. Cross beam energy transport (CBET) is accomplished via a plasma instability. This transfer is not perfectly predictable. In addition, the higher intensity required to make up for the losses along the long path can lead to stimulated backscatter as well as the generation of suprathermal electrons. The inner beams will pass through the plasma ablated from the capsule toward the end of the pulse. Heating this plasma acts as another parasitic loss. In addition, the light passing through the turbulent blow-off can be refracted in unpredictable directions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multiscale Modeling of CO 2 Migration and Trapping in Fractured Reservoirs with Validation by Model Comparison and Real-Site Applications (Final Technical Report)

This report documents the accomplishments achieved during the project titled “Multiscale Modeling of CO 2 Migration and Trapping in Fractured Reservoirs with Validation by Model Comparison and Real-Site Applications” funded by the US Department of Energy, Office of Fossil Energy. The objectives of the project were to develop modeling capabilities for predicting CO 2 and brine migration in fractured reservoirs during geologic carbon storage and to investigate the feasibility of carbon storage in fractured reservoirs using the newly developed modeling capabilities. To achieve these objectives new mass transfer functions were developed to adapt existing dual-continuum approaches – a commonly used approach in hydrocarbon reservoir modeling – to the CO 2 -brine system. In the dual-continuum approach the undisturbed rock matrix and the fractures are modeled as two separate continua, which are coupled through mass exchange of CO 2 and brine between the two continua. Mass transfer functions represent the mass exchange between fractures and the rock matrix, without requiring highly resolved numerical grids of individual fractures, thus allowing large modeling domains relevant to geologic carbon storage questions (tens of meters to hundreds of kilometers). Mass transfer functions for the exchange of both free-phase CO 2 and brine between fractures and rock matrix were developed for both the injection phase and the post-injection phase. During the injection phase the fractures are filled with CO 2 , while the rock matrix is initially filled with brine. Therefore, a mass transfer function was developed that is based on vertical displacement of brine by CO 2 in the rock matrix due to gravitational forces. Once CO 2 injection ceases and CO 2 migrates away from the direct vicinity of the injection site, the fractures will begin to re-fill with brine, while the rock matrix will have high CO 2 saturations. For these conditions a mass transfer functions was developed that is based on spontaneous imbibition, as the pore spaces in the rock matrix are much smaller than in the fractures. For formations with low permeability or high capillary entry pressure of the rock matrix, free-phase CO 2 may not be able to enter the rock matrix. However, free-phase CO 2 at the interface between fractures and rock matrix can dissolve into brine resident in the rock matrix. Therefore, additional mass transfer functions were developed based on aqueous-phase diffusion of CO 2 through the rock matrix. These mass transfer functions are based on existing solutions for diffusion and were expanded to be applicable to a variety of block shapes relevant to dual-continuum modeling. Implementations of the newly developed mass transfer functions showed that they give accurate results when compared to highly resolved models. In addition to developing new transfer functions, a new vertically-integrated dual-continuum modeling framework was developed. This development takes advantage of fast vertical segregation of CO 2 and brine in the fractures due to the density difference between CO 2 and brine and the high permeability of the fractures. Therefore, a vertically-integrated vertical equilibrium model for the fractures is coupled with a vertically-distributed model for the rock matrix. The newly developed approach significantly reduces the computational effort, but is only accurate for fracture permeabilities of 100 mD or higher. Finally, the newly developed models were applied to real and hypothetical geologic carbon storage sites to investigate the feasibility of using fractured reservoirs for carbon storage. The general results were, that storage of CO 2 – both as free-phase and dissolved – is feasible, but that storage capacities are likely lower than in unfractured systems. While the high permeability of the fractures reduces the injection pressure (or permits higher injection rates), fast migration of CO 2 may limit the overall storage volume, because CO 2 may escape laterally with not enough of the rock matrix being swept for effective storage.

20 FOSSIL-FUELED POWER PLANTS↗

Preliminary Analysis of Species Partitioning in the DWPF Melter - Sludge Batch 7b

The goal of the work described in this report is to provide additional off-gas carryover data from the Defense Waste Processing Facility (DWPF) melter to support the estimation of off-gas carryover rates from the Low Activity Waste (LAW) melter at the Hanford Site. Of particular interest is to determine the volatility of technetium-99 (Tc-99), which is a key isotope of concern for disposal and is a significant challenge to retain in glass in the Hanford waste treatment flowsheet. Since the DWPF melter is not equipped with routine offgas monitoring and/or sampling capabilities, the off-gas carryover data had to be extracted using component mass balance of the measured feed and pour stream (PS) glass compositions under time-averaged melter operating conditions over the duration of one canister-filling cycle, typically 20-25 hours. The method was first tested using the Sludge Batch 6 (SB6) PS sample data taken while Canister #3472 was being filled on 12/20/2010, approximately three months after the (argon) glass bubblers were installed. The method was tested again using the SB7a PS sample data taken while Canister #3619 was being filled on 8/5/2011 with the glass bubblers in operation.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Modeling Hydrodynamic Instabilities and Mix in National Ignition Facility Hohlraums

Radiation hydrodynamics modeling of hohlraum experiments on the National Ignition Facility has shown a chronic overprediction of hohlraum performance compared to experiment. Depending on the particular experiment, up to several hundred kilojoules of laser energy must be artificially subtracted from the simulations in order to match experimentally measured quantities like the time of peak neutron production (bang time). Fine-scale hydrodynamic mixing of the high-Z hohlraum wall with the low-Z hohlraum gas, an effect omitted in current hohlraum simulations, has long been hypothesized as possible source of some of this degradation. A preliminary computational assessment using a spherical hohlraum model has been made. It shows that the wall-gas interface is both unstable during certain phases of the experiment and has the potential for turbulent mixing of the gas and wall material, depending on the hohlraum gas fill and pulse shape. The same simulations also indicate that mixing at this interface can delay the simulated bang time by up to 140 ps with high gas fills but appears negligible for low gas fills.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Microphysics Studies of Hohlraum Dynamics (LDRD 17-ERD-060 Final Report)

The goal of this project was demonstrate the feasibility and utility of modeling the effects of kinetic physics in gas-filled hohlraums using high-fidelity simulation tools. We have successfully achieved this goal. In our work, we have found that hybrid particle-in-cell simulations (Thoma et al. 2017, Higginson et al. 2019a) utilizing quasi-neutral, fluid electrons and kinetic ions are powerful tools to understand the kinetic physics of interpenetration and mix relevant to low-gas-fill (near-vacuum) hohlraums. We have benchmarked the code, Chicago, against OMEGA laser experiments designed to mimic conditions found in hohlraums; the code was able to reproduce data from x-ray imaging and Thomson Spectra resolved over both time and space. The experiments highlighted kinetic-ion effects such increased penetration depth that were only able to be resolved with a kinetic model. We have run preliminary simulations with Chicago of full-scale hohlraum experiments that include laser raytracing, equation-of-state and ionization, radiation transport, and accurate experimental geometries. Continuation of this work will predictive hohlraum models at low-gas-fill and identify the better regimes of hohlraum design with the potential for improved coupling and capsule symmetry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗