Search NASA⌕ Search

SEARCH · Search NASA

Results for “equilibrium constant”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Microwave-driven heterogeneous catalysis for activation of dinitrogen to ammonia under atmospheric pressure

This paper presents an innovative approach to producing energy-dense, carbon–neutral liquid ammonia as a means for carrying energy. This approach synergistically integrates microwave reaction chemistry with novel heterogeneous catalysis that decouples dinitrogen activation from high-temperature and high-pressure reactions, altering reaction pathways and increasing ammonia formation rate. Results presented here demonstrate that ammonia synthesis can be conducted at 280 ? and ambient pressure to achieve ~1 mmol ammonia/g cat/h over supported ruthenium catalyst systems utilizing microwave irradiation. It is further shown that adding promoter ions such as potassium, cerium, and barium significantly improves the ammonia production rate over undoped ruthenium-based catalysts. This effect could be attributed to enhanced dielectric loss processes that lead to stronger microwave absorption by the catalyst. Measurement of the equilibrium constant under microwave conditions showed a higher ammonia yield than under thermal equilibrium conditions for both the iron- and ruthenium-based catalysts. Finally, this study also illustrates the advantages of using a variable-frequency microwave reactor for ambient-pressure ammonia synthesis. Mechanistically, investigators believed that the oscillating electric fields of the radiation can couple with adsorbed nitrogen on the surface and accelerate its dissociation. Since dinitrogen dissociation on the surface is rate limiting, this effectively accelerates the reaction.

Hu, Jianli (John)↗
NASA NTRS1975

Proton transfer reactions from H3/+/ ions to N2, O2, and CO molecules

The rate constants for proton transfer from H3(+) ions to N2, O2, and CO have been measured as function of hydrogen-buffer-gas partial pressure. The rate constant for proton transfer from H3(+) to N2 shows a very large pressure dependence, increasing from 1.0 by 10 to the -9th power cu cm/s at low H2 partial pressures to 1.7 by 10 to the -9th power cu cm/s at high H2 partial pressures. The rate constants for proton transfer from H3(+) to O2 and CO are constant with partial pressure of H2. The roles of excess vibrational energy in H3(+) ions and of equilibrium between forward and back reactions are discussed. Back reaction is observed only for the reaction of H3(+) ions with O2, and an equilibrium constant of 2.0 (plus or minus 0.4) at 298 K has been determined. From these data, the proton affinity of O2 is deduced to be 0.47 (plus or minus 0.11) kcal/mole higher than that of H2.

Kim, J. K.↗
NASA NTRS1997

Kinetics and Thermochemistry of the Cl((sup 2)P(sub J)) + C2Cl4 Association Reaction

A laser flash photolysis-resonance fluorescence technique has been employed to study the kinetics of the Cl(sup 2)P(sub j) + C2Cl4 association reaction as a function of temperature (231-390 K) and pressure (3-700 Torr) in nitrogen buffer gas. The reaction is found to be in the falloff regime between third and second order over the range of conditions investigated, although the second-order limit is approached at the highest pressures and lowest temperatures. At temperatures below 300 K, the association reaction is found to be irreversible on the experimental time scale of approximately 20 m-s. The kinetic data at T is less than 300 K have been employed to obtain falloff parameters in a convenient format for atmospheric modeling. At temperatures above 330 K, reversible addition is observed, thus allowing equilibrium constants for C2Cl5 formation and dissociation to be determined. Second- and third-law analyses of the equilibrium data lead to the following thermochemical parameters for the association reaction: Delta-H(298) = -18.1 +/- 1.3 kcal/mol, Delta-H(0) = -17.6 +/- 1.3 kcal/mol, and Delta-S(298) = -27.7 +/- 3.0 cal/mol.K. In conjunction with the well-known heats of formation of Cl((sup 2)P(sub j)) and C2Cl4 the above Delta-H values lead to the following heats of formation for C2Cl5, at 298 and 0 K: Delta-H(f,298) = 8.0 +/- 1.3 kcal/mol and Delta-H(f,0) = 8.1 +/- 1.5 kcal/mol. The kinetic and thermochemical parameters reported above are compared with other reported values, and the significance of reported association rate coefficients for understanding tropospheric chlorine chemistry is discussed.

Nicovich, J. M.↗
NASA NTRS1997

Kinetics and Thermochemistry of the Br((sup 2)P3/2) + NO2 Association Reaction

