We consider the singular Riemann problem for the rectilinear isentropic compressible Euler equations with discontinuous flux,more specifically,for pressureless flow on the left and polytropic flow on the right separat...We consider the singular Riemann problem for the rectilinear isentropic compressible Euler equations with discontinuous flux,more specifically,for pressureless flow on the left and polytropic flow on the right separated by a discontinuity x=x(t).We prove that this problem admits global Radon measure solutions for all kinds of initial data.The over-compressing condition on the discontinuity x=x(t)is not enough to ensure the uniqueness of the solution.However,there is a unique piecewise smooth solution if one proposes a slip condition on the right-side of the curve x=x(t)+0,in addition to the full adhesion condition on its left-side.As an application,we study a free piston problem with the piston in a tube surrounded initially by uniform pressureless flow and a polytropic gas.In particular,we obtain the existence of a piecewise smooth solution for the motion of the piston between a vacuum and a polytropic gas.This indicates that the singular Riemann problem looks like a control problem in the sense that one could adjust the condition on the discontinuity of the flux to obtain the desired flow field.展开更多
This paper studies the non-homogeneous generalized Riemann-Hilbert(RH)problems involving two unknown functions.Using the uniformization theorem,such problems are transformed into the case of homogeneous type.By the th...This paper studies the non-homogeneous generalized Riemann-Hilbert(RH)problems involving two unknown functions.Using the uniformization theorem,such problems are transformed into the case of homogeneous type.By the theory of classical boundary value problems,we adopt a novel method to obtain the sectionally analytic solutions of problems in strip domains,and analyze the conditions of solvability and properties of solutions in various domains.展开更多
In this article, we first introduce the general linear elliptic complex equation of first order with certain conditions, and then propose discontinuous Riemann-Hilbert problem and some kinds of modified well-posed-nes...In this article, we first introduce the general linear elliptic complex equation of first order with certain conditions, and then propose discontinuous Riemann-Hilbert problem and some kinds of modified well-posed-ness for the complex equation. Then we verify the equivalence of three kinds of well-posed-ness. The discontinuous boundary value problem possesses many applications in mechanics and physics etc.展开更多
Several approximate methods have been used to find approximate solutions of elliptic systems of first order equations. One common method is the Newton imbedding approach, i.e. the parameter extension method. In this a...Several approximate methods have been used to find approximate solutions of elliptic systems of first order equations. One common method is the Newton imbedding approach, i.e. the parameter extension method. In this article, we discuss approximate solutions to discontinuous Riemann-Hilbert boundary value problems, which have various applications in mechanics and physics. We first formulate the discontinuous Riemann-Hilbert problem for elliptic systems of first order complex equations in multiply connected domains and its modified well-posedness, then use the parameter extensional method to find approximate solutions to the modified boundary value problem for elliptic complex systems of first order equations, and then provide the error estimate of approximate solutions for the discontinuous boundary value problem.展开更多
A class of discontinuous penalty functions was proposed to solve constrained minimization problems with the integral approach to global optimization, m-mean value and v-variance optimality conditions of a constrained ...A class of discontinuous penalty functions was proposed to solve constrained minimization problems with the integral approach to global optimization, m-mean value and v-variance optimality conditions of a constrained and penalized minimization problem were investigated. A nonsequential algorithm was proposed. Numerical examples were given to illustrate the effectiveness of the algorithm.展开更多
In this paper,a fully discrete stability analysis is carried out for the direct discontinuous Galerkin(DDG)methods coupled with Runge-Kutta-type implicit-explicit time marching,for solving one-dimensional linear conve...In this paper,a fully discrete stability analysis is carried out for the direct discontinuous Galerkin(DDG)methods coupled with Runge-Kutta-type implicit-explicit time marching,for solving one-dimensional linear convection-diffusion problems.In the spatial discretization,both the original DDG methods and the refined DDG methods with interface corrections are considered.In the time discretization,the convection term is treated explicitly and the diffusion term implicitly.By the energy method,we show that the corresponding fully discrete schemes are unconditionally stable,in the sense that the time-stepis only required to be upper bounded by a constant which is independent of the mesh size h.Opti-mal error estimate is also obtained by the aid of a special global projection.Numerical experiments are given to verify the stability and accuracy of the proposed schemes.展开更多
In this paper,the minimal dissipation local discontinuous Galerkin method is studied to solve the elliptic interface problems in two-dimensional domains.The interface may be arbitrary smooth curves.It is shown that th...In this paper,the minimal dissipation local discontinuous Galerkin method is studied to solve the elliptic interface problems in two-dimensional domains.The interface may be arbitrary smooth curves.It is shown that the error estimates in L;-norm for the solution and the flux are O(h;|log h|)and O(h|log h|;),respectively.In numerical experiments,the successive substitution iterative methods are used to solve the LDG schemes.Numerical results verify the efficiency and accuracy of the method.展开更多
In this paper nonlinear boundary value problems for discontinuous delayed differen- tial equations are considered. Some existence and boundedness results of solutions are obtained via the method of upper and lower sol...In this paper nonlinear boundary value problems for discontinuous delayed differen- tial equations are considered. Some existence and boundedness results of solutions are obtained via the method of upper and lower solutions, which may be discontinuous. Our analysis can be applied to those phenomena from physics and control theory which have been successfully described by delayed differential equations with discontinuous functions, such as electric, pneu- matic, and hydraulic networks. It is also an important step to precede the design of control signals when finite-time or practical stability control are concerned.展开更多
In this paper,we consider the inverse resonance problems for the discontinuous and non-selfadjoint Sturm-Liouville problem.We prove the uniqueness theorem and provide a reconstructive algorithm for the potential by us...In this paper,we consider the inverse resonance problems for the discontinuous and non-selfadjoint Sturm-Liouville problem.We prove the uniqueness theorem and provide a reconstructive algorithm for the potential by using the Cauchy data and Weyl function.展开更多
We extend the fnite element method introduced by Lakkis and Pryer(SIAM J.Sci.Comput.33(2):786–801,2011)to approximate the solution of second-order elliptic problems in nonvariational form to incorporate the discontin...We extend the fnite element method introduced by Lakkis and Pryer(SIAM J.Sci.Comput.33(2):786–801,2011)to approximate the solution of second-order elliptic problems in nonvariational form to incorporate the discontinuous Galerkin(DG)framework.This is done by viewing the“fnite element Hessian”as an auxiliary variable in the formulation.Representing the fnite element Hessian in a discontinuous setting yields a linear system of the same size and having the same sparsity pattern of the compact DG methods for variational elliptic problems.Furthermore,the system matrix is very easy to assemble;thus,this approach greatly reduces the computational complexity of the discretisation compared to the continuous approach.We conduct a stability and consistency analysis making use of the unifed frameworkset out in Arnold et al.(SIAM J.Numer.Anal.39(5):1749–1779,2001/2002).We also give an a posteriori analysis of the method in the case where the problem has a strong solution.The analysis applies to any consistent representation of the fnite element Hessian,and thus is applicable to the previous works making use of continuous Galerkin approximations.Numerical evidence is presented showing that the method works well also in a more general setting.展开更多
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.展开更多
In this study, the boundary-value problem with eigenvalue parameter generated by the differential equation with discontinuous coefficients and boundary conditions which contains not only endpoints of the considered in...In this study, the boundary-value problem with eigenvalue parameter generated by the differential equation with discontinuous coefficients and boundary conditions which contains not only endpoints of the considered interval, but also point of discontinuity and linear functionals is investigated. So, the problem is not pure boundary-value. The authors single out a class of linear functionals and find simple algebraic conditions on coefficients, which garantee the existence of infinit number eigenvalues. Also the asymptotic formulas for eigenvalues are found.展开更多
Based on a 4 x 4 matrix spectral problem, an AKNS soliton hierarchy with six potentials is generated. Associated with this spectral problem, a kind of Riemann-Hilbert problems is formulated for a six-component system ...Based on a 4 x 4 matrix spectral problem, an AKNS soliton hierarchy with six potentials is generated. Associated with this spectral problem, a kind of Riemann-Hilbert problems is formulated for a six-component system of mKdV equations in the resulting AKNS hierarchy. Soliton solutions to the considered system of coupled mKdV equations are computed, through a reduced Riemann-Hilbert problem where an identity jump matrix is taken.展开更多
The purpose of this paper is to prove existence of minimisers of the functional where Ω is an open set of the Heisenberg group Hn, K runs over all closed sets of Hn, u varies in C_H^1(Ω\ K), α,β> 0,q≥1, g ∈ ...The purpose of this paper is to prove existence of minimisers of the functional where Ω is an open set of the Heisenberg group Hn, K runs over all closed sets of Hn, u varies in C_H^1(Ω\ K), α,β> 0,q≥1, g ∈ Lq(Ω) ∩ L∞(Ω) and f : R2n→R is a convex function satisfying some structure conditions (H1)(H2)(H3) (see below).展开更多
This paper studies the nonlinear mixed problem for a class of symmetric hyperbolic systems with the boundary condition satisfying the dissipative condition about discontinuous data in higher dimension spaces, establis...This paper studies the nonlinear mixed problem for a class of symmetric hyperbolic systems with the boundary condition satisfying the dissipative condition about discontinuous data in higher dimension spaces, establishes the local existence theorem by using the method of a prior estimates, and obtains the structure of singularities of the solutions of such problems.展开更多
The paper aims at establishing Riemann-Hilbert problems and presenting soliton solutions for nonlocal reverse-time nonlinear Schrodinger(NLS) hierarchies associated with higher-order matrix spectral problems.The Sokho...The paper aims at establishing Riemann-Hilbert problems and presenting soliton solutions for nonlocal reverse-time nonlinear Schrodinger(NLS) hierarchies associated with higher-order matrix spectral problems.The Sokhotski-Plemelj formula is used to transform the Riemann-Hilbert problems into Gelfand-Levitan-Marchenko type integral equations.A new formulation of solutions to special Riemann-Hilbert problems with the identity jump matrix,corresponding to the reflectionless inverse scattering transforms,is proposed and applied to construction of soliton solutions to each system in the considered nonlocal reversetime NLS hierarchies.展开更多
In this paper, we study Riemann boundary value problems on the Curve of Parabola. We characterized the functions which are intergrable on the Curve of Parabola. We also study the asymptotic behaviors of Cauchy-type in...In this paper, we study Riemann boundary value problems on the Curve of Parabola. We characterized the functions which are intergrable on the Curve of Parabola. We also study the asymptotic behaviors of Cauchy-type integral and Cauchy principal value integral on the Curve of Parabola at infinity. At the end, we discuss the Riemann boundary value problems for sectionally holomorphic functions with the Curve of Parabola as their jump curve and obtain the explicit form.展开更多
In this paper, the mixed initial-boundary value problem for general first order quasi- linear hyperbolic systems with nonlinear boundary conditions in the domain D = {(t, x) | t ≥ 0, x ≥0} is considered. A suffic...In this paper, the mixed initial-boundary value problem for general first order quasi- linear hyperbolic systems with nonlinear boundary conditions in the domain D = {(t, x) | t ≥ 0, x ≥0} is considered. A sufficient condition to guarantee the existence and uniqueness of global weakly discontinuous solution is given.展开更多
This paper is concerned with the initial-boundary value problem of scalar conservation laws with weak discontinuous flux, whose initial data are a function with two pieces of constant and whose boundary data are a con...This paper is concerned with the initial-boundary value problem of scalar conservation laws with weak discontinuous flux, whose initial data are a function with two pieces of constant and whose boundary data are a constant function. Under the condition that the flux function has a finite number of weak discontinuous points, by using the structure of weak entropy solution of the corresponding initial value problem and the boundary entropy condition developed by Bardos-Leroux-Nedelec, we give a construction method to the global weak entropy solution for this initial-boundary value problem, and by investigating the interaction of elementary waves and the boundary, we clarify the geometric structure and the behavior of boundary for the weak entropy solution.展开更多
In this article, we introduce a coupled approach of local discontinuous Calerkin and standard finite element method for solving convection diffusion problems. The whole domain is divided into two disjoint subdomains. ...In this article, we introduce a coupled approach of local discontinuous Calerkin and standard finite element method for solving convection diffusion problems. The whole domain is divided into two disjoint subdomains. The discontinuous Galerkin method is adopted in the subdomain where the solution varies rapidly, while the standard finite element method is used in the other subdomain due to its lower computational cost. The stability and a priori error estimate are established. We prove that the coupled method has O(ε1/2 + h1/2)hk) convergence rate in an associated norm, where ε is the diffusion coefficient, h is the mesh size and k is the degree of polynomial. The numerical results verify our theoretical results. Moreover, 2k-order superconvergence of the numerical traces at the nodes, and the optimal convergence of the errors under L2 norm are observed numerically on the uniform mesh. The numerical results also indicate that the coupled method has the same convergence order and almost the same errors as the purely LDG method.展开更多
基金supported by the National Natural Science Foundation of China(11871218,12071298)in part by the Science and Technology Commission of Shanghai Municipality(21JC1402500,22DZ2229014)。
文摘We consider the singular Riemann problem for the rectilinear isentropic compressible Euler equations with discontinuous flux,more specifically,for pressureless flow on the left and polytropic flow on the right separated by a discontinuity x=x(t).We prove that this problem admits global Radon measure solutions for all kinds of initial data.The over-compressing condition on the discontinuity x=x(t)is not enough to ensure the uniqueness of the solution.However,there is a unique piecewise smooth solution if one proposes a slip condition on the right-side of the curve x=x(t)+0,in addition to the full adhesion condition on its left-side.As an application,we study a free piston problem with the piston in a tube surrounded initially by uniform pressureless flow and a polytropic gas.In particular,we obtain the existence of a piecewise smooth solution for the motion of the piston between a vacuum and a polytropic gas.This indicates that the singular Riemann problem looks like a control problem in the sense that one could adjust the condition on the discontinuity of the flux to obtain the desired flow field.
基金Supported by National Natural Science Foundation of China(Grant No.11971015).
文摘This paper studies the non-homogeneous generalized Riemann-Hilbert(RH)problems involving two unknown functions.Using the uniformization theorem,such problems are transformed into the case of homogeneous type.By the theory of classical boundary value problems,we adopt a novel method to obtain the sectionally analytic solutions of problems in strip domains,and analyze the conditions of solvability and properties of solutions in various domains.
文摘In this article, we first introduce the general linear elliptic complex equation of first order with certain conditions, and then propose discontinuous Riemann-Hilbert problem and some kinds of modified well-posed-ness for the complex equation. Then we verify the equivalence of three kinds of well-posed-ness. The discontinuous boundary value problem possesses many applications in mechanics and physics etc.
文摘Several approximate methods have been used to find approximate solutions of elliptic systems of first order equations. One common method is the Newton imbedding approach, i.e. the parameter extension method. In this article, we discuss approximate solutions to discontinuous Riemann-Hilbert boundary value problems, which have various applications in mechanics and physics. We first formulate the discontinuous Riemann-Hilbert problem for elliptic systems of first order complex equations in multiply connected domains and its modified well-posedness, then use the parameter extensional method to find approximate solutions to the modified boundary value problem for elliptic complex systems of first order equations, and then provide the error estimate of approximate solutions for the discontinuous boundary value problem.
文摘A class of discontinuous penalty functions was proposed to solve constrained minimization problems with the integral approach to global optimization, m-mean value and v-variance optimality conditions of a constrained and penalized minimization problem were investigated. A nonsequential algorithm was proposed. Numerical examples were given to illustrate the effectiveness of the algorithm.
基金the NSFC grant 11871428the Nature Science Research Program for Colleges and Universities of Jiangsu Province grant 20KJB110011Qiang Zhang:Research supported by the NSFC grant 11671199。
文摘In this paper,a fully discrete stability analysis is carried out for the direct discontinuous Galerkin(DDG)methods coupled with Runge-Kutta-type implicit-explicit time marching,for solving one-dimensional linear convection-diffusion problems.In the spatial discretization,both the original DDG methods and the refined DDG methods with interface corrections are considered.In the time discretization,the convection term is treated explicitly and the diffusion term implicitly.By the energy method,we show that the corresponding fully discrete schemes are unconditionally stable,in the sense that the time-stepis only required to be upper bounded by a constant which is independent of the mesh size h.Opti-mal error estimate is also obtained by the aid of a special global projection.Numerical experiments are given to verify the stability and accuracy of the proposed schemes.
基金Supported by National Natural Science Foundation of China(11571002,11461046)Natural Science Foundation of Jiangxi Province,China(20151BAB211013,20161ACB21005)+2 种基金Science and Technology Project of Jiangxi Provincial Department of Education,China(150172)Science Foundation of China Academy of Engineering Physics(2015B0101021)Defense Industrial Technology Development Program(B1520133015)
文摘In this paper,the minimal dissipation local discontinuous Galerkin method is studied to solve the elliptic interface problems in two-dimensional domains.The interface may be arbitrary smooth curves.It is shown that the error estimates in L;-norm for the solution and the flux are O(h;|log h|)and O(h|log h|;),respectively.In numerical experiments,the successive substitution iterative methods are used to solve the LDG schemes.Numerical results verify the efficiency and accuracy of the method.
基金Supported by the National Natural Science Foundation of China(Grants No.70703016 and No.10001024)Research Grant of the Business School of Nanjing University
文摘In this paper nonlinear boundary value problems for discontinuous delayed differen- tial equations are considered. Some existence and boundedness results of solutions are obtained via the method of upper and lower solutions, which may be discontinuous. Our analysis can be applied to those phenomena from physics and control theory which have been successfully described by delayed differential equations with discontinuous functions, such as electric, pneu- matic, and hydraulic networks. It is also an important step to precede the design of control signals when finite-time or practical stability control are concerned.
基金supported in part by the National Natural Science Foundation of China (11871031)the Natural Science Foundation of the Jiangsu Province of China(BK 20201303)supported in part by Postgraduate Research and Practice Innovation Program of Jiangsu Province (KYCX20 0245)
文摘In this paper,we consider the inverse resonance problems for the discontinuous and non-selfadjoint Sturm-Liouville problem.We prove the uniqueness theorem and provide a reconstructive algorithm for the potential by using the Cauchy data and Weyl function.
文摘We extend the fnite element method introduced by Lakkis and Pryer(SIAM J.Sci.Comput.33(2):786–801,2011)to approximate the solution of second-order elliptic problems in nonvariational form to incorporate the discontinuous Galerkin(DG)framework.This is done by viewing the“fnite element Hessian”as an auxiliary variable in the formulation.Representing the fnite element Hessian in a discontinuous setting yields a linear system of the same size and having the same sparsity pattern of the compact DG methods for variational elliptic problems.Furthermore,the system matrix is very easy to assemble;thus,this approach greatly reduces the computational complexity of the discretisation compared to the continuous approach.We conduct a stability and consistency analysis making use of the unifed frameworkset out in Arnold et al.(SIAM J.Numer.Anal.39(5):1749–1779,2001/2002).We also give an a posteriori analysis of the method in the case where the problem has a strong solution.The analysis applies to any consistent representation of the fnite element Hessian,and thus is applicable to the previous works making use of continuous Galerkin approximations.Numerical evidence is presented showing that the method works well also in a more general setting.
基金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.
文摘In this study, the boundary-value problem with eigenvalue parameter generated by the differential equation with discontinuous coefficients and boundary conditions which contains not only endpoints of the considered interval, but also point of discontinuity and linear functionals is investigated. So, the problem is not pure boundary-value. The authors single out a class of linear functionals and find simple algebraic conditions on coefficients, which garantee the existence of infinit number eigenvalues. Also the asymptotic formulas for eigenvalues are found.
基金supported in part by NSFC(11371326,11301331,and 11371086)NSF under the grant DMS-1664561+2 种基金the 111 project of China(B16002)the China state administration of foreign experts affairs system under the affiliation of North China Electric Power University,Natural Science Fund for Colleges and Universities of Jiangsu Province under the grant 17KJB110020the Distinguished Professorships by Shanghai University of Electric Power,China and North-West University,South Africa
文摘Based on a 4 x 4 matrix spectral problem, an AKNS soliton hierarchy with six potentials is generated. Associated with this spectral problem, a kind of Riemann-Hilbert problems is formulated for a six-component system of mKdV equations in the resulting AKNS hierarchy. Soliton solutions to the considered system of coupled mKdV equations are computed, through a reduced Riemann-Hilbert problem where an identity jump matrix is taken.
基金This work is supported by NNSF(10471063), Hunan NSF(03JJY4002) & Hunan Education Administration Item(03A011)
文摘The purpose of this paper is to prove existence of minimisers of the functional where Ω is an open set of the Heisenberg group Hn, K runs over all closed sets of Hn, u varies in C_H^1(Ω\ K), α,β> 0,q≥1, g ∈ Lq(Ω) ∩ L∞(Ω) and f : R2n→R is a convex function satisfying some structure conditions (H1)(H2)(H3) (see below).
文摘This paper studies the nonlinear mixed problem for a class of symmetric hyperbolic systems with the boundary condition satisfying the dissipative condition about discontinuous data in higher dimension spaces, establishes the local existence theorem by using the method of a prior estimates, and obtains the structure of singularities of the solutions of such problems.
基金supported in part by NSFC(11975145 and 11972291)the Natural Science Foundation for Colleges and Universities in Jiangsu Province(17 KJB 110020)。
文摘The paper aims at establishing Riemann-Hilbert problems and presenting soliton solutions for nonlocal reverse-time nonlinear Schrodinger(NLS) hierarchies associated with higher-order matrix spectral problems.The Sokhotski-Plemelj formula is used to transform the Riemann-Hilbert problems into Gelfand-Levitan-Marchenko type integral equations.A new formulation of solutions to special Riemann-Hilbert problems with the identity jump matrix,corresponding to the reflectionless inverse scattering transforms,is proposed and applied to construction of soliton solutions to each system in the considered nonlocal reversetime NLS hierarchies.
文摘In this paper, we study Riemann boundary value problems on the Curve of Parabola. We characterized the functions which are intergrable on the Curve of Parabola. We also study the asymptotic behaviors of Cauchy-type integral and Cauchy principal value integral on the Curve of Parabola at infinity. At the end, we discuss the Riemann boundary value problems for sectionally holomorphic functions with the Curve of Parabola as their jump curve and obtain the explicit form.
文摘In this paper, the mixed initial-boundary value problem for general first order quasi- linear hyperbolic systems with nonlinear boundary conditions in the domain D = {(t, x) | t ≥ 0, x ≥0} is considered. A sufficient condition to guarantee the existence and uniqueness of global weakly discontinuous solution is given.
文摘This paper is concerned with the initial-boundary value problem of scalar conservation laws with weak discontinuous flux, whose initial data are a function with two pieces of constant and whose boundary data are a constant function. Under the condition that the flux function has a finite number of weak discontinuous points, by using the structure of weak entropy solution of the corresponding initial value problem and the boundary entropy condition developed by Bardos-Leroux-Nedelec, we give a construction method to the global weak entropy solution for this initial-boundary value problem, and by investigating the interaction of elementary waves and the boundary, we clarify the geometric structure and the behavior of boundary for the weak entropy solution.
基金Supported by National Natural Science Foundation of China (10571046, 10571053, and 10871066)Program for New Century Excellent Talents in University (NCET-06-0712)+2 种基金Key Laboratory of Computational and Stochastic Mathematics and Its Applications, Universities of Hunan Province, Hunan Normal Universitythe Project of Scientific Research Fund of Hunan Provincial Education Department (09K025)the Key Scientific Research Topic of Jiaxing University (70110X05BL)
文摘In this article, we introduce a coupled approach of local discontinuous Calerkin and standard finite element method for solving convection diffusion problems. The whole domain is divided into two disjoint subdomains. The discontinuous Galerkin method is adopted in the subdomain where the solution varies rapidly, while the standard finite element method is used in the other subdomain due to its lower computational cost. The stability and a priori error estimate are established. We prove that the coupled method has O(ε1/2 + h1/2)hk) convergence rate in an associated norm, where ε is the diffusion coefficient, h is the mesh size and k is the degree of polynomial. The numerical results verify our theoretical results. Moreover, 2k-order superconvergence of the numerical traces at the nodes, and the optimal convergence of the errors under L2 norm are observed numerically on the uniform mesh. The numerical results also indicate that the coupled method has the same convergence order and almost the same errors as the purely LDG method.