Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “replication”

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,495 records · Page 83
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

The solubility and speciation of REE phosphate endmembers (CePO 4 and YPO 4 ) in Cl-rich aqueous fluids from 350 to 450 °C and implications for natural systems

The rare earth elements (REE) are important metals used increasingly in advanced technologies. Within the crust, the elements Ce and Y are commonly more abundant compared to other lanthanides and comprise important end-member constituents of REE-bearing minerals. Specifically, Ce is part of the light (L) REE which have larger ionic radii than the heavy (H) REE, which are grouped together with Y. These differences in ionic radius can lead to important physico-chemical trends within the lanthanide group. Despite a recent increase in experimental and thermodynamic data for the REE at high temperature and pressure, there is still a significant lack of these data at supercritical conditions. In this study we conducted batch-type experiments to measure the solubility of REE phosphates (CePO 4 and YPO 4 ) at varying starting pH (1.5–10), and salinity (0.01–1.4 mol/kg NaCl) at 350 and P sat , and from 400 to 450 °C at 700 bar. Results show that the solubility of Ce (33–0.14 ppb) is generally higher than Y (13–0.13 ppb) and that Ce complexes more strongly with both chloride and hydroxyl ligands compared to Y. The solubilities of both REE phosphates are highly pH-dependent and, to a lesser extent, depend on salinity at the studied conditions. The solubility data from this study were implemented into the GEMSFITS program to optimize the thermodynamic properties of Ce and Y hydroxyl and chloride species. The updated standard partial molal Gibbs energies of formation (Δ f G 0 T,P ) are used within the experimental temperature and pressure range to accurately predict the CePO 4 and YPO 4 solubility and Ce and Y speciation behavior. Based on the updated thermodynamic properties we also provide formation constants (log β n Cl,OH ) for Ce and Y hydroxyl and chloride species. Updated thermodynamic properties are applied to model REE-apatite dissolution and REE mobility based on the Pea Ridge iron oxide apatite deposit in Missouri, USA. The apatite dissolution model replicates natural observations including the replacement of monazite and xenotime after apatite and is an example of the utility of the new thermodynamic constants applied to supercritical crustal fluids. Furthermore, the findings of this study advance the predictive capabilities of geochemical models, our understanding of the behavior of individual REE, and permit modeling the overarching fractionation trends between LREE and HREE in supercritical crustal fluids.

58 GEOSCIENCES↗
DOE OSTI2020

Assessing the sensitivity and repeatability of permanganate oxidizable carbon as a soil health metric: An interlab comparison across soils

Soil organic matter is central to the soil health framework. Therefore, reliable indicators of Soil organic matter is central to the soil health framework. Therefore, reliable indicators of changes in soil organic matter are essential to inform land management decisions. Permanganate oxidizable carbon (PDXC), an emerging soil health indicator, has shown promise for being sensitive to soil management. However, strict standardization is required for widespread implementation in research and commercial contexts. Here, we used 36 soils-three from each of the 12 USDA soil orders-to determine the effects of sieve size and soil mass of analysis on PDXC results. Using replicated measurements across 12 labs in the US and the EU (n = 7951 samples), we quantified the relative importance of 1) variation between labs, 2) variation within labs, 3) effect soil mass, and 4) effect of soil sieve size on the repeatability of PDXC. We found a wide range of overall variability in PDXC values across labs (0.03 to 171.8%; mean = 13.4%), and much of this variability was attributable to within-lab variation (median = 6.5%) independently of soil mass or sieve size. Greater soil mass (2.5 g) decreased absolute PDXC values by a mean of 177 mg kg -1 soil and decreased analytical variability by 6.5%. For soils with organic carbon (SOC) >10%, greater soil mass (2.5 g) resulted in more frequent PDXC values above the limit of detection whereas the lower soil mass (0.75 g) resulted in PDXC values below the limit of detection for SOC contents <5%. A finer sieve size increased absolute values of PDXC by 124 mg kg -1 while decreasing the analytical variability by 1.8%. In general, soils with greater SOC contents had lower analytical variability. These results point to potential standardizations of the PDXC protocol that can decrease the variability of the metric. We recommend that the PDXC protocol be standardized to use 2.5 g for soils <10% SOC. Sieve size was a relatively small contributor to analytical variability and therefore we recommend that this decision be tailored to the study purpose. Tradeoffs associated with these standardizations can be mitigated, ultimately providing guidance on how to standardize PDXC for routine analysis.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Soil carbon change in intensive agriculture after 25 years of conservation management

