Finite volume and finite element methods for solving a one-dimensional space-fractional Boussinesq equation

In this paper, a new one-dimensional space-fractional Boussinesq equation is proposed. Two novel numerical methods with a nonlocal operator (using nodal basis functions) for the space-fractional Boussinesq equation are derived. These methods are based on the finite volume and finite element methods, respectively. Finally, some numerical results using fractional Boussinesq equation with the maximally positive skewness and the maximally negative skewness are given to demonstrate the strong potential of these approaches. The novel simulation techniques provide excellent tools for practical problems. These new numerical models can be extended to two- and three-dimensional fractional space-fractional Boussinesq equations in future research where we plan to apply these new numerical models for simulating the tidal water table fluctuations in a coastal aquifer.     

A degenerate reaction-diffusion system modeling atmospheric dispersion of pollutants

We study the global existence and approximation of the solutions to a degenerate reaction-diffusion system modeling photochemical generation and atmospheric dispersion of pollutants.     

Finite time blow-up in a parabolic–elliptic Keller–Segel system with flux dependent chemotactic coefficient

A parabolic–elliptic Keller–Segel system with homogeneous Neumann boundary condition is considered in a radially symmetric domain where and is a -dimensional ball of radius We assert that under a condition on the initial data, radial weak solutions blow-up in finite time when      

Blow-up and decay criteria for a model of chemotaxis

For a system of equations introduced by Jäger and Luckhaus (1992) [6] as a model of chemotaxis, the questions of blow-up and global existence criteria are investigated. Specifically, for a convex region, a lower bound for the blow-up time is derived if the solution blows up, and explicit criteria to ensure non-blow-up are determined.     

Gaussian bounds for degenerate parabolic equations

Let A be a real symmetric, degenerate elliptic matrix whose degeneracy is controlled by a weight w in the A₂ or QC class. We show that there is a heat kernel Wt(x,y) associated to the parabolic equation wut=divA∇u, and Wt satisfies classic Gaussian bounds:

     

有限体积KFVS方法在二维溃坝中的应用

本文采用了基于KFVS格式的有限体积方法 (FVM)求解了控制水流运动的二维浅水方程 ,建立了二维水坝瞬间溃坝的洪水演进模型 .并应用此模型模拟了二维非对称溃坝和对称溃坝情形下坝左下角有障碍物时的洪水波演进过程 .模拟结果表明该数学模型对二维浅水运动的模拟很有效 .     

Finite volume method for solving a one-dimensional parabolic inverse problem

In this paper, finite volume method is used to solve a one-dimensional parabolic inverse problem with source term and Neumann boundary conditions for the first time. Some advantages of this approach are developing difference schemes and maintaining certain properties of the physics of the problems, especially for the treatment of the source term and the unknown boundary conditions. Numerical results show that our method is more effective.     

A posteriori error estimator for finite volume methods

In this paper, first we discuss a technique to compare finite volume method and some well-known finite element methods, namely the dual mixed methods and nonconforming primal methods, for elliptic equations. These both equivalences are exploited to give us a posteriori error estimator for finite volume methods. This estimator is explicitly given, easy to compute and asymptotically exact without any regularity of the solution in unstructured grids.     

Finite volume approximation of degenerate two‐phase flow model with unlimited air mobility

Models of two‐phase flows in porous media, used in petroleum engineering, lead to a coupled system of two equations, one elliptic and the other degenerate parabolic, with two unknowns: the saturation and the pressure. In view of applications in hydrogeology, we construct a robust finite volume scheme allowing for convergent simulations, as the ratio μ of air/liquid mobility goes to infinity. This scheme is shown to satisfy a priori estimates (the saturation is shown to remain in a fixed interval, and a discrete L2(0,T;H1(Ω)) estimate is proved for both the pressure and a function of the saturation), which are sufficient to derive the convergence of a subsequence to a weak solution of the continuous equations, as the size of the discretization tends to zero. We then show that the scheme converges to a two‐phase flow model whose limit, as the mobility of the air phase tends to infinity, is the “quasi‐Richards equation” (Eymard et al., Convergence of two phase flow to Richards model, F. Benkhaldoun, editor, Finite Volumes for Complex Applications IV, ISTE, London, 2005; Eymard et al., Discrete Cont Dynam Syst, 5 (2012) 93–113), which remains available even if the gas phase is not connected with the atmospheric pressure. Numerical examples, which show that the scheme remains robust for high values of μ, are finally given. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2013     

Numerical methods for fourth order nonlinear degenerate diffusion problems

Numerical schemes are presented for a class of fourth order diffusion problems. These problems arise in lubrication theory for thin films of viscous fluids on surfaces. The equations being in general fourth order degenerate parabolic, additional singular terms of second order may occur to model effects of gravity, molecular interactions or thermocapillarity. Furthermore, we incorporate nonlinear surface tension terms. Finally, in the case of a thin film flow driven by a surface active agent (surfactant), the coupling of the thin film equation with an evolution equation for the surfactant density has to be considered. Discretizing the arising nonlinearities in a subtle way enables us to establish discrete counterparts of the essential integral estimates found in the continuous setting. As a consequence, the resulting algorithms are efficient, and results on convergence and nonnegativity or even strict positivity of discrete solutions follow in a natural way. The paper presents a finite element and a finite volume scheme and compares both approaches. Furthermore, an overview over qualitative properties of solutions is given, and various applications show the potential of the proposed approach.     

