两相介质近场波动模拟的一种解耦有限元方法

本文将求解近场波动问题的一种解耦技术推广到两相介质,得到了一种求解两相介质近场波动问题的直接解耦方法,包括集中质量有限元模型、时域显式积分格式和局部人工边界条件. 首先应用加权残数法,并依据波动模拟的精度要求,得到了两相介质集中质量有限元模型. 然后,结合两相介质中波动的衰减特性,实现了透射边界在两相介质近场波动中的运用. 最后,通过数值实验,并与解析解对比,验证了本文方法的有效性.     

层状饱水软土地基三维非轴对称动力响应分析方法

将饱水软土地基视为两相介质、考虑水的渗流和土骨架9变形的耦合作用，用Fourier展开和Hankel积分变换分析三维非称对称饱和弹性土层波动方程，用刚度矩阵方法，建立了层状饱和软土地基三维非轴对称动力响应的解析分析方法。以数值算例对比分析了单相土介质与两相饱和土介质三维非轴对称稳态动力响应，结果表明：在饱水软土地基动力响应分析中应该考虑土体中孔隙流体的影响。     

利用走时反演二维多孔介质渗透率

综合利用直接方法和优化方法研究二维多孔介质渗透率反演问题. 问题的提法与前人不同,是由区域内各点流体渗透的走时来反演渗透率. 该反问题的求解可分成两步进行,归结为两个相应的子问题进行研究. 首先由流体走时反演二维多孔介质速度场的数值算法,然后由二维多孔介质速度场反演二维多孔介质渗透率的数值算法,最后给出数值例子. 数值结果表明了所用数值方法的正确性和有效性.     

基于横向各向同性BISQ方程的弹性波传播数值模拟

Biot流动和喷射流动是含流体多孔隙介质中流体流动的两种重要力学机制. 近年来，利用同时处理这两种力学机制的BISQ(Biot-Squirt)模型，弹性波衰减和频散的问题已被广泛研究；然而基于BISQ方程的波场数值模拟尚未见到公开的报道．本文从BISQ方程出发，利用交错网格方法对横向各向同性孔隙介质中不同频率和相界情况，以及双层介质中的弹性波传播进行数值模拟，研究了在同时考虑两种流动机制作用情况下地震波和声波的传播特性及传播过程中出现的各种波动现象．      

淤积泥砂对水平地运动作用时刚性坝面动压力的影响

文章基于固体介质，流体介质，液固两相饱和介质的时域显示有限元波动分析方法，研究了可压缩库水条件下淤积泥砂对水平地运动作用时刚性坝面上动压力的影响，研究中分析了作为两相介质处理的淤积泥砂的饱和度、厚度、渗透系数和孔隙率对坝面动压力的影响，还比较了作为流程两相介质，单相固体介质、单相重流体介质处理的不同淤积泥砂层模型间的计算结果差异。     

基于时空守恒元和解元(CE/SE)方法的孔隙介质多相流动计算

将时空守恒元/解元(CE/SE)方法推广到二维孔隙介质多相流问题的数值计算中,采用人工压缩法耦合速度和压力,同时结合杂交粒子水平集方法捕捉物质界面.提出一套完整的二维欧拉型孔隙介质非稳态多相不可压缩黏性流动计算方案.通过对溃坝和液滴在重力作用下的运动和变形问题的数值模拟,验证了方法的精度和有效性.在此基础上,提出了一个新的孔隙介质两相流物理模型——双层流体顶盖驱动方腔流.     

三维层状孔隙介质中弹性波的一种积分表达式II：正确性检验和数值模拟实验

基于前一篇文章中得到的关于三维层状孔隙介质中弹性波场的积分形式半解析解,本文通过离散波数法开展了数值模拟.将全空间均匀孔隙介质中单力点源和爆炸点源作用下弹性波场的解析解和我们的数值模拟结果进行对比,发现两者是完全一致的.而在一个两层半空间模型下的数值模拟,验证了固相位移Green函数的9组空间互易性情况.通过以上两种对比检验,验证了半解析解理论公式、数值模拟方法以及相应程序代码的正确性和可靠性.随后利用敏感度分析研究了不同的介质参数变化对爆炸点源在界面上会产生的反射波场的影响.通过垂直地震剖面模型的数值模拟,发现弹性波场能很好地反映孔隙介质物理性质的变化,同时也讨论了动力协调这一孔隙介质中的特殊现象.我们发展的基于半解析解的数值模拟方法可以为三维层状孔隙介质中弹性波传播特征的研究提供一种可供选择的有效工具和手段.     

