Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Porosity”

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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Comparison of aluminium foams prepared by different methods using X-ray tomography

Highlights: • The closed porosity of MF aluminium foam is much lower than GI aluminium foam at a similar overall porosity. • There are more micropores in cell walls of MF aluminium foams than that in GI aluminium foams. • GI method is more suitable to prepare aluminium foams with high porosity under a similar cell size condition. • The thickness difference between normal cell wall and Plateau border is greater in GI foams compared to MF foams. - Abstract: The cell structure and mechanical properties of aluminium foams prepared by melt foaming (MF) method are supposed to be better than the ones produced by gas injection (GI) method. Preparation processes of these two methods have been developed during the last recent years. Mechanical properties of aluminium foams depend strongly on characteristics of the cell structure. Therefore, it is necessary to compare cell structures of the foams fabricated with these two methods quantitatively using X-ray computed tomography and this is the purpose of the present paper. MF and GI foams have thus been prepared. True density measurements and cell structures indicate that the volume fraction of the closed cells of MF foam with high porosity is negligible. By contrast, the majority of cells of GI foam are closed. Results of cell structure analysis show that there are more micropores in cell walls of the foams prepared by MF method compared to the ones prepared by GI method, especially for the MF foams with low porosity. In addition, results of cell size study show that the cell size distribution is wide in dynamic GI foams, and there are usually big holes inside MF foams. GI method is more suitable to prepare aluminium foams with high porosity compared to MF method. Furthermore, the thickness difference between normal cell wall and Plateau border is greater in GI foams, in particular for the ones produced by the static injector.

36 MATERIALS SCIENCE↗
DOE OSTI2024

3-D Geological Modeling for Numerical Flow Simulation Studies of Gas Hydrate Reservoirs at the Kuparuk State 7-11-12 Pad in the Prudhoe Bay Unit on the Alaska North Slope

Accurate reservoir evaluation requires reliable three-dimensional (3-D) geological models. Here, this study conducted 3-D geological modeling for numerical flow simulation of the B1 sand gas hydrate reservoir at the Kuparuk State 7-11-12 pad, Prudhoe Bay Unit, Alaska North Slope. The model integrates well logs, core, and seismic data to address spatial heterogeneity in geological structures and reservoir properties. Two modeling types were performed: structural framework modeling and petrophysical property modeling. For structural framework modeling, seismic data and well log markers were used to reproduce subsurface structures characterized by a normal fault system. A volume-based modeling algorithm and stair-stepping grid were applied. The resulting 3-D model comprised 2,640,000 grid cells across 264 layers, including seven fault grids. For petrophysical property modeling, total porosity was initially modeled using sequential Gaussian simulation with collocated cokriging. To reproduce the upward coarsening of the B1 sand, upscaled log-derived total porosity and a three-dimensional (3-D) trend depicting total porosity variation were used as primary and secondary data, respectively. Gas hydrate saturation distribution was modeled similarly, with secondary data from estimated porosity distribution and seismic-derived acoustic impedance map enhancing accuracy. Results indicate higher gas hydrate saturation in the upper part of the B1 sand and areas with higher acoustic impedance. Intrinsic permeability was modeled from the total porosity and clay-bound water volume, and effective permeability was derived from the gas hydrate saturation and intrinsic permeability distributions based on the “Tokyo model”. Effective permeability distributions were influenced by the total porosity, gas hydrate saturation, and intrinsic permeability. Within the same layer, higher gas hydrate saturation leads to decreased effective permeability. In total, 100 sets of multiple scenarios were prepared, providing input data for dynamic flow simulations to evaluate the effects of lateral heterogeneity in reservoir properties and the hydraulic characteristics of faults on production behavior for preassessment before the long-term production test.

58 GEOSCIENCES↗
DOE OSTI2021

Accessibility of Pores to Methane in New Albany Shale Samples of Varying Maturity Determined Using SANS and USANS

