Two-grid virtual element discretization of semilinear elliptic problem

Virtual element method for semilinear elliptic problems on polygonal meshes

A two-grid method for semi-linear elliptic interface problems by partially penalized immersed finite element methods

Convergence analysis of virtual element method for the electric interface model on polygonal meshes with small edges

SummaryIn this paper, we establish the connections between the virtual element method (VEM) and the hourglass control techniques that have been developed since the early 1980s to stabilize underintegrated C0 Lagrange finite element methods. In the VEM, the bilinear form is decomposed into two parts: a consistent term that reproduces a given polynomial space and a correction term that provides stability. The essential ingredients of ‐continuous VEMs on polygonal and polyhedral meshes are described, which reveals that the variational approach adopted in the VEM affords a generalized and robust means to stabilize underintegrated finite elements. We focus on the heat conduction (Poisson) equation and present a virtual element approach for the isoparametric four‐node quadrilateral and eight‐node hexahedral elements. In addition, we show quantitative comparisons of the consistency and stabilization matrices in the VEM with those in the hourglass control method of Belytschko and coworkers. Numerical examples in two and three dimensions are presented for different stabilization parameters, which reveals that the method satisfies the patch test and delivers optimal rates of convergence in the L2 norm and the H1 seminorm for Poisson problems on quadrilateral, hexahedral, and arbitrary polygonal meshes. Copyright © 2015 John Wiley & Sons, Ltd.

• Weak formulation of configurational forces (CFs) in the VEM is derived. • Automated mesh refinement based on CFs improves VEM approximations. • Accurate crack tip loading analyses are demonstrated. • Internal moments in virtual elements must not be disregarded in fracture analyses. • Efficient and accurate crack growth simulations with non-convex crack tip elements. The Virtual Element Method (VEM) constitutes a generalization of the Finite Element Method, allowing for spatial discretizations with arbitrary, even non-convex polygonal elements. This offers significant advantages, particularly for problems involving moving discontinuities, local stress concentrations, or multiple phases with distinctly different dimensions. The theory of configurational forces (CF) represents a very flexible basis to quantify driving forces for virtual defect displacements within a thermodynamical framework, and at the same time provides an indicator for the local mesh quality in numerical calculations. In combination, the VEM and the CF theory thus constitute an effective basis for numerical simulations of crack growth within the context of a classical approach with strong discontinuities and sharp crack tips. Challenges to cope with are attributable to the lacking knowledge of shape functions within virtual elements, the numerical integration on arbitrary polygonal elements, and the requirement of an automated generation of polygonal meshes. The weak form of the energy-momentum balance is derived from a variational approach, uniquely separating nodal CFs acting at free surfaces from those driving the defect. Furthermore, internal degrees of freedom, which are inherent to higher-order virtual elements, have to be interpreted within the context of CFs. Unphysical, so-called spurious, CFs are exploited for an improved mesh design at crack tips, eventually generating accurate results of crack tip loading and crack path simulation.

Virtual element methods (VEM) are the latest evolution of the Mimetic Finite Difference Method and can be considered to be more close to the Finite Element approach. They combine the ductility of mimetic finite differences for dealing with rather weird element geometries with the simplicity of implementation of Finite Elements. Moreover, they make it possible to construct quite easily high-order and high-regularity approximations (and in this respect they represent a significant improvement with respect to both FE and MFD methods). In the present paper we show that, on the other hand, they can also be used to construct DG-type approximations, although numerical tests should be done to compare the behavior of DG-VEM versus DG-FEM.

Flow fields in a centrifugal fan for an indoor unit of an air-conditioner were calculated with finite element method-based large eddy simulation (LES) with the aim of predicting fan performance and aerodynamic noise in this study. The numerical simulation code employed throughout the LES was called FrontFlow/blue (FFB). We compared 10M grid [coarse grid] and 60M grid [fine grid] calculation results for investigation of influence of grid resolution. In the fine grid, the number of grid elements in blade-to-blade direction, and of region between the shroud and the bell mouth increased in particular. By calculating with the fine grid, calculated distributions of absolute velocities at blade exit reasonably agreed with experimental results. Because of this, maximum absolute velocity by fine grid near hub decreased as compared to those by coarse grid. Calculated sound pressure level by fine grid was therefore smaller than that by coarse grid, and the overestimation of sound pressure was suppressed by calculating with fine grid. This decrease of the absolute velocity was a first factor for the improvement of calculation accuracy. Moreover, number of captured streaks on the blade, hub, and shroud surfaces by fine grid increased as compared to those by coarse grid. As a result, size of streak by fine grid became smaller than that by coarse grid. Static pressure fluctuations by fine grid on the blade, hub, and shroud surfaces therefore reduced as compared to those by coarse grid. Aerodynamic noise was related to static pressure fluctuations according to Curle’s equation. This reduction of static pressure fluctuations was therefore a second factor for improvement of calculation accuracy.

