Permeability up-scaling using Haar Wavelets

In the context of flow in porous media, up-scaling is the coarsening of a geological model and it is at the core of water resources research and reservoir simulation. An ideal up-scaling procedure preserves heterogeneities at different length-scales but reduces the computational costs required by dynamic simulations. A number of up-scaling procedures have been proposed. We present a block renormalization algorithm using Haar wavelets which provide a representation of data based on averages and fluctuations. In this work, absolute permeability will be discussed for single-phase incompressible creeping flow in the Darcy regime, leading to a finite difference diffusion type equation for pressure. By transforming the terms in the flow equation, given by Darcy’s law, and assuming that the change in scale does not imply a change in governing physical principles, a new equation is obtained, identical in form to the original. Haar wavelets allow us to relate the pressures to their averages and apply the transformation to the entire equation, exploiting their orthonormal property, thus providing values for the coarse permeabilities. Focusing on the mean-field approximation leads to an up-scaling where the solution to the coarse scale problem well approximates the averaged fine scale pressure profile.     

Impact of truncation error and numerical scheme on the simulation of the early time growth of viscous fingering

The truncation error associated with different numerical schemes (first order finite volume, second order finite difference, control volume finite element) and meshes (fixed Cartesian, fixed structured triangular, fixed unstructured triangular and dynamically adapting unstructured triangular) is quantified in terms of apparent longitudinal and transverse diffusivity in tracer displacements and in terms of the early time growth rate of immiscible viscous fingers. The change in apparent numerical longitudinal diffusivity with element size agrees well with the predictions of Taylor series analysis of truncation error but the apparent, numerical transverse diffusivity is much lower than the longitudinal diffusivity in all cases. Truncation error reduces the growth rate of immiscible viscous fingers for wavenumbers greater than 1 in all cases but does not affect the growth rate of higher wavenumber fingers as much as would be seen if capillary pressure were present. The dynamically adapting mesh in the control volume finite element model gave similar levels of truncation error to much more computationally intensive fine resolution fixed meshes, confirming that these approaches have the potential to significantly reduce the computational effort required to model viscous fingering.     

Automatic,unstructured mesh generation for tidal calculations in a large domain

An automated procedure is described for the production of unstructured, finite element meshes to perform depth-integrated, hydrodynamic calculations in an ocean-scale, two-dimensional domain. Three relatively coarse meshes with nearly identical boundaries are automatically produced by basing internal size guidelines on a localized truncation error analysis that was performed using results from a highly resolved mesh.

Qualitative and quantitative comparisons of model performance are made at 150 historical tidal stations. The coarsest mesh is shown to meet or exceed the overall accuracy of the other meshes, including a highly resolved mesh that has over six times as many computational points. The automated procedure quickly and easily produces a computationally efficient and accurate finite element mesh that is reproducible. In addition, the methodology is shown to have potential for assessing the importance and accuracy of and bathymetric details and evaluating historical hydrodynamic data.     

Multiscale modeling and simulation of single-crystal MgO through an atomistic field theory

The present work is concerned with the application of an atomistic-continuum field theory (AFT) in modeling and simulation of crystalline materials. Atomistic formulation of the field theory and its finite element implementation are introduced. Single-crystal MgO under mechanical loading is modeled and simulated. With a coarse mesh, the field theory is shown to be able to simulate dynamic and nonlinear behavior of multi-atom crystalline materials without the need of additional numerical treatments. Reducing the finite element mesh to the atomic scale, i.e., the finite element size is equal to the size of the primitive unit cell, atomic-scale critical phenomena, including dislocations nucleation and motion, have been successfully reproduced.     

Consistent hybrid finite volume/element formulations: model and complex viscoelastic flows

