Search NASA⌕ Search

SEARCH · Search NASA

Results for “Volume fractions”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Deciphering the small-angle scattering of polydisperse hard spheres using deep learning

We introduce a deep learning approach for analyzing the scattering function of the polydisperse hard sphere system. We use a variational autoencoder-based neural network to learn the bidirectional mapping between the scattering function and the system parameters, including the volume fraction and polydispersity. Such that the trained model serves both as a generator that produces a scattering function from the system parameters and an inferrer that extracts system parameters from the scattering function. We first generate a scattering dataset by carrying out molecular dynamics simulations of the polydisperse hard spheres modeled by the truncated-shifted Lennard-Jones model, then analyze the scattering function dataset using singular value decomposition to confirm the feasibility of dimensional compression. Then, we split the dataset into training and testing sets and train our neural network on the training set only. Our generator model produces a scattering function with significantly higher accuracy compared to the traditional Percus–Yevick approximation and β correction, and the inferrer model can extract the volume fraction and polydispersity with much higher accuracy than traditional model functions.

Ding, Lijie [ORNL] (ORCID:0000000227454606)↗
DOE OSTI2021

Two-color pyrometry based flame to fuel surface radiative heat flux diagnostic using flamelets

In this study, a non-intrusive two-color pyrometry based radiative heat flux diagnostic is presented that accounts for the self-absorption effects. This diagnostic relies on looking up the flame temperature and soot volume fraction using the ratio of color intensities obtained from digital cameras and numerical solutions of 1D steady-state diffusion flames. Virtual two color pyrometry of computed flame show a unique one-to-one mapping of strained flames with intensity ratio; thereby allowing flames to be indexed and creation of a two-color flamelet manifold (TCFM). The TCFM is a function of two variables; the ratio of red to green intensities (flame index) and physical distance to the stoichiometric surface defined using a level-set function. The developed diagnostic is used in upward flame spread experiments where an approximate flame hull reconstruction method is applied to provide the 3D level-set function and intensity ratio for use of TCFM to provide radiative absorption properties for 3D ray tracing. Accounting for self- absorption effects provide improved soot volume fraction estimates. In conclusion, the radiative heat flux estimates are found to be in very good agreement with literature and improved results are obtained near the flame base.

42 ENGINEERING↗
DOE OSTI2022

Influence of intermolecular interactions on the infrared complex indices of refraction for binary liquid mixtures

This paper investigates the accuracy of deriving the composite optical constants of binary mixtures from only the complex indices of refraction of the neat materials. These optical constants enable the reflectance spectra of the binary mixtures to be modeled for multiple scenarios (e.g., different substrates, thicknesses, volume ratios), which is important for contact and standoff chemical detection. Using volume fractions, each mixture's complex index of refraction was approximated via three different mixing rules. To explore the impact of intermolecular interactions, these predictions are tested by experimental measurements for two representative sets of binary mixtures: (1) tributyl phosphate combined with n-dodecane, a non-polar medium, to represent mixtures which primarily interact via dispersion forces and (2) tributyl phosphate and 1-butanol to represent mixtures with polar functional groups that can also interact via dipole–dipole interactions, including hydrogen bonding. The residuals and the root-mean-square error between the experimental and calculated index values are computed and demonstrate that for miscible liquids in which the average geometry of the cross-interactions can be considered isotropic (e.g., dispersion), the refractive indices of the mixtures can be modeled using composite n and k values derived from volume fractions of the neat liquids. Conversely, in spectral regions where the geometry of the cross-interactions is more restricted and anisotropic (e.g., hydrogen bonding), the calculated n and k values vary from the measured values. The impact of these interactions on the reflectance spectra are then compared by modeling a thin film of the binary mixtures on an aluminum substrate using both the measured and the mathematically computed indices of refraction.

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

Bulk high-temperature superconductivity in pressurized tetragonal La 2 PrNi 2 O 7