Global existence and decay properties for a degenerate Keller-Segel model with a power factor in drift term

The following degenerate parabolic system modelling chemotaxis is considered:

(KS)     

Coexistence and optimal control problems for a degenerate predator-prey model

In this paper we present a predator-prey mathematical model for two biological populations which dislike crowding. The model consists of a system of two degenerate parabolic equations with nonlocal terms and drifts. We provide conditions on the system ensuring the periodic coexistence, namely the existence of two non-trivial non-negative periodic solutions representing the densities of the two populations. We assume that the predator species is harvested if its density exceeds a given threshold. A minimization problem for a cost functional associated with this process and with some other significant parameters of the model is also considered.     

多孔介质中可压缩可混溶驱动问题的有限体积元法

有界区域上多孔介质中可压缩可混溶驱动问题由两个非线性抛物型方程耦合而成：压力方程和饱和度方程均是抛物型方程．运用有限体积元法对两个方程进行数值分析,给出了全离散有限体积元格式,并通过详细的理论分析,得到了近似解与原问题真解的最优H^1模误差估计。     

A simple finite volume method for the shallow water equations

We present a new finite volume method for the numerical solution of shallow water equations for either flat or non-flat topography. The method is simple, accurate and avoids the solution of Riemann problems during the time integration process. The proposed approach consists of a predictor stage and a corrector stage. The predictor stage uses the method of characteristics to reconstruct the numerical fluxes, whereas the corrector stage recovers the conservation equations. The proposed finite volume method is well balanced, conservative, non-oscillatory and suitable for shallow water equations for which Riemann problems are difficult to solve. The proposed finite volume method is verified against several benchmark tests and shows good agreement with analytical solutions.     

Congruences for degenerate number sequences

The degenerate Stirling numbers and degenerate Eulerian polynomials are intimately connected to the arithmetic of generalized factorials. In this article, we show that these numbers and similar sequences may in fact be expressed as p-adic integrals of generalized factorials. As an application of this identification we deduce systems of congruences which are analogues and generalizations of the Kummer congruences for the ordinary Bernoulli numbers.     

Superconvergent Nyström and degenerate kernel methods for eigenvalue problems

The Nyström and degenerate kernel methods, based on projections at Gauss points onto the space of (discontinuous) piecewise polynomials of degree ?r-1, for the approximate solution of eigenvalue problems for an integral operator with a smooth kernel, exhibit order 2r. We propose new superconvergent Nyström and degenerate kernel methods that improve this convergence order to 4r for eigenvalue approximation and to 3r for spectral subspace approximation in the case where the kernel is sufficiently smooth. Moreover for a simple eigenvalue, we show that by using an iteration technique, an eigenvector approximation of order 4r can be obtained. The methods introduced here are similar to that studied by Kulkarni in [10] and exhibit the same convergence orders, so a comparison with these methods is worked out in detail. Also, the error terms are analyzed and the obtained methods are numerically tested. Finally, these methods are extended to the case of discontinuous kernel along the diagonal and superconvergence results are also obtained.     

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.     

An enhanced finite volume method to model 2D linear elastic structures

This paper details the evaluation and enhancement of the vertex-centred finite volume method for the purpose of modelling linear elastic structures undergoing bending. A matrix-free edge-based finite volume procedure is discussed and compared with the traditional isoparametric finite element method via application to a number of test-cases. It is demonstrated that the standard finite volume approach exhibits similar disadvantages to the linear Q4 finite element formulation when modelling bending. An enhanced finite volume approach is proposed to circumvent this and a rigorous error analysis conducted. It is demonstrated that the developed finite volume method is superior to both standard finite volume and Q4 finite element methods, and provides a practical alternative to the analysis of bending-dominated solid mechanics problems.     

Two-grid methods for characteristic finite volume element solution of semilinear convection-diffusion equations

Two-grid methods for characteristic finite volume element solutions are presented for a kind of semilinear convection-dominated diffusion equations. The methods are based on the method of characteristics, two-grid method and the finite volume element method. The nonsymmetric and nonlinear iterations are only executed on the coarse grid (with grid size H). And the fine-grid solution (with grid size h) can be obtained by a single symmetric and linear step. It is proved that the coarse grid can be much coarser than the fine grid. The two-grid methods achieve asymptotically optimal approximation as long as the mesh sizes satisfy H = O(h^1/3).     

Finite volume analysis of adaptive beams with piezoelectric sensors and actuators

This paper presents a finite volume (FV) formulation for the free vibration analysis and active vibration control of the smart beams with piezoelectric sensors and actuators. The governing equations based on Timoshenko beam theory are discretized using the finite volume method. For the purpose of forced vibration control of beam structures, the negative velocity feedback controller is designed for the single-input, single-output system. To achieve the best effect, the piezoelectric sensors and actuators are coupled with the host structure in different positions and then the performance of the designed control system is evaluated for each position. In the test examples, first the shear locking free feature of the present formulation is demonstrated. This has been performed by doing static and natural frequency analysis of some reference models. Then, the capability of the proposed method for the prediction of uncontrolled forced vibration response and active vibration control of a beam structure is studied.