岩体裂隙网络非稳定渗流分析与数值模拟

针对裂隙岩体的非稳定渗流问题,通过将Darcy定理扩展到包含干区的整个裂隙网络区域,并令潜在溢出边界条件为Signorini型互补边界条件,将湿区上的非稳定渗流问题转化为全域上的一个新的初边值问题。为降低试探函数选取的难度,建立与定义在整个裂隙网络区域上的偏微分方程（PDE）提法等价的抛物型变分不等式（PVI）提法,并给出裂隙网络非稳定渗流分析的有限元数值分析格式和迭代算法,与砂槽模型试验数据的对比分析,验证其有效性。最后,将文中发展的计算方法应用到含复杂裂隙网络的边坡非稳定渗流分析,计算结果很好地反映出边坡内部自由面随库水降落的变化规律,并能准确地描述裂隙网络内部渗流运动特征及流量分布的不均匀性。     

基于离散裂隙网络模型的裂隙水渗流计算

离散裂隙网络模型(Discrete Fracture Network（DFN）)是研究裂隙水渗流最为有效的手段之一。文章根据裂隙几何参数和水力参数的统计分布,利用Monte Carlo随机模拟技术生成二维裂隙网络,基于图论无向图的邻接矩阵判断裂隙网络的连通,利用递归算法提取出裂隙网络的主干网或优势流路径。基于立方定律和渗流连续性方程,利用数值解析法建立了二维裂隙网络渗流模型,分析不同边界条件下裂隙网络中的流体流动。结果表明,该方法可以模拟区域宏观水力梯度和边界条件下,裂隙网络水力梯度方向总的流量,以及节点的水位、节点间的流量和流动方向的变化特征,为区域岩溶裂隙水渗流计算提供了一种实用、可行的方法。      

基于区域分解算法的地下水耦合模型及其应用

根据岩体中不同区域裂隙发育规模的差异分为连续区域和离散区域,在不同区域中使用不同的地下水运动数学模型,应用区域分解算法来解决这类问题。其中连续区域采用了等效连续介质模型,离散区域采用了随机裂隙网络模型,通过区域公共边界上水位和流量连续的条件将两模型耦合求解。将基于区域分解算法的耦合模型应用于锦屏水电站坝址区三维渗流场的模拟中,通过钻孔观测水位和计算水位的对比发现,该方法是有效的,能够应用于实际工程。     

剖面二维裂隙网络渗流计算方法

本文针对剖面裂隙网络渗流问题,采用全区域不变网络分析方法,并引用初流量概念,化非线性分析为线性分析,成功的解决了岩体剖面上有自由面的裂隙网络渗流计算问题。     

改进的有自由面渗流问题的Bathe算法

建议了一个新的有自由面渗流问题的变分不等式提法,该提法通过将潜在出渗面上的边界条件提为Signorini型条件,从而从理论上消除了出渗点的奇性,解决了出渗点的定位问题。同时在离散求解时,通过引进依赖于网格参数的连续型Heaviside函数,克服了Bathe算法中所固有的网格依赖性,提高了这类方法的数值稳定性。     

Signorini型变分不等式方法在实际工程渗流问题中的应用

对Signorini型变分不等式方法定义的罚Heaviside函数进行了改进,建立了自适应罚Heaviside函数,从而使Signorini型变分不等式方法的数值稳定性进一步得到改善,网格依赖性进一步得到克服。通过将Signorini型变分不等式方法应用于深圳市公明水库坝区渗流分析,论证了Signorini型变分不等式方法在复杂强非线性三维渗流问题中的适用性。     

基于变分不等式的地下洞室渗流边界模拟

Signorini型变分不等式在求解有出渗点的渗流自由面问题时,消除了出渗点的奇性,克服了网格的依赖性。在迭代求解过程中多采用约束迭代法,这种数学约束比较严格,对于自由面穿过的单元计算不容易收敛,会造成结果在两种解中震荡。笔者在变分不等式的基础上修改了迭代公式,对数学约束进行了修改,建立了变带宽的迭代方法。通过修改迭代算法提高了Signorini型变分不等式方法的数值稳定性,同时减少了迭代时间。地下厂房开挖后地下水会从洞室的边墙渗出,临界出渗点的确定对分析渗漏量和排水孔效果起到关键作用。通过对工程中开挖边界和排水孔边界的渗流计算模拟分析,证明了改进迭代算法后的Signorini型变分不等式在复杂非线性强的三维渗流计算中收敛性较好。     

岩体裂隙网络渗流变水温影响分析