The Ruddlesden–Popper (R–P) bilayer nickelate, La 3 Ni 2 O 7 , was recently found to show signatures of high-temperature superconductivity (HTSC) at pressures above 14 GPa . Subsequent investigations achieved zero resistance in single-crystalline and polycrystalline samples under hydrostatic pressure conditions. Yet, obviousdiamagnetic signals, the other hallmark of superconductors, are still lacking owing to the flamentary nature with low superconducting volume fraction. The presence of a new 1313 polymorph and competing R–P phases obscured proper identification of the phase for HTSC. Thus, achieving bulk HTSC and identifying the phase at play are the most prominent tasks. Here we address these issues in the praseodymium (Pr)-doped La 2 PrNi 2 O 7 polycrystalline samples. We find that substitutions of Pr for La efectively inhibit the intergrowth of diferent R–P phases, resulting in a nearly pure bilayer structure. For La 2 PrNi 2 O 7 , pressure-induced orthorhombic to tetragonalstructural transition takes place at P c ≈ 11 GPa, above which HTSC emerges gradually on further compression. The superconducting transition temperatures at 18–20 GPa reach $T$ $^{onset}_{c}$ = $82.5$ $K$ and $T$ $^{zero}_{c}$ = $60$ $K$, which are the highest values, to our knowledge, among known nickelate superconductors. Importantly, bulk HTSC was testified by detecting clear diamagnetic signals below about 75 K with appreciable superconducting shielding volume fractions at a pressure of above 15 GPa. Further, our results not only resolve the existing controversies but also provide directions for exploring bulk HTSC in the bilayer nickelates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Fluctuations and power-law scaling of dry, frictionless granular rheology near the hard-particle limit

Here the flow of frictionless granular particles is studied with stress-controlled discrete element modeling simulations for systems varying in size from 300 to 100,000 particles. The volume fraction and shear-stress ratio μ are relatively insensitive to system size for a wide range of inertial numbers I. Second-order effects in strain rate, such as normal stress differences, require large system sizes to accurately extract meaningful results, notably a nonmonotonic dependence in the first normal stress difference with strain rate. The rheological response represented by the μ(I) scalar model works well at describing the lower-order aspects of the rheology, except near the quasistatic limit of these stress-controlled flows. The pressure is varied over five decades, and a pressure dependence of the coordination number is observed, which is not captured by the inertial number. Large fluctuations observed for small systems N≤1000 near the quasistatic limit can lead to the arrest of flow resulting in challenges to fitting the data to rheological relationships. The inertial number is also insufficient for capturing the pressure-dependent behavior of property fluctuations. Fluctuations in the flow and microstructural properties are measured in both the quasistatic and inertial regimes, including shear stress, pressure, strain rate, normal stress differences, volume fraction, coordination number, and contact fabric anisotropy. The fluctuations in flow properties scale self-similarly with pressure and system size. A transition in the scaling of fluctuations of stress properties and contact fabric anisotropy are measured and proposed as a quantitative identification of the transition from inertial to quasistatic flow.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Soot-particle core-shell and fractal structures from small-angle X-ray scattering measurements in a flame

We have characterized soot particles measured in situ in a laminar co-flow ethylene-air diffusion flame using small-angle X-ray scattering (SAXS). The analysis includes temperature measurements made with coherent anti-Stokes Raman spectroscopy (CARS) and complements soot volume-fraction and maturity measurements made with laser-induced incandescence (LII). We compared the results of fits to the SAXS measurements using a unified model and a fractal core-shell model. Power-law parameters yielded by the unified model indicate that aggregates of primary particles are in the mass-fractal regime, whereas the primary particles are in the surface-fractal regime in the middle of the flame. Higher and lower in the flame, the primary-particle power-law parameter approaches 4, suggesting smooth primary particles. Furthermore, these trends are consistent with fits using the fractal core-shell model, which indicate that particles have an established core-shell structure in the middle of the flame and are internally homogeneous at higher and lower heights in the flame. Primary-particle size distributions derived using the fractal core-shell model demonstrate excellent agreement with distributions inferred from transmission electron microscopy (TEM) images in the middle of the flame. Higher in the flame, a second small mode appears in the size distributions, suggesting particle fragmentation during oxidation. Surface oxidation would explain (1) aggregate fragmentation and (2) loss of core-shell structure leading to smoother primary-particle surfaces by removal of carbon overlayers. SAXS measurements are much more sensitive to incipient and young soot particles than LII and demonstrate significant volume fraction from particles low in the flame where the LII signal is negligible.

42 ENGINEERING↗
DOE OSTI2023

Optimizing Phase Change Composite Thermal Energy Storage Using the Thermal Ragone Framework