The accessibility of pores to methane has been investigated in Devonian New Albany Shale Formation early-mature (Ro = 0.50%) to post-mature (Ro = 1.40%) samples. A Marcellus Shale Formation sample was included to expand the maturation range to Ro 2.50%. These are organic matter-rich rocks with total organic carbon (TOC) values of 3.4 to 14.4% and porosity values of 2.19 to 6.88%. Contrast matching small-angle neutron scattering (SANS) and ultra-small angle neutron scattering (USANS) techniques were used to generate porosity-related data before and after pressure cycling under hydrostatic (in a vacuum and at 500 bar of deuterated methane) and uniaxial stress (0 to ca. 350 bar) conditions. Our results showed that the accessible porosity was small for the samples studied, ranging from zero to 2.9%. No correlation between the accessible porosity and TOC or mineralogical composition was revealed, and the most likely explanation for porosity variation was related to the thermal transformation of organic matter and hydrocarbon generation. Pressure caused improvements in accessible porosity for most samples, except the oil window sample (Ro = 0.84%). Our data show that densification of methane occurs in nanopores, generally starting at diameters smaller than 20 nm, and that the distribution of methane density is affected by pressure cycling.

03 NATURAL GAS↗
DOE OSTI2024

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↗
DOE OSTI2024

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↗
NASA NTRS1998

Hydroxyapatite fiber reinforced poly(alpha-hydroxy ester) foams for bone regeneration

A process has been developed to manufacture biodegradable composite foams of poly(DL-lactic-co-glycolic acid) (PLGA) and hydroxyapatite short fibers for use in bone regeneration. The processing technique allows the manufacture of three-dimensional foam scaffolds and involves the formation of a composite material consisting of a porogen material (either gelatin microspheres or salt particles) and hydroxyapatite short fibers embedded in a PLGA matrix. After the porogen is leached out, an open-cell composite foam remains which has a pore size and morphology defined by the porogen. By changing the weight fraction of the leachable component it was possible to produce composite foams with controlled porosities ranging from 0.47 +/- 0.02 to 0.85 +/- 0.01 (n = 3). Up to a polymer:fiber ratio of 7:6, short hydroxyapatite fibers served to reinforce low-porosity PLGA foams manufactured using gelatin microspheres as a porogen. Foams with a compressive yield strength up to 2.82 +/- 0.63 MPa (n = 3) and a porosity of 0.47 +/- 0.02 (n = 3) were manufactured using a polymer:fiber weight ratio of 7:6. In contrast, high-porosity composite foams (up to 0.81 +/- 0.02, n = 3) suitable for cell seeding were not reinforced by the introduction of increasing quantities of hydroxyapatite short fibers. We were therefore able to manufacture high-porosity foams which may be seeded with cells but which have minimal compressive yield strength, or low porosity foams with enhanced osteoconductivity and compressive yield strength.

NASA Discipline Cell Biology↗
DOE OSTI2020

Proposed Analytical Methods for Determining Filter Media Properties

High Efficiency Particulate Air (HEPA) filters, commonly used in nuclear filtration applications, are an integral part of the waste management processes in nuclear plants. HEPA filters are 99.97% efficient filtration devices characterized by their high resistance to air flow, or pressure drop. From theoretical models, the initial pressure drop across a clean filter is proven to be a function of the filter media fiber diameters and porosity of the media. These filter properties are relatively difficult to obtain with traditional manual methods; therefore, there is a need to develop analytical methods to find the fiber diameter and porosity to determine the pressure drop. Typically, these values can be found by substituting a calculated representative value based upon measured values, such as equivalent fiber diameter based upon a clean pressure drop. However, these values may also be directly recorded and measured by incorporating Scanning Electron Microscope (SEM) image analysis and other measurement methods to analyze the filter media and determine its physical properties without the need for media testing. DiameterJ, an open source Java plug-in, used with ImageJ, can process images of filter media taken by an SEM to find statistical data such as mean fiber diameter and porosity. To produce the raw data, SEM images of two filter media type samples are taken and segmented in DiameterJ using the traditional and statistical region merging segmentation algorithms. Manual segmentation is necessary after the initial segmentation by the algorithms as the images tend to be too complex for the algorithms to output with the necessary accuracy. However, complications exist in the manual segmentation process as these methods can be time intensive and prone to the individual bias of the user. This in turn can skew the final mean fiber diameter result, and lead to either an over or under prediction of the pressure drop. It was also discovered that the porosity data produced by DiameterJ is inaccurate, as the SEM analyzes the three-dimensional filter media by projecting its geometry onto a plane and analyzing it as a two-dimensional binary image. Thus, the porosity is artificially inflated through the segmentation process, rendering this result incorrect. Alternatively, density determination, gravimetric analysis, and thickness testing of the filter media is collectively used to determine the filter fiber porosity. Together, the SEM image analysis and analytical lab methods produce results through direct measurements which allow for the prediction of the initial pressure drop from clean filter media without the need for prior media testing to collect pressure data. By improving upon this proposed analytical method in the future, there is potential to streamline the process of finding these filter properties into a more direct methodology for determining the pressure drop across HEPA filters.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Oxidation and Associated Pore Structure Modification During Experimental Alteration of Granite

