COARSE-MESH-ACCURACY IMPROVEMENT OF BILINEAR Q4-PLANE ELEMENT BY THE COMBINED HYBRID FINITE ELEMENT METHOD

IntroductionByadding‘nodeless’incompatiblebubblemodesandpreservinggeometriccharacteristicofthevariationalprincipleinmechanics,thecombinedhybridmethod[1～4 ]remarkablelyenhancedcoarse_mesh_accuracyofconventionalquadrilateralelementsoflowerorder.ThequadrilateralplaneelasticelementCH(0_1 )proposedinRef.[3 ]isasuccessfulexample.Followingthegeometricpointofviewinmechanicscombinedwithmathematicalanalysis,anovelexpressionofthecombinedhybridvariationalprincipleisintroducedtoclarifyitsintrinsicmecha…     

Effects of element order and interface reconstruction in FEM/volume‐of‐fluid incompressible flow simulation

In previous studies, the moment‐of‐fluid interface reconstruction method showed dramatic accuracy improvements in static and pure advection tests over existing methods, but this did not translate into an equivalent improvement in volume‐tracked multimaterial incompressible flow simulation using low‐order finite elements. In this work, the combined effects of the spatial discretization and interface reconstruction in flow simulation are examined. The mixed finite element pairs, Q₁Q₀ (with pressure stabilization) and Q₂P _{? 1} are compared. Material order‐dependent and material order‐independent first and second‐order accurate interface reconstruction methods are used. The Q₂P _{? 1} elements show significant improvements in computed flow solution accuracy for single material flows but show reduced convergence using element‐average piecewise constant density and viscosity in volume‐tracked simulations. In general, a refined Q₁Q₀ grid, with better material interface resolution, provided an accuracy similar to the Q₂P _{? 1} element grid with a comparable number of degrees of freedom. Moment‐of‐fluid shows more benefit from the higher‐order accurate flow simulation than the LVIRA, Youngs', and power diagram interface reconstruction methods, especially on unstructured grids, but does not recover the dramatic accuracy improvements it has shown in advection tests. Published 2012. This article is a US Government work and is in the public domain in the USA.     

Iterative stabilization of the bilinear velocity–constant pressure element

Some finite element approximations of incompressible flows, such as those obtained with the bilinear velocity–constant pressure element (Q₁?P₀), are well known to be unstable in pressure while providing reasonable results for the velocity. We shall see that there exists a subspace of piecewise constant pressures that leads to a stable approximation. The main drawback associated with this subspace is the necessity of assembling groups of elements, the so-called ‘macro-elements’, which increases dramatically the bandwidth of the system. We study a variant of Uzawa's method which enables us to work in the desired subspace without increasing the bandwidth of the system. Numerical results show that this method is efficient and can be made to work at a low extra cost. The method can easily be generalized to other problems and is very attractive in three-dimensional cases.     

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.     

Stability and approximability of the 𝒫1–𝒫0 element for Stokes equations

In this paper we study the stability and approximability of the ??₁–??₀ element (continuous piecewise linear for the velocity and piecewise constant for the pressure on triangles) for Stokes equations. Although this element is unstable for all meshes, it provides optimal approximations for the velocity and the pressure in many cases. We establish a relation between the stabilities of the ??₁–??₀ element (bilinear/constant on quadrilaterals) and the ??₁–??₀ element. We apply many stability results on the ??₁–??₀ element to the analysis of the ??₁–??₀ element. We prove that the element has the optimal order of approximations for the velocity and the pressure on a variety of mesh families. As a byproduct, we also obtain a basis of divergence‐free piecewise linear functions on a mesh family on squares. Numerical tests are provided to support the theory and to show the efficiency of the newly discovered, truly divergence‐free, ??₁ finite element spaces in computation. Copyright © 2006 John Wiley & Sons, Ltd.     

Stability of some low‐order approximations for the Stokes problem

