The Dual-Reciprocity Boundary Element Analysis for Hydraulically Fractured Shale Gas Reservoirs Considering Diffusion and Sorption Kinetics

This work presents a new application of boundary element method (BEM) to model fluid transport in unconventional shale gas reservoirs with discrete hydraulic fractures considering diffusion, sorption kinetics and sorbed-phase surface diffusion. The fluid transport model consists of two governing partial differential equations (PDEs) written in terms of effective diffusivities for free and sorbed gases, respectively. Boundary integral formulations are analytically derived using the fundamental solution of the Laplace equation for the governing PDEs and Green’s second identity. The domain integrals arising due to the time-dependent function and nonlinear terms are transformed into boundary integrals employing the dual-reciprocity method. This transformation retains the domain-integral-free, boundary-integral-only character of standard BEM approaches. In the proposed solution, the free- and sorbed-gas flow in the shale matrix is solved simultaneously after coupling the fracture flow equation of free gas. Well production performance under the effect of relaxation phenomenon due to delayed responses of sorbed gas under nonequilibrium sorption condition is rigorously captured by imposing the zero-flux condition at fracture–matrix interface for the sorbed-gas transport equation. The validity of proposed solution is verified using several case studies through comparison against a commercial finite-element numerical simulator.

     

An iteration method for integral equations arising from gravity flows with free surface

Gravity potential flows with free surface still present considerable difficulties in non-linear mathematical problems. Previous researchers using analytic function theory could consider only simple geometrical boundaries. Analyzing curvilinear solid boundaries by means of analytic theory is a difficult problem that has not been solved. In this paper, using Muskhelishvili's singular integral equation theory, we turn the gravity flow problem into the Riemann-Hilbert problem. Taking the length of the streamline of the boundary as the independent variable and the velocity potential of the boundary as the function to be determined, we avoid the difficulty that the angle of the curved fixed part is unknown. Following the difference method and the finite element method, we develop a new numerical method that is suitable for complex solid boundaries and overcome the difficulties encountered in applying analytic function theory. Under known discharge, the convergence and stability of the method have been proved and an estimation of error has been obtained. The method has been successfully applied to the calculation of the flow past spillway buckets. The calculated values agree well with the measured results.     

基于边界元方法的气动弹性降阶模型及仿真

传统气动弹性的时域计算耗费了大量时间,为了提高计算效率,本文发展了基于边界元方法的降阶模型技术。首先基于边界元方法建立非定常流场的求解模型,结合特征值分析技术建立了非定常气动力的低阶模型;然后,利用边界元方法建立了气动网格和结构网格之间的信息转换矩阵;最后将非定常气动力降阶模型和结构动力学方程联合,建立了气动弹性系统的低阶状态空间模型。将所发展的降阶模型方法应用于NACA0012翼型的非定常气动力求解中,结果表明降阶模型可以在保证原系统计算精度的同时提高了计算效率;将降阶模型技术应用到三维机翼的气动弹性响应计算中,在系统阶数仅为12阶的情况下可以得到与原系统一致的极限环响应,说明降阶模型技术在求解气动弹性问题中的巨大优势。     

Solution of the Navier–Stokes equations in velocity–vorticity form using a Eulerian–Lagrangian boundary element method

This paper describes the Eulerian–Lagrangian boundary element model for the solution of incompressible viscous flow problems using velocity–vorticity variables. A Eulerian–Lagrangian boundary element method (ELBEM) is proposed by the combination of the Eulerian–Lagrangian method and the boundary element method (BEM). ELBEM overcomes the limitation of the traditional BEM, which is incapable of dealing with the arbitrary velocity field in advection‐dominated flow problems. The present ELBEM model involves the solution of the vorticity transport equation for vorticity whose solenoidal vorticity components are obtained iteratively by solving velocity Poisson equations involving the velocity and vorticity components. The velocity Poisson equations are solved using a boundary integral scheme and the vorticity transport equation is solved using the ELBEM. Here the results of two‐dimensional Navier–Stokes problems with low–medium Reynolds numbers in a typical cavity flow are presented and compared with a series solution and other numerical models. The ELBEM model has been found to be feasible and satisfactory. Copyright © 2000 John Wiley & Sons, Ltd.     

常航速船波势及船波阻力的边界元法