Weathering plays a crucial role in a number of environmental processes, and the microstructure and evolution of multi-scale pore space is a critically important factor in the weathering of geological formations. In igneous rocks the infiltration of meteoric water into initially relatively dry material can cause the rock to begin to disaggregate, increasing porosity and surface area, and allowing further disaggregation and weathering to occur. These processes, in turn, allow biota to colonize the rock, further enhancing the weathering rate. In some igneous rocks this may be driven by oxidation of primary minerals. One such mineral, biotite, has been repeatedly mentioned as a cause of cracking during oxidation. However, the scale-dependence of the processes by which this occurs are poorly understood. We cannot, therefore, accurately extrapolate laboratory reaction rates to the field in predictive numerical models. In order to better understand the effects of oxidation and test the hypothesis that fracture and disaggregation are initiated by swelling of oxidizing biotites, we reacted granite cores in a selenic acid-rich aqueous solution at 200°C for up to 438 days. Elevated temperatures and selenic acid were used to provide relatively fast reaction rates and highly oxidizing conditions in sealed reaction vessels. These experiments were analyzed using a combination of imaging, X-ray diffraction, Mössbauer spectroscopy, and small- and ultra-small angle neutron scattering to interrogate porosity and microfracture formation. The experimental results show little observable biotite swelling, but significant transport and growth of iron oxides and/or clays along grain boundaries throughout the sample. Significant increases in porosity were also observed at the sample rim, likely associated with feldspar alteration. Fractures and transport were observed throughout the core, suggesting that stresses due to crystallization pressures caused by the growing iron phases may be the initiating factors in granite weathering, possibly followed by biotite swelling after sufficient permeability is achieved.

Anovitz, Lawrence M.↗
DOE OSTI2017

Bayesian Inversion of Seismic and Electromagnetic Data for Marine Gas Reservoir Characterization Using Multi-chain Markov Chain Monte Carlo Sampling

In this study we developed an efficient Bayesian inversion framework for interpreting marine seismic amplitude versus angle (AVA) and controlled source electromagnetic (CSEM) data for marine reservoir characterization. The framework uses a multi-chain Markov-chain Monte Carlo (MCMC) sampler, which is a hybrid of DiffeRential Evolution Adaptive Metropolis (DREAM) and Adaptive Metropolis (AM) samplers. The inversion framework is tested by estimating reservoir-fluid saturations and porosity based on marine seismic and CSEM data. The multi-chain MCMC is scalable in terms of the number of chains, and is useful for computationally demanding Bayesian model calibration in scientific and engineering problems. As a demonstration, the approach is used to efficiently and accurately estimate the porosity and saturations in a representative layered synthetic reservoir. The results indicate that the seismic AVA and CSEM joint inversion provides better estimation of reservoir saturations than the seismic AVA-only inversion, especially for the parameters in deep layers. The performance of the inversion approach for various levels of noise in observational data was evaluated – reasonable estimates can be obtained with noise levels up to 25%. Sampling efficiency due to the use of multiple chains was also checked and was found to have almost linear scalability.

Ren, Huiying↗
DOE OSTI2018

Glass-bonded iodosodalite waste form for immobilization of 129I

