Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “permeability models”

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 415 records · Page 23

Optimum dimensions of power solenoids for magnetic suspension

Design optimization of power solenoids for controllable and stabilizable magnetic suspensions with force compensation in a wind tunnel is shown. It is assumed that the model of a levitating body is a sphere of ferromagnetic material with constant magnetic permeability. This sphere, with a radius much smaller than its distance from the solenoid above, is to be maintained in position on the solenoid axis by balance of the vertical electromagnetic force and the force of gravitation. The necessary vertical (axial) force generated by the solenoid is expressed as a function of relevant system dimensions, solenoid design parameters, and physical properties of the body. Three families of curves are obtained which depict the solenoid power for a given force as a function of the solenoid length with either outside radius or inside radius as a variable parameter and as a function of the outside radius with inside radius as a variable parameter. The curves indicate the optimum solenoid length and outside radius, for minimum power, corresponding to a given outside radius and inside radius, respectively.

Kaznacheyev, B. A.↗

Processing and Characterization of PETI Composites Fabricated by High Temperature VARTM

The use of composites as primary structures on aerospace vehicles has increased dramatically over the past decade. As these advanced structures increase in size and complexity, their production costs have grown significantly. A major contributor to these manufacturing costs is the requirement of elevated processing pressures, during the thermal cure, to create fully consolidated composites. For certain composite parts, high temperature vacuum assisted resin transfer molding (HT-VARTM) can offer reduced fabrication costs compared to conventional autoclave techniques. The process has been successfully used with phenylethynyl terminated imide (PETI) resins developed by NASA LaRC. In the current study, two PETI resins, LARC(TradeMark) PETI-330 and LARC(TradeMark) PETI-8 have been used to make test specimens using HT-VARTM. Based on previous work at NASA LaRC, larger panels with a quasi-isotropic lay-up were fabricated. The resultant composite specimens exhibited void contents <3% by volume depending on the type of carbon fabric preform used. Mechanical properties of the panels were determined at both room and elevated temperatures. These included open-hole compressive (OHC) and short beam shear (SBS) properties. Limited process modeling efforts were carried out including infusion times, composite panel size limitations and fabric permeability characterization. Work has also been carried out to develop new PETI based resins specifically geared towards HT-VARTM. The results of this work are presented herein.

Ghose, Sayata↗

Global Geo-processed Data of Aquifer Properties by 0.5° Grid, Country and Water Basins