Changes in soil organic carbon (SOC) and nitrogen (SON) are strongly affected by land management but few long-term comparative studies have surveyed changes throughout the whole soil profile. We quantified 25-year SOC and SON changes to 1 m in 10 replicate ecosystems at an Upper Midwest, USA site. We compared four annual cropping systems in maize (Zea mays)-soybean (Glycine max)-winter wheat (Triticum aestivum) rotations, each managed differently (Conventional, No-till, Reduced input, and Biologically based); in three managed perennial systems (hybrid Poplar (Populus × euramericana), Alfalfa (Medicago sativa), and Conifer (Pinus spp.); and in three successional systems (Early, Mid- and Late succession undergoing a gradual change in species composition and structure over time). Both Reduced input and Biologically based systems included winter cover crops. Neither SOC nor SON changed significantly in the Conventional or Late successional systems over 25 years. All other systems gained SOC and SON to different degrees. SOC accrual was fastest in the Early successional system (0.8 ± 0.1 Mg C ha –1 y –1 ) followed by Alfalfa and Conifer (avg. 0.7 ± 0.1 Mg C ha –1 y –1 ), Poplar, Reduced input, and Biologically based systems (avg. 0.4 ± 0.1 Mg C ha –1 y –1 ), and Mid-successional and No-till systems (0.3 and 0.2 Mg C ha –1 y –1 , respectively). Over the most recent 12 years, rates of SOC accrual slowed in all systems except Reduced input and Mid-successional. There was no evidence of SOC loss at depth in any system, including No-till. Rates of SON accrual ranged from 64.7 to 0.8 kg N ha –1 y –1 in the order Alfalfa ≥ Early successional > Reduced input and Biologically based ≥ Poplar > No-till and Conifer > Mid-successional systems. Pyrogenic C levels in the Conventional, Early, and Late successional systems were similar despite 17 years of annual burning in the Early successional system (~ 15 % of SOC to 50 cm, on average, and ~40 % of SOC from 50 to 100 cm). Results underscore the importance of cover crops, perennial crops, and no-till options for sequestering whole profile C in intensively managed croplands.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Effect of pore fluid chemistry on the mechanical behavior of a divalent compacted bentonite, an experimental and constitutive study

Ongoing research in isolating high-level nuclear waste and spent fuel has highlighted compacted bentonite as a suitable material for engineered barrier systems in deep geological repositories due to its extraordinary swelling and retention properties. This research focuses on the chemo-mechanical behavior of compacted bentonite exposed to different pore fluids with different concentrations and loading conditions. The study involves swelling pressure and compressibility experiments along with mineralogy analysis employing X-ray diffraction (XRD) and Cation exchange. The tests were conducted on BCV (a Mg/Ca- bentonite) compacted at a dry density of 1.48 ± .02 Mg/m 3 . An advanced chemical-mechanical constitutive model for unsaturated highly expansive clays was adopted to simulate the material response and better understand its behavior. The model is able to account for the main phenomena at both macro and microstructural levels and the interactions between them. The model successfully replicated experimental observations. The XRD analyses support the macroscopic observation, indicating that salinity impacts crystalline swelling as demonstrated by the reduction of basal spacing from 19.27 Å to 15.68 Å when the osmotic suction increases from 0 MPa to 33 MPa. The results suggested that the osmotic pressure generated by the concentration in the pore fluids promotes a reduction in swelling pressures, swelling strains, and crystalline swelling of clay minerals. Also, it affects the pre-consolidation stress and the compressibility of the compacted samples. In conclusion, it was also observed that both solution type and solution concentration impact the clay swelling pressure.

Chemo-mechanical constitutive model↗
DOE OSTI2025

Surrogate models for development of unconventional shale reservoirs by an integrated numerical approach of hydraulic fracturing, flow and geomechanics, and machine learning

We develop well-completion surrogate models by taking an integrated workflow of hydraulic fracturing, flow, geomechanics, and machine learning simulation. There are three steps in the proposed workflow. First, history-matching processes are conducted with the field data including pumping and production data for characterization. Second, full-physics simulation is performed with various parameters of the field development (e.g., cluster spacing, clusters per stage, pumping rates and times, amount of proppant, and well spacing) to generate multiple simulation results by changing the parameters of the completion design with well-known hydraulic fracturing, reservoir, geomechanics simulators to calculate fracture geometry, reservoir depressurization, induced stress changes. The workflow is demonstrated over a field in the Southern Midland Basin. Here, we take two completion scenarios: a single well case followed by a multi-well case. Finally, a Long Short-Term Memory (LSTM) machine learning algorithm is employed to create surrogate models that can replicate the full-physics simulation results. Furthermore, results show that the trained models applied in the single well and multi-well cases for a particular geological system can provide good accuracy close to those provided by full-physics simulations. Specifically, the site-specific surrogate models can predict fracture parameters (length, height, and surface area) and cumulative production accurately with computational efficiency, suggesting our proposed workflow can be used as a pragmatic tool for expediting the well completion optimization process.

Geomechanics↗
DOE OSTI2017

A plasma deflagration accelerator as a platform for laboratory astrophysics

The replication of astrophysical flows in the laboratory is critical for isolating particular phenomena and dynamics that appear in complex, highly-coupled natural systems. Particularly, plasma jets are observed in astrophysical contexts at a variety of scales, typically at high magnetic Reynolds number and driven by internal currents. In this paper, we present detailed measurements of the plasma parameters within deflagration-produced plasma jets, the scaling of these parameters against both machine operating conditions and the corresponding astrophysical phenomena. Using optical and spectroscopic diagnostics, including Schlieren cinematography, we demonstrate the production of current-driven plasma jets of »100 km/s and magnetic Reynolds numbers of »100, and discuss the dynamics of their acceleration into vacuum. Finally, the results of this study will contribute to the reproduction of various types of astrophysical jets in the laboratory and indicate the ability to further probe active research areas such as jet collimation, stability, and interaction.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Progress toward NIF opacity measurements

The opacities of iron and other mid-Z elements help to regulate the transport of energy in the sun. Recent experiments on the Sandia National Laboratories Z machine have shown large discrepancies between the measured and calculated opacities of iron at certain solar conditions. To replicate these opacity measurements a platform is being developed on the National Ignition Facility to measure the opacities of iron and other elements at the same conditions as in the Z experiments. The NIF platform consists of a hohlraum to heat the opacity sample to the desired conditions, a separate backlighter to radiograph the sample, and a spectrometer to give the spectrally resolved opacity. Not only must the opacity be measured but the temperature and density of the sample must also be accurately determined. This platform has now produced its first iron transmission measurements. Overall, these measurements will be presented along with plans for future measurements and details on how the measurements will be improved.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Shock wave formation in the thermosphere by an earthgrazing fireball: Empirical evidence for volatile-enhanced hydrodynamic shielding

Hydrodynamic shielding is a theoretically well-established but observationally elusive and experimentally difficult-to-replicate phenomenon with implications that extend far beyond meteor physics. Rare earthgrazing meteoroids with infrasound signatures that penetrate to the ground can be used to probe hydrodynamic shielding that leads to strong shock formation at high altitude. Here, we report the first coordinated optical and multi-station infrasound observations of a centimeter-scale earthgrazing fireball that generated sustained cylindrical line shock at thermospheric altitudes near 92 km. The event was recorded by numerous optical stations and three infrasound arrays, allowing trajectory reconstruction, ablation behavior, acoustic source localization, and shock characteristics. Optical observations indicate early mechanical erosion and ablation/evaporation at exceptionally low dynamic pressure, consistent with a cometary or a porous, volatile-bearing CM chondritic object. Independent infrasound detections localize shock generation to multiple points along a 164 km trajectory segment near perigee. Weak-shock modeling yields a consistent blast radius of ∼30 m, implying an acoustic-equivalent source size far exceeding the physical dimensions of the ∼45 g nucleus. We demonstrate that classical gas dynamics and ablation-driven hydrodynamic shielding alone cannot account for these observations under ambient thermospheric conditions. We show that volatile release provides the additional flow-field density enhancement required to amplify hydrodynamic shielding, reduce the effective local Knudsen number, and sustain a shock envelope capable of radiating detectable infrasound. Furthermore, these results demonstrate that small, volatile-rich meteoroids can transiently establish continuum-like flow in rarefied environments.

Astrometry↗
DOE OSTI2020

After the hurricane: Validating a resilience assessment methodology

With increasing utility grid outages in the United States, there is growing interest in assessing risk and developing mitigation strategies to reduce the impact of grid outages. Working with the U.S. Air Force, the U.S. Department of Energy's National Renewable Energy Laboratory (NREL) developed a replicable energy resilience assessment methodology and investment decision tool to: (1) identify and score hazards and vulnerabilities at the site level; (2) analyze risks to energy infrastructure; and (3) identify and prioritize energy resilience investments. This work improves on existing resilience assessment methodologies and tools by combining a bottom-up, all-hazards assessment methodology with top-down geographic information system mapping capabilities to provide an innovative, dynamic tool for identifying and prioritizing actionable solutions. This process combines probabilistic forecasting with an iterative approach for continuously updating and reassessing risks to address temporal dynamism. Relationships among systems are modeled and visualized to estimate the effectiveness of resilience actions across multiple interdependent systems and inform financial priorities through cost-difficulty-impact trade-offs. The approach is validated in a case study at Tyndall Air Force Base (AFB) in Florida, which experienced a Category 5 hurricane in 2018. The risks and mitigation strategies identified pre-hurricane are compared with post-hurricane, realized impacts. The assessment effectively identifies risks and actions to increase site energy resilience, but the methodology can be enhanced though greater consideration of the interdependencies between the energy system and related systems like transportation, communication, and food/water systems, which impact the recovery of the energy system and the base.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2017

Jet stability and wall impingement flow field in a thermal striping experiment

We present velocity and temperature field measurements for a 0.9 x 0.9 x 1.7 m glass tank in which two air jets at Re=10000 mix and impinge upon the lid at ambient temperature and pressure. Flow patterns are characterized across a 350 x 200 mm plane located 3 mm below the lid for two inlet geometries: 1) “extended”, in which inlet channels protrude above the tank base, and 2) “flush”, a flat base without protrusions. This minor geometry variation produced distinct changes in the lid flow field, appearing as three stagnant regions for the extended case and only one for flush. The dichotomy is attributed to system stability characteristics: jets are stable in the extended case and unstable for flush. In a separate set of nonisothermal tests, the impingement temperature field was measured for inlet temperature mismatches of 4 oC and jets near Re=10000. A 50 m-long fiber optic distributed temperature sensor positioned 2 mm below the lid measured at 1350 locations. Like the velocity fields, the temperature fields differ for the two inlet geometries: good thermal mixing for the flush case and subdued mixing for the extended case. Simulations with the spectral element code Nek5000 replicated the observed stability dichotomy, duplicating the number of stagnant regions observed in the experiment and matching their locations within ±10 mm. Simulation data suggests that flush case instability is due to interactions between jets and wall flows at the bottom of the tank. The clear flow dichotomy exhibited by this two-jet setup presents an unambiguous case to test the ability of CFD tools to predict subtle flow field changes driven by minor modifications in geometry in the context of thermal striping.

