非均匀介质弹性波动方程的不规则网格有限差分方法

从弹性波动方程出发，提出了一种新的空间不规则网格有限差分方法，并用于求解非均匀各向异性介质中的弹性波正演问题。这种方法简单易行，对于复杂几何结构，例如低速层、套管井和非平面界面等，在较细的不规则网格上进行离散，计算时间和占用内存更少。与多重网格差分方法相比，该方法不需要粗、细网格之间的插值，所有网格差分计算在同一次空间迭代中完成。具有复杂几何交界面的模型计算，包括地下透镜体、套管井眼等，在确定弹性常数和密度后，用不规则网格的差分方法更易实现。该方法使用了Higdon吸收边界条件解决人工边界反射问题，引入了新的稳定性条件和网格频散条件，很好地消除了非物理散射波。理论模型的效值计算表明，该方法具有良好的稳定性和计算精度，在模拟非均匀介质弹性波传播时，比相同精度的规则网格有限差分方法计算速度更快。该方法易于推广到非结构网格和三维问题中。     

透射边界条件在波动谱元模拟中的实现:一维波动

多次透射公式(multi-transmitting formula,MTF)是一种具有普适性的局部人工边界条件,但其在近场波动数值模拟中一般与有限元法结合.由于波动谱元模拟的数值格式与有限元格式有极大的不同,传统的MTF在谱元离散格式中无法直接实现.为了使物理概念清楚、精度可控的多次透射人工边界条件能够适应波动谱元模拟的需求,首先指出多次透射边界与谱元离散格式结合的基本问题,并分析了空间内插和时间内插两种方案的可行性.然后从空间内插角度出发,提出基于拉格朗日多项式插值模式的MTF谱元格式,并采用一种简单内插方法实现高阶MTF.最后通过一维波动数值试验检验这些MTF谱元格式的精度,并讨论其数值稳定性.结果表明:对于一、二阶MTF,几种格式的精度相当;对于三、四阶MTF,基于谱单元位移模式插值的格式精度最高.相反,随着插值多项式阶次的升高,不同MTF格式的稳定临界值逐步降低,但是所有格式均在人工波速大大超过物理波速时才可能发生失稳.     

一种基于增量径向基函数插值的流场重构方法

由于流场参数重构中, 用于重构的基网格单元的物理参数波动量相对于均值较小, 径向基函数（RBF） 直接插值方法重构会产生较大的数值振荡, 论文提出了一种增量RBF 插值方法, 并用于有限体积的流场重构步, 明显改善了插值格式的收敛性和稳定性. 算例首先通过简单的一维模型说明该方法的有效性, 当目标函数波动量相对于均值为小量时, 增量RBF 插值能够抑制数值振荡; 进一步通过二维亚音速、跨音速定常无黏算例、静止圆柱绕流非定常算例以及超音速前台阶算例来说明该方法在典型流场数值求解中的通用性和有效性. 研究表明增量RBF 重构方法可陡峭地捕捉激波间断, 可有效改善流场求解的收敛性和稳定性, 数值耗散小, 计算效率高.      

水下声波传播过程的时域间断Galerkin有限元方法

本文构建了声压波动方程的改进时域间断Galerkin有限元方法.传统时域连续有限元方法在计算高梯度、强间断特征水中声波传播问题时往往会出现虚假数值振荡现象,这些数值振荡会影响正常波动的计算精度.为了解决这一问题,本文通过引入人工阻尼的方式构建了改进的时域间断Galerkin有限元方法,并针对具有高梯度、强间断特征的多障碍物复杂边界和层合液体介质声传播问题进行了计算.计算结果表明,与传统时域连续方法如N ew mark方法计算结果对比,所发展方法能较好地消除高梯度和强间断声压力波传播过程中虚假的数值振荡,具有较高的计算精度.问题的求解为进一步流固声耦合问题的研究奠定了基础.     

Seismic propagation simulation in complex media with non-rectangular irregular-grid finite-difference