This repository of global hydrogeologic datasets contains aquifer properties on 0.5° scale, including depth to groundwater (Fan et al., 2013), aquifer thickness (de Graaf et al., 2015), WHYMap aquifer classes (Richts et al., 2011), recharge (Döll and Fiedler, 2008; Gleeson et al., 2016), lakes (Messager et al., 2016), porosity and permeability (Gleeson et al., 2014), digitized and geo-processed from their respective sources. Globally gridded aquifer properties could be used independently to estimate global groundwater availability or used as critical inputs to the superwell model to simulate groundwater extraction and provide estimates of pumped volumes and unit costs under user-specific scenarios. Key resources related to this data are: Niazi, H., Ferencz, S. B., Graham, N. T., Yoon, J., Wild, T. B., Hejazi, M., Watson, D. J., & Vernon, C. R. (2025). Long-term hydro-economic analysis tool for evaluating global groundwater cost and supply: Superwell v1.1. Geoscientific Model Development, 18(5), 1737-1767. https://doi.org/10.5194/gmd-18-1737-2025 superwell model repository which uses this data to simulate groundwater extraction and provides estimates of the global extractable volumes and unit-costs ($/km3) of accessible groundwater production under user-specified extraction scenarios. Repository Overview Main output: aquifer_properties_rec.csv contains all processed outputs, including aquifer properties like porosity, permeability, recharge, lake areas, aquifer thickness, and depth to groundwater. shapefiles.zip: contains all digitized GIS databases and shapefile for all aquifer properties prep_inputs.R and prep_inputs_recharge_lakes.R: R scripts that process the shapefiles to produce the aquifer_properties_rec.csv file plot_inputs.R: R script for plotting the maps and conducting preliminary analysis on the available groundwater volume basin_to_country_mapping.csv, basin_country_region_mapping.csv and continent_county_mapping.csv provide the mapping between continents, 32 energy-economic macro regions, countries, and water basins for post-processing aquifer_properties_rec.csv Maps: Each map visualizes the spatial distribution of one of the aquifer properties across the globe map_in_Porosity.png map_in_Permeability.png map_in_Aquifer_thickness.png map_in_Depth_to_water.png map_in_Recharge.png map_in_Grid_area_km.png map_in_Lake_area_km.png map_in_WHYClass.png Sample inputs sample_inputs.py: this script samples inputs from the aquifer_properties_rec dataset, ensuring the sampled and original inputs maintain the same distributions sampled_data_100.csv contains 100 sampled data points and sampled_data_100.png compares their distributions Dataset Overview The main outputs are consolidated in a comprehensive aquifer_properties_rec.csv file and include the following fields: GridCellID: Unique identifier for each (roughly 0.5°) grid cell Continent: Continent name Country: Country name GCAM_basin_ID: Identifier for GCAM hydrologic basin Basin_long_name: Full name of the basin WHYClass: Hydrogeologic classification based on WHYMap aquifer classes (Richts et al., 2011) Porosity: Soil porosity (%) (Gleeson et al., 2014) Permeability: Soil permeability (in square meters; Gleeson et al., 2014) Aquifer_thickness: Thickness of the aquifer (in meters; de Graaf et al., 2015) Depth_to_water: Depth to groundwater (in meters; Fan et al., 2013) Recharge: long-term annual averaged recharge rates (in m/yr; Döll and Fiedler, 2008; Gleeson et al., 2016) Grid_area: Area of the grid cell (in square meters) Lakes_area: Area of inland lakes (in square meters; Messager et al., 2016) Key References The datasets are digitized versions of global hydrogeologic properties from the following key literature sources: Depth to Groundwater: Fan, Y., Li, H., & Miguez-Macho, G. (2013). Global Patterns of Groundwater Table Depth. Science, 339(6122), 940-943. https://doi.org/10.1126/science.1229881 Aquifer Thickness: de Graaf, I. E. M., Sutanudjaja, E. H., van Beek, L. P. H., & Bierkens, M. F. P. (2015). A high-resolution global-scale groundwater model. Hydrol. Earth Syst. Sci., 19(2), 823-837. https://doi.org/10.5194/hess-19-823-2015 Porosity and Permeability: Gleeson, T., Moosdorf, N., Hartmann, J., & van Beek, L. P. H. (2014). A glimpse beneath earth's surface: GLobal HYdrogeology MaPS (GLHYMPS) of permeability and porosity. Geophysical Research Letters, 41(11), 3891-3898. https://doi.org/10.1002/2014GL059856 Aquifer classes: Richts, A., Struckmeier, W. F., & Zaepke, M. (2011). WHYMAP and the Groundwater Resources Map of the World 1:25,000,000. In J. A. A. Jones (Ed.), Sustaining Groundwater Resources: A Critical Element in the Global Water Crisis (pp. 159-173). Springer Netherlands. https://doi.org/10.1007/978-90-481-3426-7_10 Recharge: Döll, P., & Fiedler, K. (2008). Global-scale modeling of groundwater recharge. Hydrol. Earth Syst. Sci., 12(3), 863-885. https://doi.org/10.5194/hess-12-863-2008; Gleeson, T., Befus, K. M., Jasechko, S., Luijendijk, E., & Cardenas, M. B. (2016). The global volume and distribution of modern groundwater. Nature Geoscience, 9(2), 161-167. https://doi.org/10.1038/ngeo2590 Inland Lakes: Messager, M. L., Lehner, B., Grill, G., Nedeva, I., & Schmitt, O. (2016). Estimating the volume and age of water stored in global lakes using a geo-statistical approach. Nature Communications, 7(1), 13603. https://doi.org/10.1038/ncomms13603 Cite as Niazi, H., Watson, D., Hejazi, M., Yonkofski, C., Ferencz, S., Vernon, C., Graham, N., Wild, T., & Yoon, J. (2024). Global Geo-processed Data of Aquifer Properties by 0.5° Grid, Country and Water Basins. MultiSector Dynamics-Living, Intuitive, Value-adding, Environment. https://doi.org/10.57931/2484226 Contact Reach out to Hassan Niazi or open an issue in the superwell repository for questions or suggestions.

aquifer thickness↗

DECOVALEX-2023: Task C Final Report

