Simulation of unsteady fluid filtration caused by the exploitation of underground resources

The land subsidence near an oil-producing region is simulated. The elastic filtration regime (the Terzaghi model) is used. The numerical simulation is performed by the discretization of the boundary value problem in spatial variables with the aid of the finite element method and by the finite difference discretization in time. The numerical algorithm in use is implemented as a software package. The problem of pumping a fluid from a five-layer medium is solved as an illustrative example.     

Acoustic sounding of perforated wellbores by short waves

The reflection and transmission of harmonic waves and waves of finite duration through the boundary of the perforated zone of a cased wellbore filled with a fluid are studied. A model of the plane fluid flow in the well (in a quasi-one-dimensional approximation) and filtration absorption of the fluid by the porous medium surrounding the well is proposed. The effect of the quality of the perforation (length of perforation tunnels) on the evolution of the waves reflected from the boundary of the perforated zone of the well is studied.     

Multi-fluid model of suspension filtration in a porous medium

In the framework of a three-fluid approach, a new model of suspension filtration in a porous medium is constructed with account for the formation of a dense packing of trapped particles with finite permeability and porosity. The following three continua are considered: the carrier fluid, the suspended particles, and the deposited particles. For a one-dimensional transient flow of suspension, a system of equations for the concentrations of the suspended and deposited particles, the suspension velocity, and the pressure is constructed. Two cases of the flow in a porous medium are considered: plane and radial. Numerical solution is found using a finite-difference method. Numerical calculations are shown to be in agreement with an analytical solution for the simplest case of filtration with a constant velocity and constant porosity and permeability. A comparison is performed with the classic filtration models for five sets of experimental data on the contamination of a porous sample. It is shown that near the inlet boundary, where an intense deposition of particles takes place, the new model describes the concentration profile of the deposited particles more accurately than the classical model.     

Finite analytic boundary condition for a higher‐order finite analytic scheme

A higher‐order finite analytic scheme based on one‐dimensional finite analytic solutions is used to discretize three‐dimensional equations governing turbulent incompressible free surface flow. In order to preserve the accuracy of the numerical scheme, a new, finite analytic boundary condition is proposed for an accurate numerical solution of the partial differential equation. This condition has higher‐order accuracy. Thus, the same order of accuracy is used for the boundary. Boundary conditions were formulated and derived for fluid inflow, outflow, impermeable surfaces and symmetry planes. The derived boundary conditions are treated implicitly and updated with the solution of the problem. The basic idea for the derivation of boundary conditions was to use the discretized form of the governing equations for the fluid flow simplified on the boundaries and flow information. To illustrate the influence of the higher‐order effects at the boundaries, another, lower‐order finite analytic boundary condition, is suggested. The simulations are performed to demonstrate the validity of the present scheme and boundary conditions for a Wigley hull advancing in calm water. Copyright © 2005 John Wiley & Sons, Ltd.     

Natural convection in porous media—dual reciprocity boundary element method solution of the Darcy model

This paper describes the solution of a steady state natural convection problem in porous media by the dual reciprocity boundary element method (DRBEM). The boundary element method (BEM) for the coupled set of mass, momentum, and energy equations in two dimensions is structured by the fundamental solution of the Laplace equation. The dual reciprocity method is based on augmented scaled thin plate splines. Numerical examples include convergence studies with different mesh size, uniform and non‐uniform mesh arrangement, and constant and linear boundary field discretizations for differentially heated rectangular cavity problems at filtration with Rayleigh numbers of Ra*=25, 50, and 100 and aspect ratios of A=1/2, 1, and 2. The solution is assessed by comparison with reference results of the fine mesh finite volume method (FVM). Copyright © 2000 John Wiley & Sons, Ltd.     

Reservoir fluid flow through vertical fracture of thick stratum with underlying fluid interface

