Numerical modelling of geoelectric structures using the galerkin and finite elements methods

Summary The Galerkin and finite elements methods are used to model numerically H- and E-polarized electromagnetic fields in inhomogeneous media. The error of approximation is studied and an example of numerical results is given.     

Construction of fast dynamos using unsteady flows and maps in three dimensions

Abstract

Unsteady fast dynamos are constructed using a sequence of rapid movements separated by periods of diffusion of the magnetic field. Motivated by the physical mechanism described by Ze'dovich and by Soward, we examine the effect of single-mode Beltrami waves applied sequentially, and show that they can be approximated by a simple “stretch-fold-shear” (SFS) map of the unit cube onto itself. In the SFS map, the field points in a fixed direction and diffusion is easily computed. The numerical results indicate that fast dynamo action occurs for sufficiently large shear, and that the process is primarily a coherent feature of the larger magnetic scales. Similar results are obtained for smooth flows. Some preliminary analysis of the SFS map using a decomposition method is described.     

Higher and lowest order mixed finite element approximation of subsurface flow problems with solutions of low regularity

In this work we study mixed finite element approximations of Richards’ equation for simulating variably saturated subsurface flow and simultaneous reactive solute transport. Whereas higher order schemes have proved their ability to approximate reliably reactive solute transport (cf., e.g. [Bause M, Knabner P. Numerical simulation of contaminant biodegradation by higher order methods and adaptive time stepping. Comput Visual Sci 7;2004:61–78]), the Raviart–Thomas mixed finite element method (RT₀) with a first order accurate flux approximation is popular for computing the underlying water flow field (cf. [Bause M, Knabner P. Computation of variably saturated subsurface flow by adaptive mixed hybrid finite element methods. Adv Water Resour 27;2004:565–581, Farthing MW, Kees CE, Miller CT. Mixed finite element methods and higher order temporal approximations for variably saturated groundwater flow. Adv Water Resour 26;2003:373–394, Starke G. Least-squares mixed finite element solution of variably saturated subsurface flow problems. SIAM J Sci Comput 21;2000:1869–1885, Younes A, Mosé R, Ackerer P, Chavent G. A new formulation of the mixed finite element method for solving elliptic and parabolic PDE with triangular elements. J Comp Phys 149;1999:148–167, Woodward CS, Dawson CN. Analysis of expanded mixed finite element methods for a nonlinear parabolic equation modeling flow into variably saturated porous media. SIAM J Numer Anal 37;2000:701–724]). This combination might be non-optimal. Higher order techniques could increase the accuracy of the flow field calculation and thereby improve the prediction of the solute transport. Here, we analyse the application of the Brezzi-Douglas-Marini element (BDM₁) with a second order accurate flux approximation to elliptic, parabolic and degenerate problems whose solutions lack the regularity that is assumed in optimal order error analyses. For the flow field calculation a superiority of the BDM₁ approach to the RT₀ one is observed, which however is less significant for the accompanying solute transport.     

Numerical modeling of the penetration area of reversed submerged jet by the methods of finite elements and finite volumes

Water jet flow has many usages in the field of management and water resource operation that can be applied in mixing, dilution and aerification. The current study has calculated the flow velocity, length and height of jet penetration area (jet and main flow are of opposite directions) by the use of methods of finite element (FEM) and finite volume (FVM) and k–ε model. In order to evaluate and verify this turbulent model, the results of the numerical model have been compared with the experimental results. This model has been studied for consideration of various jet flow velocities and thicknesses. The conclusions have indicated that the length and the height of the penetration area have linear relationship with jet flow velocity; therefore, as the jet flow velocity increases, the length of jet penetration increases as well. The comparison of the results of numerical method with the experimental data have demonstrated that the FVM holds less convergence time and better results compared with FEM.     

Temperature-driven fluid flow in porous media using a mixed finite element method and a finite volume method

We present a numerical scheme for the computation of conservative fluid velocity, pressure and temperature fields in a porous medium. For the velocity and pressure we use the primal–dual mixed finite element method of Trujillo and Thomas while for the temperature we use a cell-centered finite volume method. The motivation for this choice of discretization is to compute accurate conservative quantities. Since the variant of the mixed finite element method we use is not commonly used, the numerical schemes are presented in detail. We sketch the computational details and present numerical experiments that justify the accuracy predicted by the theory.     