The Full-scale Emplacement (FE) heater experiment at the Mont Terri Underground Rock Laboratory (URL) was designed and conducted by Nagra to replicate an emplacement tunnel of Nagra’s reference repository design at 1:1 scale. Alongside testing the technical feasibility of constructing disposal tunnels, emplacing waste containers in the tunnels and then backfilling them, the main goals of the FE experiment are (1) to obtain a better understanding of the coupled effects of induced thermo-hydro-mechanical (THM) processes that may occur and (2) to validate existing coupled THM models (Müller et al., 2017). A key aspect of ensuring safety for repositories located in low-permeability rock involves minimizing any damage to the rock itself, thereby preserving its integrity and promoting a stable environment Amongst a number of processes that could damage the rock is the increase in pore pressure due to thermal loading caused by heat emitted from the waste. To reduce the potential damage of the rock, it is important to analyse the evolution of heat over time due to the heat load of the containers and assess possible consequences by coupled THM models. The aim of Task C of DECOVALEX-2023 was to build 3D numerical models of the FE experiment, focussing in particular on the heating induced pore pressure change in the Opalinus Clay. Data from a large number of sensors were available from the FE experiment for model comparison. These sensors measured temperature and relative humidity in the bentonite around the heaters, and temperature, pressure and displacement/strain in the surrounding Opalinus clay. Data were available from the start of excavation (April 2012) up to August 2020 for most sensors (more than 5 years from the start of heating in December 2014). To fulfil the overall aim of the task, the work was broken down into a number of steps, starting with simpler models to build confidence in each team’s approach and then moving to more complex models that better represent the FE experiment. Step 0 consisted of 2D benchmark models, gradually increasing the number of processes that are represented from thermal (T) only models in Step 0a, to coupled thermal hydraulic (TH) models in Step 0b with a representation of changing porosity, to coupled thermo-hydro-mechanical (THM) models in Step 0c, where porosity changes are calculated by the mechanical model. A detailed specification of processes, parameters, initial and boundary conditions was provided for this step, with the ambition that all teams would work towards close agreement in their model results, thus building confidence in the model implementations. vi It was not straightforward to achieve agreement between the teams, so additional steps (Step 0b2, 0b3, 0c2, 0c3) were added along with derivation of some analytical solutions against which the models could be compared. The reasons for the differences between teams were investigated and found to be caused primarily by different conceptual model assumptions (including temperature dependence of the thermal expansion of water), different model formulations (including porosity evolution) and differences in modelled domain sizes, boundary conditions and grid discretisation. This demonstrates that comparisons between multiple modelling teams and/or comparison with analytical results and experimental data are highly beneficial in providing an indication of uncertainty in model predictions. At the conclusion of Step 0, almost all teams had achieved a close agreement in model results and those that had not achieved an agreement knew the reason for this. Step 1 moved from 2D models to 3D models of the FE experiment without adding technical features like shotcrete or EDZ, and only considering the heating phase. Initially the 3D model was tightly specified to continue to build confidence in the model implementations (Step 1a). The results of Step 1a were compared to the data from the FE-experiment without the teams seeing the data. The teams were then provided with a sub-set of the data from the FE-experiment and invited to consider how best to use the large dataset for model comparison (Step 1b). Teams were then asked to use the data provided to calibrate their models, only changing material property values rather than adding features or processes to their models (Step 1c). In Step 1, teams were asked to only model the heating phase of the experiment, so pressure in the Opalinus Clay was reported as change in pressure since the initial conditions were specified rather than modelled. The change from 2D to 3D models was accompanied by an increase in the dispersion of results between the teams. Some of this was resolved during the task, but some remained and is potentially due to model discretisation. Calibration of parameters was useful in improving the fit of the models to the data but the remaining differences indicated that the models were missing features or processes. In Step 2, the teams were asked to update their models with additional features and processes as well as calibrating parameters to try and improve the fit of the models to the data. Teams were encouraged to represent ventilation of the open FE tunnel prior to backfilling with heaters and bentonite and in Step 2, the absolute pressure in the Opalinus Clay was compared between the teams. Teams took different approaches, but there was consideration of adding shotcrete and an EDZ into the model, representing stress change during excavation and different approaches to modelling ventilation of the FE tunnel. Overall, the documented results showed a very good agreement for temperature. The results for porewater pressure evolution showed a significant improvement for most teams compared to Step 1c with a good agreement to the measurements for several teams whereas some teams overpredicted the pressure increase and others overpredicted the drainage effect especially for the sensors close to the heater. Step 3 was an opportunity for teams to use the models developed in Step 1 and Step 2 to make predictions about the temperature and pressure changes that will be expected at the FE experiment over the next few years in light of the planned changes in thermal output of the heaters.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Prediction of Long-Term Geochemical Change in Bentonite Based on the Interpretative THMC Model of the FEBEX In Situ Test

Since nuclear energy is crucial in the decarbonization of the energy supply, one hurdle to remove is the handling of high-level radioactive waste (HLW). Disposal of HLW in a deep geological repository has long been deemed a viable permanent option. In the design of a deep geological repository, compacted bentonite is the most commonly proposed buffer material. Predicting the long-term chemical evolution in bentonite, which is important for the safety assessment of a repository, has been challenging because of the complex coupled processes. Models for large-scale tests and predictions based on such models have been some of the best practices for such purposes. An 18-year-long in situ test with two dismantling events provided a unique set of chemical data that allowed for studying chemical changes in bentonite. In this paper, we first developed coupled thermal, hydrological, mechanical, and chemical (THMC) models to interpret the geochemical data collected in the in situ test and then extended the THMC model to 200 years to make long-term prediction of the geochemical evolution of bentonite. The interpretive coupled THMC model shows that the geochemical profiles were strongly affected by THM processes such as evaporation/condensation, porosity change caused by swelling, permeability change, and the shape of concentration profiles for major cations were largely controlled by transport processes, but concentration levels were regulated by chemical reactions, and the profiles of some species such as pH, bicarbonate, and sulfate were dominated by these reactions. The long-term THMC model showed that heating prolongs the time that bentonite becomes fully saturated in the area close to the heater/canister; however, once the bentonite becomes fully saturated, high concentrations of ions in bentonite near the heater, which was observed in the field test, will disappear; illitization continues for 50 years but will not proceed further.

