有限差分法精确计算地震走时

快速精确计算走时在层析成像反演中起重要作用，本文采用有限差分法来计算。首先通过变换，将程函方程化成一守恒形式；然后再差分离散，得到一阶和二阶精度的差分格式，并证明了差分格式条件稳定。模型试算表明：本文算法简明稳定，能快速高精度地计算地震走时，是一种值得注意和采用的走时计算方法。     

间断水波在梯级水库中的传播

基于矢通量分裂得到了一维浅水方程组的隐式守恒有限差分格式,并对有底坡、有摩阻的梯级水库同时溃坝和相继溃坝的溃坝波传播进行了数值模拟。对数值结果作了分析。结果合理,方法有效可行。     

计算地球流体力学的回顾、进展及展望

简要回顾了计算地球流体力学的发展历史,概括介绍了计算地球流体力学的研究进展及最新发展方向。针对线性发展方程的计算稳定性问题,介绍了判定线性发展方程初值问题的稳定性判别条件：CLF条件。对分析线性偏微分方程差分格式计算稳定性判据的Fourier方法、启发性稳定性分析方法,也做了简要的介绍。同时,重点介绍了判定线性发展方程初边值问题计算稳定性的GKS理论。在非线性发展方程的计算稳定性方面,重点介绍的主要内容包括：非线性发展方程的计算紊乱现象和计算不稳定的原因;克服非线性发展方程计算不稳定的方法;瞬时平方守恒型差分格式的构造;隐式和显式完全平方守恒格式的设计;强迫耗散非线性发展方程的计算稳定性问题。在计算地球流体力学的近期进展方面,重点介绍了非线性发展方程的计算稳定性与初值的关系,强迫耗散非线性发展方程显式准完全平方守恒格式的构造。对计算地球流体力学需要进一步研究的问题,也做了简要介绍。这些研究工作的介绍,无疑对推动计算地球流体力学的研究和大气海洋模式的研制具有一定的指导意义。     

扩展型Boussinesq水波方程的混合求解格式

为高效求解扩展型Boussinesq水波方程,建立了基于有限差分和有限体积方法的混合数值格式。将一维控制方程写为守恒形式,方程中通量部分采用有限体积方法求解,剩余部分采用有限差分方法求解。其中,有限体积方法采用Godunov类高分辨率格式,并结合HLL(Harten-Lax and van Leer)式黎曼问题近似解求界面数值通量,黎曼问题界面左右变量通过高精度状态插值方法(MUSCL)构筑。有限差分方法则采用具有二阶精度的中心差分公式进行。采用具有TVD(Total Variation Diminishing)性质的三阶龙格-库塔多步积分法进行时间积分。对数值模式进行了验证,数值结果同解析解或实验数据吻合良好。     

Crank-Nicolson差分格式及其稳定性研究

本文以自己独特的方式,构造了一维和二维抛物型方程的Ｃｒａｎｋ－Ｎｉｃｏｌｓｏｎ差分格式。本文不仅详细地给出了离散误差的表达式,而且论证了它们的稳定性。该差分格式具有精度高,稳定性好,计算量和存储量都比较小的特点,是一个很理想,便于应用的差分格式。     

解二维扩散方程的DuFort—Frankel差分格式

Ｄｕ Ｆｏｒｔ－Ｆｒａｎｋｅｌ差分格式是对Ｒｉｃｈａｒｄｓｏｎ格式进行修正得到的差分格式。本文将它从一维推广到二维,给出了二维Ｄｕ Ｆｏｒｔ－Ｆｒａｎｋｅｌ差分格式相容性所满足的条件,并严格论证了它的绝对稳定性。     

求解渗流问题的高精度差分格式

在求解渗流问题的传统差分格式中,只有Crank Nicolson格式具有对时间t的二阶精度。本文在导数超收敛点概念的基础上,提出一种求解渗流问题的三阶精度差分格式,并将其与显式差分格式叠加形成组合差分格式以改善格式的稳定、收敛条件。算例计算结果表明,该组合格式具有精度高,稳定收敛限制宽松,易于编程等优点      

二维Crank—Nicholson差分格式及其稳定性

利用积分插值法,把二维Ｃｒａｎｋ－Ｎｉｃｈｏｌｓｏｎ差分格式,由常系数推广到变系数情形。不仅导出了差分格式的截断误差,而且还应用能量估计法,详细论证了差分格式的绝对稳定性,该差分格式是一个精度高,稳定性好,便于应用的差分格式。     

马斯京根-康吉洪水演算方法的稳定性分析

基于运动波差分解的数值扩散,在一定条件下可以模拟扩散波物理扩散的概念,利用有限差分方法,建立了马斯京根-康吉洪水演算方法.对于该差分格式中的权重系数X,国内外长期以来没有统一的定论.本文据此利用Von Neumann稳定性分析方法,对该差分格式进行了较为深入的分析,且得到了该差分格式的稳定性条件.     

