STOCHASTIC FINITE ELEMENT ANALYSIS OF GROUNDWATER FLOW USING THE FIRST- ORDER RELIABILITY METHOD
|
|
- Rose Matthews
- 5 years ago
- Views:
Transcription
1 Abstract STOCHASTIC FINITE ELEMENT ANALYSIS OF GROUNDWATER FLOW USING THE FIRST- ORDER RELIABILITY METHOD Jeffrey D. Cawlfield Dept. of Geological Engineering University of Missouri-Rolla Rolla, MO and Nicholas Sitar Dept. of Civil Engineering University of California, Berkeley Berkeley, CA An approach for stochastic analysis of groundwater flow based on the first-order reliability method is presented. The method can fully utilize any level of probabilistic information from the minimum of second-moments to the complete full joint distribution. It is well suited for problems in which the statistical information is incomplete, as is common when considering groundwater flow. The results of a first-order reliability analysis include an estimate of the probability of exceeding a specified performance criteria and the sensitivity of the stochastic solution to changes in the uncertain variables and their moments. First-order reliability results from an analysis of onedimensional flow with uncertain hydraulic conductivity compare well with results from a previously presented Monte Carlo analysis. In our two-dimensional stochastic finite element model, the uncertain variables are the element hydraulic conductivities and the values of constant head or point flux boundary conditions. Example two-dimensional applications are presented for confined flow between two constant head boundaries. The effect of spatial correlation between element hydraulic conductivities on the probability estimate and the sensitivity information is investigated. Introduction In recent years much attention has been given to the incorporation of uncertainty into the analysis of subsurface flow in order to explicitly account for porous media variability. Sources of uncertainty which arise during analysis of subsurface flow were listed by Tang -191-
2 [1984] as follows: inherent spatial variability, random measurement errors, systematic sampling bias, calibration errors and estimation errors from insufficient sampling. In addition, the solution model may contribute uncertainty due to model simplifications and assumptions. The analyst must attempt to account for each source of uncertainty, and a solution model should incorporate as much of the resulting statistical information as possible. In this paper we assume that the statistical characteristics of the uncertain parameters are available from an appropriate analysis and we present a method for utilizing that statistical information in a stochastic model using the first-order reliability method. The first-order reliability method is an attractive stochastic modeling technique for a number of reasons: (1) as a minimum, the method requires only second-moment information concerning the uncertain variables; however, additional information (such as marginal distributions) can be incorporated into the solution, (2) the method directly yields an estimate of probability associated with a particular uncertain event; for example, exceedance of fluid head at a given location in a flow domain, (3) the first-order reliability calculations directly yield sensitivity information, (4) the method can be used with numerical solution models, (5) the method is efficient even for very low probability events. After a brief review of previously used stochastic models, we present the theoretical background of the first-order reliability method and then show the basic implementation of the method to specific groundwater flow problems. In addition, example one-dimensional analytical flow solutions and two-dimensional stochastic finite element analyses are presented. A more extensive treatment can be found in Sitar et al. [1987] and Cawlfield and Sitar [1987]. Previous Stochastic Modeling Early groundwater flow models which incorporated uncertainty assumed statistical independence between hydraulic conductivity values (e.g., Warren and Price [1961], McMillan [1966] and Freeze [1975]). However, hydraulic conductivity (or transmissivity) usually exhibits spatial correlation over some distance. Gelhar [1976] incorporated spatial correlation through an application of stochastic continuum theory, where the spatial covariance function of hydraulic conductivity was related to the covariance function of the head random field through spectral analysis. The stochastic continuum approach relies upon small perturbations of the random variables and, in order to obtain closed-form functional relationships between hydraulic conductivity uncertainty and head uncertainty, the method is usually only applied to tabular aquifers of infinite extent. Additional work using the stochastic continuum approach has been carried out using spectral analysis (e.g, Bakr et al. [1978], Gutjhar et al. [1978], Gutjahr and Gelhar [1981] and Mizel at al. [1982]). Other solution methods have also been used to analyze the stochastic continuum (e.g., a Green's function approach by Dagan [1982], and classical small perturbation theory by Tang and Pinder [1977]). Smith and Freeze [1979a and 1979b] used a nearest-neighbor technique to model spatial correlation in conjunction with Monte Carlo simulation. Their approach accounted for more complex boundary conditions and problem geometry than were considered with -192-
3 stochastic continuum theory. Others have used geostatistical techniques to incorporate spatial variability into Monte Carlo analyses (e.g., Delhomme [1979], Clifton and Neuman [1982] and Clifton [1985]). These previous methods of stochastic modeling have concentrated on analysis of the uncertain moments, usually the mean and variance of the random head. Although these moment calculations are important, decisions are usually based on a quantitative evaluation of risk. For use in formal decision analysis, results must be expressed as a measure of probability, and the first-order reliability method directly yields such measures by estimating the probability associated with a particular flow event. The method also directly provides information concerning the sensitivity of the stochastic solution to changes in the uncertain variables and their moments. Therefore, the first-order reliability method provides an attractive alternative to the available stochastic continuum and Monte Carlo techniques. First-Order Reliability Theory Introduction The first-order reliability method has been used in structural stochastic analysis for some time. The method can, in principle, incorporate any amount of probabilistic information, from the minimum requirement of second moments to the full joint probability distribution function (PDF). In groundwater flow problems the statistical information is typically incomplete, but often the marginal distributions are known or can be reasonably assumed (i.e., hydraulic conductivity is commonly assumed to be lognormally distributed). Therefore, the first-order reliability method is particularly well suited for stochastic groundwater flow analysis. The method is based on two main concepts which will be discussed in some detail: (1) formulation of a performance function which describes the flow behavior of interest in terms of the random variables, and (2) transformation of the problem into standard normal space, where an estimate of probability is obtained. The performance function is symbolized as Z=g(X), where X is a vector of the n uncertain parameters in a problem. For example, X may contain hydraulic conductivities, boundary conditions and problem geometry. For convenience, the limit state is defined as Z=g(X) = 0. The limit state is, in general, an n-dimensional hypersurface which divides the performance into two regions: a region where Z=g(X) < 0, and a region where Z=g(X) > 0. We will refer to the region where Z <. 0 as the region of exceedance, and the objective is to estimate the probability of exceeding an arbitrary target value (symbolized as P[Z <_ 0]). To illustrate these concepts, consider two-dimensional space and a function q = f(x X 2 ). The objective is to evaluate the probability that q is greater than some value of interest, say q t. The performance function is formulated as follows: Z = g(x) = q t - q = q t - f(x X 2 ) Figure la illustrates this performance function, showing the limit state surface as a line for this case. The probability that q J> q t is equivalent to P[Z <. 0], which is obtained by integrating the joint PDF in the region where g(x) <. 0 as follows: P[Z <. 0] = F Z (0) = / f X (x)dx -193-
4 where Fz(0) is the cumulative distribution function (CDF) for Z evaluated at 0, fx(x)dx is the joint PDF of X and the integral is taken over the region where g(x) <. 0. Figure 2 gives a graphical interpretation of Z and its relation to the limit state. Z is a margin Figure la. Two-dimensional example performance function, Z = g(x) = q t - f(x X 2 )
5 G(y)>0 HYPERPLANE^ TANGENT AT y* Figure lb. Standard space representation of the performance function of Figure la, illustrating the first-order approximation for P[Z <_ 0]. /3, a and y* are the reliability index, vector of direction cosines and design point as described in the text
6 Figure 2. Illustration of the margin against exceedance, Z, and probability of exceedance, p. fz(z) is the (unknown) PDF of Z
7 against exceedance (i.e., a safety margin), and the probability content which we seek to evaluate is shown cross-hatched in Figure 2. The exact evaluation of P[Z <_ 0] requires an n-fold integration of the joint PDF. This exact evaluation is difficult for practical problems because the joint PDF is unknown, the n-fold integration is over a generally complex region where Z = g(x) <. 0 and the integration is virtually impossible for n greater than about 2 or 3. Reliability Index Because of the difficulties in evaluating the n-fold probability integral, approximate methods for evaluating the probability content defined by Z <. 0 have been developed based on a measure of reliability known as the reliability index. Within the context of second moment information, Cornell [1969] defined the reliability index Q3) as follows: P = Mz / z where /^ andct Z are the mean and standard deviation of the performance function, Z, respectively. Figure 2 illustrates the concept of Cornell's reliability index, which is shown to represent the distance between the mean performance and the limit state (measured in units of the standard deviation of Z). Therefore, /3 can be thought of as a measure of the margin against exceedance (or safety margin), and a high /3 indicates a low probability of exceedance. If the performance function is linear, the moments JJL^ and o~ z may be computed exactly in terms of the means and covariances of the random variables X. However, when g(x) is a nonlinear function of X, the full joint PDF is required to exactly evaluate ^ and CT Z. Because the full joint PDF is seldom known, early approaches used a first-order mean centered Taylor series approximation as follows [Ang and Cornell, 1974]: Mz *> g(m) a z 2 = 2 S ôg/3xj dg/ôxj pjj CTJCTJ i J where the double summation is over all combinations of i and j, and the partial derivatives are evaluated at the point M (the vector of mean values of X). The approximation requires only second moment information concerning the X random vari- ables, and the performance function must be differentiable with respect to each X variable. The reliability index obtained using these approximations is known as a mean value first-order second moment reliability index (/?MVFOSM)- The MVFOSM index has been shown to have a major shortcoming: the value of J^MVFOSM ' s not invariant with respect to mutually consistent formulations of Z = g(x) [Hasofer and Lind, 1974]. Consider, for example, a performance function Z = qt-x,x 2, which could be equivalently formulated as Z = (qt / X,)-X 2. In both cases, Z will be zero, less than zero and greater than zero for the same values of X X 2 and qj, so the performance functions are mutually consistent with respect to the limit state. However, the value of /^MVFOSM ' S not necessarily equal for these two performance functions. This result is unsatisfactory, because consistency requires that both the linearization point and a unique performance function must then be specified for each problem. Hasofer and Lind [1974] showed that the invariance problem can be overcome by linearization on the limit state surface (rather than at the mean point), because mutually consistent performance functions all have identical surfaces defined by Z = 0. The conceptualization of the reliability index was illustrated in Figure 2 as a measure of the distance between mean performance and the limit state. Following that logic, Hasofer and Lind [1974] defined an invariant reliability index in terms of the minimum distance -197-
8 between the limit state and the mean point in standard space. The standard space consists of uncorrelated variâtes with standard deviations equal to one and means equal to zero (therefore, the mean point (M) in standard space is conveniently the origin). Figure lb is the standard space representation of the original space shown in Figure la. The standard space exhibits rotational symmetry; that is (refer to Figure lb), for all hyperplanes of equal distance from the origin, the probability content P is constant within a second moment representation (i.e., upper Tchebycheff bound), regardless of the orientation of the hyperplane [Der Kiureghian and Liu, 1986]. Evaluation of the reliability index using the Hasofer-Lind definition requires that the original random variables (X) be transformed into uncorrelated standard variâtes (Y). This transformation is conveniently expressed as Y = L" 1 D" 1 (X - M) where D is a diagonal matrix of standard deviations of X, L is a lower-triangular decomposition of the correlation matrix R = [ Pij ] and M is the vector of means of X. The transformation matrix L" 1 may be evaluated through Cholesky decomposition. The point on the limit-state surface nearest to the origin in standard space is known as the design point (y*), and its counterpart in the original space is symbolized as x*. The reliability index, /3, can be expressed in terms of the design point as fi = a* y* where a* is the unit normal at the design point directed towards the region where Z <_ 0. An appropriate optimization algorithm for determining the design point is given by Madsen et al. [1986], as Yi+i = [<*m + G(yi)/ V y G(y i ) ]a i T here g(x)= g(dly + M) = G(Y), V y G(y) signifies the gradient of G(y) with respect to the vector yj and V v G(y) signifies the magnitude of the gradient vector. The design point (y*) can be thought of as the most likely set of circumstances leading to exceedance. In standard normal space the contours of probability are concentric spheres, and the probability decays with distance from the origin. Therefore, the design point, which is the closest point to the origin on the limit state surface, has the highest likelihood of occurrence of any point in the region where Z <_ 0. Without complete joint PDF information, the actual form of probability decay is unknown and, strictly speaking, the design point is not necessarily the most likely set of circumstances leading to exceedance. However, such an interpretation is plausible as a general concept. Probability Estimate With only second moment information, in a strict sense, P[Z <. 0] cannot be estimated. However, an ad hoc approximation given by P[Z <. 0] = *(- /3) is often used, where *( ) is the standard cumulative normal probability function. This approximation assumes that Z is normally distributed. Now consider the case where the joint PDF is known. Then it is convenient to -198-
9 transform the X random variables into standard normal variâtes Y. The transformation Y = T(X) is nonlinear for nonnormal X, and is given in terms of the joint PDF of X [Hohenbichler and Rackwitz, 1981; Der Kiureghian and Liu, 1986]. Standard normal space possesses rotational symmetry in a stricter sense than was mentioned previously for standard space. Referring to Figure lb, the probability P defined by the hyperplane is given exactly by *(- js). Therefore, for a general limit state surface, a first-order approximation of P[Z <. 0] is defined by P[Z < 0] - *(- j8) This approximation is obtained by replacing the actual limit state surface with a tangent hyperplane at the design point y*. The approximation works well because, in standard normal space, the probability density decays exponentially as a function of distance squared from the origin; therefore, most of the probability content in the region Z <. 0 is concentrated in a small area near the design point, where the tangent hyperplane is a very good approximation to the limit state surface (Figure lb). The approximation is good in most practical problems, but more refined approximations that account for the curvature of the limit state surface near the design point are discussed by Madsen et al. [1986] and Der Kiureghian et al. [1987]. In groundwater flow problems, the analyst usually operates with information that falls between the two extremes just discussed; that is, between only second moment information (where the probability estimate is ad hoc) and full joint distribution information (where the properties of standard normal space can be used to obtain an exact first-order estimate of probability). Usually, the groundwater analyst will be able to assign marginal distributions to most, if not all, the uncertain variables in the problem (i.e., lognormal hydraulic conductivities). For such problems with incomplete probability information, Der Kiureghian and Liu [1986] have developed a general joint distribution model which is consistent with all available statistical information and which greatly simplifies the transformation to standard normal space. The joint PDF according to this model is given by fxw = 4> n (z,r0) [fxiw- fxn(x n )l / [*W - <P^n)] where zj = *"'(Fxi(xj)), < ( ) is the standard normal PDF and * n (z, RQ) is the n-dimensional normal density of zero means, unit standard deviations and correlation matrix RQ- The correlation matrix RQ is obtained using the integral relation [Der Kiureghian and Liu, 1986]: Pij = / / [ (xi-mi) / o-j] [(xj-mj) / o-j] 0 2 (z;,zj,poij)dzidzj where the double integration is carried out from - o to + oo. For each pair of variables Xj and X; with known marginal distributions and correlation coefficient pjj, the integral relation can be solved for the coefficient pqij, which is an element of Roij- Liu and Der Kiureghian [1986] have developed semi-empirical formulas which allow convenient calculation of selected marginal distribution types. The transformation to the standard normal space for the general joint distribution model is given by -199-
10 Y = LQ" 1 ^(FxiOq)) *" 1 (F X n(%)) where LQ is the lower-triangular decomposition of RQ. This transformation is convenient, because it requires only marginal distributions for each X. The reliability index is evaluated as before, by iteratively solving for the minimum distance between the limit state surface and the origin in standard space. If marginal distribution information is known, the general model is used and the transformation is carried out at each trial linearization point. The reliability index resulting from the calculations can then be used to give a first-order estimate of P[Z <. 0]. The accuracy of the estimate is a function of the curvature of the limit-state surface near the design point. In addition, the probability estimate assumes that the general joint PDF model is appropriate for the problem. Sensitivity Information As mentioned previously, the first-order reliability method directly yields sensitivity information. The most elementary measure of sensitivity is the partial derivative of 0 with respect to the values of the random variables at the design point (y*). Such a sensitivity measure gives an indication of how /3 will change when a particular random variable is slightly perturbed about its design point value. The sensitivity measure is given by a* = y*/1 y* = V y * j8 where V v */3 signifies a vector of gradients of /3 with respect to the components of the vector y*. The values of a* result directly from solving for the design point. For practical applications, a sensitivity measure in the space of the original random variables (X) is required. Using the chain rule of differentiation, Der Kiureghian and Liu [1985] suggest a scaled sensitivity measure in original space given by 7 = (V x */3) D / (V X *J3)D The gamma vector indicates the relative importance of each X variable, and the use of the standard deviation matrix D as a scaling factor provides a measure of the sensitivity of /3 with respect to equally likely changes in the random variables X. Sensitivity vectors can also be derived to measure sensitivity with respect to deterministic parameters in the performance function. If the performance function is given by g(x) = g(x, -n) where T) is a vector of deterministic parameters, Madsen et al. [1986] have shown that \ P = (V-n (g(x*,ti)) / V y G(y*) In groundwater flow applications, where either a target head or target flux is usually used as a deterministic parameter in the performance function, this formulation can be used to measure the sensitivity of )8 to the target value. Madsen et al. [1986] also give a formulation for evaluating the sensitivity of /3 with respect to distribution parameters for the random variables X ( i.e., means and variances). Suppose fx(x) = fx(x, 6), where 0 is a vector of distribution parameters. Then -200-
11 V 0 js = a* J T0 where JJQ is a matrix containing the partial derivatives of the transformation T with respect to the parameters 6, evaluated at the design point. Example Applications One Dimensional Problems With Analytical Flow Solutions Sitar et al. [1987] used the first-order reliability method to evaluate the uncertain head in a one-dimensional flow region; the problem had been previously analyzed by Freeze [1975] using Monte Carlo simulation. The flow region is shown in Figure 3, and it is divided into ten equally dimensioned blocks which are assumed to have hydraulic conductivities (K) described by independent lognormal random variables with means of Y equal to 1.0 logjo cm/sec and standard deviations of Y equal to 0.5 log^o cm/sec (where Y=logioK). An analytical solution for the fluid potential ($) in any section of the flow region was given by Freeze [1975]. For comparison with the Monte Carlo results of Freeze [1975], the first-order reliability performance function is formulated to compare a target fluid potential at a specific location in the flow region with the fluid potential calculated using the analytical solution (cp), as follows: z = g(x) = <p t -4> An estimate of P[Z <_ 0] = P[ <p > 4>t ] * s obtained using the first-order reliability analysis, and a PDF for the fluid potential can be developed from the individual results over a large range of target values. Consider the PDF for <p at x=45.0 cm in the flow region. Figure 4 compares the firstorder reliability PDF to a normal curve fit to Monte Carlo results of Freeze [1975]. The curves agree quite well, although there is some deviation at either end of the PDF. This discrepancy is an artifact of fitting a normal curve to the Monte Carlo results, whereas the curve fit to the first-order reliability results simply connects the individually calculated PDF points. Had a normal curve been fit to the first-order reliability PDF points, the two curves would have been practically identical. Figure 5 shows a comparison of PDF forms for <p near the constant head boundary at x=85.0 cm in the flow region (near the high potential boundary), where a skewed PDF is obtained. In this case, the two curves are nearly identical. We should note that using the first-order reliability technique to develop a complete PDF form is somewhat inefficient, because the first-order reliability analysis yields only one estimate of probability for each calculation. However, for an analysis of a specific probability associated with one (or a few) target values, the first-order reliability approach will usually be much more efficient than Monte Carlo simulation, particularly for low probability events. Additionally, the first-order reliability evaluation provides sensitivity information, which is inefficient to obtain using Monte Carlo simulation. Stochastic Finite Element Analysis of Two-Dimensional Flow We have recently implemented the first-order reliability method in conjunction with a finite element model of two-dimensional steady state groundwater flow [Cawlfield and Sitar, 1987]. The performance function is formulated either in terms of nodal heads (hjj or flux at constant head nodes (qj) as follows: -201-
12 x=o s y > i=i K. i=2 K 2 i=3 K 3 i=m^ K H X=L i=m *L K m Figure 3. One-dimensional flow region of Freeze [1975]. L = cm
13 First-Order Reliability Monte Carlo FLUID POTENTIAL (cm) 100 Figure 4. Comparison of First-Order Reliability PDF and Monte Carlo PDF for <p at x 45.0 cm
14 0.05 i i PDF 0.0Z First -Order ~Monte Carlo Reliab Llity /\\ A i k h *<t> m L FLUID POTENTIAL (cm) 100 Figure 5. Comparison of First-Order Reliability PDF and Monte Carlo PDF for <p at x 85.0 cm
15 g(x) = h t + 2 a k h k or g(x) = q t + S a^ where h t and q t are the target values, a k or a\ are deterministic œnstants and the summations are over all values k=l to m or 1=1 to p (m is the number of nodes and p is the number of constant head nodes). The solution then requires evaluation of the partial derivatives 3g/dXj for all random variables Xj. These can be explicitly calculated by direct differentiation of the numerical equations resulting from the finite element discretization. For confined steady state flow, the finite element discretization results in the following general matrix equation: Kh = q where K is the global conductance matrix, h is the vector of nodal heads and q is the vector of fluxes. If the random variables of concern are the element hydraulic conductivities, direct differentiation of the two types of performance functions given previously leads to the following: dg/dkj = -a T K" 1 (ÔK/dKi) h for the first type of performance function and dg/diq = a T 3Q/dKi for the second type, where each entry of the ôq/dk; vector is given by octy 3Kj = SKjj / ÔK; h + Ky (-K^ÔK / dkjh) In the last equation, Kj; is a row vector with j taking on values from 1 to n and n is the number of nodes. IC* and ôk/dkj can be explicitly calculated from the global conductance matrix. The inverse matrix must be calculated only once during each iteration. The ôk/dk{ matrix must be calculated once for each uncertain Kj during each iteration, but the computational requirements are relatively minor because only a few terms in the conductance matrix are involved. In a similar fashion, explicit calculations of the partial derivative of the performance function with respect to other random variables can be formulated (for example, the vector of random variables (X) might include constant head or flux boundary values). The full derivations are not given here, but it should be apparent that explicit calculation of these partial derivatives increases the efficiency of the computations, as opposed to using, for example, a finite divided difference approximation of the derivatives. In our model, each element hydraulic conductivity, values of constant boundary heads and values of externally applied point fluxes may be considered random variables. Problem geometry is assumed to be deterministic. Any of the random variables may be correlated to each other, and the model is capable of generating correlation coefficients between element hydraulic conductivities, based on an assumed autocorrelation functional form. The autocorrelation functional form may be assigned as a point function, in which case the element properties are derived using the spatial averaging techniques of Vanmarke [1983]. In the examples discussed in this section, we have assumed that the point hydraulic conductivity spatial correlation can be described by an exponential decay function, and the x and y direction spatial correlation functions are independent (separa
16 ble). Consider the finite element mesh shown in Figure 6, containing 66 nodes and 50 elements, with confined flow between two constant head boundaries. Suppose we are interested in estimating the probability that the head at node 33 ($33) is greater than some target head (< t ). The performance function is then formulated as Z=g(X)=0 t Assume the constant head boundaries are deterministic and fixed, so that the only uncertain parameters are the element hydraulic conductivities, which are lognormally distributed and have the following point statistics: yx K = and org; = In addition, the K field is isotropic, with isotropic exponential autocorrelation (i.e., the x and y autocorrelation functional forms are identical) as follows: p = e (-1 d I / 0.58) where d is the distance of separation between two K values and 8 is the scale of fluctuation. The scale of fluctuation is equal to two times the correlation integral scale, which is also often used to express correlation strength. The steady state head at node 33 for a deterministic analysis of confined flow in a uniform porous media would be Now consider two target heads at node 33 for our stochastic analysis: the relatively likely event of <p=110.0, and the much lower probability event of <p= The probabilities estimated using the first-order reliability method and the given statistical information are shown for these two cases in Figures 7 and 8 as a function of scale of fluctuation. First note that P[<p > ] is much higher than P[(p 140.0], as expected. Perhaps the most important result illustrated by this example is the sensitivity of the probability estimate to the scale of fluctuation. It appears that, for this type of flow problem, the scale of fluctuation is extremely important with respect to probability estimates, but only if the scale of fluctuation is less than about one-half the distance between the two constant head boundaries. The gamma sensitivity results are shown in Figures 9 and 10 for the two target heads. The sensitivity results are not surprising, as it is expected that the stochastic solution would be most sensitive to the K values in elements immediately surrounding node 33. For the low probability event (<p l = 140.0), the solution is very sensitive to the K values in the elements immediately to the left of node 33, and also surprisingly sensitive to the element K values along the right constant head boundary. To demonstrate potential applications to more complex problems, consider a layered flow region with a low conductivity material surrounded by a high conductivity material, such as shown in Figure 11. Both materials have coefficients of variation equal to 1.0, and autocorrelation occurs only between element K values in the same material (with isotropic exponential autocorrelation as before). If the same performance function is used to evaluate probabilities associated with the head at node 33, we find that the probability estimates are not affected too much by the introduction of the low conductivity layers (i.e., the probability estimates are within about 10% of the estimates for the one material case considered previously). However, the sensitivity information is much different for the two material case, as shown for the gamma sensitivity in Figure 12. Note that the solution is now most sensitive to values of K in the elements along the bottom row of the mesh (in the high conductivity zone closest to node 33). These results illustrate the potential utility of the first-order reliability method for more complex analyses
17 NO FLOW BOUNDARY = = Q NO FLOW BOUNDARY Figure 6. Finite element mesh for example stochastic finite element analysis of confined flow region
18 o d O l-oo 6 < Ld cc CD Q < LU MA- =OA- = , =6y O o- oo SCALE OF FLUCTUATION Figure 7. P[033 > 110.0] at node 33 as a function of scale of fluctuation (S)
19 O" o O" < to o* ^=0,. = , =6, < o^ SCALE OF FLUCTUATION Figure 8. P[</>33 > 140.0] at node 33 as a function of scale of fluctuation (8)
20 GAMMA SENSITIVITIES , J Figure 9. Gamma sensitivity for P[#33 > 110.0]
21 GAMMA SENSITIVITIES Figure 10. Gamma sensitivity for P[#33 > 140.0]
22 .»» «' t. S ff * a ««* * H Mat. 1: * * '.ME it. 2: K X =K 0 1 K x -K yy \L K -= lxlc ^ 6 x=8 y = * ". * *» * * %.,>!% = t * -6, 1x10 a * v -.. *. * *. ««««a o Figure 11. Two-dimensional confined flow region with two materials
23 GAMMA SENSITIVITIES , f Figure 12. Gamma sensitivity for P[033 > 110.0] for two material case
24 Conclusions The first-order reliability method provides an attractive technique for stochastic analysis because it is general and flexible, it is capable of incorporating a wide range of statistical information (including marginal PDF forms), it provides an estimate of probability associated with a particular subsurface flow event and it directly yields information concerning the sensitivity of the stochastic solution to changes in the random variables and their moments. The method can be used with either analytical or numerical solutions of the direct flow problem, and the performance function can be generally formulated such that the probability associated with a wide variety of flow events may be estimated. Results from application of the first-order reliability analysis to one-dimensional problems with analytical flow solutions indicate that the first-order reliability estimates of probability compare very well to Monte Carlo simulation estimates. Our results using a stochastic finite element model of two-dimensional steady state flow illustrate the potential utility of the first-order reliability method in the analysis and interpretation of more complex groundwater flow problems. ACKNOWLEDGMENTS The work presented in this report has been supported in part by the National Science Foundation under PYI Award No. CEE The opinions, findings and conclusions expressed in this report are those of the authors and do not necessarily reflect the views of the National Science Foundation. The authors wish to express their appreciation to Professors Armen Der Kiureghian and T. N. Narasimhan for their helpful review and suggestions throughout the development of this work
25 REFERENCES Ang, Alfredo H.-S. & Cornell, A. C. (1974) Reliability bases of structural safety and design. Journal of the Structural Division, ASCE 100(ST9): Bakr, A. A., Gelhar, L. W., Gutjahr, A. L. & MacMillan, J. R. (1978) Stochastic Analysis of Spatial Variability in Subsurface Flow, 1. Comparison of One- and Three-Dimensional Flows. Water Resour. Res. 14(2): Cawlfield, Jeffrey D. & Sitar, Nicholas (1987) Application of First-Order Reliability to Stochastic Finite Element Analysis of Groundwater Flow. Report No. UCB/GT/87-01, Department of Civil Engineering, University of California, Berkeley. Clifton, Peter M. (1985) Groundwater Travel Time Uncertainty Analysis-Sensitivity of Results to Model Geometry, and Correlations and Cross-Correlations Among Input Parameters. RHO-BW-ST-70-P, Rockwell Hanford Operations, Richland, Washington. Clifton, Peter M. & Newman, S. P. (1982) Effects of Kriging and Inverse Modeling on Conditional Simulation of the Avra Valley Aquifer in Southern Arizona. Water Resour. Res. 18(4): Cornell, G A. (1969) A Probability Based Structural Code. J. of American Concrete Institute 66(12): Dagan, G. (1982) Stochastic Modeling of Groundwater Flow by Unconditional and Conditional Probabilities, 1. Conditional Simulation and the Direct Problem. Water Resour. Res. 18(4): Delhomme, J. P. (1979) Spatial Variability and Uncertainty in Groundwater Flow Parameters: A Geostatistical Approach. Water Resour. Res. 15(2): Der Kiureghian, Armen & Liu, Pei-Ling (1985) Structural Reliability Under Incomplete Probability Information, Report No. UCB/SESM/85-01, Department of Civil Engineering, University of California, Berkeley. Der Kiureghian, Armen & Liu, Pei-Ling (1986) Structural Reliability Under Incomplete Probability Information. Journal of Engineering Mechanics, ASCE 112(1): Der Kiureghian, Armen, Lin, H-Z. & Hwang, S.-J. (1987) Second-Order Reliability Approximations. Journal of Engineering Mechanics, ASCE 113(8): Freeze, R. A. (1975) A Stochastic-Conceptual Analysis of One-Dimensional Groundwater Flow in Non-Uniform Porous Media. Water Resour. Res. 11: Gelhar, L. W. (1976) Effects of Hydraulic Conductivity Variations in Groundwater Flow. Paper presented at 2nd IAHR International Symposium on Stochastic Hydraulics, Lund, Sweden. Gutjahr, A. L., Gelhard, L. W., Bakr, A. A. & MacMillan, J. R. (1978) Stochastic Analysis of Spatial Variability in Subsurface Flows, 2. Evaluation and Application. Water Resour. Res. 14(5): Gutjahr, A. L. & Gelhar, L. W. (1981) Stochastic Models of Subsurface Flow: Infinite Versus Finite Domains and Stationarity. Water Resour. Res. 17(2): Hasofer, A. M. & Lind, N. (1974) An Exact Invariant First-Order Reliability Format. Journal of Engineering Mechanics, ASCE 100. Hohenbilcher, M. & Rackwitz, R. (1981) Non-Normal Dependent Vectors in Structural Safety. Journal of Engineering Mechanics, ASCE 107(EM6): Liu, Pei-Ling & Der Kiureghian, Armen (1986) Multivariate Distribution Models With Prescribed Marginals and Covariances. Probabilistic Engineering Mechanics 1(2): McMillan, W. D. (1966) Theoretical Analysis of Groundwater Basin Operations. Water Resources Center Contribution 114, University of California. Madsen, H. O., Krenk, S. & Lind, N. C (1986) Methods of Sn-uctural Safety. Englewood Cliffs, NJ: Prentice-Hall Inc. Mizell, S. A., Gelhar, L. W. & Gutjahr, A. L. (1982) Stochastic Analysis of Spatial Variability in Two-Dimensional Steady Groundwater Flow Assuming Stationary and -215-
26 Non-Stationary Heads. Water Resour. Res. 18(4): Sitar, N., Cawlfield, J. D. & Der Kiureghian, Armen (1987) First-Order Reliability Approach to Stochastic Analysis of Subsurface Flow and Contaminant Transport. Water Resour. Res. 23(5): Smith, L. & Freeze, R. A. (1979a) Stochastic Analysis of Steady State Groundwater Flow in a Bounded Domain, 1. One-Dimensional Simulations. Water Resour. Res. 15(3): Smith, L., & Freeze, R. A. (1979b) Stochastic Analysis of Steady State Groundwater Flow in a Bounded Domain, 2. Two-Dimensional Simulations. Water Resour. Res. 15(6): Tang, W. H. (1984) Principles of Probabilistic Characterization of Soil Properties, in Probabilistic Characterization of Soil Properties: Bridge Between Theory and Practice. Proc. of ASCE Symposium (D.S. Bowles and H.-Y. Ko, eds.), Tang, D. H. & Pinder, G. F. (1977) Simulation of Groundwater Flow and Mass Transport. Advan. Water Res. 1(1): Vanmarke, E. H. (1979) Random Fields. Cambridge, MA: MIT Press. Warren, J. E. & Price, H. S. (1961) Flow in Heterogeneous Porous Media. Society of Petroleum Eng. Journal 1:
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 informationStructural Reliability
Structural Reliability Thuong Van DANG May 28, 2018 1 / 41 2 / 41 Introduction to Structural Reliability Concept of Limit State and Reliability Review of Probability Theory First Order Second Moment Method
More informationEFFICIENT MODELS FOR WIND TURBINE EXTREME LOADS USING INVERSE RELIABILITY
Published in Proceedings of the L00 (Response of Structures to Extreme Loading) Conference, Toronto, August 00. EFFICIENT MODELS FOR WIND TURBINE ETREME LOADS USING INVERSE RELIABILITY K. Saranyasoontorn
More informationA Simple Third-Moment Method for Structural Reliability
A Simple Third-Moment Method for Structural Reliability Yan-Gang Zhao* 1, Zhao-Hui Lu 2 and Tetsuro Ono 3 1 Associate Professor, Nagoya Institute of Technology, Japan 2 Graduate Student, Nagoya Institute
More informationTHIRD-MOMENT STANDARDIZATION FOR STRUCTURAL RELIABILITY ANALYSIS
THIRD-MOMENT STANDARDIZATION FOR STRUCTURAL RELIABILITY ANALYSIS By Yan-Gang Zhao and Tetsuro Ono ABSTRACT: First- and second-order reliability methods are generally considered to be among the most useful
More informationCreated by Erik Kostandyan, v4 January 15, 2017
MATLAB Functions for the First, Second and Inverse First Order Reliability Methods Copyrighted by Erik Kostandyan, Contact: erik.kostandyan.reliability@gmail.com Contents Description... References:...
More informationcomponent risk analysis
273: Urban Systems Modeling Lec. 3 component risk analysis instructor: Matteo Pozzi 273: Urban Systems Modeling Lec. 3 component reliability outline risk analysis for components uncertain demand and uncertain
More informationTransactions on Modelling and Simulation vol 18, 1997 WIT Press, ISSN X
Stochastic groundwater flow using boundary elements C. Z. Karakostas* and G. D. Manolis** ^Institute ofengineering Seismology and Earthquake Engineering, P. O. Box 53, Thessaloniki GR 55102, Greece **Department
More informationEFFICIENT SHAPE OPTIMIZATION USING POLYNOMIAL CHAOS EXPANSION AND LOCAL SENSITIVITIES
9 th ASCE Specialty Conference on Probabilistic Mechanics and Structural Reliability EFFICIENT SHAPE OPTIMIZATION USING POLYNOMIAL CHAOS EXPANSION AND LOCAL SENSITIVITIES Nam H. Kim and Haoyu Wang University
More informationSensitivity and Reliability Analysis of Nonlinear Frame Structures
Sensitivity and Reliability Analysis of Nonlinear Frame Structures Michael H. Scott Associate Professor School of Civil and Construction Engineering Applied Mathematics and Computation Seminar April 8,
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 informationBasics of Uncertainty Analysis
Basics of Uncertainty Analysis Chapter Six Basics of Uncertainty Analysis 6.1 Introduction As shown in Fig. 6.1, analysis models are used to predict the performances or behaviors of a product under design.
More informationFourth-Moment Standardization for Structural Reliability Assessment
Fourth-Moment Standardization for Structural Reliability Assessment Yan-Gang Zhao, M.ASCE 1 ; and Zhao-Hui Lu Abstract: In structural reliability analysis, the uncertainties related to resistance and load
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 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 informationPRACTICAL FIRST-ORDER RELIABILITY COMPUTATIONS USING SPREADSHEET
PRACTICAL FIRST-ORDER RELIABILITY COMPUTATIONS USING SPREADSHEET B. K. Low Associate Professor Geotechnical Research Centre School of Civil & Environ. Engineering Nanyang Technological University Republic
More informationAnalysis of Regression and Bayesian Predictive Uncertainty Measures
Analysis of and Predictive Uncertainty Measures Dan Lu, Mary C. Hill, Ming Ye Florida State University, dl7f@fsu.edu, mye@fsu.edu, Tallahassee, FL, USA U.S. Geological Survey, mchill@usgs.gov, Boulder,
More informationStatistícal Methods for Spatial Data Analysis
Texts in Statistícal Science Statistícal Methods for Spatial Data Analysis V- Oliver Schabenberger Carol A. Gotway PCT CHAPMAN & K Contents Preface xv 1 Introduction 1 1.1 The Need for Spatial Analysis
More informationMultivariate Distribution Models
Multivariate Distribution Models Model Description While the probability distribution for an individual random variable is called marginal, the probability distribution for multiple random variables is
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 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 informationNew Developments in Tail-Equivalent Linearization method for Nonlinear Stochastic Dynamics
New Developments in Tail-Equivalent Linearization method for Nonlinear Stochastic Dynamics Armen Der Kiureghian President, American University of Armenia Taisei Professor of Civil Engineering Emeritus
More informationAcceptable Ergodic Fluctuations and Simulation of Skewed Distributions
Acceptable Ergodic Fluctuations and Simulation of Skewed Distributions Oy Leuangthong, Jason McLennan and Clayton V. Deutsch Centre for Computational Geostatistics Department of Civil & Environmental Engineering
More information4. Distributions of Functions of Random Variables
4. Distributions of Functions of Random Variables Setup: Consider as given the joint distribution of X 1,..., X n (i.e. consider as given f X1,...,X n and F X1,...,X n ) Consider k functions g 1 : R n
More informationAnalysis of structural reliability under parameter uncertainties
Probabilistic Engineering Mechanics 23 (2008) 351 358 www.elsevier.com/locate/probengmech Analysis of structural reliability under parameter uncertainties Armen Der Kiureghian Department of Civil and Environmental
More informationReliability analysis of geotechnical risks
Reliability analysis of geotechnical risks Lazhar Belabed*, Hacene Benyaghla* * Department of Civil Engineering and Hydraulics, University of Guelma, Algeria Abstract The evaluation of safety or reliability
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 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 informationReduction of Random Variables in Structural Reliability Analysis
Reduction of Random Variables in Structural Reliability Analysis S. Adhikari and R. S. Langley Department of Engineering University of Cambridge Trumpington Street Cambridge CB2 1PZ (U.K.) February 21,
More informationScale of Fluctuation for Geotechnical Probabilistic Analysis
834 Geotechnical Safety and Risk V T. Schweckendiek et al. (Eds.) 015 The authors and IOS Press. This article is published online with Open Access by IOS Press and distributed under the terms of the Creative
More informationOn the Optimal Scaling of the Modified Metropolis-Hastings algorithm
On the Optimal Scaling of the Modified Metropolis-Hastings algorithm K. M. Zuev & J. L. Beck Division of Engineering and Applied Science California Institute of Technology, MC 4-44, Pasadena, CA 925, USA
More informationReliability-based calibration of design seismic response spectra and structural acceptance criteria
Reliability-based calibration of design seismic response spectra and structural acceptance criteria C. Loth & J. W. Baker Department of Civil and Environmental Engineering Stanford University ABSTRACT:
More informationEnvironmental Contours for Determination of Seismic Design Response Spectra
Environmental Contours for Determination of Seismic Design Response Spectra Christophe Loth Modeler, Risk Management Solutions, Newark, USA Jack W. Baker Associate Professor, Dept. of Civil and Env. Eng.,
More informationRisk Assessment of Highway Bridges: A Reliability-based Approach
Risk Assessment of Highway Bridges: A Reliability-based Approach by Reynaldo M. Jr., PhD Indiana University-Purdue University Fort Wayne pablor@ipfw.edu Abstract: Many countries are currently experiencing
More informationA general procedure for rst/second-order reliability method (FORM/SORM)
Structural Safety 21 (1999) 95±112 www.elsevier.nl/locate/strusafe A general procedure for rst/second-order reliability method (FORM/SORM) Yan-Gang Zhao*, Tetsuro Ono Department of Architecture, Nagoya
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 informationHowever, reliability analysis is not limited to calculation of the probability of failure.
Probabilistic Analysis probabilistic analysis methods, including the first and second-order reliability methods, Monte Carlo simulation, Importance sampling, Latin Hypercube sampling, and stochastic expansions
More informationMultivariate Distributions
IEOR E4602: Quantitative Risk Management Spring 2016 c 2016 by Martin Haugh Multivariate Distributions We will study multivariate distributions in these notes, focusing 1 in particular on multivariate
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 informationReliability of Acceptance Criteria in Nonlinear Response History Analysis of Tall Buildings
Reliability of Acceptance Criteria in Nonlinear Response History Analysis of Tall Buildings M.M. Talaat, PhD, PE Senior Staff - Simpson Gumpertz & Heger Inc Adjunct Assistant Professor - Cairo University
More informationA probabilistic method to predict fatigue crack initiation
International Journal of Fracture (2006) 137:9 17 DOI 10.1007/s10704-005-3074-0 Springer 2006 A probabilistic method to predict fatigue crack initiation SALIL. S. KULKARNI, L. SUN, B. MORAN, S. KRISHNASWAMY
More informationFigure The different sources of statistical uncertainties
3.4 Uncertainties in the treatment of statistical data and influences on the structural reliability assessment * 3.4. Introduction Structural reliability methods allow us to account for the uncertain part
More informationSUGGESTED ANALYTICAL MODEL FOR LIVE LOADS
SUGGESTED ANALYTICAL MODEL FOR LIVE LOADS A. B. Khalil, Ph. D. Assistance Professor Cairo University M. A. Ahmed, Ph.D. Assistance Professor Cairo University M. A. El-Reedy, Ph. D. Structure Engineer Email:elreedyma@yahoo.com
More informationThe Influence of Piers on the Risk of Flooding Upstream from a Bridge
Archives of Hydro-Engineering and Environmental Mechanics Vol. 52 (2005), No. 4, pp. 287 301 IBW PAN, ISSN 1231 3726 The Influence of Piers on the Risk of Flooding Upstream from a Bridge Marek Sowiński
More informationResearch Collection. Basics of structural reliability and links with structural design codes FBH Herbsttagung November 22nd, 2013.
Research Collection Presentation Basics of structural reliability and links with structural design codes FBH Herbsttagung November 22nd, 2013 Author(s): Sudret, Bruno Publication Date: 2013 Permanent Link:
More informationReliability Analysis Methods
Reliability Analysis Methods Emilio Bastidas-Arteaga, Abdel-Hamid Soubra To cite this version: Emilio Bastidas-Arteaga, Abdel-Hamid Soubra. Reliability Analysis Methods. Michael A. Hicks; Cristina Jommi.
More informationThe Generalized Likelihood Uncertainty Estimation methodology
CHAPTER 4 The Generalized Likelihood Uncertainty Estimation methodology Calibration and uncertainty estimation based upon a statistical framework is aimed at finding an optimal set of models, parameters
More informationIntroduction. Spatial Processes & Spatial Patterns
Introduction Spatial data: set of geo-referenced attribute measurements: each measurement is associated with a location (point) or an entity (area/region/object) in geographical (or other) space; the domain
More informationDEEP LEARNING CHAPTER 3 PROBABILITY & INFORMATION THEORY
DEEP LEARNING CHAPTER 3 PROBABILITY & INFORMATION THEORY OUTLINE 3.1 Why Probability? 3.2 Random Variables 3.3 Probability Distributions 3.4 Marginal Probability 3.5 Conditional Probability 3.6 The Chain
More informationGuideline for Offshore Structural Reliability Analysis - General 3. RELIABILITY ANALYSIS 38
FEBRUARY 20, 1995 3. RELIABILITY ANALYSIS 38 3.1 General 38 3.1.1 Variables 38 3.1.2 Events 39 3.1.3 Event Probability 41 3.1.4 The Reliability Index 41 3.1.5 The Design Point 42 3.1.6 Transformation of
More informationMULTIVARIATE STOCHASTIC SIMULATION WITH SUBJECTIVE MULTIVARIATE NORMAL DISTRIBUTIONS1. Peter J. Ince2 and Joseph Buongiorno3
MULTIVARIATE STOCHASTIC SIMULATION WITH SUBJECTIVE MULTIVARIATE NORMAL DISTRIBUTIONS1 Peter J. Ince2 and Joseph Buongiorno3 Abstract.-In many applications of Monte Carlo simulation in forestry or forest
More information7 Geostatistics. Figure 7.1 Focus of geostatistics
7 Geostatistics 7.1 Introduction Geostatistics is the part of statistics that is concerned with geo-referenced data, i.e. data that are linked to spatial coordinates. To describe the spatial variation
More informationUNIFORMLY MOST POWERFUL CYCLIC PERMUTATION INVARIANT DETECTION FOR DISCRETE-TIME SIGNALS
UNIFORMLY MOST POWERFUL CYCLIC PERMUTATION INVARIANT DETECTION FOR DISCRETE-TIME SIGNALS F. C. Nicolls and G. de Jager Department of Electrical Engineering, University of Cape Town Rondebosch 77, South
More informationA Polynomial Chaos Approach to Robust Multiobjective Optimization
A Polynomial Chaos Approach to Robust Multiobjective Optimization Silvia Poles 1, Alberto Lovison 2 1 EnginSoft S.p.A., Optimization Consulting Via Giambellino, 7 35129 Padova, Italy s.poles@enginsoft.it
More informationTIME-VARIANT RELIABILITY ANALYSIS FOR SERIES SYSTEMS WITH LOG-NORMAL VECTOR RESPONSE
TIME-VARIANT RELIABILITY ANALYSIS FOR SERIES SYSTEMS WITH LOG-NORMAL VECTOR RESONSE Sayan Gupta, ieter van Gelder, Mahesh andey 2 Department of Civil Engineering, Technical University of Delft, The Netherlands
More information= _(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
WATER RESOURCES RESEARCH, VOL. 35, NO. 7, PAGES 2273-2277, JULY 999 A general method for obtaining analytical expressions for the first-order velocity covariance in heterogeneous porous media Kuo-Chin
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 informationStochastic Processes. M. Sami Fadali Professor of Electrical Engineering University of Nevada, Reno
Stochastic Processes M. Sami Fadali Professor of Electrical Engineering University of Nevada, Reno 1 Outline Stochastic (random) processes. Autocorrelation. Crosscorrelation. Spectral density function.
More information22 : Hilbert Space Embeddings of Distributions
10-708: Probabilistic Graphical Models 10-708, Spring 2014 22 : Hilbert Space Embeddings of Distributions Lecturer: Eric P. Xing Scribes: Sujay Kumar Jauhar and Zhiguang Huo 1 Introduction and Motivation
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 informationSTA 4273H: Statistical Machine Learning
STA 4273H: Statistical Machine Learning Russ Salakhutdinov Department of Statistics! rsalakhu@utstat.toronto.edu! http://www.utstat.utoronto.ca/~rsalakhu/ Sidney Smith Hall, Room 6002 Lecture 7 Approximate
More informationCALIBRATED RESPONSE SPECTRA FOR COLLAPSE ASSESSMENT UNDER MULTIVARIATE HAZARD AND STRUCTURAL RESPONSE UNCERTAINTIES
10NCEE Tenth U.S. National Conference on Earthquake Engineering Frontiers of Earthquake Engineering July 21-25, 2014 Anchorage, Alaska CALIBRATED RESPONSE SPECTRA FOR COLLAPSE ASSESSMENT UNDER MULTIVARIATE
More informationSafety Envelope for Load Tolerance and Its Application to Fatigue Reliability Design
Safety Envelope for Load Tolerance and Its Application to Fatigue Reliability Design Haoyu Wang * and Nam H. Kim University of Florida, Gainesville, FL 32611 Yoon-Jun Kim Caterpillar Inc., Peoria, IL 61656
More informationENSC327 Communications Systems 19: Random Processes. Jie Liang School of Engineering Science Simon Fraser University
ENSC327 Communications Systems 19: Random Processes Jie Liang School of Engineering Science Simon Fraser University 1 Outline Random processes Stationary random processes Autocorrelation of random processes
More informationA Note on Bootstraps and Robustness. Tony Lancaster, Brown University, December 2003.
A Note on Bootstraps and Robustness Tony Lancaster, Brown University, December 2003. In this note we consider several versions of the bootstrap and argue that it is helpful in explaining and thinking about
More informationPreliminary Statistics Lecture 2: Probability Theory (Outline) prelimsoas.webs.com
1 School of Oriental and African Studies September 2015 Department of Economics Preliminary Statistics Lecture 2: Probability Theory (Outline) prelimsoas.webs.com Gujarati D. Basic Econometrics, Appendix
More informationStatistical Rock Physics
Statistical - Introduction Book review 3.1-3.3 Min Sun March. 13, 2009 Outline. What is Statistical. Why we need Statistical. How Statistical works Statistical Rock physics Information theory Statistics
More informationRESPONSE SURFACE METHODS FOR STOCHASTIC STRUCTURAL OPTIMIZATION
Meccanica dei Materiali e delle Strutture Vol. VI (2016), no.1, pp. 99-106 ISSN: 2035-679X Dipartimento di Ingegneria Civile, Ambientale, Aerospaziale, Dei Materiali DICAM RESPONSE SURFACE METHODS FOR
More informationNonparametric Bayesian Methods (Gaussian Processes)
[70240413 Statistical Machine Learning, Spring, 2015] Nonparametric Bayesian Methods (Gaussian Processes) Jun Zhu dcszj@mail.tsinghua.edu.cn http://bigml.cs.tsinghua.edu.cn/~jun State Key Lab of Intelligent
More informationA conceptual interpretation of the renewal theorem with applications
Risk, Reliability and Societal Safety Aven & Vinnem (eds) 2007 Taylor & Francis Group, London, ISBN 978-0-415-44786-7 A conceptual interpretation of the renewal theorem with applications J.A.M. van der
More informationarxiv:hep-ex/ v1 17 Sep 1999
Propagation of Errors for Matrix Inversion M. Lefebvre 1, R.K. Keeler 1, R. Sobie 1,2, J. White 1,3 arxiv:hep-ex/9909031v1 17 Sep 1999 Abstract A formula is given for the propagation of errors during matrix
More informationLessons in Estimation Theory for Signal Processing, Communications, and Control
Lessons in Estimation Theory for Signal Processing, Communications, and Control Jerry M. Mendel Department of Electrical Engineering University of Southern California Los Angeles, California PRENTICE HALL
More informationA Moving Kriging Interpolation Response Surface Method for Structural Reliability Analysis
Copyright 2013 Tech Science Press CMES, vol.93, no.6, pp.469-488, 2013 A Moving Kriging Interpolation Response Surface Method for Structural Reliability Analysis W. Zhao 1,2, J.K. Liu 3, X.Y. Li 2, Q.W.
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 informationL03. PROBABILITY REVIEW II COVARIANCE PROJECTION. NA568 Mobile Robotics: Methods & Algorithms
L03. PROBABILITY REVIEW II COVARIANCE PROJECTION NA568 Mobile Robotics: Methods & Algorithms Today s Agenda State Representation and Uncertainty Multivariate Gaussian Covariance Projection Probabilistic
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 informationUsing Method of Moments in Schedule Risk Analysis
Raymond P. Covert MCR, LLC rcovert@mcri.com BACKGROUND A program schedule is a critical tool of program management. Program schedules based on discrete estimates of time lack the necessary information
More informationShannon meets Wiener II: On MMSE estimation in successive decoding schemes
Shannon meets Wiener II: On MMSE estimation in successive decoding schemes G. David Forney, Jr. MIT Cambridge, MA 0239 USA forneyd@comcast.net Abstract We continue to discuss why MMSE estimation arises
More informationApplications of Randomized Methods for Decomposing and Simulating from Large Covariance Matrices
Applications of Randomized Methods for Decomposing and Simulating from Large Covariance Matrices Vahid Dehdari and Clayton V. Deutsch Geostatistical modeling involves many variables and many locations.
More informationSTA 4273H: Statistical Machine Learning
STA 4273H: Statistical Machine Learning Russ Salakhutdinov Department of Computer Science! Department of Statistical Sciences! rsalakhu@cs.toronto.edu! h0p://www.cs.utoronto.ca/~rsalakhu/ Lecture 7 Approximate
More informationSTA414/2104. Lecture 11: Gaussian Processes. Department of Statistics
STA414/2104 Lecture 11: Gaussian Processes Department of Statistics www.utstat.utoronto.ca Delivered by Mark Ebden with thanks to Russ Salakhutdinov Outline Gaussian Processes Exam review Course evaluations
More informationGroundwater Resources Management under. Uncertainty Harald Kunstmann*, Wolfgang Kinzelbach* & Gerrit van Tender**
Groundwater Resources Management under Uncertainty Harald Kunstmann*, Wolfgang Kinzelbach* & Gerrit van Tender**, C/f &OP3 Zwnc/?, Email: kunstmann@ihw. baum. ethz. ch Abstract Groundwater resources and
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 informationUncertainty Quantification and Validation Using RAVEN. A. Alfonsi, C. Rabiti. Risk-Informed Safety Margin Characterization. https://lwrs.inl.
Risk-Informed Safety Margin Characterization Uncertainty Quantification and Validation Using RAVEN https://lwrs.inl.gov A. Alfonsi, C. Rabiti North Carolina State University, Raleigh 06/28/2017 Assumptions
More informationModule 8. Lecture 5: Reliability analysis
Lecture 5: Reliability analysis Reliability It is defined as the probability of non-failure, p s, at which the resistance of the system exceeds the load; where P() denotes the probability. The failure
More informationCMCE Computational Methods in Civil Engineering
Comp. Meth. Civil Eng. 2010, Vol. 1 No. 1 pp. 55-64 Copyright by the University of Guilan, Printed in I.R. Iran CMCE Computational Methods in Civil Engineering New method for estimation of the scale of
More informationReview (probability, linear algebra) CE-717 : Machine Learning Sharif University of Technology
Review (probability, linear algebra) CE-717 : Machine Learning Sharif University of Technology M. Soleymani Fall 2012 Some slides have been adopted from Prof. H.R. Rabiee s and also Prof. R. Gutierrez-Osuna
More informationDiscriminative Direction for Kernel Classifiers
Discriminative Direction for Kernel Classifiers Polina Golland Artificial Intelligence Lab Massachusetts Institute of Technology Cambridge, MA 02139 polina@ai.mit.edu Abstract In many scientific and engineering
More informationInformation Updating in Infrastructure Systems
University of Washington, Seattle Civil and Environmental Engineering April 10, 2008 Information Updating in Infrastructure Systems Subject to Multiple Hazards Daniel Straub University of California, Berkeley
More informationRelevance Vector Machines for Earthquake Response Spectra
2012 2011 American American Transactions Transactions on on Engineering Engineering & Applied Applied Sciences Sciences. American Transactions on Engineering & Applied Sciences http://tuengr.com/ateas
More informationStochastic methods for aquifer protection and management
Stochastic methods for aquifer protection and management Dr Adrian Butler Department of Civil & Environmental Engineering G-WADI 07 International Workshop on Groundwater Modeling for Arid and Semi-arid
More informationUSING GEOSTATISTICS TO DESCRIBE COMPLEX A PRIORI INFORMATION FOR INVERSE PROBLEMS THOMAS M. HANSEN 1,2, KLAUS MOSEGAARD 2 and KNUD S.
USING GEOSTATISTICS TO DESCRIBE COMPLEX A PRIORI INFORMATION FOR INVERSE PROBLEMS THOMAS M. HANSEN 1,2, KLAUS MOSEGAARD 2 and KNUD S. CORDUA 1 1 Institute of Geography & Geology, University of Copenhagen,
More informationReview of Probability Theory
Review of Probability Theory Arian Maleki and Tom Do Stanford University Probability theory is the study of uncertainty Through this class, we will be relying on concepts from probability theory for deriving
More informationMAXIMUM ENTROPY-BASED UNCERTAINTY MODELING AT THE FINITE ELEMENT LEVEL. Pengchao Song and Marc P. Mignolet
MAXIMUM ENTROPY-BASED UNCERTAINTY MODELING AT THE FINITE ELEMENT LEVEL Pengchao Song and Marc P. Mignolet SEMTE, Faculties of Mechanical and Aerospace Engineering, Arizona State University, 51 E. Tyler
More informationStatistics for scientists and engineers
Statistics for scientists and engineers February 0, 006 Contents Introduction. Motivation - why study statistics?................................... Examples..................................................3
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 informationApplication of Copulas as a New Geostatistical Tool
Application of Copulas as a New Geostatistical Tool Presented by Jing Li Supervisors András Bardossy, Sjoerd Van der zee, Insa Neuweiler Universität Stuttgart Institut für Wasserbau Lehrstuhl für Hydrologie
More informationA6523 Modeling, Inference, and Mining Jim Cordes, Cornell University
A6523 Modeling, Inference, and Mining Jim Cordes, Cornell University Lecture 19 Modeling Topics plan: Modeling (linear/non- linear least squares) Bayesian inference Bayesian approaches to spectral esbmabon;
More informationSimulation of random fields in structural design. Lectures. Sebastian Wolff
Lectures Simulation of random fields in structural design Sebastian Wolff presented at the 11th International Probabilistic Workshop, Brno 2013 Source: www.dynardo.de/en/library Novák and Vořechovský:
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 information