Search NASA⌕ Search

SEARCH · Search NASA

Results for “LOGARITHM”

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 811 records · Page 45

Strain-Based Assessment of Shale Caprock during Cyclic Underground Hydrogen Storage

Successful large-scale underground hydrogen storage (UHS) in depleted gas reservoirs depends on the integrity of the overlying caprock to prevent hydrogen loss during cyclic injection and depletion. Prior studies on crushed Marcellus shale, a potential caprock, indicate that cyclic hydrogen injection and depletion induces microstructural changes, increasing porosity and permeability. However, the extent of these changes in intact shale remains unclear. This study presents a strain-based experimental approach to quantify volumetric strain evolution in intact Marcellus shale matrix under unconstrained stress conditions. A quadrant-shaped shale sample without visible fractures underwent eight hydrostatic pore-pressure cycles (injection to 1500 psi and depletion to 500 psi in 250 psi steps). Linear strain gauges measured strain in three orthogonal directions. Results indicate progressive plastic strain accumulation, leading to an ∼12% increase in matrix porosity after eight cycles, with an estimated 19% increase after 30 cycles. This porosity increase follows a logarithmic trend, suggesting a diminishing effect in later cycles. Additionally, permeability and diffusive mass flux are projected to rise by ∼70% over 30 cycles, enhancing hydrogen migration risk. The shale matrix also exhibited mechanical stiffening over successive cycles, limiting large-scale deformation but not preventing porosity enhancement. A new parameter, α, was introduced to characterize shale sensitivity to cyclic loading, aiding UHS caprock assessments. These findings underscore the necessity of incorporating cyclic loading effects in UHS site selection and operational strategies to ensure long-term storage integrity.

08 HYDROGEN↗

Reduced Chemistry for Butanol Isomers at Engine-Relevant Conditions

Butanol has received significant research attention as a second-generation biofuel in the past few years. In the present study, skeletal mechanisms for four butanol isomers were generated from two widely accepted, well-validated detailed chemical kinetic models for the butanol isomers. The detailed models were reduced using a two-stage approach consisting of the directed relation graph with error propagation and sensitivity analysis. During the reduction process, issues were encountered with pressure-dependent reactions formulated using the logarithmic pressure interpolation approach; these issues are discussed and recommendations are made to avoid ambiguity in its future implementation in mechanism development. The performance of the skeletal mechanisms generated here was compared with that of detailed mechanisms in simulations of autoignition delay times, laminar flame speeds, and perfectly stirred reactor temperature response curves and extinction residence times, over a wide range of pressures, temperatures, and equivalence ratios. Good agreement was observed between the detailed and skeletal mechanisms, demonstrating the adequacy of the resulting reduced chemistry for all the butanol isomers in predicting global combustion phenomena. In addition, the skeletal mechanisms closely predicted the time-histories of fuel mass fractions in homogeneous compression-ignition engine simulations. The performance of each butanol isomer was additionally compared with that of a gasoline surrogate with an antiknock index of 87 in a homogeneous compression-ignition engine simulation. The gasoline surrogate was consumed faster than any of the butanol isomers, with tert -butanol exhibiting the slowest fuel consumption rate. As a result, while n -butanol and isobutanol displayed the most similar consumption profiles relative to the gasoline surrogate, the two literature chemical kinetic models predicted different orderings.

09 BIOMASS FUELS↗

Binding Affinity Prediction by Pairwise Function Based on Neural Network

In this paper, we present a new approach to estimate the binding affinity from given three-dimensional poses of protein–ligand complexes. In this scheme, every protein–ligand atom pair makes an additive free-energy contribution. The sum of these pairwise contributions then gives the total binding free energy or the logarithm of the dissociation constant. The pairwise contribution is calculated by a function implemented via a neural network that takes the properties of the two atoms and their distance as input. The pairwise function is trained using a portion of the PDBbind 2018 data set. The model achieves good accuracy for affinity predictions when evaluated with PDBbind 2018 and with the CASF-2016 benchmark, comparing favorably to many scoring functions such as that of AutoDock Vina. The framework here may be extended to incorporate other factors to further improve its accuracy and power.

59 BASIC BIOLOGICAL SCIENCES↗

