Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SATURATION”

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 325 records · Page 18

Rocky Worlds Limited to ∼1.8 Earth Radii by Atmospheric Escape during a Star’s Extreme UV Saturation

Recent observations and analysis of low-mass (<10 M {sub ⊕}) exoplanets have found that rocky planets only have radii up to 1.5–2 R {sub ⊕}. Two general hypotheses exist for the cause of the dichotomy between rocky and gas-enveloped planets (or possible water worlds): either low-mass planets do not necessarily form thick atmospheres of a few wt.%, or the thick atmospheres on these planets easily escape, driven by X-ray and extreme ultraviolet (XUV) emissions from young parent stars. Here, we show that a cutoff between rocky and gas-enveloped planets due to hydrodynamic escape is most likely to occur at a mean radius of 1.76 ± 0.38 (2 σ ) R {sub ⊕} around Sun-like stars. We examine the limit in rocky planet radii predicted by hydrodynamic escape across a wide range of possible model inputs, using 10,000 parameter combinations drawn randomly from plausible parameter ranges. We find a cutoff between rocky and gas-enveloped planets that agrees with the observed cutoff. The large cross-section available for XUV absorption in the extremely distended primitive atmospheres of low-mass planets results in complete loss of atmospheres during the ∼100 Myr phase of stellar XUV saturation. In contrast, more-massive planets have less-distended atmospheres and less escape, and so retain thick atmospheres through XUV saturation—and then indefinitely as the XUV and escape fluxes drop over time. The agreement between our model and exoplanet data leads us to conclude that hydrodynamic escape plausibly explains the observed upper limit on rocky planet size and few planets (a “valley”, or “radius gap”) in the 1.5–2 R {sub ⊕} range.

79 ASTRONOMY AND ASTROPHYSICS↗

Capillary pressure-saturation relations in quartz and carbonate sands: Limitations for correlating capillary and wettability influences on air, oil, and supercritical CO 2 trapping

Capillary pressure (P c )-saturation (S w ) relations are essential for predicting equilibrium and flow of immiscible fluid pairs in soils and deeper geologic formations. In systems that are difficult to measure, behavior is often estimated based on capillary scaling of easily measured P c –S w relations (e.g., air-water, and oil-water), yet the reliability of such approximations needs to be examined. In this study, 17 sets of brine drainage and imbibition curves were measured with air-brine, decane-brine, and supercritical (sc) CO 2 -brine in homogeneous quartz and carbonate sands, using porous plate systems under ambient (0.1 MPa, 23°C) and reservoir (12.0 MPa, 45°C) conditions. Comparisons between these measurements showed significant differences in residual nonwetting phase saturation, S nw,r . Through applying capillary scaling, changes in interfacial properties were indicated, particularly wettability. With respect to the residual trapping of the nonwetting phases, S nwr, CO2 > S nwr, decane > S nwr, air . Decane-brine and scCO 2 -brine Pc–Sw curves deviated significantly from predictions assuming hydrophilic interactions. Moreover, neither the scaled capillary behavior nor S nw,r for scCO 2 -brine were well represented by decane-brine, apparently because of differences in wettability and viscosities, indicating limitations for using decane (and other organic liquids) as a surrogate fluid in studies intended to apply to geological carbon sequestration. Thus, challenges remain in applying scaling for predicting capillary trapping and multiphase displacement processes across such diverse fields as vadose zone hydrology, enhanced oil recovery, and geologic carbon sequestration.

58 GEOSCIENCES↗

Water Saturation Relations and Their Diffusion-Limited Equilibration in Gas Shale: Implications for Gas Flow in Unconventional Reservoirs

