Time–space fractional governing equations of transient groundwater flow in confined aquifers: Numerical investigation          下载免费PDF全文

Tongbi Tu  Ali Ercan  M. Levent Kavvas 《水文研究》2018,32(10):1406-1419

In order to model non‐Fickian transport behaviour in groundwater aquifers, various forms of the time–space fractional advection–dispersion equation have been developed and used by several researchers in the last decade. The solute transport in groundwater aquifers in fractional time–space takes place by means of an underlying groundwater flow field. However, the governing equations for such groundwater flow in fractional time–space are yet to be developed in a comprehensive framework. In this study, a finite difference numerical scheme based on Caputo fractional derivative is proposed to investigate the properties of a newly developed time–space fractional governing equations of transient groundwater flow in confined aquifers in terms of the time–space fractional mass conservation equation and the time–space fractional water flux equation. Here, we apply these time–space fractional governing equations numerically to transient groundwater flow in a confined aquifer for different boundary conditions to explore their behaviour in modelling groundwater flow in fractional time–space. The numerical results demonstrate that the proposed time–space fractional governing equation for groundwater flow in confined aquifers may provide a new perspective on modelling groundwater flow and on interpreting the dynamics of groundwater level fluctuations. Additionally, the numerical results may imply that the newly derived fractional groundwater governing equation may help explain the observed heavy‐tailed solute transport behaviour in groundwater flow by incorporating nonlocal or long‐range dependence of the underlying groundwater flow field.  相似文献   

有关边界层问题中一类奇异积分方程新的存在性结果     

石丽莉  张高勋 《成都信息工程学院学报》2006,21(6):891-896

对边界层理论新结果中出现的一类奇异积分方程w（t）=∫1 t（1-s）（λ＋λs＋s）/w（s）ds＋（1-t）∫t 0s/w（s）ds,t∈（0,1）进行讨论，并得出了上述方程在λ∈（-1/2,0）上正解存在性的新结果。  相似文献   

关于Diophantine方程的解     

乐茂华 《广东海洋大学学报》2004,24(4):47-48

对于正整数a ,设S(a)是Smarandache函数。证明了 :方程S(1·2 ) +S(2·3) +… +S(x(x +1) ) =S(x(x +1) (x +2 ) /3)仅有正整数解x =1。  相似文献   

用比电导法研究活性炭的吸附特性     

陈以良 《广东海洋大学学报》1992,(1)

用比电导法研究两种活性炭自稀水溶液中吸附强电解质硫酸铬和硫酸铜的吸附平衡特性。结果表明,在本文的研究条件下,两种活性炭对硫酸铬和硫酸铜都有吸附作用。同时,活性炭吸附硫酸铜的吸附平衡特性可以用 Freundlich 方程式来描述。研究的结果对固体在溶液中的吸附基础理论以及处理工业废水具有一定的意义。  相似文献   

Numerical manifold method for dynamic nonlinear analysis of saturated porous media     

H. W. Zhang  L. Zhou 《国际地质力学数值与分析法杂志》2006,30(9):927-951

It is well known that the Babuska–Brezzi stability criterion or the Zienkiewicz–Taylor patch test precludes the use of the finite elements with the same low order of interpolation for displacement and pore pressure in the nearly incompressible and undrained cases, unless some stabilization techniques are introduced for dynamic analysis of saturated porous medium where coupling occurs between the displacement of solid skeleton and pore pressure. The numerical manifold method (NMM), where the interpolation of displacement and pressure can be determined independently in an element for the solution of u–p formulation, is derived based on triangular mesh for the requirement of high accurate calculations from practical applications in the dynamic analysis of saturated porous materials. The matrices of equilibrium equations for the second‐order displacement and the first‐order pressure manifold method are given in detail for program coding. By close comparison with widely used finite element method, the NMM presents good stability for the coupling problems, particularly in the nearly incompressible and undrained cases. Numerical examples are given to illustrate the validity and stability of the manifold element developed. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

The similitude equation in magnetotelluric inversion     

S. E. Dosso  D. W. Oldenburg 《Geophysical Journal International》1991,106(2):507-509

  相似文献   

Numerical prediction of blast‐induced stress wave from large‐scale underground explosion     

Chengqing Wu  Yong Lu  Hong Hao 《国际地质力学数值与分析法杂志》2004,28(1):93-109

