Diagrammatic solutions for hydraulic head moments in 1-D and 2-D bounded domains

We present a diagrammatic method for solving stochastic 1-D and 2-D steady-state flow equations in bounded domains. The diagrammatic method results in explicit solutions for the moments of the hydraulic head. This avoids certain numerical constraints encountered in realization-based methods. The diagrammatic technique also allows for the consideration of finite domains or large fluctuations, and is not restricted by distributional assumptions. The results of the method for 1-D and 2-D finite domains are compared with those obtained through a realization-based approach. Mean and variance of head are well reproduced for all log-conductivity variances inputted, including those larger than one. The diagrammatic results also compare favorably to hydraulic head moments derived by standard analytic methods requiring a linearized form of the flow equation.     

1-D, 2-D, and 3-D analytical solutions of unsaturated flow in groundwater

An excellent tool for checking numerical models of unsaturated flow in groundwater is analytical solutions. However, because of the highly nonlinear nature of the governing partial differential equation, only a limited number of analytical solutions are available. This paper first gives some simple 1-D solutions. Next, by use of a transformation, the nonlinear partial differential equation is converted to a linear one for a specific form of the moisture content vs. pressure head and relative hydraulic conductivity vs. pressure head curves. This allows both 2-D and 3-D solutions to be derived, which is done in this paper. Finally, computations from a finite element computer program are compared with results from one of the analytical solutions to illustrate the use of the derived equations.     

Stochastic perturbation analysis of groundwater flow. Spatially variable soils,semi-infinite domains and large fluctuations

As is well known, a complete stochastic solution of the stochastic differential equation governing saturated groundwater flow leads to an infinite hierarchy of equations in terms of higher-order moments. Perturbation techniques are commonly used to close this hierarchy, using power-series expansions. These methods are applied by truncating the series after a finite number of terms, and products of random gradients of conductivity and head potential are neglected. Uncertainty regarding the number or terms required to yield a sufficiently accurate result is a significant drawback with the application of power series-based perturbation methods for such problems. Low-order series truncation may be incapable of representing fundamental characteristics of flow and can lead to physically unreasonable and inaccurate solutions of the stochastic flow equation. To support this argument, one-dimensional, steady-state, saturated groundwater flow is examined, for the case of a spatially distributed hydraulic conductivity field. An ordinary power-series perturbation method is used to approximate the mean head, using second-order statistics to characterize the conductivity field. Then an interactive perturbation approach is introduced, which yields improved results compared to low-order, power-series perturbation methods for situations where strong interactions exist between terms in such approximations. The interactive perturbation concept is further developed using Feynman-type diagrams and graph theory, which reduce the original stochastic flow problem to a closed set of equations for the mean and the covariance functions. Both theoretical and practical advantages of diagrammatic solutions are discussed; these include the study of bounded domains and large fluctuations.     

A HORIZONTAL 2-D HYDRAULIC NUMERICAL MODEL AND IT’S APPLICATIONS TO FLOOD FORECAST

1 INTRODUCTION There were many flood disasters in China in recent years. When the water level in a river is very high, weak parts of its dike may be destroyed resulted in the submersion of the protected land and properties. It is of significance for decision-makers to exactly predict the processes of flood propagation during flood control. There are many modes of dike bursting, such as seepage destroying by overflow on top of dike caused by dike body sinking induced by piping and soil fl…     

Numerical approximation of head and flux covariances in three dimensions using mixed finite elements

A numerical method is developed for accurately approximating head and flux covariances and cross-covariances in finite two- and three-dimensional domains using the mixed finite element method. The method is useful for determining head and flux covariances for non-stationary flow fields, for example those induced by injection or extraction wells, impermeable subsurface barriers, or non-stationary hydraulic conductivity fields. Because the numerical approximations to the flux covariances are obtained directly from the solution to the coupled problem rather than having to differentiate head covariances, the approximations are in general more accurate than those obtained from conventional finite difference or finite element methods. Results for uniform flow example problems are consistent with results from previously published finite domain analyses and demonstrate that head variances and covariances are quite sensitive to boundary conditions and the size of the bounded domain. Flux variances and covariances are less sensitive to boundary conditions and domain size. Results comparing approximations from lower-order Raviart–Thomas–Nedelec and higher order Brezzi–Douglas–Marini[9] finite element spaces indicate that higher order element space improve the estimate of the flux covariances, but do not significantly affect the estimate of the head covariances.     

