Finite volume element method for analysis of unsteady reaction-diffusion problems

A finite volume element method is developed for analyzing unsteady scalar reaction-diffusion problems in two dimensions. The method combines the concepts that are employed in the finite volume and the finite element method together. The finite volume method is used to discretize the unsteady reaction-diffusion equation, while the finite element method is applied to estimate the gradient quantities at cell faces. Robustness and efficiency of the combined method have been evaluated on uniform rectangular grids by using available numerical solutions of the two-dimensional reaction-diffusion problems. The numerical solutions demonstrate that the combined method is stable and can provide accurate solution without spurious oscillation along the high-gradient boundary layers.

An explicit finite volume element method for solving characteristic level set equation on triangular grids

Level set methods are widely used for predicting evolutions of complex free surface topologies,such as the crystal and crack growth,bubbles and droplets deformation,spilling and breaking waves,and two-phase flow phenomena.This paper presents a characteristic level set equation which is derived from the two-dimensional level set equation by using the characteristic-based scheme.An explicit finite volume element method is developed to discretize the equation on triangular grids.Several examples are presented to demonstrate the performance of the proposed method for calculating interface evolutions in time.The proposed level set method is also coupled with the Navier-Stokes equations for two-phase immiscible incompressible flow analysis with surface tension.The Rayleigh-Taylor instability problem is used to test and evaluate the effectiveness of the proposed scheme.

Alternating Direction Finite Volume Element Methods for Three-Dimensional Parabolic Equations

This paper presents alternating direction finite volume element methods for three-dimensional parabolic partial differential equations and gives four computational schemes, one is analogous to Douglas finite difference scheme with second-order splitting error, the other two schemes have third-order splitting error, and the last one is an extended LOD scheme. The L2 norm and H1 semi-norm error estimates are obtained for the first scheme and second one, respectively. Finally, two numerical examples are provided to illustrate the efficiency and accuracy of the methods.

A FINITE VOLUME ELEMENT METHOD FOR THERMAL CONVECTION PROBLEMS

Consider the finite volume element method for the thermal convection problem with the infinite Prandtl number. The author uses a conforming piecewise linear function on a fine triangulation for velocity and temperature, and a piecewise constant function on a coarse triangulation for pressure. For general triangulation the optimal order H1 norm error estimates are given.

Finite Volume Element Predictor-corrector Method for a Class of Nonlinear Parabolic Systems

A finite volume element predictor-corrector method for a class of nonlinear parabolic system of equations is presented and analyzed. Suboptimal L^2 error estimate for the finite volume element predictor-corrector method is derived. A numerical experiment shows that the numerical results are consistent with theoretical analysis.

Biquartic Finite Volume Element Method Based on Lobatto-Guass Structure

In this paper, a biquartic finite volume element method based on Lobatto- Guass structure is presented for variable coefficient elliptic equation on rectangular partition. Not only the optimal H1 and L2 error estimates but also some superconvergent properties are available and could be proved for this method. The numerical results obtained by this finite volume element scheme confirm the validity of the theoretical analysis and the effectiveness of this method.

Quadratic Finite Volume Element Schemes over Triangular Meshes for a Nonlinear Time-Fractional Rayleigh-Stokes Problem

In this article,we study a 2D nonlinear time-fractional Rayleigh-Stokes problem,which has an anomalous subdiffusion term,on triangular meshes by quadratic finite volume element schemes.Time-fractional derivative,defined by Caputo fractional derivative,is discretized through L2−1σformula,and a two step scheme is used to approximate the time first-order derivative at time tn−α/2,where the nonlinear term is approximated by using a matching linearized difference scheme.A family of quadratic finite volume element schemes with two parameters are proposed for the spatial discretization,where the range of values for two parameters areβ1∈(0,1/2),β2∈(0,2/3).For testing the precision of numerical algorithms,we calculate some numerical examples which have known exact solution or unknown exact solution by several kinds of quadratic finite volume element schemes,and contrast with the results of an existing quadratic finite element scheme by drawing diversified comparison plots and showing the detailed data of L2 error results and convergence orders.Numerical results indicate that,L2 error estimate of one scheme with parameters β_(1)=(3−√3)/6,β2=(6+√3−√21+6√3)/9 is O(h^(3)+△t^(2)),and L^(2) error estimates of other schemes are O(h^(2)+△t^(2)),where h and △t denote the spatial and temporal discretization parameters,respectively.