Zheng, Liange↗

Using onset times from frequent geophysical surveys to understand reservoir fluid flow and to estimate flow properties

Seismic characterization of flow properties is a difficult prospect, due to the indirect relationship between permeability and seismic velocity and attenuation. At best, seismic time-lapse changes can detect the effects of saturation and pressure changes in a reservoir due to fluid flow. Even in this case, the interpretation of the seismic observation in terms of the state of the reservoir depends intimately on the properties of the rock physics model, which are usually poorly known. The onset time, the calendar time of geophysical changes, provides an alternative datum for characterizing properties such as reservoir permeability. The main advantage of an onset time is that it is sensitive to the flow properties of the reservoir yet insensitive to the details of the rock physics model. Two examples of the utility of onset times are discussed: The use of travel time shifts induced by the injection of carbon dioxide between two wells, and the time-lapse time shifts for elastic waves propagating through a reservoir undergoing enhanced oil recovery. In both examples, the onset times are mapped into permeability estimates using a trajectory-based approach akin to seismic tomography.

58 GEOSCIENCES↗

Geomechanical and flow implications with continued bioconversion of coal to methane: Experiments and modeling

Microbial conversion of coal to methane is a promising technology during transition of coal-based energy from conventional coal mining to natural gas recovery from coal. Significant research advances have been made towards engineering ideal microbial communities and nutrients for bio-stimulation of coal. However, actual field applications require geomechanical and flow behavior characterization of coal during the bioconversion process as well as gas production over the life of the created “biogenic gas reservoirs”. This work presents the results of an experimental investigation to estimate gas/methane production using bioconversion of coal and analyses of variations in bulk modulus, strain and permeability. Bulk modulus of coal demonstrated time-dependent behavior with continued bioconversion of solid coal. The phenomenon was modeled using mass balancing in a closed environment and the logistic equation based biogenic gas production with time. The modeled gas production and changes in modulus of coal showed excellent agreement with the corresponding experimental results. Next, numerical simulation of biogenic conversion of coal in a constant stress in-situ condition, replicating the Huff ‘n Puff method under field conditions, was carried out. The results showed that bioconversion can lead to decrease in effective stress, increase in permeability and decrease in the modulus of coal with time. Repeated recharge of coal with nutrients and recovering the produced gas, that is, sequential cycles of Huff n’ Puff, showed further increase in permeability and decrease in coal strength, potentially leading to coal failure in-situ, further increasing the permeability, thus enhancing the prospect of field application of the technology.

01 COAL, LIGNITE, AND PEAT↗

Estimates for the Aerodynamic Coefficients of Ringsail and Disk-Gap-Band Parachutes Operating on Mars

The creation of flight dynamics simulations is often part of the design, development, and flight operations of missions that include entry, descent, and landing (EDL) in their concept of operation. Such flight dynamics simulations include models for many of the EDL sub-systems. For EDL systems that include a parachute, models of the parachute’s aerodynamic characteristics must be created for inclusion in the flight dynamics simulations. Models are presented for the aerodynamic coefficients (i.e., C(D),C(Tot),C(T),C(N), and C(m,SLCP)) of Supersonic Ringsail (SSRS) and Disk-Gap-Band (DGB) parachutes as functions of total porosity, λ(T), Mach number, M, and total angle of attack, λ(T) (when necessary). The source aerodynamic coefficients data used for creating these models were obtained from a wind tunnel test of subscale parachutes. In this wind tunnel test, subscale parachutes of both parachute types were fabricated from two different fabrics with very different permeabilities. By varying the fabric permeability, while maintaining the parachute geometry constant, it was possible to vary λ(T). The fabric permeability test data necessary for the calculation of λ(T) were obtained from samples of the same fabrics used to fabricate the subscale parachutes. Although the models for the aerodynamic coefficients are simple polynomial functions of λ(T) and M, they are capable of producing good reproductions of the source data. The (λ(T), M) domains over which these models are applicable are clearly defined. The models and their domains are applicable to flight operations on Mars. The models presented in this report are not unique – other models could be created from the source data. To make the creation of alternative numerical models, the aerodynamic coefficients source data are included in this report.

Cruz, Juan R.↗

Computational Modeling and Evolutionary Implications of Biochemical Reactions in Bacterial Microcompartments