Explicit algebraic influence coefficients: two-dimensional aquifer model / Coefficients d'influence algébriques explicites: modèle d'aquifère bidimensionnel

Abstract

This work extends the algebraic expression of influence coefficients developed for one-dimensional aquifer models to a two-dimensional (2-D) case. First, the partial differential equation governing the flow in a 2-D semi-confined aquifer is discretized using a finite difference scheme. This results in a system of discrete equations presented in the form of water balance equations associated with a network of interconnected compartments centred on the grid nodes. The foregoing system is transformed into a series of uncoupled 1-D equations stated in terms of some generalized hydraulic head for which they are also solved. Second, the original hydraulic head is recovered from the generalized one via an appropriate linear transformation. Whence, the algebraic expression making the hydraulic head explicit versus sources and boundary conditions is derived. This discrete expression, mapped onto its continuous counterpart, helps to deduce an algebraic form of the inter-compartment influence coefficients. Finally, a comparison with the analytical Green function is carried out.     

Numerical method of moments for solute transport in a nonstationary flow field conditioned on hydraulic conductivity and head measurements

The unconditional stochastic studies on groundwater flow and solute transport in a nonstationary conductivity field show that the standard deviations of the hydraulic head and solute flux are very large in comparison with their mean values (Zhang et al. in Water Resour Res 36:2107–2120, 2000; Wu et al. in J Hydrol 275:208–228, 2003; Hu et al. in Adv Water Resour 26:513–531, 2003). In this study, we develop a numerical method of moments conditioning on measurements of hydraulic conductivity and head to reduce the variances of the head and the solute flux. A Lagrangian perturbation method is applied to develop the framework for solute transport in a nonstationary flow field. Since analytically derived moments equations are too complicated to solve analytically, a numerical finite difference method is implemented to obtain the solutions. Instead of using an unconditional conductivity field as an input to calculate groundwater velocity, we combine a geostatistical method and a method of moment for flow to conditionally simulate the distributions of head and velocity based on the measurements of hydraulic conductivity and head at some points. The developed theory is applied in several case studies to investigate the influences of the measurements of hydraulic conductivity and/or the hydraulic head on the variances of the predictive head and the solute flux in nonstationary flow fields. The study results show that the conditional calculation will significantly reduce the head variance. Since the hydraulic head measurement points are treated as the interior boundary (Dirichlet boundary) conditions, conditioning on both the hydraulic conductivity and the head measurements is much better than conditioning only on conductivity measurements for reduction of head variance. However, for solute flux, variance reduction by the conditional study is not so significant.     

Sensitivity and moment analyses of head in variably saturated regimes

A numerical approach for approximating statistical moments of hydraulic heads of variably saturated flows in multi-dimensional porous media is developed. The approximation relies on a first-order Taylor series expansion of a finite element flow model and an adjoint state numerical method for variably saturated flows to evaluate sensitivities. This approach can be employed to analyze uncertainties associated with predictions of head of steady-state or transient flows in variably saturated porous media, with any type of boundary and initial conditions. Limitations of stochastic analytical methods such as spectral/perturbation approaches and the time-consuming Monte Carlo simulation technique are thus alleviated. An example is given to demonstrate the utility of the approach and to investigate the temporal evolution of head variances in a variably saturated flow regime. Results show that the fluctuation of the water table can have significant impacts on the propagation of the head variance.     

Mass-conservative finite volume methods on 2-D unstructured grids for the Richards’ equation

