scholarly journals A model for discrete fracture-clay rock interaction incorporating electrostatic effects on transport

Author(s):  
Carl I. Steefel ◽  
Christophe Tournassat

AbstractA model based on the code CrunchClay is presented for a fracture-clay matrix system that takes electrostatic effects on transport into account. The electrostatic effects on transport include those associated with the development of a diffusion potential as captured by the Nernst-Planck equation, and the formation of a diffuse layer bordering negatively charged clay particles within which partial anion exclusion occurs. The model is based on a dual continuum formulation that accounts for diffuse layer and bulk water pore space, providing a more flexible framework than is found in the classical mean electrostatic potential models. The diffuse layer model is obtained by volume averaging ion concentrations in the Poisson-Boltzmann equation, but also includes the treatment of longitudinal transport within this continuum. The calculation of transport within the bulk and diffuse layer porosity is based on a new formulation for the Nernst-Planck equation that considers averaging of diffusion coefficients and accumulation factors at grid cell interfaces. Equations for function residuals and the associated Jacobian matrix are presented such that the system of nonlinear differential-algebraic equations can be solved with Newton’s method. As an example, we consider a 2D system with a single discrete fracture within which flow and advective transport occurs that is coupled to diffusion in the clay-rich matrix. The simulation results demonstrate the lack of retardation for anions (e.g., 36Cl−) of the contaminant plume within the fracture flow system because they are largely excluded from the charged clay rock, while the migration of cations (e.g., 90Sr++) is more strongly attenuated. The diffusive loss of divalent cations in particular from the fracture is accentuated by their accumulation in the diffuse layer within the clay-rich matrix.

2020 ◽  
Author(s):  
Gilles Mpembele ◽  
Jonathan Kimball

<div>The analysis of power system dynamics is usually conducted using traditional models based on the standard nonlinear differential algebraic equations (DAEs). In general, solutions to these equations can be obtained using numerical methods such as the Monte Carlo simulations. The use of methods based on the Stochastic Hybrid System (SHS) framework for power systems subject to stochastic behavior is relatively new. These methods have been successfully applied to power systems subjected to</div><div>stochastic inputs. This study discusses a class of SHSs referred to as Markov Jump Linear Systems (MJLSs), in which the entire dynamic system is jumping between distinct operating points, with different local small-signal dynamics. The numerical application is based on the analysis of the IEEE 37-bus power system switching between grid-tied and standalone operating modes. The Ordinary Differential Equations (ODEs) representing the evolution of the conditional moments are derived and a matrix representation of the system is developed. Results are compared to the averaged Monte Carlo simulation. The MJLS approach was found to have a key advantage of being far less computational expensive.</div>


Author(s):  
Achim Ilchmann ◽  
Jonas Kirchhoff

AbstractWe investigate genericity of various controllability and stabilizability concepts of linear, time-invariant differential-algebraic systems. Based on well-known algebraic characterizations of these concepts (see the survey article by Berger and Reis (in: Ilchmann A, Reis T (eds) Surveys in differential-algebraic equations I, Differential-Algebraic Equations Forum, Springer, Berlin, pp 1–61. 10.1007/978-3-642-34928-7_1)), we use tools from algebraic geometry to characterize genericity of controllability and stabilizability in terms of matrix formats.


Author(s):  
J. P. Meijaard ◽  
V. van der Wijk

Some thoughts about different ways of formulating the equations of motion of a four-bar mechanism are communicated. Four analytic methods to derive the equations of motion are compared. In the first method, Lagrange’s equations in the traditional form are used, and in a second method, the principle of virtual work is used, which leads to equivalent equations. In the third method, the loop is opened, principal points and a principal vector linkage are introduced, and the equations are formulated in terms of these principal vectors, which leads, with the introduced reaction forces, to a system of differential-algebraic equations. In the fourth method, equivalent masses are introduced, which leads to a simpler system of principal points and principal vectors. By considering the links as pseudorigid bodies that can have a uniform planar dilatation, a compact form of the equations of motion is obtained. The conditions for dynamic force balance become almost trivial. Also the equations for the resulting reaction moment are considered for all four methods.


Sign in / Sign up

Export Citation Format

Share Document