This paper presents a finite-difference (FD) method with spatially non-rectangular irregular grids to simulate the elastic wave propagation. Staggered irregular grid finite difference operators with a second-order time and spatial accuracy are used to approximate the velocity-stress elastic wave equations. This method is very simple and the cost of computing time is not much. Complicated geometries like curved thin layers, cased borehole and nonplanar interfaces may be treated with nonrectangular irregular grids in a more flexible way. Unlike the multi-grid scheme, this method requires no interpolation between the fine and coarse grids and all grids are computed at the same spatial iteration. Compared with the rectangular irregular grid FD, the spurious diffractions from “staircase” interfaces can easily be eliminated without using finer grids. Dispersion and stability conditions of the proposed method can be established in a similar form as for the rectangular irregular grid scheme. The Higdon‘s absorbing boundary condition is adopted to eliminate boundary reflections. Numerical simulations show that this method has satisfactory stability and accuracy in simulating wave propagation near rough solid-fluid interfaces. The computation costs are less than those using a regular grid and rectangular grid FD method.     

A fully conservative high‐order upwind multi‐moment method using moments in both upwind and downwind cells

We propose a fully conservative high‐order upwind multi‐moment method for the conservation equation. The proposed method is based on a third‐order polynomial interpolation function and semi‐Lagrangian formulation and is a variant of the constrained interpolation profile conservative semi‐Lagrangian scheme with third‐order polynomial function method. The third‐order interpolation function is constructed based on three constraints in the upwind cell (two boundary values and a cell average) and a constraint in the downwind cell (a cell center value). The proposed method shows fourth‐order accuracy in a benchmark problem (sine wave propagation). We also propose a less oscillatory formulation of the proposed method. The less oscillatory formulation can minimize numerical oscillations. These methods were validated through scalar transport problems, and compressible flow problems (shock tube and 2D explosion problems). Copyright © 2016 John Wiley & Sons, Ltd.     

The Analysis of Dynamic Stress Intensity Factor for Semi-Circular Surface Crack Using Time-Domain BEM Formulation

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｔｒａｎｓｉｅｎｔａｎａｌｙｓｉｓｆｏｒ 3＿Ｄｅｌａｓｔｏｄｙｎａｍｉｃｃｒａｃｋｐｒｏｂｌｅｍｓｉｓａｎｉｍｐｏｒｔａｎｔｃｏｎｔｅｎｔｉｎｆｒａｃｔｕｒｅｍｅｃｈａｎｉｃｓ.Ｂｅｃａｕｓｅｏｆｔｈｅｄｉｆｆｉｃｕｌｔｉｅｓｏｆｔｈｅｐｒｏｂｌｅｍｓ,ｔｈｅｃａｌｃｕｌａｔｉｏｎｏｆＤＳＩＦｆｏｒ 3＿Ｄｃｒａｃｋｉｓａｐｒｏｂｌｅｍｕｎｓｏｌｖｅｄｃｏｍｐｌｅｔｅｌｙｕｐｔｏｔｈｅｐｒｅｓｅｎｔ.Ｃｈｅｎ[1]ｆｉｒｓｔｌｙｕｓｅｄｔｈｅｆｉｎｉｔｅｄｉｆｆｅｒｅｎｃｅ…     

An immersed smoothed point interpolation method (IS‐PIM) for fluid‐structure interaction problems

An immersed smoothed point interpolation method using 3‐node triangular background cells is proposed to solve 2D fluid‐structure interaction problems for solids with large deformation/displacement placed in incompressible viscous fluid. In the framework of immersed‐type method, the governing equations can be decomposed into 3 parts on the basis of the fictitious fluid assumption. The incompressible Navier‐Stokes equations are solved using the semi‐implicit characteristic‐based split scheme, and solids are simulated using the newly developed edge‐based smoothed point interpolation method. The fictitious fluid domain can be used to calculate the coupling force. The numerical results show that immersed smoothed point interpolation method can avoid remeshing for moving solid based on immersed operation and simulate the contact phenomenon without an additional treatment between the solid and the fluid boundary. The influence from information transfer between solid domain and fluid domain on fluid‐structure interaction problems has been investigated. The numerical results show that the proposed interpolation schemes will generally improve the accuracy for simulating both fluid flows and solid structures.     

A Control Volume Radial Basis Function Techniques for the Numerical Simulation of Saturated Flows in Semi-confined Aquifer

