Search NASA⌕ Search

SEARCH · Search NASA

Results for “Dynamic Timing Analysis”

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 937 records · Page 52

Development of High-Fidelity and Efficient Modeling Capabilities for Enabling Co-Optimization of Fuels and Multi-Mode Engines

The objective of this closely coordinated research project is to develop improved physical sub-models and innovative numerical algorithms to enable accurate and efficient predictions of multi-mode combustion to support EERE’s Co-Optima program. To this end, a research pro-gram is proposed that seeks • to develop accurate submodels for predicting multi-mode combustion regimes, wall-heat transfer, non-equilibrium plasma ignition, and combustion-mode transition, • to develop dynamic load-rebalancing algorithms and efficient semi-implicit time-integration schemes for enabling exascale computing, • to validate the computational submodels against experiments and direct numerical simulation (DNS) data in conjunction with sensitivity analysis, and • to support the direct transition of the resulting sub-models into other engine-simulation tools through the development of self-contained modules. All sub-models will be integrated into the open-source exascale engine-simulation platform Nek5000.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Transient Catalytic Reaction Analysis Through Signal Defragmentation

The Temporal Analysis of Products (TAP) pulse response technique provides valuable insights into catalytic function and reaction kinetics. However, complex fragmentation patterns in the TAP mass spectrometry signals can complicate precise quantification, particularly when analyzing transient gas flux data typical of TAP experiments. This work demonstrates a standard defragmentation method that deconvolves transient TAP signals while maintaining the temporal resolution of the experiment. First, the integrals of calibration gas fluxes are used to determine the fingerprint fragmentation pattern and construct a fragmentation matrix. This matrix is then used to defragment experimental flux data at each recorded time point via a non-negative least squares regression. The effectiveness of this method is demonstrated using virtual data and control experiments with a TAP reactor system. The defragmentation is then applied to the more complex propane dehydrogenation reaction on a chromia/alumina catalyst, which can contain up to ten significant gas species in the reactor outlet. Initial propane pulsing reveals an induction period during which propane is fully oxidized to CO2, followed by partial reduction to CO. Afterwards, there is a transition in chemistries towards coking and propylene production. Our example illustrates a practical method for the accurate determination of the time-dependent reactant/product concentrations and rates for a thorough analysis of the propane dehydrogenation kinetics. This approach can be broadly applied to any transient mass spectrometry experiment for a better understanding of catalyst-reaction dynamics.

36 - MATERIALS SCIENCE↗

EVOLUTION OF NEAR-SURFACE FLOWS INFERRED FROM HIGH-RESOLUTION RING-DIAGRAM ANALYSIS

Ring-diagram analysis of acoustic waves observed at the photosphere can provide a relatively robust determination of the sub-surface flows at a particular time under a particular region. The depth of penetration of the waves is related to the size of the region, hence the depth extent of the measured flows is inversely proportional to the spatial resolution. Most ring-diagram analysis has focused on regions of extent ∼15° (180 Mm) or more in order to provide reasonable mode sets for inversions. Helioseismic and Magnetic Imager (HMI) data analysis also provides a set of ring fit parameters on a scale three times smaller. These provide flow estimates for the outer 1% (7 Mm) of the Sun only, with very limited depth resolution, but with spatial resolution adequate to map structures potentially associated with the belts and regions of magnetic activity. There are a number of systematic effects affecting the determination of flows from a local helioseismic analysis of regions over different parts of the observable disk, and not all of them are well understood. In this study we characterize those systematic effects with higher spatial resolution so that they may be accounted for more effectively in mapping the temporal and spatial evolution of the flows. Leaving open the question of the mean structure of the global meridional circulation and the differential rotation, we describe the near-surface flow anomalies in time and latitude corresponding to the torsional oscillation pattern in differential rotation and analogous patterns in the meridional cell structure as observed by the Solar Dynamics Observatory/HMI.

79 ASTRONOMY AND ASTROPHYSICS↗

Equilibration and analysis of first-principles molecular dynamics simulations of water