Running a fine grid model with 107 - 109 of cells is possible using a supercomputer with 103 - 106 of CPUs but may not be always cost-effective. The most cost-effective way is to use a coarse grid model that is much smaller but with static/dynamic profiles very close to the fine grid model. This paper proposes a new layer optimization and upscaling method with the aim for creating a consistent coarse grid model. Unlike the industry's existing layer optimization and upscaling methods, the proposed method performs layer optimization and upscaling fully integrated with the Lorenz coefficient and curves (LCC). Coarse grid layers and their permeabilities are created by minimizing the difference between fine and coarse grid LCCs. The process consists of static and dynamic optimizations. The former is measured by LCC while the latter by pressure, GOR, and water-cut. A new LCC-based permeability upscaling method is developed to preserve the fine grid multiphase flow behaviors. A satisfactory coarse grid model is achieved when both static and dynamic criteria are met. The proposed method has been successfully applied to a giant carbonate oil field in the Caspian Sea that consists of a matrix dominated platform and a fracture/karst dominated rim. Due to the field's complex geology and high H2S content (15%), a dual porosity, dual permeability compositional model has been created to model compositional sour crude flow within and between the matrix and fracture/karst features. The reservoir drive mechanisms are fluid expansion, miscible gas injection and aquifer drive. The reservoir is undersaturated and has an abnormally high initial reservoir pressure. The fine-grid static model contains 104 million cells (370×225×625×2) and the optimized upscaled coarse-grid dynamic model has 8.3 million cells (370×225×50×2). The upscaled model can be run efficiently on the company's existing HPC infrastructure with a maximum of 64 CPUs. Excellent matches of the Lorenz coefficient maps for reservoir total/zones and Lorenz curves at all wells between the fine and coarse grid models have been achieved. Matches on the dynamic variables, e.g., pressure, gas breakthrough time, and GOR growth, in all producers are within the defined acceptable tolerances. The high quality of the static and dynamic matches between the coarse- and fine-grid models confirms that the reservoir properties of the coarse-grid model is very close to the fine-grid model and can be used a base model for history matching and uncertainty analysis.

Conforming virtual element approximations of the two-dimensional Stokes problem

Multiple solutions for resonant semilinear elliptic problems in [formula omitted

The virtual element method (VEM) was developed not too long ago, starting with the paper (Beirão-da-Veiga et al. in SIAM J Numer Anal 51:794–812, 2013) related to elasticity in solid mechanics. The virtual element method allows to revisit the construction of different elements, however has so far not applied to space-time formulations for one-dimensional structural elements like strings, trusses and beams. Here we study several VEM elements suitable for space-time formulations that are build upon the Hamilton’s principle. It will be shown that these elements can be easily incorporated in classical finite element codes since they have the same number of unknowns. Furthermore, we show that the property of VEM to deal with non-conforming meshes is of special interest for holistic space time formulation: VEM formulations allow locally varying time discretizations (time increments) in a natural and efficient way.

Seismic data are routinely used to infer in situ properties of earth materials on many scales, ranging from global studies to investigations of surficial geological formations. While inversion and imaging algorithms utilizing these data have improved steadily, there are remaining challenges that make detailed measurements of the properties of some geologic materials very difficult. For example, the determination of the concentration and orientation of fracture systems is prohibitively expensive to simulate on the fine grid and, thus, some type of coarse-grid simulations are needed. In this paper, we describe a new multiscale finite element algorithm for simulating seismic wave propagation in heterogeneous media. This method solves the wave equation on a coarse grid using multiscale basis functions and a global coupling mechanism to relate information between fine and coarse grids. Using a mixed formulation of the wave equation and staggered discontinuous basis functions, the proposed multiscale methods have the following properties. • The total wave energy is conserved. • Mass matrix is diagonal on a coarse grid and explicit energy-preserving time discretization does not require solving a linear system at each time step. • Multiscale basis functions can accurately capture the subgrid variations of the solution and the time stepping is performed on a coarse grid. We discuss various subgrid capturing mechanisms and present some preliminary numerical results.

In this article, we propose new gradient recovery schemes for the BEM‐based Finite Element Method (BEM‐based FEM) and Virtual Element Method (VEM). Supporting general polytopal meshes, the BEM‐based FEM and VEM are highly flexible and efficient tools for the numerical solution of boundary value problems in two and three dimensions. We construct the recovered gradient from the gradient of the finite element approximation via local averaging. For the BEM‐based FEM, we show that, under certain requirements on the mesh, superconvergence of the recovered gradient is achieved, which means that it converges to the true gradient at a higher rate than the untreated gradient. Moreover, we propose a simple and very efficient a posteriori error estimator, which measures the difference between the unprocessed and recovered gradient as an error indicator. Since the BEM‐based FEM and VEM are specifically suited for adaptive refinement, the resulting adaptive algorithms perform very well in numerical examples.

Parallel block preconditioners for virtual element discretizations of the time-dependent Maxwell equations

This paper presents a multidomain (MD) computational approach to multiscale (MS) compositional reservoir simulation amenable to parallel processing. The main objective is to minimize the upscaling process and preserve the pertinent flow characteristics of the fine-scale geological features of the reservoir. The preservation of the fine-scale reservoir characteristics yields higher-accuracy solution compared to an upscaled model. The computing algorithm involves breaking down the fine-grid (FG) computing mesh into an underlying coarse-grid (CG) subdomain. Unlike other MS techniques, the CG subdomain only serves as an assemblage of the FG cells and not as a computational CG overlaying the FG cells. The flow equations are then restructured into a global pressure equation which is solved on the FG using a large time step without solving for the compositions, saturations, etc. The convergence tolerance for the global pressure solution is soft (no iterations) to reduce computing time. The pressure solution is then used to calculate the molar flux of components only at the FG boundaries of the CG subdomain cells. Finally, each CG subdomain subject to its boundary conditions is solved independently using a smaller time step to compute cell pressures, compositions, and phase saturations. This procedure is amenable to parallel processing. This MS scheme was applied to time domain, in which the global pressure solution is obtained with a larger time step than in the local FG solution. This means that the molar flux boundary condition at the CG subdomain boundaries is not updated every time step, thus, buying additional computational advantage. This new technique was tested on several volatile oil and gas-condensate systems and excellent agreement was obtained with FG simulation with about 60% reduction in computing time using a single processor. Much greater time reduction is expected in a multi processing environment.