In this work, a Control Volume Radial Basis Function technique (CV-RBF) is adapted to solve ground water flow in the saturated zone of the semi-confined aquifer. The CV-RBF method differs from classical CV methods in the way that the flux at the cell surfaces is computed. A local RBF interpolation of the field variable is performed at the centres of the cell being integrated and its neighbours. This interpolation is then used to reconstruct the solution and its gradient in the integration points which support the flux computation. In addition, it is required that such interpolation satisfies the governing equation in a certain number of points placed around the cell centres. In this way, the local interpolations become equivalent to local boundary-value problems. The CV-RBF method is combined with a local remeshing technique in order to track the phreatic surface, where the gradients required to satisfy the kinematic condition are computed by the same local RBF interpolations used for the flux computation. The proposed numerical approach is validated in a series of three-dimensional groundwater flow problems where the operations of recharging and extracting water from a semi-confined aquifer are modelled.     

Numerical modeling of wave interactions with coastal structures by a constrained interpolation profile/immersed boundary method

A high‐order difference method based multiphase model is proposed to simulate nonlinear interactions between water wave and submerged coastal structures. The model is based on the Navier–Stokes equations using a constrained interpolation profile (CIP) method for the flow solver, and employs an immersed boundary method (IBM) for the treatment of wave–structure interactions. A more accurate interface capturing scheme, the volume of fluid/weighed line interface calculation (VOF/WLIC) scheme, is adopted as the interface capturing method. A series of computations are performed to verify the application of the model for simulations of fluid interaction with various structures. These problems include flow over a fixed cylinder, water entry of a circular cylinder and solitary waves passing various submerged coastal structures. Computations are compared with the available analytical, experimental and other numerical results and good agreement is obtained. The results of this study demonstrate the accuracy and applications of the proposed model to simulate the nonlinear flow phenomena and capture the complex free surface flow. Copyright © 2015 John Wiley & Sons, Ltd.     

透射边界条件在波动谱元模拟中的实现:二维波动

将邢浩洁和李鸿晶提出的多次透射公式(multi-transmitting formula,MTF)的谱元格式应用于均匀介质中线弹性SH波动问题的谱元模拟.假定紧邻人工边界的一层谱单元为具有直线边界的四边形单元,以保证每个人工边界节点都唯一对应一条指向内域的离散网格线.人工边界节点在某时刻的位移由该离散网格线上的节点在前若干时刻的位移确定,按照MTF谱元格式进行计算.通过平面波以一定角度传播的外源问题算例和点源脉冲自由扩散的内源问题算例,验证了方法的可行性以及对实际复杂波动问题的适用性.通过不同类型初值问题算例,在时域内分析了插值多项式阶次、人工波速和透射阶次三个参数对反射误差的影响.结果表明:插值多项式阶次较高的格式会表现出更好的精度,但总体上对反射误差的影响较小;人工波速对反射误差具有显著影响,当人工波速小于介质物理波速时反射误差较大,而当人工波速等于或稍大于介质物理波速时反射误差处于较低水平;透射阶次对反射误差具有决定性影响,表现在不失稳的情形下提高透射阶次能够迅速降低反射误差,但内源问题从三阶MTF开始出现飘移失稳,外源问题从二阶MTF开始出现轻微的飘移失稳.     

A RECONSTRUCTION METHOD FOR FINITE VOLUME FLOW FIELD SOLVING BASED ON INCREMENTAL RADIAL BASIS FUNCTIONS

由于流场参数重构中, 用于重构的基网格单元的物理参数波动量相对于均值较小, 径向基函数（RBF） 直接插值方法重构会产生较大的数值振荡, 论文提出了一种增量RBF 插值方法, 并用于有限体积的流场重构步, 明显改善了插值格式的收敛性和稳定性. 算例首先通过简单的一维模型说明该方法的有效性, 当目标函数波动量相对于均值为小量时, 增量RBF 插值能够抑制数值振荡; 进一步通过二维亚音速、跨音速定常无黏算例、静止圆柱绕流非定常算例以及超音速前台阶算例来说明该方法在典型流场数值求解中的通用性和有效性. 研究表明增量RBF 重构方法可陡峭地捕捉激波间断, 可有效改善流场求解的收敛性和稳定性, 数值耗散小, 计算效率高.     

Improvement of moving particle semi‐implicit method for simulation of progressive water waves

