On the non-existence of -uniform finite difference methods on uniform meshes for semilinear two-point boundary value problems

In this paper fitted finite difference methods on a uniform mesh with internodal spacing , are considered for a singularly perturbed semilinear two-point boundary value problem. It is proved that a scheme of this type with a frozen fitting factor cannot converge -uniformly in the maximum norm to the solution of the differential equation as the mesh spacing goes to zero. Numerical experiments are presented which show that the same result is true for a number of schemes with variable fitting factors.

     

Existence of generalized homoclinic solutions of a coupled KdV-type Boussinesq system under a small perturbation

This paper considers the coupled KdV-type Boussinesq system with a small perturbation $u_{xx}=6cv-6u-6uv+\varepsilon f(\varepsilon,u,u_{x},v,v_{x}),$ $ v_{xx}=6cu-6v-3u^{2}+\varepsilon g(\varepsilon,u,u_{x},v,v_{x}),$ where $c=1+\mu$, $\mu>0$ and $\varepsilon$ are small parameters. The linear operator has a pair of real eigenvalues and a pair of purely imaginary eigenvalues. We first change this system into an equivalent system with dimension 4, and then show that its dominant system has a homoclinic solution and the whole system has a periodic solution if the perturbation functions $g$ and $h$ satisfy some conditions. By using the contraction mapping theorem, the perturbation theorem, and the reversibility, we theoretically prove that this homoclinic solution, when higher order terms are added, will persist and exponentially approach to the obtained periodic solution (called generalized homoclinic solution) for small $\varepsilon$ and $\mu>0$.     

A nonconforming combination of the finite element and volume methods with an anisotropic mesh refinement for a singularly perturbed convection-diffusion equation

In this paper we formulate and analyze a discretization method for a 2D linear singularly perturbed convection-diffusion problem with a singular perturbation parameter . The method is based on a nonconforming combination of the conventional Galerkin piecewise linear triangular finite element method and an exponentially fitted finite volume method, and on a mixture of triangular and rectangular elements. It is shown that the method is stable with respect to a semi-discrete energy norm and the approximation error in the semi-discrete energy norm is bounded by with independent of the mesh parameter , the diffusion coefficient and the exact solution of the problem.

     

Unconditional Superconvergent Analysis of Quasi-Wilson Element for Benjamin-Bona-Mahoney Equation

This article aims to study the unconditional superconvergent behavior of nonconforming quadrilateral quasi-Wilson element for nonlinear Benjamin Bona Mahoney (BBM) equation. For the generalized rectangular meshes including rectangular mesh, deformed rectangular mesh and piecewise deformed rectangular mesh, by use of the special character of this element, that is, the conforming part (bilinear element) has high accuracy estimates on the generalized rectangular meshes and the consistency error can reach order $O(h^2)$, one order higher than its interpolation error, the superconvergent estimates with respect to mesh size $h$ are obtained in the broken $H^1$-norm for the semi-/ fully-discrete schemes. A striking ingredient is that the restrictions between mesh size $h$ and time step $\tau$ required in the previous works are removed. Finally, some numerical results are provided to confirm the theoretical analysis.     

Trace Averaging Domain Decomposition Method with Nonconforming Finite Elements

1.IntroductionDomaindec0mpositionreferstonumericaJmethodsf0robtainingsoluti0nsofsci-entificandengineeringproblemsbycombiningsoluti0nstoproblemspo8ed0nphysica1subdomains,or,moregeneraJly,byc0mbiningsoluti0nst0appropriatelyconstructedsubproblems.IthasbeenasubjectofintenseinterestreceDtlybecause0fitssultabil-ityforimplementationonhighperformancecomputerarchitectures.Somepapersarelistedinthereferencesherein,wlilchindicatethatmuchprogresshasbeenmadeinthestudyofnonoverlaPdomaindecompositionmethods…     

Two-level iteration penalty and variational multiscale method for steady incompressible flows

In this paper, we study two-level iteration penalty and variational multiscale method for the approximation of steady Navier-Stokes equations at high Reynolds number. Comparing with classical penalty method, this new method does not require very small penalty parameter $\varepsilon$. Moreover, two-level mesh method can save a large amount of CPU time. The error estimates in $H^1$ norm for velocity and in $L^2$ norm for pressure are derived. Finally, two numerical experiments are shown to support the efficiency of this new method.     

关于一类实二次域的类数与基本单位的上界

设素数P≡1(mod4),k,ε分别表示实二次域Q(p~(1/2))类数和基本单位．本文改进了类数h和基本单位ε的上界,证明了:hlogeε<1/4(p~(1/2) 6)log(2ep~(1/2)),并得到了几个重要的推论．     

A subdivision collocation method for solving two point boundary value problems of order three