Two‐level low‐order finite element approximations are considered for the inhomogeneous Stokes equations. The elements introduced are attractive because of their simplicity and computational efficiency. In this paper, the stability of a Q₁(h)–Q₁(2h) approximation is analysed for general geometries. Using the macroelement technique, we prove the stability condition for both two‐ and three‐dimensional problems. As a result, optimal rates of convergence are found for the velocity and pressure approximations. Numerical results for three test problems are presented. We observe that for the computed examples, the accuracy of the two‐level bilinear approximation is compared favourably with some standard finite elements. Copyright © 2007 John Wiley & Sons, Ltd.     

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.     

A triangular quasi-conforming finite element for transient dynamic analysis

A type of 3 node triangular element is constructed by the Quasi-conforming method, which may be used to solve the equation of a type of inverse problem of wave propagation after Laplace transformation Δu−A ² u=0. The strains in the element are approximated by an exponential function and the string-net function between neighbouring elements is approximated by one dimensional general solution of the equation. Furthermore the strain field satisfies the equation, and therefore in the derivation of the element formulation, no shape function is needed. In this sense, it is a kind of hybrid element. Compared with the ordinary linear triangular element, the new one features higher precision with coarse meshes. Some numerical tests are presented. The project is supported by the National Natural Science Foundation of China.     

On the selective local stabilization of the mixed Q1–P0 element

In this paper we study the stability and performance of the quadrilateral finite element ??₁–??₀ (bili‐ near/constant) for the Stokes equations. We set up a framework to show the stability of the element for a wide range of meshes with macroelement patches. We apply the new theory to show the stability of ??₁–??₀ elements on some previously studied meshes and on some newly suggested meshes. Nevertheless such earlier and newly suggested meshes are not effective in practice, compared to the traditional unstable meshes for the ??₁–??₀ element. The new theory leads naturally to a general idea in treating instability of square ??₁–??₀ elements by the local stabilization on macroelement patches of larger, but fixed sizes. The good performance of the traditional ??₁–??₀ square elements with filtering can be kept in some cases after the local stabilization. Some numerical tests are provided to support the theory and to show the performance of stabilized ??₁–??₀ elements. Copyright © 2007 John Wiley & Sons, Ltd.     

Properties of an Al/(Al<Subscript>2</Subscript>O<Subscript>3</Subscript>+TiB<Subscript>2</Subscript>+ZrB<Subscript>2</Subscript>) hybrid composite manufactured by powder metallurgy and hot pressing

The properties and microstructure of an Al/(Al₂O₃ + TiB₂ + ZrB₂) hybrid composite made by using hot pressing of aluminum combined with different amounts of TiB₂, ZrB₂, and Al₂O₃ powders are studied. The mechanical properties of the composites are investigated on the basis of microhardness and compression tests. The results show that the microstructure of the composites is uniform and the particles are well distributed in the matrix.     

Gradient‐based nodal limiters for artificial diffusion operators in finite element schemes for transport equations

This paper presents new linearity‐preserving nodal limiters for enforcing discrete maximum principles in continuous (linear or bilinear) finite element approximations to transport problems with steep fronts. In the process of algebraic flux correction, the oscillatory antidiffusive part of a high‐order base discretization is decomposed into a set of internodal fluxes and constrained to be local extremum dim inishing. The proposed nodal limiter functions are designed to be continuous and satisfy the principle of linearity preservation that implies the preservation of second‐order accuracy in smooth regions. The use of limited nodal gradients makes it possible to circumvent angle conditions and guarantee that the discrete maximum principle holds on arbitrary meshes. A numerical study is performed for linear convection and anisotropic diffusion problems on uniform and distorted meshes in two space dimensions. Copyright © 2017 John Wiley & Sons, Ltd.     

A finite element method for the one‐dimensional extended Boussinesq equations

