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.展开更多
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.展开更多
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 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 conditio...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, 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.展开更多
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.展开更多
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.展开更多
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 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.
基金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.
基金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.
基金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.
基金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.
基金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.
文摘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.
基金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.