First-principles molecular dynamics (FPMD) simulations based on density functional theory are becoming increasingly popular for the description of liquids. In view of the high computational cost of these simulations, the choice of an appropriate equilibration protocol is critical. In this study, we assess two methods of estimation of equilibration times using a large dataset of first-principles molecular dynamics simulations of water. The Gelman-Rubin potential scale reduction factor [A. Gelman and D. B. Rubin, Stat. Sci. 7, 457 (1992)] and the marginal standard error rule heuristic proposed by White [Simulation 69, 323 (1997)] are evaluated on a set of 32 independent 64-molecule simulations of 58 ps each, amounting to a combined cumulative time of 1.85 ns. The availability of multiple independent simulations also allows for an estimation of the variance of averaged quantities, both within MD runs and between runs. We analyze atomic trajectories, focusing on correlations of the Kohn-Sham energy, pair correlation functions, number of hydrogen bonds, and diffusion coefficient. The observed variability across samples provides a measure of the uncertainty associated with these quantities, thus facilitating meaningful comparisons of different approximations used in the simulations. We find that the computed diffusion coefficient and average number of hydrogen bonds are affected by a significant uncertainty in spite of the large size of the dataset used. Lastly, a comparison with classical simulations using the TIP4P/2005 model confirms that the variability of the diffusivity is also observed after long equilibration times. Complete atomic trajectories and simulation output files are available online for further analysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effect of Biomass Water Dynamics in Cosmic-Ray Neutron Sensor Observations: A Long-Term Analysis of Maize–Soybean Rotation in Nebraska

Precise soil water content (SWC) measurement is crucial for effective water resource management. This study utilizes the Cosmic-Ray Neutron Sensor (CRNS) for area-averaged SWC measurements, emphasizing the need to consider all hydrogen sources, including time-variable plant biomass and water content. Near Mead, Nebraska, three field sites (CSP1, CSP2, and CSP3) growing a maize–soybean rotation were monitored for 5 (CSP1 and CSP2) and 13 (CSP3) years. Data collection included destructive biomass water equivalent (BWE) biweekly sampling, epithermal neutron counts, atmospheric meteorological variables, and point-scale SWC from a sparse time domain reflectometry (TDR) network (four locations and five depths). In 2023, dense gravimetric SWC surveys were collected eight (CSP1 and CSP2) and nine (CSP3) times over the growing season (April to October). The N0 parameter exhibited a linear relationship with BWE, suggesting that a straightforward vegetation correction factor may be suitable (fb). Results from the 2023 gravimetric surveys and long-term TDR data indicated a neutron count rate reduction of about 1% for every 1 kg m−2 (or mm of water) increase in BWE. This reduction factor aligns with existing shorter-term row crop studies but nearly doubles the value previously reported for forests. This long-term study contributes insights into the vegetation correction factor for CRNS, helping resolve a long-standing issue within the CRNS community.

Chemistry↗

Modified Eigen-Decomposition-based Interval Analysis (MEDIA) for Power System Dynamic State Estimation

The Bayesian approach has been used for the dynamic state estimation (DSE) of a power system. However, due to the complexity of noise resources, it is difficult to quantify measurement and process noise using probability density functions (PDFs). To overcome the difficulty, the authors of this paper propose a modified eigen-decomposition-based interval analysis (MEDIA) method, which employs bounds instead of PDFs to quantify the noise, and uses the eigen decomposition method to reduce the negative impact of the overestimation problem. Using the simulation data generated from IEEE 16-machine and IEEE 10-machine systems, it is shown that the proposed MEDIA method can estimate the hard boundaries of dynamic states in real time. Furthermore, comparison with the forward-backward propagation method and the extended set-membership filter also shows that the proposed MEDIA method performs better by providing narrower boundaries in the DSE.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Transient analysis of advancing contact angle measurements on polished rock surfaces