A new finite element method for Nwogu's (O. Nwogu, ASCE J. Waterw., Port, Coast., Ocean Eng., 119 , 618–638 (1993)) one‐dimensional extended Boussinesq equations is presented using a linear element spatial discretisation method coupled with a sophisticated adaptive time integration package. The accuracy of the scheme is compared to that of an existing finite difference method (G. Wei and J.T. Kirby, ASCE J. Waterw., Port, Coast., Ocean Eng., 121 , 251–261 (1995)) by considering the truncation error at a node. Numerical tests with solitary and regular waves propagating in variable depth environments are compared with theoretical and experimental data. The accuracy of the results confirms the analytical prediction and shows that the new approach competes well with existing finite difference methods. The finite element formulation is shown to enable the method to be extended to irregular meshes in one dimension and has the potential to allow for extension to the important practical case of unstructured triangular meshes in two dimensions. This latter case is discussed. Copyright © 1999 John Wiley & Sons, Ltd.     

Accuracy and convergence of element-by-element iterative solvers for incompressible fluid flows using penalty finite element model

The ability of two types of Conjugate Gradient like iterative solvers (GMRES and ORTHOMIN) to resolve large-scale phenomena as a function of mesh density and convergence tolerance limit is investigated. The flow of an incompressible fluid inside a sudden expansion channel is analysed using three meshes of 400, 1600 and 6400 bilinear elements. The iterative solvers utilize the element-by-element data structure of the finite element technique to store and maintain the data at the element level. Both the mesh density and the penalty parameter are found to influence the choice of the convergence tolerance limit needed to obtain accurate results. An empirical relationship between the element size, the penalty parameter, and the convergence tolerance is presented. This relationship can be used to predict the proper choice of the convergence tolerance for a given penalty parameter and element size.     

基于非结构网格的TTGC有限元格式的实现及在超声速流动中的应用

基于非结构网格系统,实现了时空三阶精度的TTGC有限元格式,并在三阶TTGC格式上发展了基于人工粘性的激波捕捉技术。在非结构网格下,采用这种方法对若干典型的超声速流动问题(SOD激波管、马赫数为3的前台阶流动以及马赫数为8的高超声速圆柱流动)进行了验证计算。结果表明,TTGC格式分辨率高,在粗糙网格下能够准确的模拟超声速流场中的激波、接触间断等复杂流动现象,并且能有效的控制间断附近的数值色散现象。与传统的有限体积方法相比,本文实现的TTGC有限元格式在模拟超声速流动问题方面具有格式精度高、数值耗散小等优点。     

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.     

On the quadrilateral Q2–P1 element for the Stokes problem

The Q₂ – P₁ approximation is one of the most popular Stokes elements. Two possible choices are given for the definition of the pressure space: one can either use a global pressure approximation (that is on each quadrilateral the finite element space is spanned by 1 and by the global co‐ordinates x and y) or a local approach (consisting in generating the local space by means of the constants and the local curvilinear co‐ordinates on each quadrilateral ξ and η). The former choice is known to provide optimal error estimates on general meshes. This has been shown, as it is standard, by proving a discrete inf–sup condition. In the present paper we check that the latter approach satisfies the inf–sup condition as well. However, recent results on quadrilateral finite elements bring to light a lack in the approximation properties for the space coming out from the local pressure approach. Numerical results actually show that the second choice (local or mapped pressure approximation) is suboptimally convergent. Copyright © 2002 John Wiley & Sons, Ltd.     

Leray-α Regularization of the Smagorinsky-Closed Filtered Equations for Turbulent Jets at High Reynolds Numbers