Thermodynamics–Structure–Dynamics Correlations and Nonuniversal Effects in the Elastically Collective Activated Hopping Theory of Glass-Forming Liquids

We employ the microscopic Elastically Collective Nonlinear Langevin Equation (ECNLE) theory of activated dynamics in combination with crystal-avoiding simulations to study four inter-related questions for metastable monodisperse hard sphere fluids. The first is how significantly improved integral equation theory structural input (Modified-Verlet (MV) closure) changes the dynamical predictions of ECNLE theory. The main consequence is a modest enhancement of the importance of the collective elastic barrier relative to its local cage contribution, which increases the alpha relaxation time and fragility relative to prior results based on the Percus–Yevick closure. Second, ECNLE-MV theory predictions for the alpha time and self-diffusion constant in the metastable regime are quantitatively compared to our new simulations. The small adjustment of a numerical prefactor that enters the collective elastic barrier leads to quantitative agreement over three decades. Third, using the more accurate MV structural input, ECNLE theory is shown to predict thermodynamics–structure–dynamics “correlations” based on various long and short wavelength scalar properties all related to static two-point collective density fluctuations. The logarithm of the alpha relaxation time scales as a power law with these scalar metrics with an exponent that is significantly lower in the less dense noncooperative activated regime compared to the very dense highly cooperative regime. However, the discovered correlation of activated relaxation with a thermodynamic property (dimensionless compressibility) is not causal in ECNLE theory, but rather reflects a strong connection between the local structural quantities that quantify kinetic constraints in the theory with the amplitude of long wavelength density fluctuations. Fourth, the consequences of chemically specific nonuniversalities associated with the onset condition and relative importance of collective elasticity are studied. In conclusion, the predicted thermodynamics–structure–dynamics correlations are found to be robust, albeit with nontrivial shifts of the onset condition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Surface Phases and Surface Freezing in an Ionic Liquid

Room temperature ionic liquids (RTILs), a novel class of liquid salts, are intensively studied for their basic science and numerous emerging applications. When undercooled, RTILs comprising long alkyl chains often exhibit liquid crystal (LC) bulk phases. However, only one molecular-resolution experimental structure study was published for their LC surface phases. With this being said, we measured the temperature evolution of another LC surface phase, using surface specific Å-resolution X-ray methods. This phase’s existence range, 90 °C, much exceeds the corresponding bulk phase’s 3 °C. Its thickness, L , confirms the theory-predicted logarithmic temperature dependence, with an amplitude equaling the bulk correlation length. Surprisingly, at L’s divergence temperature, a ~20 Å thick, hexagonally packed, crystalline monolayer forms at, and fully covers, the sample’s surface. Consequently, it is identified as a surface-frozen Langmuir–Gibbs film and fundamentally differs from the only reported RTIL surface crystal, a Coulomb-dominated, four-layer, island phase, covering only 5%–15% of the surface.

36 MATERIALS SCIENCE↗

Optically Generated Free-Carrier Collection from an All Single-Walled Carbon Nanotube Active Layer

Semiconducting single-walled carbon nanotubes' (SWCNTs) broad absorption range and all-carbon composition make them attractive materials for light harvesting. We report photoinduced charge transfer from both multichiral and single-chirality SWCNT films into atomically flat SnO 2 and TiO 2 crystals. Higher-energy second excitonic SWCNT transitions produce more photocurrent, demonstrating carrier injection rates are competitive with fast hot-exciton relaxation processes. A logarithmic relationship exists between photoinduced electron-transfer driving force and photocarrier collection efficiency, becoming more efficient with smaller diameter SWCNTs. Photocurrents are generated from both conventional sensitization and in the opposite direction with the semiconductor under accumulation and acting as an ohmic contact with only the p-type nanotubes. Finally, we demonstrate that SWCNT surfactant choice and concentration play a large role in photon conversion efficiency and present methods of maximizing photocurrent yields.

36 MATERIALS SCIENCE↗

Overlap Concentration in Salt-Free Polyelectrolyte Solutions

