基于Lamb问题的谱元法和有限元法模拟精度比较

在地震动数值模拟方法中,谱元法和有限元法是应用较广泛的两种方法。基于经典的Lamb问题模型,首先推导给出地表竖向位移的解析解答。然后分别利用常用的四阶谱元法和线性有限元法,模拟了地表脉冲力源作用下模型的位移响应。考虑有意义的最短波长内的采样点个数及单元高宽比的变化,对比了两种方法的模拟精度;结果表明:对于谱元法,观测点与波源之间需至少包含两个网格,在此条件下,最短波长内包含一个网格（最短波长内5个采样点）时,数值解与解析解的误差小于1%,已达很高的精度;对于有限元法,最短波长内需包含10个网格时才能达到这一精度。此外,在满足网格尺寸要求的前提下,单元水平向与垂直向尺寸的比值在1∶1到5∶1的范围内时,谱元法和有限元法的模拟精度均变化不大。因此,单位波长内采样点个数相同时,谱元法的模拟精度比有限元法高的多,同时,在一定范围内两种方法的模拟结果对于宽高比的变化不敏感。     

有限元三角网格中波动的频散分析及数值仿真

波动问题有限元离散后会引起数值误差, 数值频散的本质就是数值误差传播引起的非物理解. 数值频散不仅没有实际意义, 而且还会影响对真实波动现象的认识. 为厘清有限元三角网格中波动数值频散的影响因素, 本文推导了集中质量矩阵和一致质量矩阵的频散函数, 同时给出了组合质量矩阵的频散函数, 并对不同质量矩阵的数值频散进行了对比研究. 理论分析和数值计算结果表明: 有限元三角网格中波动的数值频散受网格布局、 波传播方向、 单元网格纵横比以及质量矩阵的影响; 一致质量矩阵的数值频散比集中质量矩阵更易受到波传播方向的影响; 不合理的三角网格单元会对数值相速度(数值频散)产生不良影响; 正三角网格中波动的数值频散几乎不受波传播方向的影响; 一致质量矩阵与集中质量矩阵的线性组合能够有效地压制数值频散.      

Vectorized simulation of groundwater flow and contaminant transport using analytic element method and random walk particle tracking

Groundwater contaminant transport processes are usually simulated by the finite difference (FDM) or finite element methods (FEM). However, they are susceptible to numerical dispersion for advection‐dominated transport. In this study, a numerical dispersion‐free coupled flow and transport model is developed by combining the analytic element method (AEM) with random walk particle tracking (RWPT). As AEM produces continuous velocity distribution over the entire aquifer domain, it is more suitable for RWPT than FDM/finite element methods. Using the AEM solutions, RWPT tracks all the particles in a vectorized manner, thereby improving the computational efficiency. The present model performs a convolution integral of the response of an impulse contaminant injection to generate concentration distributions due to a permanent contaminant source. The RWPT model is validated with an available analytical solution and compared to an FDM solution, the RWPT model more accurately replicates the analytical solution. Further, the coupled AEM‐RWPT model has been applied to simulate the flow and transport in hypothetical and field aquifer problems. The results are compared with the FDM solutions and found to be satisfactory. The results demonstrate the efficacy of the proposed method.     

Application of multiscale finite element method in the uncertainty qualification of large-scale groundwater flow

In this article, we discuss the application of multiscale finite element method (MsFEM) to groundwater flow in heterogeneous porous media. We investigate the ability of MsFEM in qualifying the flow uncertainty. Monte Carlo simulation is employed to implement the stochastic analysis, and MsFEM is used to avoid a full resolution to the spatial variable conductivity field. Large-scale flow with high variability is investigated by inspecting the single realization as well as the probability distribution functions of head and velocity. The numerical results show that the performance of MsFEM depends on the ratio between the correlation length and the coarse element size. An accurate prediction to the velocity requires a much lower ratio than the head. The MsFEM has different convergence rates for the head and the velocity, while the convergence rates do not deteriorate as the variance grows.     

A penalty finite difference model for Navier-Stokes flow problem in sedimentation basins

An explicit finite difference solution procedure for the 2D, unsteady, turbulent and incompressible Navier-Stokes flow problem is developed. The innovative feature of the numerical model is the use of a penalty function technique. To illustrate the effectiveness of the model, two numerically predicted flow patterns for laboratory size sedimentation tanks are compared to actual experimental results. The predicted flow pattern shows a satisfactory match with experimental data.     

