改进时间分数阶模型模拟非Fick溶质运移

通过将经典时间分数阶对流-弥散方程的等待时间分布函数的尾部修改为指数型,推导出了改进时间分数阶对流-弥散方程,并提出有效的时空算子分裂数值求解方法。对两个理想算例和一个实际算例进行计算,结果表明,改进的时间分数阶对流-弥散方程继承了时间分数阶对流-弥散方程能模拟穿透曲线幂率型拖尾分布的优点,还可模拟穿透曲线尾部由幂率型转换到指数型的过程;特征时间λ、分数阶指数γ和两相容量比例系数β共同决定了运移行为。改进的新模型可以区分非均质介质中流动相和非流动相中的溶质浓度, 更细微地模拟非Fick溶质运移行为。     

基于常Q模型的分数阶粘弹介质数值模拟方法

基于常Q模型的解耦分数阶拉普拉斯算子粘滞波动方程,可以分开模拟振幅衰减和相位错动。但该方程拉普拉斯算子的阶数是随空间变化的,因此数值求解存在一定困难。这里基于截断的泰勒展开,经过一系列近似,推导出拉普拉斯算子的阶数与空间无关的解耦分数阶粘滞弹性波动方程。采用中心差分计算时间导数,使用交错网格伪谱法计算空间导数。数值算例表明,新的方程在处理非均匀介质时具有精度高,计算简便的优点。     

二维输运方程高精度数值模拟

采用剖开算子法,把二维输运问题剖分为两个子初值问题(对流分步、扩散分步)。在任意三角形网格中,分别对不同性质的算子采用各自适合的算法,即采用特征线法求解对流分步,采用半隐式有限元法求解扩散分步。重点探讨了对流插值问题,给出了一种完全对称三次插值模式,有效地减少了数值阻尼。为了克服高阶插值数值震荡问题,计算中保证了函数及其一阶偏导数连续。算例表明,数值方法模拟结果与精确解吻合较好。该算法在求解输运方程(包括纯对流输运方程)时,既能有效减少数值阻尼,也能保证计算中不出现数值震荡。     

分数阶时间导数计算方法在含黏滞流体黏弹双相VTI介质波场模拟中的应用

相对于整数阶导数,分数阶微分算子可以更简洁地描述具有历史依赖性和空间全域相关性的复杂力学和物理过程。但是对分数阶波动方程进行数值模拟,计算量和存储量均较大,尤其对长时间或大计算域的模拟更是如此。文中给出了3种计算方法：全局记忆法、短时记忆法、自适应记忆法,并将这3种方法应用于含黏滞流体黏弹双相VTI （横向各向同性）介质分数阶波传播方程正演。通过对比3种方法的模拟精度、计算时间及占用内存发现：虽然短时记忆法可以通过设置短时记忆长度来调整计算时间与所占内存,但是短时记忆长度越短,精度越差;而自适应记忆法在保证精度的前提下,是短时记忆法与全局记忆法在计算时间与占用内存两方面的折衷。最后对各方法的利弊进行总结,为后续正演模拟及新的分数阶数值算法开发提供方法上的参考。在正演过程中,不仅要使所建模型更贴近实际地下介质,还需对选取的数值算法在计算时间、计算存储量和精度之间进行利弊权衡,以得到一个比较合理的数值算法。     

对流弥散方程不同有限元处理方法的比较分析

Galerkin有限元在处理含第二类边界条件的对流弥散方程时,针对对流项和弥散项有两种不同的格林积分变换,所得数值结果的精度也不同。一种方法是把对流和弥散项整体考虑实施格林积分转换(降低微分阶数,由二阶降成一阶),应用边界条件,得出变分方程;另一种处理方法是只对弥散项实施积分变换,应用边界条件,得出变分方程。以一维问题为参考,对两种方法的数值结果与解析解进行比较分析。     

求解对流占优地下水水质模型的有限元高精度广义迎风单元均衡格式及应用