Bacterial microcompartments (BMCs) are protein-encapsulated compartments found across at least 23 bacterial phyla. BMCs contain a variety of metabolic processes that share the commonality of toxic or volatile intermediates, oxygen-sensitive enzymes and cofactors, or increased substrate concentration for magnified reaction rates. These compartmentalized reactions have been computationally modeled to explore the encapsulated dynamics, ask evolutionary-based questions, and develop a more systematic understanding required for the engineering of novel BMCs. Many crucial aspects of these systems remain unknown or unmeasured, such as substrate permeabilities across the protein shell, feasibility of pH gradients, and transport rates of associated substrates into the cell. This review explores existing BMC models, dominated in the literature by cyanobacterial carboxysomes, and highlights potentially important areas for exploration.

bacterial microcompartments↗

Thermo-Poro-Mechanical Modeling of RTV Intumescence

Room temperature vulcanizing (RTV) silicone is a high-temperature adhesive used as a gap-filler between heatshield tiles in numerous entry missions. Its propensity to intumesce, or swell upon exposure to heat, is a well-known effect that needs to be carefully quantified during design. At tile interfaces of charring ablators, intumescence, combined with differential recession, could cause the gap filler to protrude past the ablator outer mold line, forming a “fence”. Fencing can in turn cause transition to turbulence of the flow wetting the heat shield, leading to augmented surface heating. Recent experiments conducted at the Plasmatron X facility, the high enthalpy wind tunnel of the Center for Hypersonics and Entry Systems Studies, have shown prominent fencing of RTV gap fillers in PICA, under both nitrogen and air plasmas. Similar observations are well known in the arcjet literature. Further experiments under controlled environment, performed using in situ X-ray micro-computed tomography (micro-CT) at the Advanced Light Source (ALS), have shown heating rate-dependent swelling and shrinkage of RTV during pyrolysis. To simulate RTV intumescence, a novel model was introduced in the Porous Materials Analysis Toolbox based on OpenFOAM, PATO, to account for pore-pressure buildup within both closed- and open-pores. The governing equation for the thermo-poro-mechanical response were developed, assuming linear elasticity for the charring silicone. A new multi-pyrolysis model that tracks non-monotonic advancement of material properties with pyrolysis was proposed. This model addresses the limitations of state-of-the-art ablator models to capture the different stages of thermal degradation and coupled thermomechanics. Swelling of RTV was simulated using the new thermo-poro-mechanical model and compared against in situ micro-CT data. Results showed good agreement in intumescence height and temperature profiles at all heating rates, indicating that the key factor contributing to RTV swelling is the internal pressure build-up within closed- and open-pores. As RTV is cured into a soft (rubbery) compound with low-porosity and permeability, initial temperature increase and pyrolysis gas production cause a significant increase of internal pressure, causing a pronounced volume growth. As thermal degradation progresses, rigidization of the silicone occurs due to char hardening which counteract volume shrinkage after gas pressure relief. Overall, our model shows that accounting for changes in properties such porosity, permeability and key thermomechanical coefficients is crucial for capturing the RTV volume change during ablation and enable a predictive capability for heatshield tile interface response. A plan for future calibration of thermomechanical properties evolution during degradation is discussed, as a key next step to close the new model.

RTV↗

Thermo-Poro-Mechanical Modeling of RTV Intumescence

Room temperature vulcanizing (RTV) silicone is a high-temperature adhesive used as a gap-filler between heatshield tiles in numerous entry missions. Its propensity to intumesce, or swell upon exposure to heat, is a well-known effect that needs to be carefully quantified during design. At tile interfaces of charring ablators, intumescence, combined with differential recession, could cause the gap filler to protrude past the ablator outer mold line, forming a “fence”. Fencing can in turn cause transition to turbulence of the flow wetting the heat shield, leading to augmented surface heating. Recent experiments conducted at the Plasmatron X facility, the high enthalpy wind tunnel of the Center for Hypersonics and Entry Systems Studies, have shown prominent fencing of RTV gap fillers in PICA, under both nitrogen and air plasmas. Similar observations are well known in the arcjet literature. Further experiments under controlled environment, performed using in situ X-ray micro-computed tomography (micro-CT) at the Advanced Light Source (ALS), have shown heating rate-dependent swelling and shrinkage of RTV during pyrolysis. To simulate RTV intumescence, a novel model was introduced in the Porous Materials Analysis Toolbox based on OpenFOAM, PATO, to account for pore-pressure buildup within both closed- and open-pores. The governing equation for the thermo-poro-mechanical response were developed, assuming linear elasticity for the charring silicone. A new multi-pyrolysis model that tracks non-monotonic advancement of material properties with pyrolysis was proposed. This model addresses the limitations of state-of-the-art ablator models to capture the different stages of thermal degradation and coupled thermomechanics. Swelling of RTV was simulated using the new thermo-poro-mechanical model and compared against in situ micro-CT data. Results showed good agreement in intumescence height and temperature profiles at all heating rates, indicating that the key factor contributing to RTV swelling is the internal pressure build-up within closed- and open-pores. As RTV is cured into a soft (rubbery) compound with low-porosity and permeability, initial temperature increase and pyrolysis gas production cause a significant increase of internal pressure, causing a pronounced volume growth. As thermal degradation progresses, rigidization of the silicone occurs due to char hardening which counteract volume shrinkage after gas pressure relief. Overall, our model shows that accounting for changes in properties such porosity, permeability and key thermomechanical coefficients is crucial for capturing the RTV volume change during ablation and enable a predictive capability for heatshield tile interface response. A plan for future calibration of thermomechanical properties evolution during degradation is discussed, as a key next step to close the new model.

