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.     

Non-linear finite element dynamic analysis of axisymmetric solids

The subject of this paper is the finite element analysis of axisymmetric solids subjected to axisymmetric static and dynamic loading, and taking into account material as well as geometric non-linearities. A general Lagrangian formulation forms the basis for the incremental equations of motion which are solved using direct integration methods. Solution accuracy is improved by applying equilibrium correction loads at each step. Finite element discretization is achieved through the use of quadrilateral plane stress and axisymmetric elements with incompatible modes added for improvement of the element flexural characteristics. Several numerical examples are presented to demonstrate the effectiveness of the developed computer program.     

Time domain dynamic analysis of piles under impact loading

In this paper, time domain dynamic analysis of piles under impact loading is presented. For this purpose a hybrid boundary element technique is implemented. Linear beam column finite elements are used to model the piles and resulting governing equations are solved using an implicit integration scheme. The continuum is assumed to be elastic and an efficient step-by-step time integration scheme is implemented by using an approximate half space integral formulation. By enforcing displacement equilibrium conditions at each time step, a system of equations is generated which yields the solution. Results of this time domain formulation under linear material behavior are compared with Laplace domain results to validate the methods.     

A finite element-finite difference alternating direction algorithm for three-dimensional groundwater transport

While technically feasible, three-dimensional finite element groundwater transport simulation has not found widespread application because of the considerable computational burden inherent in the approach. An operator splitting algorithm which treats the horizontal plane using finite elements in the first step and the vertical dimension using finite differences in the second step provides considerable savings in both computer memory requirements and computational effort. The algorithm reduces the problem of an (N × M) array of nodes from one in which NM equations must be solved simultaneously, to one of solving N equations M times in the first step and M equations N times in the second step. Applications to type problems and field situations indicate that the method is robust and accurate.     

THREE-DIMENSIONAL NUMERICAL MODELING OF COHESIVE SEDIMENT TRANSPORT BY TIDAL CURRENT IN PEARL RIVER ESTUARV

IINTRODUCTIONEstUariesareprominentcoastalfeatUres.Estuariesareofgreateconomicssignificancetomankind.Attheseareas,manyharborsandwaterchannelshavetobebuiltforeconomicpurposes.ThedesignandconstrUctionofcoastalstrUctUresinestUariesrequireknowledgeofhydrodynamicsaswellassedimenttransportinsuchregions.ThenatUreofestuariesiscontrolledbyvariouscoastalhydrodynamicprocesses.Undertheactionofhydrodynamics,sedimentdepositionsorerosionswilloccurinestuariesornearcoastalstrUCtures.Tomaintainnavigati…     

钢框架结构的非线性静力抗震可靠性分析:Ⅰ有限元反应及其灵敏度

作为一种随机有限元方法,有限元可靠度方法通过有限元反应灵敏度分析将结构可靠度分析的近似解析方法与结构确定性分析的有限元方法结合起来,可以有效地处理结构反应是基本随机变量的隐式函数这一难题,因此成为大型复杂结构可靠度分析的有效工具。采用基于位移的非线性纤维梁柱单元对钢框架结构进行有限元建模,推导了单元与截面的基本方程。针对基于位移的纤维梁柱单元,采用逐级递进方式,分别推导了整体级、单元级、截面级和材料级的有限元反应灵敏度直接微分表达式。研究结果可为采用FORM或SORM等近似解析法进行钢框架结构的非线性静力有限元可靠性分析提供算法支持和编程依据。     

Parallel iterative solution for h and p approximations of the shallow water equations

A p finite element scheme and parallel iterative solver are introduced for a modified form of the shallow water equations. The governing equations are the three-dimensional shallow water equations. After a harmonic decomposition in time and rearrangement, the resulting equations are a complex Helmholz problem for surface elevation, and a complex momentum equation for the horizontal velocity. Both equations are nonlinear and the resulting system is solved using the Picard iteration combined with a preconditioned biconjugate gradient (PBCG) method for the linearized subproblems. A subdomain-based parallel preconditioner is developed which uses incomplete LU factorization with thresholding (ILUT) methods within subdomains, overlapping ILUT factorizations for subdomain boundaries and under-relaxed iteration for the resulting block system. The method builds on techniques successfully applied to linear elements by introducing ordering and condensation techniques to handle uniform p refinement. The combined methods show good performance for a range of p (element order), h (element size), and N (number of processors). Performance and scalability results are presented for a field scale problem where up to 512 processors are used.     

A new finite element technique for the solution of two-phase flow through porous media