The accuracy and consistency of a new cell‐vertex hybrid finite element/volume scheme are investigated for viscoelastic flows. Finite element (FE) discretization is employed for the momentum and continuity equation, with finite volume (FV) applied to the constitutive law for stress. Here, the interest is to explore the consequences of utilizing conventional cell‐vertex methodology for an Oldroyd‐B model and to demonstrate resulting drawbacks in the presence of complex source terms on structured and unstructured grids. Alternative strategies worthy of consideration are presented. It is demonstrated how high‐order accuracy may be achieved in steady state by respecting consistency in the formulation. Both FE and FV spatial discretizations are embedded in the scheme, with FV triangular sub‐cells referenced within parent triangular finite elements. Both model and complex flow problems are selected to quantify and assess accuracy, appealing to analysis and experimental validation. The test problem is that of steady sink flow, a pure extensional flow, which reflects some of the numerical difficulties involved in solving more generalized viscoelastic flows, where both source and flux terms may contribute equally to stress propagation. In addition, a complex transient filament‐stretching flow is chosen to compute the evolution of stress fields within liquid bridges. Shortcomings of the various stress upwinding schemes are discussed in this context, whilst dealing with such free‐surface type problems. Here, stress fluctuation distribution alone is advocated, and a Lax‐scheme is found to deliver accuracy and stability to the computational results, comparing well with the literature. Copyright © 2004 John Wiley & Sons, Ltd.     

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.     

Two-level stabilized finite element method for Stokes eigenvalue problem

A two-level stabilized finite element method for the Stokes eigenvalue problem based on the local Gauss integration is considered.This method involves solving a Stokes eigenvalue problem on a coarse mesh with mesh size H and a Stokes problem on a fine mesh with mesh size h = O(H 2),which can still maintain the asymptotically optimal accuracy.It provides an approximate solution with the convergence rate of the same order as the usual stabilized finite element solution,which involves solving a Stokes eigenvalue problem on a fine mesh with mesh size h.Hence,the two-level stabilized finite element method can save a large amount of computational time.Moreover,numerical tests confirm the theoretical results of the present method.     

Upwind adaptive finite element investigations of the two-dimensional reactive interaction of supersonic gaseous jets

We describe an upwind finite element method aimed at numerically simulating the two-dimensional transonic flow of a reactive gaseous mixture. The method uses in particular a triangular finite element mesh, with an adaptive procedure based on mesh refinement by triangle division, and an upwind non-oscillatory scheme based on an approximate Riemann solver for the evaluation of the convective terms for all species. Results concerning the reactive interaction of two supersonic gaseous jets are presented.     

Experimental and numerical analysis of fretting crack formation based on 3D X-FEM frictional contact fatigue crack model

Nowadays, numerical simulation of 3D fatigue crack growth is easily handled using the eXtended Finite Element Method coupled with level set techniques. The finite element mesh does not need to conform to the crack geometry. Most difficulties associated to complex mesh generation around the crack and the re-meshing steps during the possible propagation are hence avoided. A 3D two-scale frictional contact fatigue crack model developed within the X-FEM framework is presented in this article. It allows the use of a refined discretization of the crack interface independent from the underlying finite element mesh and adapted to the frictional contact crack scale. A stabilized three-field weak formulation is also proposed to avoid possible oscillations in the local solution linked to the LBB condition when tangential slip is occurring. Two basic three-dimensional numerical examples are presented. They aim at illustrating the capacities and the high level of accuracy of the proposed X-FEM model. Stress intensity factors are computed along the crack front. Finally an experimental 3D ball/plate fretting fatigue test with running conditions inducing crack nucleation and propagation is modeled. 3D crack shapes defined from actual experimental ones and fretting loading cycle are considered. This latter numerical simulation demonstrates the model ability to deal with challenging actual complex problems and the possibility to achieve tribological fatigue prediction at a design stage based on the fatigue crack modeling.     

Developing a unified FVE‐ALE approach to solve unsteady fluid flow with moving boundaries

