On the moments of time to first passage of the linear oscillator

The moments of time to first passage of the response process of the linear oscillator excited by a stationary wide-band process are obtained numerically for the symmetric two-sided barrier problem. A recursive set of partial differential equations governing the moments is solved using a Petrov-Galerkin finite element method. Stable, accurate solutions for the first few moments are calculated over a wide range of oscillator damping.     

Lanczos method for dynamic analysis of damped structural systems

In this paper we extend the Lanczos algorithm for the dynamic analysis of structures⁷ to systems with general matrix coefficients. The equations of dynamic equilibrium are first transformed to a system of first order differential equations. Then the unsymmetric Lanczos method is used to generate two sets of vectors. These vectors are used in a method of weighted residuals to reduce the equations of motion to a small unsymmetric tridiagonal system. The algorithm is further simplified for systems of equations with symmetric matrices. By appropriate choice of the starting vectors we obtain an implementation of the Lanczos method that is remarkably close to that in Reference 7, but generalized to the case with indefinite matrix coefficients. This simplification eliminates one of the sets of vectors generated by the unsymmetric Lanczos method and results in a symmetric tridiagonal, but indefinite, system. We identify the difficulties that may arise when this implementation is applied to problems with symmetric indefinite matrices such as vibration of structures with velocity feedback control forces which lead to symmetric damping matrices. This approach is used to evaluate the vibration response of a damped beam problem and a space mast structure with symmetric damping matrix arising from velocity feedback control forces. In both problems, accurate solutions were obtained with as few as 20 Lanczos vectors.     

Comparison of Galerkin-type discretization techniques for two-phase flow in heterogeneous porous media

A comparison of Standard Galerkin, Petrov-Galerkin, and Fully-Upwind Galerkin methods for the simulation of two-phase flow in heterogeneous porous media is presented. On the basis of the coupled pressure-saturation equations, a generalized formulation for all three finite element methods is derived and analysed. For flow in homogeneous media, the Petrov-Galerkin method gives excellent results. But this method fails miserably for problems with heterogeneous media. This is because it is not able to capture correctly processes that take place at interfaces when, for instance, the capillary pressure-saturation relationship after Brooks and Corey is assumed. The Fully-Upwind Galerkin method is superior to the Petrov-Galerkin approach because it is able to give correct results for flow in homogeneous and heterogeneous media for the two models of van Genuchten and Brooks-Corey. The widely used formulation which is correct for the homogeneous case cannot be used for heterogeneous media. Instead the straightforward approach of gradp_c in combination with a chord-slope technique must be utilized.     

An approximate Riemann solver for sensitivity equations with discontinuous solutions

The sensitivity of a model output (called a variable) to a parameter can be defined as the partial derivative of the variable with respect to the parameter. When the governing equations are not differentiable with respect to this parameter, problems arise in the numerical solution of the sensitivity equations, such as locally infinite values or instability. An approximate Riemann solver is thus proposed for direct sensitivity calculation for hyperbolic systems of conservation laws in the presence of discontinuous solutions. The proposed approach uses an extra source term in the form of a Dirac function to restore sensitivity balance across the shocks. It is valid for systems such as the Euler equations for gas dynamics or the shallow water equations for free surface flow. The method is first detailed and its application to the shallow water equations is proposed, with some test cases such as dike- or dam-break problems with or without source terms. An application to a two-dimensional flow problem illustrates the superiority of direct sensitivity calculation over the classical empirical approach.     

对黄河流域水资源保护的探索与实践

黄河水资源贫乏,水污染严重,水资源供需矛盾日益尖锐,黄河流域水资源保护面临着诸多的困难和压力.新《水法》的颁布,给黄河水资源保护带来新的机遇和挑战,依法建立法制化、社会化的流域与区域相结合的水资源保护管理体系,完善法规,建立健全执法体系,以水功能区划管理为重心,实施入河总量控制,建立和完善重大水污染事件快速反应机制,加强水质监测机制与技术创新,加强水资源保护前期和科研工作,建设水资源保护信息管理系统,采取多部门联合治污等项措施,是做好黄河流域水资源保护工作的重要前提和任务.     

COMPUTING WELL HYDRAULICS SOLUTIONS

Abstract. Type curves, and drawdown analysis curves of other sorts, have a well-established position in the ground-water profession. Although a plethora of solutions to well hydraulics problems exists, many of these arc not available in useful type curve formats. lpurthcr, if one wants to perform parameter estimation using many of these solutions, a computer-based form of the solution is needed. Those who have programmed analytical solutions know the difficulties that will be encountered in coding even the most innocuous-looking result. In this paper, a computer-based method of type curve generation is prcscnted based on Stehfest's method of inverting the Laplace transform solution of ground-water flow problems. The latter are usually easier to obtain than a complete closed-form solution. The method is very robust for a large class of well hydraulics problems. A companion diskette with FORTRAN source and executable example program is available from the author.     