42 ENGINEERING↗
DOE OSTI2023

Development of experimental and computational frameworks to predict subcooled flow boiling in the LANL Isotope Production Facility

Cooling is crucial to maintain the integrity of target systems in isotope production facilities. At Los Alamos National Laboratory (LANL)’s Isotope Production Facility (IPF), multiple encapsulated targets are stacked and irradiated in tandem with a 100 MeV, ~250μA proton beam. To facilitate effective heat removal, these stacked targets are separated and cooled via a series of water channels. At these beam currents, this high-energy proton beam heats the target system, likely initiating subcooled flow boiling in the cooling channels. However, in-beam monitoring of the IPF target system is not possible due to the extreme radiation environment, and the necessarily significant shielding. To better understand high-power target performance, we developed ex-situ experimental and computational frameworks to predict the behavior of subcooled flow boiling at IPF. Subcooled flow boiling experiments on Inconel 625 samples under IPF conditions (2 bar pressure, 10 GPM flow rate (i.e., 2249 kg/m 2 /s), 85 K subcooling) revealed that IPF's average operating power is at the early stage of boiling with a heat transfer coefficient of 48,000 W/m 2 /s. The proposed modeling framework enables us to predict a complete boiling curve, i.e., single-phase heat transfer, onset of nucleate boiling, two-phase heat transfer, and critical heat flux (CHF), with specification of input boiling parameters up to intermediate heat flux levels. The estimated CHF under IPF conditions is 5.2 MW/m 2 . Experimental data under reduced conditions (2 bar pressure, 1.5 GPM flow rate (i.e., 337 kg/m 2 /s), 45 K subcooling) served as validation cases for the computational modeling. This computational model can be further extended to more complicated systems replicating the real IPF configuration, for instance, to study void distribution as a function of the incident proton beam profile and coolant velocity profile of multiple cooling channels. Finally, the proposed experimental and computational frameworks provide a means to better understand cooling systems in the isotope production facilities at different accelerators, where in-beam monitoring of the cooling process is not available.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2017