The thermal Ragone framework can be applied to the design of thermal energy storage heat exchangers to identifying relationships between their power requirements and energy capabilities. This presentation focuses on the design of planar thermal energy storage heat exchangers with phase change materials and thermal conductivity additives. Finite-difference models were used to optimize thermal storage heat exchanger designs focusing on key design parameters of conductivity additive volume fraction, spacing between tubes, and the phase transition temperatures. Design trends were developed for optimal designs based on varying thermal power requirements. To aid in future device design, models of reduced complexity were developed with the ability to predict optimal designs with computational speeds up to 5,000 times faster than finite-difference models.

ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATION,↗
DOE OSTI2024

Automated Manufacturing of Grid Stiffened Panels with Radically Reduced Tooling

Grid stiffened continuous fiber reinforced composite panels are an attractive option for creating lightweight structures due to the tailorability for various applications and the resulting high specific properties. However, the panel stiffeners and stiffener intersections result in high tooling complexity and correspondingly high cost of implementation. These factors have limited the impact of such structures in the composites industry. Previous research has demonstrated the ability to produce high quality, high aspect ratio beams, representative of individual grid stiffeners, using E-glass/PET comingled tow via direct digital manufacturing. Further, prior preliminary efforts have demonstrated the potential to use the same approach to manufacture grid intersections that have continuous fiber in both directions. To expand on the previous efforts in grid stiffeners produced by direct digital manufacture with radically reduced tooling requirements, this effort compares two methods of providing positioning and consolidation, nozzle vs. roller. Both processes are based on a commingled yarn feedstock. The extrusion through a nozzle has been shown to enable grid intersection control through local variations in applied consolidation and serves as the baseline process. However, this approach requires a continuous placement path to create the complete grid stiffened panel as no mechanism for cutting and restarting has been implemented. Alternatively, a newly developed placement head incorporating cut and refeed, mounted to a 6-axis robot, offers the potential of improved path placement efficiency. The two techniques are used to produce similar grid composite stiffeners to evaluate the effectiveness of producing the grid intersections. Rate of deposition of the two end effectors are compared and the quality of the associated grid stiffeners, and intersections, are determined through measurement of geometry, fiber volume fraction and void fraction.

Engineering↗
DOE OSTI2020

Numerical predictions of cryogenic hydrogen vertical jets

Here, a comparison of Computational Fluid Dynamics (CFD) predictions with measurements is presented for cryo-compressed hydrogen vertical jets. The stagnation conditions of the experiments are characteristic of unintended leaks from pipe systems that connect cryogenic hydrogen storage tanks and could be encountered at a fuel cell refueling station. Jets with pressure up to 5 bar and temperatures just above the saturation liquid temperature were examined. Comparisons are made to the centerline mass fraction and temperature decay rates, the radial profiles of mass fraction and the contours of volume fraction. Two notional nozzle approaches are tested to model the under-expanded jet that was formed in the tests with pressures above 2 bar. In both approaches the mass and momentum balance from the throat to the notional nozzle are solved, while the temperature at the notional nozzle was assumed equal to the nozzle temperature in the first approach and was calculated by an energy balance in the second approach. The two approaches gave identical results. Satisfactory agreement with the measurements was found in terms of centerline mass fraction and temperature. However, for test with 3 and 4 bar release the concentration was overpredicted. Furthermore, a wider radial spread was observed in the predictions possibly revealing higher degree of diffusion using the k-ε turbulence model. An integral model for cryogenic jets was also developed and provided good results. Finally, a test simulation was performed with an ambient temperature jet and compared to the cold jet showing that warm jets decay faster than cold jets.

08 HYDROGEN↗
DOE OSTI2023

Comparative Life Cycle Assessment of Injection Molded and Big Area Additive Manufactured NdFeB Bonded Permanent Magnets