用二相介质有限元方法对共和7.0级地震前兆的数值模拟

利用固液二相介质平面应变问题的有限元程序对共和７．０级地震的前兆现象的时空分布进行了数值模拟．研究中不仅考虑了固相的非线性、硬化及膨胀等效应，也充分考虑了水对前兆的影响，因此模拟结果与实际的前兆分布特征符合较好．     

VTI介质qP波方程高精度有限差分算子

波动方程有限差分法是一种使用广泛的地震波数值模拟方法.但是有限差分法本身固有存在着数值频散问题，会降低地震波场模拟的精度与分辨率.为了克服常规有限差分算子的数值频散，本文针对VTI介质地震波数值模拟问题，构造了频率-空间域qP波波动方程高精度有限差分优化算子，根据最优化理论中高斯-牛顿法确定了高精度有限差分算子的优化系数.利用常规差分算子和高精度优化差分算子对归一化相速度的频散关系精度进行了对比分析，并对均匀各向同性介质和均匀VTI介质中的qP波地震波场进行了有限差分数值模拟，通过频散关系精度分析和波场数值模拟结果表明：有限差分优化算子具有较高的波场数值模拟精度，有效压制了传统有限差分算子数值模拟中的数值频散现象，提高了有限差分算子精度，为VTI介质频率-空间域qP波正演模拟奠定了基础.     

错格实数傅里叶变换微分算子及其在非均匀介质波动传播研究中的应用

提出一种新的数值微分运算方法，即错格实数傅里叶变换微 分法. 该方法的运算速度 比错格复数傅里叶变换数值微分解法快0.33倍；因为该微分算法在整个微分运算过程中保留 了奈奎斯特分量，使得它比普通分格的实数傅里叶变换数值微分算法的精度高，稳定性好. 将该方法和Cagniard De Hoop解析法在求解半无限空间地震波动的问题中进行比较，结果 表明，新微分法的精度和解析方法的精度相同. 在非均匀介质中的地震波传播数值模拟的结 果表明，该方法是一种研究非均匀介质中地震波传播问题的有效的数值微分方法.     

On different numerical inverse Laplace methods for solute transport problems

Numerical inversion is required when Laplace transform cannot be inverted analytically by manipulating tabled formulas of special cases. However, the numerical inverse Laplace transform is generally an ill-posed problem, and there is no universal method which works well for all problems. In this study, we selected seven commonly used numerical inverse Laplace transform methods to evaluate their performance for dealing with solute transport in the subsurface under uniform or radial flow condition. Such seven methods included the Stehfest, the de Hoog, the Honig–Hirdes, the Talbot, the Weeks, the Simon and the Zakian methods. We specifically investigated the optimal free parameters of each method, including the number of terms used in the summation and the numerical tolerance. This study revealed that some commonly recommended values of the free parameters in previous studies did not work very well, especially for the advection-dominated problems. Instead, we recommended new values of the free parameters for some methods after testing their robustness. For the radial dispersion, the de Hoog, the Talbot, and the Simon methods worked very well, regardless of the dispersion-dominated or advection-dominated situations. The Weeks method can be used to solve the dispersion-dominated problems, but not the advection-dominated problems. The Stehfest, the Honig–Hirdes, and the Zakian methods were recommended for the dispersion-dominated problems. The Zakian method was efficient, while the de Hoog method was time-consuming under radial flow condition. Under the uniform flow condition, all the methods could present somewhat similar results when the free parameters were given proper values for dispersion-dominated problems; while only the Simon method, the Weeks method, and the de Hoog method worked well for advection-dominated problems.     

Practical implementation of the fractional flow approach to multi-phase flow simulation

Fractional flow formulations of the multi-phase flow equations exhibit several attractive attributes for numerical simulations. The governing equations are a saturation equation having an advection diffusion form, for which characteristic methods are suited, and a global pressure equation whose form is elliptic. The fractional flow approach to the governing equations is compared with other approaches and the implication of equation form for numerical methods discussed. The fractional flow equations are solved with a modified method of characteristics for the saturation equation and a finite element method for the pressure equation. An iterative algorithm for determination of the general boundary conditions is implemented. Comparisons are made with a numerical method based on the two-pressure formulation of the governing equations. While the fractional flow approach is attractive for model problems, the performance of numerical methods based on these equations is relatively poor when the method is applied to general boundary conditions. We expect similar difficulties with the fractional flow approach for more general problems involving heterogenous material properties and multiple spatial dimensions.     