Precise simulation of the propagation of surface water waves, especially when involving breaking wave, takes a significant place in computational fluid dynamics. Because of the strong nonlinear properties, the treatment of large surface deformation of free surface flow has always been a challenging work in the development of numerical models. In this paper, the moving particle semi‐implicit (MPS) method, an entirely Lagrangian method, is modified to simulate wave motion in a 2‐D numerical wave flume preferably. In terms of consecutive pressure distribution, a new and simple free surface detection criterion is proposed to enhance the free surface recognition in the MPS method. In addition, a revised gradient model is deduced to diminish the effect of nonuniform particle distribution and then to reduce the numerical wave attenuation occurring in the original MPS model. The applicability and stability of the improved MPS method are firstly demonstrated by the calculation of hydrostatic problem. It is revealed that these modifications are effective to suppress the pressure oscillation, weaken the local particle clustering, and boost the stability of numerical algorithm. It is then applied to investigate the propagation of progressive waves on a flat bed and the wave breaking on a mild slope. Comparisons with the analytical solutions and experimental results indicate that the improved MPS model can give better results about the profiles and heights of surface waves in contrast with the previous MPS models. Copyright © 2017 John Wiley & Sons, Ltd.     

A dissipation‐free numerical method to solve one‐dimensional hyperbolic flow equations

In this paper, a numerical method to capture the shock wave propagation in 1‐dimensional fluid flow problems with 0 numerical dissipation is presented. Instead of using a traditional discrete grid, the new numerical method is built on a range‐discrete grid, which is obtained by a direct subdivision of values around the shock area. The range discrete grid consists of 2 types: continuous points and shock points. Numerical solution is achieved by tracking characteristics and shocks for the movements of continuous and shock points, respectively. Shocks can be generated or eliminated when triggering entropy conditions in a marking step. The method is conservative and total variation diminishing. We apply this new method to several examples, including solving Burgers equation for aerodynamics, Buckley‐Leverett equation for fractional flow in porous media, and the classical traffic flow. The solutions were verified against analytical solutions under simple conditions. Comparisons with several other traditional methods showed that the new method achieves a higher accuracy in capturing the shock while using much less grid number. The new method can serve as a fast tool to assess the shock wave propagation in various flow problems with good accuracy.     

Simulation seismic wave propagation in topographic structures using asymmetric staggered grids

A new 3 D finite- difference ( FD ) method of spatially asymmetric staggered grids was presented to simulate elastic wave propagation in topographic structures. The method approximated the first-order elastic wave equations by irregular grids finite difference operator with second-order time precise and fourth-order spatial precise. Additional introduced finite difference formula solved the asymmetric problem arisen in non-uniform staggered grid scheme, The method had no interpolation between the fine and coarse grids. All grids were computed at the same spatial iteration. Complicated geometrical structures like rough submarine interface, fault and nonplanar interfaces were treated with fine irregular grids. Theoretical analysis and numerical simulations show that this method saves considerable memory and computing time, at the same time, has satisfactory stability and accuracy.     

7阶WENO-S格式的计算效率研究

WENO-S格式是一类适合于含间断问题数值模拟的加权本质无振荡格式.这类格式的光滑因子满足对单频波为常数,这使得其近似色散关系与线性基底格式一致,并且具有良好的小尺度波动模拟能力．计算效率是数值方法性能指标的一个重要方面.由于WENO-S格式的光滑因子在各子模板上的计算公式除下标不同外形式一致,在计算线性对流方程相邻数值通量时,部分光滑因子完全相同．为此提出一种消除WENO-S格式冗余光滑因子计算的方法.该方法要求一条网格线上用于重构或插值的量可以表示为一个序列.基于此要求分析其对于几种不同物理问题的可行性和使用方法.以7阶WENO-S格式为例介绍了格式性质和去除冗余光滑因子计算的方法.该方法中预先计算和存储一条网格线上的所有光滑因子,在网格点较多的情况下,光滑因子计算次数约为原7阶WENO-S格式的1/4.对一维对流问题、球面波传播问题、二维旋转问题、二维小扰动传播问题及一维和二维无黏流动问题进行了数值模拟.结果表明该格式对多种流动结构具有良好的捕捉能力,并且同时具有良好的计算效率,去除冗余计算后又降低了约20%的计算时间.     

On the efficient application of weighted essentially nonoscillatory scheme