非均匀介质热蠕变流动的数值求解

针对非均匀介质中热蠕变流动问题,给出了有限单元方法与网格－粒子方法联合求解新技术,即有限单元方法求解欧拉网格节点上的未知量,分布于单元内部作为物质成分标记的粒子反映变形过程.有限元法求解动量方程和连续性方程时引入了速度场和压力场等阶插值的压力场稳定的Petrov Galerkin方法,求解能量方程时采用了流线迎风Petrov Galerkin方法,网格－粒子算法中采用双线性插值与有限单元插值函数对应.有限单元计算与网格－粒子计算相对独立,两种方法计算的数据通过有限单元节点传递.同时,实现了三角形单元的算法和程序,解决了复杂结构条件下不规则网格计算的问题.通过经典方腔热对流问题验证了程序,给出了不规则形态块体沉降算例,并分析了数值解的稳定性.     

Frequency‐domain waveform modelling and inversion for coupled media using a symmetric impedance matrix

To simulate the seismic signals that are obtained in a marine environment, a coupled system of both acoustic and elastic wave equations is solved. The acoustic wave equation for the fluid region simulates the pressure field while minimizing the number of degrees of freedom of the impedance matrix, and the elastic wave equation for the solid region simulates several elastic events, such as shear waves and surface waves. Moreover, by combining this coupled approach with the waveform inversion technique, the elastic properties of the earth can be inverted using the pressure data obtained from the acoustic region. However, in contrast to the pure acoustic and elastic cases, the complex impedance matrix for the coupled media does not have a symmetric form because of the boundary (continuity) condition at the interface between the acoustic and elastic elements. In this study, we propose a manipulation scheme that makes the complex impedance matrix for acoustic–elastic coupled media to take a symmetric form. Using the proposed symmetric matrix, forward and backward wavefields are identical to those generated by the conventional approach; thus, we do not lose any accuracy in the waveform inversion results. However, to solve the modified symmetric matrix, LDLT factorization is used instead of LU factorization for a matrix of the same size; this method can mitigate issues related to severe memory insufficiency and long computation times, particularly for large‐scale problems.     

On the primary variable switching technique for simulating unsaturated–saturated flows

Primary variable switching appears as a promising numerical technique for variably saturated flows. While the standard pressure-based form of the Richards equation can suffer from poor mass balance accuracy, the mixed form with its improved conservative properties can possess convergence difficulties for dry initial conditions. On the other hand, variable switching can overcome most of the stated numerical problems. The paper deals with variable switching for finite elements in two and three dimensions. The technique is incorporated in both an adaptive error-controlled predictor–corrector one-step Newton (PCOSN) iteration strategy and a target-based full Newton (TBFN) iteration scheme. Both schemes provide different behaviors with respect to accuracy and solution effort. Additionally, a simplified upstream weighting technique is used. Compared with conventional approaches the primary variable switching technique represents a fast and robust strategy for unsaturated problems with dry initial conditions. The impact of the primary variable switching technique is studied over a wide range of mostly 2D and partly difficult-to-solve problems (infiltration, drainage, perched water table, capillary barrier), where comparable results are available. It is shown that the TBFN iteration is an effective but error-prone procedure. TBFN sacrifices temporal accuracy in favor of accelerated convergence if aggressive time step sizes are chosen.     

Equilibrium-type and link-type lattice Boltzmann models for generic advection and anisotropic-dispersion equation

We extend lattice Boltzmann (LB) methods to advection and anisotropic-dispersion equations (AADE). LB methods are advocated for the exactness of their conservation laws, the handling of different length and time scales for flow/transport problems, their locality and extreme simplicity. Their extension to anisotropic collision operators (L-model) and anisotropic equilibrium distributions (E-model) allows to apply them to generic diffusion forms. The AADE in a conventional form can be solved by the L-model. Based on a link-type collision operator, the L-model specifies the coefficients of the symmetric diffusion tensor as linear combination of its eigenvalue functions. For any type of collision operator, the E-model constructs the coefficients of the transformed diffusion tensors from linear combinations of the relevant equilibrium projections. The model is able to eliminate the second order tensor of its numerical diffusion. Both models rely on mass conserving equilibrium functions and may enhance the accuracy and stability of the isotropic convection–diffusion LB models.The link basis is introduced as an alternative to a polynomial collision basis. They coincide for one particular eigenvalue configuration, the two-relaxation-time (TRT) collision operator, suitable for both mass and momentum conservation laws. TRT operator is equivalent to the BGK collision in simplicity but the additional collision freedom relates it to multiple-relaxation-times (MRT) models. “Optimal convection” and “optimal diffusion” eigenvalue solutions for the TRT E-model allow to remove next order corrections to AADE. Numerical results confirm the Chapman–Enskog and dispersion analysis.     