Immobilization of radioiodine (e.g., 129I, 131I) is an important need for current and future nuclear fuel cycles. For the current work, iodosodalite [Na8(AlSiO4)6I2] was synthesized hydrothermally from metakaolin, NaI, and NaOH. Following hydrothermal treatment, dried unwashed powders were used to make glass-bonded iodosodalite waste forms by heating pressed pellets at 650, 750, or 850 °C with two different types of sodium borosilicate glass binders, i.e., NBS-4 and SA-800. These heat-treated specimens were characterized with X-ray diffraction, Fourier-transform infrared spectroscopy, scanning electron microscopy, energy dispersive spectroscopy, thermal analysis, porosity and density measurements, neutron activation analysis, and inductively-coupled plasma mass spectrometry. The pellets mixed with 10 mass% of NBS-4 or SA-800 and heat-treated at 750 °C contained relatively high percentage iodine retention (~44-47 % of the maximum iodine loading) with relatively low porosities, while other pellets with higher percentages iodine retention either contained higher porosity or were not completely sintered. ASTM C1308 chemical durability tests of monolithic specimens showed a large initial release of Na, Al, Si, and I on the first day, possibly from water-soluble salt crystals or non-durable amorphous phases. Release rates of Na and Si were higher than for Al and I, probably due to a poorly durable Na-Si-O phase from the glass bonding matrix. The cumulative normalized release of iodine was 12.5 g m-2 for the first 10 1-d exchanges, suggestive of coherent dissolution. The average release rate from 10-24 days during the 7-d exchange intervals was 0.2336 g m-2 d-1.

sodalite, iodosodalite, waste forms, iodine↗
DOE OSTI2024

Chemical composition based machine learning model to predict defect formation in additive manufacturing

With a goal of exploiting additive manufacturing to improve the manufacturing of existing reactor materials, we developed a chemical composition-based machine learning model to predict the printability of any given alloy in laser powder bed fusion (L-PBF) using experimental data from peer-reviewed literature. We defined printability as the ability to avoid defects like cracking, balling, porosity, and lack of fusion, that are caused by thermal stresses (during solidification or liquation), molten pool disintegration into disconnected small beads or lack of heat input respectively. Our models predict the tendency of balling defect formation and porosity percentage for a given composition, under a given set of processing conditions. To predict the likelihood of balling defect, three models: a random forest classifier, a gradient boost regressor and a neural network were trained on a dataset containing both traditional alloys and high entropy alloys. The neural network model showed the highest accuracy of 92.3 % in predicting the balling defect formation. A random forest regressor, gradient boost regressor and neural network were trained and tested on a dataset of various alloys to predict porosity. The random forest regressor showed the best predictions with an R 2 score of 0.97. The models also revealed the relative importance of the input descriptors on defect-formation tendency. Of particular significance was the identification of carbon as an important element in determining the occurrence of balling and percent porosity in alloys like steel, as well as being moderately important to the percentage porosity in other alloys as well as steel. Manganese was also identified as a key descriptor for the percentage of porosity in steel and other alloys. Manganese’s low thermal conductivity and consistent presence in the dataset is the likely cause for its contribution. Carbon’s role is attributable to its relatively high specific heat and high melting temperature. In conclusion, our model serves as a swift, chemistry-based tool to design experiments and find modified compositions better suited for additive manufacturing.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Morphology-Driven Control of Metabolite Selectivity Using Nanostructure-Initiator Mass Spectrometry

Nanostructure-initiator mass spectrometry (NIMS) is a laser desorption/ionization analysis technique based on the vaporization of a nanostructure-trapped liquid "initiator" phase. Here in this work, we report an intriguing relationship between NIMS surface morphology and analyte selectivity. Scanning electron microscopy and spectroscopic ellipsometry were used to characterize the surface morphologies of a series of NIMS substrates generated by anodic electrochemical etching. Mass spectrometry imaging was applied to compare NIMS sensitivity of these various surfaces toward the analysis of diverse analytes. The porosity of NIMS surfaces was found to increase linearly with etching time where the pore size ranged from 4 to 12 nm with corresponding porosities estimated to be 7-70%. Surface morphology was found to significantly and selectively alter NIMS sensitivity. The small molecule (<2k Da) sensitivity was found to increase with increased porosity, whereas low porosity had the highest sensitivity for the largest molecules examined. Estimation of molecular sizes showed that this transition occurs when the pore size is <3× the maximum of molecular dimensions. While the origins of selectivity are unclear, increased signal from small molecules with increased surface area is consistent with a surface area restructuring-driven desorption/ionization process where signal intensity increases with porosity. In contrast, large molecules show highest signal for the low-porosity and small-pore-size surfaces. We attribute this to strong interactions between the initiator-coated pore structures and large molecules that hinder desorption/ionization by trapping large molecules. Lastly, this finding may enable us to design NIMS surfaces with increased specificity to molecules of interest.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2014

