Search NASA⌕ Search

SEARCH · Search NASA

Results for “curved surfaces”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Band structure and Fermi surface nesting in LaSb 2

Here we use high-resolution angle-resolved photoemission spectroscopy (ARPES) and density functional theory (DFT) to investigate the electronic structure of the charge density wave (CDW) system LaSb 2 . This compound is among an interesting group of materials that manifests both a CDW transition and lower temperature superconductivity. We find the DFT calculations to be in good agreement with our ARPES data. The Fermi surface of LaSb 2 consists of two small hole pockets close to Γ and four larger pockets near the Brillouin zone boundary. The key features of the Fermi surface do not vary significantly with temperature. A saddle point is present at -0.19 eV below the Fermi level at Γ. Saddle points in the band structure have more pronounced effects on materials properties when they are located closer to the Fermi level, making doped LaSb 2 compounds a potential interesting subject of future research. Multiple peaks are present in the generalized, electronic susceptibility calculations, indicating the presence of three possible nesting vectors. We were not able to detect any signatures of the CDW transition at 355 K down to the lowest temperature of 7 K, pointing to the subtle nature of this transition. This is unusual, given that such a high transition temperature is expected to be associated with the presence of a large CDW gap. This is confirmed through investigation of the Fermi surface and through analysis of momentum distribution curves. It is possible that changes are subtle and occur below current sensitivity of our measurements.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

A curvature-based mechanism for the spontaneous emission of electromagnetic radiation at curved plasma–vacuum interfaces

In this work, we derive from first principles the dispersion relation for electromagnetic radiation that is spontaneously emitted at a 2D curved plasma–vacuum interface in the absence of both background and external electric and magnetic fields. This plasma geometry is motivated by plasmas generated through hypervelocity impact projectiles—an area that is of importance to satellite safety in the aerospace community, and may be of importance for the developing study of projectile impact-driven inertial confinement fusion. We find spontaneous radiation occurs when thermal surface waves propagate along the edge of a circular plasma, with a frequency lower than the plasma frequency and spatial decay proportional to 1/r. This is distinct from more well-known mechanisms of radiation, in which the radiation frequency is greater than or equal to the plasma frequency. This theory is validated against 2D, electromagnetic particle-in-cell simulations and shows reasonable agreement. Finally, we discuss the effects of plasma geometry and attribute the radiation generation mechanism to the acceleration of charged surface waves as they traverse azimuthally along a curved path.

Lau, Raymond (ORCID:0000000299318101)↗
DOE OSTI2021

A highly efficient atomically thin curved PdIr bimetallene electrocatalyst

Abstract The multi-metallene with an ultrahigh surface area has great potential in precise tuning of surface heterogeneous d-electronic correlation by surface strain effect for the distinctive surface electronic structure, which is a brand new class of promising 2D electrocatalyst for sustainable energy device application. However, achieving such an atomically thin multi-metallene still presents a great challenge. Herein, we present a new synthetic method for an atomic-level palladium-iridium (PdIr) bimetallene with an average thickness of only ∼1.0 nm for achieving superior catalysis in the hydrogen evolution reaction (HER) and the formic acid oxidation reaction (FAOR). The curved PdIr bimetallene presents a top-ranked high electrochemical active area of 127.5 ± 10.8 m2 gPd+Ir−1 in the reported noble alloy materials, and exhibits a very low overpotential, ultrahigh activity and improved stability for HER and FAOR. DFT calculation reveals that the PdIr bimetallene herein has a unique lattice tangential strain, which can induce surface distortion while concurrently creating a variety of concave-convex featured micro-active regions formed by variously coordinated Pd sites agglomeration. Such a strong strain effect correlates the abnormal on-site active 4d10-t2g-orbital Coulomb correlation potential and directly elevates orbital-electronegativity exposure within these active regions, resulting in a preeminent barrier-free energetic path for significant enhancement of FAOR and HER catalytic performance.

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

High Carrier Lifetimes in Epitaxial Germanium-Tin/Al(In)As Heterostructures with Variable Tin Compositions

