This paper presents a new method to estimate the height of the atmospheric boundary layer(ABL) by using COSMIC radio occultation bending angle(BA) data. Using the numerical differentiation method combined with the reg...This paper presents a new method to estimate the height of the atmospheric boundary layer(ABL) by using COSMIC radio occultation bending angle(BA) data. Using the numerical differentiation method combined with the regularization technique, the first derivative of BA profiles is retrieved, and the height at which the first derivative of BA has the global minimum is defined to be the ABL height. To reflect the reliability of estimated ABL heights, the sharpness parameter is introduced, according to the relative minimum of the BA derivative. Then, it is applied to four months of COSMIC BA data(January, April, July, and October in 2008), and the ABL heights estimated are compared with two kinds of ABL heights from COSMIC products and with the heights determined by the finite difference method upon the refractivity data. For sharp ABL tops(large sharpness parameters), there is little difference between the ABL heights determined by different methods, i.e.,the uncertainties are small; whereas, for non-sharp ABL tops(small sharpness parameters), big differences exist in the ABL heights obtained by different methods, which means large uncertainties for different methods. In addition, the new method can detect thin ABLs and provide a reference ABL height in the cases eliminated by other methods. Thus, the application of the numerical differentiation method combined with the regularization technique to COSMIC BA data is an appropriate choice and has further application value.展开更多
A fast explicit finite difference method (FEFDM),derived from the differential equations of one-dimensional steady pipe flow,was presented for calculation of wellhead injection pressure.Recalculation with a traditiona...A fast explicit finite difference method (FEFDM),derived from the differential equations of one-dimensional steady pipe flow,was presented for calculation of wellhead injection pressure.Recalculation with a traditional numerical method of the same equations corroborates well the reliability and rate of FEFDM.Moreover,a flow rate estimate method was developed for the project whose injection rate has not been clearly determined.A wellhead pressure regime determined by this method was successfully applied to the trial injection operations in Shihezi formation of Shenhua CCS Project,which is a good practice verification of FEFDM.At last,this method was used to evaluate the effect of friction and acceleration terms on the flow equation on the wellhead pressure.The result shows that for deep wellbore,the friction term can be omitted when flow rate is low and in a wide range of velocity the acceleration term can always be deleted.It is also shown that with flow rate increasing,the friction term can no longer be neglected.展开更多
In this paper, the pseudo-spectral approximations for a class of the Kdv-Burgers type equation is presented. Convergence and stability of the approximation have been proved by Sobolev's inequalities and the bounde...In this paper, the pseudo-spectral approximations for a class of the Kdv-Burgers type equation is presented. Convergence and stability of the approximation have been proved by Sobolev's inequalities and the bounded extensive method of the nonlinear function. Finally, the numerical examples are proposed.展开更多
Mathematics is very important for the engineering and scientist but to make understand the mathematics is very difficult if without proper tools and suitable measurement. A numerical method is one of the algorithms wh...Mathematics is very important for the engineering and scientist but to make understand the mathematics is very difficult if without proper tools and suitable measurement. A numerical method is one of the algorithms which involved with computer programming. In this paper, Scilab is used to carter the problems related the mathematical models such as Matrices, operation with ODE＇s and solving the Integration.展开更多
In order to simulate the airflow in anhydrous case and the water-air flow in groundwater case, a numerical model of airflow in soil was developed. For the nonlinearity of the governing partial differential equation, t...In order to simulate the airflow in anhydrous case and the water-air flow in groundwater case, a numerical model of airflow in soil was developed. For the nonlinearity of the governing partial differential equation, the corresponding discretization and linearization methods were given. Due to the mass transfer between air-phase and water-phase, phase states of the model elements were constantly changing. Thus, parameters of the model were divided into primary ones and secondary ones, and the primary variables changing with phase states and the secondary variables can be obtained by their functional relationship with the primary variables. Additionally, the special definite condition of this numerical model was illustrated. Two examples were given to simulate the airflow in soil whether there was groundwater or not, and the effectiveness of the numerical model is verified by comparing the results of simulation with that of exoeriment.展开更多
The artificial boundary method is one of the most popular and effective numerical methods tor solving partial differential equations on unbounded domains, with more than thirty years development. The artificiM boundar...The artificial boundary method is one of the most popular and effective numerical methods tor solving partial differential equations on unbounded domains, with more than thirty years development. The artificiM boundary method has reached maturity in recent years. It has been applied to various problems in scientific and engineering computations, and the theoretical issues such as the convergence and error estimates of the artificial boundary method have been solved gradually. This paper reviews the development and discusses different forms of the artificial boundary method.展开更多
A robust nonlinear control method is presented for spacecraft precise formation flying.With the constraint forces and consid-ering nonlinearity and perturbations,the problem of the formation keeping is changed to the ...A robust nonlinear control method is presented for spacecraft precise formation flying.With the constraint forces and consid-ering nonlinearity and perturbations,the problem of the formation keeping is changed to the Lagrange systems with the holonomic constraints and the differential algebraic equations (DAE).The nonlinear control laws are developed by solving the DAE.Because the traditional numerical solving methods of DAE are very sensitive to the various errors and the resulting con-trol laws are not robust in engineering application,the robust control law designed method is further developed by designing the correct coefficients to correct the errors of the formation array constraints.A numeral study simulated the robustness of this method for the various errors in the formation flying mission,including the initial errors of spacecraft formation,the reference satellite orbit determination errors,the relative perturbation forces model errors,and so on.展开更多
The main goal of this article is to discuss the numerical solution to a nonlinear wave equation associated with the first of the celebrated Painlevd transcendent ordinary differential equations. In order to solve nume...The main goal of this article is to discuss the numerical solution to a nonlinear wave equation associated with the first of the celebrated Painlevd transcendent ordinary differential equations. In order to solve numerically the above equation, whose solutions blow up in finite time, the authors advocate a numerical methodology based on the Strang＇s symmetrized operator-splitting scheme. With this approach, one can decouple nonlinearity and differential operators, leading to the alternate solution at every time step of the equation as follows： （i） The first Painlevd ordinary differential equation, （ii） a linear wave equation with a constant coefficient. Assuming that the space dimension is two, the authors consider a fully discrete variant of the above scheme, where the space-time discretization of the linear wave equation sub-steps is achieved via a Galerkin/finite element space approximation combined with a second order accurate centered time discretization scheme. To handle the nonlinear sub-steps, a second order accurate centered explicit time discretization scheme with adaptively variable time step is used, in order to follow accurately the fast dynamic of the solution before it blows up. The results of numerical experiments are presented for different coefficients and boundary conditions. They show that the above methodology is robust and describes fairly accurately the evolution of a rather ＂violent＂ phenomenon.展开更多
基金supported by the National Natural Science Foundation of China (Grant No. 41475021)
文摘This paper presents a new method to estimate the height of the atmospheric boundary layer(ABL) by using COSMIC radio occultation bending angle(BA) data. Using the numerical differentiation method combined with the regularization technique, the first derivative of BA profiles is retrieved, and the height at which the first derivative of BA has the global minimum is defined to be the ABL height. To reflect the reliability of estimated ABL heights, the sharpness parameter is introduced, according to the relative minimum of the BA derivative. Then, it is applied to four months of COSMIC BA data(January, April, July, and October in 2008), and the ABL heights estimated are compared with two kinds of ABL heights from COSMIC products and with the heights determined by the finite difference method upon the refractivity data. For sharp ABL tops(large sharpness parameters), there is little difference between the ABL heights determined by different methods, i.e.,the uncertainties are small; whereas, for non-sharp ABL tops(small sharpness parameters), big differences exist in the ABL heights obtained by different methods, which means large uncertainties for different methods. In addition, the new method can detect thin ABLs and provide a reference ABL height in the cases eliminated by other methods. Thus, the application of the numerical differentiation method combined with the regularization technique to COSMIC BA data is an appropriate choice and has further application value.
基金Project(Z110803)supported by the State Key Laboratory of Geomechanics and Geotechnical Engineering,ChinaProject(2008AA062303)supported by the National High Technology Research and Development Program of China
文摘A fast explicit finite difference method (FEFDM),derived from the differential equations of one-dimensional steady pipe flow,was presented for calculation of wellhead injection pressure.Recalculation with a traditional numerical method of the same equations corroborates well the reliability and rate of FEFDM.Moreover,a flow rate estimate method was developed for the project whose injection rate has not been clearly determined.A wellhead pressure regime determined by this method was successfully applied to the trial injection operations in Shihezi formation of Shenhua CCS Project,which is a good practice verification of FEFDM.At last,this method was used to evaluate the effect of friction and acceleration terms on the flow equation on the wellhead pressure.The result shows that for deep wellbore,the friction term can be omitted when flow rate is low and in a wide range of velocity the acceleration term can always be deleted.It is also shown that with flow rate increasing,the friction term can no longer be neglected.
基金Supported by the Natural Science Foundation of Henan Educational Committee(2003110005)Supported by the Natural Science Foundation of Henan University(XK02069)
文摘In this paper, the pseudo-spectral approximations for a class of the Kdv-Burgers type equation is presented. Convergence and stability of the approximation have been proved by Sobolev's inequalities and the bounded extensive method of the nonlinear function. Finally, the numerical examples are proposed.
文摘Mathematics is very important for the engineering and scientist but to make understand the mathematics is very difficult if without proper tools and suitable measurement. A numerical method is one of the algorithms which involved with computer programming. In this paper, Scilab is used to carter the problems related the mathematical models such as Matrices, operation with ODE＇s and solving the Integration.
基金Project(Y5080022) supported by the Natural Science Foundation of Zhejiang Province,ChinaProject(RC1202) supported by Scientific and Technological Program of Water Resources Department of Zhejiang Province in 2012,ChinaProject(Y201224384) supported by Scientific Research Program of Education Department of Zhejiang Province in 2012,China
文摘In order to simulate the airflow in anhydrous case and the water-air flow in groundwater case, a numerical model of airflow in soil was developed. For the nonlinearity of the governing partial differential equation, the corresponding discretization and linearization methods were given. Due to the mass transfer between air-phase and water-phase, phase states of the model elements were constantly changing. Thus, parameters of the model were divided into primary ones and secondary ones, and the primary variables changing with phase states and the secondary variables can be obtained by their functional relationship with the primary variables. Additionally, the special definite condition of this numerical model was illustrated. Two examples were given to simulate the airflow in soil whether there was groundwater or not, and the effectiveness of the numerical model is verified by comparing the results of simulation with that of exoeriment.
基金supported by National National Science Foundation of China(Grant No.10971116)FRG of Hong Kong Baptist University(Grant No.FRG1/11-12/051)
文摘The artificial boundary method is one of the most popular and effective numerical methods tor solving partial differential equations on unbounded domains, with more than thirty years development. The artificiM boundary method has reached maturity in recent years. It has been applied to various problems in scientific and engineering computations, and the theoretical issues such as the convergence and error estimates of the artificial boundary method have been solved gradually. This paper reviews the development and discusses different forms of the artificial boundary method.
基金supported by the China Postdoctoral Foundation (Grant Nos. 20080440217, 200902666)
文摘A robust nonlinear control method is presented for spacecraft precise formation flying.With the constraint forces and consid-ering nonlinearity and perturbations,the problem of the formation keeping is changed to the Lagrange systems with the holonomic constraints and the differential algebraic equations (DAE).The nonlinear control laws are developed by solving the DAE.Because the traditional numerical solving methods of DAE are very sensitive to the various errors and the resulting con-trol laws are not robust in engineering application,the robust control law designed method is further developed by designing the correct coefficients to correct the errors of the formation array constraints.A numeral study simulated the robustness of this method for the various errors in the formation flying mission,including the initial errors of spacecraft formation,the reference satellite orbit determination errors,the relative perturbation forces model errors,and so on.
文摘The main goal of this article is to discuss the numerical solution to a nonlinear wave equation associated with the first of the celebrated Painlevd transcendent ordinary differential equations. In order to solve numerically the above equation, whose solutions blow up in finite time, the authors advocate a numerical methodology based on the Strang＇s symmetrized operator-splitting scheme. With this approach, one can decouple nonlinearity and differential operators, leading to the alternate solution at every time step of the equation as follows： （i） The first Painlevd ordinary differential equation, （ii） a linear wave equation with a constant coefficient. Assuming that the space dimension is two, the authors consider a fully discrete variant of the above scheme, where the space-time discretization of the linear wave equation sub-steps is achieved via a Galerkin/finite element space approximation combined with a second order accurate centered time discretization scheme. To handle the nonlinear sub-steps, a second order accurate centered explicit time discretization scheme with adaptively variable time step is used, in order to follow accurately the fast dynamic of the solution before it blows up. The results of numerical experiments are presented for different coefficients and boundary conditions. They show that the above methodology is robust and describes fairly accurately the evolution of a rather ＂violent＂ phenomenon.