In this study, an arbitrary Lagrangian–Eulerian (ALE) approach is incorporated with a mixed finite‐volume–element (FVE) method to establish a novel moving boundary method for simulating unsteady incompressible flow on non‐stationary meshes. The method collects the advantages of both finite‐volume and finite‐element (FE) methods as well as the ALE approach in a unified algorithm. In this regard, the convection terms are treated at the cell faces using a physical‐influence upwinding scheme, while the diffusion terms are treated using bilinear FE shape functions. On the other hand, the performance of ALE approach is improved by using the Laplace method to improve the hybrid grids, involving triangular and quadrilateral elements, either partially or entirely. The use of hybrid FE grids facilitates this achievement. To show the robustness of the unified algorithm, we examine both the first‐ and the second‐order temporal stencils. The accuracy and performance of the extended method are evaluated via simulating the unsteady flow fields around a fixed cylinder, a transversely oscillating cylinder, and in a channel with an indented wall. The numerical results presented demonstrate significant accuracy benefits for the new hybrid method on coarse meshes and where large time steps are taken. Of importance, the current method yields the second‐order temporal accuracy when the second‐order stencil is used to discretize the unsteady terms. Copyright © 2009 John Wiley & Sons, Ltd.     

An h-adaptive solution of the spherical blast wave problem

Shock waves and contact discontinuities usually appear in compressible flows, requiring a fine mesh in order to achieve an acceptable accuracy of the numerical solution. The usage of a mesh adaptation strategy is convenient as uniform refinement of the whole mesh becomes prohibitive in three-dimensional (3D) problems. An unsteady h-adaptive strategy for unstructured finite element meshes is introduced. Non-conformity of the refined mesh and a bounded decrease in the geometrical quality of the elements are some features of the refinement algorithm. A 3D extension of the well-known refinement constraint for 2D meshes is used to enforce a smooth size transition among neighbour elements with different levels of refinement. A density-based gradient indicator is used to track discontinuities. The solution procedure is partially parallelised, i.e. the inviscid flow equations are solved in parallel with a finite element SUPG formulation with shock capturing terms while the adaptation of the mesh is sequentially performed. Results are presented for a spherical blast wave driven by a point-like explosion with an initial pressure jump of 10⁵ atmospheres. The adapted solution is compared to that computed on a fixed mesh. Also, the results provided by the theory of self-similar solutions are considered for the analysis. In this particular problem, adapting the mesh to the solution accounts for approximately 4% of the total simulation time and the refinement algorithm scales almost linearly with the size of the problem.     

Using Vorticity as an Indicator for the Generation of Optimal Coarse Grid Distribution

An improved vorticity-based gridding technique is presented and applied to create optimal non-uniform Cartesian coarse grid for numerical simulation of two-phase flow. The optimal coarse grid distribution (OCGD) is obtained in a manner to capture variations in both permeability and fluid velocity of the fine grid using a single physical quantity called “vorticity”. Only single-phase flow simulation on the fine grid is required to extract the vorticity. Based on the fine-scale vorticity information, several coarse grid models are generated for a given fine grid model. Then the vorticity map preservation error is used to predict how well each coarse grid model reproduces the fine-scale simulation results. The coarse grid model which best preserves the fine-scale vorticity, i.e. has the minimum vorticity map preservation error is recognized as an OCGD. The performance of vorticity-based optimal coarse grid is evaluated for two highly heterogeneous 2D formations. It is also shown that two-phase flow parameters such as mobility ratio have only minor impact on the performance of the predicted OCGD.     

Nonconforming stabilized combined finite element method for Reissner-Mindlin plate

Based on combination of two variational principles, a nonconforming stabilized finite element method is presented for the Reissner-Mindlin plates. The method is convergent when the finite element space is energy-compatible. Error estimates are derived. In particular, three finite element spaces are applied in the computation. Numerical results show that the method is insensitive to the mesh distortion and has better performence than the MITC4 and DKQ methods. With properly chosen parameters, high accuracy can be obtained at coarse meshes.     

Upscaling of Channel Systems in Two Dimensions Using Flow-Based Grids

