This paper is devoted to the Polynomial Preserving Recovery (PPR) based a posteriori error analysis for the second-order elliptic non-symmetric eigenvalue problem. An asymptotically exact a posteriori error estimator ...This paper is devoted to the Polynomial Preserving Recovery (PPR) based a posteriori error analysis for the second-order elliptic non-symmetric eigenvalue problem. An asymptotically exact a posteriori error estimator is proposed for solving the convection-dominated non-symmetric eigenvalue problem with non-smooth eigenfunctions or multiple eigenvalues. Numerical examples confirm our theoretical analysis.展开更多
During the use of robotics in applications such as antiterrorism or combat,a motion-constrained pursuer vehicle,such as a Dubins unmanned surface vehicle(USV),must get close enough(within a prescribed zero or positive...During the use of robotics in applications such as antiterrorism or combat,a motion-constrained pursuer vehicle,such as a Dubins unmanned surface vehicle(USV),must get close enough(within a prescribed zero or positive distance)to a moving target as quickly as possible,resulting in the extended minimum-time intercept problem(EMTIP).Existing research has primarily focused on the zero-distance intercept problem,MTIP,establishing the necessary or sufficient conditions for MTIP optimality,and utilizing analytic algorithms,such as root-finding algorithms,to calculate the optimal solutions.However,these approaches depend heavily on the properties of the analytic algorithm,making them inapplicable when problem settings change,such as in the case of a positive effective range or complicated target motions outside uniform rectilinear motion.In this study,an approach employing a high-accuracy and quality-guaranteed mixed-integer piecewise-linear program(QG-PWL)is proposed for the EMTIP.This program can accommodate different effective interception ranges and complicated target motions(variable velocity or complicated trajectories).The high accuracy and quality guarantees of QG-PWL originate from elegant strategies such as piecewise linearization and other developed operation strategies.The approximate error in the intercept path length is proved to be bounded to h^(2)/(4√2),where h is the piecewise length.展开更多
Errors inevitably exist in numerical weather prediction （NWP） due to imperfect numeric and physical parameterizations. To eliminate these errors, by considering NWP as an inverse problem, an unknown term in the pred...Errors inevitably exist in numerical weather prediction （NWP） due to imperfect numeric and physical parameterizations. To eliminate these errors, by considering NWP as an inverse problem, an unknown term in the prediction equations can be estimated inversely by using the past data, which are presumed to represent the imperfection of the NWP model （model error, denoted as ME）. In this first paper of a two-part series, an iteration method for obtaining the MEs in past intervals is presented, and the results from testing its convergence in idealized experiments are reported. Moreover, two batches of iteration tests were applied in the global forecast system of the Global and Regional Assimilation and Prediction System （GRAPES-GFS） for July-August 2009 and January-February 2010. The datasets associated with the initial conditions and sea surface temperature （SST） were both based on NCEP （National Centers for Environmental Prediction） FNL （final） data. The results showed that 6th h forecast errors were reduced to 10% of their original value after a 20-step iteration. Then, off-line forecast error corrections were estimated linearly based on the 2-month mean MEs and compared with forecast errors. The estimated error corrections agreed well with the forecast errors, but the linear growth rate of the estimation was steeper than the forecast error. The advantage of this iteration method is that the MEs can provide the foundation for online correction. A larger proportion of the forecast errors can be expected to be canceled out by properly introducing the model error correction into GRAPES-GFS.展开更多
Model error is one of the key factors restricting the accuracy of numerical weather prediction （NWP）. Considering the continuous evolution of the atmosphere, the observed data （ignoring the measurement error） can ...Model error is one of the key factors restricting the accuracy of numerical weather prediction （NWP）. Considering the continuous evolution of the atmosphere, the observed data （ignoring the measurement error） can be viewed as a series of solutions of an accurate model governing the actual atmosphere. Model error is represented as an unknown term in the accurate model, thus NWP can be considered as an inverse problem to uncover the unknown error term. The inverse problem models can absorb long periods of observed data to generate model error correction procedures. They thus resolve the deficiency and faultiness of the NWP schemes employing only the initial-time data. In this study we construct two inverse problem models to estimate and extrapolate the time-varying and spatial-varying model errors in both the historical and forecast periods by using recent observations and analogue phenomena of the atmosphere. Numerical experiment on Burgers＇ equation has illustrated the substantial forecast improvement using inverse problem algorithms. The proposed inverse problem methods of suppressing NWP errors will be useful in future high accuracy applications of NWP.展开更多
This paper deals with a-posteriori error estimates for piecewise linear finite element approximations of parabolic problems in two space dimensions. The analysis extends previous results for elliptic problems to the p...This paper deals with a-posteriori error estimates for piecewise linear finite element approximations of parabolic problems in two space dimensions. The analysis extends previous results for elliptic problems to the parabolic context.展开更多
In this paper,we investigate a streamline diffusion finite element approxi- mation scheme for the constrained optimal control problem governed by linear con- vection dominated diffusion equations.We prove the existenc...In this paper,we investigate a streamline diffusion finite element approxi- mation scheme for the constrained optimal control problem governed by linear con- vection dominated diffusion equations.We prove the existence and uniqueness of the discretized scheme.Then a priori and a posteriori error estimates are derived for the state,the co-state and the control.Three numerical examples are presented to illustrate our theoretical results.展开更多
An online systematic error correction is presented and examined as a technique to improve the accuracy of real-time numerical weather prediction, based on the dataset of model errors （MEs） in past intervals. Given t...An online systematic error correction is presented and examined as a technique to improve the accuracy of real-time numerical weather prediction, based on the dataset of model errors （MEs） in past intervals. Given the analyses, the ME in each interval （6 h） between two analyses can be iteratively obtained by introducing an unknown tendency term into the prediction equation, shown in Part I of this two-paper series. In this part, after analyzing the 5-year （2001-2005） GRAPES- GFS （Global Forecast System of the Global and Regional Assimilation and Prediction System） error patterns and evolution, a systematic model error correction is given based on the least-squares approach by firstly using the past MEs. To test the correction, we applied the approach in GRAPES-GFS for July 2009 and January 2010. The datasets associated with the initial condition and SST used in this study were based on NCEP （National Centers for Environmental Prediction） FNL （final） data. The results indicated that the Northern Hemispheric systematically underestimated equator-to-pole geopotential gradient and westerly wind of GRAPES-GFS were largely enhanced, and the biases of temperature and wind in the tropics were strongly reduced. Therefore, the correction results in a more skillful forecast with lower mean bias and root-mean-square error and higher anomaly correlation coefficient.展开更多
In this paper, we consider the following problem:The quadratic spline collocation, with uniform mesh and the mid-knot points are taken as the collocation points for this problem is considered. With some assumptions, w...In this paper, we consider the following problem:The quadratic spline collocation, with uniform mesh and the mid-knot points are taken as the collocation points for this problem is considered. With some assumptions, we have proved that the solution of the quadratic spline collocation for the nonlinear problem can be written as a series expansions in integer powers of the mesh-size parameter. This gives us a construction method for using Richardson’s extrapolation. When we have a set of approximate solution with different mesh-size parameter a solution with high accuracy can he obtained by Richardson’s extrapolation.展开更多
In this paper, a posteriori error estimates were derived for piecewise linear finite element approximations to parabolic obstacle problems. The instrumental ingredient was introduced as a new interpolation operator wh...In this paper, a posteriori error estimates were derived for piecewise linear finite element approximations to parabolic obstacle problems. The instrumental ingredient was introduced as a new interpolation operator which has optimal approximation properties and preserves positivity. With the help of the interpolation operator the upper and lower bounds were obtained.展开更多
One of the classical approaches in the analysis of a variational inequality problem is to transform it into an equivalent optimization problem via the notion of gap function. The gap functions are useful tools in deri...One of the classical approaches in the analysis of a variational inequality problem is to transform it into an equivalent optimization problem via the notion of gap function. The gap functions are useful tools in deriving the error bounds which provide an estimated distance between a specific point and the exact solution of variational inequality problem. In this paper, we follow a similar approach for set-valued vector quasi variational inequality problems and define the gap functions based on scalarization scheme as well as the one with no scalar parameter. The error bounds results are obtained under fixed point symmetric and locally α-Holder assumptions on the set-valued map describing the domain of solution space of a set-valued vector quasi variational inequality problem.展开更多
In this paper,we consider the global error bound for the generalized complementarity problem(GCP)with analytic functions.Based on the new technique,we establish computable global error bound under milder conditions,wh...In this paper,we consider the global error bound for the generalized complementarity problem(GCP)with analytic functions.Based on the new technique,we establish computable global error bound under milder conditions,which refines the previously known results.展开更多
In this paper, we propose the nonconforming virtual element method (NCVEM) discretization for the pointwise control constraint optimal control problem governed by elliptic equations. Based on the NCVEM approximation o...In this paper, we propose the nonconforming virtual element method (NCVEM) discretization for the pointwise control constraint optimal control problem governed by elliptic equations. Based on the NCVEM approximation of state equation and the variational discretization of control variables, we construct a virtual element discrete scheme. For the state, adjoint state and control variable, we obtain the corresponding prior estimate in H<sup>1</sup> and L<sup>2</sup> norms. Finally, some numerical experiments are carried out to support the theoretical results.展开更多
In this article, a reduced mixed finite element （MFE） formulation based on proper orthogonal decomposition （POD） for the non-stationary conduction-convection problems is presented. Also the error estimates between...In this article, a reduced mixed finite element （MFE） formulation based on proper orthogonal decomposition （POD） for the non-stationary conduction-convection problems is presented. Also the error estimates between the reduced MFE solutions based on POD and usual MFE solutions are derived. It is shown by numerical examples that the results of numerical computation are consistent with theoretical conclusions. Moreover, it is shown that the reduced MFE formulation based on POD is feasible and efficient in finding numerical solutions for the non-stationary conduction-convection problems.展开更多
The initial value error and the imperfect numerical model are usually considered as error sources of numerical weather prediction （NWP）. By using past multi-time observations and model output, this study proposes a ...The initial value error and the imperfect numerical model are usually considered as error sources of numerical weather prediction （NWP）. By using past multi-time observations and model output, this study proposes a method to estimate imperfect numerical model error. This method can be inversely estimated through expressing the model error as a Lagrange interpolation polynomial, while the coefficients of polyno- mial are determined by past model performance. However, for practical application in the full NWP model, it is necessary to determine the following criteria： （1） the length of past data sufficient for estimation of the model errors, （2） a proper method of estimating the term ＂model integration with the exact solution＂ when solving the inverse problem, and （3） the extent to which this scheme is sensitive to the observational errors. In this study, such issues are resolved using a simple linear model, and an advection diffusion model is applied to discuss the sensitivity of the method to an artificial error source. The results indicate that the forecast errors can be largely reduced using the proposed method if the proper length of past data is chosen. To address the three problems, it is determined that （1） a few data limited by the order of the corrector can be used, （2） trapezoidal approximation can be employed to estimate the ＂term＂ in this study; however, a more accurate method should be explored for an operational NWP model, and （3） the correction is sensitive to observational error.展开更多
In this paper, we consider the Cauchy problem for the Laplace equation, which is severely ill-posed in the sense that the solution does not depend continuously on the data. A modified Tikhonov regularization method is...In this paper, we consider the Cauchy problem for the Laplace equation, which is severely ill-posed in the sense that the solution does not depend continuously on the data. A modified Tikhonov regularization method is proposed to solve this problem. An error estimate for the a priori parameter choice between the exact solution and its regularized approximation is obtained. Moreover, an a posteriori parameter choice rule is proposed and a stable error estimate is also obtained. Numerical examples illustrate the validity and effectiveness of this method.展开更多
In this paper, a high accuracy finite volume element method is presented for two-point boundary value problem of second order ordinary differential equation, which differs from the high order generalized difference me...In this paper, a high accuracy finite volume element method is presented for two-point boundary value problem of second order ordinary differential equation, which differs from the high order generalized difference methods. It is proved that the method has optimal order error estimate O(h3) in H1 norm. Finally, two examples show that the method is effective.展开更多
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.展开更多
基金Supported by the National Natural Science Foundation of China (Grant Nos.1236108412001130)。
文摘This paper is devoted to the Polynomial Preserving Recovery (PPR) based a posteriori error analysis for the second-order elliptic non-symmetric eigenvalue problem. An asymptotically exact a posteriori error estimator is proposed for solving the convection-dominated non-symmetric eigenvalue problem with non-smooth eigenfunctions or multiple eigenvalues. Numerical examples confirm our theoretical analysis.
基金supported by the National Natural Sci‐ence Foundation of China(Grant No.62306325)。
文摘During the use of robotics in applications such as antiterrorism or combat,a motion-constrained pursuer vehicle,such as a Dubins unmanned surface vehicle(USV),must get close enough(within a prescribed zero or positive distance)to a moving target as quickly as possible,resulting in the extended minimum-time intercept problem(EMTIP).Existing research has primarily focused on the zero-distance intercept problem,MTIP,establishing the necessary or sufficient conditions for MTIP optimality,and utilizing analytic algorithms,such as root-finding algorithms,to calculate the optimal solutions.However,these approaches depend heavily on the properties of the analytic algorithm,making them inapplicable when problem settings change,such as in the case of a positive effective range or complicated target motions outside uniform rectilinear motion.In this study,an approach employing a high-accuracy and quality-guaranteed mixed-integer piecewise-linear program(QG-PWL)is proposed for the EMTIP.This program can accommodate different effective interception ranges and complicated target motions(variable velocity or complicated trajectories).The high accuracy and quality guarantees of QG-PWL originate from elegant strategies such as piecewise linearization and other developed operation strategies.The approximate error in the intercept path length is proved to be bounded to h^(2)/(4√2),where h is the piecewise length.
基金funded by the National Natural Science Foundation Science Fund for Youth (Grant No.41405095)the Key Projects in the National Science and Technology Pillar Program during the Twelfth Fiveyear Plan Period (Grant No.2012BAC22B02)the National Natural Science Foundation Science Fund for Creative Research Groups (Grant No.41221064)
文摘Errors inevitably exist in numerical weather prediction （NWP） due to imperfect numeric and physical parameterizations. To eliminate these errors, by considering NWP as an inverse problem, an unknown term in the prediction equations can be estimated inversely by using the past data, which are presumed to represent the imperfection of the NWP model （model error, denoted as ME）. In this first paper of a two-part series, an iteration method for obtaining the MEs in past intervals is presented, and the results from testing its convergence in idealized experiments are reported. Moreover, two batches of iteration tests were applied in the global forecast system of the Global and Regional Assimilation and Prediction System （GRAPES-GFS） for July-August 2009 and January-February 2010. The datasets associated with the initial conditions and sea surface temperature （SST） were both based on NCEP （National Centers for Environmental Prediction） FNL （final） data. The results showed that 6th h forecast errors were reduced to 10% of their original value after a 20-step iteration. Then, off-line forecast error corrections were estimated linearly based on the 2-month mean MEs and compared with forecast errors. The estimated error corrections agreed well with the forecast errors, but the linear growth rate of the estimation was steeper than the forecast error. The advantage of this iteration method is that the MEs can provide the foundation for online correction. A larger proportion of the forecast errors can be expected to be canceled out by properly introducing the model error correction into GRAPES-GFS.
基金Project supported by the Special Scientific Research Project for Public Interest(Grant No.GYHY201206009)the Fundamental Research Funds for the Central Universities,China(Grant Nos.lzujbky-2012-13 and lzujbky-2013-11)the National Basic Research Program of China(Grant Nos.2012CB955902 and 2013CB430204)
文摘Model error is one of the key factors restricting the accuracy of numerical weather prediction （NWP）. Considering the continuous evolution of the atmosphere, the observed data （ignoring the measurement error） can be viewed as a series of solutions of an accurate model governing the actual atmosphere. Model error is represented as an unknown term in the accurate model, thus NWP can be considered as an inverse problem to uncover the unknown error term. The inverse problem models can absorb long periods of observed data to generate model error correction procedures. They thus resolve the deficiency and faultiness of the NWP schemes employing only the initial-time data. In this study we construct two inverse problem models to estimate and extrapolate the time-varying and spatial-varying model errors in both the historical and forecast periods by using recent observations and analogue phenomena of the atmosphere. Numerical experiment on Burgers＇ equation has illustrated the substantial forecast improvement using inverse problem algorithms. The proposed inverse problem methods of suppressing NWP errors will be useful in future high accuracy applications of NWP.
文摘This paper deals with a-posteriori error estimates for piecewise linear finite element approximations of parabolic problems in two space dimensions. The analysis extends previous results for elliptic problems to the parabolic context.
基金supported by the National Basic Research Program under the Grant 2005CB321701the National Natural Science Foundation of China under the Grants 60474027 and 10771211.
文摘In this paper,we investigate a streamline diffusion finite element approxi- mation scheme for the constrained optimal control problem governed by linear con- vection dominated diffusion equations.We prove the existence and uniqueness of the discretized scheme.Then a priori and a posteriori error estimates are derived for the state,the co-state and the control.Three numerical examples are presented to illustrate our theoretical results.
基金funded by the National Natural Science Foundation Science Fund for Youth (Grant No.41405095)the Key Projects in the National Science and Technology Pillar Program during the Twelfth Fiveyear Plan Period (Grant No.2012BAC22B02)the National Natural Science Foundation Science Fund for Creative Research Groups (Grant No.41221064)
文摘An online systematic error correction is presented and examined as a technique to improve the accuracy of real-time numerical weather prediction, based on the dataset of model errors （MEs） in past intervals. Given the analyses, the ME in each interval （6 h） between two analyses can be iteratively obtained by introducing an unknown tendency term into the prediction equation, shown in Part I of this two-paper series. In this part, after analyzing the 5-year （2001-2005） GRAPES- GFS （Global Forecast System of the Global and Regional Assimilation and Prediction System） error patterns and evolution, a systematic model error correction is given based on the least-squares approach by firstly using the past MEs. To test the correction, we applied the approach in GRAPES-GFS for July 2009 and January 2010. The datasets associated with the initial condition and SST used in this study were based on NCEP （National Centers for Environmental Prediction） FNL （final） data. The results indicated that the Northern Hemispheric systematically underestimated equator-to-pole geopotential gradient and westerly wind of GRAPES-GFS were largely enhanced, and the biases of temperature and wind in the tropics were strongly reduced. Therefore, the correction results in a more skillful forecast with lower mean bias and root-mean-square error and higher anomaly correlation coefficient.
基金The Project was supported by National Natural Science Foundation of China
文摘In this paper, we consider the following problem:The quadratic spline collocation, with uniform mesh and the mid-knot points are taken as the collocation points for this problem is considered. With some assumptions, we have proved that the solution of the quadratic spline collocation for the nonlinear problem can be written as a series expansions in integer powers of the mesh-size parameter. This gives us a construction method for using Richardson’s extrapolation. When we have a set of approximate solution with different mesh-size parameter a solution with high accuracy can he obtained by Richardson’s extrapolation.
基金Project supported by National Natural Science Foundation ofChina (Grant No .10471089)
文摘In this paper, a posteriori error estimates were derived for piecewise linear finite element approximations to parabolic obstacle problems. The instrumental ingredient was introduced as a new interpolation operator which has optimal approximation properties and preserves positivity. With the help of the interpolation operator the upper and lower bounds were obtained.
文摘One of the classical approaches in the analysis of a variational inequality problem is to transform it into an equivalent optimization problem via the notion of gap function. The gap functions are useful tools in deriving the error bounds which provide an estimated distance between a specific point and the exact solution of variational inequality problem. In this paper, we follow a similar approach for set-valued vector quasi variational inequality problems and define the gap functions based on scalarization scheme as well as the one with no scalar parameter. The error bounds results are obtained under fixed point symmetric and locally α-Holder assumptions on the set-valued map describing the domain of solution space of a set-valued vector quasi variational inequality problem.
基金supported by National Natural Science Foundation of China(Nos.11171180 and 11101303)Specialized Research Fund for the Doctoral Program of Chinese Higher Education(No.20113705110002)Shandong Provincial Natural Science Foundation(Nos.ZR2010AL005 and ZR2011FL017)
文摘In this paper,we consider the global error bound for the generalized complementarity problem(GCP)with analytic functions.Based on the new technique,we establish computable global error bound under milder conditions,which refines the previously known results.
文摘In this paper, we propose the nonconforming virtual element method (NCVEM) discretization for the pointwise control constraint optimal control problem governed by elliptic equations. Based on the NCVEM approximation of state equation and the variational discretization of control variables, we construct a virtual element discrete scheme. For the state, adjoint state and control variable, we obtain the corresponding prior estimate in H<sup>1</sup> and L<sup>2</sup> norms. Finally, some numerical experiments are carried out to support the theoretical results.
基金supported by the National Science Foundation of China (10871022 11061009+5 种基金 40821092)the National Basic Research Program (2010CB428403 2009CB421407 2010CB951001)Natural Science Foundation of Hebei Province (A2010001663)Chinese Universities Scientific Fund (2009-2-05)
文摘In this article, a reduced mixed finite element （MFE） formulation based on proper orthogonal decomposition （POD） for the non-stationary conduction-convection problems is presented. Also the error estimates between the reduced MFE solutions based on POD and usual MFE solutions are derived. It is shown by numerical examples that the results of numerical computation are consistent with theoretical conclusions. Moreover, it is shown that the reduced MFE formulation based on POD is feasible and efficient in finding numerical solutions for the non-stationary conduction-convection problems.
基金funded by the Special Scientific Research Project for Public Interest (GYHY201206009)the National Key Technologies Research and Development Program (Grant No. 2012BAC22B02)+2 种基金the National Natural Science Foundation Science Fund for Creative Research Groups (Grant No.41221064)the Special Scientific Research Project for Public Interest (Grant No. GYHY201006013)the National Natural Science Foundation of China (Grant No. 41105070 )
文摘The initial value error and the imperfect numerical model are usually considered as error sources of numerical weather prediction （NWP）. By using past multi-time observations and model output, this study proposes a method to estimate imperfect numerical model error. This method can be inversely estimated through expressing the model error as a Lagrange interpolation polynomial, while the coefficients of polyno- mial are determined by past model performance. However, for practical application in the full NWP model, it is necessary to determine the following criteria： （1） the length of past data sufficient for estimation of the model errors, （2） a proper method of estimating the term ＂model integration with the exact solution＂ when solving the inverse problem, and （3） the extent to which this scheme is sensitive to the observational errors. In this study, such issues are resolved using a simple linear model, and an advection diffusion model is applied to discuss the sensitivity of the method to an artificial error source. The results indicate that the forecast errors can be largely reduced using the proposed method if the proper length of past data is chosen. To address the three problems, it is determined that （1） a few data limited by the order of the corrector can be used, （2） trapezoidal approximation can be employed to estimate the ＂term＂ in this study; however, a more accurate method should be explored for an operational NWP model, and （3） the correction is sensitive to observational error.
基金supported by the National Natural Science Foundation of China(1117113611261032)+2 种基金the Distinguished Young Scholars Fund of Lan Zhou University of Technology(Q201015)the basic scientific research business expenses of Gansu province collegethe Natural Science Foundation of Gansu province(1310RJYA021)
文摘In this paper, we consider the Cauchy problem for the Laplace equation, which is severely ill-posed in the sense that the solution does not depend continuously on the data. A modified Tikhonov regularization method is proposed to solve this problem. An error estimate for the a priori parameter choice between the exact solution and its regularized approximation is obtained. Moreover, an a posteriori parameter choice rule is proposed and a stable error estimate is also obtained. Numerical examples illustrate the validity and effectiveness of this method.
基金heprojectissupportedbyNNSFofChina (No .1 9972 0 39) .
文摘In this paper, a high accuracy finite volume element method is presented for two-point boundary value problem of second order ordinary differential equation, which differs from the high order generalized difference methods. It is proved that the method has optimal order error estimate O(h3) in H1 norm. Finally, two examples show that the method is effective.
基金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.