The solution to the 2-D time-dependent unsaturated flow equation is numerically approximated by a second-order accurate cell-centered finite-volume discretization on unstructured grids. The approximation method is based on a vertex-centered Least Squares linear reconstruction of the solution gradients at mesh edges.A Taylor series development in time of the water content dependent variable in a finite-difference framework guarantees that the proposed finite volume method is mass conservative. A Picard iterative scheme solves at each time step the resulting non-linear algebraic problem. The performance of the method is assessed on five different test cases and implementing four distinct soil constitutive relationships. The first test case deals with a column infiltration problem. It shows the capability of providing a mass-conservative behavior. The second test case verifies the numerical approximation by comparison with an analytical mixed saturated–unsaturated solution. In this case, the water drains from a fully saturated portion of a 1-D column. The third and fourth test cases illustrate the performance of the approximation scheme on sharp soil heterogeneities on 1-D and 2-D multi-layered infiltration problems. The 2-D case shows the passage of an abrupt infiltration front across a curved interface between two layers. Finally, the fifth test case compares the numerical results with an analytical solution that is developed for a 2-D heterogeneous soil with a source term representing plant roots. This last test case illustrates the formal second-order accuracy of the method in the numerical approximation of the pressure head.     

2-D SIMULATION OF CHANNEL FLOWS WITH MOVEABLE BED

1 INTRODUCTIONFor many hydraulic engineering problems, the analysis of flow and bed level variations in openchannels is a fundamental prerequisite. forcal methOds fOr alluvial rivers are well develoPednowadays as far as onediInensional descriPtions are concemed. A cOmPrhensive analysis of Ihe wellknown models is Presented by Habersack(l998). HOwever, for a number of Problems such as channelwidening, flow pattem close to sPuds and etc. a more deailed knowledge of the bed level behavio…     

A generalized transformation approach for simulating steady-state variably-saturated subsurface flow

A numerical method is proposed to accurately and efficiently compute a direct steady-state solution of the nonlinear Richards equation. In the proposed method, the Kirchhoff integral transformation and a complementary transformation are applied to the governing equation in order to separate the nonlinear hyperbolic characteristic from the linear parabolic part. The separation allows the transformed governing equation to be applied to partially- to fully-saturated systems with arbitrary constitutive relations between primary (pressure head) and secondary variables (relative permeability). The transformed governing equation is then discretized with control volume finite difference/finite element approximations, followed by inverse transformation. The approach is compared to analytical and other numerical approaches for variably-saturated flow in 1-D and 3-D domains. The results clearly demonstrate that the approach is not only more computationally efficient but also more accurate than traditional numerical solutions. The approach is also applied to an example flow problem involving a regional-scale variably-saturated heterogeneous system, where the vadose zone is up to 1 km thick. The performance, stability, and effectiveness of the transform approach is exemplified for this complex heterogeneous example, which is typical of many problems encountered in the field. It is shown that computational performance can be enhanced by several orders of magnitude with the described integral transformation approach.     

Finite element-fractional steps solution of the 3-D coastal circulation model

The mathematical model for the nearly horizontal circulation due to wind, tides and density gradients in 3-D coastal areas is solved by a combined use of the method of finite elements and the integration in fractional steps. The discretisation of the flow domain is achieved through a system of 1-D finite elements over the depth, z, and 2-D finite elements in x?y space. The differential operators of the momentum equations in x and y, are split and integrated separately in z and x?y dimensions. The method is an extension of a previously presented approach combining finite differences and expansion in series. The application refers to the wind induced circulation in the 3-D coastal basin of Thessaloniki Bay.     

An example of solute spreading in nonstationary, bounded geological formations