利用满足Laplace方程,线性化自由面条件及无穷远处条件的Havelock兴波源涵数,建立了关于常航速稳态船波势函数的边界积分方程.针对这个积分方程,建立了相应的数值计算方法,编制了一般三维问题的边界元法计算机程序,可用来计算全潜和半潜物体的稳态绕流场及船舶兴波阻力.     

A net inflow method for incompressible viscous flow with moving free surface

A new finite element procedure called the net inflow method has been developed to simulate time-dependent incompressible viscous flow including moving free surfaces and inertial effects. As a fixed mesh approach with triangular element, the net inflow method can be used to analyse the free surface flow in both regular and irregular domains. Most of the empty elements are excluded from the computational domain, which is adjusted successively to cover the entire region occupied by the liquid. The volume of liquid in a control volume is updated by integrating the net inflow of liquid during each iteration. No additional kinetic equation or material marker needs to be considered. The pressure on the free surface and in the liquid region can be solved explicitly with the continuity equation or implicitly by using the penalty function method. The radial planar free surface flow near a 2D point source and the dam-breaking problem on either a dry bed or a still liquid have been analysed and presented in this paper. The predictions agree very well with available analytical solutions, experimental measurements and/or other numerical results.     

A marker-and-cell approach to viscoelastic free surface flows using the PTT model

This work is concerned with the development of a numerical method capable of simulating two-dimensional viscoelastic free surface flows governed by the non-linear constitutive equation PTT (Phan-Thien–Tanner). In particular, we are interested in flows possessing moving free surfaces. The fluid is modelled by a marker-and-cell type method and employs an accurate representation of the fluid surface. Boundary conditions are described in detail and the full free surface stress conditions are considered. The PTT equation is solved by a high order method which requires the calculation of the extra-stress tensor on the mesh contour. The equations describing the numerical technique are solved by the finite difference method on a staggered grid. In order to validate the numerical method fully developed flow in a two-dimensional channel was simulated and the numerical solutions were compared with known analytic solutions. Convergence results were obtained throughout by using mesh refinement. To demonstrate that complex free surface flows using the PTT model can be computed, extrudate swell and a jet flowing onto a rigid plate were simulated.     

Numerical simulations of inviscid and viscous free‐surface flows with application to nonlinear water waves

The purpose of the present study is to establish a numerical model appropriate for solving inviscid/viscous free‐surface flows related to nonlinear water wave propagation. The viscous model presented herein is based on the Navier–Stokes equations, and the free‐surface is calculated through an arbitrary Lagrangian–Eulerian streamfunction‐vorticity formulation. The streamfunction field is governed by the Poisson equation, and the vorticity is obtained on the basis of the vorticity transport equation. For computing the inviscid flow the Laplace streamfunction equation is used. These equations together with the respective (appropriate) fully nonlinear free‐surface boundary conditions are solved using a finite difference method. To demonstrate the model feasibility, in the present study we first simulate collision processes of two solitary waves of different amplitudes, and compute the phenomenon of overtaking of such solitary waves. The developed model is subsequently applied to calculate (both inviscid and the viscous) flow field, as induced by passing of a solitary wave over submerged rectangular structures and rigid ripple beds. Our study provides a reasonably good understanding of the behavior of (inviscid/viscous) free‐surface flows, within the framework of streamfunction‐vorticity formulation. The successful simulation of the above‐mentioned test cases seems to suggest that the arbitrary Lagrangian–Eulerian/streamfunction‐vorticity formulation is a potentially powerful approach, capable of effectively solving the fully nonlinear inviscid/viscous free‐surface flow interactions. Copyright © 2009 John Wiley & Sons, Ltd.     

可压缩流场中气泡脉动数值模拟

在应用边界元方法对气泡动力学的研究中, 绝大多数模型是建立在不压缩势流理论基础之上, 针对可压缩流场中气泡运动特性的研究很少. 从波动方程出发, 分别在气泡运动前期和后期对波动方程进行简化, 得到气泡运动局部和全局简化方程, 采用双渐进方法对简化方程进行匹配, 提出了考虑流场可压缩性的非球状气泡运动模型. 该模型的计算结果与Prospertti 等的解析结果吻合很好, 气泡脉动最大半径和内部最大压力随气泡脉动逐渐减小. 基于该模型对比了自由场中药包爆炸考虑可压缩性与不考虑可压缩性的计算结果, 发现考虑可压缩性气泡射流速度较小, 随后基于该模型计算了刚性边界下气泡的运动特性.     

A coupled boundary element–finite difference model of surface wave motion over a wall turbulent flow