An exact solution is obtained for the problem of steady-state filtration of a heavy dense incompressible fluid in a thin, infinitely deep, inclined reservoir having a crack of given depth along the reservoir rise. The region of filtration of the lighter liquid (oil) has an impermeable upper boundary in the form of a horizontal fault line. Below the filtration region there is a free boundary, below which lies the region of stationary fluid (bottom water). The interface of the fluids, the fissure profile, and the reservoir fluid flow rate are determined from the solution of the problem on the basis of the given parameters (permeability of the reservoir and of the material filling the fissure, viscosity of the filtering fluid, specific weight of the upper and lower fluids, depth of the fissure, pressure differential between a point at the fissure and a point at the interface of the fluids). In the case when the thin reservoir is a vertical filtering layer, the considered flow is interpreted as the motion of the reservoir fluid through a vertical fissure of a thick reservoir (half-space) in the presence of an underlying fluid interface. The problem is solved in finite form with the aid of known analytic functions using integrals of the Cauchy type. The fundamental solution is first found of the special problem of flow with a point singularity. The fundamental solution is also of independent importance as an extension of the solution of certain known problems [1–4].     

Numerical methods for solving singular integral equations obtained by fracture mechanical analysis of cracked wedge简

The numerical solutions to the singular integral equations obtained by the fracture mechanical analyses of a cracked wedge under three different conditions are considered. The three considered conditions are: (i) a radial crack on a wedge with a non-finite radius under the traction-traction boundary condition, (ii) a radial crack on a wedge with a finite radius under the traction-traction boundary condition, and (iii) a radial crack on a finite radius wedge under the traction-displacement boundary condition. According to the boundary conditions, the extracted singular integral equations have different forms. Numerical methods are used to solve the obtained coupled singular integral equations, where the Gauss-Legendre and the Gauss-Chebyshev polynomials are used to approximate the responses of the singular integral equations. The results are presented in figures and compared with those obtained by the analytical response. The results show that the obtained Gauss-Chebyshev polynomial response is closer to the analytical response.     

A numerical method for coupling three dimensional nonlinear finite elements with elastic boundary elements

This paper systematically deals with the following three problems: (1) Some numerical schemes in coupling FEM- and BEM: including condensation of the boundary integral equation, symmetrization of equivalent stiffness matrix and treatment of traction discontinuity, (2) Coupling of elastoplastic finite elements to elastic boundary elements, (3) Coupling of elasto-viscoplastic finite elements to elastic boundary elements and numerical stability condition.     

离心压缩机叶轮三维有限元强度分析系统

根据离心压缩机叶轮形状的特点,通过有效的网格自动划分和边界条件的处理,采用三维有限元方法进行了强度计算与分析,从而建立了离心叶轮三维有限元强度分析系统。对于叶轮与主轴之间的过盈配合,求出配合处由于变形而产生的弹性支力,以载荷形式施加于配合之处,由于计算时不加入主轴,从而减少了一,通过等厚度圆尖力的计算验证了有限元程序的可靠性,最二,对一个实际的三元流叶轮进行了计算,通过编制的后处理软件,对影响叶轮     

Compact finite difference-Fourier spectral method for three-dimensional incompressible Navier-Stokes equations

A new compact finite difference-Fourier spectral hybrid method for solving the three dimensional incompressible Navier-Stokes equations is developed in the present paper. The fifth-order upwind compact finite difference schemes for the nonlinear convection terms in the physical space, and the sixth-order center compact schemes for the derivatives in spectral space are described, respectively. The fourth-order compact schemes in a single nine-point cell for solving the Helmholtz equations satisfied by the velocities and pressure in spectral space is derived and its preconditioned conjugate gradient iteration method is studied. The treatment of pressure boundary conditions and the three dimensional non-reflecting outflow boundary conditions are presented. Application to the vortex dislocation evolution in a three dimensional wake is also reported. The project supported by the National Natural Science Foundation of China     

Analytical solutions for finite cylindrical dynamic cavity expansion in compressible elastic-plastic materials

Analytical solutions for the dynamic cylindrical cavity expansion in a com-pressible elastic-plastic cylinder with a finite radius are developed by taking into account of the effect of lateral free boundary, which are different from the traditional cavity expan-sion models for targets with infinite dimensions. The finite cylindrical cavity expansion process begins with an elastic-plastic stage followed by a plastic stage. The elastic-plastic stage ends and the plastic stage starts when the plastic wave front reaches the lateral free boundary. Approximate solutions of radial stress on cavity wall are derived by using the Von-Mise yield criterion and Forrestal’s similarity transformation method. The effects of the lateral free boundary and finite radius on the radial stress on the cavity wall are discussed, and comparisons are also conducted with the finite cylindrical cavity expansion in incompressible elastic-plastic materials. Numerical results show that the lateral free boundary has significant influence on the cavity expansion process and the radial stress on the cavity wall of metal cylinder with a finite radius.     