We analyze the impact of a linear trend in the mean log-conductivity on the transport of a conservative tracer in a bounded domain. The effects of such a linear trend on solute transport were analyzed in depth for unbounded domains (Rubin and Seong, Water Resour Res 30(11):2901–2911, 1994; Indelman and Rubin, Water Resour Res 31(5):1257–1265, 1995; Water Resour Res 32(5):1257–1265, 1996), whereas studies concerning this special case of medium nonstationarity in finite domains usually focus on head or flow statistics (Guadagnini et al., Stoch Environ Res Risk Assess, 17:394–407, 2003). In this study both ensemble and effective plume moments are provided for an instantaneous release of a solute through a linear source normal to the mean flow direction, by taking into account different sizes of the source. The analysis involving a steady velocity field spatially nonstationary is developed by using the stochastic finite element method. Results show that ensemble moments are affected by increasing trends both parallel and normal to the mean flow direction, but the impact on effective plume moments is very different. A parallel trend does not seem to influence the effective second moments; while a normal trend, although modifies the transverse effective moment only weakly, strongly increases the longitudinal one, especially for large initial sizes of the source. Furthermore, the increase of the particle displacement variance produced by a parallel trend in the finite domain disagrees with the results obtained in an unbounded domain, due to the boundary conditions here considered making both head and velocity moments nonstationary and nonsymmetric.     

Eulerian spatial moments for solute transport in three-dimensional heterogeneous,dual-permeability media

A Eulerian analytical method is developed for nonreactive solute transport in heterogeneous, dual-permeability media where the hydraulic conductivities in fracture and matrix domains are both assumed to be stochastic processes. The analytical solution for the mean concentration is given explicitly in Fourier and Laplace transforms. Instead of using the fast fourier transform method to numerically invert the solution to real space (Hu et al., 2002), we apply the general relationship between spatial moments and concentration (Naff, 1990; Hu et al., 1997) to obtain the analytical solutions for the spatial moments up to the second for a pulse input of the solute. Owing to its accuracy and efficiency, the analytical method can be used to check the semi-analytical and Monte Carlo numerical methods before they are applied to more complicated studies. The analytical method can be also used during screening studies to identify the most significant transport parameters for further analysis. In this study, the analytical results have been compared with those obtained from the semi-analytical method (Hu et al., 2002) and the comparison shows that the semi-analytical method is robust. It is clearly shown from the analytical solution that the three factors, local dispersion, conductivity variation in each domain and velocity convection flow difference in the two domains, play different roles on the solute plume spreading in longitudinal and transverse directions. The calculation results also indicate that when the log-conductivity variance in matrix is 10 times less than its counterpart in fractures, it will hardly influence the solute transport, whether the conductivity field is matrix is treated as a homogeneous or random field.     

Experimental and Numerical Investigation of the Flow Hydraulic in Gradual Transition Open Channels

Channel expansions are common in both natural and artificial open channels. With increasing cross-sectional dimensions in an expansion, the flow decelerates. Due to separation of flow and subsequent eddy formation, a significant head loss is occurred along the transition. This study presents the results of experimental investigations on subcritical flow along the expansive transition of rectangular to trapezoidal channels. Also, a numerical simulation was developed using the finite volume method with Reynolds Stress turbulent model. Water surface profiles and velocity distributions of flow through the transition were measured experimentally and compared with the numerical results. Also, hydraulic efficiency of the transition and coefficient of energy head loss were calculated. The results show that with increasing the upstream Froude number, hydraulic efficiency of the transition and coefficient of energy head loss are decreased and increased, respectively. The results also showed the ability of numerical simulation for simulating the flow separation zones and bed shear stress along the transition for different inlet discharges and inflow Froude numbers.     

构造热演化的数值模拟方法

对模拟构造热演化的二维和三维非稳态热传导-平流微分方程的数值求解进行了研究.在空间和时间上分别应用有限单元法和差分法进行离散化处理;推导了非稳态热传导-平流有限元方程.在求解技术上发展了“迎风”权函数方法,将有限元法、有限差分法和“迎风”权函数法结合起来,使用变网格、变结点方法对地质体的运动介质、运动边界、多热源的非线性、非稳态的构造热演化进行了模拟.并对各种热源体的处理、初始条件和边界条件的处理进行了分析.     

用边界单元法模拟二维地形对大地电磁场的影响