Coupled hydrodynamic response of concrete gravity dams using finite and infinite elements

This paper presents the application of the finite element method for analysing the two-dimensional response of reservoir-dam systems subjected to horizontal ground motion. The interaction between the dam and the reservoir as well as the compressibility of water has been taken into account. The complete system has been considered to be composed of two substructures, namely the reservoir and the dam. To take into account the large extent of the reservoir, it has been idealized using specially developed infinite elements coupled with standard finite elements while the dam is represented using finite elements alone. Structural damping of the dam and radiation damping in the fluid phase have been accounted for in the analysis. It is concluded that the effect of radiation damping is considerable at high frequencies of excitation. The coupled response of the system is significantly large at and near the fundamental natural frequency of the system in comparison to the uncoupled responses. The method is computationally quite economical, capable of taking into account the arbitrary geometry of the system and is recommended for practical application. Further applications and extensions of the approach to three dimensional analyses are possible.     

Numerical simulation of groundwater flow patterns using flux as upper boundary

Since the 1960s, most of the studies on groundwater flow systems by analytical and numerical modelling have been based on given‐head upper boundaries. The disadvantage of the given‐head approach is that the recharge into and discharge from a basin vary with changes in hydraulic conductivity and/or basin geometry. Consequently, flow patterns simulated with given‐head boundaries but with different hydraulic conductivities and/or basin geometry may not reflect the effects of these variables. We conducted, therefore, numerical simulations of groundwater flow in theoretical drainage basins using flux as the upper boundary and realistically positioned fluid‐potential sinks while changing the infiltration intensity, hydraulic conductivities, and geometric configuration of the basin. The simulated results demonstrate that these variables are dominant factors controlling the flow pattern in a laterally closed drainage basin. The ratio of infiltration intensity to hydraulic conductivity (R_ic) has been shown to be an integrated pattern‐parameter in a basin with a given geometric configuration and possible fluid‐potential‐sink distribution. Successively, the changes in flow patterns induced by stepwise reductions in R_ic are identical, regardless of whether the reductions are due to a decrease in infiltration intensity or an increase in hydraulic conductivity. The calculated examples show five sequential flow patterns containing (i) only local, (ii) local–intermediate, (iii) local–intermediate–regional, (iv) local–regional, and (v) just regional flow systems. The R_ic was found to determine also whether a particular sink is active or not as a site of discharge. Flux upper boundary is preferable for numerical simulation when discussing the flow patterns affected by a change of infiltration, the hydraulic conductivity, or the geometry of a basin. Copyright © 2012 John Wiley & Sons, Ltd.     

Separation of PP- and PS-wave reflected seismic data using two-dimensional finite offset common-reflection-surface traveltime approximation

Recently, the interest in PS-converted waves has increased for several applications, such as sub-basalt layer imaging, impedance estimates and amplitude-versus-offset analysis. In this study, we consider the problem of separation of PP- and PS-waves from pre-stacked multicomponent seismic data in two-dimensional isotropic medium. We aim to demonstrate that the finite-offset common-reflection-surface traveltime approximation is a good alternative for separating PP- and PS-converted waves in common-offset and common shot configurations by considering a two-dimensional isotropic medium. The five parameters of the finite-offset common-reflection-surface are firstly estimated through the inversion methodology called very fast simulated annealing, which estimates all parameters simultaneously. Next, the emergence angle, one of the inverted parameters, is used to build an analytical separation function of PP and PS reflection separation based on the wave polarization equations. Once the PP- and PS-converted waves were separated, the sections are stacked to increase the signal-to-noise ratio using the special curves derived from finite-offset common-reflection-surface approximation. We applied this methodology to a synthetic dataset from simple-layered to complex-structured media. The numerical results showed that the inverted parameters of the finite offset common-reflection-surface and the separation function yield good results for separating PP- and PS-converted waves in noisy common-offset and common shot gathers.     

Symmetric modes of vibration for A circular plate using finite dynamic elements

The dynamic element method has been shown previously to provide a computational advantage over the ordinary finite element method for various beam elements. The Taylor expansions are computed here for the dynamic shape functions (two terms) and dynamic stiffness matrix (four terms) for the axisymmetric vibrations of an annular plate element. The complicated matrices which result are made more tractable by expressing them as power series in powers of the aspect ratio. The percentage error in the natural frequencies is then calculated using both the two- and the three-term dynamic stiffness matrix, demonstrating the increased accuracy for a given number of elements.     