Large volumes of water are used for hydraulic fracturing of low permeability shale reservoirs to stimulate gas production, with most of the water remaining unrecovered and distributed in a poorly understood manner within stimulated regions. Because water partitioning into shale pores controls gas release, we measured the water saturation dependence on relative humidity (rh) and capillary pressure (P c ) for imbibition (adsorption) as well as drainage (desorption) on samples of Woodford Shale. Experiments and modeling of water vapor adsorption into shale laminae at rh = 0.31 demonstrated that long times are needed to characterize equilibrium in larger (5 mm thick) pieces of shales, and yielded effective diffusion coefficients from 9 × 10 -9 to 3 × 10 -8 m 2 s -1 , similar in magnitude to the literature values for typical low porosity and low permeability rocks. Most of the experiments, conducted at 50°C on crushed shale grains in order to facilitate rapid equilibration, showed significant saturation hysteresis, and that very large P c (~1 MPa) are required to drain the shales. These results quantify the severity of the water blocking problem, and suggest that gas production from unconventional reservoirs is largely associated with stimulated regions that have had little or no exposure to injected water. Finally, gravity drainage of water from fractures residing above horizontal wells reconciles gas production in the presence of largely unrecovered injected water, and is discussed in the broader context of unsaturated flow in fractures.

58 GEOSCIENCES↗

An analytical approach to ascertain saturation‐excess versus infiltration‐excess overland flow in urban and reference landscapes

Abstract Uncontrolled overland flow drives flooding, erosion, and contaminant transport, with the severity of these outcomes often amplified in urban areas. In pervious media such as urban soils, overland flow is initiated via either infiltration‐excess (where precipitation rate exceeds infiltration capacity) or saturation‐excess (when precipitation volume exceeds soil profile storage) mechanisms. These processes call for different management strategies, making it important for municipalities to discern between them. In this study, we derived a generalized one‐dimensional model that distinguishes between infiltration‐excess overland flow (IEOF) and saturation‐excess overland flow (SEOF) using Green–Ampt infiltration concepts. Next, we applied this model to estimate overland flow generation from pervious areas in 11 U.S. cities. We used rainfall forcing that represented low‐ and high‐intensity events and compared responses among measured urban versus predevelopment reference soil hydraulic properties. The derivation showed that the propensity for IEOF versus SEOF is related to the equivalence between two nondimensional ratios: (a) precipitation rate to depth‐weighted hydraulic conductivity and (b) depth of soil profile restrictive layer to soil capillary potential. Across all cities, reference soil profiles were associated with greater IEOF for the high‐intensity set of storms, and urbanized soil profiles tended towards production of SEOF during the lower intensity set of storms. Urban soils produced more cumulative overland flow as a fraction of cumulative precipitation than did reference soils, particularly under conditions associated with SEOF. These results will assist cities in identifying the type and extent of interventions needed to manage storm water produced from pervious areas.

Stewart, Ryan D.↗

Estimating uncertainty in ambient and saturation nutrient uptake metrics from nutrient pulse releases in stream ecosystems

Nutrient spiraling is an important ecosystem process characterizing nutrient transport and uptake in streams. Various nutrient addition methods are used to estimate uptake metrics; however, uncertainty in the metrics is not often evaluated. A method was developed to quantify uncertainty in ambient and saturation nutrient uptake metrics estimated from saturating pulse nutrient additions (Tracer Additions for Spiraling Curve Characterization; TASCC). Using a Monte Carlo (MC) approach, the 95% confidence interval (CI) was estimated for ambient uptake lengths (S w-amb ) and maximum areal uptake rates (U max ) based on 100,000 datasets generated from each of four nitrogen and five phosphorous TASCC experiments conducted seasonally in a forest stream in eastern Tennessee, U.S.A. Uncertainty estimates from the MC approach were compared to the CIs estimated from ordinary least squares (OLS) and non-linear least squares (NLS) models used to calculate S w-amb and U max , respectively, from the TASCC method. The CIs for Sw-amb and Umax were large, but were not consistently larger using the MC method. Despite the large CIs, significant differences (based on nonoverlapping CIs) in nutrient metrics among seasons were found with more significant differences using the OLS/NLS vs. the MC method. Lastly, we suggest that the MC approach is a robust way to estimate uncertainty, as the calculation of S w-amb and U max violates assumptions of OLS/NLS while the MC approach is free of these assumptions. The MC approach can be applied to other ecosystem metrics that are calculated from multiple parameters, providing a more robust estimate of these metrics and their associated uncertainties.