An ADI Finite Volume Element Method for a Viscous Wave Equation with Variable Coefficients

Based on rectangular partition and bilinear interpolation,we construct an alternating-direction implicit(ADI)finite volume element method,which combined the merits of finite volume element method and alternating direction implicit method to solve a viscous wave equation with variable coefficients.This paper presents a general procedure to construct the alternating-direction implicit finite volume element method and gives computational schemes.Optimal error estimate in L2 norm is obtained for the schemes.Compared with the finite volume element method of the same convergence order,our method is more effective in terms of running time with the increasing of the computing scale.Numerical experiments are presented to show the efficiency of our method and numerical results are provided to support our theoretical analysis.

A STABILIZED CRANK-NICOLSON MIXED FINITE VOLUME ELEMENT FORMULATION FOR THE NON-STATIONARY PARABOLIZED NAVIER-STOKES EQUATIONS

A time semi-discrete Crank-Nicolson （CN） formulation with second-order time accuracy for the non-stationary parabolized Navier-Stokes equations is firstly established. And then, a fully discrete stabilized CN mixed finite volume element （SCNMFVE） formu- lation based on two local Gaussian integrals and parameter-free with the second-order time accuracy is established directly from the time semi-discrete CN formulation so that it could avoid the discussion for semi-discrete SCNMFVE formulation with respect to spatial wriables and its theoretical analysis becomes very simple. Finally, the error estimates of SCNMFVE solutions are provided.

Meshfree Finite Volume Element Method for Constrained Optimal Control Problem Governed by Random Convection Diffusion Equations

In this paper,we investigate a stochastic meshfree finite volume element method for an optimal control problem governed by the convection diffusion equations with random coefficients.There are two contributions of this paper.Firstly,we establish a scheme to approximate the optimality system by using the finite volume element method in the physical space and the meshfree method in the probability space,which is competitive for high-dimensional random inputs.Secondly,the a priori error estimates are derived for the state,the co-state and the control variables.Some numerical tests are carried out to confirm the theoretical results and demonstrate the efficiency of the proposed method.

Finite Volume Element Method for Solving the Elliptic Neumann Boundary Control Problems

Solving optimization problems with partial differential equations constraints is one of the most challenging problems in the context of industrial applications. In this paper, we study the finite volume element method for solving the elliptic Neumann boundary control problems. The variational discretization approach is used to deal with the control. Numerical results demonstrate that the proposed method for control is second-order accuracy in the L2 (Γ) and L∞ (Γ) norm. For state and adjoint state, optimal convergence order in the L2 (Ω) and H1 (Ω) can also be obtained.

Maximum Norm Estimates for Finite Volume Element Method for Non-selfadjoint and Indefinite Elliptic Problems

In this paper, we establish the maximum norm estimates of the solutions of the finite volume element method （FVE） based on the P1 conforming element for the non-selfadjoint and indefinite elliptic problems.

Uniform Convergence for Finite Volume Element Method for Non-selfadjoint and Indefinite Elliptic Problems

In this paper, we prove the existence, uniqueness and uniform convergence of the solution of finite volume element method based on the P1 conforming element for non-selfadjoint and indefinite elliptic problems under minimal elliptic regularity assumption.

Mixed Finite Volume Element Method for Vibration Equations of Beam with Structural Damping

In this paper, for the initial and boundary value problem of beams with structural damping, by introducing intermediate variables, the original fourth-order problem is transformed into second-order partial differential equations, and the mixed finite volume element scheme is constructed, and the existence, uniqueness and convergence of the scheme are analyzed. Numerical examples are provided to confirm the theoretical results. In the end, we test the value of δ to observe its influence on the model.

A Quadratic Serendipity Finite Volume Element Method on Arbitrary Convex Polygonal Meshes

Based on the idea of serendipity element,we construct and analyze the first quadratic serendipity finite volume element method for arbitrary convex polygonalmeshes in this article.The explicit construction of quadratic serendipity element shape function is introduced from the linear generalized barycentric coordinates,and the quadratic serendipity element function space based on Wachspress coordinate is selected as the trial function space.Moreover,we construct a family of unified dual partitions for arbitrary convex polygonal meshes,which is crucial to finite volume element scheme,and propose a quadratic serendipity polygonal finite volume element method with fewer degrees of freedom.Finally,under certain geometric assumption conditions,the optimal H1 error estimate for the quadratic serendipity polygonal finite volume element scheme is obtained,and verified by numerical experiments.