For strongly charged polyelectrolytes in salt-free solutions, we use molecular dynamics simulations of a coarse-grained bead-spring model to calculate overlap concentrations c* and chain structure for polymers containing N = 10 to 1600 monomers. Over much of this range, we find that the end-to-end distance R* at c* increases faster than linearly with increasing N, as chains at the overlap concentration approach strongly extended conformations. This trend results in the overlap concentration c* decreasing as a stronger function of N than the classical prediction c* ~ N –2 . Finally, this stronger dependence can be fit either by a logarithmic correction to scaling or by an apparent scaling c* ~ N –m , with m > 2.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Signature of Many-Body Localization of Phonons in Strongly Disordered Superlattices

Many-body localization (MBL) has attracted significant research attention due to its immunity to thermalization, role in logarithmic entanglement entropy growth, and opportunities to reach exotic quantum orders. However, experimental realization of MBL in solid-state systems has remained challenging. Here we report evidence of a possible phonon MBL phase in disordered GaAs/AlAs superlattices. Through grazing-incidence inelastic X-ray scattering, we observe a strong deviation of the phonon population from equilibrium in disordered samples doped with ErAs nanodots at low temperature, signaling a departure from thermalization. This behavior occurs within finite phonon energy and wavevector windows, suggesting a localization-thermalization crossover. We support our observation by proposing a theoretical model for the effective phonon Hamiltonian in disordered superlattices, and showing that it can be mapped exactly to a disordered 1D Bose-Hubbard model with a known MBL phase. The large on-site repulsion qualitatively explains the experimental findings. Furthermore, our work provides mode-resolved experimental evidence of phonon localization, potentially extending the scope of MBL to disordered solid-state systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Ionic-Based Electrochemical Gas Sensors for Low-Cost, High-Sensitivity SO2 Detection

Sulfur dioxide (SO2) is a toxic gas associated with adverse health and environmental effects that necessitate reliable monitoring techniques. Here, we report the development of an all-solid-state electrochemical sensor utilizing a lithium borate (Li3BO3) solid electrolyte capable of subppm of SO2 detection. While subppm of SO2 sensing has been previously demonstrated in other solid-state electrolyte systems─such as stabilized zirconia, natrium super ionic conductors (NASICON) under mixed-potential conditions─here we establish Li3BO3 as an alternative solid electrolyte enabling equilibrium potentiometric sensing in an all-solid architecture. This sensor demonstrates a detection limit of at least 0.25 ppm, surpassing the human-olfactory threshold and meeting the rigorous requirements for industrial and personal monitoring applications. The sensing mechanism relies on the formation of Li2SO4 on the electrode surface, as evidenced by multimodal characterization techniques, including Raman spectroscopy, scanning electron microscopy (SEM), and scanning transmission electron microscopy (STEM). The strong linear correlation between the open-circuit potential (OCV) and the logarithm of SO2 concentration between 0.25 and 2 ppm indicates that the response is Nernstian in nature.

Lagunas, Francisco (ORCID:000000026377683X)↗

Uncondensed Graphitic Carbon Nitride on Reduced Graphene Oxide for Oxygen Sensing via a Photoredox Mechanism

Melon, a polymeric, uncondensed graphitic carbon nitride with a two-dimensional structure, has been coupled with reduced graphene oxide (rGO) to create an oxygen chemiresistor sensor that is active under UV photoactivation. Oxygen gas is an important sensor target in a variety of areas including industrial safety, combustion process monitoring, as well as environmental and biomedical fields. Because of the intimate electrical interface formed between melon and rGO, charge transfer of photoexcited electrons occurs between the two materials when under UV (λ = 365 nm) irradiation. A photoredox mechanism wherein oxygen is reduced on the rGO surface provides the basis for sensing oxygen gas in the concentration range 300–100,000 ppm. The sensor response was found to be logarithmically proportional to oxygen gas concentration. DFT calculations of a melon-oxidized graphene composite found that slight protonation of melon leads to charge accumulation on the rGO layer and a corresponding charge depletion on the melon layer. As a result, this work provides an example of a metal-free system for solid–gas interface sensing via a photoredox mechanism.

14 SOLAR ENERGY↗

Kinetic Isotope Effects: Interpretation and Prediction Using Degrees of Rate Control