A laser flash photolysis-resonance fluorescence technique has been employed to study the kinetics of the Br((sup 2)P3/2) + NO2 association reaction as a function of temperature (259-432 K) pressure (12.5 - 700 Torr), and buffer gas identity (He, Ar, H2, N2, CO2, CF4, SF6). The reaction is found to be in the falloff regime between third and second order over the entire range of conditions investigated. At temperatures below 350 K, the association reaction is found to be irreversible on the time scale of the experiment (approximately 30 ms). At higher temperatures reversible addition is observed, allowing equilibrium constants for BrNO2 formation and dissociation to be determined. Second- and third-law analyses of the equilibrium data are in only fair agreement and lead to the following thermochemical parameters for the association reaction: Delta-H(298) = 19.6 +/- 1.7 kcal/mol, Delta-H(0) = -18.6 +/- 2.0 kcal/mol, Delta-S(298) = 29.3 +/- 4.2 cal/mol/K, Delta-H(sub f)(sub 298)(BrNO2) = 17.0 +/-1.8 kcal/mol(uncertainties are 2 sigma estimates of absolute accuracy). The value for Delta-H(0) determined in this study has been employed to calculate k(sub 0)(sup SC), the low-pressure third-order rate coefficient in the strong collision limit, by using the method of Troe; calculated values of k(sub 0)(sup SC) are inconsistent with experimental results unless Delta-H(0) is assigned a value near the lower limit derived from analysis of the high-temperature approach to equilibrium data, i.e. delta-H(0) approximately equals -16.6 kcal/mol. A potential source of systematic error in the calculation of both k(sub 0)(sup SC) and the absolute entropy of BrNO2 results from the complete lack of knowledge of the energies and degeneracies of the electronic states of BrNO3. The procedure developed by Troe and co-workers has been employed to extrapolate experimental falloff curves to the low- and high-pressure limits. Derived values for k(sub 0)(M,298K) in units of 10(exp -31) cm(exp 6)/sq molecule/s range from 2.75 for M = He to 6.54 for M = CO2; 2 sigma uncertainties are estimated to be +/- 20%. Values for k(sub 0)(N2,T) in units of 10(exp -31) cm(exp 6)/sq molecule/s are 5.73 at 259 K, 4.61 at 298 K, and 3.21 at 346 K; the observed temperature dependence for k(sub 0)(N2,T) is consistent with the theoretical temperature dependence for Beta(sub c)k(sub 0)(sup SC). Values for k(sub infinity)(T) in units of 10(exp -11) cu cm/molecule/s are 2.86 at 259 K, 3.22 at 298 K, and 3.73 at 346 K; 2 sigma uncertainties are estimated to be a factor of 2. Approximate falloff parameters in a convenient format for atmospheric modeling are also derived.

Kreutter, K. D.↗
NASA NTRS1999

Structure of the Soot Growth Region of Laminar Premixer Methane/Oxygen Flames

Soot is a dominant feature of hydrocarbon/air flames, affecting their reaction mechanisms and structure. As a result, soot processes affect capabilities for computational combustion as well as predictions of flame radiation and pollution emissions. Motivated by these observations, the present investigation extended past work on soot growth in laminar premixed flames, seeking to evaluate model predictions of flame structure. Xu et al. report direct measurements of soot residence times, soot concentrations, soot structure, gas temperatures and gas compositions for premixed flames similar to those studied by Harris and Weiner and Ramer et al. respectively. It was found that predictions of major stable gas species concentrations based on mechanisms of Leung and Lindstedt and Frenklach and coworkers, were in good agreement with the measurements. The results were also used to evaluate the hydrogen-abstraction/carbon-addition (HACA) soot growth mechanisms of Frenklach and coworkers and Colket and Hall. It was found that these mechanisms were effective using quite reasonable correlations for the steric factors appearing in the theories. The successful evaluation of the HACA mechanism of soot growth in Refs. 1 and 2 is encouraging but one aspect of this evaluation is a concern. In particular, H-atom concentrations play a crucial role in the HACA mechanism and it was necessary to estimate these concentrations because they were not measured directly. These estimates were made assuming local thermodynamic equilibrium between H, and H based on measured temperatures and H2 concentrations and the equilibrium constant data of Kee et al.. This approach was justified by the flame structure predictions; nevertheless, direct evaluation of equilibrium estimates of H-atom concentrations in the soot growth regions of laminar premixed flames is needed to provide more convincing proof of this behavior. Thus, the objective of the present investigation was to complete new measurements of the structure of the soot growth region of laminar premixed flames and to use these results to evaluate whether H and H2 are in thermodynamic equilibrium and to extend the earlier evaluation of predictions of concentrations of major gas species.