The article reports on blending of the Leray-α regularization with the conventional Smagorinsky subgrid-scale closure as an option for large-eddy-simulation of turbulent flows at very high Reynolds number on coarse meshes. The model has been tested in the self-similar far-field region of a jet at a range of Reynolds numbers spanning over two decades (4×10³, 4×10⁴ and 4×10⁵) on two very coarse meshes of 2×10⁵ and 3×10⁴ mesh cells. The results are compared with the well-resolved DNS for $Re_D=4\times 10^3$ on 15 million cells and experimental data for higher Re numbers. While the pure Leray-α can fail badly at high Re numbers on very coarse meshes, a blending of the two strategies by adding a small amount of extra-dissipation performs well even at a huge jet Reynolds number of $Re_D=4\times 10^5$ on a very coarse mesh (2×10⁵ cells), despite the ratio of the typical mesh spacing to the Kolmogorov length exceeding 300. It is found that the main prerequisite for successful LES, both for the classic Smagorinsky and the blended Leray-α/Smagorinsky model, is to resolve the shear-length $L_s=\sqrt{\varepsilon/{\cal S}^3}$ (where ${\cal S}$ is the shear-rate modulus), defined by the constraint Δ/L _s ?<?1, where Δ is the typical mesh-cell size. For the mixed Leray-α/Smagorinsky model the regularization parameter should also be related to the shear-length rather than the local mesh size or Reynolds number, for which we propose a guide criterion α?=?0.15÷0.3 L _s .     

A Lagrangian Mixed Subgrid-Scale Model in Generalized Coordinates

This paper presents the formulation of a mixed dynamic subgrid-scale model in non-Cartesian geometries suitable for the study of complex flows. Following the approach developed by Jordan [J. Comput. Phys. 148, 322 (1999)], the variables are first transformed into a contravariant form and then filtered in the computational space. A dynamic localized mixed model, previously developed within the Cartesian framework has been entirely re-formulated for non-orthogonal meshes. The model performance was evaluated by carrying out two tests. First, a plane channel flow at Re_τ = 395 was simulated using both Cartesian and curvilinear grids; the results show that the model formulation is consistent and insensitive to grid distortion, and compares well with the reference data. Then, computations of the turbulent flow over a two-dimensional channel with a wavy wall were performed. Accurate first- and second-order statistics were obtained using relatively coarse grids. This revised version was published online in July 2006 with corrections to the Cover Date.     

STREAMLINE-BASED MATHEMATICAL MODEL FOR CO2 MISCIBLE FLOODING

According to the research theory of improved black oil simulator, a practical mathematical model for C02 miscible flooding was presented. In the model, the miscible process simulation was realized by adjusting oil/gas relative permeability and effective viscosity under the condition of miscible flow. In order to predict the production performance fast, streamline method is employed to solve this model as an alternative to traditional finite difference methods. Based on streamline distribution of steady-state flow through porous media with complex boundary confirmed with the boundary element method (BEM), an explicit total variation diminishing (TVD) method is used to solve the one-dimensional flow problem. At the same time, influences of development scheme, solvent slug size, and injection periods on CO2 drive recovery are discussed. The model has the advantages of less information need, fast calculation, and adaptation to calculate CO2 drive performance of all kinds of patterns in a random shaped porous media with assembly boundary. It can be an effective tool for early stage screening andmiscible oil field.reservoir dynamic management of the CO2 miscible oil field.     

Matched Boundary Extrapolation Solutions for CO<Subscript>2</Subscript> Well-Injection into a Saline Aquifer

The injection of supercritical CO₂ through wells into deep brine reservoirs is a topic of interest for geologic carbon sequestration. The injected CO₂ is predominantly immiscible with the brine and its low density relative to brine leads to strong buoyancy effects. The displacement of brine by CO₂ in general is a multidimensional, complex nonlinear problem that requires numerical methods to solve. The approximations of vertical equilibrium and complete gravity segregation (sharp interface) have been introduced to reduce the complexity and dimensionality of the problem. Furthermore, for the radial displacement process considered here, the problem can be formulated in terms of a similarity variable that reduces spatial and temporal dependencies to a single variable. However, the resulting ordinary differential equation is still nonlinear and exact solutions are not available. The existing analytical solutions are approximations limited to certain parameter ranges that become inaccurate over a large portion of the parameter space. Here, I use a matched boundary extrapolation method to provide much greater accuracy for analytical/semi-analytical approximations over the full parameter range.