根据水质模型的具体特点，对不同的方程采用不同方法，水流问题用有限元法；对流弥散方程先用算子分裂的方法分解为两个方程，即对流方程和弥散方程，前者用高精度广义迎风格式求解，对弥散方程则采用多单元均衡格式法求解，最后合成为高精度广义迎风均衡格式求出溶质浓度。通过对数值实验例子的计算和实验溶质迁移的模拟，可以看出在求解对流弥散定解问题时，广义迎风均衡格式克服了有限元数值波动和浓度出现负值的问题，与有限元相比有较大改进。     

基于分数阶Kelvin模型的饱和黏土一维流变固结分析

引入基于Caputo分数导数的弹壶元件修正Kelvin模型,以描述饱和黏土的一维流变本构关系。沿用Terzaghi饱和土一维固结理论的假设推导流变固结方程,引入Laplace变换和基于Fourier级数展开的Laplace数值逆变换解法进行数值求解。通过与基于整数阶导数模型解析解的对比,证明数值解法的有效性。通过对文献中一维流变固结试验结果的模拟,验证修正Kelvin模型的适用性。然后分析弹壶元件中分数导数阶数和黏滞系数对地基流变固结进程的影响。计算结果表明,在固结开始相当长的一段时间内,孔隙水压的整体消散速度要快于Terzaghi一维固结理论,但在固结后期则会慢于后者;而且在整个固结过程中,地基沉降速率都要慢于后者。总体来看,地基沉降滞后于孔压消散,并且分数导数阶数越小或黏滞系数越大,这种现象就越明显,而且沉降稳定需要的时间越长。     

非均匀孔隙介质中的弥散作用——(2)层状二维空间间断系统的预报

第一部分(Plumb和whitaker,本期)提出的闭合法对层状系统和二维空间间断系统进行了求解,其中我们所说的w和η区形成了在空间分布的孔隙介质。闭合法问题的解使我们得到确定大范围弥散张量所需要的数据和解大范围弥散方程所需要的另外一些参数。层状系统的理论结果展示的纵向弥散系数与利用室内视均质孔隙介质根据传统关系式所得到的纵向弥散系数相比要增大几个数量级。对于这种情况的研究,混合时间—空间系数看来是重要的,而且能够引起弥散系数增大的条件该系数值也增大。二维空间间断系统的结果与实验观测是一致的,而且我们发现非均匀性的规模和两个区之间的水力传导系数的差别对大范围弥散张量具有重要影响。正如层状介质情况那样,时间—空间混合导数(不对称向量)的系数的预测值表明该项可能是非均质孔隙介质的物质运移预报的一个重要的因子。     

2.5维各向异性介质中多波波场正演模拟与波场分析

完全3维弹性波数值模拟计算时间长,并且占用庞大的计算资源,这不利于在计算机配置不高的情况下进行科学研究,而二维弹性波数值模拟又达不到三维模拟的精度;同时,当模型、波场空间分布比较复杂时,传统的3维波动方程拟谱法模拟结果比较差.因此,在较高数值精度的一阶应力-速度弹性波动方程的基础上,采用傅氏变换仅计算y方向的偏导数,利用有限差分方法计算x、z方向和时间的偏导数,即利用2.5维数值模拟方法,实现在二维介质中计算三维弹性波场.最后通过数值模拟实现了在各向异性介质中多波波场的数值模拟,验证了2.5维方法是一种高精度、高效率、且能适应复杂模型的正演模拟方法,通过波场分析进一步认识了波在各向异性介质中的传播规律.     

地下水溶质运移数值模拟中减少误差的新方法

地下水中污染物运移的数值模拟方法一直是学界的研究热点问题.而如何减少与消除对流-弥散方程数值解中浓度陡锋面附近的数值振荡与数值弥散,更是研究的前沿与难点.提出了一种地下水溶质运移数值模拟中减少数值弥散的新方法.该方法的核心思想是在水动力弥散系数上加上一个数值弥散估算值,得到一个修正弥散系数,用其替代方程中有明确物理意义的水动力弥散系数进行计算.并提出了一个参数——数值弥散因子(μ_NDF),可以根据研究需要进行参数分区并适当调节该因子的大小,从而达到控制数值振荡,减小数值弥散的目的.从一维到二维的多个数值算例的模拟计算结果表明,该方法能在消除数值振荡的基础上,较好地减少数值弥散,达到满意的精度.     

