In this paper, based on the theorem of the high-order velocity energy, integration and variation principle, the high-order Hamilton＇s principle of general holonomic systems is given. Then, three-order Lagrangian equa...In this paper, based on the theorem of the high-order velocity energy, integration and variation principle, the high-order Hamilton＇s principle of general holonomic systems is given. Then, three-order Lagrangian equations and four-order Lagrangian equations are obtained from the high-order Hamilton＇s principle. Finally, the Hamilton＇s principle of high-order Lagrangian function is given.展开更多
An exact augmented Lagrangian function for the nonlinear nonconvex programming problems with inequality constraints was discussed. Under suitable hypotheses, the relationship was established between the local unconstr...An exact augmented Lagrangian function for the nonlinear nonconvex programming problems with inequality constraints was discussed. Under suitable hypotheses, the relationship was established between the local unconstrained minimizers of the augmented Lagrangian function on the space of problem variables and the local minimizers of the original constrained problem. Furthermore, under some assumptions, the relationship was also established between the global solutions of the augmented Lagrangian function on some compact subset of the space of problem variables and the global solutions of the constrained problem. Therefore, f^om the theoretical point of view, a solution of the inequality constrained problem and the corresponding values of the Lagrange multipliers can be found by the well-known method of multipliers which resort to the unconstrained minimization of the augmented Lagrangian function presented.展开更多
A multifractal model is developed to connect the Lagrangian multifractal dimensions with their Eulerian counterparts. We propose that the characteristic time scale of a Lagrangian quantity should be the Lagrangian tim...A multifractal model is developed to connect the Lagrangian multifractal dimensions with their Eulerian counterparts. We propose that the characteristic time scale of a Lagrangian quantity should be the Lagrangian time scale, and it should not be the Eulerian time scale which was widely used in previous studies on Lagrangian statistics. Using the present model, we can obtain the scaling exponents of Lagrangian velocity structure functions from the existing data or models of scaling exponents of Eulerian velocity structure functions. This model is validated by comparing its prediction with the results of experiments, direct numerical simulations, and the previous theoretical models. The comparison shows that the proposed model can better predict the scaling exponents of Lagrangian velocity structure functions, especially for orders larger than 6.展开更多
In this paper, a new augmented Lagrangian penalty function for constrained optimization problems is studied. The dual properties of the augmented Lagrangian objective penalty function for constrained optimization prob...In this paper, a new augmented Lagrangian penalty function for constrained optimization problems is studied. The dual properties of the augmented Lagrangian objective penalty function for constrained optimization problems are proved. Under some conditions, the saddle point of the augmented Lagrangian objective penalty function satisfies the first-order Karush-Kuhn-Tucker (KKT) condition. Especially, when the KKT condition holds for convex programming its saddle point exists. Based on the augmented Lagrangian objective penalty function, an algorithm is developed for finding a global solution to an inequality constrained optimization problem and its global convergence is also proved under some conditions.展开更多
Tensor robust principal component analysis（TRPCA） problem aims to separate a low-rank tensor and a sparse tensor from their sum. This problem has recently attracted considerable research attention due to its wide ra...Tensor robust principal component analysis（TRPCA） problem aims to separate a low-rank tensor and a sparse tensor from their sum. This problem has recently attracted considerable research attention due to its wide range of potential applications in computer vision and pattern recognition. In this paper, we propose a new model to deal with the TRPCA problem by an alternation minimization algorithm along with two adaptive rankadjusting strategies. For the underlying low-rank tensor, we simultaneously perform low-rank matrix factorizations to its all-mode matricizations; while for the underlying sparse tensor,a soft-threshold shrinkage scheme is applied. Our method can be used to deal with the separation between either an exact or an approximate low-rank tensor and a sparse one. We established the subsequence convergence of our algorithm in the sense that any limit point of the iterates satisfies the KKT conditions. When the iteration stops, the output will be modified by applying a high-order SVD approach to achieve an exactly low-rank final result as the accurate rank has been calculated. The numerical experiments demonstrate that our method could achieve better results than the compared methods.展开更多
A Lagrangian lattice Boltzmann method for solving Euler equations is proposed. The key step in formulating this method is the introduction of the displacement distribution function. The equilibrium distribution functi...A Lagrangian lattice Boltzmann method for solving Euler equations is proposed. The key step in formulating this method is the introduction of the displacement distribution function. The equilibrium distribution function consists of macroscopic Lagrangian variables at time steps n and n + 1. It is different from the standard lattice Boltzmann method. In this method the element, instead of each particle, is required to satisfy the basic law. The element is considered as one large particle, which results in simpler version than the corresponding Eulerian one, because the advection term disappears here. Our numerical examples successfully reproduce the classical results.展开更多
Path integral technique is discussed using Hamilton Jacobi method. The Hamilton Jacobi function of non-natural Lagrangian is obtained using separation of variables method. This function makes an important role in path...Path integral technique is discussed using Hamilton Jacobi method. The Hamilton Jacobi function of non-natural Lagrangian is obtained using separation of variables method. This function makes an important role in path integral quantization. The path integral is obtained as integration over the canonical phase space coordinates, which contains the generalized coordinate q and the generalized momentum p. One illustrative example is considered to explain the application of our formalism.展开更多
In this paper,the forecasting equations of a 2nd-order space-time differential remainder are deduced from the Navier-Stokes primitive equations and Eulerian operator by Taylor-series expansion.Here we introduce a cubi...In this paper,the forecasting equations of a 2nd-order space-time differential remainder are deduced from the Navier-Stokes primitive equations and Eulerian operator by Taylor-series expansion.Here we introduce a cubic spline numerical model(Spline Model for short),which is with a quasi-Lagrangian time-split integration scheme of fitting cubic spline/bicubic surface to all physical variable fields in the atmospheric equations on spherical discrete latitude-longitude mesh.A new algorithm of"fitting cubic spline—time step integration—fitting cubic spline—……"is developed to determine their first-and2nd-order derivatives and their upstream points for time discrete integral to the governing equations in Spline Model.And the cubic spline function and its mathematical polarities are also discussed to understand the Spline Model’s mathematical foundation of numerical analysis.It is pointed out that the Spline Model has mathematical laws of"convergence"of the cubic spline functions contracting to the original functions as well as its 1st-order and 2nd-order derivatives.The"optimality"of the 2nd-order derivative of the cubic spline functions is optimal approximation to that of the original functions.In addition,a Hermite bicubic patch is equivalent to operate on a grid for a 2nd-order derivative variable field.Besides,the slopes and curvatures of a central difference are identified respectively,with a smoothing coefficient of 1/3,three-point smoothing of that of a cubic spline.Then the slopes and curvatures of a central difference are calculated from the smoothing coefficient 1/3 and three-point smoothing of that of a cubic spline,respectively.Furthermore,a global simulation case of adiabatic,non-frictional and"incompressible"model atmosphere is shown with the quasi-Lagrangian time integration by using a global Spline Model,whose initial condition comes from the NCEP reanalysis data,along with quasi-uniform latitude-longitude grids and the so-called"shallow atmosphere"Navier-Stokes primitive equations in the spherical coordinates.The Spline Model,which adopted the Navier-Stokes primitive equations and quasi-Lagrangian time-split integration scheme,provides an initial ideal case of global atmospheric circulation.In addition,considering the essentially non-linear atmospheric motions,the Spline Model could judge reasonably well simple points of any smoothed variable field according to its fitting spline curvatures that must conform to its physical interpretation.展开更多
基金the Natural Science Foundation of Jiangxi Provincethe Foundation of Education Department of Jiangxi Province under Grant No.[2007]136
文摘In this paper, based on the theorem of the high-order velocity energy, integration and variation principle, the high-order Hamilton＇s principle of general holonomic systems is given. Then, three-order Lagrangian equations and four-order Lagrangian equations are obtained from the high-order Hamilton＇s principle. Finally, the Hamilton＇s principle of high-order Lagrangian function is given.
文摘An exact augmented Lagrangian function for the nonlinear nonconvex programming problems with inequality constraints was discussed. Under suitable hypotheses, the relationship was established between the local unconstrained minimizers of the augmented Lagrangian function on the space of problem variables and the local minimizers of the original constrained problem. Furthermore, under some assumptions, the relationship was also established between the global solutions of the augmented Lagrangian function on some compact subset of the space of problem variables and the global solutions of the constrained problem. Therefore, f^om the theoretical point of view, a solution of the inequality constrained problem and the corresponding values of the Lagrange multipliers can be found by the well-known method of multipliers which resort to the unconstrained minimization of the augmented Lagrangian function presented.
基金supported by the National Natural Science Foundation of China(11072247,11021262,and 11232011)National Natural Science Associate Foundation of China(NSAF)(U1230126)973 program of China(2013CB834100)
文摘A multifractal model is developed to connect the Lagrangian multifractal dimensions with their Eulerian counterparts. We propose that the characteristic time scale of a Lagrangian quantity should be the Lagrangian time scale, and it should not be the Eulerian time scale which was widely used in previous studies on Lagrangian statistics. Using the present model, we can obtain the scaling exponents of Lagrangian velocity structure functions from the existing data or models of scaling exponents of Eulerian velocity structure functions. This model is validated by comparing its prediction with the results of experiments, direct numerical simulations, and the previous theoretical models. The comparison shows that the proposed model can better predict the scaling exponents of Lagrangian velocity structure functions, especially for orders larger than 6.
文摘In this paper, a new augmented Lagrangian penalty function for constrained optimization problems is studied. The dual properties of the augmented Lagrangian objective penalty function for constrained optimization problems are proved. Under some conditions, the saddle point of the augmented Lagrangian objective penalty function satisfies the first-order Karush-Kuhn-Tucker (KKT) condition. Especially, when the KKT condition holds for convex programming its saddle point exists. Based on the augmented Lagrangian objective penalty function, an algorithm is developed for finding a global solution to an inequality constrained optimization problem and its global convergence is also proved under some conditions.
基金Supported by the National Natural Science Foundation of China(Grant Nos.6157209961320106008+2 种基金91230103)National Science and Technology Major Project(Grant Nos.2013ZX040050212014ZX04001011)
文摘Tensor robust principal component analysis（TRPCA） problem aims to separate a low-rank tensor and a sparse tensor from their sum. This problem has recently attracted considerable research attention due to its wide range of potential applications in computer vision and pattern recognition. In this paper, we propose a new model to deal with the TRPCA problem by an alternation minimization algorithm along with two adaptive rankadjusting strategies. For the underlying low-rank tensor, we simultaneously perform low-rank matrix factorizations to its all-mode matricizations; while for the underlying sparse tensor,a soft-threshold shrinkage scheme is applied. Our method can be used to deal with the separation between either an exact or an approximate low-rank tensor and a sparse one. We established the subsequence convergence of our algorithm in the sense that any limit point of the iterates satisfies the KKT conditions. When the iteration stops, the output will be modified by applying a high-order SVD approach to achieve an exactly low-rank final result as the accurate rank has been calculated. The numerical experiments demonstrate that our method could achieve better results than the compared methods.
基金This work is supported by National Science Foundation of PRC grant under 10 1710 5 5 and Nature Science Foundation of Shandong province granted under (Q98A0 8116 )
文摘A Lagrangian lattice Boltzmann method for solving Euler equations is proposed. The key step in formulating this method is the introduction of the displacement distribution function. The equilibrium distribution function consists of macroscopic Lagrangian variables at time steps n and n + 1. It is different from the standard lattice Boltzmann method. In this method the element, instead of each particle, is required to satisfy the basic law. The element is considered as one large particle, which results in simpler version than the corresponding Eulerian one, because the advection term disappears here. Our numerical examples successfully reproduce the classical results.
文摘Path integral technique is discussed using Hamilton Jacobi method. The Hamilton Jacobi function of non-natural Lagrangian is obtained using separation of variables method. This function makes an important role in path integral quantization. The path integral is obtained as integration over the canonical phase space coordinates, which contains the generalized coordinate q and the generalized momentum p. One illustrative example is considered to explain the application of our formalism.
文摘In this paper,the forecasting equations of a 2nd-order space-time differential remainder are deduced from the Navier-Stokes primitive equations and Eulerian operator by Taylor-series expansion.Here we introduce a cubic spline numerical model(Spline Model for short),which is with a quasi-Lagrangian time-split integration scheme of fitting cubic spline/bicubic surface to all physical variable fields in the atmospheric equations on spherical discrete latitude-longitude mesh.A new algorithm of"fitting cubic spline—time step integration—fitting cubic spline—……"is developed to determine their first-and2nd-order derivatives and their upstream points for time discrete integral to the governing equations in Spline Model.And the cubic spline function and its mathematical polarities are also discussed to understand the Spline Model’s mathematical foundation of numerical analysis.It is pointed out that the Spline Model has mathematical laws of"convergence"of the cubic spline functions contracting to the original functions as well as its 1st-order and 2nd-order derivatives.The"optimality"of the 2nd-order derivative of the cubic spline functions is optimal approximation to that of the original functions.In addition,a Hermite bicubic patch is equivalent to operate on a grid for a 2nd-order derivative variable field.Besides,the slopes and curvatures of a central difference are identified respectively,with a smoothing coefficient of 1/3,three-point smoothing of that of a cubic spline.Then the slopes and curvatures of a central difference are calculated from the smoothing coefficient 1/3 and three-point smoothing of that of a cubic spline,respectively.Furthermore,a global simulation case of adiabatic,non-frictional and"incompressible"model atmosphere is shown with the quasi-Lagrangian time integration by using a global Spline Model,whose initial condition comes from the NCEP reanalysis data,along with quasi-uniform latitude-longitude grids and the so-called"shallow atmosphere"Navier-Stokes primitive equations in the spherical coordinates.The Spline Model,which adopted the Navier-Stokes primitive equations and quasi-Lagrangian time-split integration scheme,provides an initial ideal case of global atmospheric circulation.In addition,considering the essentially non-linear atmospheric motions,the Spline Model could judge reasonably well simple points of any smoothed variable field according to its fitting spline curvatures that must conform to its physical interpretation.