Kinetic isotope effects (KIEs) have been used for decades in catalysis research as a tool for clarifying reaction mechanisms. Significant primary kinetic isotope effects have usually been interpreted as being a result of isotope substitution at a site of bond breaking (or forming) in the rate-determining step in the reaction mechanism. However, quantitative analysis of the magnitude of the KIE in complex multistep reaction mechanisms is seldom reported. We prove here that the logarithm of the rate ratio for two isotopes is the weighted average over all species in the mechanism of their standard-state free energy difference between the two isotopes, divided by RT. The weighting factor is the degree of rate control (DRC) for that species (e.g., transition state, intermediate, reactant) when the rate is measured separately for each isotope. It is instead the degree of selectivity ratio control (DSRC) when the KIE is measured as the product selectivity in a parallel competition between two isotopes within the same reactant molecule. Since only a few species have non-zero DRCs (or DSRCs) for most reactions, only these few contribute to this weighted average and the KIE. We show that this provides a simple way to interpret and quantitatively predict kinetic isotope effects that is powerful in the insights it provides, allowing one to evaluate directly which species contribute most to the KIE. By applying it to H/D KIEs in several example mechanisms, we further show that the traditional way of interpreting KIEs that focuses only on the rate-determining step can easily lead to misunderstanding of KIE and the reaction mechanism. Finally, this highlights the importance to consider the effect of isotope substitution on the energies of all species with large DRCs (i.e, those whose energies are kinetically relevant). This method also offers opportunities for quantitative validation of mechanism-based microkinetic models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Stability of Nd Hydroxyl Complexes at Near Neutral to Alkaline pH from 25 to 75 °C: Implications for Monazite Solubility in Hydrothermal Aqueous Fluids

Speciation of rare earth elements (REE) in aqueous fluids controls their mobilization during fluid-rock interaction. Thermodynamic modeling provides important insights into the factors controlling REE mobility and fractionation in aqueous fluids, in which aqueous REE hydroxyl complexes are particularly important across a broad pH range. However, the stability of these REE hydroxyl complexes remains poorly determined due to limited experimental data above ambient temperature. In this work, UV–vis spectrophotometry experiments were conducted to determine the cumulative hydrolysis constants for neodymium (Nd) from 25 to 75 °C in near-neutral to alkaline solutions. The color indicator m -cresol purple was used to determine in situ pH. Alkaline NaOH-bearing aqueous solutions were doped with varying initial Nd concentrations (0 to ∼0.155 mmol/kg), which resulted in the release of protons (H + ) and pH decrease during hydrolysis of Nd 3+ according to Nd 3+ + n H 2 O = Nd(OH) n 3– n + H + . The average OH – ligand number coordinated to Nd 3+ was found to increase from 1.0 to 1.6 at 25 °C and from 1.2 to 2.8 at 75 °C over a pH range from 6.3 to 9.0. The measured speciation shows an increased predominance of Nd(OH) 3 0 over the Nd(OH) 2 + and Nd(OH) 2+ species with increased temperature and pH. The increasing cumulative formation constants (log*β° n , n = 1 to 3) retrieved from 25 to 75 °C differ by 0.1 to 1.9 logarithmic units in comparison to existing thermodynamic databases. These updated thermodynamic data have important implications for geochemical modeling of the speciation of Nd hydroxyl complexes and the solubility of monazite as a function of pH.

58 GEOSCIENCES↗

In Situ Oxidation Studies of High-Entropy Alloy Nanoparticles

Although high-entropy alloys (HEAs) have shown tremendous potential for elevated temperature, anticorrosion, and catalysis applications, little is known on how HEA materials behave under complex service environments. Herein, we studied the high-temperature oxidation behavior of Fe 0.28 Co 0.21 Ni 0.20 Cu 0.08 Pt 0.23 HEA nanoparticles (NPs) in an atmospheric pressure dry air environment by in situ gas-cell transmission electron microscopy. It is found that the oxidation of HEA NPs is governed by Kirkendall effects with logarithmic oxidation rates rather than parabolic as predicted by Wagner's theory. Further, the HEA NPs are found to oxidize at a significantly slower rate compared to monometallic NPs. The outward diffusion of transition metals and formation of disordered oxide layer are observed in real time and confirmed through analytical energy dispersive spectroscopy, and electron energy loss spectroscopy characterizations. Localized ordered lattices are identified in the oxide, suggesting the formation of Fe 2 O 3 , CoO, NiO, and CuO crystallites in an overall disordered matrix. Hybrid Monte Carlo and molecular dynamics simulations based on first-principles energies and forces support these findings and show that the oxidation drives surface segregation of Fe, Co, Ni, and Cu, while Pt stays in the core region. The present work offers key insights into how HEA NPs behave under high-temperature oxidizing environment and sheds light on future design of highly stable alloys under complex service conditions.

