skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Time scales of relaxation dynamics during transient conditions in two-phase flow: RELAXATION DYNAMICS

Abstract

The relaxation dynamics toward a hydrostatic equilibrium after a change in phase saturation in porous media is governed by fluid reconfiguration at the pore scale. Little is known whether a hydrostatic equilibrium in which all interfaces come to rest is ever reached and which microscopic processes govern the time scales of relaxation. Here we apply fast synchrotron-based X-ray tomography (X-ray CT) to measure the slow relaxation dynamics of fluid interfaces in a glass bead pack after fast drainage of the sample. The relaxation of interfaces triggers internal redistribution of fluids, reduces the surface energy stored in the fluid interfaces, and relaxes the contact angle toward the equilibrium value while the fluid topology remains unchanged. The equilibration of capillary pressures occurs in two stages: (i) a quick relaxation within seconds in which most of the pressure drop that built up during drainage is dissipated, a process that is to fast to be captured with fast X-ray CT, and (ii) a slow relaxation with characteristic time scales of 1–4 h which manifests itself as a spontaneous imbibition process that is well described by the Washburn equation for capillary rise in porous media. The slow relaxation implies that a hydrostatic equilibrium is hardlymore » ever attained in practice when conducting two-phase experiments in which a flux boundary condition is changed from flow to no-flow. Implications for experiments with pressure boundary conditions are discussed.« less

Authors:
ORCiD logo [1]; ORCiD logo [2];  [3]; ORCiD logo [4]; ORCiD logo [3]
  1. School of Chemical, Biological and Environmental Engineering, Oregon State University, Corvallis Oregon USA; Department Soil Physics, Helmholtz-Centre for Environmental Research-UFZ, Halle Germany
  2. Shell Global Solutions International B.V., Rijswijk Netherlands
  3. School of Chemical, Biological and Environmental Engineering, Oregon State University, Corvallis Oregon USA
  4. Department Soil Physics, Helmholtz-Centre for Environmental Research-UFZ, Halle Germany; Institut für Agrar- und Ernährungswissenschaften, Martin-Luther-Universität Halle-Wittenberg, Halle Germany
Publication Date:
Research Org.:
Argonne National Lab. (ANL), Argonne, IL (United States). Advanced Photon Source (APS)
Sponsoring Org.:
NSFOTHER
OSTI Identifier:
1372226
Resource Type:
Journal Article
Resource Relation:
Journal Name: Water Resources Research; Journal Volume: 53; Journal Issue: 6
Country of Publication:
United States
Language:
ENGLISH
Subject:
54 ENVIRONMENTAL SCIENCES

Citation Formats