In this paper, the efficient application of high‐order weighted essentially nonoscillatory (WENO) reconstruction to the subsonic and transonic engineering problems is studied. On the basis of the physical considerations, two techniques are proposed to enhance the accuracy and efficiency of the WENO reconstruction. First, it is observed that the WENO scheme using characteristic variable has better accuracy and convergence speed than the scheme using primitive variable. For engineering problems with shock of moderate amplitude, on the basis of the Rankine–Hugoniot conditions, a simplified characteristic‐variable‐based WENO is developed. The simplified version significantly reduces the cost overhead without sacrificing the shock‐capturing capability. Second, in this work, it is found for viscous case that it is better to include the viscous effect. On the basis of a simple analysis, the viscous correction to the parameter ε in the WENO reconstruction is proposed. Numerical results indicate, with the proposed simplified characteristic‐variable‐based reconstruction and the viscous correction, that the nonlinear WENO interpolation is sharply activated in the region of shock jump, whereas in the shockless area, the WENO interpolation weights are tuned towards the designed optimal value for better accuracy. Compared with the original characteristic‐variable‐based WENO, the current implementation has similar accuracy and reduced cost. At the same time, compared with the primitive variable‐based WENO, better accuracy and convergence speed are obtained at marginal cost overhead. Several practical cases are calculated to demonstrate the accuracy and efficiency of the current methodology. Copyright © 2012 John Wiley & Sons, Ltd.     

Advances in viscous vortex methods—meshless spatial adaption based on radial basis function interpolation

Vortex methods have a history as old as finite differences. They have since faced difficulties stemming from the numerical complexity of the Biot–Savart law, the inconvenience of adding viscous effects in a Lagrangian formulation, and the loss of accuracy due to Lagrangian distortion of the computational elements. The first two issues have been successfully addressed, respectively, by the application of the fast multipole method, and by a variety of viscous schemes which will be briefly reviewed in this article. The standard method to deal with the third problem is the use of remeshing schemes consisting of tensor product interpolation with high‐order kernels. In this work, a numerical study of the errors due to remeshing has been performed, as well as of the errors implied in the discretization itself using vortex blobs. In addition, an alternative method of controlling Lagrangian distortion is proposed, based on ideas of radial basis function (RBF) interpolation (briefly reviewed here). This alternative is formulated grid‐free, and is shown to be more accurate than standard remeshing. In addition to high‐accuracy, RBF interpolation allows core size control, either for correcting the core spreading viscous scheme or for providing a variable resolution in the physical domain. This formulation will allow in theory the application of error estimates to produce a truly adaptive spatial refinement technique. Proof‐of‐concept is provided by calculations of the relaxation of a perturbed monopole to a tripole attractor. Copyright © 2005 John Wiley & Sons, Ltd.     

A novel finite point method for flow simulation

A novel finite point method is developed to simulate flow problems. The mashes in the traditional numerical methods are supplanted by the distribution of points in the calculation domain. A local interpolation based on the properties of Taylor series expansion is used to construct an approximation for unknown functions and their derivatives. An upwind‐dominated scheme is proposed to efficiently handle the non‐linear convection. Comparison with the finite difference solutions for the two‐dimensional driven cavity flow and the experimental results for flow around a cylinder shows that the present method is capable of satisfactorily predicting the flow separation characteristic. The present algorithm is simple and flexible for complex geometric boundary. The influence of the point distribution on computation time and accuracy of results is included. Copyright © 2002 John Wiley & Sons, Ltd.     

A simple high‐resolution advection scheme

A simple, robust, mass‐conserving numerical scheme for solving the linear advection equation is described. The scheme can estimate peak solution values accurately even in regions where spatial gradients are high. Such situations present a severe challenge to classical numerical algorithms. Attention is restricted to the case of pure advection in one and two dimensions since this is where past numerical problems have arisen. The authors' scheme is of the Godunov type and is second‐order in space and time. The required cell interface fluxes are obtained by MUSCL interpolation and the exact solution of a degenerate Riemann problem. Second‐order accuracy in time is achieved via a Runge–Kutta predictor–corrector sequence. The scheme is explicit and expressed in finite volume form for ease of implementation on a boundary‐conforming grid. Benchmark test problems in one and two dimensions are used to illustrate the high‐spatial accuracy of the method and its applicability to non‐uniform grids. Copyright © 2007 John Wiley & Sons, Ltd.