Modeling bimolecular reactions and transport in porous media
|
|
- Jasper Bell
- 5 years ago
- Views:
Transcription
1 GEOPHYSICAL RESEARCH LETTERS, VOL. 36, L02407, doi: /2008gl036381, 2009 Modeling bimolecular reactions and transport in porous media Yaniv Edery, 1 Harvey Scher, 1 and Brian Berkowitz 1 Received 21 October 2008; revised 8 December 2008; accepted 18 December 2008; published 28 January [1] We quantitatively account for the measured concentration profile from a laboratory experiment of bimolecular (A + B! C) reactive transport in a porous medium with a particle tracking (PT) model. The PT results are in contrast to the analytical solution of the continuum scale advectiondispersion-reaction equation, which results in an excess quantity of reaction product (C). The approaches differ in the treatment of the mixing zone, the fluctuations due to the low reactant concentrations, and the localized nature of the reaction. The PT can accommodate a range of transport modes with different temporal distributions. An exponential temporal distribution is equivalent to Fickian transport, which we use for the comparison to the laboratory data; a truncated power-law (TPL) temporal distribution yields a non-fickian transport characteristic of heterogeneous media. We study the influence of disorder on the mixing zone and the product concentration profiles via these contrasting transport modes. Citation: Edery, Y., H. Scher, and B. Berkowitz (2009), Modeling bimolecular reactions and transport in porous media, Geophys. Res. Lett., 36, L02407, doi: /2008gl Introduction [2] Considerable effort has been invested over the last three decades in conceptualizing and modeling reactive transport in porous media. The problem is of key importance for understanding the evolution of different chemical contaminants in groundwater. In porous media, the reactions are localized to pore-scale sites and heterogeneity can affect the nature of the transport of reactants. The high degree of fluctuations in pore-scale mixing, especially for low concentrations, can impact the larger-scale reactive transport behavior. The actual amount of mixing between reactants (and hence the amount of product) may be less than would be otherwise predicted by continuum models, due to averaging the concentrations and the mixing zone [e.g., Rashidi et al., 1996; Cao and Kitanidis, 1998; Raje and Kapoor, 2000]. Palanichamy et al. [2007] discussed the inappropriateness of assumptions made in a continuum picture of low concentrations of reactants undergoing random collisions. Methods to account for non-uniform mixing and random reaction of the reactants, in both space and time, therefore appear necessary. The objective of this study is to develop such a method, and to compare the results to available data from a laboratory experiment of bimolecular reactive transport in a porous medium and to solutions of continuum models. 1 Department of Environmental Sciences and Energy Research, Weizmann Institute of Science, Rehovot, Israel. Copyright 2009 by the American Geophysical Union /09/2008GL [3] We use a particle tracking (PT) technique to simulate the subtle features of reactions and different transport modes in porous media. The PT approach has a long history, having been first described by Smoluchowski about a century ago to treat Brownian motion. The extension of PT to account for reactions, by specifying probabilistic rules of particle interaction, has been discussed extensively [e.g., Gillespie, 1976, 1977; Fabriol et al., 1993; Sun, 1999; Cao et al., 2005; Lindenberg and Romero, 2007; Palanichamy et al., 2007; Srinivasan et al., 2007; Yuste et al., 2008]. [4] We have used a PT approach previously to demonstrate non-fickian transport in a random fracture network, modeling it within the continuous time random walk (CTRW) framework [Berkowitz and Scher, 1998]. The results of direct CTRW-PT simulations are in excellent agreement with analytical and numerical solutions of the CTRW equations [Dentz et al., 2004]. The CTRW-PT approach allows considerable flexibility in the choice of probability density functions (pdf s) governing the spatial and temporal aspects of the particle motion [Dentz et al., 2008]. The limiting choice of an exponential temporal pdf with a spatial pdf having finite moments (e.g., an exponential pdf) conveniently reproduces the solutions of the advective-dispersion equation (ADE). Changing the temporal pdf to a truncated power law (TPL) allows treatment of a full range of non-fickian behaviors [Berkowitz et al., 2006; Dentz et al., 2008]. [5] Here we focus on the bimolecular chemical reaction A + B! C, introducing into the CTRW-PT framework interactions among migrating particles. We maintain the same degree of flexibility for reaction (i.e., local equilibrium or kinetic controlled) as for transport. For our present purposes we simplify a general picture of the reacting species; we specify that the reaction time is fast compared to the transport time and focus on the type of transport and on the spatial feature of the bimolecular process, which we use to capture the local nature of the interaction in a porous medium. The radius of interaction R is a key parameter in defining the mixing zone in this picture; it is the upper limit of the volume within which A and B react. The variation of R from local interaction to uniform mixing presents an opportunity to compare the PT results to the spatial distribution of the reactant product (i.e., C) in laboratory experiments and to the analytical solution of the advection-dispersion-reaction equation (ADRE) discussed below. Few detailed comparisons of PT for reactive transport, laboratory experiments and the ADRE have been undertaken to date. The ADRE is limited to Fickian transport. Depending on the degree of porous medium heterogeneity, we consider the influence of both Fickian and non-fickian transport as well as the stochastic nature of the bimolecular chemical reaction on the effective mixing zone. We stress that consideration of PT simulations in two and three dimensions (2d, 3d) is critical. In contrast to real systems, random walks in 1d guarantee L of5
2 that all A and B particles will interact, generally over short length scales, giving rise to a distorted mixing zone. [6] A prime example of the continuum approach mentioned above is the ADRE [e.g., Raje and Kapoor, 2000; Gramling et al., 2002], which for the bimolecular reaction and in one dimension is written i /@t + v@c i i = r(c A, c B )wherec i = c i (x, t) (i = A, B) isthe concentration, v is the velocity, D is the dispersion coefficient and jr(c A, c B )j is the total rate of product creation via reaction. A common form is r(c A, c B )= Gc A c B ; this coupled system involves an equation for each reactant and one for the product (c C ) with the source Gc A c B. Many studies are based on analytical and numerical solutions of the ADRE [e.g., Najafi and Hajinezhad, 2008; Rubio et al., 2008]. The ADRE assumes that dispersion is Fickian and that the reactants undergo continuous mixing along the flow path. As such, the reactants are assumed to be well-mixed so that deterministic rate laws based on local concentrations of reacting species are employed. While Eulerian and Lagrangian descriptions of (reactive) transport should in principle yield identical results [e.g., Gardiner, 2004], the ADRE description may in the case of reactive transport require high spatial resolution to fully capture, if at all, the local-scale heterogeneous (non-uniform) reactant and product concentration distributions. An informative critique of this approach is given by Palanichamy et al. [2007]. We stress that in contrast to the continuum (ADRE or other) framework, the PT approach monitors particle positions in both space and time, and accounts naturally for both particle spreading (characterized by variance) and mixing (i.e., relative positions between particles, characterized by covariance). [7] As an example of bimolecular reactive transport in porous media, we consider the laboratory column experiment described by Gramling et al. [2002]. Colorimetric reactions between CuSO 4 and EDTA 4 were measured in a 1d flow field, in a cell (36 cm (length) 5.5 cm (width) 1.8 cm (depth)) packed with cryolite. The flow cell was saturated with EDTA 4 (denoted as B) and CuSO 4 (denoted as A) was subsequently injected as a step input. The mixing and reaction zone of the product CuEDTA 2 was measured and compared to an analytical solution of the ADRE. An earlier, similar experimental study was reported by Raje and Kapoor [2000]. 2. Methods [8] We consider a general CTRW-PT model that incorporates bimolecular reactions. The model involves numerical simulation via application of the equations of motion [Dentz et al., 2004, 2008], s (N+1) = s (N) + V (N), t (N+1) = t (N) +t (N) where (s (N), t (N) ) denotes the location of a particle in space-time after N steps. The spatial and temporal random increments V (N) and t (N) are distributed according to the joint transition displacement and time probability density y(s, t). These equations comprise a full range of CTRWs, including, e.g., general coupled space-time behavior and correlation with previous steps. A particular feature is that this model accounts for the time difference between the particle step movements and their typical reaction times. [9] Here, we use a decoupled form of the random walk, y(s, t) = F(s)y(t), wherein the spatial s, angle f and time steps are sampled from different, independent distributions, p(s), W(f) and y(t), respectively (see Dentz et al. [2008] for detailed discussions of coupled and decoupled formulations and the angular dependence W(f)). To obtain an ADE from these distributions, we assign decaying exponential distributions, exp( l s t) and exp( l t t), to p(s) and y(t), respectively. R We derive l s andr l t using the definitions (for 2d) v y = R ps ðþs 2 ds yðþtdt t and D y = R ps ðþs 3 ds, where v y and D y are the yðþtdt t generalized particle velocity and dispersion, respectively [Berkowitz et al., 2006]. In simulations of non-fickian transport, we choose a truncated power law (TPL) for the time distribution and an exponential distribution for the spatial transition. The TPL is characterized by three parameters: the power law exponent b, the characteristic transition time t 1, and the cut-off time to Fickian transport, t 2 (see, e.g., Berkowitz et al. [2006] for details). [10] To simulate bimolecular reactions, two types of particles, marked A and B are introduced into the domain. We record the position of each particle at specified times, as it migrates through the domain. A sampling time step, Dt, is specified, at which we calculate the distance between a given A particle and a B particle. If this distance is smaller than a prescribed effective reaction radius, R, then A and B are replaced by a C particle. The C particles are allowed to migrate through space by the same rules as the A and B particles, although if desired, they can be held in place ( precipitated ), or the reaction can be treated as reversible. A further modification can be introduced to effect different reaction rates, i.e., a probability to produce a C particle even if the A, B separation is R. A simple way to effect this stochastic variation in the rate is to replace the fixed R by a normal distribution of R, N(R m, s), where R m is the mean radius and s is the standard deviation. The reaction rate features of the model are quite flexible and other investigations are in progress (e.g., by relating R directly to specific reaction rate coefficients). Clearly, the choice of R will affect the degree of particle interaction and production of C particles; comparing simulations to experiments can in fact be used to estimate the effective reaction radius in a physical system (see discussion below). [11] The simulation procedure is numerically intensive, but not impractical; calculations with a C compiler, using for example 50,000 particles each of A and B, required only one hour on a Pentium 4, 2.4 GHz processor with 1 GB of RAM. With 50,000 particles, we find that discrete artifacts in the spatial distribution are virtually eliminated; the difference in peak value of the spatial distribution of C particles between 50,000 particles and 100,000 particles is less than 1%. [12] We focus on an experiment of Gramling et al. [2002]. To compare the CTRW-PT model against this experiment and the ADRE solution, the initial conditions of the experiment require uniformly distributed B particles in the bounded flow domain. From this spatial distribution, the concentration of particles per unit length (c 0 = number of particles/length of flow cell (36 cm)) is used to define the influx of A particles distributed uniformly along the inlet face. Specifically, the inflow of A particles per time step Dt is set equal to the number of particles in a length of flow cell Dx/v. Thus, the inflow spatial concentration of A particles 2of5
3 velocity. The dispersion coefficient, D, was estimated p by fitting the solution c A (x, t)/c 0 = 0.5 erfc [(x vt)/(2 ffiffiffiffiffi Dt )] to the PT simulations, which was run for a constant influx of A particles. This is an approximate solution valid for D vx in a flow domain with infinite boundaries. Note that because in general v y 6¼ v and D y 6¼ D, small variations between the results can arise. These values of v y and D y were then used in the analytical solution of the ADRE for the distribution of C particles p [Gramling et al., 2002]: c C (x, t)/c 0 = 0.5 erfc [(vt x)/(2 ffiffiffiffiffi p ffiffiffiffiffi Dt )], for x < vt; 0.5 erfc [(x vt)/(2 Dt )], for x > vt. Figure 1. Spatial distribution of particles A, B, and C after 202 time steps, with R = 0.5 cm, using the PT-ADE form, with Dt = 0.1 s; 50,000 particles each of A and B. equals the outflow spatial concentration of B particles, with C particles being formed in the mixing zone. [13] We set v y = 0.67 cm/s and D y = cm 2 /s, to match the fluid velocity and dispersion values of a specific experiment of Gramling et al. [2002], so that l s = 5.74 cm 1 and l t =1.92s 1 for the decaying exponential forms of p(s) and y(t), respectively. The angle f representing the direction of each spatial transition was chosen from a normal N(0, 1) distribution between ( p 13, p 13 ); each particle that hit the boundary of the domain was advanced at the same angle of transition. With this angular range, the PT simulations using the prescribed values of l s and l t yielded concentration profiles with approximately the same values of v and D chosen above. Employing the ADE form of our PT model with only the uniformly distributed A particles, we confirmed that the calculated x position of the particle center of mass at different times matched the center of mass 3. Results and Discussion [14] Results are shown for particle profiles c/c 0, where c denotes the PT simulation of c A,B,C, the analytical solution of the ADRE and the measured profile [Gramling et al., 2002]. A feature of the experiment is the direct measurement of a profile rather than only breakthrough curves. There are three basic comparisons among these various c: the PT c A,B,C in the Fickian (ADE) mode (Figure 1), with the experimental data and the analytical ADRE c C (Figure 2) and with the PT-TPL (non-fickian) mode (Figure 3). [15] Figure 1 shows the spatial distribution of particles A, B, and C, with R = 0.5 cm, using the ADE form. The A particles show the step inflow (c A /c 0 = 1) and then a sharp decline and forward tail. The B particles, originally distributed uniformly, leave a backward tail as they advance with the bulk fluid flow. Between these two profiles, C particles are produced in the mixing zone; c C is distributed approximately symmetrically between the forward and backward tails of the A and B particles. Increasing R leads to an increased number of C particles. The time step Dt was determined to be the value where further decrease caused no change in the production of C particles. [16] Figure 2 shows PT-ADE mode simulations, compared to the ADRE analytical solution and to the experi- Figure 2. Simulations (smoothed) with the PT-ADE, compared to the ADRE analytical solution (dotted line) and the experimental measurements (points) of Gramling et al. [2002], showing the relative concentration of C particles at time s (Dt = 0.01 s); dot-dashed line is the PT-ADE with fixed R = 0.5 cm; solid line is the PT-ADE with N(0.5, 1) values of R sampled at each time step. Inset shows experimental measurements and the PT-ADE simulation with N(0.5, 1) distributed R, allowing for a 1% deviation in the velocity in the experiment of Gramling et al. [2002]. 3of5
4 Figure 3. Spatial distribution of particles A, B, and C after 202 time steps, with R normally distributed N(0.5, 1) sampled at each time step, using the PT-TPL form, with Dt =0.1s,b = 1.6, t 1 =2L/v (see text), t 2 =10 5 s, 50,000 particles each of A and B. mental measurements (points) of Gramling et al. [2002]. The simulations are given for two cases, with a fixed R = 0.5 cm, and with normally distributed N(0.5, 1) values of R sampled at each time step. Gramling et al. [2002] showed that the analytical solution of the ADRE significantly overpredicts the experimentally measured reaction product by up to 40% of the peak value, as seen in Figure 2. The experiments of Raje and Kapoor [2000] also found that while the analytical solution of the ADE without a reactive component quantifies transport of each reactant separately, the analytical solution of the ADRE over-predicts the amount of reaction product by this amount. Both studies suggested that this error is due to the high degree of mixing at the local scale assumed in the ADRE formulation. [17] On the other hand the resemblance of the C particle concentration profiles between the PT simulations (Figure 2) and the measurements of Gramling et al. [2002] is striking. The value R = 0.5 cm can be considered representative of the mixing volume in which the reaction takes place; this value seems reasonable given the average grain size of 0.13 cm, and the relatively high fluid velocity of 0.67 cm/s. Comparing the PT simulations for the fixed R and a variable R (N(0.5,1) distributed) at each time step, it is seen that the latter condition leads to a slightly broader concentration profile. Moreover, an R can be defined above which the particles appear well mixed, i.e., the peak of the C profile increases, so that the ADRE can be recovered. The inset in Figure 2 shows a comparison of the fit of the PT-ADE against the measurements, allowing for a small degree of experimental error. Clearly, even a 1% deviation in the estimated inlet fluid velocity or initial time at which A particles enter the system not unreasonable given the practical limitations of controlling the change in fluid inflows, and the resolution of images used to estimate c C modifies the location of the predicted peak concentration and of the relative fit to the tails of the concentration profile. These results clearly contrast the ADRE assumption of uniform mixing and the particle tracking accounting of local-scale heterogeneity and nonuniform mixing, and support the arguments against the use of the ADRE given in the Introduction. In other words, the ADRE picture of mixing does not hold because the degree of reaction is strongly dependent on the local-scale spatial distribution of the reactants. [18] While the column experiment of Gramling et al. [2002] is well-matched by the ADE form of the CTRW-PT model, asymmetric spatial particle plumes that are distinct from Fickian patterns (whether in the forward and/or backward directions) are often present in both homogeneous and heterogeneous porous media [Berkowitz et al., 2006]. In such cases, transport is often better described within the CTRW framework [Berkowitz et al., 2006] using a temporal distribution with a significant power-law region. Particle tracking simulations with the TPL form were therefore considered, assuming a moderate degree of non- Fickian transport (b = 1.6); this value is typical of transport behavior in many columns [e.g., Berkowitz et al., 2006], and leads to spatial concentration profiles with heavier backward tails. [19] Figure 3 shows the spatial distribution of particles A, B, and C, with normally distributed N(0.5, 1) values of R sampled at each time step, using the TPL form (b = 1.6, t 1 =2L/v, where L = 0.13 cm (the grain size) and v = 0.67 cm/s, t 2 =10 5 s); for comparison, the corresponding C profile using the PT-ADE form is also shown. With the TPL, the backward tail of B particles is heavier than the forward tail of A particles, and as a result, the concentration profile of reaction product C particles is highly asymmetrical. In general, the spatial distribution of C particles is broader than that for the ADE form. Note that the relative concentration of A reaches a value greater than 1, which is related to the integral effect of the inflows for a TPL pdf (which will be investigated in a further study). 4. Concluding Remarks [20] The reactive transport problem considered here involves a separation between two different length scales. One is the length scale of the flow cell, at which on average, A and B particles may appear well mixed. The second length scale is of the pores where the mixing particles react. As a consequence, there is less reaction than might be expected at the large scale, such as from the analytical solution of the ADRE which assumes a greater degree of mixing. [21] The CTRW-PT model (i.e., ADE-PT, TPL-PT) enables a general representation of reactive transport, which is not bound by the assumptions of ADRE-like models that require the concentrations to be continuous and deterministic, and the transport to be Fickian. Though the CTRW-PT model does not account explicitly for physical pores in the medium, it includes the effects of local-scale dynamics via the pdf s applied to describe the particle transport. The effective distance over which particles mix and react in any specific experiment is given by the effective reaction distance R; increasing R has the effect of increasing the mixing between A and B particles, and thus of increasing the amount of product C. [22] The analysis here leads to two principal conclusions. First, the use of partial differential equation approaches to quantify reactive transport must be called into question, because of the difficulty in capturing local scale mixing and reaction. Particle tracking approaches thus appear more attractive, in spite of their potential computational limita- 4of5
5 tions for large-scale domains. Second, and related to the use of a PT approach, reactions require that reactant particles be found in close proximity in both space and time. [23] Acknowledgments. We thank the Israel Ministry of Science and Technology for financial support. References Berkowitz, B., and H. Scher (1998), Theory of anomalous chemical transport in fracture networks, Phys. Rev. E, 57, Berkowitz, B., A. Cortis, M. Dentz, and H. Scher (2006), Modeling non- Fickian transport in geological formations as a continuous time random walk, Rev. Geophys., 44, RG2003, doi: /2005rg Cao, J., and P. K. Kitanidis (1998), Pore-scale dilution of conservative solutes: An example, Water Resour. Res., 8, Cao, Y., L. Gillespie, and D. Petzold (2005), Multiscale stochastic simulation algorithm with stochastic partial equilibrium assumption for chemically reacting systems, J. Comput. Phys., 206, Dentz, M., A. Cortis, H. Scher, and B. Berkowitz (2004), Time behavior of solute transport in heterogeneous media: Transition from anomalous to normal transport, Adv. Water Resour., 27, , doi: /j. advwatres Dentz, M., H. Scher, D. Holder, and B. Berkowitz (2008), Transport behavior of coupled continuous-time random walks, Phys. Rev. E, 78, , doi: /physreve Fabriol, R., J. P. Sauty, and G. Ouzounian (1993), Coupling geochemistry with a particle tracking transport model, J. Contam. Hydrol., 13, Gardiner, C. W. (2004), Handbook of Stochastic Methods for Physics, Chemistry and the Natural Sciences, Springer Ser. Synergetics, vol. 13, 3rd ed., Springer, New York. Gillespie, D. (1976), A general method for numerically simulating the stochastic time evolution of coupled chemical reactions, J. Comput. Phys., 22, Gillespie, D. (1977), Exact stochastic simulation of coupled chemical reactions, J. Phys. Chem., 25, Gramling, C. M., C. F. Harvey, and L. C. Meigs (2002), Reactive transport in porous media: A comparison of model prediction with laboratory visualization, Environ. Sci. Technol., 36, Lindenberg, K., and A. H. Romero (2007), Numerical study of A + A! 0 and A + B! 0 reactions with inertia, J. Chem. Phys, 127(17), , doi: / Najafi, H. S., and H. Hajinezhad (2008), Solving one-dimensional advection-dispersion with reaction using some finite-diference methods, Appl. Math. Sci., 2, Palanichamy, J., T. Becker, M. Spiller, J. Köngeter, and S. Mohan (2007), Multicomponent reaction modelling using a stochastic algorithm, Comput. Visual. Sci., doi: /s00, y. Raje, D. S., and V. Kapoor (2000), Experimental study of bimolecular reaction kinetics in porous media, Environ. Sci. Technol., 34, Rashidi, M., L. Peurrung, A. F. B. Tompson, and T. J. Kulp (1996), Experimental analysis of pore-scale flow and transport in porous media, Adv. Water Resour., 3, Rubio, A. D., A. Zalts, and C. D. El Hasi (2008), Numerical solution of the advection-reaction-diffusion equation at different scales, Environ. Model. Software, 1, Srinivasan, G., D. M. Tartakovsky, B. A. Robinson, and A. B. Aceves (2007), Quantification of uncertainty in geochemical reactions, Water Resour. Res., 43, W12415, doi: /2007wr Sun, N. Z. (1999), A finite cell method for simulating the mass transport process in porous media, Water Resour. Res., 35, Yuste, S. B., J. Klafter, and K. Lindenberg (2008), Number of distinct sites visited by a subdiffusive random walker, Phys. Rev. E, 77, , doi: /physreve B. Berkowitz, Y. Edery, and H. Scher, Department of Environmental Sciences and Energy Research, Weizmann Institute of Science, Rehovot Israel. (brian.berkowitz@weizmann.ac.il; yaniv.edery@weizmann. ac.il; harvey.scher@weizmann.ac.il) 5of5
LAGRANGIAN SIMULATION OF BIMOLECULAR REACTIONS
LAGRANGIAN SIMULATION OF BIMOLECULAR REACTIONS YONG ZHANG Abstract. A fully Lagrangian approach is developed to model bimolecular reactions undergoing Fickian diffusion. The basis of the Lagrangian algorithm
More informationarxiv: v1 [physics.flu-dyn] 18 Mar 2018
APS Persistent incomplete mixing in reactive flows Alexandre M. Tartakovsky 1 and David Barajas-Solano 1 1 Pacific Northwest National Laboratory, Richland, WA 99352, USA arxiv:1803.06693v1 [physics.flu-dyn]
More informationAccepted Manuscript. Modeling Bimolecular Reactions and Transport in Porous Media Via Particle Tracking
Accepted Manuscript Modeling Bimolecular Reactions and Transport in Porous Media Via Particle Tracking Dong Ding, David A. Benson, Amir Paster, Diogo Bolster PII: S0309-1708(12)00276-X DOI: http://dx.doi.org/10.1016/j.advwatres.2012.11.001
More informationQuantification of uncertainty in geochemical reactions
WATER RESOURCES RESEARCH, VOL. 43,, doi:10.1029/2007wr006003, 2007 Quantification of uncertainty in geochemical reactions Gowri Srinivasan, 1,2 Daniel M. Tartakovsky, 1,3 Bruce A. Robinson, 4 and Alejandro
More informationRAM C. ACHARYA, ALBERT J. VALOCCHI, THOMAS W. WILLINGHAM, AND CHARLES J. WERTH
PORE-SCALE SIMULATION OF DISPERSION AND REACTION ALONG A TRANSVERSE MIXING ZONE IN TWO-DIMENSIONAL HETEROGENEOUS POROUS MEDIUM COMPOSED OF RANDOMLY DISTRIBUTED CIRCULAR AND ELLIPTICAL CYLINDERS RAM C.
More informationSOLUTE TRANSPORT. Renduo Zhang. Proceedings of Fourteenth Annual American Geophysical Union: Hydrology Days. Submitted by
THE TRANSFER FUNCTON FOR SOLUTE TRANSPORT Renduo Zhang Proceedings 1994 WWRC-94-1 1 n Proceedings of Fourteenth Annual American Geophysical Union: Hydrology Days Submitted by Renduo Zhang Department of
More informationA continuous time random walk approach to transient flow in heterogeneous porous media
WATER RESOURCES RESEARCH, VOL. 42,, doi:10.1029/2006wr005227, 2006 A continuous time random walk approach to transient flow in heterogeneous porous media Andrea Cortis 1 and Christen Knudby 2 Received
More informationLecture notes for /12.586, Modeling Environmental Complexity. D. H. Rothman, MIT September 24, Anomalous diffusion
Lecture notes for 12.086/12.586, Modeling Environmental Complexity D. H. Rothman, MIT September 24, 2014 Contents 1 Anomalous diffusion 1 1.1 Beyond the central limit theorem................ 2 1.2 Large
More informationTailing of the breakthrough curve in aquifer contaminant transport: equivalent longitudinal macrodispersivity and occurrence of anomalous transport
5 GQ7: Securing Groundwater Quality in Urban and Industrial Environments (Proc. 6th International Groundwater Quality Conference held in Fremantle, Western Australia, 2 7 December 27). IAHS Publ. no. XXX,
More informationA tempered multiscaling stable model to simulate transport in regional scale fractured media
Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 37,, doi:10.1029/2010gl043609, 2010 A tempered multiscaling stable model to simulate transport in regional scale fractured media Yong Zhang,
More informationApplication of the random walk method to simulate the transport of kinetically adsorbing solutes
Groundwater Contamination (Proceedings of the Symposium held during the Third IAHS Scientific Assembly, Baltimore, MD, May 1989), IAHS Publ. no. 185, 1989 Application of the random walk method to simulate
More informationCode-to-Code Benchmarking of the PORFLOW and GoldSim Contaminant Transport Models using a Simple 1-D Domain
Code-to-Code Benchmarking of the PORFLOW and GoldSim Contaminant Transport Models using a Simple 1-D Domain - 11191 Robert A. Hiergesell and Glenn A. Taylor Savannah River National Laboratory SRNS Bldg.
More informationNumerical Solution of the Two-Dimensional Time-Dependent Transport Equation. Khaled Ismail Hamza 1 EXTENDED ABSTRACT
Second International Conference on Saltwater Intrusion and Coastal Aquifers Monitoring, Modeling, and Management. Mérida, México, March 3-April 2 Numerical Solution of the Two-Dimensional Time-Dependent
More informationBrian Berkowitz, 1 Andrea Cortis, 2 Marco Dentz, 3 and Harvey Scher 1
Click Here for Full Article MODELING NON-FICKIAN TRANSPORT IN GEOLOGICAL FORMATIONS AS A CONTINUOUS TIME RANDOM WALK Brian Berkowitz, 1 Andrea Cortis, 2 Marco Dentz, 3 and Harvey Scher 1 Received 18 July
More informationOn the advective-diffusive transport in porous media in the presence of time-dependent velocities
GEOPHYSICAL RESEARCH LETTERS, VOL. 3, L355, doi:.29/24gl9646, 24 On the advective-diffusive transport in porous media in the presence of time-dependent velocities A. P. S. Selvadurai Department of Civil
More informationContinuous Time Random Walks for Non-Local Radial Solute Transport
Continuous Time Random Walks for Non-Local Radial Solute Transport Marco Dentz Institute of Environmental Assessment and Water Research (IDÆA), Spanish National Research Council (CSIC), 08034 Barcelona,
More informationEffective pore-scale dispersion upscaling with a correlated continuous time random walk approach
WATER RESOURCES RESEARCH, VOL. 47, W12538, doi:10.1029/2011wr010457, 2011 Effective pore-scale dispersion upscaling with a correlated continuous time random walk approach T. Le Borgne, 1 D. Bolster, 2
More informationNon-Fickian dispersion in porous media: 2. Model validation from measurements at different scales
WATER RESOURCES RESEARCH, VOL. 44, W06427, doi:10.1029/2007wr006279, 2008 Non-Fickian dispersion in porous media: 2. Model validation from measurements at different scales T. Le Borgne 1,2 and P. Gouze
More informationSubordinated advection-dispersion equation for contaminant transport
WATER RESOURCES RESEARCH, VO. 37, NO. 6, PAGES 1543 1550, JUNE 2001 Subordinated advection-dispersion equation for contaminant transport Boris Baeumer, 1 David A. Benson, 2 Mark M. Meerschaert, 3 and Stephen
More informationReactive Transport in Porous Media
Reactive Transport in Porous Media Daniel M. Tartakovsky Department of Mechanical & Aerospace Engineering University of California, San Diego Alberto Guadagnini, Politecnico di Milano, Italy Peter C. Lichtner,
More informationAn alternative approach to simulate transport based on the Master equation
Tracers and Modelling in Hydrogeology (Proceedings of (lie TraM'2000 Conference held at Liège, Belgium, May 2000). IAHS Publ. no. 262, 2000. 121 An alternative approach to simulate transport based on the
More informationEffect of spatial concentration fluctuations on effective kinetics in diffusion-reaction systems
WATER RESOURCES RESEARCH, VOL.???, XXXX, DOI:10.1029/, Effect of spatial concentration fluctuations on effective kinetics in diffusion-reaction systems A. M. Tartakovsky Pacific Northwest National Laboratory,
More informationComparison of Heat and Mass Transport at the Micro-Scale
Comparison of Heat and Mass Transport at the Micro-Scale E. Holzbecher, S. Oehlmann Georg-August Univ. Göttingen *Goldschmidtstr. 3, 37077 Göttingen, GERMANY, eholzbe@gwdg.de Abstract: Phenomena of heat
More informationAn advective-dispersive stream tube approach for the transfer of conservative-tracer data to reactive transport
WATER RESOURCES RESEARCH, VOL. 36, NO. 5, PAGES 1209 1220, MAY 2000 An advective-dispersive stream tube approach for the transfer of conservative-tracer data to reactive transport Olaf A. Cirpka and Peter
More informationBME STUDIES OF STOCHASTIC DIFFERENTIAL EQUATIONS REPRESENTING PHYSICAL LAW
7 VIII. BME STUDIES OF STOCHASTIC DIFFERENTIAL EQUATIONS REPRESENTING PHYSICAL LAW A wide variety of natural processes are described using physical laws. A physical law may be expressed by means of an
More informationA MATLAB Implementation of the. Minimum Relative Entropy Method. for Linear Inverse Problems
A MATLAB Implementation of the Minimum Relative Entropy Method for Linear Inverse Problems Roseanna M. Neupauer 1 and Brian Borchers 2 1 Department of Earth and Environmental Science 2 Department of Mathematics
More informationMixing and reactions in porous media
Mixing and reactions in porous media Pietro De Anna To cite this version: Pietro De Anna. Mixing and reactions in porous media. Hydrology. Université Rennes 1, 2012. English. HAL Id: tel-00822932
More information1D Verification Examples
1 Introduction 1D Verification Examples Software verification involves comparing the numerical solution with an analytical solution. The objective of this example is to compare the results from CTRAN/W
More informationA numerical study of passive transport through fault zones
R Available online at www.sciencedirect.com Earth and Planetary Science Letters 214 (2003) 633^643 www.elsevier.com/locate/epsl A numerical study of passive transport through fault zones G.S. O Brien a;,
More informationWATER RESOURCES RESEARCH, VOL. 44, W05410, doi: /2008wr006858, 2008
Click Here for Full Article WATER RESOURCES RESEARCH, VOL. 44, W05410, doi:10.1029/2008wr006858, 2008 Transport of conservative solutes in simulated fracture networks: 2. Ensemble solute transport and
More informationThe distortion observed in the bottom channel of Figure 1 can be predicted from the full transport equation, C t + u C. y D C. z, (1) x D C.
1 8. Shear Dispersion. The transport models and concentration field solutions developed in previous sections assume that currents are spatially uniform, i.e. u f(,y,). However, spatial gradients of velocity,
More informationLagrangian characterization of contaminant transport through multidimensional heterogeneous media with limited heterogeneity information
Lagrangian characterization of contaminant transport through multidimensional heterogeneous media with limited heterogeneit information Zhang Yong 1, David A. Benson 1, Mark M. Meerschaert 2, Eric M. LaBolle
More informationLecture 25: Large Steps and Long Waiting Times
Lecture 25: Large Steps and Long Waiting Times Scribe: Geraint Jones (and Martin Z. Bazant) Department of Economics, MIT Proofreader: Sahand Jamal Rahi Department of Physics, MIT scribed: May 10, 2005,
More informationSimple closed form formulas for predicting groundwater flow model uncertainty in complex, heterogeneous trending media
WATER RESOURCES RESEARCH, VOL. 4,, doi:0.029/2005wr00443, 2005 Simple closed form formulas for predicting groundwater flow model uncertainty in complex, heterogeneous trending media Chuen-Fa Ni and Shu-Guang
More informationSTOCHASTIC CONTINUUM ANALYSIS OF GROUNDWATER FLOW PATHS FOR SAFETY ASSESSMENT OF A RADIOACTIVE WASTE DISPOSAL FACILITY
STOCHASTIC CONTINUUM ANALYSIS OF GROUNDWATER FLOW PATHS FOR SAFETY ASSESSMENT OF A RADIOACTIVE WASTE DISPOSAL FACILITY K. Chang*, C.L. Kim, E.Y. Lee, J.W.Park, H.Y.Park, C.G. Rhee, M.J. Song Nuclear Environment
More informationChemical Hydrogeology
Physical hydrogeology: study of movement and occurrence of groundwater Chemical hydrogeology: study of chemical constituents in groundwater Chemical Hydrogeology Relevant courses General geochemistry [Donahoe]
More information11280 Electrical Resistivity Tomography Time-lapse Monitoring of Three-dimensional Synthetic Tracer Test Experiments
11280 Electrical Resistivity Tomography Time-lapse Monitoring of Three-dimensional Synthetic Tracer Test Experiments M. Camporese (University of Padova), G. Cassiani* (University of Padova), R. Deiana
More informationBrownian motion and the Central Limit Theorem
Brownian motion and the Central Limit Theorem Amir Bar January 4, 3 Based on Shang-Keng Ma, Statistical Mechanics, sections.,.7 and the course s notes section 6. Introduction In this tutorial we shall
More informationTheory of anomalous chemical transport in random fracture networks
PHYSICAL REVIEW E VOLUME 57, NUMBER 5 MAY 1998 Theory of anomalous chemical transport in random fracture networks Brian Berkowitz and Harvey Scher Department of Environmental Sciences and Energy Research,
More informationDispersion, mixing and reaction in porous media Tanguy Le Borgne Geosciences Rennes, OSUR, Université de Rennes 1, France
Dispersion, mixing and reaction in porous media Tanguy Le Borgne Geosciences Rennes, OSUR, Université de Rennes 1, France Regis Turuban, Joaquin Jimenez, Yves Méheust, Olivier Bour, Jean-Raynald de Dreuzy,
More informationCHARACTERIZATION OF HETEROGENEITIES AT THE CORE-SCALE USING THE EQUIVALENT STRATIFIED POROUS MEDIUM APPROACH
SCA006-49 /6 CHARACTERIZATION OF HETEROGENEITIES AT THE CORE-SCALE USING THE EQUIVALENT STRATIFIED POROUS MEDIUM APPROACH Mostafa FOURAR LEMTA Ecole des Mines de Nancy, Parc de Saurupt, 54 04 Nancy, France
More informationOn using the equivalent conductivity to characterize solute spreading in environments with low-permeability lenses
WATER RESOURCES RESEARCH, VOL. 38, NO. 8, 1132, 10.1029/2001WR000528, 2002 On using the equivalent conductivity to characterize solute spreading in environments with low-permeability lenses Andrew J. Guswa
More informationNumerical simulation of non-fickian transport in geological formations with multiple-scale heterogeneities
WATER RESOURCES RESEARCH, VOL. 40, W04209, doi:10.1029/2003wr002750, 2004 Numerical simulation of non-fickian transport in geological formations with multiple-scale heterogeneities Andrea Cortis Department
More informationCHAPTER 11: CHROMATOGRAPHY FROM A MOLECULAR VIEWPOINT
CHAPTER 11: CHROMATOGRAPHY FROM A MOLECULAR VIEWPOINT Contrasting approaches 1. bulk transport (e.g., c c = W ) t x + D c x goal: track concentration changes advantage: mathematical rigor (for simple models).
More informationReaction time distributions in chemical kinetics: Oscillations and other weird behaviors
Introduction The algorithm Results Summary Reaction time distributions in chemical kinetics: Oscillations and other weird behaviors Ramon Xulvi-Brunet Escuela Politécnica Nacional Outline Introduction
More informationImproved inverse modeling for flow and transport in subsurface media: Combined parameter and state estimation
GEOPHYSICAL RESEARCH LETTERS, VOL. 32, L18408, doi:10.1029/2005gl023940, 2005 Improved inverse modeling for flow and transport in subsurface media: Combined parameter and state estimation Jasper A. Vrugt,
More informationWind and turbulence experience strong gradients in vegetation. How do we deal with this? We have to predict wind and turbulence profiles through the
1 2 Wind and turbulence experience strong gradients in vegetation. How do we deal with this? We have to predict wind and turbulence profiles through the canopy. 3 Next we discuss turbulence in the canopy.
More informationDiffusion. CS/CME/BioE/Biophys/BMI 279 Nov. 15 and 20, 2016 Ron Dror
Diffusion CS/CME/BioE/Biophys/BMI 279 Nov. 15 and 20, 2016 Ron Dror 1 Outline How do molecules move around in a cell? Diffusion as a random walk (particle-based perspective) Continuum view of diffusion
More informationHandbook of Stochastic Methods
C. W. Gardiner Handbook of Stochastic Methods for Physics, Chemistry and the Natural Sciences Third Edition With 30 Figures Springer Contents 1. A Historical Introduction 1 1.1 Motivation I 1.2 Some Historical
More informationIntroduction: the plague. The Scaling Laws Of Human Travel. Introduction: SARS (severe acute respiratory syndrom)
The Scaling Laws Of Human Travel Verena Zuber Department of Statistics University of Munich 7. July 2006 Introduction: the plague highly infectious disease several pandemics have influenced the human history
More informationPseudo-seismic wavelet transformation of transient electromagnetic response in engineering geology exploration
GEOPHYSICAL RESEARCH LETTERS, VOL. 34, L645, doi:.29/27gl36, 27 Pseudo-seismic wavelet transformation of transient electromagnetic response in engineering geology exploration G. Q. Xue, Y. J. Yan, 2 and
More informationDiffusion and cellular-level simulation. CS/CME/BioE/Biophys/BMI 279 Nov. 7 and 9, 2017 Ron Dror
Diffusion and cellular-level simulation CS/CME/BioE/Biophys/BMI 279 Nov. 7 and 9, 2017 Ron Dror 1 Outline How do molecules move around in a cell? Diffusion as a random walk (particle-based perspective)
More informationLagrangian simulation of mixing and reactions in complex geochemical systems
Lagrangian simulation of mixing and reactions in complex geochemical systems Nicholas B. Engdahl 1, David A. Benson 2, and Diogo Bolster 3 1 Civil and Environmental Engineering, Washington State University
More informationLecture 13. Drunk Man Walks
Lecture 13 Drunk Man Walks H. Risken, The Fokker-Planck Equation (Springer-Verlag Berlin 1989) C. W. Gardiner, Handbook of Stochastic Methods (Springer Berlin 2004) http://topp.org/ http://topp.org/species/mako_shark
More informationLarge-Time Spatial Covariance of Concentration of Conservative Solute and Application to the Cape Cod Tracer Test
Transport in Porous Media 42: 109 132, 2001. 2001 Kluwer Academic Publishers. Printed in the Netherlands. 109 Large-Time Spatial Covariance of Concentration of Conservative Solute and Application to the
More informationApplication of the theory of micromixing to groundwater reactive transport models
WATER RESOURCES RESEARCH, VOL. 39, NO. 11, 1313, doi:10.1029/2003wr002368, 2003 Application of the theory of micromixing to groundwater reactive transport models Bruce A. Robinson and Hari S. Viswanathan
More informationWeak Ergodicity Breaking. Manchester 2016
Weak Ergodicity Breaking Eli Barkai Bar-Ilan University Burov, Froemberg, Garini, Metzler PCCP 16 (44), 24128 (2014) Akimoto, Saito Manchester 2016 Outline Experiments: anomalous diffusion of single molecules
More informationVolatility and Returns in Korean Futures Exchange Markets
Volatility and Returns in Korean Futures Exchange Markets Kyungsik Kim*,, Seong-Min Yoon and Jum Soo Choi Department of Physics, Pukyong National University, Pusan 608-737, Korea Division of Economics,
More informationChapter 2 Theory. 2.1 Continuum Mechanics of Porous Media Porous Medium Model
Chapter 2 Theory In this chapter we briefly glance at basic concepts of porous medium theory (Sect. 2.1.1) and thermal processes of multiphase media (Sect. 2.1.2). We will study the mathematical description
More informationAnomalous diffusion of volcanic earthquakes
Anomalous diffusion of volcanic earthquakes SUMIYOSHI ABE 1,2 and NORIKAZU SUZUKI 3 1 Department of Physical Engineering, Mie University, Mie 514-8507, Japan 2 Institute of Physics, Kazan Federal University,
More informationPhysical pictures of transport in heterogeneous media: Advection-dispersion, random-walk, and fractional derivative formulations
WATER RESOURCES RESEARCH, VOL. 38, NO. 1, 1191, doi:1.129/21wr13, 22 Physical pictures of transport in heterogeneous media: Advection-dispersion, random-walk, and fractional derivative formulations Brian
More informationModeling Gas Flooding in the Presence of Mixing Dr. Russell T. Johns Energy and Minerals Engineering The Pennsylvania State University
Vienna 2011 Modeling Gas Flooding in the Presence of Mixing Dr. Russell T. Johns Energy and Minerals Engineering The Pennsylvania State University Why is Mixing Important? A) Mixing of gas and oil alters
More informationContinuous and Discrete random process
Continuous and Discrete random and Discrete stochastic es. Continuous stochastic taking values in R. Many real data falls into the continuous category: Meteorological data, molecular motion, traffic data...
More informationAlgebraic Turbulence-Chemistry Interaction Model
50th AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition 09-12 January 2012, Nashville, Tennessee AIAA 2012-0183 50th AIAA Aerospace Sciences Meeting, January 9-12,
More informationLecture 3: From Random Walks to Continuum Diffusion
Lecture 3: From Random Walks to Continuum Diffusion Martin Z. Bazant Department of Mathematics, MIT February 3, 6 Overview In the previous lecture (by Prof. Yip), we discussed how individual particles
More informationA colloid-facilitated transport model with variable colloid transport properties
GEOPHYSICAL RESEARCH LETTERS, VOL. 34, L09401, doi:10.1029/2007gl029625, 2007 A colloid-facilitated transport model with variable colloid transport properties Bruce A. Robinson, 1 Andrew V. Wolfsberg,
More informationA SPATIAL STRUCTURAL DERIVATIVE MODEL FOR ULTRASLOW DIFFUSION
THERMAL SCIENCE: Year 7, Vol., Suppl., pp. S-S7 S A SPATIAL STRUCTURAL DERIVATIVE MODEL FOR ULTRASLOW DIFFUSION by Wei XU a, Wen CHEN a*, Ying-Jie LIANG a*, and Jose WEBERSZPIL b a State Key Laboratory
More informationThe Truth about diffusion (in liquids)
The Truth about diffusion (in liquids) Aleksandar Donev Courant Institute, New York University & Eric Vanden-Eijnden, Courant In honor of Berni Julian Alder LLNL, August 20th 2015 A. Donev (CIMS) Diffusion
More informationTurbulence is a ubiquitous phenomenon in environmental fluid mechanics that dramatically affects flow structure and mixing.
Turbulence is a ubiquitous phenomenon in environmental fluid mechanics that dramatically affects flow structure and mixing. Thus, it is very important to form both a conceptual understanding and a quantitative
More informationNonlinear effects of salt concentrations on evaporation from porous media
GEOPHYSICAL RESEARCH LETTERS, VOL. 39,, doi:10.1029/2011gl050763, 2012 Nonlinear effects of salt concentrations on evaporation from porous media Mansoureh Norouzi Rad 1 and Nima Shokri 1 Received 23 December
More informationTransport of Chemicals, Kinetics, and Equilibrium
Transport of Chemicals, Kinetics, and Equilibrium Module 1: Basic Concepts, Lecture Chemical Fate and Transport in the Environment, nd edition. H.F. Hemond and E.J. Fechner-Levy. Academic Press. London.
More informationPUBLICATIONS. Water Resources Research. Nickel migration and retention dynamics in natural soil columns RESEARCH ARTICLE 10.
PUBLICATIONS Water Resources Research RESEARCH ARTICLE Key Points: Retention by sorption and precipitation is needed to quantify dynamic processes Preferential pathways alter retention and overall nickel
More informationPoint Propagation Model of Pollution Sources of Heavy Metals in the Surface Soil
Int. J. Electrochem. Sci., 8 (013) 5153-5160 International Journal of ELECTROCHEMICAL SCIENCE www.electrochemsci.org Point Propagation Model of Pollution Sources of Heavy Metals in the Surface Soil Hai-ming
More informationHandbook of Stochastic Methods
Springer Series in Synergetics 13 Handbook of Stochastic Methods for Physics, Chemistry and the Natural Sciences von Crispin W Gardiner Neuausgabe Handbook of Stochastic Methods Gardiner schnell und portofrei
More informationSTOCHASTIC CHEMICAL KINETICS
STOCHASTIC CHEICAL KINETICS Dan Gillespie GillespieDT@mailaps.org Current Support: Caltech (NIGS) Caltech (NIH) University of California at Santa Barbara (NIH) Past Support: Caltech (DARPA/AFOSR, Beckman/BNC))
More informationTRACER TESTS FOR VARIOUS CARBONATE CORES USING X-RAY CT
SCA005-56 / TRACER TESTS FOR VARIOUS CARBONATE CORES USING X-RAY CT M. Fourar,, G. Konan, C. Fichen, E. Rosenberg, P. Egermann, and R. Lenormand Ecole des Mines de Nancy Laboratoire d Energétique et de
More informationEVALUATION OF CRITICAL FRACTURE SKIN POROSITY FOR CONTAMINANT MIGRATION IN FRACTURED FORMATIONS
ISSN (Online) : 2319-8753 ISSN (Print) : 2347-6710 International Journal of Innovative Research in Science, Engineering and Technology An ISO 3297: 2007 Certified Organization, Volume 2, Special Issue
More informationModeling of 1D Anomalous Diffusion In Fractured Nanoporous Media
LowPerm2015 Colorado School of Mines Low Permeability Media and Nanoporous Materials from Characterisation to Modelling: Can We Do It Better? IFPEN / Rueil-Malmaison - 9-11 June 2015 CSM Modeling of 1D
More informationComparison of Heat and Mass Transport at the Micro-Scale
Comparison of Heat and Mass Transport at the Micro-Scale Ekkehard Holzbecher, Sandra Oehlmann October 10 th, 2012 Excerpt from the Proceedings of the 2012 COMSOL Conference in Milan Heat & Mass Transfer
More informationTemporal evolution of pore geometry, fluid flow, and solute transport resulting from colloid deposition
WATER RESOURCES RESEARCH, VOL. 45,, doi:10.1029/2008wr007252, 2009 Temporal evolution of pore geometry, fluid flow, and solute transport resulting from colloid deposition Cheng Chen, 1,2 Boris L. T. Lau,
More information3-DPTV Experiments of Anomalous, Steady Transport of a Conservative Tracer in Homogeneous and Heterogeneous Porous Media
3-DPTV Experiments of Anomalous, Steady Transport of a Conservative Tracer in Homogeneous and Heterogeneous Porous Media by A. Cenedese, J.H. Cushman & and M. Moroni University of Rome La Sapienza, DITS,
More informationA breakthrough curve analysis of unstable density-driven flow and transport in homogeneous porous media
WATER RESOURCES RESEARCH, VOL. 40,, doi:10.1029/2003wr002668, 2004 A breakthrough curve analysis of unstable density-driven flow and transport in homogeneous porous media M. Wood, C. T. Simmons, and J.
More informationCO 2 storage capacity and injectivity analysis through the integrated reservoir modelling
CO 2 storage capacity and injectivity analysis through the integrated reservoir modelling Dr. Liuqi Wang Geoscience Australia CO 2 Geological Storage and Technology Training School of CAGS Beijing, P.
More informationWeak Ergodicity Breaking WCHAOS 2011
Weak Ergodicity Breaking Eli Barkai Bar-Ilan University Bel, Burov, Korabel, Margolin, Rebenshtok WCHAOS 211 Outline Single molecule experiments exhibit weak ergodicity breaking. Blinking quantum dots,
More informationSuperstatistics and temperature fluctuations. F. Sattin 1. Padova, Italy
Superstatistics and temperature fluctuations F Sattin 1 Padova, Italy Abstract Superstatistics [C Beck and EGD Cohen, Physica A 322, 267 (2003)] is a formalism aimed at describing statistical properties
More informationDisordered Structures. Part 2
Disordered Structures Part 2 Composites and mixtures Consider inhomogeneities on length scales > 10-20 Å Phase separation two (or multi-) phase mixtures Mixtures of particles of different kinds - solids,
More informationPollution. Elixir Pollution 97 (2016)
42253 Available online at www.elixirpublishers.com (Elixir International Journal) Pollution Elixir Pollution 97 (2016) 42253-42257 Analytical Solution of Temporally Dispersion of Solute through Semi- Infinite
More informationRandom walks. Marc R. Roussel Department of Chemistry and Biochemistry University of Lethbridge. March 18, 2009
Random walks Marc R. Roussel Department of Chemistry and Biochemistry University of Lethbridge March 18, 009 1 Why study random walks? Random walks have a huge number of applications in statistical mechanics.
More informationCover times of random searches
Cover times of random searches by Chupeau et al. NATURE PHYSICS www.nature.com/naturephysics 1 I. DEFINITIONS AND DERIVATION OF THE COVER TIME DISTRIBUTION We consider a random walker moving on a network
More informationTheory of fractional Lévy diffusion of cold atoms in optical lattices
Theory of fractional Lévy diffusion of cold atoms in optical lattices, Erez Aghion, David Kessler Bar-Ilan Univ. PRL, 108 230602 (2012) PRX, 4 011022 (2014) Fractional Calculus, Leibniz (1695) L Hospital:
More informationA Decomposition Method for Solving Coupled Multi-Species Reactive Transport Problems
Transport in Porous Media 37: 327 346, 1999. c 1999 Kluwer Academic Publishers. Printed in the Netherlands. 327 A Decomposition Method for Solving Coupled Multi-Species Reactive Transport Problems YUNWEI
More informationDETERMINATION OF MODEL VALID PREDICTION PERIOD USING THE BACKWARD FOKKER-PLANCK EQUATION
.4 DETERMINATION OF MODEL VALID PREDICTION PERIOD USING THE BACKWARD FOKKER-PLANCK EQUATION Peter C. Chu, Leonid M. Ivanov, and C.W. Fan Department of Oceanography Naval Postgraduate School Monterey, California.
More informationApplication of high-order spectral method for the time fractional mobile/immobile equation
Computational Methods for Differential Equations http://cmde.tabrizu.ac.ir Vol. 4, No. 4, 2016, pp. 309-322 Application of high-order spectral method for the time fractional mobile/immobile equation Hossein
More informationThe Role of Dispersion in Radionuclide Transport - Data and Modeling Requirements
S-N/99205-003 Shaw/13052-200 The Role of Dispersion in Radionuclide Transport - Data and Modeling Requirements Revision No.: 1 February 2004 Prepared for U.S. Department of Energy under Contract No. DE-AC52-03NA99205
More informationModeling of dispersed phase by Lagrangian approach in Fluent
Lappeenranta University of Technology From the SelectedWorks of Kari Myöhänen 2008 Modeling of dispersed phase by Lagrangian approach in Fluent Kari Myöhänen Available at: https://works.bepress.com/kari_myohanen/5/
More informationWhen do diffusion-limited trajectories become memoryless?
When do diffusion-limited trajectories become memoryless? Maciej Dobrzyński CWI (Center for Mathematics and Computer Science) Kruislaan 413, 1098 SJ Amsterdam, The Netherlands Abstract Stochastic description
More informationComputational modelling of reactive transport in hydrogeological systems
Water Resources Management III 239 Computational modelling of reactive transport in hydrogeological systems N. J. Kiani, M. K. Patel & C.-H. Lai School of Computing and Mathematical Sciences, University
More informationNew Fast Kalman filter method
New Fast Kalman filter method Hojat Ghorbanidehno, Hee Sun Lee 1. Introduction Data assimilation methods combine dynamical models of a system with typically noisy observations to obtain estimates of the
More informationTurbulent Flows. g u
.4 g u.3.2.1 t. 6 4 2 2 4 6 Figure 12.1: Effect of diffusion on PDF shape: solution to Eq. (12.29) for Dt =,.2,.2, 1. The dashed line is the Gaussian with the same mean () and variance (3) as the PDF at
More informationAnomalous diffusion in biology: fractional Brownian motion, Lévy flights
Anomalous diffusion in biology: fractional Brownian motion, Lévy flights Jan Korbel Faculty of Nuclear Sciences and Physical Engineering, CTU, Prague Minisymposium on fundamental aspects behind cell systems
More information