In this work,we focus on the conforming and nonconforming leap-frog virtual element methods for the generalized nonlinear Schrodinger equation,and establish their unconditional stability and optimal error estimates.By...In this work,we focus on the conforming and nonconforming leap-frog virtual element methods for the generalized nonlinear Schrodinger equation,and establish their unconditional stability and optimal error estimates.By constructing a time-discrete system,the error between the solutions of the continuous model and the numerical scheme is separated into the temporal error and the spatial error,which makes the spatial error τ-independent.The inverse inequalities in the existing conforming and new constructed nonconforming virtual element spaces are utilized to derive the L^(∞)-norm uniform boundedness of numerical solutions without any restrictions on time-space step ratio,and then unconditionally optimal error estimates of the numerical schemes are obtained naturally.What needs to be emphasized is that if we use the pre-existing nonconforming virtual elements,there is no way to derive the L^(∞)-norm uniform boundedness of the functions in the nonconforming virtual element spaces so as to be hard to get the corresponding inverse inequalities.Finally,several numerical examples are reported to confirm our theoretical results.展开更多
With the observational wind data and the Zebiak-Cane model, the impact of Madden-Iulian Oscillation (MJO) as external forcing on El Nino-Southern Oscillation (ENSO) predictability is studied. The observational dat...With the observational wind data and the Zebiak-Cane model, the impact of Madden-Iulian Oscillation (MJO) as external forcing on El Nino-Southern Oscillation (ENSO) predictability is studied. The observational data are analyzed with Continuous Wavelet Transform (CWT) and then used to extract MJO signals, which are added into the model to get a new model. After the Conditional Nonlinear Optimal Perturbation (CNOP) method has been used, the initial errors which can evolve into maximum prediction error, model errors and their join errors are gained and then the Nifio 3 indices and spatial structures of three kinds of errors are investigated. The results mainly show that the observational MJO has little impact on the maximum prediction error of ENSO events and the initial error affects much greater than model error caused by MJO forcing. These demonstrate that the initial error might be the main error source that produces uncertainty in ENSO prediction, which could provide a theoretical foundation for the adaptive data assimilation of the ENSO forecast and contribute to the ENSO target observation.展开更多
To improve the spectral characteristics of the high-order weighted compact nonlinear scheme(WCNS),optimized flux difference schemes are proposed.The disadvantages in previous optimization routines,i.e.,reducing formal...To improve the spectral characteristics of the high-order weighted compact nonlinear scheme(WCNS),optimized flux difference schemes are proposed.The disadvantages in previous optimization routines,i.e.,reducing formal orders,or extending stencil widths,are avoided in the new optimized schemes by utilizing fluxes from both cell-edges and cell-nodes.Optimizations are implemented with Fourier analysis for linear schemes and the approximate dispersion relation(ADR)for nonlinear schemes.Classical difference schemes are restored near discontinuities to suppress numerical oscillations with use of a shock sensor based on smoothness indicators.The results of several benchmark numerical tests indicate that the new optimized difference schemes outperform the classical schemes,in terms of accuracy and resolution for smooth wave and vortex,especially for long-time simulations.Using optimized schemes increases the total CPU time by less than 4%.展开更多
A new mixed scheme which combines the variation of constants and the H1-Galerkin mixed finite element method is constructed for nonlinear Sobolev equation with nonlinear con- vection term. Optimal error estimates are ...A new mixed scheme which combines the variation of constants and the H1-Galerkin mixed finite element method is constructed for nonlinear Sobolev equation with nonlinear con- vection term. Optimal error estimates are derived for both semidiscrete and fully discrete schemes. Finally, some numerical results are given to confirm the theoretical analysis of the proposed method.展开更多
In this paper, a computational approach is proposed for solving the discrete-time nonlinear optimal control problem, which is disturbed by a sequence of random noises. Because of the exact solution of such optimal con...In this paper, a computational approach is proposed for solving the discrete-time nonlinear optimal control problem, which is disturbed by a sequence of random noises. Because of the exact solution of such optimal control problem is impossible to be obtained, estimating the state dynamics is currently required. Here, it is assumed that the output can be measured from the real plant process. In our approach, the state mean propagation is applied in order to construct a linear model-based optimal control problem, where the model output is measureable. On this basis, an output error, which takes into account the differences between the real output and the model output, is defined. Then, this output error is minimized by applying the stochastic approximation approach. During the computation procedure, the stochastic gradient is established, so as the optimal solution of the model used can be updated iteratively. Once the convergence is achieved, the iterative solution approximates to the true optimal solution of the original optimal control problem, in spite of model-reality differences. For illustration, an example on a continuous stirred-tank reactor problem is studied, and the result obtained shows the applicability of the approach proposed. Hence, the efficiency of the approach proposed is highly recommended.展开更多
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.展开更多
Output measurement for nonlinear optimal control problems is an interesting issue. Because the structure of the real plant is complex, the output channel could give a significant response corresponding to the real pla...Output measurement for nonlinear optimal control problems is an interesting issue. Because the structure of the real plant is complex, the output channel could give a significant response corresponding to the real plant. In this paper, a least squares scheme, which is based on the Gauss-Newton algorithm, is proposed. The aim is to approximate the output that is measured from the real plant. In doing so, an appropriate output measurement from the model used is suggested. During the computation procedure, the control trajectory is updated iteratively by using the Gauss-Newton recursion scheme. Consequently, the output residual between the original output and the suggested output is minimized. Here, the linear model-based optimal control model is considered, so as the optimal control law is constructed. By feed backing the updated control trajectory into the dynamic system, the iterative solution of the model used could approximate to the correct optimal solution of the original optimal control problem, in spite of model-reality differences. For illustration, current converted and isothermal reaction rector problems are studied and the results are demonstrated. In conclusion, the efficiency of the approach proposed is highly presented.展开更多
In the numerical prediction of weather or climate events,the uncertainty of the initial values and/or prediction models can bring the forecast result’s uncertainty.Due to the absence of true states,studies on this pr...In the numerical prediction of weather or climate events,the uncertainty of the initial values and/or prediction models can bring the forecast result’s uncertainty.Due to the absence of true states,studies on this problem mainly focus on the three subproblems of predictability,i.e.,the lower bound of the maximum predictable time,the upper bound of the prediction error,and the lower bound of the maximum allowable initial error.Aimed at the problem of the lower bound estimation of the maximum allowable initial error,this study first illustrates the shortcoming of the existing estimation,and then presents a new estimation based on the initial observation precision and proves it theoretically.Furthermore,the new lower bound estimations of both the two-dimensional ikeda model and lorenz96 model are obtained by using the cnop(conditional nonlinear optimal perturbation)method and a pso(particle swarm optimization)algorithm,and the estimated precisions are also analyzed.Besides,the estimations yielded by the existing and new formulas are compared;the results show that the estimations produced by the existing formula are often incorrect.展开更多
In this paper,we consider the local discontinuous Galerkin method with generalized alter-nating numerical fluxes for two-dimensional nonlinear Schrödinger equations on Carte-sian meshes.The generalized fluxes not...In this paper,we consider the local discontinuous Galerkin method with generalized alter-nating numerical fluxes for two-dimensional nonlinear Schrödinger equations on Carte-sian meshes.The generalized fluxes not only lead to a smaller magnitude of the errors,but can guarantee an energy conservative property that is useful for long time simulations in resolving waves.By virtue of generalized skew-symmetry property of the discontinuous Galerkin spatial operators,two energy equations are established and stability results con-taining energy conservation of the prime variable as well as auxiliary variables are shown.To derive optimal error estimates for nonlinear Schrödinger equations,an additional energy equation is constructed and two a priori error assumptions are used.This,together with properties of some generalized Gauss-Radau projections and a suitable numerical initial condition,implies optimal order of k+1.Numerical experiments are given to demonstrate the theoretical results.展开更多
基金supported by the NSF of China(Grant Nos.11801527,11701522,11771163,12011530058,11671160,1191101330)by the China Postdoctoral Science Foundation(Grant Nos.2018M632791,2019M662506).
文摘In this work,we focus on the conforming and nonconforming leap-frog virtual element methods for the generalized nonlinear Schrodinger equation,and establish their unconditional stability and optimal error estimates.By constructing a time-discrete system,the error between the solutions of the continuous model and the numerical scheme is separated into the temporal error and the spatial error,which makes the spatial error τ-independent.The inverse inequalities in the existing conforming and new constructed nonconforming virtual element spaces are utilized to derive the L^(∞)-norm uniform boundedness of numerical solutions without any restrictions on time-space step ratio,and then unconditionally optimal error estimates of the numerical schemes are obtained naturally.What needs to be emphasized is that if we use the pre-existing nonconforming virtual elements,there is no way to derive the L^(∞)-norm uniform boundedness of the functions in the nonconforming virtual element spaces so as to be hard to get the corresponding inverse inequalities.Finally,several numerical examples are reported to confirm our theoretical results.
基金The National Natural Science Foundation of China under contract No.41405062
文摘With the observational wind data and the Zebiak-Cane model, the impact of Madden-Iulian Oscillation (MJO) as external forcing on El Nino-Southern Oscillation (ENSO) predictability is studied. The observational data are analyzed with Continuous Wavelet Transform (CWT) and then used to extract MJO signals, which are added into the model to get a new model. After the Conditional Nonlinear Optimal Perturbation (CNOP) method has been used, the initial errors which can evolve into maximum prediction error, model errors and their join errors are gained and then the Nifio 3 indices and spatial structures of three kinds of errors are investigated. The results mainly show that the observational MJO has little impact on the maximum prediction error of ENSO events and the initial error affects much greater than model error caused by MJO forcing. These demonstrate that the initial error might be the main error source that produces uncertainty in ENSO prediction, which could provide a theoretical foundation for the adaptive data assimilation of the ENSO forecast and contribute to the ENSO target observation.
基金Project supported by the National Key Project(No.GJXM92579)the Defense Industrial Technology Development Program(No.C1520110002)the State Administration of Science,Technology and Industry for National Defence,China。
文摘To improve the spectral characteristics of the high-order weighted compact nonlinear scheme(WCNS),optimized flux difference schemes are proposed.The disadvantages in previous optimization routines,i.e.,reducing formal orders,or extending stencil widths,are avoided in the new optimized schemes by utilizing fluxes from both cell-edges and cell-nodes.Optimizations are implemented with Fourier analysis for linear schemes and the approximate dispersion relation(ADR)for nonlinear schemes.Classical difference schemes are restored near discontinuities to suppress numerical oscillations with use of a shock sensor based on smoothness indicators.The results of several benchmark numerical tests indicate that the new optimized difference schemes outperform the classical schemes,in terms of accuracy and resolution for smooth wave and vortex,especially for long-time simulations.Using optimized schemes increases the total CPU time by less than 4%.
基金Supported by National Natural Science Fund of China (11061021)Key Project of Chinese Ministry of Education (12024)+2 种基金Natural Science Fund of Inner Mongolia Autonomous Region (2012MS0108,2012MS0106,2011BS0102)Scientific Research Projection of Higher Schools of Inner Mongolia (NJZZ12011,NJZY13199)Program of Higher-level talents of Inner Mongolia University (125119,Z200901004,30105-125132)
文摘A new mixed scheme which combines the variation of constants and the H1-Galerkin mixed finite element method is constructed for nonlinear Sobolev equation with nonlinear con- vection term. Optimal error estimates are derived for both semidiscrete and fully discrete schemes. Finally, some numerical results are given to confirm the theoretical analysis of the proposed method.
文摘In this paper, a computational approach is proposed for solving the discrete-time nonlinear optimal control problem, which is disturbed by a sequence of random noises. Because of the exact solution of such optimal control problem is impossible to be obtained, estimating the state dynamics is currently required. Here, it is assumed that the output can be measured from the real plant process. In our approach, the state mean propagation is applied in order to construct a linear model-based optimal control problem, where the model output is measureable. On this basis, an output error, which takes into account the differences between the real output and the model output, is defined. Then, this output error is minimized by applying the stochastic approximation approach. During the computation procedure, the stochastic gradient is established, so as the optimal solution of the model used can be updated iteratively. Once the convergence is achieved, the iterative solution approximates to the true optimal solution of the original optimal control problem, in spite of model-reality differences. For illustration, an example on a continuous stirred-tank reactor problem is studied, and the result obtained shows the applicability of the approach proposed. Hence, the efficiency of the approach proposed is highly recommended.
基金Supported by National Natural Science Foundation of China (61135001, 61075029, 61074179, 61074155) and the Postdoctoral Science Foundation of China (20110491692)
基金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.
文摘Output measurement for nonlinear optimal control problems is an interesting issue. Because the structure of the real plant is complex, the output channel could give a significant response corresponding to the real plant. In this paper, a least squares scheme, which is based on the Gauss-Newton algorithm, is proposed. The aim is to approximate the output that is measured from the real plant. In doing so, an appropriate output measurement from the model used is suggested. During the computation procedure, the control trajectory is updated iteratively by using the Gauss-Newton recursion scheme. Consequently, the output residual between the original output and the suggested output is minimized. Here, the linear model-based optimal control model is considered, so as the optimal control law is constructed. By feed backing the updated control trajectory into the dynamic system, the iterative solution of the model used could approximate to the correct optimal solution of the original optimal control problem, in spite of model-reality differences. For illustration, current converted and isothermal reaction rector problems are studied and the results are demonstrated. In conclusion, the efficiency of the approach proposed is highly presented.
基金supported by the National Natural Science Foundation of China(Grant No.41331174)
文摘In the numerical prediction of weather or climate events,the uncertainty of the initial values and/or prediction models can bring the forecast result’s uncertainty.Due to the absence of true states,studies on this problem mainly focus on the three subproblems of predictability,i.e.,the lower bound of the maximum predictable time,the upper bound of the prediction error,and the lower bound of the maximum allowable initial error.Aimed at the problem of the lower bound estimation of the maximum allowable initial error,this study first illustrates the shortcoming of the existing estimation,and then presents a new estimation based on the initial observation precision and proves it theoretically.Furthermore,the new lower bound estimations of both the two-dimensional ikeda model and lorenz96 model are obtained by using the cnop(conditional nonlinear optimal perturbation)method and a pso(particle swarm optimization)algorithm,and the estimated precisions are also analyzed.Besides,the estimations yielded by the existing and new formulas are compared;the results show that the estimations produced by the existing formula are often incorrect.
基金the National Natural Science Foundation of China Grants U1637208 and 71773024.the National Natural Science Foundation of China Grant 11971132.
文摘In this paper,we consider the local discontinuous Galerkin method with generalized alter-nating numerical fluxes for two-dimensional nonlinear Schrödinger equations on Carte-sian meshes.The generalized fluxes not only lead to a smaller magnitude of the errors,but can guarantee an energy conservative property that is useful for long time simulations in resolving waves.By virtue of generalized skew-symmetry property of the discontinuous Galerkin spatial operators,two energy equations are established and stability results con-taining energy conservation of the prime variable as well as auxiliary variables are shown.To derive optimal error estimates for nonlinear Schrödinger equations,an additional energy equation is constructed and two a priori error assumptions are used.This,together with properties of some generalized Gauss-Radau projections and a suitable numerical initial condition,implies optimal order of k+1.Numerical experiments are given to demonstrate the theoretical results.