In this paper, we propose a method for the numerical solution of self adjoint singularly perturbed third order boundary value problems in which the highest order derivative is multiplied by a small parameter $\varepsilon$. In this method, first we introduce the derivatives of two scale relations satisfied by the subdivision schemes. After that we use these derivatives to construct the subdivision collocation method for the numerical solution of singularly perturbed boundary value problems. Convergence of the subdivision collocation method is also discussed. Numerical examples are presented to illustrate the proposed method.     

Analysis of a fully discrete finite element method for the phase field model and approximation of its sharp interface limits

We propose and analyze a fully discrete finite element scheme for the phase field model describing the solidification process in materials science. The primary goal of this paper is to establish some useful a priori error estimates for the proposed numerical method, in particular, by focusing on the dependence of the error bounds on the parameter , known as the measure of the interface thickness. Optimal order error bounds are shown for the fully discrete scheme under some reasonable constraints on the mesh size and the time step size . In particular, it is shown that all error bounds depend on only in some lower polynomial order for small . The cruxes of the analysis are to establish stability estimates for the discrete solutions, to use a spectrum estimate result of Chen, and to establish a discrete counterpart of it for a linearized phase field operator to handle the nonlinear effect. Finally, as a nontrivial byproduct, the error estimates are used to establish convergence of the solution of the fully discrete scheme to solutions of the sharp interface limits of the phase field model under different scaling in its coefficients. The sharp interface limits include the classical Stefan problem, the generalized Stefan problems with surface tension and surface kinetics, the motion by mean curvature flow, and the Hele-Shaw model.

     

A Two-Grid Finite Element Approximation for Nonlinear Time Fractional Two-Term Mixed Sub-Diffusion and Diffusion Wave Equations

In this paper, we develop a two-grid method (TGM) based on the FEM for 2D nonlinear time fractional two-term mixed sub-diffusion and diffusion wave equations. A two-grid algorithm is proposed for solving the nonlinear system, which consists of two steps: a nonlinear FE system is solved on a coarse grid, then the linearized FE system is solved on the fine grid by Newton iteration based on the coarse solution. The fully discrete numerical approximation is analyzed, where the Galerkin finite element method for the space derivatives and the finite difference scheme for the time Caputo derivative with order $\alpha\in(1,2)$ and $\alpha_{1}\in(0,1)$. Numerical stability and optimal error estimate $O(h^{r+1}+H^{2r+2}+\tau^{\min\{3-\alpha,2-\alpha_{1}\}})$ in $L^{2}$-norm are presented for two-grid scheme, where $t,$ $H$ and $h$ are the time step size, coarse grid mesh size and fine grid mesh size, respectively. Finally, numerical experiments are provided to confirm our theoretical results and effectiveness of the proposed algorithm.     

CONVERGENCE OF AN ALTERNATING A-φ SCHEME FOR QUASI-MAGNETOSTATIC EDDY CURRENT PROBLEM

We propose in this paper an alternating A-$\phi$ method for the quasi-magnetostatic eddy current problem by means of finite element approximations. Bounds for continuous and discrete error in finite time are given. And it is verified that provided the time step $\tau$ is sufficiently small, the proposed algorithm yields for finite time $T$ an error of $O(h+\tau^{1/2})$ in the $L^2$-norm for the magnetic field $H(= \mu^{-1} \nabla \times A)$, where $h$ is the mesh size, $\mu$ the magnetic permeability.     

TWO-GRID CHARACTERISTIC FINITE VOLUME METHODS FOR NONLINEAR PARABOLIC PROBLEMS＊

In this work, two-grid characteristic finite volume schemes for the nonlinear parabolic problem are considered. In our algorithms, the diffusion term is discretized by the finite volume method, while the temporal differentiation and advection terms are treated by the characteristic scheme. Under some conditions about the coefficients and exact solution, optimal error estimates for the numerical solution are obtained. Furthermore, the two- grid characteristic finite volume methods involve solving a nonlinear equation on coarse mesh with mesh size H, a large linear problem for the Oseen two-grid characteristic finite volume method on a fine mesh with mesh size h = O（H2） or a large linear problem for the Newton two-grid characteristic finite volume method on a fine mesh with mesh size h = 0（I log hll/2H3）. These methods we studied provide the same convergence rate as that of the characteristic finite volume method, which involves solving one large nonlinear problem on a fine mesh with mesh size h. Some numerical results are presented to demonstrate the efficiency of the proposed methods.     

An Improved Two-Grid Technique for the Nonlinear Time-Fractional Parabolic Equation Based on the Block-Centered Finite Difference Method

A combined scheme of the improved two-grid technique with the block-centered finite difference method is constructed and analyzed to solve the nonlinear time-fractional parabolic equation. This method is considered where the nonlinear problem is solved only on a coarse grid of size $H$ and two linear problems based on the coarse-grid solutions and one Newton iteration is considered on a fine grid of size $h$. We provide the rigorous error estimate, which demonstrates that our scheme converges with order $\mathcal{O}(\Delta t^{2-\alpha}+h^2+H^4)$ on non-uniform rectangular grid. This result indicates that the improved two-grid method can obtain asymptotically optimal approximation as long as the mesh sizes satisfy $h=\mathcal{O}(H^2).$ Finally, numerical tests confirm the theoretical results of the presented method.     