Hierarchical Pore Morphology of Cretaceous Shale: A Small-Angle Neutron Scattering and Ultrasmall-Angle Neutron Scattering Study

Shale reservoirs are becoming an increasingly important source of oil and natural gas supply and a potential candidate for CO 2 sequestration. Understanding the pore morphology in shale may provide clues to making gas extraction more efficient and cost-effective. The porosity of Cretaceous shale samples from Alberta, Canada, collected from different depths with varying mineralogical compositions, has been investigated by small- and ultrasmall-angle neutron scattering. Moreover these samples come from the Second White Specks and Belle Fourche formations, and their organic matter content ranges between 2 and 3%. The scattering length density of the shale specimens has been estimated using the chemical composition of the different mineral components. Scattering experiments reveal the presence of fractal and non-fractal pores. It has been shown that the porosity and specific surface area are dominated by the contribution from meso- and micropores. The fraction of closed porosity has been calculated by comparing the porosities estimated by He pycnometry and scattering techniques. There is no correlation between total porosity and mineral components, a strong correlation has been observed between closed porosity and major mineral components in the studied specimens.

58 GEOSCIENCES↗
DOE OSTI2023

Using small angle x-ray scattering to examine the aggregation mechanism in silica nanoparticle-based ambigels for improved optical clarity

Silica-based aerogels are a promising low-cost solution for improving the insulation efficiency of single-pane windows and reducing the energy consumption required for space heating and cooling. Two key material properties required are high porosity and small pore sizes, which lead to low thermal conductivity and high optical transparency, respectively. However, porosity and pore size are generally directly linked, where high porosity materials also have large pore sizes. This is unfavorable as large pores scatter light, resulting in reduced transmittance in the visible regime. In this work, we utilized preformed silica colloids to explore methods for reducing pore size while maintaining high porosity. Further, the use of preformed colloids allows us to isolate the effect of solution conditions on porous gel network formation by eliminating simultaneous nanoparticle growth and aggregation found when using typical sol–gel molecular-based silica precursors. Specifically, we used in situ synchrotron-based small-angle x-ray scattering during gel formation to better understand how pH, concentration, and colloid size affect particle aggregation and pore structure. Ex situ characterization of dried gels demonstrates that peak pore widths can be reduced from 15 to 13 nm, accompanied by a narrowing of the overall pore size distribution, while maintaining porosities of 70%–80%. Optical transparency is found to increase with decreasing pore sizes while low thermal conductivities ranging from 95 +/− 13 mW/m K are maintained. Mechanical performance was found to depend primarily on effective density and did not show a significant dependence on solution conditions. Overall, our results provide insights into methods to preserve high porosity in nanoparticle-based aerogels while improving optical transparency.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Direct Numerical Simulation of the Flow through a Randomly Packed Pebble Bed

The proposition for molten salt and high-temperature gas-cooled reactors has increased the focus on the dynamics and physics in randomly packed pebble beds. Research is being conducted on the validity of these designs as a possible contestant for the fourth-generation nuclear power systems. A detailed understanding of the coolant flow behavior is required in order to ensure proper cooling of the reactor core during normal and accident conditions. In order to increase the understanding of the flow through these complex geometries and enhance the accuracy of lower-fidelity modeling, high-fidelity approaches such as direct numerical simulation (DNS) can be utilized. Nek5000, a spectral-element computational fluid dynamics (CFD) code, was used to develop DNS fluid flow data. The flow domain consisted of 147 pebbles enclosed by a bounding wall. In the work presented, the Reynolds numbers ranged from 430 to 1050 based on the pebble diameter and inlet velocity. Characteristics of the flow domain such as volume averaged porosity, axial porosity, and radial porosity were studied and compared with correlations available in the literature. Friction factors from the DNS results for all Reynolds numbers were compared with correlations in the literature. The first- and second-order statistics show good agreement with the available experimental data. Turbulence length scales were analyzed in the flow. Reynolds stress anisotropy was characterized by utilizing invariant analysis. Overall, the results of the analysis in this study provide deeper understanding of the flow behavior and the effect of the wall in packed beds.