Permanent magnets are expected to play a crucial role in the realization of the clean economy. In particular, the neodymium-iron-boron (Nd2Fe14B or NdFeB) magnets, which have the highest energy density among rare earth permanent magnets, are needed for building more efficient windmill generators, electric vehicle motors, etc. Currently, near-net shape magnets can be either made through sintering and compression molding with extensive post machining or directly through injection molding. However, injection molding has a loading volume fraction limitation of 0.65 for nylon binders. A novel method of manufacturing bonded permanent magnets with loading fraction greater than 0.65 has been demonstrated using Big Area Additive Manufacturing (BAAM) printers. As energy density is directly proportional to the square of the magnet loading fraction, magnets produced using BAAM printers require less volume and magnetic material compared to that of injection molded magnets on average. A comparative LCA shows that this difference in magnetic powder consumption nearly constitutes the difference in the environmental impact categories. Even after assuming recycled magnetic input, the BAAM magnets perform better environmentally than injection molded magnets, especially in the ozone depletion category. Since BAAM printers can accommodate even higher loading fractions, at scale, BAAM printers possibly can bring about a significant decrease in rare earth mineral consumption and environmental emissions. So furthermore, single screw extrusion enables BAAM printers to have high print speeds and allow them to be economically competitive against injection molding. Therefore, BAAM printed magnets show promise in transitioning towards the clean economy.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Dimensionality-Induced Change in Topological Order in Multiferroic Oxide Superlattices

We construct ferroelectric (LuFeO 3 ) m /(LuFe 2 O 4 ) superlattices with varying index m to study the effect of confinement on topological defects. We observe a thickness-dependent transition from neutral to charged domain walls and the emergence of fractional vortices. In thin LuFeO 3 layers, the volume fraction of domain walls grows, lowering the symmetry from P6 3 cm to P3c1 before reaching the nonpolar P6 3 /mmc state, analogous to the group-subgroup sequence observed at the high-temperature ferroelectric to paraelectric transition. Our study shows how dimensional confinement stabilizes textures beyond those in bulk ferroelectric systems.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Quality Investigation of Pultruded Carbon Fiber Panels Subjected to Four-Point Flexure via Fiber Optic Sensing

Pultruded carbon fiber-reinforced composites are attractive to the wind energy industry due to the rapid production of highly aligned unidirectional composites with enhanced fiber volume fractions and increased specific strength and stiffness. However, high volume carbon fiber manufacturing remains cost-prohibitive. This study investigates the feasibility of a pultruded low-cost textile carbon fiber-reinforced epoxy composite as a promising material in spar cap production was undertaken based on mechanical response to four-point flexure loading. As spar caps are primarily subjected to flexural loading, large-span four-point flexure was considered, and coupon testing was restricted to tensile modulus and compression strength assessment. High-resolution spatial fiber optic strain sensing was utilized to determine spatial strain distribution during four-point flexure, revealing consistent strain along the length of the part and proved to be an excellent option for process manufacturing quality examination. Additionally, holes with diameters of 2.49 mm, 5.08 mm, and 1.93 mm were drilled through the thickness of full-width parts to determine the feasibility of structural health monitoring of pultruding parts internal to wind blades via fiber optic strain sensing.

Chemistry↗
DOE OSTI2024

Globally Gridded Groundwater Extraction Volumes and Costs under Six Depletion and Ponded Depth Targets

