Transient behavior of three-dimensional semiconductor device with heat conduc- tion is described by a coupled mathematical system of four quasi-linear partial differential equations with initial-boundary value conditi...Transient behavior of three-dimensional semiconductor device with heat conduc- tion is described by a coupled mathematical system of four quasi-linear partial differential equations with initial-boundary value conditions. The electric potential is defined by an ellip- tic equation and it appears in the following three equations via the electric field intensity. The electron concentration and the hole concentration are determined by convection-dominated diffusion equations and the temperature is interpreted by a heat conduction equation. A mixed finite volume element approximation, keeping physical conservation law, is used to get numerical values of the electric potential and the accuracy is improved one order. Two con- centrations and the heat conduction are computed by a fractional step method combined with second-order upwind differences. This method can overcome numerical oscillation, dispersion and decreases computational complexity. Then a three-dimensional problem is solved by computing three successive one-dimensional problems where the method of speedup is used and the computational work is greatly shortened. An optimal second-order error estimate in L2 norm is derived by using prior estimate theory and other special techniques of partial differential equations. This type of mass-conservative parallel method is important and is most valuable in numerical analysis and application of semiconductor device.展开更多
In this paper,by using trapezoidal rule and the integration-by-parts formula of Malliavin calculus,we propose three new numerical schemes for solving decoupled forward-backward stochastic differential equations.We the...In this paper,by using trapezoidal rule and the integration-by-parts formula of Malliavin calculus,we propose three new numerical schemes for solving decoupled forward-backward stochastic differential equations.We theoretically prove that the schemes have second-order convergence rate.To demonstrate the effectiveness and the second-order convergence rate,numerical tests are given.展开更多
We propose new numerical schemes for decoupled forward-backward stochastic differ- ential equations （FBSDEs） with jumps, where the stochastic dynamics are driven by a d- dimensional Brownian motion and an independen...We propose new numerical schemes for decoupled forward-backward stochastic differ- ential equations （FBSDEs） with jumps, where the stochastic dynamics are driven by a d- dimensional Brownian motion and an independent compensated Poisson random measure. A semi-discrete scheme is developed for discrete time approximation, which is constituted by a classic scheme for the forward SDE [20, 28] and a novel scheme for the backward SDE. Under some reasonable regularity conditions, we prove that the semi-discrete scheme can achieve second-order convergence in approximating the FBSDEs of interest; and such convergence rate does not require jump-adapted temporal discretization. Next, to add in spatial discretization, a fully discrete scheme is developed by designing accurate quadrature rules for estimating the involved conditional mathematical expectations. Several numerical examples are given to illustrate the effectiveness and the high accuracy of the proposed schemes.展开更多
In this paper, under weak conditions, we theoretically prove the second-order convergence rate of the Crank-Nicolson scheme for solving a kind of decoupled forward-backward stochastic differential equations.
In this paper, a new discrete formulation and a type of new posteriori error estimators for the second-order element discretization for Stokes problems are presented, where pressure is approximated with piecewise fir...In this paper, a new discrete formulation and a type of new posteriori error estimators for the second-order element discretization for Stokes problems are presented, where pressure is approximated with piecewise first-degree polynomials and velocity vector field with piecewise second-degree polynomials with a cubic bubble function to be added. The estimators are the globally upper and locally lower bounds for the error of the finite element discretization. It is shown that the bubble part for this second-order element approximation is substituted for the other parts of the approximate solution.展开更多
This paper considers the so-called expected residual minimization(ERM)formulation for stochastic second-order cone complementarity problems,which is based on a new complementarity function called termwise residual com...This paper considers the so-called expected residual minimization(ERM)formulation for stochastic second-order cone complementarity problems,which is based on a new complementarity function called termwise residual complementarity function associated with second-order cone.We show that the ERM model has bounded level sets under the stochastic weak R0-property.We further derive some error bound results under either the strong monotonicity or some kind of constraint qualifications.Then,we apply the Monte Carlo approximation techniques to solve the ERM model and establish a comprehensive convergence analysis.Furthermore,we report some numerical results on a stochastic second-order cone model for optimal power flow in radial networks.展开更多
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 paper,a nonlinear time-fractional coupled diffusion system is solved by using a mixed finite element(MFE)method in space combined with L1-approximation and implicit second-order backward difference scheme in t...In this paper,a nonlinear time-fractional coupled diffusion system is solved by using a mixed finite element(MFE)method in space combined with L1-approximation and implicit second-order backward difference scheme in time.The stability for nonlinear fully discrete finite element scheme is analyzed and a priori error estimates are derived.Finally,some numerical tests are shown to verify our theoretical analysis.展开更多
基金supported by National Natural Science Foundation of China(11101244,11271231)National Tackling Key Problems Program(20050200069)Doctorate Foundation of the Ministry of Education of China(20030422047)
文摘Transient behavior of three-dimensional semiconductor device with heat conduc- tion is described by a coupled mathematical system of four quasi-linear partial differential equations with initial-boundary value conditions. The electric potential is defined by an ellip- tic equation and it appears in the following three equations via the electric field intensity. The electron concentration and the hole concentration are determined by convection-dominated diffusion equations and the temperature is interpreted by a heat conduction equation. A mixed finite volume element approximation, keeping physical conservation law, is used to get numerical values of the electric potential and the accuracy is improved one order. Two con- centrations and the heat conduction are computed by a fractional step method combined with second-order upwind differences. This method can overcome numerical oscillation, dispersion and decreases computational complexity. Then a three-dimensional problem is solved by computing three successive one-dimensional problems where the method of speedup is used and the computational work is greatly shortened. An optimal second-order error estimate in L2 norm is derived by using prior estimate theory and other special techniques of partial differential equations. This type of mass-conservative parallel method is important and is most valuable in numerical analysis and application of semiconductor device.
基金supported by National Natural Science Foundation of China (Grant Nos. 91130003 and 11171189)Natural Science Foundation of Shandong Province (Grant No. ZR2011AZ002)
文摘In this paper,by using trapezoidal rule and the integration-by-parts formula of Malliavin calculus,we propose three new numerical schemes for solving decoupled forward-backward stochastic differential equations.We theoretically prove that the schemes have second-order convergence rate.To demonstrate the effectiveness and the second-order convergence rate,numerical tests are given.
基金The authors would like to thank the referees for their valuable comments, which improved much of the quality of the paper. This work is partially support- ed by the National Natural Science Foundations of China under grant numbers 91130003,11171189 and 11571206 by Natural Science Foundation of Shandong Province under grant number ZR2011AZ002+1 种基金 by the U.S. Department of Energy, Office of Science, Office of Ad- vanced Scientific Computing Research, Applied Mathematics program under contract number ERKJE45 and by the Laboratory Directed Research and Development program at the Oak Ridge National Laboratory, which is operated by UT-Battelle, LLC, for the U.S. Department of Energy under Contract DE-AC05-00OR22725.
文摘We propose new numerical schemes for decoupled forward-backward stochastic differ- ential equations （FBSDEs） with jumps, where the stochastic dynamics are driven by a d- dimensional Brownian motion and an independent compensated Poisson random measure. A semi-discrete scheme is developed for discrete time approximation, which is constituted by a classic scheme for the forward SDE [20, 28] and a novel scheme for the backward SDE. Under some reasonable regularity conditions, we prove that the semi-discrete scheme can achieve second-order convergence in approximating the FBSDEs of interest; and such convergence rate does not require jump-adapted temporal discretization. Next, to add in spatial discretization, a fully discrete scheme is developed by designing accurate quadrature rules for estimating the involved conditional mathematical expectations. Several numerical examples are given to illustrate the effectiveness and the high accuracy of the proposed schemes.
文摘In this paper, under weak conditions, we theoretically prove the second-order convergence rate of the Crank-Nicolson scheme for solving a kind of decoupled forward-backward stochastic differential equations.
基金the National Natural Science Foundation of China, the Evolvement Plan of Science and Technology of Beijing Educational Council,
文摘In this paper, a new discrete formulation and a type of new posteriori error estimators for the second-order element discretization for Stokes problems are presented, where pressure is approximated with piecewise first-degree polynomials and velocity vector field with piecewise second-degree polynomials with a cubic bubble function to be added. The estimators are the globally upper and locally lower bounds for the error of the finite element discretization. It is shown that the bubble part for this second-order element approximation is substituted for the other parts of the approximate solution.
基金This work was supported in part by the National Natural Science Foundation of China(Nos.71831008,11671250,11431004 and 11601458)Humanity and Social Science Foundation of Ministry of Education of China(No.15YJA630034)+2 种基金Shandong Province Natural Science Fund(No.ZR2014AM012)Higher Educational Science and Technology Program of Shandong Province(No.J13LI09)Scientific Research of Young Scholar of Qufu Normal University(No.XKJ201315).
文摘This paper considers the so-called expected residual minimization(ERM)formulation for stochastic second-order cone complementarity problems,which is based on a new complementarity function called termwise residual complementarity function associated with second-order cone.We show that the ERM model has bounded level sets under the stochastic weak R0-property.We further derive some error bound results under either the strong monotonicity or some kind of constraint qualifications.Then,we apply the Monte Carlo approximation techniques to solve the ERM model and establish a comprehensive convergence analysis.Furthermore,we report some numerical results on a stochastic second-order cone model for optimal power flow in radial networks.
基金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.
基金the National Natural Science Fund(11661058,11301258,11361035)the Natural Science Fund of Inner Mongolia Autonomous Region(2016MS0102,2015MS0101)+1 种基金the Scientific Research Projection of Higher Schools of Inner Mongolia(NJZZ12011)the National Undergraduate Innovative Training Project(201510126026).
文摘In this paper,a nonlinear time-fractional coupled diffusion system is solved by using a mixed finite element(MFE)method in space combined with L1-approximation and implicit second-order backward difference scheme in time.The stability for nonlinear fully discrete finite element scheme is analyzed and a priori error estimates are derived.Finally,some numerical tests are shown to verify our theoretical analysis.