Group IV-based germanium-tin (Ge1-ySny) compositional materials have recently shown great promise for infrared detection, light emission and ultra-low power transistors. High carrier lifetimes are desirable for enhancing the detection limit and efficiency of photodetectors, low threshold current density in lasers, and low tunneling barrier height by lowering defects and dislocations at the heterointerface of a source and a channel. Here, carrier lifetimes in epitaxial germanium (Ge) and variable tin (Sn) compositional Ge1-ySny materials were experimentally determined on GaAs substrates using the contactless microwave photoconductive decay (u-PCD) technique at an excitation wavelength of 1500 nm. Sharp (2 x 2) reflection high energy electron diffraction patterns and low surface roughness were observed from the surface of the Ge0.97Sn0.03 epilayer. X-ray rocking curves from Ge0.97Sn0.03 and Ge0.94Sn0.06 layers demonstrated the pseudomorphic and lattice-matched growth on AlAs and In0.12Al0.88As buffers, respectively, further substantiated by reciprocal space maps and abrupt heterointerfaces evident from the presence of Pendellosung oscillations. High effective carrier lifetimes of 150 ns to 450 ns were measured for Ge1-ySny epilayers as a function of Sn composition, surface roughness, growth temperature, and layer thickness. The observed increase in the carrier lifetime with an increasing Ge layer thickness and a reducing surface roughness, by incorporating Sn, were explained. The enhancement of the carrier lifetime with an increasing Sn concentration was achieved by controlling the defects with lattice-matched Ge0.94Sn0.06/In0.12Al0.88As heterointerfaces or the pseudomorphic growth of Ge0.94Sn0.06 on GaAs. Therefore, our monolithic integration of variable Sn alloy compositional Ge1-ySny materials with high carrier lifetimes opens avenues to realize electronic and optoelectronic devices.

carrier lifetimes↗
DOE OSTI2020

Minimal-Area Metrics on the Swiss Cross and Punctured Torus

The closed string theory minimal-area problem asks for the conformal metric of least area on a Riemann surface with the condition that all non-contractible closed curves have length greater than or equal to an arbitrary constant that can be set to one. Through every point in such a metric there is a geodesic that saturates the length condition, and the saturating geodesics in a given homotopy class form a band. The extremal metric is unknown when bands of geodesics cross, as it happens for surfaces of non-zero genus. Here, we use recently proposed convex programs to numerically find the minimal-area metric on the square torus with a square boundary, for various sizes of the boundary. For large enough boundary the problem is equivalent to the “Swiss cross” challenge posed by Strebel. We find that the metric is positively curved in the two-band region and flat in the single-band regions. For small boundary the metric develops a third band of geodesics wrapping around it, and has both regions of positive and of negative curvature. This surface can be completed to provide the minimal-area metric on a once-punctured torus, representing a closed-string tadpole diagram.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Multicolor Variability of Young Stars in the Lagoon Nebula: Driving Causes and Intrinsic Timescales

Space observatories have provided unprecedented depictions of the many variability behaviors typical of low-mass, young stars. However, those studies have so far largely omitted more massive objects (~2 M {sub ⊙} to 4–5 M {sub ⊙}) and were limited by the absence of simultaneous, multiwavelength information. We present a new study of young star variability in the ~1–2 Myr old, massive Lagoon Nebula region. Our sample encompasses 278 young, late B to K-type stars, monitored with Kepler/K2. Auxiliary u, g, r, i, Hα time-series photometry, simultaneous with K2, was acquired at the Paranal Observatory. We employed this comprehensive data set and archival infrared photometry to determine individual stellar parameters, assess the presence of circumstellar disks, and tie the variability behaviors to inner disk dynamics. We found significant mass-dependent trends in variability properties, with B/A stars displaying substantially reduced levels of variability compared to G/K stars for any light-curve morphology. These properties suggest different magnetic field structures at the surface of early-type and later-type stars. We also detected a dearth of some disk-driven variability behaviors, particularly dippers, among stars earlier than G. This indicates that their higher surface temperatures and more chaotic magnetic fields prevent the formation and survival of inner disk dust structures corotating with the star. Finally, we examined the characteristic variability timescales within each light curve and determined that the day-to-week timescales are predominant over the K2 time series. These reflect distinct processes and locations in the inner disk environment, from intense accretion triggered by instabilities in the innermost disk regions to variable accretion efficiency in the outer magnetosphere.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Binder-Coated Carbon Cloth Electrodes for All-Vanadium Redox Flow Batteries