Xu, F.↗
DOE OSTI2022

Autoprocessing and oxyanion loop reorganization upon GC373 and nirmatrelvir binding of monomeric SARS-CoV-2 main protease catalytic domain

The monomeric catalytic domain (residues 1–199) of SARS-CoV-2 main protease (MPro 1-199 ) fused to 25 amino acids of its flanking nsp4 region mediates its autoprocessing at the nsp4-MPro 1-199 junction. We report the catalytic activity and the dissociation constants of MPro 1-199 and its analogs with the covalent inhibitors GC373 and nirmatrelvir (NMV), and the estimated monomer-dimer equilibrium constants of these complexes. Mass spectrometry indicates the presence of the accumulated adduct of NMV bound to MPro WT and MPro 1-199 and not of GC373. A room temperature crystal structure reveals a native-like fold of the catalytic domain with an unwound oxyanion loop (E state). In contrast, the structure of a covalent complex of the catalytic domain-GC373 or NMV shows an oxyanion loop conformation (E* state) resembling the full-length mature dimer. These results suggest that the E-E* equilibrium modulates autoprocessing of the main protease when converting from a monomeric polyprotein precursor to the mature dimer.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1991

Thermodynamics of acetylene van der Waals dimerization

Integrated band intensities of the 620/cm absorption in (C2H2)2 are measured by FTIR spectroscopy at constant acetylene pressure between 198 and 273 K. These data, in conjunction with ab initio results for (C2H2)2, are used for the statistical evaluation of the equilibrium constant Kp(T) for acetylene-cluster dimerization. The present results are used to clarify the role of molecular clusters in chemical systems at or near equilibrium, in particular in Titan's stratosphere.

Colussi, A. J.↗
NASA NTRS1990

Modification of the PARC Navier-Stokes Code to predict rocket engine nozzle performance

The PARC2D Navier-Stokes Code was modified to compute the performance parameters for rocket engine nozzles. The perfect gas code was applied to the SSME engine nozzle for inviscid, laminar and turbulent flow. Inviscid computations compare well with Rocketdyne computations. Performance degradation due to the boundary layers is very reasonable. Application of the code to nontraditional nozzle geometries and to low Reynolds nozzles is demonstrated. Modification of the code for equilibrium H2/O2 chemistry is described. Thermodynamic and equilibrium constants are determined from statistical mechanics and the transport properties from exact kinetic theory, using collison integrals determined from appropriate intermolecular potentials. The equilibrium code was used to compute the SSME flowfield. Modifications of the flowfield due to the change of composition are described.

Collins, Frank G.↗
DOE OSTI2023

Tumor‐pH‐value responsive non‐peripheral substituted phthalocyanines: Synthesis, investigation of photophysical and photochemical properties

Cancer is one of the diseases with the highest mortality rate worldwide. Although PDT has recently produced encouraging outcomes, there are still many areas that need to be improved. The first of these is the negative consequences faced by patients treated with the PDT method when exposed to sunlight. For this reason, a new PDT method has been developed in recent years, and it is aimed at using photosensitizer molecules that can be active in acidic conditions. Since the pH values of tumor tissues are more acidic than normal tissues, preparing molecules that act effectively in acidic conditions will allow for more effective results in treating cancer with PDT. In this context, within the scope of this study, 3‐(4‐propionylphenoxy)phthalonitrile ( 1 ) and its non‐peripheral tetra‐substituted phthalocyanine derivatives [( 2 ), ( 3 ), and ( 4 )] were prepared. With these phthalocyanine derivatives, the novel compounds ( 5 ), ( 6 ), and ( 7 ) were synthesized for the first time. The aggregation tendencies of newly synthesized phthalocyanines ( 5–7 ) were investigated in solvent media. The effects of pH changes upon UV–Vis and fluorescence spectra were performed. The electronic and emission spectra of synthesized phthalocyanine derivatives are highly sensitive to pH changes. Formation constant (Log K) values of mono‐ and di‐protonated phthalocyanine forms were calculated by the Henderson–Hasselback equation. The mono‐ and di‐protonated species' equilibrium constants (logK 1 and logK 2 ) were calculated as ~5.0. This value may be promising for pH‐sensitizing photosensitizers. Also, the photophysical and photochemical properties of synthesized metallophthalocyanine derivatives ( 2 ) and ( 5 ) were studied at different pH values. The singlet oxygen quantum yield of ( 2 ) and ( 5 ) was calculated to be 0.78 and 0.81 in DMSO, respectively. When pH = 6.4, that is, tumor‐pH‐values, this value for ( 5 ) has increased to 0.92. The newly synthesized phthalocyanines are suitable photosensitizers for PDT applications, especially with high singlet oxygen quantum yield at pH 6.4.