三维非线性有限元与弹性边界元耦合数值方法

本文系统地讨论了以下三个问题:(1) 有限元与边界元耦合中的几个数值问题,其中包括:边界积分方程的凝聚、等效刚度矩阵的对称化及面力不连续的处理;(2) 弹塑性有限元与弹性边界元的耦合;(3) 弹粘塑性有限元与弹性边界元的耦合及数值计算稳定性条件。     

A Dorodnitsyn finite element formulation for laminar boundary layer flow

The Dorodnitsyn boundary later formulation is given a finite element interpretation and found to generate very accurate and economical solutions when combined with an implicit, non-iterative marching scheme in the downstream direction. The algorithm is of order (Δ²u, Δx) whether linear or quadratic elements are used across the boundary layer. Solutions are compared with a Dorodnitsyn spectral formulation and a conventional finite difference formulation for three Falkner-Skan pressure gradient cases and the flow over a circular cylinder. With quadratic elements the Dorodnitsyn finite element formulation is approximately five times more efficient than the conventional finite difference formulation.     

Steady filtration problems with seawater intrusion: macro‐hybrid penalized finite element approximations

Macro‐hybrid penalized finite element approximations are studied for steady filtration problems with seawater intrusion. On the basis of nonoverlapping domain decompositions with vertical interfaces, sections of coastal aquifers are decomposed into subsystems with simpler geometries and small scales, interconnected via transmission conditions of pressure and flux continuity. Corresponding local penalized formulations are derived from the global penalized variational formulation of the two‐free boundary flow problem, with continuity transmission conditions modelled variationally in a dual sense. Then, macro‐hybrid finite element approximations are derived for the system, defined on independent subdomain grids. Parallel relaxation penalty‐duality algorithms are proposed from fixed‐point problem characterizations. Numerical experiments exemplify the macro‐hybrid penalized theory, showing a good agreement with previous primal conforming penalized finite element approximations (Comput. Methods Appl. Mech. Engng. 2000; 190 :609–624). Copyright © 2005 John Wiley & Sons, Ltd.     

A structured tri-tree search method for generation of optimal unstructured finite element grids in two and three dimensions

A new method for generating finite element grids in two and three dimensions is developed. The method is based on a new search tree structure. The search tree is built upon triangles in two dimensions and tetrahedra in three dimensions. The density of elements can be varied throughout the computational domain. Efficient search algorithms for finding points in space and for finding the boundary of the domain have been developed. The speed of the grid algorithm will permit adaptive gridding during computation. The grid algorithm is generally applicable to both hydrodynamic as well as aerodynamic finite element computations. The technique has been used with success for gridding the North Sea-Skagerrak area.     

The boundary value problems for thin coats

We formulate and solve three boundary value problems for the coats of finite and infinite domains in the plane.     

Free axisymmetric vibration of transversely isotropic laminated circular plates

Based on the three dimensional elastic theory, the state equation of the axisymmetric free vibration of transversely isotropic circular plates is established. Taking the advantage of finite Hankel transform, two exact solutions are derived for two boundary conditions, i. e. the rigid-slipping boundary and elastic simply supported boundary. Finally, numerical results are presented and compared with those of FEM. The project is supported by the National Natural Science Foundation of China and the Zhejiang Provincial Natural Science Foundation.     

基于近场动力学的梁结构振动特性分析方法研究

基于欧拉梁理论推导了两自由度梁的常规态型近场动力学(Peridynamics,PD)模型,并提出一种新的自由边界条件施加方法,对不同边界条件的PD梁进行了模态分析,与局部梁的有限元结果进行了对比,验证了模型的收敛性,分析了PD非局部参数对固有频率的影响.结果表明,当近场作用域内物质点密度较小时,PD梁模型非局部性较弱,...     

特征值问题的边界形状灵敏度

研究连续系统振动特征值问题的边界形状灵敏度满足什么方程和边界条件,如何离散化作近似计算结果表明：如果采用相同的有限单元剖分模式,边界形状灵敏度方程和特征值问题方程具有相同的系数矩阵,但前者是非齐次方程,后者是齐次方程;前者需要施加非齐次边界条件,后者施加齐次边界条件。     

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.