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 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.展开更多
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.展开更多
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 article, on the basis of two-level discretizations and multiscale finite element method, two kinds of finite element algorithms for steady Navier-Stokes problem are presented and discussed. The main technique ...In this article, on the basis of two-level discretizations and multiscale finite element method, two kinds of finite element algorithms for steady Navier-Stokes problem are presented and discussed. The main technique is first to use a standard finite element discretization on a coarse mesh to approximate low frequencies, then to apply the simple and Newton scheme to linearize discretizations on a fine grid. At this process, multiscale finite element method as a stabilized method deals with the lowest equal-order finite element pairs not satisfying the inf-sup condition. Under the uniqueness condition, error analyses for both algorithms are given. Numerical results are reported to demonstrate the effectiveness of the simple and Newton scheme.展开更多
Consider the finite volume element method for the thermal convection problem with the infinite Prandtl number. The author uses a conforming piecewise linear function on a fine triangulation for velocity and temperatur...Consider the finite volume element method for the thermal convection problem with the infinite Prandtl number. The author uses a conforming piecewise linear function on a fine triangulation for velocity and temperature, and a piecewise constant function on a coarse triangulation for pressure. For general triangulation the optimal order H1 norm error estimates are given.展开更多
In this paper, the Crank-Nicolson/Newton scheme for solving numerically second- order nonlinear parabolic problem is proposed. The standard Galerkin finite element method based on P2 conforming elements is used to the...In this paper, the Crank-Nicolson/Newton scheme for solving numerically second- order nonlinear parabolic problem is proposed. The standard Galerkin finite element method based on P2 conforming elements is used to the spatial discretization of the problem and the Crank-Nieolson/Newton scheme is applied to the time discretization of the resulted finite element equations. Moreover, assuming the appropriate regularity of the exact solution and the finite element solution, we obtain optimal error estimates of the fully discrete Crank- Nicolson/Newton scheme of nonlinear parabolic problem. Finally, numerical experiments are presented to show the efficient performance of the proposed scheme.展开更多
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.展开更多
We investigate numerically the dynamical reactions of a moving wheelset model to real measured track irregularities.The background is to examine whether the dynamics are suitable as the input to the inverse problem:de...We investigate numerically the dynamical reactions of a moving wheelset model to real measured track irregularities.The background is to examine whether the dynamics are suitable as the input to the inverse problem:determine the true track geometry from measured wheelset dynamical reactions.It is known that the method works well for the vertical position of the rails but the computed lateral position is often flawed.We find that the lateral motion of the wheelset often may differ from the track geometry.The cases are investigated closely but the reasons remain unknown.While the wheelset dynamics reflect the larger(>4-6 mm)aperiodic track disturbances and single large disturbances quite well,this does not seem to be the case for general smaller or periodic track irregularities or sections behind single large disturbances.The resulting dynamics of a wheelset to lateral track irregularities are in general not sufficiently accurate to be used as the basis for a description of the track irregularities.展开更多
基金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 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.
基金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.
基金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.
文摘In this article, on the basis of two-level discretizations and multiscale finite element method, two kinds of finite element algorithms for steady Navier-Stokes problem are presented and discussed. The main technique is first to use a standard finite element discretization on a coarse mesh to approximate low frequencies, then to apply the simple and Newton scheme to linearize discretizations on a fine grid. At this process, multiscale finite element method as a stabilized method deals with the lowest equal-order finite element pairs not satisfying the inf-sup condition. Under the uniqueness condition, error analyses for both algorithms are given. Numerical results are reported to demonstrate the effectiveness of the simple and Newton scheme.
基金Supported by National Natural Science Foundation of China(10071044)the Research Fund of Doctoral Program of High Education by State Education Ministry of China.
文摘Consider the finite volume element method for the thermal convection problem with the infinite Prandtl number. The author uses a conforming piecewise linear function on a fine triangulation for velocity and temperature, and a piecewise constant function on a coarse triangulation for pressure. For general triangulation the optimal order H1 norm error estimates are given.
基金in part supported by the Distinguished Young Scholars Fund of Xinjiang Province(2013711010)NCET-13-0988the NSF of China(11271313,11271298,61163027,and 11362021)
文摘In this paper, the Crank-Nicolson/Newton scheme for solving numerically second- order nonlinear parabolic problem is proposed. The standard Galerkin finite element method based on P2 conforming elements is used to the spatial discretization of the problem and the Crank-Nieolson/Newton scheme is applied to the time discretization of the resulted finite element equations. Moreover, assuming the appropriate regularity of the exact solution and the finite element solution, we obtain optimal error estimates of the fully discrete Crank- Nicolson/Newton scheme of nonlinear parabolic problem. Finally, numerical experiments are presented to show the efficient performance of the proposed scheme.
基金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.
文摘We investigate numerically the dynamical reactions of a moving wheelset model to real measured track irregularities.The background is to examine whether the dynamics are suitable as the input to the inverse problem:determine the true track geometry from measured wheelset dynamical reactions.It is known that the method works well for the vertical position of the rails but the computed lateral position is often flawed.We find that the lateral motion of the wheelset often may differ from the track geometry.The cases are investigated closely but the reasons remain unknown.While the wheelset dynamics reflect the larger(>4-6 mm)aperiodic track disturbances and single large disturbances quite well,this does not seem to be the case for general smaller or periodic track irregularities or sections behind single large disturbances.The resulting dynamics of a wheelset to lateral track irregularities are in general not sufficiently accurate to be used as the basis for a description of the track irregularities.