On Numerical Solution of Quasilinear Boundary Value Problems with Two Small Parameters

We consider the singular perturbation problem $$-\varepsilon^2u"+\mu b(x,u)u'+c(x,u)=0,u(0),u(1)$$ given with two small parameters $\varepsilon$ and $\mu$ , $\mu =\varepsilon^{1+p},p>0$. The problem is solved numerically by using finite difference schemes on the mesh which is dense in the boundary layers. The convergence uniform in $\varepsilon$ is proved in the discrete $L^1$ norm. Some convergence results are given in the maximum norm as well.     

The Finite Element Analysis of the Controlled-Source Electromagnetic Induction Problems by Fractional-Step Projection Method

This paper provides a convergence analysis of a fractional-step projection method for the controlled-source electromagnetic induction problems in heterogenous electrically conduting media by means of finite element approximations. Error estimates in finite time are given. And it is verified that provided the time step $\tau$ is sufficiently small, the proposed algorithm yields for finite time $T$ an error of $\mathcal{O}(h^s+\tau)$) in the $L^2$-norm for the magnetic field $\boldsymbol{H},$ where $h$ is the mesh size and $1/2 < s≤1$.     

ERROR ANALYSIS OF A STABILIZED FINITE ELEMENT METHOD FOR THE GENERALIZED STOKES PROBLEM

This paper is devoted to the establishment of sharper $a$ $priori$stability and error estimates of a stabilized finite element method proposed by Barrenechea and Valentin for solving the generalized Stokes problem, which involves a viscosity $\nu$ and a reaction constant $\sigma$. With the establishment of sharper stability estimates and the help of $ad$ $hoc$finite element projections, we can explicitly establish the dependence of error bounds of velocity and pressure on the viscosity $\nu$, the reaction constant $\sigma$, and the mesh size $h$. Our analysis reveals that the viscosity $\nu$ and the reaction constant $\sigma$ respectively act in the numerator position and the denominator position in the error estimates of velocity and pressure in standard norms without any weights. Consequently, the stabilization method is indeed suitable for the generalized Stokes problem with a small viscosity $\nu$ and a large reaction constant $\sigma$. The sharper error estimates agree very well with the numerical results.     

Numerical Approximations of the Spectral Discretization of Flame Front Model

In this paper, we consider the numerical solution of the flame front equation, which is one of the most fundamental equations for modeling combustion theory. A schema combining a finite difference approach in the time direction and a spectral method for the space discretization is proposed. We give a detailed analysis for the proposed schema by providing some stability and error estimates in a particular case. For the general case, although we are unable to provide a rigorous proof for the stability, some numerical experiments are carried out to verify the efficiency of the schema. Our numerical results show that the stable solution manifolds have a simple structure when $\beta$ is small, while they become more complex as the bifurcation parameter $\beta$ increases. At last numerical experiments were performed to support the claim the solution of flame front equation preserves the same structure as K-S equation.     

Symmetry-breaking bifurcation of analytic solutions to free boundary problems: An application to a model of tumor growth

In this paper we develop a general technique for establishing analyticity of solutions of partial differential equations which depend on a parameter . The technique is worked out primarily for a free boundary problem describing a model of a stationary tumor. We prove the existence of infinitely many branches of symmetry-breaking solutions which bifurcate from any given radially symmetric steady state; these asymmetric solutions are analytic jointly in the spatial variables and in .

     

A Two-Grid Algorithm of Fully Discrete Galerkin Finite Element Methods for a Nonlinear Hyperbolic Equation

A two-grid finite element approximation is studied in the fully discrete scheme obtained by discretizing in both space and time for a nonlinear hyperbolic equation. The main idea of two-grid methods is to use a coarse-grid space ($S_H$) to produce a rough approximation for the solution of nonlinear hyperbolic problems and then use it as the initial guess on the fine-grid space ($S_h$). Error estimates are presented in $H^1$-norm, which show that two-grid methods can achieve the optimal convergence order as long as the two different girds satisfy $h$ = $\mathcal{O}$($H^2$). With the proposed techniques, we can obtain the same accuracy as standard finite element methods, and also save lots of time in calculation. Theoretical analyses and numerical examples are presented to confirm the methods.     

随机变延迟微分方程平衡方法的均方收敛性与稳定性

针对一类变延迟微分方程,应用全隐式方法—平衡方法,研究了其收敛性和稳定性.结果表明平衡方法以$\frac{1}{2}\gamma,\gamma\in（0,1]$阶收敛到精确解;并且强平衡方法和弱平衡方法都能保持解析解的均方稳定性;进一步数值实验验证了算法理论分析的正确性,并且表明全隐式的平衡方法比显式方法—Euler方法具有更好的稳定性.