54 ENVIRONMENTAL SCIENCES↗

First saturation correction in high energy proton-nucleus collisions. Part II. Single inclusive semi-hard gluon production

Exploiting recently obtained analytic solutions of classical Yang-Mills equations for higher order perturbations in the field of the dilute object (proton), we derive the complete first saturation correction to the single inclusive semi-hard gluon production in high energy proton-nucleus collisions by applying the Lehmann-Symanzik-Zimmermann reduction formula. We thus finalize the program started by Balitsky (see ref. [1]) and independently by Chirilli, Kovchegov and Wertepny (see ref. [2]) albeit using a very different approach to carry out our calculations. We extracted the functional dependence of gluon spectrum on the color charge densities of the colliding objects; thus our results can be used to evaluate complete first saturation correction to the double/multiple inclusive gluon productions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Back-to-back inclusive dijets in DIS at small x : Sudakov suppression and gluon saturation at NLO

Back-to-back dijet cross-sections in deeply inelastic scattering (DIS) at small xBj are suppressed by many-body multiple scattering and screening effects arising from gluon saturation at high parton densities. They are similarly sensitive in these kinematics to large Sudakov logarithms from soft gluon radiation. Uncovering novel physics in this DIS channel therefore requires understanding the interplay of the two phenomena. In this work, we compute the small $x_\text{Bj}$ inclusive dijet DIS cross-section in back-to-back kinematics at next-to-leading order (NLO) in the Color Glass Condensate effective field theory (CGC EFT). Our result includes, for the first time, all real and virtual NLO contributions to the impact factor. These include all Sudakov double and single logarithm contributions, as well as all other finite $\mathcal{O}(α_s)$ terms that contribute at this order. We demonstrate explicitly that resummations of small x and Sudakov logarithms can be performed simultaneously in the CGC EFT. This requires that the JIMWLK kernel for small x evolution of the Weizsäcker-Williams (WW) gluon distribution satisfies a kinematic constraint imposed by lifetime ordering of successive gluon emissions; the corresponding modifications to the kernel, corresponding to resummations of large double transverse logarithms, are precisely of the type required to stabilize JIMWLK evolution beyond leading logarithmicaccuracy. We compute the azimuthal harmonics of the NLO back-to-back distributions and show their sensitivity to both the unpolarized and linearly polarized WW gluon distributions. Finally, we discuss how TMD factorization is broken by an emergent saturation scale at small $x$.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Divisible Rigid Groups. II. Stability, Saturation, and Elementary Submodels

A group G is said to be rigid if it contains a normal series G = G{sub 1} > G{sub 2} > . . . > G{sub m} > G{sub m+1} = 1, whose quotients G{sub i}/G{sub i+1} are Abelian and, treated as right ℤ[G/G{sub i}]- modules, are torsion-free. A rigid group G is divisible if elements of the quotient G{sub i}/G{sub i+1} are divisible by nonzero elements of the ring ℤ[G/G{sub i}]. Every rigid group is embedded in a divisible one. Previously, it was stated that the theory 𝔗{sub m} of divisible m-rigid groups is complete. Here, it is proved that this theory is ω-stable. Furthermore, we describe saturated models, study elementary submodels of an arbitrary model, and find a representation for a countable saturated model in the form of a limit group in the Fraïssé system of all finitely generated m-rigid groups. Also, it is proved that the theory 𝔗{sub m} admits quantifier elimination down to a Boolean combination of ∀∃-formulas.

97 MATHEMATICS AND COMPUTING↗

An Experimental Method to Measure Gaseous Diffusivity in Tight and Partially Saturated Porous Media via Continuously Monitored Mass Spectrometry