We report that contact angle measurements for gas-liquid-rock systems are important for modeling multi-phase flow and transport in the subsurface. These data are needed in applications such as the extraction of oil and gas resources, geologic sequestration of carbon dioxide, contaminant fate and transport, and aquifer recharge through the vadose zone. Contact angles are frequently measured with the sessile drop method. Previous research has largely ignored the dynamic behavior of sessile droplets on geologic materials. This study investigates the dynamic behavior of sessile water droplets on prepared rock surfaces in the presence of air. Droplet diameter and advancing contact angle were determined at 0.5 s intervals for ~90 s on flat polished disks of Burlington limestone, Crossville sandstone, Mancos shale, Sierra White granite, Vermilion Bay granite, and Westerly granite using a Krüss DSA 30 Drop Shape Analyzer. The droplet diameter and advancing contact angle data sets were nonlinearly regressed against time using two different two-parameter models. The median coefficients of determination for the fits were 0.85 and 0.96, respectively. The resulting parameter estimates were used to compute the apparent equilibrium contact angle, θ e , for each disk following droplet diameter stabilization. Estimates of θ e ranged from 37.2° for Mancos shale to 75.6° for Burlington limestone. Analysis of variance indicated statistically significant differences in θ e between the rock types at the 95% confidence level. The variability of θe on the polished rock surfaces, as quantified by the coefficient of variation (CV) for θ e , varied between ~3 and ~9%; there were no significant differences in CV between the rock types. Neutron radiography indicated changes droplet morphology over time were due to the spontaneous imbibition of water into the rock matrix. Finally, the transient analysis employed in this study permits a more meaningful estimate of the equilibrium contact angle for rocks than taking the initial value or averaging over time as is frequently done.

58 GEOSCIENCES↗

Solvent Effect on the Photoinduced Structural Change of a Phosphorescent Molecular Butterfly

Abstract Photoinduced structural changes (PSC) is one of the fundamental excited‐state dynamic processes, and yet often very challenging to distinguish from competing electronic excited‐state relaxation channels having similar or even comparable timescales. Here, we report a detailed study on the PSC of a pyrazolate bridged platinum(II) binuclear complex, BFPtPZ (C^NPt(μ‐pz′) 2 PtC^N, C^N=2‐(2,4‐difluorophenyl)pyridine, pz′=pyrazolate), a molecular butterfly, using time‐correlated single photon counting measurements at different wavelengths and sample temperatures. Analysis of the results obtained using dichloromethane (DCM) and ethylene carbonate (EC) as solvents enabled us to reveal an unexpected, strong solvent effect on the PSC processes. We show that a rapid PSC process with a characteristic timescale of 323 ps is observed in DCM, which leads to an excitation equilibrium between the ligand center/metal‐to‐ligand charge transfer ( 3 LC/MLCT) and metal‐metal‐to‐ligand charge transfer ( 3 MMLCT) triplet states. The subsequent relaxation from these electronic states to the ground state takes place in several nanoseconds. In contrast, the corresponding PSC process in EC appears slow at all temperatures studied in our experiments and showed no sign of such excitation equilibrium. The observed solvent effect is found to arise from distinct solvent properties including their viscosities and polarities as well as the peculiar electronic excited‐states of the butterfly‐like molecules with charge transfer character.

Ma, Ying‐Zhong↗

The effects of phase decorrelation on the dynamics of the bump-on-tail instability

The dynamics of the bump-on-tail instability has been studied. The novelty of the work is the analysis of how the bump-on-tail dynamics is affected by an extrinsic stochastisation of the phase of the wave-particle interaction; here referred to as phase decorrelation. For this purpose, a nonlinear Monte Carlo model has been developed. When the characteristic time scale for macroscopic phase decorrelation becomes shorter than time scales of nonlinear wave-particle dynamics, the system may be described quasilinearly, with the phase decorrelation being replaced by a quasilinear diffusion coefficient in particle energy. A purely quasilinear Monte Carlo model, which is typically less computationally demanding than the fully nonlinear description due to the reduced dimensionality of phase space, has been developed for comparison. In this paper, parameter regimes, where the nonlinear and the quasilinear descriptions quantitatively agree on a macroscopic level, have been investigated, using combined theoretical and numerical analyses. Qualitative effects on the macroscopic dynamics by the presence of phase decorrelation and/or by structures of the energy distribution function in the proximity of the wave-particle resonance are also studied.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Thermodynamics of optical Bloch equations

Abstract Optical Bloch equations (OBE) describe the coherent exchange of energy between a quantum bit (qubit) and a quasi-resonant driving field in the presence of a thermal bath. Despite it being an ubiquitous process in quantum technologies, a sound thermodynamic analysis is still missing. We hereby provide such an analysis, by deriving the relevant framework from first principles. We start from a complete microscopic description of the qubit-bath system where definitions of heat, work and entropy production are unambiguous. We trace out the bath and coarse-grain the resulting expressions in time, using a methodology similar to the derivation of the dynamical master equation, to derive closed expressions for the first and second law in terms of system properties. Long coarse graining times yield the Floquet Master equation and its already known thermodynamic description. Short coarse-graining times yield instead the OBE and a novel thermodynamic framework which explicitly depends on quantum coherences in the qubit’s energy basis which produce quantum signatures in the heat and entropy production flows. This allows us to characterize a genuinely quantum non-equilibrium situation, where the coherences created by the driving field are continuously erased by the bath. Our findings can be readily extended to larger open quantum systems. They carry the seeds for future thermodynamic analyses of quantum gates and the design of quantum engines in the strong coherent driving regime.