A new upstream weighting finite element technique is developed for improved solution of the two-phase immiscible flow equations. Unlike the upstream weighting technique used by previous investigators, the new technique does not employ finite difference concepts to achieve the required upstream weighting of relative permeabilities or mobilities. Instead, upstream weighting is achieved by (1) representing the relative permeabilities or mobilities as continuous functions expressed in terms of the shape functions and nodal values (2) using asymmetric weighting functions to weight the spatial terms in the flow equations. These weighting functions are constructed such that they are dependent on the flow direction along each side of an element.In conjunction with the proposed technique, two solution schemes for treating the resulting set of non-linear algebraic equations are presented. These are the fully-implicit chord slope incremental solution scheme and the Newton-Raphson solution scheme. Both schemes allow the use of large time steps without being unstable.The proposed numerical technique is applied to two problems (1) the one-dimensional Buckley-Leverett problem (2) the two-dimensional five-spot well flow problem. Results indicate that this technique is superior to not only earlier finite element schemes but also five-point upstream finite difference formulae.     

考虑流-固耦合影响的挡水结构自振特性

对流体-固体之间的动力相互作用的问题，采用流体有限元法研究结构自振特性的影响。建立了流固耦合运动方程体系，并对耦合方程组的求解过程进行了推导与阐述。编制了流体有限元自动生成程序，使得流固耦合问题从建立模型到分析能够利用计算机自动实现。最后以挡水结构作为算例，分别就不同水位工况和水体采用不同长度模型对结构振动的影响进行了分析。     

库水模型及库水深度对拱坝动力响应的影响

基于接触非线性有限元模型,以锦屏一级拱坝为例,库水分别采用附加质量模型、可压缩流体有限元模型、不可压缩流体有限元模型计算了正常蓄水位及运行低水位时坝体的动力响应,结果表明:库水模型对拱坝动力响应有较大影响,随库水深度的增大,各模型计算结果差异增大;相比于流体可压缩模型,采用不可压缩流体模型所得动力响应普遍偏大;运行低水位工况,由于静水压力减小导致拱效应减弱,从而降低了拱坝的整体性,因此运行低水位工况各缝开度普遍高于正常蓄水位工况,且其拉应力范围较大,因此,运行低水位工况将对抗震设计起控制作用。     

Time-stepping schemes for finite element tidal model computations

A number of temporal procedures for solving the long-wave surface water equations using the finite element method in space are presented and analyzed. The analysis determines the stability of the schemes and the error in wave amplitude and phase that can be expected. The computational efficiency of the various methods is also discussed. The results ofthis analysis indicate types of errors that might be manifested in finite element surface water modeling using the different schemes.     

A model for tides and currents in the English Channel and southern North Sea

The amplitude and phase of 11 tidal constituents for the English Channel and southern North Sea are calculated using a frequency domain, finite element model. The governing equations — the shallow water equations — are modifed such that sea level is calculated using an elliptic equation of the Helmholz type followed by a back-calculation of velocity using the primitive momentum equations. Triangular elements with linear basis functions are used. The modified form of the governing equations provides stable solutions with little numerical noise. In this field-scale test problem, the model was able to produce the details of the structure of 11 tidal constituents including O₁, K₁, M₂, S₂, N₂, K₂, M₄, MS₄, MN₄, M₆, and 2MS₆.     

频率域激电有限元数值模拟

在三维地电条件下,应用有限元方法模拟了不同频率谐变电流激发下的地表激电响应,首先根据麦克斯韦方程推导出谐变电磁场复电位的波动方程,利用广义变分原理得到复电位波动方程及其边界条件下边值问题等价的变分问题.尔后将区域剖分离散化,在单元中进行插值,得到一组线性方程组,通过合理存储刚度矩阵、解方程组,得到各个结点上的复电位值,最后得到表征频率域激电响应的幅频率等参数.模拟结果表明该方法合理正确.     

Stability and accuracy of two-dimensional kinematic wave overland flow modeling

A two-dimensional finite element based overland flow model was developed and used to study the accuracy and stability of three numerical schemes and watershed parameter aggregation error. The conventional consistent finite element scheme results in oscillations for certain time step ranges. The lumped and the upwind finite element schemes are tested as alternatives to the consistent scheme. The upwind scheme did not improve on the stability or the accuracy of the solution, while the lumped scheme provided stable and accurate solutions for time steps twice the size of time steps needed for the consistent scheme. A new accuracy based dynamic time step estimate for the two-dimensional overland flow kinematic wave solution is developed for the lumped scheme. The newly developed dynamic time step estimates are functions of the mesh size, and time of concentration of the watershed hydrograph. Due to lack of analytical solutions, the time step was developed by comparing numerical solutions of various levels of discretization to a reference solution using a very fine mesh and a very small time step. The time step criteria were tested on a different set of problems and proved to be adequate for accurate and stable solutions. A sensitivity analysis for the watershed slope, Manning’s roughness coefficient and excess rainfall rate was conducted in order to test the effect of parameter aggregation on the stability and accuracy of the solution. The results of this analysis show that aggregation of the slope data resulted in the highest error. The roughness coefficient had a smaller effect on the solution while the rainfall intensity did not show any significant effect on the flow rate solution for the range of rainfall intensity used. This work pioneers the challenge of providing guidelines for accurate and stable numerical solutions of the two-dimensional kinematic wave equations for overland flow.     