Detection of radioxenon and radioargon produced by underground nuclear explosions is one of the primary methods by which the Comprehensive Nuclear-Test–Ban Treaty (CTBT) monitors for nuclear activities. However, transport of these noble gases to the surface via barometric pumping is a complex process relying on advective and diffusive processes in a fractured porous medium to bring detectable levels to the surface. To better understand this process, experimental measurements of noble gas and chemical surrogate diffusivity in relevant lithologies are necessary. However, measurement of noble gas diffusivity in tight or partially saturated porous media is challenging due to the transparent nature of noble gases, the lengthy diffusion times, and difficulty maintaining consistent water saturation. Here, the quasi-steady-state Ney–Armistead method is modified to accommodate continuous gas sampling via effusive flow to a mass spectrometer. An analytical solution accounting for the cumulative sampling losses and induced advective flow is then derived. Experimental results appear in good agreement with the proposed theory, suggesting the presence of retained groundwater reduces the effective diffusivity of the gas tracers by 10–1000 times. Finally, by using a mass spectrometer, the method described herein is applicable to a broad range of gas species and porous media.

47 OTHER INSTRUMENTATION↗

Observer-based H{sub ∞} control for nonlinear Markovian jump systems with time-delay and input saturation

This paper focus on the problem of observer-based H{sub ∞} control for uncertain nonlinear Markovian jump systems with time-delay and actuator saturation. A delay-dependent sufficient condition is proposed which guarantees that the uncertain nonlinear Markovian jump system with time-delay and input saturation is stochastically stable via Lyapunov theory and linear matrix inequality (LMI) approach. Then, with this condition, the estimation of stability region and the design method of observer-based H{sub ∞} controller are given by solving LMIs and convex optimization problems. Finally, numerical examples are exploited to illustrate the effectiveness of the proposed method.

97 MATHEMATICS AND COMPUTING↗

On the importance of phase saturation heterogeneity in the analysis of laboratory studies of hydrate dissociation

Methane hydrates (MHs) have been considered as the future source of energy because of its vast resource volume and high energy density. Energy recovery from MH-bearing sediments has attracted intensifying research activities. Fundamentally, heat transfer, fluid flow through porous media, and the kinetics of hydrate reaction are the three key processes controlling the behavior of MH dissociation and the associated fluid production. Earlier studies have suggested that heterogeneous spatial distribution of SH is inevitable in MH-bearing samples synthesized in laboratory. In this paper, we extend our study to analyze numerically the simulation results from the two realizations of the samples (homogeneous and heterogeneous) to identify differences in the fluid production and to determine if they are sufficiently different. Additionally, we conduct a sensitivity analysis and a statistical analysis on the key transport and kinetic rate parameters that could affect hydrate dissociation and fluid production in the context of a heterogeneous hydrate-bearing sample, in an effort to provide insights that could lead to improved designs for laboratory experiments and (possibly) field applications. Our results suggest that the approximation of an artificial hydrate-bearing core with heterogeneous phase saturations by an assumption of uniform phase saturation distributions results in practically similar fluid production profile except for the very early stage with maximum 20.0% deviation in the water production. From the sensitivity and statistical analysis, we determine that gas production depends strongly on the kinetic rate constant, K d 0 and the composite thermal conductivity of the hydrate-bearing sediments, λ θ ; while, water production is very sensitive to K d 0 and the absolute permeability of the sandy medium, k. Understanding the effect of phase heterogeneity and the relative importance of key parameters on the production behavior of hydrate-bearing sediments could provide basis for novel production technologies that lead to enhanced gas production and energy efficiency in the energy recovery process.

42 ENGINEERING↗

The performance of a surface-applied corrosion inhibitor for the carbon steel in saturated Ca(OH){sub 2} solutions

In the present work, the performance of an amino alcohol based surface applied inhibitor was studied by the electrochemical techniques in saturated Ca(OH){sub 2} solutions. The surface morphology of the carbon steel was observed by scanning electron microscope, and the energy diffraction spectrum was also tested. Results showed that the inhibitor used in this work demonstrated obvious inhibition efficiency on the carbon steel in saturated Ca(OH){sub 2} solutions. The inhibition mechanism of the inhibitor lies in the quick adsorption of the active component on carbon steel surface.

36 MATERIALS SCIENCE↗

Experimental and modeling study of the autoignition behavior of a saturated heterocycle: Pyrrolidine