Using a shifted Grünwald-Letnikov scheme for the Caputo derivative to study anomalous solute transport in porous medium

We propose an extension of the shifted Grünwald-Letnikov method to solve fractional partial differential equations in the Caputo sense with arbitrary fractional order derivative α and with an advective term. The method uses the relation between Caputo and Riemann-Liouville definitions, the shifted Grünwald-Letnikov, and the traditional backward and forward finite difference method. The stability of the method is investigated for the implicit and explicit scheme with homogeneous boundary conditions, and a stability criterion is found for the advective-dispersive equation. An application of the method is used to solve contaminant diffusion and advective-dispersive problems. The numerical solution for the fractional diffusion and fractional advection-dispersion is compared with their respective analytical solutions for different time and space grid refinements. The diffusion simulation exhibited a good fit between the analytical and numerical solutions, with the explicit scheme going from stable to unstable as the time and space refinement changes. The fractional advection-dispersion application produced small deviations from the analytical solution. These deviations, however, are analogous to the numerical dispersions encountered in conventional finite difference solutions of the advection-dispersion equation. The new method is also compared with the traditional L2 method. Notably, an example that involves asymmetrical fractional conditions, a fractional diffusivity that depends on time, and a source term show how the methods compare. Overall, this study assesses the quality and easiness of use of the numerical method.     

Two-dimensional modelling of contaminant transport through saturated porous media using the radial point interpolation method (RPIM)

A new numerical tool is presented which models the two-dimensional contaminant transport through saturated porous media using a meshfree method called the radial point interpolation method (RPIM) with polynomial reproduction. In RPIM, an approximate solution is constructed entirely in terms of a set of nodes and no characterisation of the interrelationship of the nodes is needed. An advection-dispersion equation with sorption is considered to illustrate the applicability of the RPIM. The Galerkin weak form of the governing equation is formulated using two-dimensional meshfree shape functions constructed using thin plate spline radial basis functions. A computer program is developed for the implementation of the RPIM procedure. Three numerical examples are presented and the results are compared with those obtained from the analytical solution and finite element method. The experimental results are also used to validate the approach. The proposed RPIM has generated results with no oscillations and they are insensitive to Peclet constraints.     

Inverse modeling of natural tracer transport in a granite massif with lumped-parameter and physically based models: case study of a tunnel in Czechia

This study deals with numerical modelling of hydraulic and transport phenomena in granite of the Bohemian massif in Bedrichov, Czechia (Czech Republic). Natural tracers represented by stable isotopes δ¹⁸O and δ²H were collected at the tunnel outflow points and nearby catchment and their concentrations were monitored for seven years. The study compared transport simulations by a two-dimensional (2D) physically based model (advection-dispersion) developed in Flow123d software and a simpler lumped-parameter model, calculated with FLOWPC. Both variants were calibrated with UCODE software, either fitting the concentration data alone, or including the tunnel inflow rates in the case of the 2D model calibration (either in separate steps or within a single optimization problem). Since each of the models describes the tracer transport with different parameters, the models were compared based on the mean transit time as a postprocessed quantity. Besides this, two different options for processing the recharge data (input for both models) were evaluated. Calibration and data interpretation were possible for three of the four observed places in the tunnel, thus determining the depth limit of applicability of the stable isotopes. The estimates for discharge sampling at 25–35 m depth based on inverse modelling provide reasonable values of mean transit time (20–40 months) for the lumped parameter models, little revising the results of previous studies at the site. The resulting transport parameters of the advection-dispersion model (porosity and dispersivity) are in accordance with the hydrogeological structures present at the sampling sites.

     

一维均质与非均质土柱中溶质迁移的分数微分对流-弥散模拟