Elouard, Cyril↗

On the optical properties of carbon nanotubes. Part I. A general formula for the dynamical optical conductivity

This paper is the first one in a series of two articles in which we revisit the optical properties of single-walled carbon nanotubes (SWNTs). Produced by rolling up a graphene sheet, SWNTs owe their intriguing properties to their cylindrical quasi-one-dimensional (quasi-1D) structure (the ratio length/radius is experimentally of order of 10{sup 3}). We model SWNT by circular cylinders of small diameters on the surface of which the conduction electron gas is confined by the electric field generated by the fixed carbon ions. The pair-interaction potential considered is the 3D Coulomb potential restricted to the cylinder. To reflect the quasi-1D structure, we introduce a 1D effective many-body Hamiltonian which is the starting-point of our analysis. To investigate the optical properties, we consider a perturbation by a uniform time-dependent electric field modeling an incident light beam along the longitudinal direction. By using Kubo’s method, we derive within the linear response theory an asymptotic expansion in the low-temperature regime for the dynamical optical conductivity at fixed density of particles. The leading term only involves the eigenvalues and associated eigenfunctions of the (unperturbed) 1D effective many-body Hamiltonian and allows us to account for the sharp peaks observed in the optical absorption spectrum of SWNT.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Risk Importance Ranking of Fire Data Parameters to Enhance Fire PRA Model Realism

Fire is historically and analytically a significant contributor to nuclear power plant risk. The level of fire risk and the methods, tools and data for modeling this risk is highly debated by experts. One area of debate is the input data used in fire modeling and how to deal with this data’s high uncertainty. This report outlines initial work performed for determining the key parameters causing this uncertainty and how it propagates into nuclear power plant models. This research paves the way for the development of methods to reduce fire data uncertainty used in modeling. The Nuclear Regulatory Commission has mandated that nuclear power plants perform fire risk modeling. However, there are several issues with the current risk modeling implementation that affect the results. Approved modeling methods can be overly conservative and often do not match plant experience. Also, the data used in the modeling can have high uncertainties and is influenced by expert judgement. To evaluate input data uncertainty, researchers performed an initial review of several fire experiments done at Sandia National Laboratories. Uncertainties for fire data can come from many sources, such as experiment design constraints, environmental conditions, or other plant-specific aspects. There are many different significant and insignificant parameters driving the uncertainty. Additionally, the uncertainty of the different input data used in the fire modeling could have a significant or insignificant effect on the entire plant risk. A four-step methodology was developed to perform Integrated Probabilistic Risk Assessment Importance Ranking. A demonstration case using these steps was set up and three of the four steps were completed in fiscal year (FY) 2019 and the fourth step done FY 2020. These steps are: 1. The qualitative analysis of potential sources was conducted with the following items identified for the demonstration. • Maximum heat release rate • Time to maximum heat release rate • Duration of max heat release rate • Time to decay • Thermal conductivity of concrete • Specific heat of concrete • Density of concrete • Cable jacket thickness 2. A quantitative characterization of dominant sources of uncertainty was performed. A list of distributions and determined values of the dominant sources is shown in Appendix A. 3. A quantitative screening of the potential sources of uncertainty using Morris Elementary Effects Analysis was completed. An experimental model using the physics-based fire modeling tool Fire Dynamics Simulator was developed and coupled with the Risk Analysis Virtual Environment. The Morris analysis identified at least two parameters that can be eliminated as significant contributors (specific heat of concrete and cable jacket thickness). 4. Global importance measure (Global IM) analysis to generate a comprehensive ranking based on their influence on the plant risk. In this research, a moment-independent Global IM is used since it can address (a) uncertainty in the input parameters of the fire model, (b) uncertainty in the risk outputs, and (c) non-linearity and interactions among input parameters in the fire model, more accurately than the correlation-based and variance-based global methods. The observations from the research showed that, depending on the initial and boundary conditions of the fire scenarios, fire-induced damage could have a very small probability and could be dominated by the tail of the uncertainty distribution; hence, the accuracy of the correlation-based and variance-based methods is questionable. The moment-independent Global IM analysis in this research provides a better understanding of how experimental uncertainty data affects industry’s plant models and where improvements in that data will have the largest benefit for improving fire modeling accuracy in causing core damage. Among the five unscreened parameters obtained from the Morris EE analysis, the Global IM analysis results for the case study indicated that max heat release rate and fire location are the most important parameters. The report also outlines benefits of using a unified computational platform that integrates the underlying simulations (e.g., a fire progression model), quantitative screening (using the Morris EE method), and the Global IM analysis. A unified platform can (i) facilitate the ranking of input parameters considering multiple key fire scenarios simultaneously, rather than considering one scenario at a time, (ii) contribute to more explicit and accurate treatment of dependencies at multiple levels of Fire PRA, (iii) facilitate the sampling-based uncertainty quantification for Fire PRA, and (iv) help generating both “industry-wide” and "plant-specific" ranking of uncertainty sources in Fire PRA. Future research should be done to include additional parameters such as detection/suppression or cable fire spread. Adding a PRA software such as SAPHIRE to the RAVEN platform would help with plant model integration and improve treatment of fire-induced dependency. The I-PRA risk importance ranking methodology offered in this report can provide valuable information for efficiently (a) enhancing the realism of Fire PRA for existing plants and (b) supporting the development of Dynamic Fire PRA for advanced reactors and new plants.