An effective numerical technique is presented to model turbulent motion of a standing surface wave in a tank. The equations of motion for turbulent boundary layers at the solid surfaces are coupled with the potential flow in the bulk of the fluid, and a mixed BEM–finite difference technique is used to model the wave motion and the corresponding boundary layer flow. A mixing‐length theory is used for turbulence modelling. The model results are in good agreement with previous physical and numerical experiments. Although the technique is presented for a standing surface wave, it can be easily applied to other free surface problems. Copyright © 2005 John Wiley & Sons, Ltd.     

Development of the interfacial air layer in the non-aerated region of high-velocity spillway flows. Instabilities growth,entrapped air and influence on the self-aeration onset

Self-aeration is traditionally explained by the water turbulent boundary layer outer edge intersection with the free surface. This paper presents a discussion on the commonly accepted hypothesis behind the computation of the critical point of self-aeration in spillway flows and a new formulation is proposed based on the existence of a developing air flow over the free surface. Upstream of the inception point of self-aeration, some surface roughening has been often reported in previous studies which consequently implies some entrapped air transport and air–water flows coupling. Such air flow is proven in this study by presenting measured air velocities and computing the air boundary layer thickness for a 1V:2H smooth chute flow. Additionally, the growth rate of free surface waves has been analysed by means of Ultrasonic Sensors measurements, obtaining also the entrapped air concentration. High-speed camera imaging has been used for qualitative study of the flow perturbations.     

The Analysis of Dynamic Stress Intensity Factor for Semi-Circular Surface Crack Using Time-Domain BEM Formulation

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｔｒａｎｓｉｅｎｔａｎａｌｙｓｉｓｆｏｒ 3＿Ｄｅｌａｓｔｏｄｙｎａｍｉｃｃｒａｃｋｐｒｏｂｌｅｍｓｉｓａｎｉｍｐｏｒｔａｎｔｃｏｎｔｅｎｔｉｎｆｒａｃｔｕｒｅｍｅｃｈａｎｉｃｓ.Ｂｅｃａｕｓｅｏｆｔｈｅｄｉｆｆｉｃｕｌｔｉｅｓｏｆｔｈｅｐｒｏｂｌｅｍｓ,ｔｈｅｃａｌｃｕｌａｔｉｏｎｏｆＤＳＩＦｆｏｒ 3＿Ｄｃｒａｃｋｉｓａｐｒｏｂｌｅｍｕｎｓｏｌｖｅｄｃｏｍｐｌｅｔｅｌｙｕｐｔｏｔｈｅｐｒｅｓｅｎｔ.Ｃｈｅｎ[1]ｆｉｒｓｔｌｙｕｓｅｄｔｈｅｆｉｎｉｔｅｄｉｆｆｅｒｅｎｃｅ…     

A finite element model for free surface flows on fixed meshes

In this paper, we present a finite element model for free surface flows on fixed meshes. The main novelty of the approach, compared with typical fixed mesh finite element models for such flows, is that we take advantage of the particularities of free surface flow, instead of considering it a particular case of two‐phase flow. The fact that a given free surface implies a known boundary condition on the interface, allows us to solve the Navier–Stokes equations on the fluid domain uncoupled from the solution on the rest of the finite element mesh. This, together with the use of enhanced integration allows us to model low Froude number flows accurately, something that is not possible with typical two‐phase flow models applied to free surface flow. Copyright © 2007 John Wiley & Sons, Ltd.     

二维复合材料结构模态分析的边界元法

将边界元方法用于分析二维复合材料结构的自由振动模态，利用特解处理体积力（惯性力）仅需静态基本解就可求解问题，对一各向同性悬臂梁，用该法得到的结果与用有限元或各向性边界元法得到的结果符合得很好，但该法可解各向异性问题，对层状复合材料简支梁，用该法得到了数值结果与用一维层状复合材料梁的理论解的比较表明，当结构的长厚比大于２０时，二者符合得很好，当结构的长厚比小于２０时，一维层状材料梁的理论将产生很大的     

Numerical simulation of high‐Reynolds number flow around circular cylinders by a three‐step FEM–BEM model

