A Mixed Finite Volume Element Method for Nonlinear Time Fractional Fourth-Order Reaction–Diffusion Models

<abstract><p>In this paper, a two-grid mixed finite volume element (MFVE) algorithm is presented for the nonlinear time fractional reaction-diffusion equations, where the Caputo fractional derivative is approximated by the classical $ L1 $-formula. The coarse and fine grids (containing the primal and dual grids) are constructed for the space domain, then a nonlinear MFVE scheme on the coarse grid and a linearized MFVE scheme on the fine grid are given. By using the Browder fixed point theorem and the matrix theory, the existence and uniqueness for the nonlinear and linearized MFVE schemes are obtained, respectively. Furthermore, the stability results and optimal error estimates are derived in detailed. Finally, some numerical results are given to verify the feasibility and effectiveness of the proposed algorithm.</p></abstract>

A positive semi-definite problem of three-dimensional incompressible contamination treatment from nuclear waste in porous media is discussed in this paper. The mathematical model is defined by a nonlinear initial-boundary system consisting of partial differential equations. Four important equations (an elliptic equation, two convection-diffusion equations and a heat conductor equation) determine the physical features. Considering the physical natures and computational efficiency, the authors introduce the conservative mixed finite volume element, upwind approximation and multistep difference to solve this system. The pressure and Darcy velocity are computed by a mixed finite volume element. The concentrations and temperature are solved by a combination of upwind approximation, multistep difference and mixed finite volume element. A multistep difference is used for approximating the partial derivative with respect to time. Mixed finite volume element and upwind differences are given for solving the convection-diffusions equations. Numerical dispersion and nonphysical oscillations could be eliminated, and the computational efficiency is improved by using a large time step. Furthermore, a conservative law is preserved and error estimates in $L^2$-norm is obtained. Finally, two numerical experiments are given to show the efficiency and possible applications.

Herein, we mainly focus on developing a new two-grid Crank-Nicolson (CN) mixed finite element (MFE) (TGCNMFE) method for the generalized nonlinear time fractional fourth-order reaction diffusion equation. To do so, by introducing an auxiliary function, the nonlinear time fractional fourth-order reaction diffusion equation is first split into two second-order nonlinear equations. Thereafter, a new time semi-discrete mixed CN (TSDMCN) scheme is constructed through discretizing the time derivative and time fractional derivative by the CN difference quotient, and the existence, steadiness, and errors of the TSDMCN solutions are analysed. Next, a new TGCNMFE method is developed through using two-grid MFE technique to discretize the spacial variables, and the existence, steadiness, and error estimations for the TGCNMFE solutions are discussed. Lastly, the correctness of theory results and the superiority of the TGCNMFE method are verified by some numerical experiments.

Mixed finite volume element-upwind mixed volume element of compressible two-phase displacement and its numerical analysis

A one-dimensional nonlinear fractional filtration equation with the Riemann–Liouville time-fractional derivative is proposed for modeling fluid flow through a porous medium. This equation is derived under an assumption that the fluid has a fractional equation of state in which the fluid density depends on the time-fractional derivative of pressure. The obtained equation belongs to the diffusion-wave type of equations. A case when the order of fractional differentiation is close to an integer number is considered, and a small parameter is introduced into the fractional filtration equation under consideration. An expansion of the Riemann–Liouville time-fractional derivative into the series with respect to this small parameter is obtained. Using this expansion, a first-order approximation of the derived fractional filtration equation is performed. Next, the problem of approximate Lie point symmetry group classification for this approximate nonlinear filtration equation with a small parameter is studied. It is shown that approximate symmetry groups admitted by different realizations of the approximate filtration equation have much more dimensions than the corresponding exact Lie point symmetry groups admitted by unperturbed fractional diffusion-wave equations. Obtained classification results permit to construct approximate invariant solutions for the considered nonlinear time-fractional filtration equations.

Heterogeneous geological features, such as fractures, faults, shale laminae, and cross-stratified beds, can exert significant effects on the flow performance. Accurate simulation of these complexities requires two model functions, i.e. abilities to accurately compute both pressure and velocity and to efficiently handle permeability tensor. In this paper, a flux continuous model incorporating full tensor permeability was developed for two -phase flow of immiscible and incompressible fluids. Flow systems were split into two subsystems of configurations, one for elliptic equation of global pressure and the other for parabolic equation of saturation, and solve both sequentially by means of the mixed finite volume element (MFVE) method. Numerical examples of simulating flow in sand-shale system and fractured reservoirs are presented to prove the effectiveness and robustness of the model demonstrating particular suitability for handling high heterogeneities and accurately resolving steep pressure gradients. The MFVE method can accurately calculate the flow variables, and generate more realistic streamlines than the standard finite difference method for a discontinuous permeability field. Applications to naturally fractured reservoirs validate the use of full-tensor effective permeability for accurate and effective simulation.

