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 superconvergence in the finite element method is a phenomenon in which the fi-nite element approximation converges to the exact solution at a rate higher than the optimal order error estimate. Wang proposed and an...The superconvergence in the finite element method is a phenomenon in which the fi-nite element approximation converges to the exact solution at a rate higher than the optimal order error estimate. Wang proposed and analyzed superconvergence of the conforming finite element method by L2-projections. However, since the conforming finite element method (CFEM) requires a strong continuity, it is not easy to construct such finite elements for the complex partial differential equations. Thus, the nonconforming finite element method (NCFEM) is more appealing computationally due to better stability and flexibility properties compared to CFEM. The objective of this paper is to establish a general superconvergence result for the nonconforming finite element approximations for second-order elliptic problems by L2-projection methods by applying the idea presented in Wang. MATLAB codes are published at https://github.com/annaleeharris/Superconvergence-NCFEM for anyone to use and to study. The results of numerical experiments show great promise for the robustness, reliability, flexibility and accuracy of superconvergence in NCFEM by L2- projections.展开更多
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.展开更多
Let n≥2 and let L be a second-order elliptic operator of divergence form with coefficients consisting of both an elliptic symmetric part and a BMO anti-symmetric part in ℝ^(n).In this article,we consider the weighted...Let n≥2 and let L be a second-order elliptic operator of divergence form with coefficients consisting of both an elliptic symmetric part and a BMO anti-symmetric part in ℝ^(n).In this article,we consider the weighted Kato square root problem for L.More precisely,we prove that the square root L^(1/2)satisfies the weighted L^(p)estimates||L^(1/2)(f)||L_(ω)^p(R^(n))≤C||■f||L_(ω)^p(R^(n);R^(n))for any p∈(1,∞)andω∈Ap(ℝ^(n))(the class of Muckenhoupt weights),and that||■f||L_(ω)^p(R^(n);R^(n))≤C||L^(1/2)(f)||L_(ω)^p(R^(n))for any p∈(1,2+ε)andω∈Ap(ℝ^(n))∩RH_(2+ε/p),(R^(n))(the class of reverse Hölder weights),whereε∈(0,∞)is a constant depending only on n and the operator L,and where(2+ε/p)'denotes the Hölder conjugate exponent of 2+ε/p.Moreover,for any given q∈(2,∞),we give a sufficient condition to obtain that||■f||L_(ω)^p(R^(n);R^(n))≤C||L^(1/2)(f)||L_(ω)^p(R^(n))for any p∈(1,q)andω∈A_(p)(R^(n))∩pRH_(q/p),(R^(n)).As an application,we prove that when the coefficient matrix A that appears in L satisfies the small BMO condition,the Riesz transform∇L^(−1/2)is bounded on L_(ω)^(p)(ℝ^(n))for any given p∈(1,∞)andω∈Ap(ℝ^(n)).Furthermore,applications to the weighted L^(2)-regularity problem with the Dirichlet or the Neumann boundary condition are also given.展开更多
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.展开更多
This paper provides a proof for the uniform convergence rate (independently of the number of mesh levels) for the nonnested V-cycle multigrid method for nonsymmetric and indefinite second-order elliptic problems.
A class of nonlinear boundary value problems(BVP) for the second_order E 2 class elliptic systems in general form is discussed. By introducing a kind of transformation,this kind of BVP is reduced to a class of genera...A class of nonlinear boundary value problems(BVP) for the second_order E 2 class elliptic systems in general form is discussed. By introducing a kind of transformation,this kind of BVP is reduced to a class of generalized nonlinear Riemann_Hilbert BVP. And then some singular integral operators are introduced to establish the equivalent nonlinear singular integral equations. The solvability is proved under some suitable hypotheses by means of the properties of singular integral operators and the function theoretic methods.展开更多
In this paper the least-squares mixed finite element is considered for solving second-order elliptic problems in two dimensional domains. The primary solution u and the flux σ are approximated using finite element sp...In this paper the least-squares mixed finite element is considered for solving second-order elliptic problems in two dimensional domains. The primary solution u and the flux σ are approximated using finite element spaces consisting of piecewise polynomials of degree k and r respectively. Based on interpolation operators and an auxiliary projection, superconvergent H1-error estimates of both the primary solution approximation uh and the flux approximation σh are obtained under the standard quasi-uniform assumption on finite element partition. The superconvergence indicates an accuracy of O(hr+2) for the least-squares mixed finite element approximation if Raviart-Thomas or Brezzi-Douglas-Fortin-Marini elements of order r are employed with optimal error estimate of O(hr+1).展开更多
Diagonalized Chebyshev rational spectral methods for solving second-order elliptic problems on the half/whole line are proposed.Some Sobolev bi-orthogonal rational basis functions are constructed which lead to the dia...Diagonalized Chebyshev rational spectral methods for solving second-order elliptic problems on the half/whole line are proposed.Some Sobolev bi-orthogonal rational basis functions are constructed which lead to the diagonalization of discrete systems.Accordingly,both the exact solutions and the approximate solutions can be represented as infinite and truncated Fourier-like Chebyshev rational series.Numerical results demonstrate the effectiveness of the suggested approaches.展开更多
In this paper,a new hybridized mixed formulation of weak Galerkin method is studied for a second order elliptic problem.This method is designed by approximate some operators with discontinuous piecewise polynomials in...In this paper,a new hybridized mixed formulation of weak Galerkin method is studied for a second order elliptic problem.This method is designed by approximate some operators with discontinuous piecewise polynomials in a shape regular finite element partition.Some discrete inequalities are presented on discontinuous spaces and optimal order error estimations are established.Some numerical results are reported to show super convergence and confirm the theory of the mixed weak Galerkin method.展开更多
The singularly perturbed elliptic equation boundary value problem with turning point is considered. Using the method of multiple scales and the comparison theorem, the asymptotic behavior of solution for the boundary ...The singularly perturbed elliptic equation boundary value problem with turning point is considered. Using the method of multiple scales and the comparison theorem, the asymptotic behavior of solution for the boundary value problem is studied.展开更多
A class of nonlocal boundary value probl em s for elliptic systems in the unbounded domains are considered. Under suitable c onditions, the existence of solution and the comparison theorem for the boundary value prob...A class of nonlocal boundary value probl em s for elliptic systems in the unbounded domains are considered. Under suitable c onditions, the existence of solution and the comparison theorem for the boundary value problems are studied.展开更多
The geometric characteristics of fractures within a rock mass can be inferred by the data sampling from boreholes or exposed surfaces.Recently,the universal elliptical disc(UED)model was developed to represent natural...The geometric characteristics of fractures within a rock mass can be inferred by the data sampling from boreholes or exposed surfaces.Recently,the universal elliptical disc(UED)model was developed to represent natural fractures,where the fracture is assumed to be an elliptical disc and the fracture orientation,rotation angle,length of the long axis and ratio of short-long axis lengths are considered as variables.This paper aims to estimate the fracture size-and azimuth-related parameters in the UED model based on the trace information from sampling windows.The stereological relationship between the trace length,size-and azimuth-related parameters of the UED model was established,and the formulae of the mean value and standard deviation of trace length were proposed.The proposed formulae were validated via the Monte Carlo simulations with less than 5%of error rate between the calculated and true values.With respect to the estimation of the size-and azimuth-related parameters using the trace length,an optimization method was developed based on the pre-assumed size and azimuth distribution forms.A hypothetical case study was designed to illustrate and verify the parameter estimation method,where three combinations of the sampling windows were used to estimate the parameters,and the results showed that the estimated values could agree well with the true values.Furthermore,a hypothetical three-dimensional(3D)elliptical fracture network was constructed,and the circular disc,non-UED and UED models were used to represent it.The simulated trace information from different models was compared,and the results clearly illustrated the superiority of the proposed UED model over the existing circular disc and non-UED models。展开更多
The nonlocal boundary value problems for nonlinear elliptic systems in the unbounded domain are considered. Under suitable conditions the existence of solution and comparison theorem for the boundary value problems ar...The nonlocal boundary value problems for nonlinear elliptic systems in the unbounded domain are considered. Under suitable conditions the existence of solution and comparison theorem for the boundary value problems are studied.展开更多
Firstly, the Riemann boundary value problem for a kind of degenerate elliptic sys- tem of the first order equations in R4 is proposed. Then, with the help of the one-to-one correspondence between the theory of Cliffor...Firstly, the Riemann boundary value problem for a kind of degenerate elliptic sys- tem of the first order equations in R4 is proposed. Then, with the help of the one-to-one correspondence between the theory of Clifford valued generalized regular functions and that of the degenerate elliptic system's solution, the boundary value problem as stated above is trans- formed into a boundary value problem related to the generalized regular functions in Clifford analysis. Moreover, the solution of the Riemann boundary value problem for the degenerate elliptic system is explicitly described by using a kind of singular integral operator. Finally, the conditions for the existence of solutions of the oblique derivative problem for another kind of degenerate elliptic system of the first order equations in R4 are derived.展开更多
In this paper a singular perturbation of boundary value problem for elliptic partial differential equations of higher order is considered by using the differential inequalities. The uniformly valid asymptotic expansio...In this paper a singular perturbation of boundary value problem for elliptic partial differential equations of higher order is considered by using the differential inequalities. The uniformly valid asymptotic expansion in entire region is obtained.展开更多
Some superconvergence results of generalized difference solution for elliptic boundary value problem are given. It is shown that optimal points of the stresses for generalized difference method are the same as that fo...Some superconvergence results of generalized difference solution for elliptic boundary value problem are given. It is shown that optimal points of the stresses for generalized difference method are the same as that for finite element method.展开更多
This paper deals with the Neumann problem for a class of semilinear elliptic equations -△u + u =|u|2*-2u+ μ|u|q-2u in Ω, au/ar= |u|(?)*-2u on aΩ, where 2 = 2N/N-2, s=2(N-1)/N-2, 1 <q<2,N(?)3,μ>γ denotes...This paper deals with the Neumann problem for a class of semilinear elliptic equations -△u + u =|u|2*-2u+ μ|u|q-2u in Ω, au/ar= |u|(?)*-2u on aΩ, where 2 = 2N/N-2, s=2(N-1)/N-2, 1 <q<2,N(?)3,μ>γ denotes the unit outward normal to boundary aΩ. By vaxiational method and dual fountain theorem, the existence of infinitely many solutions with negative energy is proved.展开更多
In this paper, the first boundary value problem for quasilinear equation of the form △A(u,x)+∑i=1^m δb^i(u,x)/δxi+c(u,x)=0,Au(u,x) ≥0is studied. By using the compensated compactness theory, some results...In this paper, the first boundary value problem for quasilinear equation of the form △A(u,x)+∑i=1^m δb^i(u,x)/δxi+c(u,x)=0,Au(u,x) ≥0is studied. By using the compensated compactness theory, some results on the existence of weak solution are established. In addition, under certain condition the uniqueness of solution is proved.展开更多
In this paper, authors discuss the numerical methods of general discontinuous boundary value problems for elliptic complex equations of first order, They first give the well posedness of general discontinuous boundary...In this paper, authors discuss the numerical methods of general discontinuous boundary value problems for elliptic complex equations of first order, They first give the well posedness of general discontinuous boundary value problems, reduce the discontinuous boundary value problems to a variation problem, and then find the numerical solutions of above problem by the finite element method. Finally authors give some error-estimates of the foregoing numerical solutions.展开更多
基金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 superconvergence in the finite element method is a phenomenon in which the fi-nite element approximation converges to the exact solution at a rate higher than the optimal order error estimate. Wang proposed and analyzed superconvergence of the conforming finite element method by L2-projections. However, since the conforming finite element method (CFEM) requires a strong continuity, it is not easy to construct such finite elements for the complex partial differential equations. Thus, the nonconforming finite element method (NCFEM) is more appealing computationally due to better stability and flexibility properties compared to CFEM. The objective of this paper is to establish a general superconvergence result for the nonconforming finite element approximations for second-order elliptic problems by L2-projection methods by applying the idea presented in Wang. MATLAB codes are published at https://github.com/annaleeharris/Superconvergence-NCFEM for anyone to use and to study. The results of numerical experiments show great promise for the robustness, reliability, flexibility and accuracy of superconvergence in NCFEM by L2- projections.
文摘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 the Key Project of Gansu Provincial National Science Foundation(23JRRA1022)the National Natural Science Foundation of China(12071431)+1 种基金the Fundamental Research Funds for the Central Universities(lzujbky-2021-ey18)the Innovative Groups of Basic Research in Gansu Province(22JR5RA391).
文摘Let n≥2 and let L be a second-order elliptic operator of divergence form with coefficients consisting of both an elliptic symmetric part and a BMO anti-symmetric part in ℝ^(n).In this article,we consider the weighted Kato square root problem for L.More precisely,we prove that the square root L^(1/2)satisfies the weighted L^(p)estimates||L^(1/2)(f)||L_(ω)^p(R^(n))≤C||■f||L_(ω)^p(R^(n);R^(n))for any p∈(1,∞)andω∈Ap(ℝ^(n))(the class of Muckenhoupt weights),and that||■f||L_(ω)^p(R^(n);R^(n))≤C||L^(1/2)(f)||L_(ω)^p(R^(n))for any p∈(1,2+ε)andω∈Ap(ℝ^(n))∩RH_(2+ε/p),(R^(n))(the class of reverse Hölder weights),whereε∈(0,∞)is a constant depending only on n and the operator L,and where(2+ε/p)'denotes the Hölder conjugate exponent of 2+ε/p.Moreover,for any given q∈(2,∞),we give a sufficient condition to obtain that||■f||L_(ω)^p(R^(n);R^(n))≤C||L^(1/2)(f)||L_(ω)^p(R^(n))for any p∈(1,q)andω∈A_(p)(R^(n))∩pRH_(q/p),(R^(n)).As an application,we prove that when the coefficient matrix A that appears in L satisfies the small BMO condition,the Riesz transform∇L^(−1/2)is bounded on L_(ω)^(p)(ℝ^(n))for any given p∈(1,∞)andω∈Ap(ℝ^(n)).Furthermore,applications to the weighted L^(2)-regularity problem with the Dirichlet or the Neumann boundary condition are also given.
基金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.
文摘This paper provides a proof for the uniform convergence rate (independently of the number of mesh levels) for the nonnested V-cycle multigrid method for nonsymmetric and indefinite second-order elliptic problems.
文摘A class of nonlinear boundary value problems(BVP) for the second_order E 2 class elliptic systems in general form is discussed. By introducing a kind of transformation,this kind of BVP is reduced to a class of generalized nonlinear Riemann_Hilbert BVP. And then some singular integral operators are introduced to establish the equivalent nonlinear singular integral equations. The solvability is proved under some suitable hypotheses by means of the properties of singular integral operators and the function theoretic methods.
基金Supported by National Science Foundation of China and the Foundation of China State Education Commission and the Special Funds for Major State Basic Research Projects.
文摘In this paper the least-squares mixed finite element is considered for solving second-order elliptic problems in two dimensional domains. The primary solution u and the flux σ are approximated using finite element spaces consisting of piecewise polynomials of degree k and r respectively. Based on interpolation operators and an auxiliary projection, superconvergent H1-error estimates of both the primary solution approximation uh and the flux approximation σh are obtained under the standard quasi-uniform assumption on finite element partition. The superconvergence indicates an accuracy of O(hr+2) for the least-squares mixed finite element approximation if Raviart-Thomas or Brezzi-Douglas-Fortin-Marini elements of order r are employed with optimal error estimate of O(hr+1).
基金This work was supported in part by National Natural Science Foun-dation of China(Nos.11571238 and 11601332).
文摘Diagonalized Chebyshev rational spectral methods for solving second-order elliptic problems on the half/whole line are proposed.Some Sobolev bi-orthogonal rational basis functions are constructed which lead to the diagonalization of discrete systems.Accordingly,both the exact solutions and the approximate solutions can be represented as infinite and truncated Fourier-like Chebyshev rational series.Numerical results demonstrate the effectiveness of the suggested approaches.
文摘In this paper,a new hybridized mixed formulation of weak Galerkin method is studied for a second order elliptic problem.This method is designed by approximate some operators with discontinuous piecewise polynomials in a shape regular finite element partition.Some discrete inequalities are presented on discontinuous spaces and optimal order error estimations are established.Some numerical results are reported to show super convergence and confirm the theory of the mixed weak Galerkin method.
文摘The singularly perturbed elliptic equation boundary value problem with turning point is considered. Using the method of multiple scales and the comparison theorem, the asymptotic behavior of solution for the boundary value problem is studied.
文摘A class of nonlocal boundary value probl em s for elliptic systems in the unbounded domains are considered. Under suitable c onditions, the existence of solution and the comparison theorem for the boundary value problems are studied.
基金funded by National Natural Science Foundation of China(Grant No.41972264)Zhejiang Provincial Natural Science Foundation of China(Grant No.LR22E080002)the Observation and Research Station of Geohazards in Zhejiang,Ministry of Natural Resources,China(Grant No.ZJDZGCZ-2021).
文摘The geometric characteristics of fractures within a rock mass can be inferred by the data sampling from boreholes or exposed surfaces.Recently,the universal elliptical disc(UED)model was developed to represent natural fractures,where the fracture is assumed to be an elliptical disc and the fracture orientation,rotation angle,length of the long axis and ratio of short-long axis lengths are considered as variables.This paper aims to estimate the fracture size-and azimuth-related parameters in the UED model based on the trace information from sampling windows.The stereological relationship between the trace length,size-and azimuth-related parameters of the UED model was established,and the formulae of the mean value and standard deviation of trace length were proposed.The proposed formulae were validated via the Monte Carlo simulations with less than 5%of error rate between the calculated and true values.With respect to the estimation of the size-and azimuth-related parameters using the trace length,an optimization method was developed based on the pre-assumed size and azimuth distribution forms.A hypothetical case study was designed to illustrate and verify the parameter estimation method,where three combinations of the sampling windows were used to estimate the parameters,and the results showed that the estimated values could agree well with the true values.Furthermore,a hypothetical three-dimensional(3D)elliptical fracture network was constructed,and the circular disc,non-UED and UED models were used to represent it.The simulated trace information from different models was compared,and the results clearly illustrated the superiority of the proposed UED model over the existing circular disc and non-UED models。
基金The project supported by the National Natural Science Foundation of China
文摘The nonlocal boundary value problems for nonlinear elliptic systems in the unbounded domain are considered. Under suitable conditions the existence of solution and comparison theorem for the boundary value problems are studied.
基金Supported by the National Science Foundation of China(11401162,11571089,11401159,11301136)the Natural Science Foundation of Hebei Province(A2015205012,A2016205218,A2014205069,A2014208158)Hebei Normal University Dr.Fund(L2015B03)
文摘Firstly, the Riemann boundary value problem for a kind of degenerate elliptic sys- tem of the first order equations in R4 is proposed. Then, with the help of the one-to-one correspondence between the theory of Clifford valued generalized regular functions and that of the degenerate elliptic system's solution, the boundary value problem as stated above is trans- formed into a boundary value problem related to the generalized regular functions in Clifford analysis. Moreover, the solution of the Riemann boundary value problem for the degenerate elliptic system is explicitly described by using a kind of singular integral operator. Finally, the conditions for the existence of solutions of the oblique derivative problem for another kind of degenerate elliptic system of the first order equations in R4 are derived.
文摘In this paper a singular perturbation of boundary value problem for elliptic partial differential equations of higher order is considered by using the differential inequalities. The uniformly valid asymptotic expansion in entire region is obtained.
基金This work is supported by the Foundatiorl of Zhongshan University Advanced Research Centre
文摘Some superconvergence results of generalized difference solution for elliptic boundary value problem are given. It is shown that optimal points of the stresses for generalized difference method are the same as that for finite element method.
文摘This paper deals with the Neumann problem for a class of semilinear elliptic equations -△u + u =|u|2*-2u+ μ|u|q-2u in Ω, au/ar= |u|(?)*-2u on aΩ, where 2 = 2N/N-2, s=2(N-1)/N-2, 1 <q<2,N(?)3,μ>γ denotes the unit outward normal to boundary aΩ. By vaxiational method and dual fountain theorem, the existence of infinitely many solutions with negative energy is proved.
文摘In this paper, the first boundary value problem for quasilinear equation of the form △A(u,x)+∑i=1^m δb^i(u,x)/δxi+c(u,x)=0,Au(u,x) ≥0is studied. By using the compensated compactness theory, some results on the existence of weak solution are established. In addition, under certain condition the uniqueness of solution is proved.
文摘In this paper, authors discuss the numerical methods of general discontinuous boundary value problems for elliptic complex equations of first order, They first give the well posedness of general discontinuous boundary value problems, reduce the discontinuous boundary value problems to a variation problem, and then find the numerical solutions of above problem by the finite element method. Finally authors give some error-estimates of the foregoing numerical solutions.