Yildiz, Mustafa Alper↗
DOE OSTI2014

Gas Flow Tightly Coupled to Elastoplastic Geomechanics for Tight- and Shale-Gas Reservoirs: Material Failure and Enhanced Permeability

We investigate coupled flow and geomechanics in gas production from extremely low permeability reservoirs such as tight and shale gas reservoirs, using dynamic porosity and permeability during numerical simulation. In particular, we take the intrinsic permeability as a step function of the status of material failure, and the permeability is updated every time step. We consider gas reservoirs with the vertical and horizontal primary fractures, employing the single and dynamic double porosity (dual continuum) models. We modify the multiple porosity constitutive relations for modeling the double porous continua for flow and geomechanics. The numerical results indicate that production of gas causes redistribution of the effective stress fields, increasing the effective shear stress and resulting in plasticity. Shear failure occurs not only near the fracture tips but also away from the primary fractures, which indicates generation of secondary fractures. These secondary fractures increase the permeability significantly, and change the flow pattern, which in turn causes a change in distribution of geomechanical variables. From various numerical tests, we find that shear failure is enhanced by a large pressure drop at the production well, high Biot's coefficient, low frictional and dilation angles. Smaller spacing between the horizontal wells also contributes to faster secondary fracturing. When the dynamic double porosity model is used, we observe a faster evolution of the enhanced permeability areas than that obtained from the single porosity model, mainly due to a higher permeability of the fractures in the double porosity model. These complicated physics for stress sensitive reservoirs cannot properly be captured by the uncoupled or flow-only simulation, and thus tightly coupled flow and geomechanical models are highly recommended to accurately describe the reservoir behavior during gas production in tight and shale gas reservoirs and to smartly design production scenarios.

02 PETROLEUM↗
DOE OSTI2014

Modeling CO 2 -Brine-Rock Interaction Including Mercury and H 2 S Impurities in the Context of CO 2 Geologic Storage

This study uses modeling and simulation approaches to investigate the impacts on injectivity of trace amounts of mercury (Hg) in a carbon dioxide (CO 2 ) stream injected for geologic carbon sequestration in a sandstone reservoir at ~2.5 km depth. At the range of Hg concentrations expected (7-190 ppbV, or ~ 0.06-1.6 mg/std.m 3 CO 2 ), the total volumetric plugging that could occur due to complete condensation of Hg, or due to complete precipitation of Hg as cinnabar, results in a very small porosity change. In addition, Hg concentration much higher than the concentrations considered here would be required for Hg condensation to even occur. Concentration of aqueous Hg by water evaporation into CO 2 is also unlikely because the higher volatility of Hg relative to H 2 O at reservoir conditions prevents the Hg concentration from increasing in groundwater as dry CO 2 sweeps through, volatilizing both H 2 O and Hg. Using a model-derived aqueous solution to represent the formation water, batch reactive geochemical modeling show that the reaction of the formation water with the CO 2 -Hg mixture causes the pH to drop to about 4.7 and then become buffered near 5.2 upon reaction with the sediments, with a negligible net volume change from mineral dissolution and precipitation. Cinnabar (HgS(s)) is found to be thermodynamically stable as soon as the Hg-bearing CO 2 reacts with the formation water which contains small amounts of dissolved sulfide. Liquid mercury (Hg(l)) is not found to be thermodynamically stable at any point during the simulation. Two-dimensional radial reactive transport simulations of CO 2 injection at a rate of 14.8 kg/s into a 400 m-thick formation at isothermal conditions of 106°C and average pressure near 215 bar, with varying amounts of Hg and H 2 S trace gases, show generally that porosity changes only by about ±0.05% (absolute, i.e., new porosity = initial porosity ±0.0005) with Hg predicted to readily precipitate from the CO 2 as cinnabar in a zone mostly matching the single-phase CO 2 plume. The precipitation of minerals other than cinnabar, however, dominates the evolution of porosity. Main reactions include the replacement of primarily Fe-chlorite by siderite, of calcite by dolomite, and of K-feldspar by muscovite. Chalcedony is also predicted to precipitate from the dissolution of feldspars and quartz. Although the range of predicted porosity change is quite small, the amount of dissolution and precipitation predicted for these individual minerals is not negligible. These reactive transport simulations assume that Hg gas behaves ideally. To examine effects of non-ideality on these simulations, approximate calculations of the fugacity coefficient of Hg in CO 2 were made. Results suggest that Hg condensation could be significantly overestimated when assuming ideal gas behavior, making our simulation results conservative with respect to impacts on injectivity. The effect of pressure on Henry’s constant for Hg is estimated to yield Hg solubilities about 10% lower than when this effect is not considered, a change that is considered too small to affect the conclusions of this report. Although all results in this study are based on relatively mature data and modeling approaches, in the absence of experimental data and more detailed site-specific information, it is not possible to fully validate the results and conclusions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Geodyn Material Library: Pseudocap models for dry porous tocks