A methodology for the gridding and upscaling of geological systems characterized by channeling is presented. The overall approach entails the use of a flow-based gridding procedure for the generation of variably refined grids capable of resolving the channel geometry, a specialized full-tensor upscaling method to capture the effects of permeability connectivity, and the use of a flux-continuous finite volume method applicable to full tensor permeability fields and non-orthogonal grids. The gridding and upscaling procedures are described in detail and then applied to several two-dimensional systems. Significant improvement in the accuracy of the coarse scale models, relative to that obtained using uniform Cartesian coarse scale models, is achieved in all cases. It is shown that, for some systems, improvement results from the use of the flow-based grid, while in other cases the improvement is mainly due to the new upscaling method.     

应变局部化分析的嵌入强间断多尺度有限元法

固体材料的应变局部化行为是导致结构破坏失效的重要因素之一,开展相关数值模拟分析对于结构安全性评估具有重要意义.然而由于材料的非均质和多尺度特性,采用传统数值方法进行求解时通常需要从最小特征尺度离散求解的结构,这将大幅度增加计算规模和成本.针对这一问题,本文提出了一种基于嵌入强间断模型的多尺度有限元方法.该方法从粗细两个尺度离散求解模型,首先在细尺度单元上引入嵌入强间断模型来描述单元间断特性,所附加的跳跃位移自由度则通过凝聚技术进行消除,从而保持细尺度单元刚度阵维度不变.其次,提出了一种增强多节点粗单元技术,其可根据局部化带与粗单元边界相交情况自适应动态地增加粗节点,新构造的增强数值基函数可以捕捉细尺度间断特性,完成物理信息从细单元到粗单元的准确传递以及宏观响应的快速分析;再次,在细尺度解的计算中,将细尺度解分解为降尺度解与单胞局部摄动解,从而消除弹塑性分析时单胞内部的不平衡力.最后,通过两个典型算例分析,并与完全采用细单元的嵌入有限元结果进行对比,验证了所提出算法的正确性与有效性.     

A NEW HYBRID QUADRILATERAL FINITE ELEMENT FOR MINDLIN PLATE

ＡＮＥＷＨＹＢＲＩＤＱＵＡＤＲＩＬＡＴＥＲＡＬＦＩＮＩＴＥＥＬＥＭＥＮＴＦＯＲＭＩＮＤＬＩＮＰＬＡＴＥＣｈｉｎＹｉ（秦奕）（ＴｉａｎｊｉｎＡｒｃｈｉｔｅｃｔｕｒａｌＤｅｓｉｇｎＩｎｓｔｉｔｕｔｅ，Ｔｉａｎｊｉｎ）ＺｈａｎｇＪｉｎｇ－ｙｕ（张敬宇）（Ｉ...     

Unstructured lattice Boltzmann method in three dimensions

Over the last decade, the lattice Boltzmann method (LBM) has evolved into a valuable alternative to continuum computational fluid dynamics (CFD) methods for the numerical simulation of several complex fluid‐dynamic problems. Recent advances in lattice Boltzmann research have considerably extended the capability of LBM to handle complex geometries. Among these, a particularly remarkable option is represented by cell‐vertex finite‐volume formulations which permit LBM to operate on fully unstructured grids. The two‐dimensional implementation of unstructured LBM, based on the use of triangular elements, has shown capability of tolerating significant grid distortions without suffering any appreciable numerical viscosity effects, to second‐order in the mesh size. In this work, we present the first three‐dimensional generalization of the unstructured lattice Boltzmann technique (ULBE as unstructured lattice Boltzmann equation), in which geometrical flexibility is achieved by coarse‐graining the lattice Boltzmann equation in differential form, using tetrahedrical grids. This 3D extension is demonstrated for the case of 3D pipe flow and moderate Reynolds numbers flow past a sphere. The results provide evidence that the ULBE has significant potential for the accurate calculation of flows in complex 3D geometries. Copyright © 2005 John Wiley & Sons, Ltd.     

A parallel monolithic algorithm for the numerical simulation of large‐scale fluid structure interaction problems