This repository contains simulated outputs from superwell – a hydro-economic tool for long-term assessment of groundwater cost and supply – providing globally gridded groundwater extractable volumes and associated unit costs ($/km³) for accessible groundwater production, based on a variety of user-defined depletion and ponded depth scenarios. Key model documentation: 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 Find the source code of the superwell model on GitHub: https://github.com/JGCRI/superwell Repository Overview Main output: superwell_outputs.7z contains 6 files (4.5 GB) named as superwell_py_deep_all_0.*PD_0.*DL.csv. These files present superwell outputs of global groundwater extraction volumes and cost estimates on a 0.5° scale for six scenarios with different Ponded Depth (PD; 0.3 and 0.6 m) and Depletion Limit (DL; 5%, 25%, and 40% of available volume) targets over the entire pumping lifetime of a grid cell superwell_py_deep_all_0.3PD_0.25DL_sample_100.csv contains superwell outputs for 100 data points sampled to match the global inputs' distribution superwell_py_deep_all_0.3PD_0.25DL_Grid_72548.csv contains superwell output for a single grid cell concept_v5.png provides an overview of the superwell workflow Outputs Description year_number: year of pumping depletion_limit: set depletion limit (DL) as a volume fraction of total available groundwater Mappings: continent, country, gcam_basin_id, Basin_long_name, grid_id: geographic identifiers and basin information Inputs: grid_area (km²): area of the grid cell whyclass: hydrogeological classification of the aquifer permeability (m/day), porosity (%), total_thickness (m), depth_to_water (m): aquifer properties. The geo-processed input data has been published separately: https://doi.org/10.57931/2307831 Model outputs: orig_aqfr_sat_thickness (m), aqfr_sat_thickness (m): original and remaining/instantaneous saturated thickness of the aquifer hydraulic_conductivity (m/day), transmissivity (m²/day): hydraulic properties of the aquifer radius_of_influence (m), areal_extent (km²): well radius and area of influence from the center of the well number_of_wells (-): number of wells in a grid cell determined by a ratio of well area and grid area max_drawdown (m), drawdown (m), drawdown_interference (m): well and aquifer drawdown during extraction total_head (m): total lift for the groundwater (depth to water plus drawdown) total_well_length (m): total depth of wells drilled well_yield (m³/day): pumping rate or well yield power (kW), energy (kWh): power and energy required for pumping groundwater Volume Outputs: volume_produced_perwell (m³), cumulative_vol_produced_perwell (m³): production volume metrics per well volume_produced_allwells (m³), cumulative_vol_produced_allwells (m³): aggregate extraction volumes for all wells in a grid cell available_volume (m³): available groundwater in storage for the grid cell as determined by aquifer properties depleted_vol_fraction: fraction of total volume pumped over available volumes in a grid cell (same as depletion limit) Cost Outputs: well_installation_cost ($): well installation cost based on the hydrogeological complexity of the aquifer annual_capital_cost, maintenance_cost, nonenergy_cost ($): nonenergy costs energy_cost_rate ($/kWh): electricity rate energy_cost ($): energy cost of pumping groundwater total_cost_perwell ($), total_cost_allwells ($): total annual energy and non-energy cost for each and all wells in a grid cell a unit_cost ($/m³), unit_cost_per_km3 ($/km³), unit_cost_per_acreft ($/acre-ft): total cost of pumping a unit of groundwater, indicated for different spatial units Key Resources Model documentation: 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. Input data: 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 superwell source code: https://github.com/JGCRI/superwell Cite as Niazi, H., Ferencz, S., Yoon, J., Graham, N., Wild, T., Hejazi, M., Watson, D., & Vernon, C. (2024). Globally Gridded Groundwater Extraction Volumes and Costs under Six Depletion and Ponded Depth Targets. MSD-LIVE Data repository. https://doi.org/10.57931/2307832 Contact Reach out to Hassan Niazi or Stephen Ferencz or open an issue in the superwell repository for questions or suggestions.

Earth Systems↗
DOE OSTI2021

Origins of Lithium/Sodium Reverse Permeability Selectivity in 12-Crown-4-Functionalized Polymer Membranes

Direct lithium extraction via membrane separations has been fundamentally limited by lack of monovalent ion selectivity exhibited by conventional polymeric membranes, particularly between sodium and lithium ions. Recently, a 12- Crown-4-functionalized polynorbornene membrane was shown to have the largest lithium/sodium permeability selectivity observed in a fully aqueous system to date. Using atomistic molecular dynamics simulations, we reveal that this selectivity is due to strong interactions between sodium ions and 12-Crown-4 moieties, which reduce sodium ion diffusivity while leaving lithium ion mobility relatively unaffected. Furthermore, the ion diffusivities in the membrane, when scaled by their respective solution diffusivities and free ion fractions, can be collapsed to an almost universal relationship depending on solvent volume fraction.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Percolation and nanosecond fluctuators in V 2 O 3 films within the metal–insulator transition

Vanadium sesquioxide (V 2 O 3 ) exhibits a metal–insulator transition (MIT) at 160 K between a low temperature, monoclinic, antiferromagnetic Mott insulator and a high temperature, rhombohedral, paramagnetic, metallic phase. In thin films, a percolative transition takes place over a finite temperature range of phase coexistence. We study the fluctuating dynamics of this percolative MIT by measuring voltage noise spectra at both low frequencies (up to 100 kHz) and radio frequencies (between 10 MHz and 1 GHz). Noise intensity quadratic in bias is observed in the MIT region, as expected for resistive fluctuations probed nonperturbatively by the current. The low frequency noise resembles flicker-type 1/f β noise, often taking on the form of Lorentzian noise dominated by a small number of fluctuators as the volume fraction of the insulating phase dominates. Radio frequency noise intensity also quadratic in the bias current allows the identification of resistance fluctuations with lifetimes below 1 ns, approaching timescales seen in non-equilibrium pump–probe studies of the transition. We find quantitative consistency with a model for fluctuations in the percolative fraction. The thermodynamics of the MIT suggests that dominant fluctuations are ones that alter small volumes affecting the connectivity of domain boundaries. This noise serves as a sensitive and nonperturbative probe for the dynamics of switching phenomena in this system.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Insights into Controls of Mineralogy and Pore Structure on the Density of Methane Adsorption Phase in Shales under Supercritical Conditions