This report describes the second edition of the Pseudocap Strength models for porous rocks implemented in GEODYN material library. The first model was developed in 2007 and calibrated for concrete. Then, the model parameters were calibrated based on triaxial tests reported for limestones and sandstones of various porosities. In these models some key parameters were chosen as functions of the reference porosity,Φ. Since then multiple modifications were implemented in the model, therefore, it has been recalibrated for some common porous materials such as limestones, sandstones, alluvium, tuffs and granite. Two types of models are described in this repot. The first type (called Pseudocap Model or PM) is for rocks from a specific location. Parameters were calibrated for several specific geologic materials. The second type (called Generic Pseudocap Model or GPM) is useful for the sites where only basic information (rock type, porosity) is available. Generic models include built-in correlations between porosities and other mechanical properties observed for certain rock types. The models of both types were validated by comparing not only to quasi-static triaxial tests for these materials but also to shock Hugoniot data and spherical explosion data for some materials. All models were derived in the frame of isotropic plasticity. They are designed to be used in explicit finite element/difference codes. Tangent stiffness tensor is not provided but can be calculated numerically for the model to be used in implicit finite element codes. For an isotropic material the stress can be decomposed into volumetric and deviatoric parts. The volumetric part is modeled using an Equation of state (EOS) which calculates the pressure and the bulk sound speed as functions of the internal specific energy and density. Here a simple, Mie-Gruneisen EOS is presented, but tabulated EOS (LEOS) provided by the library can be used as well. The stress is limited by the yield surface which depends on three invariants of the stress tensor and specific internal energy. In addition, to capture the strain-rate dependence a simple multiplier is used for the yield surface which depends on the equivalent plastic strain rate. The failure surface (the ultimate yield, Y f , defined later) is chosen in the Hoek-Brown form, commonly used in rock mechanics. It includes measurable parameters such as Unconfined Compressive Strength (UCS) as well as scale parameters characterizing the quality of the rock such as GSI (Geologic Strength Index). Thus, even though the model is calibrated for small samples it offers a way to extrapolate the strength to the field scale using geological characterization of the rock mass. The model captures effects of brittle-ductile transition in rocks by introducing a cap multiplier to the yield function. The rate of dilatancy (bulking) is proportional to the slope of the yield surface affected by the cap. Therefore, it takes place only at low confinements when the pressure is less than the brittle-ductile transition pressure, P BD . On the contrary, the porous compaction takes place at pressures higher than P BD . The cap moves as the porosity is compacted or new porosity is generated due to dilatancy. The porous compaction is modeled using an evolution equation which includes deviatoric stress so that the onset of compaction corresponds to the cap surface. The model captures effects of shear-enhanced compaction which is an important for porous rocks. Section 2 describes the modeling framework and Section 3 presents the model calibration procedure. Section 4 compares experimental data for various rocks versus model predictions. The model parameters used for this comparison are given in Appendix. The files with material constants are available with the latest GEODYN material library distribution.

58 GEOSCIENCES↗