频率域线源大地电磁法有限元正演模拟

本文介绍了频率域线源大地电磁法有限元正演模拟的研究结果. 在外边界上统一应用适合于人工源的一阶吸收边界条件来形成边值问题,可减小基于平面波假设造成的人为截断边界的影响. 程序编辑中设计了两个二维数组分别存储总体系数矩阵的非零元素和在总体结点编号中的位置,使内存占用量减少,且物理意义明确,方便用高斯-赛德尔等迭代法解有限元方程时调用. 采用视δ函数模拟线源,提高了解方程组的稳定性. 最后通过对1个简单模型和1个复杂模型的模拟,证明所用的方法对异常体能够有明显的反映,说明了该方法的可靠性和有效性.     

A 3-D finite element conjugate gradient model of subsurface flow with automatic mesh generation

The 3-D flow modelling of groundwater systems of realistic size generally requires a big effort for the preparation of the input data as well as large computational costs. A numerical finite element model (MAITHREE) is developed for the efficient analysis of the steady and unsteady behaviour of natural confined 3-D basins. Starting from an initial triangular grid the code automatically generates a set of tetrahedral elements in each of the geologic units or subunits specified by the user in the vertical profile. The original element incidences list is then rearranged in order to provide conforming 3-D elements throughout the domain. The model is designed with a view to saving much of the labour involved in setting a 3-D grid and to providing flexibility as well as economical convenience through a high computational efficiency. The latter task is achieved by the aid of a solver based on the modified conjugate gradient (MCG) method which has proved to be an excellent technique for the solution of large linear finite element sets of sparse 3-D subsurface equations. Some examples derived from both hypothetical and real-world situations are discussed to illustrate the innovative features of MAITHREE and its computational performance.     

Coupled groundwater flow and transport: 1. Verification of variable density flow and transport models

This work examines variable density flow and corresponding solute transport in groundwater systems. Fluid dynamics of salty solutions with significant density variations are of increasing interest in many problems of subsurface hydrology. The mathematical model comprises a set of non-linear, coupled, partial differential equations to be solved for pressure/hydraulic head and mass fraction/concentration of the solute component. The governing equations and underlying assumptions are developed and discussed. The equation of solute mass conservation is formulated in terms of mass fraction and mass concentration. Different levels of the approximation of density variations in the mass balance equations are used for convection problems (e.g. the Boussinesq approximation and its extension, fully density approximation). The impact of these simplifications is studied by use of numerical modelling.Numerical models for nonlinear problems, such as density-driven convection, must be carefully verified in a particular series of tests. Standard benchmarks for proving variable density flow models are the Henry, Elder, and salt dome (HYDROCOIN level 1 case 5) problems. We studied these benchmarks using two finite element simulators - ROCKFLOW, which was developed at the Institute of Fluid Mechanics and Computer Applications in Civil Engineering and FEFLOW, which was developed at the Institute for Water Resources Planning and Systems Research Ltd. Although both simulators are based on the Galerkin finite element method, they differ in many approximation details such as temporal discretization (Crank-Nicolson vs predictor-corrector schemes), spatial discretization (triangular and quadrilateral elements), finite element basis functions (linear, bilinear, biquadratic), iteration schemes (Newton, Picard) and solvers (direct, iterative). The numerical analysis illustrates discretization effects and defects arising from the different levels of the density of approximation. We contribute new results for the salt dome problem, for which inconsistent findings exist in literature. Applications of the verified numerical models to more complex problems, such as thermohaline and three-dimensional convection systems, will be presented in the second part of this paper.     

大地电磁二维对称各向异性介质的有限元数值模拟

假定垂直轴为二维对称各向异性介质主轴之一，构造走向与另一主轴方向成任一夹角，用伽勒金（Ｇａｌｅｒｋｉｎ）有限元法和矩形网格，优化地合成总体刚度矩阵，使二维各向异性介质的基本方程形成有限元代数方程组，求出各节点场值，并利用ＭＯＭ法求出辅助场，进而求出张量阻抗等响应函数．同时还对有关文献的计算模型进行了数值模拟检验     

Fully coupled finite element model for dynamics of partially saturated soils

Understanding the response of partially saturated earth structures under various static and dynamic loads is important for the design and construction of economical and safe geotechnical engineering structures. In this study, the numerical approach is used to understand the dynamics of partially saturated soils. The mathematical equations governing the dynamics of partially saturated soils are derived based on the theory of mixtures and implemented within a finite element framework. The stress–strain behavior of the soil is represented by an elasto-plastic constitutive model for unsaturated soil based on bounding surface concept and the moisture-suction behavior is modeled using van Genuchten model. Fully coupled finite element simulations are performed to study the response of partially saturated soil embankment under earthquake loading and validated with centrifuge test results available in the literature. The predicted displacement responses are in good agreement with the measured responses. The pore water pressure, pore air pressure, matric suction, the degree of saturation in various elements and the response of the embankment under different initial moisture content are also discussed.