本文用边界单元法模拟二维地形对大地电磁场的影响.首先用格林公式将二维大地电磁的边值问题转变成积分方程,然后用边界单元法解积分方程,得到地形上的大地电磁场和它的法向导数,由此可计算电阻率.与有限单元法相比,本方法剖分后的地形与实际地形的拟合程度高,向计算机输入原始数据的工作简单,可以在微机上计算地形引起的视电阻率.本方法适合于在野外生产现场进行大地电磁法的地形改正.计算表明,地形对H_x型波的影响比对E_x型波严重得多;随着周期的增长,地形对H_x型波的影响变得稳定.     

LE TRITIUM DANS LES PRÉCIPITATIONS DE L'ASIE DU SUD-EST ET LE PHÉNOMÉNE DE LA MOUSSON EN 1966 ET EN 1967

Abstract

The problem of non-steady flow of water in a soil-plant system can be described by adding a sink term to the continuity equation for soil water flow. In this paper the sink term is defined in two different ways. Firstly it is considered to be dependent on the hydraulic conductivity of the soil, on the difference in pressure head between the soil and the root-soil interface and some root effectiveness function. Secondly the sink is taken to be a prescribed function of the soil water content. The partial differential equation applying to the first problem is solved by both a finite difference (FD 1) and a finite element (FE 1) technique, that applying to the second problem by a finite difference approach (FD 2). The purpose of this paper is to verify the numerical models against field measurements, to compare the results obtained by the three numerical methods and to show how the finite element method can be applied to complex but realistic two-dimensional flow situations. Two examples are given. The first concerns one-dimensional flow and it compares numerical results with those obtained experimentally in the field from water balance studies on red cabbage (Brassica oleracea L. ‘Rode Herfst’) grown on a clay soil in the presence of a water table. The second example describes two-dimensional flow in a complex field situation in the Netherlands where flow takes place under cropped field conditions through five anisotropic layers. Water is supplied to the system by infiltration from two unlined ditches and is withdrawn from the system by evapotranspiration and by leakage to an underlying pumped aquifer.     

Analysis of biodegradation in a 3-D heterogeneous porous medium using nonlinear stochastic finite element method

Probabilistic analysis by Monte Carlo Simulation method (MCSM) is a computationally prohibitive task for a reactive solute transport involving coupled PDEs with nonlinear source/sink terms in 3-D heterogeneous porous media. The perturbation based stochastic finite element method (SFEM) is an attractive alternative method to MCSM as it is computationally efficient and accurate. In the present study SFEM is developed for solving nonlinear reactive solute transport problem in a 3-D heterogeneous medium. Here the solution of the biodegradation problem involving a single solute by a single class of microorganisms coupled with dynamic microbial growth is attempted using this method. The SFEM here produces a second-order accurate solution for the mean and a first-order accurate solution for the standard deviation of concentrations. In this study both the physical parameters (hydraulic conductivity, porosity, dispersivity and diffusion coefficient) and the biological parameters (maximum substrate utilization rate and the coefficient of cell decay) are considered as spatially varying random fields. A comparison between the MCSM and SFEM for the mean and standard deviation of concentration is made for 1-D and 3-D problem. The effects of heterogeneity on the degradation of substrate and growth of biomass concentrations for a range of variances of input parameters are discussed for both 1-D and 3-D problems.     

From the pore scale to the lab scale: 3-D lab experiment and numerical simulation of drainage in heterogeneous porous media

A well-controlled 3-D experiment with pre-defined block heterogeneities is conducted, where neutron tomography is used to map 3-D water distribution after two successive drainage steps. The material and hydraulic properties of the two sands are first measured in the laboratory with multistep outflow experiments. Additionally, the pore structure of the sands is acquired by means of image analysis of synchrotron tomography data and the structure is used for pore-scale simulation of one- and two-phase flow with Lattice-Boltzmann methods. This gives us another set of material and hydraulic parameters of the sands. The two sets of hydraulic properties (from the lab scale and from the pore scale) are then used in numerical simulations of the 3-D experiment.