97 MATHEMATICS AND COMPUTING↗

Multiple parameter dynamic photoresponse microscopy for data-intensive optoelectronic measurements of van der Waals heterostructures

Quantum devices made from van der Waals (vdW) heterostructures of two dimensional (2D) materials may herald a new frontier in designer materials that exhibit novel electronic properties and unusual electronic phases. However, due to the complexity of layered atomic structures and the physics that emerges, experimental realization of devices with tailored physical properties will require comprehensive measurements across a large domain of material and device parameters. Such multi-parameter measurements require new strategies that combine data-intensive techniques—often applied in astronomy and high energy physics—with the experimental tools of solid state physics and materials science. We discuss the challenges of comprehensive experimental science and present a technique, called Multi-Parameter Dynamic Photoresponse Microscopy (MPDPM), which utilizes ultrafast lasers, diffraction limited scanning beam optics, and hardware automation to characterize the photoresponse of 2D heterostructures in a time efficient manner. Using comprehensive methods on vdW heterostructures results in large and complicated data sets; in the case of MPDPM, we measure a large set of images requiring advanced image analysis to extract the underlying physics. We discuss how to approach such data sets in general and in the specific case of a graphene–boron nitride–graphite heterostructure photocell.

36 MATERIALS SCIENCE↗

Evaluation of AI-enabled Digital Documented Safety Analysis

The National Reactor Innovation Center (NRIC) is leading a transformative initiative to accelerate advanced reactor deployment by fundamentally reimagining how nuclear safety basis documentation is developed, reviewed, and maintained. Traditional Documented Safety Analysis (DSA) processes for DOE-authorized facilities rely on static, document-centric workflows that consume significant time and resources, exemplified by recent major licensing efforts requiring hundreds of thousands of staff hours and millions of pages of documentation review. These conventional approaches create barriers to the rapid, cost-effective deployment of advanced reactors that America's future energy needs demand. NRIC's DOE Authorization Digital Transformation Project addresses these challenges through an innovative framework that integrates artificial intelligence (AI), digital engineering, and systems-based data management into a cohesive digital ecosystem. This white paper presents NRIC's methodology for evaluating AI-enabled document generation capabilities within this broader digital infrastructure, using the Demonstration of Microreactor Experiments (DOME) facility as a pilot case study. The evaluation will assess an AI tool's ability to generate a Preliminary Documented Safety Analysis (PDSA) through progressive integration stages—from standalone document processing to full digital thread connectivity—while maintaining rigorous verification, validation, and regulatory acceptance standards. By establishing dynamic, traceable connections between design data and safety documentation, NRIC's approach has the potential to reduce both document development time and regulatory review cycles by as much as 50%, while simultaneously improving accuracy, consistency, and traceability. This initiative represents a critical step toward establishing reusable digital infrastructure that reactor developers can leverage to accelerate their path from concept to commercial operation, directly supporting NRIC's mission to demonstrate and deploy advanced nuclear energy technologies.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Bryan Mound InSAR Analysis U.S. Strategic petroleum Reserve.