Genc Bilgicli, Hayriye↗
DOE OSTI2021

Electrochemical Modeling of Iodide Oxidation in Metal-Halide Molten Salts

The oxidation of iodide in NaI-AlBr 3 , NaI-AlCl 3 , and NaI-GaCl 3 molten salts was analyzed using simulation software to extract relevant kinetic parameters. The experimental oxidation potentials were ordered AlCl 3 < AlBr 3 < GaCl 3 , with higher oxidation potentials correlating with softer Lewis acidity of the metal halide. An iodide oxidation and metal halide speciation model was developed and simulated to fit the electrochemical response, enabling determination of electrochemical charge transfer parameters and chemical equilibrium constants. NaI-AlBr 3 displayed the fastest electron transfer rates yet showed the lowest current densities. All salts revealed smaller than expected current densities, explained by equilibrium between various species, where some are not electrochemically active at the studied potentials. These equilibrium reactions are due to the various metal halide species, controlling the reactant concentration of iodide and the resultant current. We hypothesize the electrochemically active iodide species, present as a metal halide monomer (MX 3 I − ), is decreased dramatically from the expected concentration, sequestered as a more stable metal halide dimer species (M 2 X 6 I − ) with a higher oxidation potential. Traditional Tafel analysis of the experimental data supports the validity of the simulations. These results increase understanding of iodide oxidation in low-temperature Lewis acidic molten salts and inform task-specific molten salt design.

Lee, Rose Y. (ORCID:0000000226587220)↗
DOE OSTI2023

On the Thermodynamic Condition for Adsorption Azeotropes

Adsorption azeotropy is a common phenomenon in mixed-gas adsorption equilibria. Since the first literature report of binary adsorption azeotropes in 1933, numerous studies investigated the thermodynamic conditions for adsorption azeotropes but failed to reach definitive conclusions. Based on the generalized Langmuir isotherm model for multicomponent adsorption equilibria which takes into account the vacant site as part of the adsorbed phase, this study presents the thermodynamic condition for adsorption azeotropes as derived from the generalized Langmuir isotherm to be the equality of the ratios of adsorbed phase activity coefficient γi and adsorption equilibrium constant $K$$^{o}_{i}$ for the two adsorbates in the binary 1–2 adsorption system, i.e., γ 1 /$K$$^{o}_{1}$ = γ 2 /$K$$^{o}_{2}$. This adsorption azeotropic condition is analogous to the vapor–liquid equilibrium azeotropic condition, i.e., the equality of the products of liquid phase activity coefficient and saturation vapor pressure Pisat for the two components, i.e., γ 1 $P$$^{sat}_{1}$ = γ 2 $P$$^{sat}_{2}$. We validated the thermodynamic condition for adsorption azeotropes with 14 azeotrope-forming adsorption systems in this study. As a result, we investigated the effects of the pressure, temperature, and adsorbed phase nonideality on azeotrope formation.

42 ENGINEERING↗
NASA NTRS1976

Absolute rate of the reaction of Cl(p-2) with molecular hydrogen from 200 - 500 K

Rate constants for the reaction of atomic chlorine with hydrogen are measured from 200 - 500 K using the flash photolysis-resonance fluorescence technique. The results are compared with previous work and are discussed with particular reference to the equilibrium constant for the reaction and to relative rate data for chlorine atom reactions. Theoretical calculations, using the BEBO method with tunneling, give excellent agreement with experiment.

Whytock, D. A.↗
DOE OSTI2022

Re-evaluation of rate constants for the reaction N 2 H 4 (+ M) $^→_←$ NH 2 + NH 2 (+ M)