silicone intumescence↗

Global Geo-processed Data of Aquifer Properties by 0.5° Grid, Country and Water Basins

This repository of global hydrogeologic datasets contains aquifer properties on 0.5° scale, including depth to groundwater (Fan et al., 2013), aquifer thickness (de Graaf et al., 2015), WHYMap aquifer classes (Richts et al., 2011), porosity and permeability (Gleeson et al., 2014), digitized and geo-processed from their respective sources. Globally gridded aquifer properties could be used independently to estimate global groundwater availability or used as critical inputs to the superwell model to simulate groundwater extraction and provide estimates of pumped volumes and unit costs under user-specific scenarios. Key resources related to this data are: Niazi, H., Ferencz, S., Graham, N., Yoon, J., Wild, T., Hejazi, M., Watson, D., & Vernon, C. (2024; In-prep). Long-term Hydro-economic Assessment Tool for Evaluating Global Groundwater Cost and Supply: Superwell v1. Geoscientific Model Development. superwell model repository which uses this data to simulate groundwater extraction and provides estimates of the global extractable volumes and unit-costs ($/km3) of accessible groundwater production under user-specified extraction scenarios. Repository Overview Main output: aquifer_properties.csv contains all processed outputs, including aquifer properties like porosity, permeability, aquifer thickness, and depth to groundwater. shapefiles.zip: contains all digitized GIS databases and shapefile for all aquifer properties prep_inputs.R: R script that processes the shapefiles to produce the aquifer_properties file plot_inputs.R: R script for plotting the maps and conducting preliminary analysis on the available groundwater volume basin_to_country_mapping.csv, basin_country_region_mapping.csv and continent_county_mapping.csv provide the mapping between continents, 32 energy-economic macro regions, countries, and water basins for post-processing aquifer_properties.csv Maps: Each map visualizes the spatial distribution of one of the aquifer properties across the globe map_in_Porosity.png map_in_Permeability.png map_in_Aquifer_thickness.png map_in_Depth_to_water.png map_in_Grid_area_km.png map_in_WHYClass.png Sample inputs sample_inputs.py: this script samples inputs from the aquifer_properties dataset, ensuring the sampled and original inputs maintain the same distributions sampled_data_100.csv contains 100 sampled data points and sampled_data_100.png compares their distributions Dataset Overview The main outputs are consolidated in a comprehensive aquifer_properties.csv file and include the following fields: GridCellID: Unique identifier for each (roughly 0.5°) grid cell Continent: Continent name Country: Country name GCAM_basin_ID: Identifier for GCAM hydrologic basin Basin_long_name: Full name of the basin WHYClass: Hydrogeologic classification based on WHYMap aquifer classes (Richts et al., 2011) Porosity: Soil porosity (%) (Gleeson et al., 2014) Permeability: Soil permeability (in square meters; Gleeson et al., 2014) Aquifer_thickness: Thickness of the aquifer (in meters; de Graaf et al., 2015) Depth_to_water: Depth to groundwater (in meters; Fan et al., 2013) Grid_area: Area of the grid cell (in square meters) Key References The datasets are digitized versions of global hydrogeologic properties from the following key literature sources: Depth to Groundwater: Fan, Y., Li, H., & Miguez-Macho, G. (2013). Global Patterns of Groundwater Table Depth. Science, 339(6122), 940-943. https://doi.org/10.1126/science.1229881 Aquifer Thickness: de Graaf, I. E. M., Sutanudjaja, E. H., van Beek, L. P. H., & Bierkens, M. F. P. (2015). A high-resolution global-scale groundwater model. Hydrol. Earth Syst. Sci., 19(2), 823-837. https://doi.org/10.5194/hess-19-823-2015 Porosity and Permeability: Gleeson, T., Moosdorf, N., Hartmann, J., & van Beek, L. P. H. (2014). A glimpse beneath earth's surface: GLobal HYdrogeology MaPS (GLHYMPS) of permeability and porosity. Geophysical Research Letters, 41(11), 3891-3898. https://doi.org/10.1002/2014GL059856 Aquifer classes: Richts, A., Struckmeier, W. F., & Zaepke, M. (2011). WHYMAP and the Groundwater Resources Map of the World 1:25,000,000. In J. A. A. Jones (Ed.), Sustaining Groundwater Resources: A Critical Element in the Global Water Crisis (pp. 159-173). Springer Netherlands. https://doi.org/10.1007/978-90-481-3426-7_10 Cite as Niazi, H., Watson, D., Hejazi, M., Yonkofski, C., Ferencz, S., Vernon, C., Graham, N., Wild, T., & Yoon, J. (2024). Global Geo-processed Data of Aquifer Properties by 0.5° Grid, Country and Water Basins. MSD-LIVE Data Repository. https://doi.org/10.57931/2307831 Contact Reach out to Hassan Niazi or open an issue in superwell repository for questions or suggestions.