An innovative computational model, developed to simulate high‐Reynolds number flow past circular cylinders in two‐dimensional incompressible viscous flows in external flow fields is described in this paper. The model, based on transient Navier–Stokes equations, can solve the infinite boundary value problems by extracting the boundary effects on a specified finite computational domain, using the projection method. The pressure is assumed to be zero at infinite boundary and the external flow field is simulated using a direct boundary element method (BEM) by solving a pressure Poisson equation. A three‐step finite element method (FEM) is used to solve the momentum equations of the flow. The present model is applied to simulate high‐Reynolds number flow past a single circular cylinder and flow past two cylinders in which one acts as a control cylinder. The simulation results are compared with experimental data and other numerical models and are found to be feasible and satisfactory. Copyright © 2001 John Wiley & Sons, Ltd.     

一种考虑薄壁散射效应的声学计算模型

采用薄壁边界元/FW-H理论混合方法建立了考虑薄壁声学散射效应的数值计算模型.这种声学计算模型可以预测存在薄壁如风扇机匣、蜗壳等条件下的声波的传播及散射问题.计算模型的建立主要包含噪声源的计算和声源的传播两方面:首先建立FW-H的频域方程,并采用计算流体力学方法计算流场,通过流场数据计算气动噪声源;然后采用薄壁面边界元法计算固壁对声波的散射,并计算声波在固壁散射后的声场分布.数值计算结果和实验结果及经典的叶轮机管道风扇噪声理论进行了对比,结果表明,这种计算模型与理论计算结果及实验结果吻合较好,可以准确的预测机匣壁的散射效应对声源传播的影响.     

Boundary element method for thermoelastic analysis of three-dimensional transversely isotropic solids

Thermal effects are well known to manifest themselves as additional volume integral terms in the direct formulation of the boundary integral equation (BIE) for linear elastic solids when using the boundary element method (BEM). This domain integral has been successfully transformed in an exact manner to surface ones only in isotropy and in 2D anisotropy, thereby restoring the BEM as a truly boundary solution technique. The difficulties with extending it to 3D general anisotropic solids lie in the mathematical complexity of the Green’s function and its derivatives for such materials. These quantities are required items in the BEM formulation. In this paper, the exact, analytical transformation of the volume integral associated with thermal effects to surface ones is achieved for a transversely isotropic material using a similar approach which the authors have previously employed for the same task in BEM for 2D general anisotropy. A numerical scheme, however, needs to be employed to evaluate some of the new terms introduced in the surface integrals that arise from this process here. The mathematical soundness of the formulation is demonstrated by a few examples; the numerical results obtained are checked by alternative means, including those obtained from the commercial FEM code, ANSYS.     

Numerical computations of two-dimensional solitary waves generated by moving disturbances

Two-dimensional solitary waves generated by disturbances moving near the critical speed in shallow water are computed by a time-stepping procedure combined with a desingularized boundary integral method for irrotational flow. The fully non-linear kinematic and dynamic free-surface boundary conditions and the exact rigid body surface condition are employed. Three types of moving disturbances are considered: a pressure on the free surface, a change in bottom topography and a submerged cylinder. The results for the free surface pressure are compared to the results computed using a lower-dimensional model, i.e. the forced Korteweg–de Vries (fKdV) equation. The fully non-linear model predicts the upstream runaway solitons for all three types of disturbances moving near the critical speed. The predictions agree with those by the fKdV equation for a weak pressure disturbance. For a strong disturbance, the fully non-linear model predicts larger solitons than the fKdV equation. The fully non-linear calculations show that a free surface pressure generates significantly larger waves than that for a bottom bump with an identical non-dimensional forcing function in the fKdV equation. These waves can be very steep and break either upstream or downstream of the disturbance.     

An adaptive numerical scheme for solving incompressible 2‐phase and free‐surface flows

In this paper, we present a numerical scheme for solving 2‐phase or free‐surface flows. Here, the interface/free surface is modeled using the level‐set formulation, and the underlying mesh is adapted at each iteration of the flow solver. This adaptation allows us to obtain a precise approximation for the interface/free‐surface location. In addition, it enables us to solve the time‐discretized fluid equation only in the fluid domain in the case of free‐surface problems. Fluids here are considered incompressible. Therefore, their motion is described by the incompressible Navier‐Stokes equation, which is temporally discretized using the method of characteristics and is solved at each time iteration by a first‐order Lagrange‐Galerkin method. The level‐set function representing the interface/free surface satisfies an advection equation that is also solved using the method of characteristics. The algorithm is completed by some intermediate steps like the construction of a convenient initial level‐set function (redistancing) as well as the construction of a convenient flow for the level‐set advection equation. Numerical results are presented for both bifluid and free‐surface problems.