Vanadium redox flow batteries (VRFBs) are a promising solution for integrating intermittent renewable energy sources into the existing power grid. However, enhancing the electrochemical performance of VRFBs is critical for their widespread adoption in grid-scale energy storage. This study investigates the impact of adding a porous binder to a carbon-cloth electrode, with a focus on optimizing thermal activation conditions. The electrochemical performance of the binder-coated electrodes compared to uncoated electrodes is evaluated through electrochemical impedance spectroscopy, polarization curve measurements, and charge-discharge cycling. The surface morphology and structural integrity of the binder-coated electrodes at each activation stage are examined using various material characterization techniques to assess the effects of thermal activation. The results are benchmarked against the experiments using non-coated electrodes to determine the performance improvements offered by the binder coating. Notably, the study reveals that binder-coated electrodes exhibit significantly lower resistance and improved efficiency compared to their uncoated counterparts, with optimal activation conditions enhancing performance metrics crucial for VRFB applications. These findings provide valuable insights for further optimizing electrode design and activation strategies, advancing the development of more efficient VRFB systems for large-scale energy storage.

Caiado, Ashley A.↗
DOE OSTI2020

MOCVD Growth and Characterization of ZnSnN 2

Since the first few reports of its synthesis in 2011 and 20121-3, interest in ZnSnN 2 for diverse applications, including photovoltaics1,3,4, visible light photocatalysis5, and optoelectronics, has been growing. Novel nitride-based LED device designs promise high gains in efficiency at emission wavelengths into the green, amber and even red.6-8 Their realization requires the precise insertion of an ultrathin layer of ZnGeN 2 or ZnSnN 2 into the active region of an InGaAs-GaAs LED structure. While thin films of ZnSnN 2 have been synthesized by RF sputtering1,9 and by MBE 2,10,11, and thin films of ZnGeN 2 , ZnSiN 2 and their alloys by MOCVD, the synthesis of ZnSnN 2 thin films by MOCVD has not yet been reported. It is a challenging goal, given the narrow temperature window between efficient activation of ammonia and the thermal decomposition of ZnSnN 2 at near-atmospheric pressures. We report here the results of studies of ZnSnN 2 films grown on r- and c-plane sapphire substrates, GaN templates on sapphire, and In 0.15 Ga 0.85 N templates on GaN-on-sapphire, at temperatures from 500 °C to 560 °C, using tetramethyltin (TMSn), diethylzinc (DEZn), germane and ammonia. Chamber pressures were varied between 200 and 500 torr. In addition, the ratios of the cation precursors and the total cation-to-anion precursor ratios were varied. Single crystal growth, as evidenced by 2θ-ω x-ray diffraction scans, was achieved on all of these substrates. Wurtzitic (0002) ZnSnN 2 diffraction were observed at 2θ ~ 32.90 for the samples grown on c-sapphire, GaN (0001) and InGaN (0001) templates, and wurtzitic (11-20) peak at 2θ ~ 54.80 for the films grown on r-sapphire. Higher pressures and higher temperatures promoted better crystallinity and more continuous films, as revealed by scanning electron microscopy (SEM) images and 2θ-ω x-ray diffraction peak widths and intensities. Crystal quality was observed to be highly dependent on growth temperature. For example, in one series a 10 °C change in growth temperature resulting in a factor of 15 increase in diffraction peak intensity. SEM images show that films grown on the GaN templates had better uniformity in nucleation than those grown on sapphire. The presence of satellite peaks around the Bragg peak in the XRD 2θ-ω scan for a film grown on a GaN template, and the absence of satellite peaks for the films grown simultaneously on c- and r-sapphire, also indicate superior quality for the film grown on GaN. HAADF-STEM images of films grown on these templates showed smooth interfaces at the film-substrate boundary and show growth rates as high as 167 nm/h, considerably higher than reported previously for MBE growth.9,10 Electron diffraction patterns show single crystallinity. Rocking curves show FWHM as low as 0.100. RMS surfaces roughnesses, measured by atomic force microscopy, are of the order of 3 nm for 1µm x 1µm areas. Current work is focusing on optimizing growth on InGaN templates. These templates are promising for their relatively lower lattice mismatch with ZnSnN 2 . Their use is a next step toward producing LED device structures.

Jayatunga, Benthara Hewage Dinushi↗
DOE OSTI2020

Effects of surface runoff and infiltration partition methods on hydrological modeling: A comparison of four schemes in two watersheds in the Northeastern US