青藏高原东缘中下地壳流与地壳变形

地壳缩短导致青藏隆升造山是普遍的认识.然而,在青藏东部,越来越多的观测数据和研究支持了中下地壳流与隆升造山的关系.目前,地壳缩短造山机制和中下地壳流造山机制仍然处于争论之中.本文建立了二维黏弹塑性有限元模型,模拟了龙门山断层带的多个地震循环的应变与变形,探讨了无与有中下地壳流情况下,地壳地表的位移、速度与变形的分布和演化;以及有中下地壳流情况下,不同流动范围、速度与黏度对模型结果的影响;并结合地形变观测数据的约束,推测了青藏东缘中下地壳流的流动状态.模拟结果显示,通过对比有和无中下地壳流的模拟结果,发现青藏东部震间的地表垂向速度在变形样式及数值上存在较大差异,即存在地壳流的地表垂向抬升速率显著大于无地壳流;震间在龙门山断层西侧附近产生的垂向凸状隆起随中下地壳流的速度、黏度及通道长度的变化而变化.此外,本文研究结果对青藏其他地区可能存在的地壳流的研究也具有一定的参考意义.     

On the mechanisms of the saltating motion of bedload

Saltation of sediment particles is an important pattern of bedload transport.Based on force analysis for sediment particles,a Lagrangian model was proposed for the saltating motion of bedload in river flows,which was then solved with numerical method.Simulation results on the saltating trajectories neglecting particle rotation and turbulence effects compare fairly well with experimental observations.The mean values of the saltation parameters (saltation height,length and velocity) also agree well with the previous experimental data.Based on the numerical results,regression equations for the dimensionless saltation height,length and velocity were presented.Using the numerically achieved characteristics of the sediment saltation,we also obtained mathematical expression for the sediment transport rate.The studies in this paper are significant for its contribution to mechanism of the bedload motion and the computation of sediment transport rate.     

Application of the finite element method to gravity data case study: Western Turkey

Gravity anomalies always include the total effects (combination of the structures which have different densities and depths) of the study area and beyond. And the well-known non-uniqueness of potential field modelling may lead to very different interpretation results.The finite element method (FEM), which has been used in potential field interpretation for decades, makes complex problems to be solved easily and accurately. The first step of FEM is to identify the elements and then to decide on the boundary of the solution space. In this step, the solution space is divided into elements. After determination of the geometrical structure of the solution space, the most suitable elements should be chosen for this geometrical structure. The agreement between the geometry and the elements is quite important for the convergence to the best possible solution.In this work, the methods of trend analysis, filtering, analytical continuation and FEM were applied to a theoretical model and to gravity data from western Turkey to produce the regional and the residual anomalies. The results were compared and it was found that the FEM produced more accurate results than other methods did.     

考虑海底地形的三维频率域可控源电磁响应有限体积法模拟

为了提高海洋可控源电磁的地形响应模拟精度,探讨了三维交错网格剖分有限体积法求解频率域电磁场解的数值方法.在似稳场近似条件下,推导了有限差分法和有限体积法低频电磁场控制方程的离散表达式.通过对离散表达式的分析比较表明,有限体积法较之有限差分具有更高的计算精度,同时也具有与有限差分相当的计算效率.对含地形界面网格单元电导率采用加权平均处理,通过与二维有限元程序的计算结果对比,验证了该方法可以有效提高有限体积法对地形变化的模拟精度.应用有限体积法计算了一个三维带地形储层模型的可控源电磁响应,分析表明地形变化对电场分量影响明显,磁场分量对地形变化不敏感.     

修正辛格式有限元法的地震波场模拟

三角网格有限元法具有网格剖分的灵活性,能有效模拟地震波在复杂介质中的传播.但传统有限元法用于地震波场模拟时计算效率较低,消耗较大计算资源.本文采用改进的核矩阵存储(IKMS)策略以提高有限元法的计算效率,该方法不用组合总体刚度矩阵,且相比于常规有限元法节省成倍的内存.对于时间离散,将有限元离散后的地震波运动方程变换至Hamilton体系,在显式二阶辛Runge-Kutta-Nystr9m(RKN)格式的基础之上加入额外空间离散算子构造修正辛差分格式,通过Taylor展开式得到具有四阶时间精度时间格式,且辛系数全为正数.本文从理论上分析了时空改进方法相比传统辛-有限元方法在频散压制、稳定性提升等方面的优势.数值算例进一步证实本方法具有内存消耗少、稳定性强和数值频散弱等优点.     