Application of the iterative solution method with Schur complement reduction to mixed finite elements based in a tetrahedral discretization

The iterative solution method for mixed finite element methods is applied to a 3-D domain partitioned with tetrahedral elements. For the particular discretization technique of first partitioning the domain with hexahedral cells, and then subsequently partitioning cells with five tetrahedral elements, a Schur complement decomposition is devised wherein the actual number of equations solved is reduced by 80%. Although this Schur complement reduction requires a fair amount of computational overhead, its application within the iterative solution method can reduce overall solution time by about 44%, depending on closure criterion and other factors.     

Numerical determination of residence time and age in a partially mixed estuary using three-dimensional hydrodynamic model

This paper documents a numerical modeling study to calculate the residence time and age of dissolved substances in a partially mixed estuary. A three-dimensional, time-dependent hydrodynamic model was established and applied to the Danshuei River estuarine system and adjacent coastal sea in Taiwan. The model showed good agreement with observations of surface elevation, tidal currents and salinity made in 2002. The model was then applied to calculate the residence time and age distribution response to different freshwater discharges with and without density-induced circulations in the Danshuei River estuarine system. Regression analysis of model results reveals that an exponential equation can be used to correlate the residence time to change of freshwater input. The simulated results show it takes approximately 10, 4.5, and 3 days, respectively, for a water parcel that has entered the headwaters of the estuary to be transported out of the estuary under low, mean, and high flow conditions with density-induced circulation. The calculated age with density-induced circulation is less than that without density-induced circulation. The age of the surface layer is less than that at the bottom layer. Overall the study shows that freshwater discharges are the important factors in controlling the transport of dissolved substances in the Danshuei River estuarine system.     

Stochastic analysis of flux and head moments in a heterogeneous aquifer system

This study investigates the behavior of flux and head in a strongly heterogeneous three-dimensional aquifer system. The analyses relied on data from 520 slug tests together with 38,000 one-foot core intervals lithological data from the site of the General Separations Area in central Savannah River Site, South Carolina, USA. The skewness in the hydraulic conductivity histograms supported the geologic information for the top two aquifers, but revealed stronger clay content, than was reported for the bottom aquifer. The log-normal distribution model described adequately the hydraulic conductivity measurements for all three aquifers although, other distributions described equally well the bottom aquifer measurements. No apparent anisotropy on the horizontal plane was found for the three aquifers, but ratios of horizontal to vertical correlation lengths between 33 and 75 indicated a strong stratification at the site. Three-dimensional Monte Carlo stochastic simulations utilized a grid with larger elements than the support volume of measurements, but of sub-REV (representative elementary volume) dimensions. This necessitated, on one hand, the use of upscaled hydraulic conductivity expressions, but on the other hand did not allow for the use of anisotropic effective hydraulic conductivity expressions (Sarris and Paleologos in J Stoch Environ Res Risk Assess 18: 188–197, 2004). Flux mean and standard deviations components were evaluated on three vertical cross-sections. The mean and variance of the horizontal flux component normal to a no-flow boundary tended to zero at approximately two to three integral scales from that boundary. Close to a prescribed head boundary both the mean and variance of the horizontal flux component normal to the boundary increased from a stable value attained at a distance of about five integral scales from that boundary. The velocity field 〈q_x〉 was found to be mildly anisotropic in the top two aquifers, becoming highly anisotropic in the bottom aquifer; 〈q_y〉 was anisotropic in all three aquifers with directions of high continuity normal to those of the 〈q_x〉 field; finally, 〈q_z〉 was highly anisotropic in all three aquifers, with higher continuity along the east–west direction. The mean head field was found to be continuous, despite the high heterogeneity of the underlying hydraulic conductivity field. Directions of high continuity were in alignment with field boundaries and mean flow direction. Conditioning did not influence significantly the expected value of the flux terms, with more pronounced being the effect on the standard deviation of the flux vector components. Conditioning reduced the standard deviations of the horizontal flux components by as much as 50% in the bottom aquifer. Variability in the head cross-sections was affected only marginally, with an average 10% reduction in the respective standard deviation. Finally, the location of the conditioning data did not appear to have a significant effect on the surrounding area, with uniform reduction in standard deviations.     