This paper presents a numerical model for predicting the dynamic response of rock mass subjected to large‐scale underground explosion. The model is calibrated against data obtained from large‐scale field tests. The Hugoniot equation of state for rock mass is adopted to calculate the pressure as a function of mass density. A piecewise linear Drucker–Prager strength criterion including the strain rate effect is employed to model the rock mass behaviour subjected to blast loading. A double scalar damage model accounting for both the compression and tension damage is introduced to simulate the damage zone around the charge chamber caused by blast loading. The model is incorporated into Autodyn3D through its user subroutines. The numerical model is then used to predict the dynamic response of rock mass, in terms of the peak particle velocity (PPV) and peak particle acceleration (PPA) attenuation laws, the damage zone, the particle velocity time histories and their frequency contents for large‐scale underground explosion tests. The computed results are found in good agreement with the field measured data; hence, the proposed model is proven to be adequate for simulating the dynamic response of rock mass subjected to large‐scale underground explosion. Extended numerical analyses indicate that, apart from the charge loading density, the stress wave intensity is also affected, but to a lesser extent, by the charge weight and the charge chamber geometry for large‐scale underground explosions. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

非线性Rossby波所满足的三阶Zakharov方程与低频振荡     

罗德海 《成都信息工程学院学报》1993,(4)

本文在基本气流具有水平切变的情况下,利用摄动法导出了非线性Rossby波所满足的三阶Zakharov方程,然后,考虑了基流具有弱切变的情况,通过使用三阶Zakharov方程研究了Rossby波列的第一类不稳定性问题。结果表明:通过非线性作用,大气中的Rossby波列可产生调制不稳定。同时,本文对这种不稳定的区域,增长率和周期进行了详细的计算,并讨论了波振幅、波数、纬度和基流切变对它们的影响,指出Rossby波列的调制不稳定可以激发30～60天的低频振荡。  相似文献   

Viscoelastodynamics of A Stratified, Compressible Planet: Incremental Field Equations and Short- and Long-Time Asymptotes     

Detlef Wolf 《Geophysical Journal International》1991,104(2):401-417

  相似文献   

Restricted Three-Body Problem: An Approximation of its General Solution – Part One – the Manifold of Symmetric Periodic Solutions     

Kostantinos Papadakis  Constantine Goudas 《Astrophysics and Space Science》2006,305(2):99-124

This paper gives the results of a programme attempting to exploit ‘la seule bréche’ (Poincaré, 1892, p. 82) of non-integrable systems, namely to develop an approximate general solution for the three out of its four component-solutions of the planar restricted three-body problem. This is accomplished by computing a large number of families of ‘solutions précieuses’ (periodic solutions) covering densely the space of initial conditions of this problem. More specifically, we calculated numerically and only for μ = 0.4, all families of symmetric periodic solutions (1st component of the general solution) existing in the domain D:(x ₀ ∊ [−2,2],C ∊ [−2,5]) of the (x ₀, C) space and consisting of symmetric solutions re-entering after 1 up to 50 revolutions (see graph in Fig. 4). Then we tested the parts of the domain D that is void of such families and established that they belong to the category of escape motions (2nd component of the general solution). The approximation of the 3rd component (asymmetric solutions) we shall present in a future publication. The 4th component of the general solution of the problem, namely the one consisting of the bounded non-periodic solutions, is considered as approximated by those of the 1st or the 2nd component on account of the `Last Geometric Theorem of Poincaré' (Birkhoff, 1913). The results obtained provoked interest to repeat the same work inside the larger closed domain D:(x ₀ ∊ [−6,2], C ∊ [−5,5]) and the results are presented in Fig. 15. A test run of the programme developed led to reproduction of the results presented by Hénon (1965) with better accuracy and many additional families not included in the sited paper. Pointer directions construed from the main body of results led to the definition of useful concepts of the basic family of order n, n = 1, 2,… and the completeness criterion of the solution inside a compact sub-domain of the (x ₀, C) space. The same results inspired the ‘partition theorem’, which conjectures the possibility of partitioning an initial conditions domain D into a finite set of sub-domains D i that fulfill the completeness criterion and allow complete approximation of the general solution of this problem by computing a relatively small number of family curves. The numerical results of this project include a large number of families that were computed in detail covering their natural termination, the morphology, and stability of their member solutions. Zooming into sub-domains of D permitted clear presentation of the families of symmetric solutions contained in them. Such zooming was made for various values of the parameter N, which defines the re-entrance revolutions number, which was selected to be from 50 to 500. The areas generating escape solutions have being investigated. In Appendix A we present families of symmetric solutions terminating at asymptotic solutions, and in Appendix B the morphology of large period symmetric solutions though examples of orbits that re-enter after from 8 to 500 revolutions. The paper concludes that approximations of the general solution of the planar restricted problem is possible and presents such approximations, only for some sub-domains that fulfill the completeness criterion, on the basis of sufficiently large number of families.  相似文献