The ratio according to which the precipitation/snowmelt is divided into surface and subsurface water is of primary importance to watershed scale hydrological modeling. The present study compared four different surface runoff and infiltration partition methods in the Richards-equation-based SWAT model (RSWAT), i.e., daily curve number method (DC-RSWAT), hourly curve number method (HC-RSWAT), Green-Ampt method (GA-RSWAT), and effective infiltration capacity method (EIC-RSWAT). These four versions of RSWAT and a daily-curve-number-based SWAT (DC-SWAT) were used to simulate daily flow rate (from 2001 to 2015) at the outlets of the Tuckahoe Creek Watershed (TCW) and Greensboro Watershed (GW) of the Chesapeake Bay in Maryland. Parameter sensitivity analysis, calibration, and prediction uncertainty analysis were conducted. The results indicate that surface runoff and infiltration partition methods have pronounced impacts on model performance. Comparative analysis showed that DC-RSWAT outperformed other versions of RSWAT for both calibration (from 2001 to2007) and validation (from 2008 to 2015) periods in the two test watersheds. In general, DC-RSWAT achieved more accurate deterministic flow simulation and reliable uncertainty estimation than other versions. As such, daily-curve-number method is recommended for surface runoff and infiltration partition in RSWAT for simulating stream flow.

hydrology, modeling, streamflow↗
DOE OSTI2020

X-Ray Afterglows from the Gamma-Ray Burst “Large-angle” Emission

In this work, we derive basic analytical results for the timing and decay of the gamma-ray burst (GRB) counterpart and delayed afterglow light curves for a brief emission episode from a relativistic surface endowed with angular structure, consisting of a uniform core of size ${\theta }_{c}$ (Lorentz factor ${{\rm{\Gamma }}}_{c}$ and surface emissivity ${i}_{\nu ^{\prime} }^{{\prime} }$ are angle independent) and an axially symmetric power-law envelope (${\rm{\Gamma }}\sim {\theta }^{-g}$). In this "large-angle emission" model, radiation produced during the prompt emission phase (GRB) at angles $\theta \gt {\theta }_{c}$ arrives at the observer well after the burst (delayed emission). The dynamical time range of the very fast decaying GRB "tail" and of the flat afterglow "plateau" and the morphology of the GRB counterpart/afterglow are all determined by two parameters: the core's parameter ${{\rm{\Gamma }}}_{c}{\theta }_{c}$ and the envelope's Lorentz factor index g, leading to three types of light curves that display three post-GRB phases (type 1: tail, plateau/slow decay, post-plateau/normal decay), two post-GRB phases (type 2: tail and fast decay), or just one (type 3: normal decay). We show how X-ray light-curve features can be used to determine core and envelope dynamical and spectral parameters. Testing of the large-angle emission model is done using the Swift/XRT X-ray emission of two afterglows of type 1 (GRB 060607A, GRB 061121), one of type 2 (GRB 061110A), and one of type 3 (GRB 061007). We find that the X-ray afterglows with plateaus require an envelope Lorentz factor ${\rm{\Gamma }}\sim {\theta }^{-2}$ and a comoving-frame emissivity ${i}_{\nu ^{\prime} }^{{\prime} }\sim {\theta }^{2}$; thus, for a typical afterglow spectrum ${F}_{\nu }\sim {\nu }^{-1}$, the lab-frame energy release is uniform over the emitting surface.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

The SN Ia runaway LP 398-9: detection of circumstellar material and surface rotation

ABSTRACT A promising progenitor scenario for Type Ia supernovae (SNeIa) is the thermonuclear detonation of a white dwarf in a close binary system with another white dwarf. After the primary star explodes, the surviving donor can be spontaneously released as a hypervelocity runaway. One such runaway donor candidate is LP 398-9, whose orbital trajectory traces back ≈105 yr to a known supernova remnant. Here, we report the discovery of carbon-rich circumstellar material around LP 398-9, revealed by a strong infrared excess and analysed with follow-up spectroscopy. The circumstellar material is most plausibly composed of inflated layers from the star itself, mechanically and radioactively heated by the past companion’s supernova. We also detect a 15.4 h periodic signal in the UV and optical light curves of LP 398-9, which we interpret as surface rotation. The rotation rate is consistent with theoretical predictions from this supernova mechanism, and the brightness variations could originate from surface inhomogeneity deposited by the supernova itself. Our observations strengthen the case for this double-degenerate SNIa progenitor channel, and motivate the search for more runaway SNIa donors.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Estimation of the Thermal-Fluids and Thermal-Structural Performance of Helium-Cooled Modular Finger-Type Divertors