烃类垂向微渗漏过程的积木块模型及差分格式

这里提出以代表质量守恒的反应对流扩散方程作为主控方程的烃类垂向微渗漏方程组的差分格式,即双向一维分裂校正差分格式,并建立地层积木块模型对该格式的边界进行讨论。差分格式是预估～校正差分格式的一种改进形式,它融合了Crank-Nicolson格式、交替方向隐格式、预估～校正差分格式的特点,具有二阶差分精度,且无条件稳定。由于差分格式将每一步都归结为求解三对角线方程组,因此适合并行运算。数值实验表明,应用差分格式的数值模拟结果符合烃类垂向微渗漏过程的理论模型,可作为烃类垂向微渗漏过程分析的计算方法。     

Optimal difference schemes for Maxwell’s equations in solving forward problems of electromagnetic soundings

In this paper, the solution of two-dimensional Maxwell’s equations is considered using the Laguerre transform. Optimal parameters of the difference schemes for the equations are obtained and presented. Numerical values of these optimal parameters are given. Second-order difference schemes with the optimal parameters provide an accuracy of the solution of the equations that is comparable to the accuracy of the solution using fourth-order schemes. It is shown that, when using the Laguerre transform, the number of optimal parameters can be reduced compared to the Fourier transform. This reduction leads to a simplification of the difference scheme and a reduction in the amount of computation, i.e., to efficiency of the algorithm.     

Unified thermo-compositional-mechanical framework for reservoir simulation

We present a reservoir simulation framework for coupled thermal-compositional-mechanics processes. We use finite-volume methods to discretize the mass and energy conservation equations and finite-element methods for the mechanics problem. We use the first-order backward Euler for time. We solve the resulting set of nonlinear algebraic equations using fully implicit (FI) and sequential-implicit (SI) solution schemes. The FI approach is attractive for general-purpose simulation due to its unconditional stability. However, the FI method requires the development of a complex thermo-compositional-mechanics framework for the nonlinear problems of interest, and that includes the construction of the full Jacobian matrix for the coupled multi-physics discrete system of equations. On the other hand, SI-based solution schemes allow for relatively fast development because different simulation modules can be coupled more easily. The challenge with SI schemes is that the nonlinear convergence rate depends strongly on the coupling strength across the physical mechanisms and on the details of the sequential updating strategy across the different physics modules. The flexible automatic differentiation-based framework described here allows for detailed assessment of the robustness and computational efficiency of different coupling schemes for a wide range of multi-physics subsurface problems.     

Well-balanced central WENO schemes for the sediment transport model in shallow water

Sediment transport model in shallow water admits steady-state solutions in which the non-zero flux gradient is exactly balanced by the source term. In this paper, we develop high-order well-balanced central weighted essentially non-oscillatory schemes for the sediment transport model. In order to maintain the well-balanced property, we first reformulate the governing equations by an equivalent form and propose a novel reconstruction procedure for the Runge-Kutta flux. Rigorous theoretical analysis as well as extensive numerical examples all suggest that the present schemes preserve the well-balanced property. Moreover, the resulting schemes keep genuine high-order accuracy for general solutions.     

二维分数阶对流-弥散方程的数值解

对二维时间分数阶对流-弥散方程和二维空间分数阶对流-弥散方程分别建立了差分格式,实现了对其的数值求解。针对理想算例进行计算求解,分析了时间和空间分数阶阶数取不同值时的扩散变化规律,验证了各自所描述的时间相关性与空间相关性。同时与传统的二维整数阶对流-弥散方程的求解结果作了对比。当时间和空间分数阶阶数α与γ分别取整数时,二维时间分数阶对流-弥散方程和二维空间分数阶对流-弥散方程都与传统二维整数阶对流-弥散方程的计算结果相同,说明提出的对二维分数阶对流-弥散方程的数值求解方法是可行的。其结果对地下水溶质运移的进一步研究提供了有效的手段。     

On the application of the maximum entropy meshfree method for elastoplastic geotechnical analysis

In this study, the Maximum Entropy Meshfree (MEM) method is employed for analysing geotechnical problems involving material nonlinearity, assuming small strains. The efficiency of the MEM method is evaluated through several solution schemes for the global governing equations as well as the local constitutive equations. The conventional implicit approach involving the Newton-Raphson method and an explicit adaptive dynamic relaxation technique are employed for solving the governing equations, while local constitutive equations are solved numerically as well as analytically. Two- and three-dimensional numerical experiments are performed to study the efficiency of different configurations of the solution scheme, which leads to some important conclusions about application of the MEM method in geotechnical problems.     

Automatic substepping integration of the subloading tij model with stress path dependent hardening

