This paper deals with superlinear fourth-order elliptic problem under Navier boundary condition. By using the mountain pass theorem and suitable truncation, a multiplicity result is established for all λ〉 0 and some...This paper deals with superlinear fourth-order elliptic problem under Navier boundary condition. By using the mountain pass theorem and suitable truncation, a multiplicity result is established for all λ〉 0 and some previous result is extended.展开更多
The nonconforming Crouzeix-Raviart type linear triangular finite element approximate to second-order elliptic problems is studied on anisotropic general triangular meshes in 2D satisfying the maximal angle condition a...The nonconforming Crouzeix-Raviart type linear triangular finite element approximate to second-order elliptic problems is studied on anisotropic general triangular meshes in 2D satisfying the maximal angle condition and the coordinate system condition. The optimal-order error estimates of the broken energy norm and L2-norm are obtained.展开更多
The author proves that there exist three solutions u(0), u(1) and u(2) in the following problem [GRAPHICS] where some conditions are imposed on Q and f. Here, 0 < u(0) < u(1), u(2) changes sign.
In this paper,we will study the nonlocal and nonvariational elliptic problem{−(1+a||u||_(q)^(αq))Δu=|u|^(p−1)u+h(x,u,∇_(u))inΩ,u=0 on∂Ω,(0.1)(1)where a>0,α>0,1<q<2^(∗),p∈(0,2^(∗)−1)∖{1}andΩis a boun...In this paper,we will study the nonlocal and nonvariational elliptic problem{−(1+a||u||_(q)^(αq))Δu=|u|^(p−1)u+h(x,u,∇_(u))inΩ,u=0 on∂Ω,(0.1)(1)where a>0,α>0,1<q<2^(∗),p∈(0,2^(∗)−1)∖{1}andΩis a bounded smooth domain in R^(N)(N≥2).Under suitable assumptions about h(x,u,∇u),we obtain\emph{a priori}estimates of positive solutions for the problem(0.1).Furthermore,we establish the existence of positive solutions by making use of these estimates and of the method of continuity.展开更多
A two-grid partition of unity method for second order elliptic problems is proposed and analyzed. The standard two-grid method is a local and parallel method usually leading to a discontinuous solution in the entire c...A two-grid partition of unity method for second order elliptic problems is proposed and analyzed. The standard two-grid method is a local and parallel method usually leading to a discontinuous solution in the entire computational domain. Partition of unity method is employed to glue all the local solutions together to get the global continuous one, which is optimal in HI-norm. Furthermore, it is shown that the L^2 error can be improved by using the coarse grid correction. Numerical experiments are reported to support the theoretical results.展开更多
The purpose of this paper is to solve nonselfadjoint elliptic problems with rapidly oscillatory coefficients. A two-order and two-scale approximate solution expression for nonselfadjoint elliptic problems is considere...The purpose of this paper is to solve nonselfadjoint elliptic problems with rapidly oscillatory coefficients. A two-order and two-scale approximate solution expression for nonselfadjoint elliptic problems is considered, and the error estimation of the twoorder and two-scale approximate solution is derived. The numerical result shows that the presented approximation solution is effective.展开更多
This paper is concerned with convergence and superconvergence properties of the local discontinuous Galerkin(LDG)method for two-dimensional semilinear second-order elliptic problems of the form−Δu=f(x,y,u)on Cartesia...This paper is concerned with convergence and superconvergence properties of the local discontinuous Galerkin(LDG)method for two-dimensional semilinear second-order elliptic problems of the form−Δu=f(x,y,u)on Cartesian grids.By introducing special GaussRadau projections and using duality arguments,we obtain,under some suitable choice of numerical fuxes,the optimal convergence order in L2-norm of O(h^(p+1))for the LDG solution and its gradient,when tensor product polynomials of degree at most p and grid size h are used.Moreover,we prove that the LDG solutions are superconvergent with an order p+2 toward particular Gauss-Radau projections of the exact solutions.Finally,we show that the error between the gradient of the LDG solution and the gradient of a special Gauss-Radau projection of the exact solution achieves(p+1)-th order superconvergence.Some numerical experiments are performed to illustrate the theoretical results.展开更多
This article extends a recently developed superconvergence result for weak Galerkin(WG)approximations for modeling partial differential equations from constant coefficients to variable coefficients.This superconvergen...This article extends a recently developed superconvergence result for weak Galerkin(WG)approximations for modeling partial differential equations from constant coefficients to variable coefficients.This superconvergence features a rate that is two orders higher than the optimal-order error estimates in the usual energy and L^(2)norms.The extension from constant to variable coefficients for the modeling equations is highly non-trivial.The underlying technical analysis is based on a sequence of projections and decompositions.Numerical results confirm the superconvergence theory for second-order elliptic problems with variable coefficients.展开更多
This paper considers the regularity of solutions to mixed boundary value problems in small-angle regions for elliptic equations. By constructing a specific barrier function, we proved that under the assumption of suff...This paper considers the regularity of solutions to mixed boundary value problems in small-angle regions for elliptic equations. By constructing a specific barrier function, we proved that under the assumption of sufficient regularity of boundary conditions and coefficients, as long as the angle is sufficiently small, the regularity of the solution to the mixed boundary value problem of the second-order elliptic equation can reach any order.展开更多
In this paper,we study the following critical elliptic problem with a variable exponent:{-Δu=u^(p+ϵa(x))inΩ,u>0 inΩ,u=0 on∂Ω,where a(x)2∈C^(2)(Ω),p=N+2/N-2,∈>0,andΩis a smooth bounded domain in R^(N)(N&g...In this paper,we study the following critical elliptic problem with a variable exponent:{-Δu=u^(p+ϵa(x))inΩ,u>0 inΩ,u=0 on∂Ω,where a(x)2∈C^(2)(Ω),p=N+2/N-2,∈>0,andΩis a smooth bounded domain in R^(N)(N>4).We show that for∈small enough,there exists a family of bubble solutions concentrating at the negative stable critical point of the function a(x).This is a new perturbation to the critical elliptic equation in contrast to the usual subcritical or supercritical perturbation,and gives the first existence result for the critical elliptic problem with a variable exponent.展开更多
The multigrid V-cycle methods for adaptive finite element discretizations of two-dimensional elliptic problems with discontinuous coefficients are considered.Under the conditions that the coefficient is quasi-monotone...The multigrid V-cycle methods for adaptive finite element discretizations of two-dimensional elliptic problems with discontinuous coefficients are considered.Under the conditions that the coefficient is quasi-monotone up to a constant and the meshes are locally refined by using the newest vertex bisection algorithm,some uniform convergence results are proved for the standard multigrid V-cycle algorithm with Gauss-Seidel relaxations performed only on new nodes and their immediate neighbours.The multigrid V-cycle algorithm uses O(N)operations per iteration and is optimal.展开更多
In this paper, we obtain some existence results for a class of singular semilinear elliptic problems where we improve some earlier results of Zhijun Zhang. We show the existence of entire positive solutions without th...In this paper, we obtain some existence results for a class of singular semilinear elliptic problems where we improve some earlier results of Zhijun Zhang. We show the existence of entire positive solutions without the monotonic condition imposed in Zhang’s paper. The main point of our technique is to choose an approximating sequence and prove its convergence. The desired compactness can be obtained by the Sobolev embedding theorems.展开更多
In this paper, we discuss the finite volume element method of P1-nonconforming quadrilateral element for elliptic problems and obtain optimal error estimates for general quadrilateral partition. An optimal cascadic mu...In this paper, we discuss the finite volume element method of P1-nonconforming quadrilateral element for elliptic problems and obtain optimal error estimates for general quadrilateral partition. An optimal cascadic multigrid algorithm is proposed to solve the non-symmetric large-scale system resulting from such discretization. Numerical experiments are reported to support our theoretical results.展开更多
We give an existence result of entropy and renormalized solutions for strongly nonlinear elliptic equations in the framework of Sobolev spaces with variable exponents of the type: -div (a(x, u,▽u)+φ(u))+g(...We give an existence result of entropy and renormalized solutions for strongly nonlinear elliptic equations in the framework of Sobolev spaces with variable exponents of the type: -div (a(x, u,▽u)+φ(u))+g(x, u,▽u)=μ, where the right-hand side belongs to L^1(Ω)+W^-1,p'(x)(Ω), -div(a(x, u,▽u)) is a Leray-Lions operator defined from W^-1,p'(x)(Ω) into its dual and φ∈C^0(R,R^N). The function g(x, u,▽u) is a non linear lower order term with natural growth with respect to |▽u| satisfying the sign condition, that is, g(x, u,▽u)u ≥ 0.展开更多
In this paper, the method of non-conforming mixed finite element for second order elliptic problems is discussed and a format of real optimal order for the lowest order error estimate.
The maximum principle is a basic qualitative property of the solution of second-order elliptic boundary value problems.The preservation of the qualitative characteristics,such as the maximum principle,in discrete mode...The maximum principle is a basic qualitative property of the solution of second-order elliptic boundary value problems.The preservation of the qualitative characteristics,such as the maximum principle,in discrete model is one of the key requirements.It is well known that standard linear finite element solution does not satisfy maximum principle on general triangular meshes in 2D.In this paper we consider how to enforce discrete maximum principle for linear finite element solutions for the linear second-order self-adjoint elliptic equation.First approach is based on repair technique,which is a posteriori correction of the discrete solution.Second method is based on constrained optimization.Numerical tests that include anisotropic cases demonstrate how our method works for problems for which the standard finite element methods produce numerical solutions that violate the discrete maximum principle.展开更多
Suppose Ω belong to R^N(N≥3) is a smooth bounded domain,ξi∈Ω,0〈ai〈√μ,μ:=((N-1)/2)^2,0≤μi〈(√μ-ai)^2,ai〈bi〈ai+1 and pi:=2N/N-2(1+ai-bi)are the weighted critical Hardy-Sobolev exponents, i ...Suppose Ω belong to R^N(N≥3) is a smooth bounded domain,ξi∈Ω,0〈ai〈√μ,μ:=((N-1)/2)^2,0≤μi〈(√μ-ai)^2,ai〈bi〈ai+1 and pi:=2N/N-2(1+ai-bi)are the weighted critical Hardy-Sobolev exponents, i = 1, 2,..., k, k ≥ 2. We deal with the conditions that ensure the existence of positive solutions to the multi-singular and multi-critical elliptic problem ∑i=1^k(-div(|x-ξi|^-2ai△↓u)-μiu/|x-ξi|^2(1+ai)-u^pi-1/|x-ξi|^bipi)=0with Dirichlet boundary condition, which involves the weighted Hardy inequality and the weighted Hardy-Sobolev inequality. The results depend crucially on the parameters ai, bi and #i, i -- 1, 2,..., k.展开更多
A second order isoparametric finite element method (IPFEM) is proposed for elliptic interface problems. It yields better accuracy than some existing second-order methods, when the coefficients or the flux across the...A second order isoparametric finite element method (IPFEM) is proposed for elliptic interface problems. It yields better accuracy than some existing second-order methods, when the coefficients or the flux across the immersed curved interface is discontinuous. Based on an initial Cartesian mesh, a mesh optimization strategy is presented by employing curved boundary elements at the interface, and an incomplete quadratic finite element space is constructed on the optimized mesh. It turns out that the number of curved boundary elements is far less than that of the straight one, and the total degree of freedom is almost the same as the uniform Cartesian mesh. Numerical examples with simple and complicated geometrical interfaces demonstrate the efficiency of the proposed method.展开更多
This paper is devoted to study of an iterative procedure for domain decomposition method of second order elliptic problem with mixed boundary conditions (i.e., Dirichlet condition on a part of boundary and Neumann con...This paper is devoted to study of an iterative procedure for domain decomposition method of second order elliptic problem with mixed boundary conditions (i.e., Dirichlet condition on a part of boundary and Neumann condition on the another part of boundary). For the pure Dirichlet problem, Marini and Quarteroni [3], [4] considered a similar approach, which is extended to more complex problem in this paper.展开更多
基金The 985 Program of Jilin Universitythe Science Research Foundation for Excellent Young Teachers of College of Mathematics at Jilin University
文摘This paper deals with superlinear fourth-order elliptic problem under Navier boundary condition. By using the mountain pass theorem and suitable truncation, a multiplicity result is established for all λ〉 0 and some previous result is extended.
基金supported by the National Natural Science Foundation of China (No. 10971203)
文摘The nonconforming Crouzeix-Raviart type linear triangular finite element approximate to second-order elliptic problems is studied on anisotropic general triangular meshes in 2D satisfying the maximal angle condition and the coordinate system condition. The optimal-order error estimates of the broken energy norm and L2-norm are obtained.
文摘The author proves that there exist three solutions u(0), u(1) and u(2) in the following problem [GRAPHICS] where some conditions are imposed on Q and f. Here, 0 < u(0) < u(1), u(2) changes sign.
基金supported by National Natural Science Foundation of China(11801167)Hunan Provincial Natural Science Foundation of China(2019JJ50387).
文摘In this paper,we will study the nonlocal and nonvariational elliptic problem{−(1+a||u||_(q)^(αq))Δu=|u|^(p−1)u+h(x,u,∇_(u))inΩ,u=0 on∂Ω,(0.1)(1)where a>0,α>0,1<q<2^(∗),p∈(0,2^(∗)−1)∖{1}andΩis a bounded smooth domain in R^(N)(N≥2).Under suitable assumptions about h(x,u,∇u),we obtain\emph{a priori}estimates of positive solutions for the problem(0.1).Furthermore,we establish the existence of positive solutions by making use of these estimates and of the method of continuity.
基金Project supported by the National Natural Science Foundation of China(No.40074031)the Science Foundation of the Science and Technology Commission of Shanghai Municipalitythe Program for Young Excellent Talents in Tongji University(No.2007kj008)
文摘A two-grid partition of unity method for second order elliptic problems is proposed and analyzed. The standard two-grid method is a local and parallel method usually leading to a discontinuous solution in the entire computational domain. Partition of unity method is employed to glue all the local solutions together to get the global continuous one, which is optimal in HI-norm. Furthermore, it is shown that the L^2 error can be improved by using the coarse grid correction. Numerical experiments are reported to support the theoretical results.
基金Project supported by the National Natural Science Foundation of China(No.10590353)the Science Research Project of National University of Defense Technology(No.JC09-02-05)
文摘The purpose of this paper is to solve nonselfadjoint elliptic problems with rapidly oscillatory coefficients. A two-order and two-scale approximate solution expression for nonselfadjoint elliptic problems is considered, and the error estimation of the twoorder and two-scale approximate solution is derived. The numerical result shows that the presented approximation solution is effective.
基金This research was supported by the NASA Nebraska Space Grant(Federal Grant/Award Number 80NSSC20M0112).
文摘This paper is concerned with convergence and superconvergence properties of the local discontinuous Galerkin(LDG)method for two-dimensional semilinear second-order elliptic problems of the form−Δu=f(x,y,u)on Cartesian grids.By introducing special GaussRadau projections and using duality arguments,we obtain,under some suitable choice of numerical fuxes,the optimal convergence order in L2-norm of O(h^(p+1))for the LDG solution and its gradient,when tensor product polynomials of degree at most p and grid size h are used.Moreover,we prove that the LDG solutions are superconvergent with an order p+2 toward particular Gauss-Radau projections of the exact solutions.Finally,we show that the error between the gradient of the LDG solution and the gradient of a special Gauss-Radau projection of the exact solution achieves(p+1)-th order superconvergence.Some numerical experiments are performed to illustrate the theoretical results.
基金supported by U.S.National Science Foundation IR/D program while working at U.S.National Science Foundationsupported by U.S.National Science Foundation(Grant No.DMS-1620016)+1 种基金supported by Zhejiang Provincial Natural Science Foundation of China(Grant No.LY23A010005)National Natural Science Foundation of China(Grant No.12071184)。
文摘This article extends a recently developed superconvergence result for weak Galerkin(WG)approximations for modeling partial differential equations from constant coefficients to variable coefficients.This superconvergence features a rate that is two orders higher than the optimal-order error estimates in the usual energy and L^(2)norms.The extension from constant to variable coefficients for the modeling equations is highly non-trivial.The underlying technical analysis is based on a sequence of projections and decompositions.Numerical results confirm the superconvergence theory for second-order elliptic problems with variable coefficients.
文摘This paper considers the regularity of solutions to mixed boundary value problems in small-angle regions for elliptic equations. By constructing a specific barrier function, we proved that under the assumption of sufficient regularity of boundary conditions and coefficients, as long as the angle is sufficiently small, the regularity of the solution to the mixed boundary value problem of the second-order elliptic equation can reach any order.
基金supported by National Natural Science Foundation of China(Grant No.11971147)supported by National Natural Science Foundation of China(Grant Nos.11831009 and 12171183)the Fundamental Research Funds for the Central Universities(Grant No.KJ02072020-0319).
文摘In this paper,we study the following critical elliptic problem with a variable exponent:{-Δu=u^(p+ϵa(x))inΩ,u>0 inΩ,u=0 on∂Ω,where a(x)2∈C^(2)(Ω),p=N+2/N-2,∈>0,andΩis a smooth bounded domain in R^(N)(N>4).We show that for∈small enough,there exists a family of bubble solutions concentrating at the negative stable critical point of the function a(x).This is a new perturbation to the critical elliptic equation in contrast to the usual subcritical or supercritical perturbation,and gives the first existence result for the critical elliptic problem with a variable exponent.
基金supported by the NSF of China (Grant Nos.12171238,12261160361)supported in part by the China NSF for Distinguished Young Scholars (Grant No.11725106)by the China NSF major project (Grant No.11831016).
文摘The multigrid V-cycle methods for adaptive finite element discretizations of two-dimensional elliptic problems with discontinuous coefficients are considered.Under the conditions that the coefficient is quasi-monotone up to a constant and the meshes are locally refined by using the newest vertex bisection algorithm,some uniform convergence results are proved for the standard multigrid V-cycle algorithm with Gauss-Seidel relaxations performed only on new nodes and their immediate neighbours.The multigrid V-cycle algorithm uses O(N)operations per iteration and is optimal.
基金supported in part by NSF(Youth) of Shandong Province and NNSF of China
文摘In this paper, we obtain some existence results for a class of singular semilinear elliptic problems where we improve some earlier results of Zhijun Zhang. We show the existence of entire positive solutions without the monotonic condition imposed in Zhang’s paper. The main point of our technique is to choose an approximating sequence and prove its convergence. The desired compactness can be obtained by the Sobolev embedding theorems.
文摘In this paper, we discuss the finite volume element method of P1-nonconforming quadrilateral element for elliptic problems and obtain optimal error estimates for general quadrilateral partition. An optimal cascadic multigrid algorithm is proposed to solve the non-symmetric large-scale system resulting from such discretization. Numerical experiments are reported to support our theoretical results.
文摘We give an existence result of entropy and renormalized solutions for strongly nonlinear elliptic equations in the framework of Sobolev spaces with variable exponents of the type: -div (a(x, u,▽u)+φ(u))+g(x, u,▽u)=μ, where the right-hand side belongs to L^1(Ω)+W^-1,p'(x)(Ω), -div(a(x, u,▽u)) is a Leray-Lions operator defined from W^-1,p'(x)(Ω) into its dual and φ∈C^0(R,R^N). The function g(x, u,▽u) is a non linear lower order term with natural growth with respect to |▽u| satisfying the sign condition, that is, g(x, u,▽u)u ≥ 0.
基金Project supported by the Cultivating Foundation of Youthful Backbone of Science and Technologyof Beijing, the National Science
文摘In this paper, the method of non-conforming mixed finite element for second order elliptic problems is discussed and a format of real optimal order for the lowest order error estimate.
基金the National Nuclear Security Administration of the U.S.Department of Energy at Los Alamos National Laboratory under Contract No.DE-AC52-06NA25396the DOE Office of Science Advanced Scientific Computing Research(ASCR)Program in Applied Mathematics Research.The first author has been supported in part by the Czech Ministry of Education projects MSM 6840770022 and LC06052(Necas Center for Mathematical Modeling).
文摘The maximum principle is a basic qualitative property of the solution of second-order elliptic boundary value problems.The preservation of the qualitative characteristics,such as the maximum principle,in discrete model is one of the key requirements.It is well known that standard linear finite element solution does not satisfy maximum principle on general triangular meshes in 2D.In this paper we consider how to enforce discrete maximum principle for linear finite element solutions for the linear second-order self-adjoint elliptic equation.First approach is based on repair technique,which is a posteriori correction of the discrete solution.Second method is based on constrained optimization.Numerical tests that include anisotropic cases demonstrate how our method works for problems for which the standard finite element methods produce numerical solutions that violate the discrete maximum principle.
基金supported partly by the National Natural Science Foundation of China (10771219)the Science Foundation of the SEAC of China (07ZN03)
文摘Suppose Ω belong to R^N(N≥3) is a smooth bounded domain,ξi∈Ω,0〈ai〈√μ,μ:=((N-1)/2)^2,0≤μi〈(√μ-ai)^2,ai〈bi〈ai+1 and pi:=2N/N-2(1+ai-bi)are the weighted critical Hardy-Sobolev exponents, i = 1, 2,..., k, k ≥ 2. We deal with the conditions that ensure the existence of positive solutions to the multi-singular and multi-critical elliptic problem ∑i=1^k(-div(|x-ξi|^-2ai△↓u)-μiu/|x-ξi|^2(1+ai)-u^pi-1/|x-ξi|^bipi)=0with Dirichlet boundary condition, which involves the weighted Hardy inequality and the weighted Hardy-Sobolev inequality. The results depend crucially on the parameters ai, bi and #i, i -- 1, 2,..., k.
基金Supported by the National Natural Science Foundation of China (11071216 and 11101361)
文摘A second order isoparametric finite element method (IPFEM) is proposed for elliptic interface problems. It yields better accuracy than some existing second-order methods, when the coefficients or the flux across the immersed curved interface is discontinuous. Based on an initial Cartesian mesh, a mesh optimization strategy is presented by employing curved boundary elements at the interface, and an incomplete quadratic finite element space is constructed on the optimized mesh. It turns out that the number of curved boundary elements is far less than that of the straight one, and the total degree of freedom is almost the same as the uniform Cartesian mesh. Numerical examples with simple and complicated geometrical interfaces demonstrate the efficiency of the proposed method.
文摘This paper is devoted to study of an iterative procedure for domain decomposition method of second order elliptic problem with mixed boundary conditions (i.e., Dirichlet condition on a part of boundary and Neumann condition on the another part of boundary). For the pure Dirichlet problem, Marini and Quarteroni [3], [4] considered a similar approach, which is extended to more complex problem in this paper.