aquifer thickness↗

Advances in Computational Approaches for Estimating Passive Permeability in Drug Discovery

Passive permeation of cellular membranes is a key feature of many therapeutics. The relevance of passive permeability spans all biological systems as they all employ biomembranes for compartmentalization. A variety of computational techniques are currently utilized and under active development to facilitate the characterization of passive permeability. These methods include lipophilicity relations, molecular dynamics simulations, and machine learning, which vary in accuracy, complexity, and computational cost. This review briefly introduces the underlying theories, such as the prominent inhomogeneous solubility diffusion model, and covers a number of recent applications. Various machine-learning applications, which have demonstrated good potential for high-volume, data-driven permeability predictions, are also discussed. Due to the confluence of novel computational methods and next-generation exascale computers, we anticipate an exciting future for computationally driven permeability predictions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Estimates for the Aerodynamic Coefficients of Ringsail and Disk-Gap-Band Parachutes Operating on Mars

Models are presented for the aerodynamic coefficients of Supersonic Ringsail and Disk-Gap-Band parachutes as functions of total porosity, Lambda(sub t), Mach number, M, and total angle of attack, Alpha(sub t) (when necessary). The source aerodynamic coefficients data used for creating these models were obtained from a wind tunnel test of subscale parachutes. In this wind tunnel test, subscale parachutes of both parachute types were fabricated from two different fabrics with very different permeabilities. By varying the fabric permeability, while maintaining the parachute geometry constant, it was possible to vary Alpha(sub t). The fabric permeability test data necessary for the calculation of Alpha(sub t) were obtained from samples of the same fabrics used to fabricate the subscale parachutes. Although the models for the aerodynamic coefficients are simple polynomial functions of Alpha(sub t) and M, they are capable of producing good reproductions of the source data. The (Alpha(sub t), M) domains over which these models are applicable are clearly defined. The models are applicable to flight operations on Mars.

Cruz, Juan R.↗

DECOVALEX-2023: Task E Final Report

This is the Task E final report for DECOVALEX-2023. Task E is focused on understanding thermal, two-phase hydrological, and mechanical (TH2M) processes, especially related to predicting brine migration in the excavation damaged zone around a heated excavation in salt. Salt is attractive as a disposal medium for radioactive waste because it is self-healing and is essentially impermeable and essentially non-porous in the far field (away from excavations). Investigation of the short-term (days to years) near-field (centimeters to tens of meters) behavior of salt is important for radioactive waste disposal because this early period strongly controls the amount of brine in a salt repository. Brine leads to corrosion of waste forms and waste packages, and possible dissolution of radionuclides with brine transport being a potential transport vector to the accessible environment. The main test case used in Task E is the ongoing Brine Availability Test in Salt (BATS) heater test located underground at the Waste Isolation Pilot Plant (WIPP) near Carlsbad, New Mexico, USA. The Task was divided into a series of Steps. Step 0 was an introduction to processes in salt, that included matching historical unheated brine inflow data from boreholes at WIPP and matching temperature observations during BATS heater test 1a. Step 1 included validation of models against a thermo-poroelastic analytical solution relevant to heated boreholes in salt, and two-phase flow around an excavation in salt. Step 2 required all the individual components covered in steps 0 and 1 to come together to match observed brine inflow behavior during the BATS 1a heater test. There were a range of approaches from the teams, from mechanistic to prescriptive. Given the uncertainties in the problem, some teams used one- or two-dimensional models of the processes, while other teams included more geometrical complexity in three-dimensional models. The key learning points from Task E have been: • Heat conduction through salt typically requires non-linear thermal conductivity (as a function of temperature), but most models do a good job matching observations, given appropriate adjustments to the applied power and some thermocouple locations. • Thermal pressurization requires coupled thermal-hydrological-mechanical (THM) responses that consider the thermal expansion of the fluid and solid phases. • Initialization of two-phase flow models around a borehole or excavation in salt are more realistically represented as “wetting up”, rather than “drying down” (i.e., the initial state after excavation is mostly dry, rather than mostly wet). • The BATS 1a heater test includes a significant release of brine after the end of heating, which requires a large increase in permeability to recreate. Task E has been a great learning experience for all the teams involved, and feedback from the modeling teams has led to changes in the design of follow-on BATS experiments, which are now ongoing underground at WIPP. There was a balance throughout the task between freedom to model phenomena how each team saw fit, and prescriptiveness in problem design to bring the modeling teams closer together to allow attribution of smaller differences between models to different modeling choices. The modeling approaches seem to go through two phases: an early phase of discovery or testing, and a later phase of refinement and improvement. In future modeling efforts, different field data could be used (e.g., BATS 2) and more time should be included in the processes for teams to make multiple model refinement or even significant changes to their conceptual model or setup, based on lessons learned from the modeling exercise.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Estimation of stress and stress-induced permeability change in a geological nuclear waste repository in a thermo-hydrologically coupled simulation