Modern Computational Technologies for Calculating Open Flow Dynamics

General approaches to solving the problems of open flow dynamics are discussed. Most attention is concentrated on the selection of numerical methods, construction of special computational grids, numerical experiments, and analysis of the adequacy of the obtained solution to the field observation results. The results of effective solution of some important practical problems are given to illustrate the potentialities of computational technologies.     

Numerical issues in plasticity models for granular materials

Friction plays a fundamental role in the mechanics of granular materials. Two problems are considered: (i) heap formation and (ii) granular flow. Both problems admit closely related mathematical models. In each case, analytical and numerical difficulties are discussed. Efficient and reliable numerical methods are proposed and implemented. The results are illustrated by several computational experiments.     

A spatially and temporally adaptive solution of Richards’ equation

Efficient, robust simulation of groundwater flow in the unsaturated zone remains computationally expensive, especially for problems characterized by sharp fronts in both space and time. Standard approaches that employ uniform spatial and temporal discretizations for the numerical solution of these problems lead to inefficient and expensive simulations. In this work, we solve Richards’ equation using adaptive methods in both space and time. Spatial adaption is based upon a coarse grid solve and a gradient error indicator using a fixed-order approximation. Temporal adaption is accomplished using variable order, variable step size approximations based upon the backward difference formulas up to fifth order. Since the advantages of similar adaptive methods in time are now established, we evaluate our method by comparison with a uniform spatial discretization that is adaptive in time for four different one-dimensional test problems. The numerical results demonstrate that the proposed method provides a robust and efficient alternative to standard approaches for simulating variably saturated flow in one spatial dimension.     

Saltwater intrusion modeling in heterogeneous confined aquifers using two-relaxation-time lattice Boltzmann method

This study develops a lattice Boltzmann method (LBM) with a two-relaxation-time collision operator (LTRT) to solve saltwater intrusion problems. A directional-speed-of-sound (DSS) technique is introduced to take into account the hydraulic conductivity heterogeneity and discontinuity, as well as the velocity-dependent dispersion coefficient. The forcing terms in the LTRT model are customized in order to recover the density-dependent groundwater flow and mass transport equations. Using the LTRT with the squared DSS achieves at least second-order accuracy. The LTRT results are verified with Henry’s analytical solution as well as compared with several numerical examples and modified Henry problems that consider heterogeneous hydraulic conductivity and velocity-dependent dispersion. The numerical results show good agreement with the Henry analytical solution and with the numerical solutions obtained by other numerical methods.     

关于解地球物理中病态方程的若干问题

讨论了三种重要的数值方法：１．吉洪诺夫正则化方法和适应正则化方法；２．预条件子方法；３．数值相关性技术．首先对这几种方法的新进展进行了评述，其次提出了几种合理选择正则子的方法并且建立了一个适应正则化算法．根据地球物理中病态方程的特性，推广了预条件子的概念并且提出了选择预条件子的若干方法．同时讨论了这些方法在地球物理中的应用．     

A mass conservative scheme for simulating shallow flows over variable topographies using unstructured grid

Most available numerical methods face problems, in the presence of variable topographies, due to the imbalance between the source and flux terms. Treatments for this problem generally work well for structured grids, but most of them are not directly applicable for unstructured grids. On the other hand, despite of their good performance for discontinuous flows, most available numerical schemes (such as HLL flux and ENO schemes) induce a high level of numerical diffusion in simulating recirculating flows. A numerical method for simulating shallow recirculating flows over a variable topography on unstructured grids is presented. This mass conservative approach can simulate different flow conditions including recirculating, transcritical and discontinuous flows over variable topographies without upwinding of source terms and with a low level of numerical diffusion. Different numerical tests cases are presented to show the performance of the scheme for some challenging problems.     

Optimal weighting in the finite difference solution of the convection-dispersion equation

Oscillation and numerical dispersion limit the reliability of numerical solutions of the convection-dispersion equation when finite difference methods are used. To eliminate oscillation and reduce the numerical dispersion, an optimal upstream weighting with finite differences is proposed. The optimal values of upstream weighting coefficients numerically obtained are a function of the mesh Peclet number used. The accuracy of the proposed numerical method is tested against two classical problems for which analytical solutions exist. The comparison of the numerical results obtained with different numerical schemes and those obtained by the analytical solutions demonstrates the possibility of a real gain in precision using the proposed optimal weighting method. This gain in precision is verified by interpreting a tracer experiment performed in a laboratory column.