Modeling of a thermal energy storage system based on coupled metal hydrides (magnesium iron – sodium alanate) for concentrating solar power plants

Concentrating solar power plants represent low cost and efficient solutions for renewable electricity production only if adequate thermal energy storage systems are included. Metal hydride thermal energy storage systems have demonstrated the potential to achieve very high volumetric energy densities, high exergetic efficiencies, and low costs. The current work analyzes the technical feasibility and the performance of a storage system based on the high temperature Mg 2 FeH 6 hydride coupled with the low temperature Na 3 AlH 6 hydride. To accomplish this, a detailed transport model has been set up and the coupled metal hydride system has been simulated based on a laboratory scale experimental configuration. Proper kinetics expressions have been developed and included in the model to replicate the absorption and desorption process in the high temperature and low temperature hydride materials. The system showed adequate hydrogen transfer between the two metal hydrides, with almost complete charging and discharging, during both thermal energy storage and thermal energy release. The system operating temperatures varied from 450°C to 500°C, with hydrogen pressures between 30 bar and 70 bar. This makes the thermal energy storage system a suitable candidate for pairing with a solar driven steam power plant. The model results, obtained for the selected experimental configuration, showed an actual thermal energy storage system volumetric energy density of about 132 kWh/m 3 , which is more than 5 times the U.S. Department of Energy SunShot target (25 kWh/m 3 ).