This paper discusses the numerical integration of the subloading tij model. This is an elastoplastic model with stress path dependent hardening, which can predict the behaviour of normally consolidated clays or loose sands, as well as of over-consolidated clays or dense sands, with a small number of material parameters. Three features distinguish the subloading tij model from the conventional ones: (a) the use of a modified stress space given by tensor t_ij; (b) the split of the plastic strain increments in two components leading to a stress path dependent hardening; and (c) the use of two yield surfaces (subloading yield surface and normal yield surface). This last feature is based on the concept of sub-yielding stress states and adds an extra internal strain-like hardening variable, related to the relative density state, which demands its own evolution law. The three characteristics above greatly improve the prediction capabilities of the model, with respect to those of the well-known Cam clay model, at the cost of only two additional parameters. Nonetheless, the numerical integration of the constitutive equations of subloading tij model is a bit challenging, mainly due to the stress path dependent hardening. In order to integrate the equations of subloading tij model in the same way as for any conventional model, the authors reformulated its equations in a simpler and direct manner. Here, these equations are integrated using multi-step explicit schemes, such as modified-Euler and Runge–Kutta–Dormand–Price, with automatic error control. Simple forward-Euler scheme is also used for the sake of comparison. The results show that the modified-Euler scheme is more accurate as well as faster than the other schemes analysed over a wide range of error tolerance. Besides, the automatic feature of these schemes is a great convenience for the users of numerical codes.     

Finite element methods for geothermal reservoir simulation

Two finite element algorithms suitable for long term simulation of geothermal reservoirs are presented. Both methods use a diagonal mass matrix and a Newton iteration scheme. The first scheme solves the 2N unsymmetric algebraic equations resulting from the finite element discretization of the equations governing the flow of heat and mass in porous media by using a banded equation solver. The second method, suitable for problems in which the transmissibility terms are small compared to the accumulation terms, reduces the set of N equations for the Newton corrections to a symmetric system. Comparison with finite difference schemes indicates that the proposed algorithms are competitive with existing methods.     

A fully-coupled discontinuous Galerkin method for two-phase flow in porous media with discontinuous capillary pressure

In this paper, we formulate and test numerically a fully-coupled discontinuous Galerkin (DG) method for incompressible two-phase flow with discontinuous capillary pressure. The spatial discretization uses the symmetric interior penalty DG formulation with weighted averages and is based on a wetting-phase potential/capillary potential formulation of the two-phase flow system. After discretizing in time with diagonally implicit Runge-Kutta schemes, the resulting systems of nonlinear algebraic equations are solved with Newton’s method and the arising systems of linear equations are solved efficiently and in parallel with an algebraic multigrid method. The new scheme is investigated for various test problems from the literature and is also compared to a cell-centered finite volume scheme in terms of accuracy and time to solution. We find that the method is accurate, robust, and efficient. In particular, no postprocessing of the DG velocity field is necessary in contrast to results reported by several authors for decoupled schemes. Moreover, the solver scales well in parallel and three-dimensional problems with up to nearly 100 million degrees of freedom per time step have been computed on 1,000 processors.     

An efficient optimum interpolation scheme for objective analysis over Indian region

In the optimum interpolation scheme, the weights for the observations are computed by solving a set of linear equations for every grid point. As the number of observations increases particularly over data-rich regions, the matrix dimension increases and the computer time required to solve these equations to determine weights increases considerably. In order to reduce the computer time for computing the weights, Tanguay and Robert suggested schemes in which the gaussian function representing the autocorrelation function has been approximated by a second-order and also by a fourth-order Taylor series expansion. This resulted in the solution of matrices of order 4 or 9 respectively to obtain weighting functions irrespective of the number of observations used in the analysis. In the present study, the analyses of mean sea level pressure and geopotential height at 700 mbar level have been carried out for five days using the above two schemes and the regular OI scheme. The analyses are found to be similar in all the three cases suggesting that a lot of computer time could be saved without sacrificing the analysis accuracy by using the modified scheme in which the second-order approximation is utilized.     

辅助微分方程完全匹配层在声波方程数值模拟中的应用

在地震波数值模拟中,为提高算法精度,需要使用高阶时间更新格式,而普通的非分裂完全匹配层（PML）吸收边界局限于低阶时间格式。辅助微分方程完全匹配层（ADE PML）是一种可以适应任意阶时间格式的非分裂完全匹配层技术,且可以直接应用复频移拉伸算子以提高PML在高角度入射时的效果。作者将ADE PML应用于声波方程四阶Runge Kutta时间格式的数值模拟中,对其吸收效能进行了检验。数值模拟表明,复频移ADE PML在高角度入射时表现优于非复频移ADE PML。另外,不同辅助变量更新格式的吸收效果存在微小差异,显格式下计算结果与解析解吻合较好。长时间能量衰减计算表明ADE PML可以稳定至2 × 10⁵时间步。