Mesh Conditions of the Preserving-Maximum-Principle Linear Finite Volume Element Method for Anisotropic Diffusion-Convection-Reaction Equations

We develop mesh conditions for linear finite volume element approximations of anisotropic diffusionconvectionreaction problems to satisfy the discrete maximum principle.We obtain the sufficient conditions to gurantee the both upper and lower bounds of the numerical solution when each angle of arbitrary triangle is O(∥q∥_∞h+∥g∥_∞h~2)-acute and h is small enough,where h denotes the mesh size,q and g are coefficients of the convection and reaction terms,respectively.To deal with the convection-dominated problems,we use the upwind triangle technique.For such scheme,the mesh condition can be sharper to O(∥g∥_∞h~2)-acute.Some numerical examples are presented to demonstrate the theoretical results.

Improvement of predictions of petrophysical transport behavior using three-dimensional finite volume element model with micro-CT images

Due to the intricate structure of porous media, the macroscopic petrophysical transport properties such as the permeability and the saturation used for the reservoir prediction also show a very complex nature and are difficult to obtain. Thus, a better understanding of the influence of the rock structure on the petrophysical transport properties is important. In this paper, we present a universal finite volume element modeling approach to reconstruct the three dimensional pore models from the micro-CT images based on the commercial software Mimics and ICEM, prior to the pore network model based on some basic assumptions. Moreover, tetra finite volume elements are piled up to realize the geometry reconstruction and the meshing process. Compared with the former methods, this process avoids the tremendously large storage requirement for the reconstructed porous geometry and the failures of meshing these complex polygon geometries, and at the same time improves the predictions of petrophysical transport behaviors. The model is tested on two Berea sandstones, four sandstone samples, two carbonate samples, and one Synthetic Silica. Single- and two phase flow simulations are conducted based on the Navier-Stokes equations in the Fluent software. Good agreements are obtained on both the network structures and predicted single- and two- phase transport properties against benchmark experimental data.

Two-Grid Immersed Finite Volume Element Methods for Semi-Linear Elliptic Interface Problems with Non-Homogeneous Jump Conditions

In this paper,we propose an immersed finite volume element method for solving the semi-linear elliptic interface problems with non-homogeneous jump conditions.Furthermore,two-grid techniques are used to improve the computational efficiency.In this way,we only need to solve a non-linear system on the coarse grid,and a linear system on the fine grid.Numerical results illustrate that the proposed method can solve the semi-linear elliptic interface problems efficiently.Approximate secondorder accuracy for the solution in the L¥norm can be obtained for the considered examples.

A Symmetric Characteristic Finite Volume Element Scheme for Nonlinear Convection-Diffusion Problems

In this paper, we implement alternating direction strategy and construct a symmetric FVE scheme for nonlinear convection-diffusion problems. Comparing to general FVE methods, our method has two advantages. First, the coefficient matrices of the discrete schemes will be symmetric even for nonlinear problems. Second, since the solution of the algebraic equations at each time step can be inverted into the solution of several one-dimensional problems, the amount of computation work is smaller. We prove the optimal H1-norm error estimates of order O（△t2 ＋ h） and present some numerical examples at the end of the paper.

Reduced-order proper orthogonal decomposition extrapolating finite volume element format for two-dimensional hyperbolic equations

This paper is concerned with establishing a reduced-order extrapolating fi- nite volume element （FVE） format based on proper orthogonal decomposition （POD） for two-dimensional （2D） hyperbolic equations. For this purpose, a semi discrete variational format relative time and a fully discrete FVE format for the 2D hyperbolic equations are built, and a set of snapshots from the very few FVE solutions are extracted on the first very short time interval. Then, the POD basis from the snapshots is formulated, and the reduced-order POD extrapolating FVE format containing very few degrees of freedom but holding sufficiently high accuracy is built. Next, the error estimates of the reduced-order solutions and the algorithm procedure for solving the reduced-order for- mat are furnished. Finally, a numerical example is shown to confirm the correctness of theoretical conclusions. This means that the format is efficient and feasible to solve the 2D hyperbolic equations.