14 SOLAR ENERGY↗
DOE OSTI2024

Spatial decay of electromagnetic waves from hypervelocity impact plasmas

Due to the high kinetic energy of impactors in space, hypervelocity impacts on spacecraft can ablate spacecraft material and produce plasmas. These plasmas have been observed in ground-based experiments to produce electromagnetic pulses (EMPs) and are suspected to be a cause of spacecraft anomalies. The initial EMP that is produced after plasma formation is replicated using a 2D electromagnetic plasma simulation, consisting of a discontinuous Galerkin, particle-in-cell (DG-PIC) framework. The spontaneous formation of EMPs from an initial quasineutral plasma is consistent with previous works and show behavior similar to that of waves emitted by antennas. The EMP emits from the plasma with a frequency less than the plasma frequency associated with the point of maximum density, showing the importance of the density gradient, and the far-field region of the “plasma antenna” is estimated. From the proportionality between plasma frequency and EMP frequency, there is a direct correlation between plasma density and radiation. Finally, as the plasma density decreases, the EMP frequency decreases, which suggests that several observable phenomena, such as optical flashes, radio frequency emission, and microwave emission, all result from the same plasma-related mechanism within a single hypervelocity impact.

42 ENGINEERING↗
DOE OSTI2022

A physics-informed multi-agents model to predict thermo-oxidative/hydrolytic aging of elastomers

This paper introduces a novel physics-informed multi-agents constitutive model to propose prediction in quasi-static constitutive behavior of cross-linked elastomer and the loss of mechanical performance during environmental aging. The presented model is used to simulate the effect of single-mechanism chemical aging (i.e. thermal-inducedor hydrolytic aging) on the behavior of the material in this hybrid framework. Those environmental single-mechanism damages change the polymer matrix over time due to massive chain scission, chain formations, and changing the arrangement of molecules in the polymer matrix. Here we propose a data-driven super-constrained machine-learned engine to represent damage in the polymer matrix and capture the changes in material behavior, including its inelastic features such as Mullins effect and permanent set in the course of aging. We have simplified the 3D stress–strain tensor mapping problem into a small number of super-constrained 1D mapping problems by means of a sequential order reduction. An assembly of multiple replicated conditional neural-network learning-agents (L-agents) is trained to systematically simplify the high-dimensional mapping problem into multiple 1D problems, each represented by a different type of agent. Our hybrid framework is designed to capture the effect of deformation history, aging time, and aging temperature. The model is validated with respect to a comprehensive set of experiments specifically designed to benchmark model capabilities and also against available data in the literature. Thermodynamic consistency and frame independency have been verified. Besides acceptable predictive abilities, a significant reduction of computational cost to predict behavior at multiple states of deformation is the most significant feature of this model.

42 ENGINEERING↗
DOE OSTI2023

A physics-driven Σ-Y atomization model for heavy-duty engine simulations