Experiments are conducted in both rapid compression machine (RCM) and shock tube (ST) to better quantify autoignition behavior (e.g., ignition delay, heat release) and understand heteroatomic effects in heterocyclic compounds, which are important reference components for the combustion of biomass-derived liquid fuels. These tests focus on the nitrogen-containing, five-membered saturated ring, pyrrolidine, at diluted conditions covering pressures of 20 and 50 bar, temperatures of 720–1450 K and a range of stoichiometries (ϕ = 0.5–2). A chemical kinetic model is developed and coupled to an existing combustion kinetics framework describing key nitrogen containing intermediates (e.g. pyrrole, ammonia and NOx). H-abstraction reactions by OH, H, CH 3 and HO 2 , are determined using ab-initio transition state theory methods, while analogies to cyclopentane are adopted for many other reactions, such as ring-opening. Here, the autoignition measurements reveal the lack of negative temperature coefficient (NTC) behavior and low-temperature chemistry for pyrrolidine, as opposed to its saturated hydrocarbon analogue, cyclopentane. Interestingly, at the lowest temperatures (T < 750 K), the reactivity of cyclopentane is greater than pyrrolidine, while at higher temperatures, pyrrolidine becomes more reactive. Agreement between the experimental measurements and the model is good, and it is found that H-abstraction reactions by HO 2 and ensuing chemistry play key roles in controlling the reactivity of this cyclic amine. Most of the fluxes, i.e., >70 %, are predicted to move through 1- or 2-pyrroline (C 4 H 7 N) and then the cyclic C 4 H 6 N radical, at both lower and higher temperatures, to form either CH 2 CHCHCHNH via ring-opening or pyrrole via β-scission. It appears that the ring opens more easily at lower temperature whereas the C–H β-scission dominates at higher temperature and lower pressure, such that the reaction of the fuel radical intermediate carrying an unpaired electron on the nitrogen atom with HO 2 is the next most notable in promoting oxidation. When comparing pyrrolidine and cyclopentane, which exhibits distinct pathways in different temperature regimes, the pyrrolidine pathways and sensitivity analysis align more closely to the high temperature case of cyclopentane where the important role of HO 2 radicals is seen to provide chain branching through HO 2 reaction with the fuel, accompanied by H 2 O 2 formation and decomposition to OH. The formation of 5-membered diene rings and ring opening reactions are also found to be highly relevant. Of particular note, it is found that there is little influence of small molecule nitrogen-chemistry, e.g., NH 2 , HCN, NO/NO 2 on the reactivity of the pyrrolidine mixtures investigated here where no recirculated combustion gases are included.

Autoignition↗

Stochastic DC optimal power flow with reserve saturation

We propose an optimization framework for stochastic optimal power flow with uncertain loads and renewable generator capacity. Our model follows previous work in assuming that generator outputs respond to load imbalances according to an affine control policy, but introduces a model of saturation of generator reserves by assuming that when a generator’s target level hits its limit, it abandons the affine policy and produces at that limit. This is a particularly interesting feature in models where wind power plants, which have uncertain upper generation limits, are scheduled to provide reserves to balance load fluctuations. Here, the resulting model is a nonsmooth nonconvex two-stage stochastic program, and we use a stochastic approximation method to find stationary solutions to a smooth approximation. Computational results on 6-bus and 118-bus test instances demonstrate that by considering the effects of saturation, our model can yield solutions with lower expected generation costs (at the same target line violation probability level) than those obtained from a model that enforces the affine policy to stay within generator limits with high probability.

17 WIND ENERGY↗

Machine learning and deep learning for mineralogy interpretation and CO 2 saturation estimation in geological carbon Storage: A case study in the Illinois Basin