Reinvestigation on mixing length in an open channel turbulent flow

The present study proposes a model on vertical distribution of streamwise velocity in an open channel turbulent flow through a newly proposed mixing length, which is derived for both clear water and sediment-laden turbulent flows. The analysis is based on a theoretical consideration which explores the effect of density stratification on the streamwise velocity profile. The derivation of mixing length makes use of the diffusion equation where both the sediment diffusivity and momentum diffusivity are taken as a function of height from the channel bed. The damping factor present in the mixing length of sediment-fluid mixture contains velocity and concentration gradients. This factor is capable of describing the dip-phenomenon of velocity distribution. From the existing experimental data of velocity, the mixing length data are calculated. The pattern shows that mixing length increases from bed to the dip-position, having a larger value at dip-position and then decreases up to the water surface with a zero value thereat. The present model agrees well with these data sets and this behavior cannot be described by any other existing model. Finally, the proposed mixing length model is applied to find the velocity distribution in wide and narrow open channels. The derived velocity distribution is compared with laboratory channel data of velocity, and the comparison shows good agreement.     

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.     

Numerical and Physical Modelling of the Effect of a Dam on Powder Avalanche Motion: Comparison with Previous Approaches

As part of the European project CADZIE, the interactionbetween powder avalanches and dams has been studied. First, we simulated a powder avalancheusing a heavy salt solution in a water tank. These experiments showed that (1) the ratio betweenthe maximum velocity parallel to the slope and the front velocity is greater than 1 (up to 1.5),and that (2) the velocity normal to the slope can rise to 74% of the front velocity. The presence ofa dam leads to a reduction of front velocity. The ratio between the maximum horizontal velocityand the front velocity increases with the dam height and reaches a maximum value of 1.9. Aswe did not have access to the density distribution in our experimental set-up, we calibrateda numerical model thanks to the experiments and we deduced the dynamic pressure. Our resultsshowed that, for an avalanche without obstacle, AVAER (Rapin, 1995) and Beghin and Closet'shypothesis (1990) lead to an overestimation of the dynamic pressure in the lower partof the flow as well as in the upper part, and that the use of front velocity for estimating the lossof kinetic energy leads to an overestimation of the dam's effectiveness.     

Adaptive finite element simulation of Stokes flow in porous media

The Stokes problem describes flow of an incompressible constant-viscosity fluid when the Reynolds number is small so that inertial and transient-time effects are negligible. The numerical solution of the Stokes problem requires special care, since classical finite element discretization schemes, such as piecewise linear interpolation for both the velocity and the pressure, fail to perform. Even when an appropriate scheme is adopted, the grid must be selected so that the error is as small as possible. Much of the challenge in solving Stokes problems is how to account for complex geometry and to capture important features such as flow separation. This paper applies adaptive mesh techniques, using a posteriori error estimates, in the finite element solution of the Stokes equations that model flow at pore scales. Different selected numerical test cases associated with various porous geometrics are presented and discussed to demonstrate the accuracy and efficiency of our methodology.     

Modelling estuarine and coastal flows using an unstructured triangular finite volume algorithm

Details are given of the development and application of a numerical model for predicting free-surface flows in estuarine and coastal basins using the finite volume method. Both second- and third-order accurate and oscillation free explicit numerical schemes have been used to solve the shallow water equations. The model deploys an unstructured triangular mesh and incorporates two types of mesh layouts, namely the ‘cell centred’ and ‘mesh vertex’ layouts, and provides a powerful mesh generator in which a user can adjust the mesh-size distribution interactively to create a desirable mesh. The quality of mesh has been shown to have a major impact on the overall performance of the numerical model.The model has been applied to simulate two-dimensional dam break flows for which transient water level distributions measured within a laboratory flume were available. In total 12 model runs were undertaken to test the model for various flow conditions. These conditions include: (1) different bed slopes (ranging from zero to 0.8%), (2) different upstream and downstream water level conditions, and (3) initially wet and dry bed conditions, downstream of the dam. Detailed comparisons have been made between model predicted and measured water levels and good agreement achieved between both sets of results. The model was then used to predict water level and velocity distributions in a real estuary, i.e. the Ribble Estuary, where the bed level varies rapidly at certain locations. In order to model the whole estuary, a 1-D numerical model has also been used to model the upper part of the estuary and this model was linked dynamically to the 2-D model. Findings from this application are given in detail.     