The U.S. Strategic Petroleum Reserve (SPR) is a stockpile of emergency crude oil to be tapped into if a disruption in the nation's oil supply occurs. The SPR is comprised of four salt dome sites. Subsidence surveys have been conducted either annually or biennially at all four sites over the life of the program. Monitoring of surface behavior is a first line defense to detecting possible subsurface cavern integrity issues. Over the life of the Bryan Mound site, subsidence rates over abandoned Cavern 3 have continuously been the highest at the site. In an effort to try and understand the subsurface dynamics, specifically over Bryan Mound Cavern 3, historic interferometric synthetic aperture radar (InSAR) data was acquired and processed by TRE Altamira. InSAR involves the processing of multiple satellite synthetic aperture radar scenes acquired across the same location of the Earth's surface at different times to map surface deformation. The analysis of the data has the ability to detect millimeters of motion spanning days, months, year and decades, across specific sites. The intent in regards to the Bryan Mound site was (1) to confirm the higher subsidence rates recorded over abandoned Cavern 3 indicated by land survey and (2) understand the regional surface behavior. This report describes the InSAR analysis results, how those results compare to the historical collection of land survey data, and what additional information the data has provided towards understanding the response recorded at the surface.

02 PETROLEUM↗

Visual HPC Workflows for the Analysis of System Dynamics Models

Visual analytics supported by high performance computing (HPC) accelerates and enhances the discovery, exploration, and analysis of causal patterns in complex system dynamics (SD) models. We present a suite of visualization-assisted ensemble-based techniques for hypothesis generation and testing, and for sensitivity analysis. By employing HPC to provide parallel, on-demand simulation of SD models, one can “steer” an ensemble of simulated scenarios in real time as one first formulates and then informally tests those hypotheses: this provides rapid feedback for analysts to refine their understanding of the causal relationships emergent from a model. Such understandings can be followed and augmented by rigorous application of statistical methods, namely global variance-based sensitivity analysis, Monte-Carlo filtering, adaptive regional sensitivity analysis, and self-organized maps: here timely computation relies on HPC, while effective presentation emphasizes high-dimensional multivariate data visualization. Immersive visualization in virtual 3D environments provides an excellent adjunct to the traditional 2D graphics typically used for SD models, as it generates an embodied understanding of model behavior and facilitates an active, collaborative critique of model structure and output. Finally, we summarize prospects for HPC-enabled visual analytics applied to SD modeling.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Physical Origins of the Transient Absorption Spectra and Dynamics in Thin-Film Semiconductors: The Case of BiVO 4

Transient absorption (TA) spectroscopy is uniquely suited for understanding kinetic processes initiated by light over vast ranges of time. In combination with white light probes, the recorded differential absorption spectra can contain spectroscopic signatures characteristic of specific charge carrier population densities. However, disentangling the often-complex and convoluted spectra is made challenging without robust analysis methods relating the underlying physical mechanisms to the spectral components. In this work, we address the origin of the transient spectra of a model system of emerging solar energy harvesting materials using a monoclinic BiVO 4 thin film. Using ground-state optical properties of the semiconductor, we find the main derivative-like spectral response to be related to shifting and broadening of oscillators, rather than specific carrier-related transitions. However, by using the Drude optical model of free carriers, we also identify the transient response related to free-hole density. Importantly, sample heating from the optical pumping, which begins at ~10 ps and plateaus by ~200 ps, dominates the overall spectral response at longer times. On the basis of a physical model of the spectral response, a kinetic model is developed that describes the pump power dependence of the free-hole density, as well as the temporal evolution of the spectral changes associated with shifting and broadening of oscillators. First-principles density functional theory calculations are used to rationalize experimental measurements. Here, this comprehensive approach to analyzing and modeling the TA spectra offers a generalizable basis for understanding the complex pump-probe data, reveals thermal heating artifacts that are frequently erroneously assigned to long-lived photocarriers, and offers a path to eliminating ambiguity in analysis of photocarrier dynamics in solid-state systems.

36 MATERIALS SCIENCE↗

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↗