Finite element techniques for convective-diffusive transport in porous media

Operator-splitting techniques are applied to convective-diffusive transport problems in porous media. The convection is treated by applying a modified method of characteristics to time-step along the characteristics of the convective part of the flow. The nonsymmetry in the spatial operator is addressed via a Petrov-Galerkin method which uses a test function to achieve stability through a balancing of the remaining convection, the diffusion, and any possible reaction terms. The use of time-stepping along characteristics allows the use of large time-steps in a stable but accurate fashion. If local phenomena are important, self-adaptive local grid refinement techniques can be coupled with the operator splitting.     

Use of Picard and Newton iteration for solving nonlinear ground water flow equations

This study examines the use of Picard and Newton iteration to solve the nonlinear, saturated ground water flow equation. Here, a simple three-node problem is used to demonstrate the convergence difficulties that can arise when solving the nonlinear, saturated ground water flow equation in both homogeneous and heterogeneous systems with and without nonlinear boundary conditions. For these cases, the characteristic types of convergence patterns are examined. Viewing these convergence patterns as orbits of an attractor in a dynamical system provides further insight. It is shown that the nonlinearity that arises from nonlinear head-dependent boundary conditions can cause more convergence difficulties than the nonlinearity that arises from flow in an unconfined aquifer. Furthermore, the effects of damping on both convergence and convergence rate are investigated. It is shown that no single strategy is effective for all problems and how understanding pitfalls and merits of several methods can be helpful in overcoming convergence difficulties. Results show that Picard iterations can be a simple and effective method for the solution of nonlinear, saturated ground water flow problems.     

无网格局部Petrov-Galerkin法大地电磁场二维正演模拟

有限差分法和有限单元法在大地电磁场数值模拟中已经得到了广泛的应用,但其数值结果的精度在很大程度上依赖于网格的离散程度.当模拟起伏地形、弯曲界面等复杂地电模型大地电磁场响应时,常常需要花费大量的时间以便得到较合理的离散网格.无网格局部Petrov-Galerkin法(MLPG)不同于有限差分法和有限元法,其形函数和权函数脱离了网格的束缚.本文详细推导了二维大地电磁场边值问题的弱式形式,并将其离散为局部积分域内的表达形式.通过模拟二维海洋地电模型大地电磁场响应,并与结构网格有限元结果进行对比,验证了本文算法和程序的正确性及精度.设计了一个含有弯曲界面的二维地电模型,讨论了不同离散网格对MLPG无网格法模拟结果的影响,并与结构有限元法结果进行了比较,结果表明MLPG无网格法模拟结果受离散网格影响较小.最后利用MLPG无网格法计算了两个海洋起伏地形模型的大地电磁响应,讨论了海底起伏地形对大地电磁响应的影响.     

Tunnel formation processes in the Riverina area of N.S.W., Australia

Soil and site characteristics associated with tunnel erosion are diverse. This contributes to difficulties in identifying areas of potential tunnelling and the reclamation of land affected by this form of water erosion. Soil, site, and tunnel characteristics at three locations having different morphological expressions of this erosion form were examined to determine which features were important in tunnel development in the Riverina Area of N.S.W., Australia. Soil dispersion was a dominant process at each site and is related to high levels of exchangeable sodium. Field and laboratory studies demonstrated the importance of dispersion and its effect on water movement, particle transport, and subsequent tunnel initiation. Tunnel initiation mechanisms differed between the sites but could be identified by morphological examination of the tunnels. Potential tunnel eroding sites are characterized by a dispersible layer, a gradient for water flow, and soil characteristics that promote water Sow in defined paths. Appropriate control techniques need to be formulated for each differing form of tunnel erosion.     

复杂二维/三维大地电磁的有限单元法正演模拟策略