考虑岩体裂隙渗流变水温影响,推导单裂隙变水温水流近似解析解和有限元解,在此基础上分别建立裂隙二维网络变水温渗流数值求解方程,分别对应裂隙网络变水温渗流分析的近似解析法和子结构法。分析变温水流运动规律发现：（1）单裂隙内水流水头与水力坡降成非线性关系,当水流由高温区向低温区流动时,水头分布曲线为凸曲线,此时按线性渗流简化水头整体偏小;当由低温区向高温区流动时,水头分布曲线为凹曲线,此时按线性渗流简化水头整体偏大。（2）单裂隙内,高水温处水力坡降小,低水温处水力坡降大;裂隙平均水温越高,流速越快;裂隙网络内存在与裂隙宽度相似的温度偏流效应,即交叉节点水流有偏向水流温度高的裂隙流动的趋势。在温度较高和温度梯度较大的区域,应该考虑水流温度变化对渗流场的影响。     

混合似变分不等式解的一个四步迭代算法

利用不动点和预解方程这一技巧,给出一个求解混合似变分不等式的四步迭代算法。在算子T伪单调连续的条件下,即可证明新提出的算法的收敛性,并且所得到的结果可以看作是对先前求解变分不等式算法的推广和改进。     

岩体离散裂隙网络的非饱和渗流数值分析

针对裂隙岩体的非饱和渗流问题,基于离散裂隙网络模型并结合非饱和Darcy定律、Richards方程、非饱和本构模型以及Signorini型饱和-非饱和互补溢出边界,提出了离散裂隙网络非饱和渗流问题的数学模型。采用有限单元法建立了裂隙网络非饱和渗流模型的数值求解格式和对应的迭代算法。通过与矩形坝稳定渗流、一维竖直裂隙非饱和入渗以及室内二维瞬态排水渗流的试验、数值及理论结果对比分析,验证了文中算法的有效性;根据流量等效原则,指出了裂隙网络模型应用于求解连续介质非饱和渗流问题的有效性。验证了该算法对于求解裂隙边坡降雨入渗问题的可靠性,揭示了降雨入渗过程裂隙网络流量分布的非均匀性及裂隙产状对降雨入渗流动具有重要的控制作用。     

排水砂槽渗流试验与Signorini型变分不等式方法验证

排水是岩土体及工程构筑物渗流控制的主要措施之一,其实质是通过在渗流域内形成潜在溢出边界或低水位边界而实现渗流控制的。在排水渗控条件下,渗流场往往具有强烈的边界非线性特征,Signorini型变分不等式方法从理论上为稳定和非稳定排水渗流问题提供了有效的分析方法,但其实际效果还需要得到试验的验证。通过开展含5个排水廊道的排水砂槽模型试验,研究了复杂排水条件下砂土渗流的基本规律,并通过试验数据与数值计算成果的对比分析,论证了Signorini型变分不等式方法的有效性和正确性。试验结果表明,在稳定渗流条件下,排水砂槽上游侧3个排水廊道对渗流控制起主导作用,而下游侧2个廊道则失去排水功能,数值计算与试验成果吻合较好;在非稳定渗流条件下,受测压管精度、砂样均匀性和毛细效应的影响,数值计算与试验成果存在一定偏差,但也较好地揭示了复杂排水条件下砂槽中的非稳定渗流特征。排水砂槽试验结果验证了Signorini型变分不等式方法的有效性和正确性,为复杂排水条件下岩土体及工程构筑物的渗控结构优化设计提供了有效的分析手段。     

A Numerical Framework for Wall Dissolution Modeling

Scallops and flutes are common dissolution rock forms encountered in karst caves and surface streams. Their evolution is only partially understood and no numerical model that simulates their formation has been presented. This work at least partially fills the gap by introducing a numerical approach to simulate the evolution of different initial forms of soluble surfaces embedded in a turbulent fluid. The aim is to analyze wall dissolution phenomena from basic principles and to identify stable profiles. The analysis is based on a finite volume moving boundary method. The underlying mathematical model is a \(k-\epsilon \) turbulent model for fluid flow coupled with turbulent scalar transport. The rock wall is treated as a moving boundary, where the normal wall retreat velocity is proportional to the under-saturation of the boundary fluid cells with respect to the mineral comprising the wall. As the flow time scale is several orders of magnitude smaller than the dissolution time scale, stationary flow field, concentration field and wall propagation velocity are calculated for each iteration. The boundary at all points is then moved by distracting minimal velocity along the entire boundary from the actual velocity at a certain location, and then normalized to the maximum allowed shift, which is equal to half the height of the boundary cell. In this way only deformation of the initial wall is calculated. The method was applied to several different initial profiles. During the evolution, the profiles progressively converged towards stable forms. In this work, a framework is proposed for a computation of the moving boundary problem related to slow dissolution of a soluble surface.     

结点虚流量法理论基础阐释及改进算法