The atomization of a liquid jet is a multi-scale and multi-physics problem of interest for many engineering applications. Particularly, it drives the fuel-oxidizer mixing and dictates the efficiency of combustion engines. High-fidelity multi-phase simulations remain challenging due to the excessive computational cost required to capture all the atomization scales. Thus, atomization models are necessary to represent sub-grid liquid structures. Here, in this work, a modified Σ-Y model in the context of Eulerian-Lagrangian Spray Atomization (ELSA) is used to transport the surface area density of the spray. The model's predictive performance is assessed under various operating conditions relevant to heavy-duty engines using the Engine Combustion Network (ECN) Spray C and Spray D research-grade injectors. The classic droplet collision formulation of the Σ-Y model alone does not replicate the response of the measured spray surface area to changes in injector, ambient pressure, and injection pressure, requiring individual tuning of the model's parameters. Instead, a transition between dense and dilute spray breakup mechanisms is proposed in terms of the average droplet spacing. The collision breakup mechanism represents the dilute spray, whereas the droplet size in the dense spray is driven by a competition between the integral scale of turbulence and the balance between the turbulent kinetic energy and the surface energy of the droplets. Such an approach minimizes the requirements for model tuning. Moreover, the role of the constants of a compressible standard k-ε RANS model is assessed in the injector's internal flow and external spray simulation framework, and an updated set is proposed. The results are validated against x-ray radiography and Ultra-Small Angle X-ray Scattering (USAXS) data and highlight the predictive capabilities of the proposed physics-driven Σ-Y model, which is compatible with engine simulation turnaround times.

42 ENGINEERING↗
DOE OSTI2021

A probabilistic creep model incorporating test condition, initial damage, and material property uncertainty

Uncertainty is prevalent in the creep resistance of alloys, where at elevated temperature and low pressure, rupture can range across logarithmic decades. In this study, a probabilistic continuum-damage-mechanics (CDM)-based model is derived to capture the uncertainty of creep resistance. To meet this objective, creep data for alloy 304 Stainless Steel is gathered. A constitutive model, “Sinh”, is calibrated deterministically to determine the statistical variability of the material properties. Three sources of uncertainty are injected into the model: test condition (stress and temperature), initial damage, and material properties. Probabilistic simulations are carried out by (a) calibrating probability distribution functions (pdfs) for each source of uncertainty (b) randomly sampling the pdfs using Monte Carlo methods and (c) executing simulations to replicate the uncertain creep behavior. A sensitivity analysis is performed to evaluate the relative effect of each source of uncertainty. In full probabilistic simulations, the cumulative uncertainty of creep behavior is evaluated. The probabilistic model accurately predicts the creep deformation and rupture of the available experiments. The probabilistic model is validated for interpolation but lacks extrapolation ability. Several future works are proposed to further improve the model.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Validation of strongly coupled geomechanics and gas hydrate reservoir simulation with multiscale laboratory tests

In this work, we validate a coupled flow-geomechanics simulator for gas hydrate deposits, named T+M AM , performing two meter-scale laboratory experiments of gas hydrates for production by depressurization, replicating the gas hydrate deposit in the Ulleung Basin, East Sea, South Korea. The first experiment with a sand-only specimen is a 1D 1 m-scale depressurization test based on the excess gas method, which represents the grain coating hydrate growth. On the other hand, the second is a 3D 1.5 m-scale test with the excess water method for a sand-mud alternating layer system, representing the pore filling hydrate growth. We measure production and displacement at the top with different depressurization levels. In particular, the 3D test exhibits high coupling strength of substantial deformation induced by incompressibility of water and high deformability of the specimen. For validation, we match pressure, flow rate, and displacement between the experimental data and numerical results. Thus, we identify that T+M AM is a reliable simulator, which can be applied to fields in both permafrost and deep oceanic hydrate deposits of strongly coupled flow and geomechanics systems. This validation also implies that other coupled simulators based on the same coupling formulation as T+M AM can be validated when individual flow and geomechanics simulators are stable and reliable.

02 PETROLEUM↗
DOE OSTI2021

Small Animal IMRT Using 3D-Printed Compensators

Preclinical radiation replicating clinical intensity modulated radiation therapy (IMRT) techniques can provide data translatable to clinical practice. For this work, treatment plans were created for oxygen-guided dose-painting in small animals using inverse-planned IMRT. Spatially varying beam intensities were achieved using 3-dimensional (3D)-printed compensators.

3D PRINTING↗