Carbon capture and storage (CCS) is a promising approach to simultaneously maintaining energy security and reducing carbon dioxide (CO 2 ) emissions under the current energy portfolio that is dominated by fossil fuel energy. Pre-injection formation characterization and post-injection CO 2 monitoring are two critical tasks to guarantee storage efficiency in CCS. The CCS projects in the Illinois Basin, the first large-scale CO 2 injection into saline aquifers in the United States, employed conventional and the latest pulsed neutron logging (PNL) tools for mineralogy interpretation and CO 2 saturation estimation, which provide valuable references for future CCS projects. Because of the inherent fuzziness of petrophysical measurements and complex subsurface heterogeneity, interpreting well-logging data is time-consuming, and its accuracy can be user-biased. In recent years, data-driven methods have been widely used to capture the non-linear patterns between input features and interpretation results. This work applied and evaluated four commonly used machine learning (ML) models, including ridge regression (RR), random forest (RF), gradient boosting regression (GBR), support vector regression (SVR), and one deep learning (DL) model, the artificial neural network (ANN). We optimized the hyperparameters of the four ML models and the DL model using the simulated annealing algorithm and the grid search strategy, respectively. The input features of the mineralogy interpretation models were eleven conventional well-logging parameters, and the label data (i.e., ground truth) were the porosity and volumetric fractions of six minerals, including quartz, feldspar, dolomite, calcite, clay, and iron minerals. The results demonstrated that the GBR and RF models were superior in predicting volumetric fractions of minerals and porosity; label data with low coefficient of variation (CV) values tended to yield better performance. For CO 2 saturation estimation, the RF was the best-performing model, followed by SVR, ANN, GBR, and RR. Furthermore, we conducted feature importance ranking using the permutation importance algorithm and found that the formation sigma and well pressure were the most important features in this study. In conclusion, the study of CCS projects in the Illinois Basin bridges the gap between the limited knowledge and understanding of geological carbon storage and the increasing demand for reliable, cost-effective, and sustainable energy solutions.

58 GEOSCIENCES↗

The effect of original and initial saturation on residual nonwetting phase capillary trapping efficiency

Injection of supercritical carbon dioxide (CO 2 ) into geological formations is a strategy for both atmospheric greenhouse gas reduction (climate change mitigation) and enhanced oil recovery. To understand CO 2 trapping efficiency, the capillary trapping behaviors that immobilize subsurface fluids were analyzed at the pore-scale using pairs of proxy fluids representing the range of in situ nonwetting and wetting fluid properties encountered in geologic storage reservoirs. The pairs of fluids were cycled through imbibition and drainage processes using a flow cell apparatus containing a sintered glass bead column. Computed x-ray microtomography (microCT) was used to identify immobilized nonwetting fluid after imbibition and drainage events. From microCT images, the trapped residual (post-secondary imbibition) nonwetting phase was spatially correlated to both the original (post-primary imbibition) and the initial (post-primary drainage) nonwetting phase; this relationship is referred to here as the original saturation dependence (SO-dependence) and initial saturation dependence (SI-dependence), respectively. Significant trends of decreasing SO- and SI-dependence with increasing wetting and nonwetting fluid phase viscosities were observed. This finding implies that the amount of CO 2 injected and ultimately trapped is dependent on the nonwetting phase (e.g. oil or gas) already present in the formation, as well as on the manner in which supercritical CO 2 is initially injected, which are factors not considered in current trapping models. To explore the potential effect that varying viscosity and interfacial tension (IFT) might have on this process, we report on a variety of fluid pairs with different viscosities and IFT.

42 ENGINEERING↗

Modeling the dynamic response of rock masses with multiple compliant fluid saturated joint sets - Part II: Continuum modeling

Part I of this two-part paper describes enhancements of a mesoscale model to include: many non-persistent fluid saturated joints with finite area, normal compliance with hysteresis and the weakening effects of angular deviation in the joint sets and of fluid pressure in the joints. The objective of this Part II is to use the results of the mesoscale model as simulated data to inspire functional forms in a continuum model for rock masses with multiple compliant fluid saturated joint sets which include these enhancements. Also, an algorithm is proposed for reducing the number of joint sets being modeled. Examples are presented here which show that the continuum model can predict similar anisotropic response to that observed in mesoscale simulations. In contrast with the mesoscale model, the continuum model can be used for large scale simulations of real applications, smoothly transitioning from coupled inelastic thermomechanical response near an explosive source to orthotropic elastic response far from the source.

58 GEOSCIENCES↗