Rate constants for the dissociation/recombination reaction N 2 H 4 (+ M) $^→_←$ NH 2 + NH 2 (+ M) are determined by a combination of quantum-chemical calculations and statistical unimolecular rate theory. Between 1100 and 2500 K, limiting low-pressure rate constants for hydrazine dissociation in the bath gas Ar of k 0 = [Ar] 6.1 × 10 20 (T/1000 K) –7.3 exp (-34490 K/T) cm 3 mol –1 s –1 , limiting high-pressure rate constants of k ∞ = 7.6 × 10 16 (T/1000 K) –1.0 exp(-33600 K/T) s –1 , and center broadening factors of the falloff curves (between 1100 and 1600 K in Ar) of F cent = 0.71 exp(- T/1460 K) + 0.29 exp(- T/21 K) + exp(-13400 K/T) were calculated. Using equilibrium constants K C = 1.7 × 10 3 (T/1000 K) –1.5 exp(-33,460 K/T) mol cm –3 , between 300 and 600 K limiting low-pressure rate constants for the reverse recombination of NH 2 radicals in the bath gas N 2 of k rec,0 = [N 2 ] 4.4 × 10 20 (T/400 K) –6.9 exp (-1630 K/T) cm 6 mol –2 s –1 , limiting high-pressure rate constants of k rec,∞ = 4.4 × 10 13 (T/1000 K) 0.44 exp(-140 K/T) cm 3 mol -1 s –1 , and center broadening factors (between 300 and 600 K in N 2 ) of F cent = exp(- T/1130 K) + exp (-10340 K/T) were obtained. A comparison with experimental results for hydrazine dissociation from the literature suggests incomplete falloff extrapolations toward k ∞ and experimental problems in the determination of k 0 at temperatures above about 1600 K. Finally, implications of the present re-evaluated rate constants for the modeling of high temperature ammonia oxidation kinetics are discussed, showing an only small influence of their precise values on the overall properties of the process.

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

Machine Learning in Environmental Chemistry: Application to Surface Complexation Modeling

Environmental chemistry – or biogeochemistry – is the scientific discipline typically invoked when examining and quantifying groundwater or surface water contamination, and nutrient cycling in the environment. Over the last three decades, there have been significant advances in mechanistic model development to describe and predict these complex biogeochemical processes. In particular, surface complexation models (SCMs) have been developed to describe the rock/soil surface reactions of metals and radionuclides, and their partitioning between various mobile species in the aqueous phase or immobile species sorbed on solid surfaces. Often represented by a simplified linear isotherm constant – Kd – in reactive transport models, these reactions play a critical role in many environmental science applications; particularly in contamination risk assessments and nuclear waste disposal performance assessments. In the past several decades, efforts by various institutions across the world have focused on developing SCMs based on datasets from laboratory measurements, including the identification of key parameters such as equilibrium constants.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Kinetics and Energetics of Electron Transfer to Dimer Radical Cations

In this work, spectra of the dimer cations naphthalene (Nap 2 •+ ) and ethene (Ethene 2 •+ ) were measured in liquid dichloromethane (DCM). The spectra peak at very different energies, 1.2 and 3.3 eV. In DCM dimerization stabilizes Nap 2 •+ by Δ$G_\text{d}°$(Nap 2 •+ ) = -218 meV relative to the monomer Nap •+ as determined from the dimerization equilibrium constant. Both dimers can transfer a positive charge to hole acceptor molecules, but for both the rate constants rise more gradually with reaction energetics than do many charge transfer reactions previously studied. A striking observation finds that the rate constant for hole transfer from the Nap 2 •+ dimer to phenanthrene is smaller by two decades than that from biphenyl•+ monomer to Nap, although both reactions have the same –Δ$\textit{G}°$=0.05 eV. A plausible interpretation for these observations is the presence of an energy of reorganization, λ(M 2 ), for the dimer that involves movement apart of the two partners in the dimer. While the dimerization equilibrium cannot be measured for Ethene 2 •+ the charge transfer data imply that both Δ$G_\text{d}°$(Ethene 2 •+ ) and λ(Ethene 2 •+ ), are considerably larger, perhaps by factors of 2-4 than for Nap 2 •+ .

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

Reaction path model of the formation of abiotic immiscible hydrocarbon fluids in subducted carbonated serpentinites, Lanzo Massif (Western Italian Alps)