Coupling of boundary and finite elements for soil-structure interaction problems

The investigation of complex soil-structure interaction problems is usually carried out with numerical solution procedures such as the finite element or the boundary element method. It must be noted, however, that the choice of one or the other of these approaches is not just a matter of preferences; depending on the type of the problem under consideration, either boundary or finite elements may be more advantageous. A considerable expansion in the computational power can be obtained, on the other hand, if one resorts to hybrid schemes which retain the main advantages of the two methods and eliminate their respective disadvantages. This paper presents results obtained with a boundary element-finite element coupling procedure, and discusses its applicability to some representative soil-structure interaction problems. The structures considered are elastic systems, such as foundations, tunnels and filled trenches (modelled by finite elements), which are coupled with homogeneous elastic halfspaces (modelled by boundary elements). The examples demonstrate the importance of using a model that includes wave radiation effects. The coupling approach is formulated entirely in the time domain so that an extension of the algorithm to non-linear analyses seems to present no further difficulties.     

Dynamic response of flexible strip-foundations by boundary and finite elements

A time domain Boundary Element-Finite method is employed to determine the dynamic response of flexible surface two-dimensional foundations under conditions of plane strain placed on an elastic soil medium and subjected either to transient external forces or to obliquely incident seismic waves. The elastic, isotropic, and homogeneous soil medium is treated by the time domain Direct Boundary Element Method, while the flexible foundation is treated by the Finite Element Method. The two methods are appropriately combined through equilibrium and compatibility considerations at the soil-foundation interface. Parametric studies examining the effect of the relative stiffness between the foundation and the soil and the spatial distribution of the dynamic disturbances on the foundation response are presented.     

On the propagation of Rayleigh waves in three dimensions in alluvial soils

Summary The propagation of Rayleigh waves in three dimensions in alluvial soils which do not behave like ordinary isotropic elastic solids have been discussed in this paper. The frequency equation has been solved for different soil constants.     

Coupling method of finite and infinite elements for strip foundation wave problems

Starting from two-dimensional wave equations and making use of Galerkin weighted residual approximations, discretized formulations for wave problems of a visco-elastic foundation have been derived. With considerations of geometrical and mechanical characteristics of a semi-infinite domain, a frequency-dependent compatible infinite element is also presented. Finally, by coupling the infinite elements with ordinary finite elements the system is used for simulation of propagating waves in a semi-infinite foundation. This model is not only suitable for simulations of complicated variations of geometrical conditions, but also for describing the unbounded behaviour of arbitrary multiple layers. Examples given indicate the model has excellent computational accuracy and feasibility for analysing the effects on foundation response due to the existence of faults or any other soft layers.     

Distribution of trace elements in subcellular fractions of three aquatic birds

This study investigated the subcellular distribution of Cu, Zn, Cd and Ag in liver of mallard (Anas platyrhynchos), spot-billed duck (A. poecilorhyncha) and great cormorant (Phalacrocorax carbo) to better understand metal fractionation and dynamics in avian liver. Most of the total Cu, Zn, Cd and Ag were present in hepatocytosol, and their concentrations increased with total hepatic levels in all the three avian species. Copper, Zn and Cd in hepatocytosol were present mostly in metallothionein fractions (MTs), and Cu and Cd contents in this fraction were positively correlated with hepatocytosol levels in all the species. Silver was observed mostly in both high-molecular weight and MT fractions in hepatocytosol in mallard and great cormorant, whereas it was present in both low-molecular weight and MT fractions in spot-billed duck, suggesting that distribution of Ag in the hepatocytosol was species-specific. The elution profile of metals in MT fractions revealed six metallic peaks in mallard and spot-billed duck, and three peaks in great cormorant, implying the presence of multiple MT isoforms in the liver of these avian species. The present study demonstrated that MTs are closely associated with metal regulation, especially Cu and Cd, in these three aquatic birds.     

Coupling of one- and two-dimensional finite elements for the computation of tidal flow in estuaries

A new method for the numerical computation of tidal flow in estuaries is demonstrated. One- and two-dimensional finite elements are directly coupled in a single model requiring the solution of only one common set of equations for all unknown parameters in each time step. The utility of the method is demonstrated with the computation of the tide in the Elbe Estuary. The computed water levels are compared with measured hydrographs at different gauges on the Elbe River.