In the present article, an efficient operational matrix based on the famous Laguerre polynomials is applied for the numerical solution of two-dimensional non-linear time fractional order reaction–diffusion equation. An operational matrix is constructed for fractional order differentiation and this operational matrix converts our proposed model into a system of non-linear algebraic equations through collocation which can be solved by using the Newton Iteration method. Assuming the surface layers are thermodynamically variant under some specified conditions, many insights and properties are deduced e.g., nonlocal diffusion equations and mass conservation of the binary species which are relevant to many engineering and physical problems. The salient features of present manuscript are finding the convergence analysis of the proposed scheme and also the validation and the exhibitions of effectiveness of the method using the order of convergence through the error analysis between the numerical solutions applying the proposed method and the analytical results for two existing problems. The prominent feature of the present article is the graphical presentations of the effect of reaction term on the behavior of solute profile of the considered model for different particular cases.

An analytic algorithm for time fractional nonlinear reaction–diffusion equation based on a new iterative method

The paper studies qualitatively new equations of moisture transfer, which generalize the Aller and Aller-Lykov equations. The generalization contributes to revealing in the original equations the specific features of the studied massifs, their structure, physical properties, processes occurring in them through the introduction of the notion of the rates of change of the fractal dimension. We have obtained solutions to the constant coefficient difference equations as a system arising when using the method of lines for the equations with a Riemann-Liouville time fractional derivative with boundary conditions of the first kind. A priori estimates are obtained that imply convergence of the obtained solutions to systems of ordinary differential equations with variable fractional coefficients. Numerical tests have been carried out to confirm theoretical results of the study.

Unstructured anisotropic mesh adaptation is known to be an efficient way to control discretization errors in computational fluid dynamics simulations. Method verification is required to provide the confidence for routine use in production analysis. The current work aims at verification of anisotropic mesh adaptation for Reynolds-averaged Navier–Stokes simulations over the ONERA M6 wing. The present verification study is performed using four different flow solvers, three different implementations of the metric field, and three mesh mechanics packages. Two of the flow solvers use stabilized finite element discretizations (FUN3D-SFE and general geometry Navier–Stokes), one uses finite volume discretization (FUN3D-FV), and the last one uses mixed finite volume and finite element discretizations (Wolf). The mesh adaptation is based on an error estimator that aims to control the quadratic error term in the linear interpolation of the Mach number. Two sets of adaptations were performed; the first one controls the interpolation error in the norm, and the second one controls the interpolation error in the norm. Convergence studies were performed on the forces and the pitching moment using all four solvers, and the results are compared with previously verified convergence studies on fixed (nonadapted) meshes. Both the forces and pitching moment on adapted meshes are found to be converging to the fine-mesh values faster than those on fixed meshes. In addition to forces and moments, the convergence of surface pressure and skin-friction coefficients at various measurement locations on the wing are also presented. Adapted-mesh surface pressure distributions agree with the fine fixed mesh pressure distributions. Adapted-mesh skin-friction distributions contain high-frequency noise with mean values approaching the fixed mesh pressure skin-friction distributions.

In this paper, we analyze the stability of the equilibrium point and Hopf bifurcation point in the three-component time-fractional differential equation, which describes the predator-prey interaction between different species. In the dynamics, the classical first-order derivative in time is modelled by either the Caputo or the Atangana-Baleanu fractional derivative of order α,0<α<1. We utilized a fractional version of the Adams-Bashforth formula to discretize these fractional derivatives in time. The results of the linear stability analysis presented are confirmed by computer simulation results.

Periodic property of the time-fractional Kundu–Mukherjee–Naskar equation

Two new finite element schemes and their analysis for modeling of wave propagation in graphene

This paper presents unconditionally stable and conformal FDTD schemes which are based on the alternating-direction implicit finite difference time domain (ADI-FDTD) method for accurate modeling of perfectly electric conducting (PEC) objects. The proposed schemes are formulated within the framework of the matrix-vector notation of the finite integration technique (FIT), which allows a systematic and consistent extension of finite difference solution of Maxwell's equations on dual grids. As possible choices of second-order convergent conformal method, we apply the partially filled cell (PFC) and the uniformly stable conformal (USC) schemes for the ADI-FDTD method. The unconditional stability and the rates of convergence of the proposed conformal ADI-FDTD (CADI-FDTD) schemes are verified by means of numerical examples of waveguide problems.

Residual-based a posteriori error estimates were derived within one unifying framework for lowest-order conforming, nonconforming, and mixed finite element schemes in Carstensen [Numer Math 100:617–637, 2005]. Therein, the key assumption is that the conforming first-order finite element space $$V_h^c$$ annulates the linear and bounded residual l written $$V_h^c \subseteq {\rm ker} \ell$$. That excludes particular nonconforming finite element methods (NCFEMs) on parallelograms in that $$V_h^c \not\subset {\rm ker} \ell$$. The present paper generalises the aforementioned theory to more general situations to deduce new a posteriori error estimates, also for mortar and discontinuous Galerkin methods. The key assumption is the existence of some bounded linear operator $$\Pi: V_h^c \rightarrow V_h^{nc}$$ with some elementary properties. It is conjectured that the more general hypothesis (H1)–(H3) can be established for all known NCFEMs. Applications on various nonstandard finite element schemes for the Laplace, Stokes, and Navier–Lame equations illustrate the presented unifying theory of a posteriori error control for NCFEM.