= _(2,r)af OG(x, a) 0p(a, y da (2) Aj(r) = (2*r)a (Oh(x)y,(y)) A (r) = -(2,r) a Ox Oxj G(x, a)p(a, y) da
|
|
- Annabel Elliott
- 5 years ago
- Views:
Transcription
1 WATER RESOURCES RESEARCH, VOL. 35, NO. 7, PAGES , JULY 999 A general method for obtaining analytical expressions for the first-order velocity covariance in heterogeneous porous media Kuo-Chin Hsu Daniel B. Stephens and Associates, Inc., Albuquerque, New Mexico Abstract. Velocity covariance of seepage flow is the essential quantity in understanding the uncertainty of groundwater flow and contaminant transport in heterogeneous porous media. Nevertheless, a closed-form analytical expression is available only for first-order analysis. These available expressions are limited to some specific correlation functions of natural log hydraulic conductivity. The present paper develops a general method for obtaining the first-order velocity covariance for any statistically isotropic log-hydraulic conductivity correlation function. The general analytical expressions are derived for both two and three dimensions. The results are demonstrated for several frequently used correlation functions.. Introduction Velocity covariance of groundwater is important in understanding the uncertainty of flow and transport in heterogeneous porous media. For example, calculating the dispersion tensor requires knowing the velocity covariance [Dagan, 989]. The first-order (in the variance 0 '2 of the naturalog hydraulic conductivity Y) velocity covariance for the exponential correlation Y has been derived in two dimensions by Rubin [990] and Hsu et al. [996] with seemingly different but actually equivalent expressions. For the same correlation function, Rubin and Dagan [992] and Zhang and Neuman [992] have derived the corresponding three-dimensional expressions. As an example of a different correlation function for Y, Hsu and Neuman [997] used a Gaussian model in an attempt to derive the first-order velocity covariance. Hsu et al. [996] and Hsu and Neuman [997] investigated higher-order effects on velocity covariance by using a semianalytical method, while Deng and Cushman [995, 998] used a numerical method to obtain the velocity covariance in the Fourier domain. The closed-form expression of velocity covariance is only available for firstorder analysis. Commonly, in previous research, the analytical solutions for velocity covariance require a specific correlation function for Y. The present study utilizes a Green's function approach to derive the general forms of the first-order velocity covariance. Expressions for both two and three dimensions are derived and are applicable to any statistically isotropic correlation function for Y. 2. First-Order Velocity Covariance Theoretical derivations of velocity covariance usually assume () infinite domain, (2) uniform mean flow, (3) steady state flow, and (4) natural log hydraulic conductivity Y as a stationary homogeneous, multivariate Gaussian random field with constant ensemble mean (Y) and variance 0 '2. With these assumptions, together with the assumption that mean flow takes place parallel to the x direction, we begin with the general expression for the first-order velocity covariance given in (A5) of Hsu et al. [996], i.e., Copyright 999 by the American Geophysical Union. Paper number 999WR /99/999 WR9007 $ /3/v4 (r) + B; (r)] () where ui, is the first-order groundwater velocity covariance, Kg is the geometric mean of hydrauliconductivity K, J is the uniform mean driving force (negative hydraulic gradient), qb is the effective porosity (treated as a constant), /3 is the Kronecker delta, d is the number of space dimensions, and r = x - y is the separation vector of points x and y. Functions.4 and B may be related to the normalized log transmissivity-head covariance and head variogram by comparing () with equation (6) of Rubin [990]. Functions.4 and B can be expressed in terms of Green's function G, via equations (All) and (A2) of Hsu et al. [996] as Aj(r) = (2*r)a (Oh(x)y,(y)) Jo '2 Oxj = _(2,r)af OG(x, a) 0p(a, y da (2) Oxj Oa, Bj (r) = (2 r)d J2o'2 (Oh(x)Oh(y)) OXj Oyk =(2)dff OG(x, Oxj a)og(y, Oyk a')02p(a, Oa Oa dada' (3) Here h is the first-order (in o-) hydraulic head fluctuation, Y' is the random fluctuation of log hydraulic conductivity, p is the log conductivity correlation function, and angle brackets represent the ensemble mean (expected value). Applying integration by parts, the integrated terms vanish as a, a' - +_% and (2) and (3) can thus be written as A (r) = -(2,r) a Ox Oxj G(x, a)p(a, y) da O2(I) (r) = -(2,r) a -- (4) OX 0Xj
2 2274 HSU: TECHNICAL NOTE Table. Expressions for Function Wi Function W Wlo W Expression p(r) - 2T(r)/r p(r)/8 + T(r)/2r- 3U(r)/2 p(r) + T(r)/r p(r) - 3T(r)/r 3p(r) + 3T(r)/2r- 9U(r)/2 p(r) - 9T(r)/2r + 5U(r)/2 p(r) + 5T(r)/2r- 35U(r)/2 p(r) - 3T(r)/2r- 5U(r)/2 p(r) - 5T(r)/2r + 9U(r)/2-3T(r)/2r + 5U(r)/2 p(r) + 3T(r)/2r- 5U(r)/2 For W and W2, T(r) - (l/r) f[ r'p(r') dr' (from (0)), and U(r) = (/r 4) fro r'3p(r') dr'. For W3-W o, T(r) = (/r 2) fro r'2p(r') dr' (from(22)), and e(r) - (/r s) f[ r'4p(r ') dr'. Bjk(r) = (2,r) a 0 4 OxiOylOxjOyk ß ff G(x, a)g(a', y)p(a, a') da da' Expanding the right-hand side of (2) and collecting terms with the same angular dependence, we find that (2) becomes (2,r) 2 A,(r) = 2 [o(r) + W,(r) cos 20] (3) where function W is shown in Table. Similarly, A 2 can be obtained in the form A2(r) = -(2 r) 2 cos ---sin 0 0 0) 0 r r - ß sin0 - + cos0 0 0) (r) (2,tr) 2 2 Wa(r) sin 20 (4) To determine the function B, we begin by writing (9) in polar coordinates and differentiating to obtain T(r) =- d = dr 7 - r - - d d2xi t d3xi t - r d dr 3 (5) where - (2*r)a Oxi Oyl 9xj Oy - G(a', y) (x, a') da' 04xit (r) = (2*r)a 02X0Xj OX k (5) ß (r) = I G(x, a)p(a, y) da (6) ß (r) = f G(a', y)ci>(x, a') da' The functions ß and ß satisfy V2 (r) = --p(r) V2 (r) = - (r) Since p depends merely on distance r = Ir], only the radial component of the Laplacian is required. For two dimensions a first integration of (8) then yields d (r) df F 0 t' (7) (8) (9) r'p(r') dr' = -T(r) (0) and the function T can be related to p through dt T d-7+¾ = () Introducing polar coordinates r - r (cos O, sin O) r and specializing to the case j =, we have 02 (r) Al(r) = --(2'n') 2 -- OXl sin 0 ci) 0 0) 2 = --(2'n') 2 COS 0 r (r) (2) Specializing to the case j = k =, after manipulation by using () and (5), we find that (5) becomes 04xXt (r) B (r) = (2,r) 2-- Oxl ( 0 r 0 r 0) 4 = (2 r) 2 COS sin 0 (r) (2 r)2i 3 p(r) + W (r) 2 cos 20 + W2(r ) cos ] 4 0 (6) where functions W and W2 are defined in Table. When this procedure is applied successively to the cases (j, k) = (, 2), (2, ), and (2, 2), we find B 2(r) = B2(r) = (2 r)2[ W (r) 4 sin20 + W2(r) sin40] (7) B22(r) = (2'rr)28-p(r)- W2(r) cos 40] (8) Substituting (3), (4), (6), (7), and (8) into (), we obtain the general two-dimensional expressions of first-order velocity covariance for any stationary isotropic function p in the form Ull(r ) = Ky2o-2 2 [3. p(r) - Wi(r) cos 20 + W2(r) cos 40] (9) KY 20'2 W,(r) sin 20 + W2(r) sin 40] (20) t,t 2(r) -- q[2 [-- K2-2 2 gd O' u22(r ) = qb 2 [ -p(r) - W2 cos 40] (2) These expressions are valid for any stationary correlation function of Y for which the resulting integrals converge. For three dimensions, applying spherical coordinates to (8), we have
3 HSU: TECHNICAL NOTE 2275 d (r) df r2 r' 2p(r' ) dr' = - T(r) (22) B3(r) = 32 {sin 2½[6W6(r) cos 0 + 2W7(r) cos 30] dt dr T + 2- = p (23) Introducing the spherical coordinates r = r (sin ½ cos 0, sin ½ sin 0, cos ½)r, so that 0 0 sin O x --= sin ½ cos 0 (24) r Or rsin½00 +/cøs½cøso cos =sin½sino + - cos½sino (25) O y r sin ½ O0 7 - o o o... sin ½ (26) oz-cøs½ r r we find that the functions A become (2 ) 3 A (r) = 4 [W3(r) + W4(r) 0½ ß (-cos 2½ + cos 20 - cos 2½ cos 20)] (27) A2(r) = (2 ) 3 W4(r)(sin 20 - cos 2½ sin 20) (28) where functions W3 and W4 can be found in Table. To derive the function B, we write (5) in spherical coordinates as T(r) = - dti) d = dr d? - d ' (30) u2(r ) = qb 2 sin W 6 ( d ) 2 d2xi 2 r + d3xi dr 3 r r d- r By manipulations similar to those used in two dimensions, the functions B can be obtained in the form Bll(F) = {3Ws(r)- 2W6(r) cos 2½ + 3W7(r) cos 4½ + cos 202W6(r) - 6Wa(r)cos 2½ + 4W7(r)cos 4½] + W7(r) cos 40(3-4 cos 2½ + cos 4½)} (3) B22(r) = (2 ) 3 [Ws(r) - 4W6(r) cos 2½ + W7(r) cos 4½ + W7(r) cos 40( cos 2½- cos 4½)] (32) B33(r ) = {W9(r) + 4W0(r) cos 2½- W?(r) cos 4½ 6 + COS 20[Wll(r) - 4W0(r) cos 2½ - W7(r) cos 4½]} B2(r) = -- {6 W6(r) sin W7(r) sin 40 - cos 2½[8Wa(r) sin W7(r) sin 40] (33) + W7(r) cos 4½(2 sin 20 + sin 40)} (34) B23(r) = - W7(r) sin 4½(3 cos 0 + cos 30)} (35) 32 {2 sin 2½[W6(r) sin 0 + Wv(r) sin 30] - Wv(r) sin 4½(sin 0 + sin 30)} (36) where the functions W (i - 5, 6,.-' ) are defined in Table. Equations (),(27)-(29), and (3)-(36) constitute the threedimensional first-order velocity covariance function with the following expressions, #d O' W 3 3 W 4 3 u(r) = qb2 p- ' -+ - Ws W6 cos 2½ + Wvcos4½+cos20 I( W6 ) + 2 cos 2½ + + cos 40(3 g3j2o '2 cos 4½ ] -4co8 2½ + cos 4½)} (37) u22(r) = qb2 [W5-4W6 cos 2½ + W7 cos 4½ A3(r) = 2 W4(r) sin 2½ cos 0 (29) + W7 cos 40( cos 2½- cos 4½)] (38) g3j2o '2 tt33(r ) = 6qb 2 [W 9 q- 4W0 cos 2½ - W7 cos 4½ + cos 20(Wll - 4W0 cos 2½ - W7 cos 4½)] (39) K2-22{ [ #d O' W cos 2½ + -cos 4½ + - sin 40(3-4 cos 2½ + cos 4½) (40) [( ) #d O' W 4 3 u 3(r) - tk2 cos 0 - -q- - W6 sin 2½ W7 sin 4½ q- cos 30(sin 2½ - sin 4½) (4) g3j2 rr 2 u23(r) = 324)2 [sin 0(2W6 sin 2½- W, sin 4½) q- W7 sin 30(2 sin 2½- sin 4½)] (42) These expressions are valid for any isotropicorrelation function p for which the resulting integrals converge. 3. Applications We apply the results derived above to commonly used isotropic correlation functions, i.e., exponential, Gaussian, spherical, and linear correlation functions IDagan, 989], in both two and three dimensions. These results are readily utilized for the nested correlation function as well. The expressions for isotro-
4 2276 HSU: TECHNICAL NOTE pic exponential, Gaussian, spherical, and linear correlation functions are, respectively, For the two-dimensional isotropic Gaussian correlation p(r) = e -an p(r) = e -br2 p(r) =- 53r +5 ( ) 2 p(r) = 0 >.. p(r) = I <. p(r) = 0 >. < (43) (44) (45) (46) 3.. Two-Dimensional Isotropic Correlation Functions Substituting (56) and (57) into (9)-(2), we recover the expressions in (70), (7), and (72) of Hsu and Neuman [997] In two dimensions, where r 2 = x 2 + x22, the isotropic derived by using Fourier transforms. exponential correlation function (43) yields, from (0), For the two-dimensional isotropic spherical correlation function, we obtain T(r) = a -- 2' [ - ( + ar)e-ar] (47) and, from Table, S ( r ) - - a - r + a- r 63 + a- - r ) e-ar + a r 6 4 (48) Therefore functions W and W2 become, on using the relation a = /X in which X is the integral scale, W = e -r/x (49) r X 2 9A 4 ( X 4A 29A 39 4 W2 = 2r 2 r e -r/x (50) Substituting (49) and (50) into (9)-(2), we obtain the firstorder veloci covarianc expressions. U (r) 3 X r = 2 e-r/x e r/x -2 ß cos20+ ( +-+4 X r ) ß e -r/x cos 40 (5) g 2r2 a 2 { -r/x --[( + -- X ) ß e -r/ cos40 } (52) K2-2 gj 2 { ( X r 2) --2 2] sin20 + [( +-+ r ) e -r/x ] sin40 } (53) u22(r) = 2 e U 2(r) = r -r/x T(r) = 2-- r (- e -br2) (54) r ) S(.) = - 2b2r + 2 e -br2 + 2b2r 4 (55) Substituting (44),(54), and (55) into functions W and W2, and relating b to the integral scale X by b -,r/4x 2, we obtain W -- I q- - F2) 4A2 e -*rr2/4a2,7.f 2 2A 7-/- 2 2A4' 7-2/- 4 ) - 7-? X 2-2 2A 7.2/, (56) W2 = e- 2/4 2 (57) T(r) = r (rlr3) - j + - g r < 2 r T(r) = l Or > 3r r 3 r U(r) = 4 0 l 4 l 3 < 3 4 r U(r) = 40 '4 [ > (s8) (59) Substituting (45), (58), and (59) into function W, W2, and (9)-(2), we then obtain the velocity covariance function. For the two-dimensional isotropic linear correlation function (46), we have T(r) = r 5-5 < 2 r T(r) = 6 r > lr r U(r) = 4 5 l [ < 4 r U(p) = 20 p4 > (60) (6) The velocity covariance function is calculated by using (46) and the above two equations Three-Dimensional Isotropic Correlation Functions For three dimensions, where r 2 = x 2 + x22 + x 23, the isotropic exponential correlation function (43), yields T, from (22), in the form T(r) = a - + a -- 2' + a' 2 e-ar (62) and, from Table, Equations (5)-(53) are given in a different but equivalent 24 form by Rubin [990, equations (9), (0), (2)] and Hsu et al. [996, equations (), (2), (3)]. U(F) a ar a r 2 a- r 3 a a ) = q- q- q- q- e -an (63)
5 HSU: TECHNICAL NOTE 2277 Substituting (43), (62), and (63) into (27)-(29) and (3)-(36), we obtain the first-order velocity covarianc expressions, which are given in a different but equivalent form by Rubin and Dagan [992, equations (A2)-(A8)] and Zhang and Neuman [992, equations ()-(4)]. For the three dimensional isotropic Gaussian correlation x/- r] ) r(r) = - 2b e-br2 q- 4b3/2r 2 erf [ ( ( 2 + 4b2r 34 ) e -br2 + 8bS/2r 3x/ 5 err [ x/- r] ( 65) U(r) = r where erf is the error function. Substituting (43) and ()-(65) into (27)-(29) and (3)-(36), we obtain the expressions for the isotropic Gaussian model. For 0-0 ø, ½ = 90 ø, and relating b to the integral scale X by b =,r/4x 2, we recover the results of (6)-(63) of Hsu and Neuman [997] derived by using the Fourier transform. For the three-dimensional isotropic spherical correlation r(r)=r (3rlr3) 5- [+ - x r [< l 3 r T(r) = 24 r 2 > lr r 3 r U(r) = - [ +63 < l 5 r U(r) = 80 r s [ > (66) (67) Equations (45), (66), and (67) are now used for the calculation of functions Wi, and the velocity covariance function (37)- (42). For the three-dimensional isotropic linear correlation function, we have r(r) = r ( 5- ) r [ < 3 r T(r) = 2 r 2 > (68) The corresponding velocity covariance function is readily calculated by using (46), (68), and (69). 4. Summary We have utilized the general expressions of the first-order velocity covariance given by Hsu et al. [996] to derive the ready-to-use expressions in both two and three dimensions. These expressions can be applied to any stationary isotropic correlation function of natural log-hydraulic conductivity for which the resulting integrals converge. We have applied the results to exponential, Gaussian, spherical, and linear functions in two and three dimensions. Some of these expressions are presented for the first time; others appear in the same or different forms given in previous research. The first-order velocity covariance expressions derived here offer a general method for obtaining analytical expressions for first-order seepage velocity covariance in heterogeneous porous media. Acknowledgments. The author thanks George L. Lamb Jr. for helpful discussion in developing the methodology. Comments from Dongxiao Zhang and an anonymous reviewer are appreciated. References Dagan, G., Flow and Transport in Porous Formations, 465 pp., Springer-Verlag, New York, 989. Deng, F.-W., and J. H. Cushman, On higher-order corrections to the flow velocity covariance tensor, Water Resour. Res., 3(7), , 995. Deng, F.-W., and J. H. Cushman, Higher-order correction to the flow velocity covariance tensor, revisited, Water Resour. Res., 34(), 03-06, 998. Hsu, K.-C., and S. P. Neuman, Second-order expressions for velocity moments in two- and three-dimensional statistically anisotropic media, Water Resour. Res., 33(4), , 997. Hsu, K.-C., D. Zhang, and S. P. Neuman, Higher-order effects on flow and transport in randomly heterogeneous media, Water Resour. Res., 32(3), , 996. Rubin, Y., Stochastic modeling of macrodispersion in heterogeneous porous media, Water Resour. Res., 26(), 33-4, 990. Rubin, Y., and G. Dagan, A note on head and velocity covariances in three-dimensional flow through heterogeneous anisotropic porous media, Water Resour. Res., 28(5), , 992. Zhang, D., and S. P. Neuman, Comment on "A note on head and velocity covariances in three-dimensional flow through heterogeneous anisotropic porous media" by Y. Rubin and G. Dagan, Water Resour. Res., 28(2), , 992. lr r U(r) =5 6 l [ < i l 5 r U(r) = 30 r s J > (69) K.-C. Hsu, Daniel B. Stephens and Associates, 6020 Academy NE, Suite 00, Albuquerque, NM (khsu@dbstephens.com) (Received July 4, 998; revised March 29, 999; accepted April, 999.)
Simple 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 informationDeterministic solution of stochastic groundwater flow equations by nonlocalfiniteelements A. Guadagnini^ & S.P. Neuman^
Deterministic solution of stochastic groundwater flow equations by nonlocalfiniteelements A. Guadagnini^ & S.P. Neuman^ Milano, Italy; ^Department of. Hydrology & Water Resources, The University ofarizona,
More informationAssessment of Hydraulic Conductivity Upscaling Techniques and. Associated Uncertainty
CMWRXVI Assessment of Hydraulic Conductivity Upscaling Techniques and Associated Uncertainty FARAG BOTROS,, 4, AHMED HASSAN 3, 4, AND GREG POHLL Division of Hydrologic Sciences, University of Nevada, Reno
More informationEnsemble Kalman filter assimilation of transient groundwater flow data via stochastic moment equations
Ensemble Kalman filter assimilation of transient groundwater flow data via stochastic moment equations Alberto Guadagnini (1,), Marco Panzeri (1), Monica Riva (1,), Shlomo P. Neuman () (1) Department of
More informationSecond-order expressions for velocity moments in twoand three-dimensional statistically anisotropic media
WATER RESOURCES RESEARCH VOL. 33 NO. 4 PAGES 625-637 APRIL 1997 Second-order expressions for velocity moments in twoand three-dimensional statistically anisotropic media Kuo-Chin Hsu and Shlomo P. Neuman
More informationHydraulic tomography: Development of a new aquifer test method
WATER RESOURCES RESEARCH, VOL. 36, NO. 8, PAGES 2095 2105, AUGUST 2000 Hydraulic tomography: Development of a new aquifer test method T.-C. Jim Yeh and Shuyun Liu Department of Hydrology and Water Resources,
More informationA note on the influence of a constant velocity boundary condition
WATER RESOURCES RESEARCH, VOL. 36, NO. 10, PAGES 3095-3101, OCTOBER 000 A note on the influence of a constant velocity boundary condition on flow and transport in heterogeneous formations S.C. Lessoff
More informationIn all of the following equations, is the coefficient of permeability in the x direction, and is the hydraulic head.
Groundwater Seepage 1 Groundwater Seepage Simplified Steady State Fluid Flow The finite element method can be used to model both steady state and transient groundwater flow, and it has been used to incorporate
More informationUse of the Fourier Laplace transform and of diagrammatical methods to interpret pumping tests in heterogeneous reservoirs
PII: S39-78(97)4-6 Advances in Water Resources, Vol. 2, pp. 58 59, 998 998 Elsevier Science Limited All rights reserved. Printed in Great Britain 39-78/98/$9. +. Use of the Fourier Laplace transform and
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 informationThe Stochastic Inverse Problem in Groundwater Hydrology: A Kalman Filter Approach
The Stochastic Inverse Problem in Groundwater Hydrology: A Kalman Filter Approach Franz Konecny Institute of Mathematics and Applied Statistics C/muerszZ?/ o/ Agrzcw/Zwm/ ^czences, Fzermo /W J#0 Ore^or
More informationInverse Modelling for Flow and Transport in Porous Media
Inverse Modelling for Flow and Transport in Porous Media Mauro Giudici 1 Dipartimento di Scienze della Terra, Sezione di Geofisica, Università degli Studi di Milano, Milano, Italy Lecture given at the
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 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 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 informationTechnical note: Analytical drawdown solution for steady-state pumping tests in two-dimensional isotropic heterogeneous aquifers
Hydrol. Earth Syst. Sci., 0, 655 667, 06 www.hydrol-earth-syst-sci.net/0/655/06/ doi:0.594/hess-0-655-06 Authors 06. CC Attribution 3.0 License. Technical note: Analytical drawdown solution for steady-state
More informationStochastic uncertainty analysis for solute transport in randomly heterogeneous media using a Karhunen-Loève-based moment equation approach
WATER RESOURCES RESEARCH, VOL. 43,, doi:10.1029/2006wr005193, 2007 Stochastic uncertainty analysis for solute transport in randomly heterogeneous media using a Karhunen-Loève-based moment equation approach
More informationType-curve estimation of statistical heterogeneity
WATER RESOURCES RESEARCH, VOL. 40,, doi:10.1029/2003wr002405, 2004 Type-curve estimation of statistical heterogeneity Shlomo P. Neuman Department of Hydrology and Water Resources, University of Arizona,
More informationA A + B. ra + A + 1. We now want to solve the Einstein equations in the following cases:
Lecture 29: Cosmology Cosmology Reading: Weinberg, Ch A metric tensor appropriate to infalling matter In general (see, eg, Weinberg, Ch ) we may write a spherically symmetric, time-dependent metric in
More informationExact transverse macro dispersion coefficients for transport in heterogeneous porous media
Stochastic Environmental Resea (2004) 18: 9 15 Ó Springer-Verlag 2004 DOI 10.1007/s00477-003-0160-6 ORIGINAL PAPER S. Attinger Æ M. Dentz Æ W. Kinzelbach Exact transverse macro dispersion coefficients
More informationUncertainty quantification for flow in highly heterogeneous porous media
695 Uncertainty quantification for flow in highly heterogeneous porous media D. Xiu and D.M. Tartakovsky a a Theoretical Division, Los Alamos National Laboratory, Mathematical Modeling and Analysis Group
More informationWATER RESOURCES RESEARCH, VOL. 24, NO. 12, PAGES , DECEMBER A Note on the Recent Natural Gradient Tracer Test at the Borden Site
WATER RESOURCES RESEARCH, VOL. 24, NO. 12, PAGES 2099-2103, DECEMBER 1988 A Note on the Recent Natural Gradient Tracer Test at the Borden Site R. L. N AFF U.S. Geological Survey, Denver, Colorado T.-C.
More informationFlow and travel time statistics in highly heterogeneous porous media
WATER RESOURCES RESEARCH, VOL. 45,, doi:1.129/28wr7168, 29 Flow and travel time statistics in highly heterogeneous porous media Hrvoje Gotovac, 1,2 Vladimir Cvetkovic, 1 and Roko Andricevic 2 Received
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 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 informationCharacterization of heterogeneous hydraulic conductivity field via Karhunen-Loève expansions and a measure-theoretic computational method
Characterization of heterogeneous hydraulic conductivity field via Karhunen-Loève expansions and a measure-theoretic computational method Jiachuan He University of Texas at Austin April 15, 2016 Jiachuan
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 informationGreen s Function for Tenuous Random Media
EECS 730, Winter 2009 c K. Sarabandi Green s Function for Tenuous Random Media In the previous section we demonstrated the application of wave theory to a tenuous medium excited by a plane wave for evaluating
More informationBasics in Geostatistics 2 Geostatistical interpolation/estimation: Kriging methods. Hans Wackernagel. MINES ParisTech.
Basics in Geostatistics 2 Geostatistical interpolation/estimation: Kriging methods Hans Wackernagel MINES ParisTech NERSC April 2013 http://hans.wackernagel.free.fr Basic concepts Geostatistics Hans Wackernagel
More informationpt)x+ x control volume
Appendix I DERIVATION OF EQUATIONS OF STEADY-STATE GROUNDWATER FLOW Consider a unit volume of saturated porous media (Figure AI. 1). In fluid mechanics, such a volume is called a control volume. The boundaries
More informationLegendre s Equation. PHYS Southern Illinois University. October 13, 2016
PHYS 500 - Southern Illinois University October 13, 2016 PHYS 500 - Southern Illinois University Legendre s Equation October 13, 2016 1 / 10 The Laplacian in Spherical Coordinates The Laplacian is given
More informationELEG 3143 Probability & Stochastic Process Ch. 6 Stochastic Process
Department of Electrical Engineering University of Arkansas ELEG 3143 Probability & Stochastic Process Ch. 6 Stochastic Process Dr. Jingxian Wu wuj@uark.edu OUTLINE 2 Definition of stochastic process (random
More informationGeostatistics for Seismic Data Integration in Earth Models
2003 Distinguished Instructor Short Course Distinguished Instructor Series, No. 6 sponsored by the Society of Exploration Geophysicists European Association of Geoscientists & Engineers SUB Gottingen 7
More informationStochastic representation of random positive-definite tensor-valued properties: application to 3D anisotropic permeability random fields
Sous la co-tutelle de : LABORATOIRE DE MODÉLISATION ET SIMULATION MULTI ÉCHELLE CNRS UPEC UNIVERSITÉ PARIS-EST CRÉTEIL UPEM UNIVERSITÉ PARIS-EST MARNE-LA-VALLÉE Stochastic representation of random positive-definite
More informationLecture 15. Theory of random processes Part III: Poisson random processes. Harrison H. Barrett University of Arizona
Lecture 15 Theory of random processes Part III: Poisson random processes Harrison H. Barrett University of Arizona 1 OUTLINE Poisson and independence Poisson and rarity; binomial selection Poisson point
More informationOne-Dimensional Steady State Infiltration in Heterogeneous Soils
WATER RESOURCES RESEARCH, VOL. 25, NO. 10, PAGES 2149-2158, OCTOBER 1989 One-Dimensional Steady State Infiltration in Heterogeneous Soils T.-C. JIM YEI Department of Hydrology and Water Resources, University
More informationObservable and measurable quantities
The Jeans Theorem The Jeans Theorem is an important result in stellar dynamics that states the importance of integrals of the motion in solving the CBE for gravitational potentials that do not change with
More informationOptical Imaging Chapter 5 Light Scattering
Optical Imaging Chapter 5 Light Scattering Gabriel Popescu University of Illinois at Urbana-Champaign Beckman Institute Quantitative Light Imaging Laboratory http://light.ece.uiuc.edu Principles of Optical
More informationEfficient geostatistical simulation for spatial uncertainty propagation
Efficient geostatistical simulation for spatial uncertainty propagation Stelios Liodakis University of the Aegean University Hill Mytilene, Greece stelioslio@geo.aegean.gr Phaedon Kyriakidis Cyprus University
More informationThis paper is part of the following report: UNCLASSIFIED
UNCLASSIFIED Defense Technical Information Center Compilation Part Notice ADPO 11594 TITLE: Strong-Property-Fluctuation Theory for Homogenization of Bianisotropic Composites DISTRIBUTION: Approved for
More informationGeostatistics in Hydrology: Kriging interpolation
Chapter Geostatistics in Hydrology: Kriging interpolation Hydrologic properties, such as rainfall, aquifer characteristics (porosity, hydraulic conductivity, transmissivity, storage coefficient, etc.),
More information18 Single vertical fractures
18 Single vertical fractures 18.1 Introduction If a well intersects a single vertical fracture, the aquifer s unsteady drawdown response to pumping differs significantly from that predicted by the Theis
More informationDiffusion of a density in a static fluid
Diffusion of a density in a static fluid u(x, y, z, t), density (M/L 3 ) of a substance (dye). Diffusion: motion of particles from places where the density is higher to places where it is lower, due to
More informationStochastic disaggregation of spatial-temporal rainfall with limited data
XXVI General Assembly of the European Geophysical Society Nice, France, 25-3 March 2 Session NP5.3/ Scaling, multifractals and nonlinearity in geophysics: Stochastic rainfall modelling: scaling and other
More informationReflection of SV- Waves from the Free Surface of a. Magneto-Thermoelastic Isotropic Elastic. Half-Space under Initial Stress
Mathematica Aeterna, Vol. 4, 4, no. 8, 877-93 Reflection of SV- Waves from the Free Surface of a Magneto-Thermoelastic Isotropic Elastic Half-Space under Initial Stress Rajneesh Kakar Faculty of Engineering
More informationPure Random process Pure Random Process or White Noise Process: is a random process {X t, t 0} which has: { σ 2 if k = 0 0 if k 0
MODULE 9: STATIONARY PROCESSES 7 Lecture 2 Autoregressive Processes 1 Moving Average Process Pure Random process Pure Random Process or White Noise Process: is a random process X t, t 0} which has: E[X
More informationElements of Multivariate Time Series Analysis
Gregory C. Reinsel Elements of Multivariate Time Series Analysis Second Edition With 14 Figures Springer Contents Preface to the Second Edition Preface to the First Edition vii ix 1. Vector Time Series
More informationWHITE NOISE APPROACH TO FEYNMAN INTEGRALS. Takeyuki Hida
J. Korean Math. Soc. 38 (21), No. 2, pp. 275 281 WHITE NOISE APPROACH TO FEYNMAN INTEGRALS Takeyuki Hida Abstract. The trajectory of a classical dynamics is detrmined by the least action principle. As
More informationB005 A NEW FAST FOURIER TRANSFORM ALGORITHM FOR FLUID FLOW SIMULATION
1 B5 A NEW FAST FOURIER TRANSFORM ALGORITHM FOR FLUID FLOW SIMULATION LUDOVIC RICARD, MICAËLE LE RAVALEC-DUPIN, BENOÎT NOETINGER AND YVES GUÉGUEN Institut Français du Pétrole, 1& 4 avenue Bois Préau, 92852
More information6. GRID DESIGN AND ACCURACY IN NUMERICAL SIMULATIONS OF VARIABLY SATURATED FLOW IN RANDOM MEDIA: REVIEW AND NUMERICAL ANALYSIS
Harter Dissertation - 1994-132 6. GRID DESIGN AND ACCURACY IN NUMERICAL SIMULATIONS OF VARIABLY SATURATED FLOW IN RANDOM MEDIA: REVIEW AND NUMERICAL ANALYSIS 6.1 Introduction Most numerical stochastic
More informationCollapse of Kaluza-Klein Bubbles. Abstract. Kaluza-Klein theory admits \bubble" congurations, in which the
UMDGR{94{089 gr{qc/940307 Collapse of Kaluza-Klein Bubbles Steven Corley and Ted Jacobson 2 Department of Physics, University of Maryland, College Park, MD 20742{4 Abstract Kaluza-Klein theory admits \bubble"
More informationA geostatistical inverse method for variably saturated flow in the vadose zone
WATER RESOURCES RESEARCH, VOL. 32, NO. 9, PAGES 2757 2766, SEPTEMBER 1996 A geostatistical inverse method for variably saturated flow in the vadose zone T.-C. Jim Yeh and Jinqi Zhang Department of Hydrology
More informationNumerical evaluation of solute dispersion and dilution in unsaturated heterogeneous media
WATER RESOURCES RESEARCH, VOL. 38, NO. 11, 1220, doi:10.1029/2001wr001262, 2002 Numerical evaluation of solute dispersion and dilution in unsaturated heterogeneous media Olaf A. Cirpka Universität Stuttgart,
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 informationAn Introduction to Partial Differential Equations in the Undergraduate Curriculum
An Introduction to Partial Differential Equations in the Undergraduate Curriculum John C. Polking LECTURE 12 Heat Transfer in the Ball 12.1. Outline of Lecture The problem The problem for radial temperatures
More informationThe Concept of Block-Effective Macrodispersion for Numerical Modeling of Contaminant Transport. Yoram Rubin
The Concept of Block-Effective Macrodispersion for Numerical Modeling of Contaminant Transport Yoram Rubin University of California at Berkeley Thanks to Alberto Bellin and Alison Lawrence Background Upscaling
More informationISSN: X (p); (e)
TORSIONA SURFACE WAVE IN SEF-REINFORCED AYER SANDWICHED BETWEEN TWO VISCO-EASIC HAF-SPACES UNDER THE INITIA STRESS Nidhi Dewangan, Sanjeev A. Sahu and Soniya Chaudhary S.G.G. Govt P.G. College, Kurud,
More informationModeling Seepage Control in Hydraulic Structures
Number 3 Volume 13 September2006 Journal of Engineering Dr. Rafa H. El-Suhaili Professor Environmental Engineering Department College of Engineering Fawaz M. Aziz Al-Joubori M.Sc Student Irrigation and
More informationMACRODISPERSION AND DISPERSIVE TRANSPORT BY UNSTEADY RIVER FLOW UNDER UNCERTAIN CONDITIONS
MACRODISPERSION AND DISPERSIVE TRANSPORT BY UNSTEADY RIVER FLOW UNDER UNCERTAIN CONDITIONS M.L. Kavvas and L.Liang UCD J.Amorocho Hydraulics Laboratory University of California, Davis, CA 95616, USA Uncertainties
More informationGaussian processes. Basic Properties VAG002-
Gaussian processes The class of Gaussian processes is one of the most widely used families of stochastic processes for modeling dependent data observed over time, or space, or time and space. The popularity
More informationPerformance Evaluation of Generalized Polynomial Chaos
Performance Evaluation of Generalized Polynomial Chaos Dongbin Xiu, Didier Lucor, C.-H. Su, and George Em Karniadakis 1 Division of Applied Mathematics, Brown University, Providence, RI 02912, USA, gk@dam.brown.edu
More informationPROBABILISTIC PREDICTIONS OF INFILTRATION INTO HETEROGENEOUS MEDIA WITH UNCERTAIN HYDRAULIC PARAMETERS
International Journal for Uncertainty Quantification, 1(1):35 47, 211 PROBABILISTIC PREDICTIONS OF INFILTRATION INTO HETEROGENEOUS MEDIA WITH UNCERTAIN HYDRAULIC PARAMETERS Peng Wang & Daniel M. Tartakovsky
More informationUpscaled flow and transport properties for heterogeneous unsaturated media
WATER RESOURCES RESEARCH, VOL. 38, NO. 5, 1053, 10.1029/2000WR000072, 2002 Upscaled flow and transport properties for heterogeneous unsaturated media Raziuddin Khaleel Fluor Federal Services, Richland,
More informationAN INTRODUCTION TO CURVILINEAR ORTHOGONAL COORDINATES
AN INTRODUCTION TO CURVILINEAR ORTHOGONAL COORDINATES Overview Throughout the first few weeks of the semester, we have studied vector calculus using almost exclusively the familiar Cartesian x,y,z coordinate
More informationGG655/CEE623 Groundwater Modeling. Aly I. El-Kadi
GG655/CEE63 Groundwater Modeling Model Theory Water Flow Aly I. El-Kadi Hydrogeology 1 Saline water in oceans = 97.% Ice caps and glaciers =.14% Groundwater = 0.61% Surface water = 0.009% Soil moisture
More informationBoundary value problems for partial differential equations
Boundary value problems for partial differential equations Henrik Schlichtkrull March 11, 213 1 Boundary value problem 2 1 Introduction This note contains a brief introduction to linear partial differential
More informationMOMENT DIFFERENTIAL EQUATIONS FOR FLOW IN HIGHLY HETEROGENEOUS POROUS MEDIA
MOMENT DIFFERENTIAL EQUATIONS FOR FLOW IN HIGHLY HETEROGENEOUS POROUS MEDIA C.L. WINTER Mathematical Modeling & Analysis Group and Center for Nonlinear Studies, Los Alamos National Laboratory, Los Alamos,
More informationBayesian methodology to stochastic capture zone determination: Conditioning on transmissivity measurements
WATER RESOURCES RESEARCH, VOL. 38, NO. 9, 1164, doi:10.1029/2001wr000950, 2002 Bayesian methodology to stochastic capture zone determination: Conditioning on transmissivity measurements L. Feyen, 1 P.
More informationA Stochastic Collocation based. for Data Assimilation
A Stochastic Collocation based Kalman Filter (SCKF) for Data Assimilation Lingzao Zeng and Dongxiao Zhang University of Southern California August 11, 2009 Los Angeles Outline Introduction SCKF Algorithm
More informationUpscaling of Richards equation for soil moisture dynamics to be utilized in mesoscale atmospheric models
Exchange Processes at the imnd Surface for a Range of Space and Time Scales (Proceedings of the Yokohama Symposium, July 1993). [AHS Publ. no. 212, 1993. 125 Upscaling of Richards equation for soil moisture
More informationMonte Carlo Simulation. CWR 6536 Stochastic Subsurface Hydrology
Monte Carlo Simulation CWR 6536 Stochastic Subsurface Hydrology Steps in Monte Carlo Simulation Create input sample space with known distribution, e.g. ensemble of all possible combinations of v, D, q,
More informationJordan Stoyanov Institute of Mathematics Bulgarian Academy of Sciences Sofia, Bulgaria
TEIE USE OF COUNTEIUWWPLES IN LEAIWING PROBABILITY AND STATISTICS Jordan Stoyanov Institute of Mathematics Bulgarian Academy of Sciences Sofia, Bulgaria Introductorv Notes Many of the lecturers here expressed
More information0 t b r 6, 20 t l nf r nt f th l t th t v t f th th lv, ntr t n t th l l l nd d p rt nt th t f ttr t n th p nt t th r f l nd d tr b t n. R v n n th r
n r t d n 20 22 0: T P bl D n, l d t z d http:.h th tr t. r pd l 0 t b r 6, 20 t l nf r nt f th l t th t v t f th th lv, ntr t n t th l l l nd d p rt nt th t f ttr t n th p nt t th r f l nd d tr b t n.
More informationP a g e 3 6 of R e p o r t P B 4 / 0 9
P a g e 3 6 of R e p o r t P B 4 / 0 9 p r o t e c t h um a n h e a l t h a n d p r o p e r t y fr om t h e d a n g e rs i n h e r e n t i n m i n i n g o p e r a t i o n s s u c h a s a q u a r r y. J
More informationMetrics and Curvature
Metrics and Curvature How to measure curvature? Metrics Euclidian/Minkowski Curved spaces General 4 dimensional space Cosmological principle Homogeneity and isotropy: evidence Robertson-Walker metrics
More informationConnection to Laplacian in spherical coordinates (Chapter 13)
Connection to Laplacian in spherical coordinates (Chapter 13) We might often encounter the Laplace equation and spherical coordinates might be the most convenient 2 u(r, θ, φ) = 0 We already saw in Chapter
More informationFundamentals of Noise
Fundamentals of Noise V.Vasudevan, Department of Electrical Engineering, Indian Institute of Technology Madras Noise in resistors Random voltage fluctuations across a resistor Mean square value in a frequency
More informationA full-bayesian approach to the groundwater inverse problem for steady state flow
WATER RESOURCES RESEARCH, VOL. 36, NO. 8, PAGES 2081 2093, AUGUST 2000 A full-bayesian approach to the groundwater inverse problem for steady state flow Allan D. Woodbury Department of Civil Engineering,
More information4 4 N v b r t, 20 xpr n f th ll f th p p l t n p pr d. H ndr d nd th nd f t v L th n n f th pr v n f V ln, r dn nd l r thr n nt pr n, h r th ff r d nd
n r t d n 20 20 0 : 0 T P bl D n, l d t z d http:.h th tr t. r pd l 4 4 N v b r t, 20 xpr n f th ll f th p p l t n p pr d. H ndr d nd th nd f t v L th n n f th pr v n f V ln, r dn nd l r thr n nt pr n,
More informationExpansion of 1/r potential in Legendre polynomials
Expansion of 1/r potential in Legendre polynomials In electrostatics and gravitation, we see scalar potentials of the form V = K d Take d = R r = R 2 2Rr cos θ + r 2 = R 1 2 r R cos θ + r R )2 Use h =
More informationLaplace s Equation in Cylindrical Coordinates and Bessel s Equation (I)
Laplace s Equation in Cylindrical Coordinates and Bessel s Equation I) 1 Solution by separation of variables Laplace s equation is a key equation in Mathematical Physics. Several phenomena involving scalar
More informationTwo-dimensional flow in a porous medium with general anisotropy
Two-dimensional flow in a porous medium with general anisotropy P.A. Tyvand & A.R.F. Storhaug Norwegian University of Life Sciences 143 Ås Norway peder.tyvand@umb.no 1 Darcy s law for flow in an isotropic
More informationA wavenumber approach to characterizing the diffuse field conditions in reverberation rooms
PROCEEDINGS of the 22 nd International Congress on Acoustics Isotropy and Diffuseness in Room Acoustics: Paper ICA2016-578 A wavenumber approach to characterizing the diffuse field conditions in reverberation
More information1 Reduced Mass Coordinates
Coulomb Potential Radial Wavefunctions R. M. Suter April 4, 205 Reduced Mass Coordinates In classical mechanics (and quantum) problems involving several particles, it is convenient to separate the motion
More informationA note on Gaussian distributions in R n
Proc. Indian Acad. Sci. (Math. Sci. Vol., No. 4, November 0, pp. 635 644. c Indian Academy of Sciences A note on Gaussian distributions in R n B G MANJUNATH and K R PARTHASARATHY Indian Statistical Institute,
More informationProcess of whiskey maturation & Current Efforts to improve Comparison Goals Modeling nonlinear Diffusion
Brian England Objective Process of whiskey maturation & Current Efforts to improve Comparison Goals Modeling nonlinear Diffusion Cylindrical Coordinates Initial and Boundary Conditions Methodologies and
More informationHomogenization Theory
Homogenization Theory Sabine Attinger Lecture: Homogenization Tuesday Wednesday Thursday August 15 August 16 August 17 Lecture Block 1 Motivation Basic Ideas Elliptic Equations Calculation of Effective
More informationHomogenization and numerical Upscaling. Unsaturated flow and two-phase flow
Homogenization and numerical Upscaling Unsaturated flow and two-phase flow Insa Neuweiler Institute of Hydromechanics, University of Stuttgart Outline Block 1: Introduction and Repetition Homogenization
More informationPROCEEDINGS. The present investigation will show that these two line elements, together
PROCEEDINGS OF THE NATIONAL ACADEMY OF SCIENCES Volume 15 April 15, 1929 Number 4 ON THE POSSIBLE LINE ELEMENTS FOR THE UNIVERSE By RICHARD C. TOLMAN NoRMAN BRIDGs LABORATORY OF PHYSICS, CIFORNLIA INS~flTrTr
More informationNote the diverse scales of eddy motion and self-similar appearance at different lengthscales of the turbulence in this water jet. Only eddies of size
L Note the diverse scales of eddy motion and self-similar appearance at different lengthscales of the turbulence in this water jet. Only eddies of size 0.01L or smaller are subject to substantial viscous
More informationMeasure-Theoretic parameter estimation and prediction for contaminant transport and coastal ocean modeling
Measure-Theoretic parameter estimation and prediction for contaminant transport and coastal ocean modeling Steven Mattis and Lindley Graham RMSWUQ 2015 July 17, 2015 Collaborators 2/45 The University of
More informationGaussian Random Fields
Gaussian Random Fields March 22, 2007 Random Fields A N dimensional random field is a set of random variables Y (x), x R N, which has a collection of distribution functions F (Y (x ) y,..., Y (x n ) y
More informationAppendix: Orthogonal Curvilinear Coordinates. We define the infinitesimal spatial displacement vector dx in a given orthogonal coordinate system with
Appendix: Orthogonal Curvilinear Coordinates Notes: Most of the material presented in this chapter is taken from Anupam G (Classical Electromagnetism in a Nutshell 2012 (Princeton: New Jersey)) Chap 2
More informationJournal of Hydrology
Journal of Hydrology 387 (21) 32 327 Contents lists available at ScienceDirect Journal of Hydrology journal homepage: www.elsevier.com/locate/jhydrol Steady-state discharge into tunnels in formations with
More informationResearch Article The Laplace Likelihood Ratio Test for Heteroscedasticity
International Mathematics and Mathematical Sciences Volume 2011, Article ID 249564, 7 pages doi:10.1155/2011/249564 Research Article The Laplace Likelihood Ratio Test for Heteroscedasticity J. Martin van
More informationIntroduction to Spatial Data and Models
Introduction to Spatial Data and Models Sudipto Banerjee 1 and Andrew O. Finley 2 1 Department of Forestry & Department of Geography, Michigan State University, Lansing Michigan, U.S.A. 2 Biostatistics,
More informationPHYS 211 Lecture 9 - Examples of 3D motion 9-1
PHYS 211 Lecture 9 - Examples of 3D motion 9-1 Lecture 9 - Examples of 3D motion Text: Fowles and Cassiday, Chap. 4 In one dimension, the equations of motion to be solved are functions of only one position
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 informationPoisson line processes. C. Lantuéjoul MinesParisTech
Poisson line processes C. Lantuéjoul MinesParisTech christian.lantuejoul@mines-paristech.fr Bertrand paradox A problem of geometrical probability A line is thrown at random on a circle. What is the probability
More informationMonte Carlo analysis of macro dispersion in 3D heterogeneous porous media
Monte Carlo analysis of macro dispersion in 3D heterogeneous porous media Arthur Dartois and Anthony Beaudoin Institute P, University of Poitiers, France NM2PourousMedia, Dubrovnik, Croatia 29 Sep - 3
More information