复杂二维和三维大地电磁模型的正演数值模拟具有一定的挑战性。对于复杂的二维和三维大地电磁正演问题,我们采用有限单元法进行求解。有限单元法最后形成一个线性方程组,系数矩阵是大型稀疏的带状对称复系数矩阵,并且其条件数远大于1,为严重病态矩阵,求解其对应方程组会遇到很多困难。不完全LU分解处理的Bi-CGSTAB迭代方法可用于该线性方程组的求解,并且具有速度快、精度高和稳定性好等优点;为了模拟无穷远边界及满足计算机的内存需求,在保证计算精度的情况下设计了非均匀网格剖分;在程序编制中,只存储有限元系数矩阵的非零元素,大大减少了正演计算的时间。通过对二维和三维模型电磁响应的计算,验证了算法的正确性。     

Coupled Aquifer-Borehole Simulation

A model coupling fluid hydraulics in a borehole with fluid flow in an aquifer is developed in this paper. Conservation of momentum is used to create a one-dimensional steady-state model of vertical flow in an open borehole combined with radially symmetric flow in an aquifer and with inflow to the well through the wellbore screen. Both laminar and turbulent wellbore conditions are treated. The influence of inflow through the wellbore screen on vertical flow in the wellbore is included, using a relation developed by Siwoń (1987) . The influence of inflow reduces the predicted vertical variation in head up to 15% compared to a calculation of head losses due to fluid acceleration and the conventional Colebrook-White formulation of friction losses in a circular pipe. The wellbore flow model is embedded into the MODFLOW-2000 ground water flow code. The nonlinear conservation of momentum equations are iteratively linearized to calculate the conductance terms for vertical flow in the wellbore. The resulting simulations agree favorably with previously published results when the model is adjusted to meet the assumptions of the previous coupled models.     

Numerical treatment of the kinematic dynamo model

Summary Since the initial equations are complicated, the treatment of the kinematic dynamo model requires the use of numerical methods. In applying them to the given problem difficulties are encountered, which are not easy to overcome. This paper deals with the analysis of the experience acquired in treating the model of a nearly symmetric dynamo. Three different methods were employed (stationary, oscillatory and general non-stationary), because a combination of several solutions will yield more comprehensive information about the model being studied. Although the results are based on the study of a single particular model, similar problems also occur in other excercises and, therefore, the conclusions have a more general validity.     

Evaluation of plan configuration irregularity effects on seismic response demands of L-shaped MRF buildings

Damage assessments after past earthquakes have frequently revealed that plan configuration irregular buildings have more severe damage due to excessive torsional responses and stress concentration than regular buildings. The plan configuration irregularities introduce major challenges in the seismic design of buildings. One such form of irregularity is the presence of re-entrant corners in the L-shaped buildings that causes stress concentration due to sudden changes in stiffness and torsional response amplification; hence causes early collapse. A constructive research into re-entrant corner and torsional irregularity problems is essentially needed greater than ever. Therefore, the focus of this study is to investigate structural seismic response demands for the class of L-shaped buildings through evaluating the plan configuration irregularity of re-entrant corners and lateral–torsion coupling effects on measured seismic response demands. The measured responses include story drift, inter-story drift, story shear force, overturning moment, torsion moment at the base and over building height, and torsional irregularity ratio. Three dimensional finite element model for nine stories symmetric buildings as reference model is developed. In addition, six L-shaped building models are formulated with gradual reduction in the plan of the reference building model. The results prove that building models with high irregularity are more vulnerable due to the stress concentration and lateral torsional coupling behavior than that with regular buildings. In addition, the related lateral shear forces in vertical resisting elements located on the periphery of the L-shaped buildings could be significantly increased in comparison with the corresponding values for a symmetric building.     

Comparison of fully coupled approaches for approximating nonlinear transport and reaction problems

The purpose of this work is to evaluate the computational efficiency of fully coupled approaches for approximating a common class of nonlinear, two-phase advective–dispersive–reactive equations. The general problem considered includes homogeneous phase chemical kinetics, equilibrium interphase mass transfer, and rate-controlled interphase mass transfer––all of which may be nonlinear. Aspects of the problem investigated include discrete mass conservative formulations, temporal discretization approaches, and nonlinear equation solution methods. Their effect on computational efficiency is investigated through a series of numerical experiments using a nondimensional model problem. The effect of problem characteristics such as large sorption capacity, strong sorption nonlinearity, fast mass transfer, fast reactions, and strong diffusion is investigated. Comparisons of solution efficiency show that the optimal approach depends upon: (1) the characteristics of the problem considered, which may be described in a nondimensional form; and (2) the accuracy achieved in the solution. Results offer general guidance for selecting solution approaches for the class of problems investigated and introduce some new solution approaches to the water resources field that may be applicable to other problems.