Over the past decade, our group has investigated the thermal-fluid performance of the helium-cooled modular divertor with multiple jets (HEMJ) and a simplified “flat” design of the HEMJ for long-pulse magnetic fusion energy (MFE) reactors. Experimental studies were performed in a helium (He) loop at the prototypical pressure of 10 MPa, nearly prototypical He temperatures and incident heat fluxes using test sections made from stainless steel and tungsten alloys. Correlations for average Nusselt numbers and pressure loss coefficients were developed from the data and are used to validate computational fluid dynamics (CFD) models. Here, this work presents updated thermal-fluids performance curves based on these correlations that estimate the maximum heat flux that can be accommodated by the plasma-facing surface and coolant pumping power requirements at prototypical operating conditions. Thermal-structural performance curves developed from ITER structural design criteria are introduced, which include protection against ductile and non-ductile failure, ratcheting fatigue, and creep fatigue. The performance design curves for these finger-type divertors demonstrate that the “flat” design, with a significantly less complicated geometry than the HEMJ, has thermal-fluid and thermal-structural performance comparable to the original HEMJ concept, and may be superior in terms of non-ductile failure criteria.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Density Measurements of Molten LiF–BeF 2 and LiF–BeF 2 –LaF 3 Salt Mixtures by Neutron Radiography

The densities of eutectic (LiF) 2 –BeF 2 and mixtures of this salt (FLiBe) with LaF 3 were measured by dilatometry and by neutron attenuation from 673 K to 1,073 K. Because LaF 3 has a limited solubility in FLiBe, it was necessary to determine the amount of LaF 3 in solution before the density could be determined. The FLiBe density determination was favorably benchmarked against the literature data. A simple comparison was not available for the LaF 3 –FLiBe mixtures, so extrapolation of published data was necessary based on analysis using the Molten Salt Thermal Properties Database-Thermochemistry, or MSTDB-TC, developed by the US Department of Energy. Solubilities for LaF 3 in FLiBe ranged from 1 to 4 mol % over 673 to 1,073 K. The salt system was heated and cooled over 24 h to evaluate potential changes in composition and hysteresis during the measurement. Changes in the meniscus were observed, and these were included in the correction for density determinations. Salt surface tension may have led to supersaturation of LaF 3 in the salt because the solubility curve was nonlinear with respect to the inverse temperature, as would be expected for an ideal system. Surface tension measurements are currently underway to test this hypothesis.

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

Correlation of coercivity, microstructure, and surface defects in HDDR-processed Nd-Fe-B powders for bonded magnet applications

Hydrogenation-Disproportionation-Desorption-Recombination (HDDR) is an effective method for producing high coercivity, anisotropic powder from either fresh Nd-Fe-B alloy or recycled Nd-Fe-B magnets for use in bonded magnets. We investigated the impact of surface defects on the coercivity of HDDR-processed Nd-Fe-B using scanning electron microscopy (SEM), magnetic measurements, and micromagnetic simulations. We observed that coercivity decreases as particle size reduces, with SEM revealing surface defects and the detachment of Nd 2 Fe 14 B grains and Nd-rich phases from the particle surface. Micromagnetic simulations indicate that demagnetization initiates at the particle surface, where these defects are most concentrated, leading to reduced coercivity. The reduction in squareness of demagnetization curve, knee point field and coercivity for smaller HDDR particles is attributed to an increased specific surface area, which exhibits reduced nucleation field and weak domain wall pinning field during magnetization reversal. By addressing the role of surface defects in coercivity degradation, this study provides insights for improving both new powder production and recycling strategies, ultimately leading to enhanced performance of bonded magnets and contributing to more sustainable practices in the rare earth supply chain. One potential strategy to enhance the performance of HDDR Nd-Fe-B materials involves reducing the fraction of fine particles ($<$35 µm) and promoting the formation of grain boundary phases on particle surfaces.