A novel parallel monolithic algorithm has been developed for the numerical simulation of large‐scale fluid structure interaction problems. The governing incompressible Navier–Stokes equations for the fluid domain are discretized using the arbitrary Lagrangian–Eulerian formulation‐based side‐centered unstructured finite volume method. The deformation of the solid domain is governed by the constitutive laws for the nonlinear Saint Venant–Kirchhoff material, and the classical Galerkin finite element method is used to discretize the governing equations in a Lagrangian frame. A special attention is given to construct an algorithm with exact total fluid volume conservation while obeying both the global and the local discrete geometric conservation law. The resulting large‐scale algebraic nonlinear equations are multiplied with an upper triangular right preconditioner that results in a scaled discrete Laplacian instead of a zero block in the original system. Then, a one‐level restricted additive Schwarz preconditioner with a block‐incomplete factorization within each partitioned sub‐domains is utilized for the modified system. The accuracy and performance of the proposed algorithm are verified for the several benchmark problems including a pressure pulse in a flexible circular tube, a flag interacting with an incompressible viscous flow, and so on. John Wiley & Sons, Ltd.     

颗粒材料多尺度分析的连接尺度方法

本文提出了耦合细尺度上基于离散颗粒集合体模型的离散单元法（DEM）和粗尺度上基于Cosserat连续体模型的有限元法（FEM）的连接尺度方法（BSM）以研究颗粒材料的力学行为。采用Cosserat连续体模型和FEM模拟的粗尺度域覆盖全域,而采用离散颗粒集合体模型的DEM模拟的细尺度域仅限于需特别关注材料微结构演变和非连续变形行为的局部区域。对这两个区域间的界面提出了适当的界面条件及其实施方案。通过采用适当的连接尺度投影算子,空间离散的粗、细尺度耦合系统多尺度运动方程具有解耦和允许分别求解、因而也允许分别采用不同时间步长对粗、细尺度计算的特点,可极大地提高BSM的计算效率。文中二维地基数值算例结果说明了所陈述方法的可应用性,以及相对基于Cosserat连续体模型的FEM和基于离散颗粒集合体模型的DEM的优越性。      

Numerical Homogenization of the Absolute Permeability Using the Conformal-Nodal and Mixed-Hybrid Finite Element Method

The modeling of hydrocarbon reservoirs and of aquifer-aquitard systems can be separated into two activities: geological modeling and fluid flow modeling. The geological model focuses on the geometry and the dimensions of the subsurface layers and faults, and on its rock types. The fluid flow model focuses on quantities like pressure, flux and dissipation, which are related to each other by rock parameters like permeability, storage coefficient, porosity and capillary pressure. The absolute permeability, which is the relevant parameter for steady single-phase flow of a fluid with constant viscosity and density, is studied here. When trying to match the geological model with the fluid flow model, it generally turns out that the spatial scale of the fluid flow model is built from units that are at least a hundred times larger in volume than the units of the geological model. To counter this mismatch in scales, the fine-scale permeabilities of the geological data model have to be upscaled' to coarse-scale permeabilities that relate the spatially averaged pressure, flux and dissipation to each other. The upscaled permeabilities may be considered as complicated averages, which are derived from the spatially averaged flow quantities in such a way that the continuity equation, Darcy's law and the dissipation equation remain valid on the coarse scale. In this paper the theory of upscaling will be presented from a physical point of view aiming at understanding, rather than mathematical rigorousness. Under the simplifying assumption of spatial periodicity of the fine-scale permeability distributions, homogenization theory can be applied. However, even then the spatial distribution of the permeability is generally so intricate that exact solutions of the homogenized permeability cannot be found. Therefore, numerical approximation methods have to be applied. To be able to estimate the approximation error, two numerical methods have been developed: one based on the conventional nodal finite element method (CN-FEM) and the other based on the mixed-hybrid finite element method (MH-FEM). CN-FEM gives an upper bound for the sum of the diagonal components of the homogenized mobility matrix, while MH-FEM gives a lower bound. Three numerical examples are presented.