36 MATERIALS SCIENCE↗

Discovery of a Wurtzite-like Cu 2 FeSnSe 4 Semiconductor Nanocrystal Polymorph and Implications for Related CuFeSe 2 Materials

I 2 -II-IV-VI 4 and I-III-VI 2 semiconductor nanocrystals have found applications in photovoltaics and other optoelectronic technologies because of their low toxicity and efficient light absorption into the near-infrared. Herein, we report the discovery of a metastable wurtzite-like polymorph of Cu 2 FeSnSe 4 , a member of the I 2 -II-IV-VI 4 family of semi-conductors containing only Earth-abundant metals. Density functional theory calculations on this metastable poly-morph of Cu 2 FeSnSe 4 indicate that it may be a superior semiconductor for solar energy and optoelectronics applications compared to the thermodynamically preferred stannite polymorph, since the former displays a sharper dispersion of energy levels near the conduction band minimum that can enhance electron mobility and suppress hot electron cooling. The experimental optical band gap was measured by the inverse logarithmic derivative method to be direct, in agreement with theory, and in the range of 1.48-1.59 eV. Mechanistic studies reveal that this metastable phase de-rives from intermediate Cu 3 Se 2 nanocrystals that serve as a structural template for the final hexagonal wurtzite-like product. We compare the chemistry of wurtzite-like Cu 2 FeSnSe 4 to the related CuFeSe 2 material system. Our experimental and computational comparisons between Cu 2 FeSnSe 4 and CuFeSe 2 help explain both the crystal chemistry of CuFeSe 2 that prevents it from forming wurtzite-like polymorphs and the essential role of Sn in stabilizing the metastable structure of Cu 2 FeSnSe 4 . Furthermore, this work provides insight into the importance of elemental composition when designing syntheses for metastable materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An Analysis of the Solution Structure and Signaling Mechanism of LovK, a Sensor Histidine Kinase Integrating Light and Redox Signals

Flavin-binding LOV domains are broadly conserved in plants, fungi, archaea, and bacteria. These {approx}100-residue photosensory modules are generally encoded within larger, multidomain proteins that control a range of blue light-dependent physiologies. The bacterium Caulobacter crescentus encodes a soluble LOV-histidine kinase, LovK, that regulates the adhesive properties of the cell. Full-length LovK is dimeric as are a series of systematically truncated LovK constructs containing only the N-terminal LOV sensory domain. Nonconserved sequence flanking the LOV domain functions to tune the signaling lifetime of the protein. Size exclusion chromatography and small-angle X-ray scattering (SAXS) demonstrate that the LOV sensor domain does not undergo a large conformational change in response to photon absorption. However, limited proteolysis identifies a sequence flanking the C-terminus of the LOV domain as a site of light-induced change in protein conformation and dynamics. On the basis of SAXS envelope reconstruction and bioinformatic prediction, we propose this dynamic region of structure is an extended C-terminal coiled coil that links the LOV domain to the histidine kinase domain. To test the hypothesis that LOV domain signaling is affected by cellular redox state in addition to light, we measured the reduction potential of the LovK FMN cofactor. The measured potential of -258 mV is congruent with the redox potential of Gram-negative cytoplasm during logarithmic growth (-260 to -280 mV). Thus, a fraction of LovK in the cytosol may be in the reduced state under typical growth conditions. Chemical reduction of the FMN cofactor of LovK attenuates the light-dependent ATPase activity of the protein in vitro, demonstrating that LovK can function as a conditional photosensor that is regulated by the oxidative state of the cellular environment.

59 BASIC BIOLOGICAL SCIENCES↗

Comparison of Energy Efficiency and Power Density in Pressure Retarded Osmosis and Reverse Electrodialysis