含自由面的无压渗流问题本质上是一类非线性自由边值问题,固定网格的结点虚流量法在全域范围内不断扣除虚域流量贡献,从而使该问题得到求解。它具有网格依赖性小、出逸点收敛快等优点,但其内在理论基础尚未被完全揭示。通过引入互补型约束条件建立了结点虚流量法和Signorini型变分不等式提法的等价性桥梁,在此基础上引入过渡区放大系数 对自由面判别准则进行优化,并以砂槽模型试验为例进行验证。对比结果表明,优化后算法数值稳定性更好,计算结果与试验数据吻合度更高。研究成果为超大规模网格的渗控结构优化设计提供了有效分析手段。     

A numerical formulation with unified unilateral boundary condition for unsaturated flow problems in porous media

This paper proposed a numerical formulation for unsaturated flow problems with nonlinear boundaries of seepage face and soil–atmosphere interface via the concept of parabolic variational inequality (PVI) method. A unified unilateral boundary condition was first proposed to represent the conditions on the seepage face and soil–atmosphere interface boundaries within the partial differential equation (PDE) formulation. A PVI formulation mathematically equivalent to the PDE formulation was then proposed, which automatically transforms the flux part of the unified unilateral boundary condition into the natural boundary condition and eliminates the singularity at seepage points. By discretizing the PVI formulation, a finite element procedure together with an iterative algorithm was suggested. An existing experiment of unsaturated flow in a layered hillside and a laboratory test of unsaturated flow through sand flume performed in this study were used to validate the proposed method, with a good agreement between the measured and computed results and a satisfactory balance of mass being maintained during the simulations. The numerical results also indicated that the problem of mesh dependence associated with unsaturated flow simulations is well addressed with the proposed numerical method. Finally, the process of unsaturated flow in a soil slope with layers of horizontal drains subjected to rainfall/evaporation was further examined. The numerical results reveal that the deployment of drains in a soil slope can significantly lower the pore water pressure around the drains, with the bottom layer drains being most effective in controlling the seepage flow.     

三维地电断面电阻率测深有限元数值模拟

用有限单元法进行了电导率分块均匀的三维点源电场电阻率测深的数值模拟.首先给出了三维构造中点源电场的边值问题、变分问题; 然后, 用有限单元法求解变分问题, 将区域剖分成六面体单元, 在单元中进行三线性函数插值, 将变分方程化为线性代数方程组; 最后解方程组, 得各节点的电位值, 进而计算出地表的视电阻率.对几例较典型的地电模型进行试算, 结果表明本方法是行之有效的      

Physical characteristics of a lava flow determined from thermal measurements at the lava’s surface

We consider the problem about determination of characteristics of a lava flow from the physical parameters measured on its surface. The problem is formulated as an inverse boundary problem for the model simulating the dynamics of a viscous heat-conducting incompressible inhomogeneous fluid, where, on the basis of additional data at one part of the model boundary, the missing conditions at another part of the boundary have to be determined, and then the characteristics of fluid in the entire model domain have to be reconstructed. The considered problem is ill-posed. We develop a numerical approach to the solution of the problem in the case of a steady-state flow. Assuming that the temperature and the heat flow are known at the upper surface of the lava, we determine the flow characteristics inside the lava. We compute model examples and show that the lava temperature and flow velocity can be determined with a high precision when the initial data are smooth or slightly noisy.     

非恒定水流计算的最优控制问题及其变分求解

为综合利用各类可用的信息源(如水流运动规律,观测资料和统计信息),在反问题的框架下构建了非恒定水流计算的变分模型.针对模型预测的可能误差来源,给出了初始条件、边界条件、物理参数三类独立参量及综合控制问题的定义.在此基础上,阐述了变分法求解这类偏微分方程最优化控制问题的基本原理及模型求解步骤.同时,为便于应用参考,选取初始条件、边界条件和物理参数作为控制向量,以实际常用的水流运动微分方程为基础,导出了一、二维非恒定水流计算的通用伴随方程.并以珠江黄浦至大虎段二维潮流计算为例,展示了变分模型的应用.     

A numerical solution to seepage problems with complex drainage systems

Seepage problems with complex drainage systems are commonly encountered in civil engineering, with strong non-linearity. A numerical solution based on the Finite Element Method combining the substructure technique with a variational inequality formulation of Signorini’s type is proposed to solve these problems. The aims of this work are to accurately characterize the boundary conditions of the drainage systems, to reduce the difficulty in mesh generation resulting from the drainage holes with small radius and dense spacing, and to eliminate the singularity at the seepage points and the resultant mesh dependency. Numerical stability and robustness of the proposed method are guaranteed by an adaptive procedure for progressively relaxing the penalized Heaviside function associated with the formulation of the discrete variational inequality. Two challenging numerical examples are presented to validate the effectiveness and robustness of the proposed method.