Schlüter, Steffen, Berg, Steffen, Li, Tianyi, Vogel, Hans-Jörg, and Wildenschild, Dorthe. Time scales of relaxation dynamics during transient conditions in two-phase flow: RELAXATION DYNAMICS. United States: N. p., 2017. Web. doi:10.1002/2016WR019815.
Schlüter, Steffen, Berg, Steffen, Li, Tianyi, Vogel, Hans-Jörg, & Wildenschild, Dorthe. Time scales of relaxation dynamics during transient conditions in two-phase flow: RELAXATION DYNAMICS. United States. doi:10.1002/2016WR019815.
Schlüter, Steffen, Berg, Steffen, Li, Tianyi, Vogel, Hans-Jörg, and Wildenschild, Dorthe. Thu . "Time scales of relaxation dynamics during transient conditions in two-phase flow: RELAXATION DYNAMICS". United States. doi:10.1002/2016WR019815.
@article{osti_1372226,
title = {Time scales of relaxation dynamics during transient conditions in two-phase flow: RELAXATION DYNAMICS},
author = {Schlüter, Steffen and Berg, Steffen and Li, Tianyi and Vogel, Hans-Jörg and Wildenschild, Dorthe},
abstractNote = {The relaxation dynamics toward a hydrostatic equilibrium after a change in phase saturation in porous media is governed by fluid reconfiguration at the pore scale. Little is known whether a hydrostatic equilibrium in which all interfaces come to rest is ever reached and which microscopic processes govern the time scales of relaxation. Here we apply fast synchrotron-based X-ray tomography (X-ray CT) to measure the slow relaxation dynamics of fluid interfaces in a glass bead pack after fast drainage of the sample. The relaxation of interfaces triggers internal redistribution of fluids, reduces the surface energy stored in the fluid interfaces, and relaxes the contact angle toward the equilibrium value while the fluid topology remains unchanged. The equilibration of capillary pressures occurs in two stages: (i) a quick relaxation within seconds in which most of the pressure drop that built up during drainage is dissipated, a process that is to fast to be captured with fast X-ray CT, and (ii) a slow relaxation with characteristic time scales of 1–4 h which manifests itself as a spontaneous imbibition process that is well described by the Washburn equation for capillary rise in porous media. The slow relaxation implies that a hydrostatic equilibrium is hardly ever attained in practice when conducting two-phase experiments in which a flux boundary condition is changed from flow to no-flow. Implications for experiments with pressure boundary conditions are discussed.},
doi = {10.1002/2016WR019815},
journal = {Water Resources Research},
number = 6,
volume = 53,
place = {United States},
year = {Thu Jun 01 00:00:00 EDT 2017},
month = {Thu Jun 01 00:00:00 EDT 2017}
}
  • We analyze a three-variable rate equation model that takes into account carrier capture and Pauli blocking in quantum dot semiconductor lasers. The exponential decay of the relaxation oscillations is analyzed from the linearized equations in terms of three key parameters that control the time scales of the laser. Depending on their relative values, we determine two distinct two-variable reductions of the rate equations in the limit of large capture rates. The first case leads to the rate equations for quantum well lasers, exhibiting relaxation oscillations dynamics. The second case corresponds to dots nearly saturated by the carriers and is characterizedmore » by the absence of relaxation oscillations.« less
  • The development of efficient theoretical methods for describing electron transfer (ET) reactions in condensed phases is important for a variety of chemical and biological applications. Previously, dynamical dielectric continuum theory was used to derive Langevin equations for a single collective solvent coordinate describing ET in a polar solvent. In this theory, the parameters are directly related to the physical properties of the system and can be determined from experimental data or explicit molecular dynamics simulations. Herein, we combine these Langevin equations with surface hopping nonadiabatic dynamics methods to calculate the rate constants for thermal ET reactions in polar solvents formore » a wide range of electronic couplings and reaction free energies. Comparison of explicit and implicit solvent calculations illustrates that the mapping from explicit to implicit solvent models is valid even for solvents exhibiting complex relaxation behavior with multiple relaxation time scales and a short-time inertial response. The rate constants calculated for implicit solvent models with a single solvent relaxation time scale corresponding to water, acetonitrile, and methanol agree well with analytical theories in the Golden rule and solvent-controlled regimes, as well as in the intermediate regime. The implicit solvent models with two relaxation time scales are in qualitative agreement with the analytical theories but quantitatively overestimate the rate constants compared to these theories. Analysis of these simulations elucidates the importance of multiple relaxation time scales and the inertial component of the solvent response, as well as potential shortcomings of the analytical theories based on single time scale solvent relaxation models. This implicit solvent approach will enable the simulation of a wide range of ET reactions via the stochastic dynamics of a single collective solvent coordinate with parameters that are relevant to experimentally accessible systems.« less
  • A method for controlling the thermal boundary conditions of non-equilibrium molecular dynamics simulations is presented. The method is simple to implement into a conventional molecular dynamics code and independent of the atomistic model employed. It works by regulating the temperature in a thermostatted boundary region by feedback control to achieve the desired temperature at the edge of an inner region where the true atomistic dynamics are retained. This is necessary to avoid intrinsic boundary effects in non-equilibrium molecular dynamics simulations. Three thermostats are investigated: the global deterministic Nose-Hoover thermostat and two local stochastic thermostats, Langevin and stadium damping. The lattermore » thermostat is introduced to avoid the adverse reflection of phonons that occurs at an abrupt interface. The method is then extended to allow atomistic/continuum models to be thermally coupled concurrently for the analysis of large steady state and transient heat conduction problems. The effectiveness of the algorithm is demonstrated for the example of heat flow down a three-dimensional atomistic rod of uniform cross-section subjected to a variety of boundary conditions.« less
  • The analysis of the field equations and the conservation laws is extended into the skin boundary region which maintains the initially static configuration of the two fluid spheres. Although this leads to a correction of the gravitational-radiation energy-loss rate from a dependence of Eapprox...cap alpha../sup -4/rho/sub 0/ /sup -4/ to ..cap alpha../sup -6/rho/sub 0/ /sup -2/, the importance of the nonlinear structure-dependent terms remain and the essential conclusion, that the quadrupole formula does not apply to this problem, is unaltered. The hydrodynamic, stress-breaking, and free-fall time scales are considered. It is shown that insofar as the quadrupole-formula comparison for free-fallmore » is concerned, only the contribution from bulk motion of the fluid spheres need be considered since tidal quadrupole deformation contributes negligibly to the quadrupole formula. With reference to our problem, it is shown that a recent derivation of the quadrupole formula for free-fall is incorrect and it is suggested that certain other derivations may be applicable only to the radiation damping of a single body.« less