Liu, Xubo B. [Ames Laboratory (AMES), Ames, IA (Un↗
DOE OSTI2024

The Dark Energy Survey Supernova Program: Light Curves and 5 Yr Data Release

We present griz photometric light curves for the full 5 yr of the Dark Energy Survey Supernova (DES-SN) program, obtained with both forced point-spread function photometry on difference images (DiffImg) performed during survey operations, and scene modelling photometry (SMP) on search images processed after the survey. This release contains 31,636 DiffImg and 19,706 high-quality SMP light curves, the latter of which contain 1635 photometrically classified SNe that pass cosmology quality cuts. This sample spans the largest redshift (z) range ever covered by a single SN survey (0.1 < z < 1.13) and is the largest single sample from a single instrument of SNe ever used for cosmological constraints. We describe in detail the improvements made to obtain the final DES-SN photometry and provide a comparison to what was used in the 3 yr DES-SN spectroscopically confirmed Type Ia SN sample. We also include a comparative analysis of the performance of the SMP photometry with respect to the real-time DiffImg forced photometry and find that SMP photometry is more precise, more accurate, and less sensitive to the host-galaxy surface brightness anomaly. The public release of the light curves and ancillary data can be found at github.com/des-science/DES-SN5YR and doi:10.5281/zenodo.12720777.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Robust Containment Queries over Collections of Rational Parametric Curves via Generalized Winding Numbers

Point containment queries for regions bound by watertight geometric surfaces, i.e., closed and without self-intersections, can be evaluated straightforwardly with a number of well-studied algorithms. When this assumption on domain geometry is not met, such methods are either unusable, or prone to misclassifications that can lead to cascading errors in downstream applications. More robust point classification schemes based on generalized winding numbers have been proposed, as they are indifferent to these imperfections. However, existing algorithms are limited to point clouds and collections of linear elements. We extend this methodology to encompass more general curved shapes with an algorithm that evaluates the winding number scalar field over unstructured collections of rational parametric curves. In particular, we evaluate the winding number for each curve independently, making the derived containment query robust to how the curves are arranged. We ensure geometric fidelity in our queries by treating each curve as equivalent to an adaptively constructed polyline that provably has the same generalized winding number at the point of interest. Our algorithm is numerically stable for points that are arbitrarily close to the model, and explicitly treats points that are coincident with curves. We demonstrate the improvements in computational performance granted by this method over conventional techniques as well as the robustness induced by its application.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Prediction of alkaline earth metal ion adsorption on goethite for various background electrolytes with the CD-MUSIC model

As water scarcity drives the use of more saline water sources, contaminant fate and transport models must capture the impact of high concentrations of alkaline earth metal ions (AEMs) and background electrolytes in these more complex waters. By utilizing macroscopic adsorption data from various electrolyte systems, a Charge Distribution – Multisite Complexation (CD-MUSIC) model, capable of incorporating electrolyte adsorption, was able to accurately simulate the adsorption behavior of alkaline earth metal ions onto goethite. The modeling effort was guided by previous spectroscopic and surface complexation modeling of alkaline earth metal adsorption and built on previous CD-MUSIC modeling that accounted for changes in crystal face contributions to the surface site density as a function of specific surface area. Further, the model was constrained to consider only two dominant surface complex species for each metal ion adsorption reaction. These two species were selected from 44 possible species through objective curve fitting of single-solute macroscopic adsorption data. While most of the alkaline earth metal surface complexes formed outer-sphere complexes at the goethite surface, an inner-sphere species was utilized for Mg 2+ . With the surface complex species and equilibrium constants obtained from this study, the calibrated model successfully predicted alkaline earth metal ion adsorption over a wide range of solution and surface conditions; the model predictions encompassed a wide range of pH (5–11), solute/solid ratio (1.37 × 10 -5 – 8.33 × 10 -4 mol -solute /g -solid ), ionic strengths (0.01 M – 0.7 M), and background electrolytes (Na + , Cs + , Rb + , Cl - , and NO 3 - ) using the same crystal face contribution methodology for site density, capacitance values, and surface acidity constants adopted for proton and cadmium adsorption in previous work (Han and Katz, 2019). Model simulations for a range of background water chemistries demonstrated the potential for Mg 2+ to reduce Cd 2+ adsorption to goethite in model seawater and oil- and gas-produced waters.

42 ENGINEERING↗
DOE OSTI2024

An integrated approach to derive relative permeability from capillary pressure

Surface tension affects all aspects of fluid flow in porous media. Through measurements of surface tension interaction under multiphase conditions, a relative permeability curve can be determined. Relative permeability is a numerical description of the interaction between two or more fluids and the porous media. It is a critical parameter for various tools that characterize subsurface multiphase flow systems, such as numerical simulation for carbon sequestration, oil and gas development, and groundwater contamination remediation. Therefore, it is critical to get a good statistical distribution of relative permeability in the porous media under study. Empirical formula for determining relative permeability from capillary pressure are already well established but do not provide the needed flexibility that is required to match laboratory-derived relative permeability curves. By expanding the existing methods for calculating relative permeability from capillary pressure data, it is possible to create both two and three-phase relative permeability curves. Mercury intrusion capillary pressure (MICP) data from the Morrow 'B' Sandstone coupled with interfacial tension and contact angle measurements were used to create a suite of relative permeability curves. Furthermore, these curves were then calibrated to a small sample of existing laboratory curves to elucidate common fitting parameters for the formation that were then used to create relative permeability curves from MICP data that does not have an associated laboratory-measured relative permeability curve.

58 GEOSCIENCES↗