Geologic disposal is a promising solution for a safe permanent isolation of accumulated high-level nuclear waste from nuclear powerplants. The behavior of host rock is highly coupled thermally, hydromechanically, and chemically. Numerical simulations of such coupled phenomena for the extremely long term (>100,000 years) and large length scale (>1 km) of geologic disposal remain to be computationally challenging. In this study, a methodology has been developed to approximate stress and stress-induced permeability change in host rock using only thermo-hydrological (TH) variables. A coupled thermo-hydromechanical (THM) simulation is carried out using TOUGH-FLAC simulator to model THM behaviors of a generic nuclear waste repository, in order to evaluate the performance of the developed methodology, which is implemented in a coupled TH simulation. Results show that stress and permeability change estimated by the developed methodology in the TH-coupled simulation match those calculated in the THM-coupled simulation over the simulated timeframe of over 10,000 years. Details about errors in stress and permeability estimates accrued by the developed methodology are discussed here. The developed methodology will help incorporate stress-induced permeability change into existing TH simulators for the long-term radionuclide transport in geologic disposal.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Electrical Properties of Lunar Environment Used for Predicting Lunar RF Propagation Characteristics

This contribution treats the lunar propagation environment as a three region medium: lunar exosphere, lunar regolith, and lunar bedrocks. Then it provides models for predicting the electromagnetic characteristics of each region. The electromagnetic characteristics could be electric characteristics represented by the complex relative permittivity, or magnetic characteristics represented by the complex relative permeability or both electric and magnetic characteristics. The lunar exosphere and the lunar bedrocks have only electric characteristics. The lunar regolith has both electric and magnetic characteristics. The complex relative permittivity models for lunar regolith are mapping of the corresponding models for Earth surface components reported in ITU-R P. 527-6. The complex relative permittivity prediction model of lunar exosphere is expressed in terms of a plasma frequency similar to the ordinary wave critical frequency in the corresponding model for the ionosphere. Based on this fascicle the following can be concluded for the frequency bands of 390 MHz and above: • The lunar exosphere can be treated as a free space, • The lunar regolith can be considered as non-magnetic, • The regolith complex relative permittivity has no temperature dependence, • The real part of the regolith complex relative permittivity depends only on regolith bulk density and it has no frequency dependence, and • The variation of regolith complex relative permittivity with regolith depth should be taken into consideration. Moreover, at frequencies of 2400 MHz and above, the lunar regolith can be treated as a uniform medium with complex relative permittivity equal to the corresponding complex relative permittivity at the regolith surface.

Electrical Permittivity↗

Field-scale estimation of soil properties from spectral induced polarization tomography

Estimates of soil properties such as Cation Exchange Capacity (CEC), water content, grain size characteristics, and permeability are important in geotechnical engineering, water resources, and agriculture. We develop a non-intrusive approach to estimate these properties in the field using spectral induced polarization (SIP) tomography. This geophysical method provides information about the frequency dependence of the complex electrical conductivity of porous media. Using 18 soil samples collected from a Bordeaux vineyard, we first conducted a laboratory study using SIP over the frequency range 10 mHz-45 kHz. The laboratory data were used to confirm the accuracy of a recently developed dynamic Stern layer petrophysical model. The results are consistent with published values from previous works using soils. A comparison was made by comparing the field complex conductivity spectra and the experimental data at two locations where core samples were obtained. The model was then used in concert with field data to image the spatial distribution of CEC, water content, permeability, and mean grain size along a vineyard transect. For clay and sandy textures found in the field, measured and estimated CEC agree rather well (from 6 to 40% discrepancy). Furthermore, our approach provides an efficient way to estimate important soil properties in a non-invasive manner, in high resolution, and over field-relevant scales of the critical zone of the Earth.

54 ENVIRONMENTAL SCIENCES↗