With gradually increasing enhancing population and industrialization, exploitation of shale gas is enhancing in the world to satisfy the growing demand for energy worldwide. To exploit shale gas efficiently, accurate assessment of shale gas in place (GIP) is necessary for determining production strategies. Unlike free and dissolved gas, adsorbed gas contributes to the shale GIP up to 85% due to the well-developed surface area and micropores. Measured excess adsorption of methane (up to over 95% of subsurface shale gas) is often corrected to absolute adsorption to obtain the actual amount of adsorption under supercritical (or geological) conditions. During the correction, adsorption phase density (APD) is critical. However, the APD of subsurface shale gas and the effects of shale properties (e.g., mineralogy and pore structure) on APD remain poorly understood. A series of high-pressure methane adsorption isotherms (HPMAI) on Caney Shales were collected and analyzed in conjunction with other United States and Chinese shales from the literature at 35–125 °C and up to 15 MPa. Here, a three-layer Ono–Kondo (OK3) model is utilized to derive the temperature- and pressure-dependent APD coupling low-pressure nitrogen adsorption isotherms (LPNAI) and HPMAI. X-ray diffraction and organic geochemistry are combined to reveal the mineralogy. Brunauer–Emmett–Teller, Barrett–Joyner–Halenda, and Horvath–Kawazoe analyses via LPNAI are used to investigate the pore structures. Results show that APD increases with organic matter (OM) proxied by total organic carbon and decreases with clay minerals and the sum of quartz and feldspar. OM dramatically contributes to the APD as multiple-layer adsorption exists, and the APD for OM could be 1.4–8.5 times that for clay minerals. Other inorganic minerals contribute less to APD. The properties and constitution of the surface area instead of the volume fraction contribute to the APD in shales. APD does not show an obvious correlation with micropore volume, likely related to the ratio of micropore volume to the total pore volume. Here, we provide a significant and comprehensive study of petrological factors that impact the APD of subsurface shale gas, which will improve the estimation of supercritical adsorption and shale GIP under reservoir conditions. Also, the findings in this work can provide applications for subsurface carbon dioxide adsorption and storage.

58 GEOSCIENCES↗
DOE OSTI2025

Analysis of the weighted shifted boundary method for the Poisson and Stokes problems

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods, and relies on reformulating the original boundary value problem over a surrogate (approximate) computational domain. Accuracy is maintained by properly shifting the location and values of the boundary conditions. This avoids integration over cut cells and the associated implementation issues. Recently, the Weighted SBM (WSBM) was proposed for the Navier-Stokes equations with free surfaces and the Stokes flow with moving boundaries. The attribute “weighted” in the name WSBM stems from the fact that its variational form is weighted with the elemental volume fraction of active fluid. The motivation for the development of the WSBM was the preservation of the volume of active fluid to a higher degree of accuracy, which in turn resulted in improved stability and robustness characteristics in moving-boundary, time-dependent simulations. In this article, we present the numerical analysis of the WSBM formulations for the Poisson and Stokes problems. We give mathematical conditions under which the bilinear forms defining the discrete variational formulations are uniformly coercive (Poisson problem) or inf-sup stable (Stokes problem). By these results, stability and optimal convergence is proven in the natural norm; L2-error estimates can also be derived.

Approximate domain boundaries↗
DOE OSTI2021

Uranium hydride corrosion. II Modeling spatially random growth using reaction and nucleation kinetics

A uranium hydride (UH 3 ) corrosion kinetics model has been developed, based on dimensionally restricted nucleation kinetics, which combines heretofore separate reaction and nucleation kinetics models into a single predictive construct. This theoretical framework accounts for the superficial and random pitting behaviour of UH 3 corrosion. The model can be used to generate the net corrosion surface area and volume fractions for a uranium specimen, generally accounting for its geometric form factor through the surface-to-volume ratio. Furthermore, the combined corrosion kinetics model enables the generation of quantities which may be directly compared to those obtained by non-destructive, non-contact evaluation methods.

36 MATERIALS SCIENCE↗