Pressure retarded osmosis (PRO) and reverse electrodialysis (RED) are emerging membrane-based technologies that can convert chemical energy in salinity gradients to useful work. The two processes have intrinsically different working principles: controlled mixing in PRO is achieved by water permeation across salt-rejecting membranes, whereas RED is driven by ion flux across charged membranes. This study compares the energy efficiency and power density performance of PRO and RED with simulated technologically available membranes for natural, anthropogenic, and engineered salinity gradients (seawater-river water, desalination brine-wastewater, and synthetic hypersaline solutions, respectively). The analysis shows that PRO can achieve both greater efficiencies (54-56%) and higher power densities (2.4-38 W/m(2)) than RED (18-38% and 0.77-1.2 W/m(2)). The superior efficiency is attributed to the ability of PRO membranes to more effectively utilize the salinity difference to drive water permeation and better suppress the detrimental leakage of salts. On the other hand, the low conductivity of currently available ion exchange membranes impedes RED ion flux and, thus, constrains the power density. Both technologies exhibit a trade-off between efficiency and power density: employing more permeable but less selective membranes can enhance the power density, but undesired entropy production due to uncontrolled mixing increases and some efficiency is sacrificed. When the concentration difference is increased (i.e., natural -> anthropogenic -> engineered salinity gradients), PRO osmotic pressure difference rises proportionally but not so for RED Nernst potential, which has logarithmic dependence on the solution concentration. Because of this inherently different characteristic, RED is unable to take advantage of larger salinity gradients, whereas PRO power density is considerably enhanced. Additionally, high solution concentrations suppress the Donnan exclusion effect of the charged RED membranes, severely reducing the permselectivity and diminishing the energy conversion efficiency. This study indicates that PRO is more suitable to extract energy from a range of salinity gradients, while significant advancements in ion exchange membranes are likely necessary for RED to be competitive with PRO.

Yip, NY↗

Surface Platinum Electrooxidation in the Presence of Oxygen

Understanding the influences of Pt surface oxides is crucial for elucidating the oxygen reduction reaction (ORR) in fuel cells, particularly with an eye toward the design of next-generation Pt-based electrocatalysts with improved activity. Although gaseous O{sub 2} is always present during the ORR, the majority of previous Pt surface oxide studies were conducted in the absence of O{sub 2} due to experimental limitations. In this study, multiple in situ techniques were applied to study the ORR on platinum in the presence of O{sub 2}. The thin channel flow Pt electrode and the electrochemical quartz microbalance (EQCM) techniques on Pt polycrystalline electrodes suggest that the influence of O{sub 2} in the electrolyte on mass change and charge transfer is negligible. The oxide formation followed a logarithm-growth behavior initiating as early as 0.2 s after potential steps. This suggests that no slow conversion of oxide species (e.g., between OH{sub ads} and O{sub ads}) takes place. Despite the negligible effect of O{sub 2} on the measured oxide coverage, steady-state X-ray absorption spectroscopy (XAS) measurements conducted on dispersed Pt nanoparticles suggested that, only when under O{sub 2}-sparging, place exchange between adsorbed oxide(s) and the Pt surface layer(s) initiated at potentials as low as 0.75 V. This is significantly lower than that observed in an O{sub 2}-free electrolyte (>1.1 V). The effects of O{sub 2} on adsorbed oxide species illustrate the complexities of applying model systems to the real world in order to arrive at a comprehensive description of both the oxygen reduction and the Pt dissolution processes.

30 DIRECT ENERGY CONVERSION↗

Long-Time Relaxation Induced by Dynamic Forcing in Geomaterials

We introduce a theoretical model and experimental evidence of the long-time relaxation process (slow dynamics) in rocks and other geomaterials following a dynamic wave excitation, at scales ranging from the laboratory to the Earth. The model is based on the slow recovery of an ensemble of grain contacts and asperities broken by a mechanical impact. It includes an Arrhenius-type equation for recovery of the metastable, broken contacts. The model provides a characteristic size of the broken contacts (order 10–9 m) and predicts that their number increases with impact amplitude. Theoretical results are in sufficient agreement with the laboratory and field data in that they suggest both the logarithmic law of recovery rate and deviations from this law.

58 GEOSCIENCES↗