分数微分对流-弥散方程(FADE)是模拟溶质迁移问题的新理论,但应用FADE来模拟溶质迁移时能否克服弥散的尺度效应尚待验证。利用长土柱实验资料结合FADE的解析解拟合推求FADE的弥散系数,并分析其与尺度之间的相关关系。研究结果表明,FADE的弥散系数具有随尺度增大而增大的现象,且均质土柱中FADE的弥散系数尺度效应小于非均质土柱中弥散系数尺度效应。在均质土柱中,弥散系数与尺度之间成指数相关关系,在非均质土柱中,弥散系数与尺度之间成幂相关关系。考虑了弥散系数分别与迁移时间和迁移距离呈线性递增两种相关关系,进而分别构建了3种考虑弥散尺度效应的FADE模型,并提出了求解的差分方法。利用上述3种考虑弥散尺度效应的FADE来模拟和预测不同空间位置处的溶质迁移过程。结果表明,对均质土柱中的溶质迁移可得到较好的模拟结果;对于非均质土柱,其模拟结果与实测结果仍然存在一定的差异。     

考虑固结效应的溶质传输研究

传统的孔隙介质水动力学采用对流-扩散方程,研究溶质在流体中的迁移。在这个过程中,孔隙介质被认为是不变形的,因而是一个稳态问题。针对二维情况下孔隙介质变形对溶质传输的影响,给出了考虑孔隙介质固结效应的溶质传输方程,并且探讨了该类问题的求解方法。     

Computational procedures for non-linear three-dimensional analysis with some advanced constitutive laws

Computational procedures for implementing some constitutive models are described and introduced in three-, and two-dimensional finite element procedures; here variable moduli, Drucker–Prager, critical state and cap models are considered. Consistent numerical schemes are presented with applications to a number of example problems. These procedures can provide successful results with advanced constitutive laws for three-dimensional analysis of a wide range of non-linear problems.     

Two-dimensional patterns of metamorphic fluid flow and isotopic resetting in layered and fractured rocks'pa

One-dimensional advection-dispersion models predict that characteristic δ¹⁸O vs. distance and δ¹⁸O vs. δ¹³C profiles should be produced during isothermal metamorphic fluid flow under equilibrium conditions. However, the patterns of isotopic resetting in rocks that have experienced fluid flow are often different from the predictions. Two-dimensional advection-dispersion simulations in systems with simple geometries suggest that such differences may be as a result of fluid channelling and need not indicate disequilibrium, high dispersivities, or polythermal flow. The patterns of isotopic resetting are a function of: (1) the permeability contrast between more permeable layers ('channels') and less permeable layers ('matrix'); (2) the width and spacing of the channels; (3) the width and spacing of discrete fractures; and (4) the orientation of the pressure gradient with respect to layering. In fractured systems, the efficiency of isotopic transport depends on the fracture aperture and the permeability of the surrounding rock. Resetting initially occurs along and immediately adjacent to the fractures, but with time isotopic resetting because of flow through the rock as a whole increases in importance. Application of the one-dimensional advection-dispersion equations to metamorphic fluid flow systems may yield incorrect estimates of fluid fluxes, intrinsic permeabilities, dispersivities, and permeability contrasts unless fluid flow occurred through zones of high permeability that were separated by relatively impermeable layers.     

Research Progress on Physical Parameterization Schemes in Numerical Weather Prediction Models

Atmospheric physics in numerical weather prediction model which predominantly determines the evolution of atmospheric processes is mainly described by physical parameterization. As a result, the development of physical parameterization has been a hot research issue in the area of numerical prediction for a long time. In this regard, the theoretical background and history of physical parameterization schemes for convection, microphysics, and planetary boundary layer, were reviewed in this study. It is suggested that the advance of physical parameterization for the model with high-resolution grid spaces should be considered as a principle issue for numerical model development in the future. Although the gird spaces in current operational numerical models generally decrease toward 10 km owing to the rapid development of high-performance computation, yet most of these schemes are designed for coarse grid spaces. Because of this kind of deficiency, the theoretical basis of these schemes inevitably faces controversy. Directions for development of physical parameterization were also suggested according to the trends of research in numerical prediction.