Fluids generated from subducted slabs participate in the cycling of deep carbon in the crust and upper mantle. In these fluids, aqueous carbon species vary in oxidation state between +IV and -IV depending on whether the fluids are oxidizing or reducing, respectively. Most studies of subduction-zone fluids have focused on oxidized carbon species. However, recent studies of natural samples have demonstrated the occurrence of deep, reducing fluids, generated in both subducted oceanic upper mantle and crustal rocks. CH 4 -H 2 -rich fluid inclusions in subducted carbonated serpentinites have demonstrated the existence of abiotic, immiscible, hydrocarbon fluids at upper mantle conditions. To investigate the formation of such immiscible hydrocarbon fluids during the evolution of subducted carbonated serpentinites, we used equilibrium constants from the Deep Earth Water model to carry out predictive chemical mass transfer modeling to simulate the alteration reactions. A novel feature of the models was the inclusion of an immiscible hydrocarbon fluid containing six components (CH 4,f , C 2 H 6,f , C 3 H 8,f , isoC 4 H 10,f , CO 2,f , H 2,f ). This feature enabled prediction of the formation of a separate immiscible fluid in equilibrium with aqueous species and minerals. We developed a predictive reaction path model of invasive H 2,f reacting with carbonated serpentinites and interstitial aqueous fluids for comparison with the natural samples from the Lanzo Massif, western Italian Alps. Over a range of temperatures and pressures, immiscible hydrocarbon fluids formed in association with altered mineral assemblages. Reaction progress caused the transformation of carbonated serpentinites and the formation of clinopyroxene, brucite, graphite, and hydrocarbon fluids, along with changes of pH, logfO 2 , and aqueous species. CH 4,f was the most abundant hydrocarbon species in all the models. The overall results at 2.0 GPa and 400 to 450 °C were consistent with the natural samples from the Lanzo Massif. Interestingly, large amounts of H 2 O formed due to oxidation of H 2 . More hydrocarbons and H 2 O formed in models with lower fluid/rock mass ratios or with more reactant H 2 . Models at different pressure and temperature conditions showed similar results with some variation in the relative stabilities of aragonite, graphite and olivine solid solution, and associated differences in mineral sequences, hydrocarbon fluids, values of aqueous species, and the final log fO 2 and pH. As a result, our models strongly support the laboratory and field evidence that reduction of carbonated serpentinites by infiltrating H 2 fluids can cause the formation of immiscible, abiotic hydrocarbon fluids in subduction zones.

58 GEOSCIENCES↗
DOE OSTI2020

Onset dynamics of air-water menisci on rock fracture surfaces

Information on menisci dynamics and equilibrium contact angles, θ e , is needed for modeling multiphase flow of fluids in geologic formations. The wettability of fracture faces is relevant to applications such as waterflooding in enhanced oil recovery and evaluation of caprock integrity for geologic carbon sequestration. We investigated the vertical capillary rise of air-water menisci on exposed fracture faces for a sandstone, a shale, and three granites. The experimental design approximated that employed in the Wilhelmy plate method. Replicate fracture faces were prepared by fracturing cylindrical rock cores using the Brazilian method and splitting the fractured cores apart along their fracture planes. Meniscus onset was imaged from a cross-sectional perspective using dynamic neutron radiography, and quantified with change point analysis. An existing model for meniscus onset fitted the experimental data very well (median R 2 = 0.96). Capillary rise was typified by t 2/3 scaling (where t is time), followed by a constant equilibrium height, z e . The average time taken for the menisci to reach z e was 0.39 s, with no differences between rock types. The menisci achieved a greater vertical extent on the fracture faces of the sedimentary rocks as compared to those of the granites. Apparent θ e values were calculated from the best fit estimates of z e using the Wilhelmy plate equation for a smooth and planar vertical surface. Despite the fractal roughness of some of the fracture faces, there was a significant linear relationship (R 2 = 0.80) between the apparent θ e values and intrinsic θ e values measured on polished surfaces of the same rocks using the sessile drop method. The regression equation indicated the apparent θ e values were consistently greater than the intrinsic values. This overestimation may be due to averaging of pixels at the air-water-solid interface, a required step in the image analysis. Further exploration of the range of applicability and measurement accuracy of this proposed new technique would be valuable.

58 GEOSCIENCES↗