Theoretical aspects of variational data assimilation (VDA) for a simple model with both global and local observational data are discussed. For the VDA problems with global observational data, the initial conditions ...Theoretical aspects of variational data assimilation (VDA) for a simple model with both global and local observational data are discussed. For the VDA problems with global observational data, the initial conditions and parameters for the model are revisited and the model itself is modified. The estimates of both error and convergence rate are theoretically made and the vahdity of the method is proved. For VDA problem with local observation data, the conventional VDA method are out of use due to the ill-posedness of the problem. In order to overcome the difficulties caused by the ill-posedness, the initial conditions and parameters of the model are modified by using the improved VDA method, and the estimates of both error and convergence rate are also made. Finally, the validity of the improved VDA method is proved through theoretical analysis and illustrated with an example, and a theoretical criterion of the regularization parameters is proposed.展开更多
In this paper we propose and analyze a second order accurate numericalscheme for the Cahn-Hilliard equation with logarithmic Flory Huggins energy potential. A modified Crank-Nicolson approximation is applied to the l...In this paper we propose and analyze a second order accurate numericalscheme for the Cahn-Hilliard equation with logarithmic Flory Huggins energy potential. A modified Crank-Nicolson approximation is applied to the logarithmic nonlinear term, while the expansive term is updated by an explicit second order AdamsBashforth extrapolation, and an alternate temporal stencil is used for the surface diffusion term. A nonlinear artificial regularization term is added in the numerical scheme,which ensures the positivity-preserving property, i.e., the numerical value of the phasevariable is always between -1 and 1 at a point-wise level. Furthermore, an unconditional energy stability of the numerical scheme is derived, leveraging the special formof the logarithmic approximation term. In addition, an optimal rate convergence estimate is provided for the proposed numerical scheme, with the help of linearizedstability analysis. A few numerical results, including both the constant-mobility andsolution-dependent mobility flows, are presented to validate the robustness of the proposed numerical scheme.展开更多
We propose a nonparametric change point estimator in the distributions of a sequence of independent observations in terms of the test statistics given by Huˇskov′a and Meintanis(2006) that are based on weighted empi...We propose a nonparametric change point estimator in the distributions of a sequence of independent observations in terms of the test statistics given by Huˇskov′a and Meintanis(2006) that are based on weighted empirical characteristic functions. The weight function ω(t; a) under consideration includes the two weight functions from Huˇskov′a and Meintanis(2006) plus the weight function used by Matteson and James(2014),where a is a tuning parameter. Under the local alternative hypothesis, we establish the consistency, convergence rate, and asymptotic distribution of this change point estimator which is the maxima of a two-side Brownian motion with a drift. Since the performance of the change point estimator depends on a in use, we thus propose an algorithm for choosing an appropriate value of a, denoted by a_s which is also justified. Our simulation study shows that the change point estimate obtained by using a_s has a satisfactory performance. We also apply our method to a real dataset.展开更多
基金Project supported by the National Natural Science Foundation of China(Nos.40675020, 50505005)
文摘Theoretical aspects of variational data assimilation (VDA) for a simple model with both global and local observational data are discussed. For the VDA problems with global observational data, the initial conditions and parameters for the model are revisited and the model itself is modified. The estimates of both error and convergence rate are theoretically made and the vahdity of the method is proved. For VDA problem with local observation data, the conventional VDA method are out of use due to the ill-posedness of the problem. In order to overcome the difficulties caused by the ill-posedness, the initial conditions and parameters of the model are modified by using the improved VDA method, and the estimates of both error and convergence rate are also made. Finally, the validity of the improved VDA method is proved through theoretical analysis and illustrated with an example, and a theoretical criterion of the regularization parameters is proposed.
基金This work is supported in part by the grants NSFC 12071090(W.Chen)NSF DMS-2012669(C.Wang)+2 种基金NSFC 11871159Guangdong Provincial Key Laboratory for Computational Science and Material Design 2019B030301001(X.Wang)NSF DMS-1719854,DMS-2012634(S.Wise).C.Wang also thanks the Key Laboratory of Mathematics for Nonlinear Sciences,Fudan University,for the support.
文摘In this paper we propose and analyze a second order accurate numericalscheme for the Cahn-Hilliard equation with logarithmic Flory Huggins energy potential. A modified Crank-Nicolson approximation is applied to the logarithmic nonlinear term, while the expansive term is updated by an explicit second order AdamsBashforth extrapolation, and an alternate temporal stencil is used for the surface diffusion term. A nonlinear artificial regularization term is added in the numerical scheme,which ensures the positivity-preserving property, i.e., the numerical value of the phasevariable is always between -1 and 1 at a point-wise level. Furthermore, an unconditional energy stability of the numerical scheme is derived, leveraging the special formof the logarithmic approximation term. In addition, an optimal rate convergence estimate is provided for the proposed numerical scheme, with the help of linearizedstability analysis. A few numerical results, including both the constant-mobility andsolution-dependent mobility flows, are presented to validate the robustness of the proposed numerical scheme.
基金supported by Natural Sciences and the Engineering Research Council of Canada (Grant No. 105557-2012)National Natural Science Foundation for Young Scientists of China (Grant No. 11201108)+1 种基金the National Statistical Research Plan Project (Grant No. 2012LZ009)the Humanities and Social Sciences Project from Ministry of Education of China (Grant No. 12YJC910007)
文摘We propose a nonparametric change point estimator in the distributions of a sequence of independent observations in terms of the test statistics given by Huˇskov′a and Meintanis(2006) that are based on weighted empirical characteristic functions. The weight function ω(t; a) under consideration includes the two weight functions from Huˇskov′a and Meintanis(2006) plus the weight function used by Matteson and James(2014),where a is a tuning parameter. Under the local alternative hypothesis, we establish the consistency, convergence rate, and asymptotic distribution of this change point estimator which is the maxima of a two-side Brownian motion with a drift. Since the performance of the change point estimator depends on a in use, we thus propose an algorithm for choosing an appropriate value of a, denoted by a_s which is also justified. Our simulation study shows that the change point estimate obtained by using a_s has a satisfactory performance. We also apply our method to a real dataset.