Inertia-gravity waves in the troposphere and lower stratosphere associated with a jet stream exit region

Radar measurements at Aberystwyth (52.4°N, 4.1°W) of winds at tropospheric and lower stratospheric heights are shown for 12–13 March 1994 in a region of highly curved flow, downstream of the jet maximum. The perturbations of horizontal velocity have comparable amplitudes in the troposphere and lower stratosphere with downward and upward phase propagation, respectively, in these two height regions. The sense of rotation with increasing height in hodographs of horizontal perturbation velocity derived for hourly intervals show downwards propagation of energy in the troposphere and upward propagation in the lower stratosphere with vertical wavelengths of 1.7 to 2.3 km. The results indicate inertia-gravity waves propagating in a direction similar to that of the jet stream but at smaller velocities. Some of the features observed contrast with those of previous observations of inertia-gravity waves propagating transverse to the jet stream. The interpretation of the hodographs to derive wave parameters has taken account of the vertical shear of the background wind transverse to the direction of wave propagation.     

Experiment and simulation of turbulent flow in local scour around a spur dyke

The turbulent flow in the local scour hole around a single non-submerged spur dyke is investigated with both experimental and numerical methods. The experiments are conducted under clear-water scour regime with an impermeable spur dyke. The scour geometry and flow velocities are measured in details with a high-resolution laser displacement meter, electro-magnetic velocimetries and PIV （Particle image velocimetry）. A 3D non-linear k-ε model is developed to simulate the complex local flow field around the scour area. The numerical model is formulated using FVM （Finite volume method） on a collocated unstructured mesh, capable of resolving complex geometries and boundaries. It is found that the simulation results are reasonably consistent with those of the experimental measurements. Based on the study results, the nature of the flow structure around a spur dyke with local scour hole is analyzed.     

Convection Regime of the Development of an Impulse–Buoyant Transient Jet Submerged in a Homogeneous Fluid

An intergral model of a transient vertical impulse–buoyant jet is suggested. The model contains a universal equation describing the propagation of the upper boundary of the convection front depending on the strength of the point source of buoyancy and momentum. The convective regime of jet propagation is considered, which includes a class of self-similar solutions corresponding to the buoyancy sources, whose strength varies with time following power and exponential laws. The obtained numerical solutions are compared with available experimental data on the profiles of vertical velocity and buoyancy on the jet axis.     

三角网格有限元法声波与弹性波模拟频散分析

本文对声波与弹性波方程进行有限元法离散,构造有限元法频散关系的一般特征值问题,分析了时间离散格式为中心差分的三角网格有限元法声波与弹性波模拟的频散特性. 比较了三种质量矩阵即分布式质量矩阵、集中质量矩阵和混合质量矩阵对有限元法频散的影响;选取四种典型三角网格,分析了混合质量矩阵有限元（MFEM）频散的方向各向异性;数值频散、方向各向异性随插值阶数的增加逐渐减弱,当空间为三阶插值时,频散主要表现为随采样率的变化而几乎无明显方向各向异性, 其频散幅值也较小. 控制其他影响因素不变的情况下,研究了不同波速比介质中弹性波的数值频散. 最后给出了三角网格MFEM的数值耗散性.     

Correlation of dynamic characteristics of a super‐tall building from full‐scale measurements and numerical analysis with various finite element models

The Di Wang Tower located in Shenzhen has 79 storeys and is about 325 m high. Field measurements have been conducted to investigate the dynamic characteristics of the super‐tall building. In parallel with the field measurements, seven finite element models have been established to model the multi‐outrigger‐braced tall building and to analyse the effects of various arrangements of outrigger belts and vertical bracings on the dynamic characteristics and responses of the Di Wang Tower under the design wind load and earthquake action. The distributions of shear forces in vertical structural components along the building height are also presented and discussed. The results from detailed modelling of group shear walls with several types of finite elements are addressed and compared to investigate various modelling assumptions. Finally, the performance of the finite element models is evaluated by correlating the natural frequencies and mode shapes from the numerical analysis with the finite element models and the field measurements. The results generated from this study are expected to be of interest to professionals and researchers involved with the design of tall buildings. Copyright © 2004 John Wiley & Sons, Ltd.