This paper discusses the two-block large-scale nonconvex optimization problem with general linear constraints.Based on the ideas of splitting and sequential quadratic optimization(SQO),a new feasible descent method fo...This paper discusses the two-block large-scale nonconvex optimization problem with general linear constraints.Based on the ideas of splitting and sequential quadratic optimization(SQO),a new feasible descent method for the discussed problem is proposed.First,we consider the problem of quadratic optimal(QO)approximation associated with the current feasible iteration point,and we split the QO into two small-scale QOs which can be solved in parallel.Second,a feasible descent direction for the problem is obtained and a new SQO-type method is proposed,namely,splitting feasible SQO(SF-SQO)method.Moreover,under suitable conditions,we analyse the global convergence,strong convergence and rate of superlinear convergence of the SF-SQO method.Finally,preliminary numerical experiments regarding the economic dispatch of a power system are carried out,and these show that the SF-SQO method is promising.展开更多
In this paper,we explore bound preserving and high-order accurate local discontinuous Galerkin(LDG)schemes to solve a class of chemotaxis models,including the classical Keller-Segel(KS)model and two other density-depe...In this paper,we explore bound preserving and high-order accurate local discontinuous Galerkin(LDG)schemes to solve a class of chemotaxis models,including the classical Keller-Segel(KS)model and two other density-dependent problems.We use the convex splitting method,the variant energy quadratization method,and the scalar auxiliary variable method coupled with the LDG method to construct first-order temporal accurate schemes based on the gradient flow structure of the models.These semi-implicit schemes are decoupled,energy stable,and can be extended to high accuracy schemes using the semi-implicit spectral deferred correction method.Many bound preserving DG discretizations are only worked on explicit time integration methods and are difficult to get high-order accuracy.To overcome these difficulties,we use the Lagrange multipliers to enforce the implicit or semi-implicit LDG schemes to satisfy the bound constraints at each time step.This bound preserving limiter results in the Karush-Kuhn-Tucker condition,which can be solved by an efficient active set semi-smooth Newton method.Various numerical experiments illustrate the high-order accuracy and the effect of bound preserving.展开更多
The objective of the paper is to develop a new algorithm for numerical solution of dynamic elastic-plastic strain hardening/softening problems. The gradient dependent model is adopted in the numerical model to overcom...The objective of the paper is to develop a new algorithm for numerical solution of dynamic elastic-plastic strain hardening/softening problems. The gradient dependent model is adopted in the numerical model to overcome the result mesh-sensitivity problem in the dynamic strain softening or strain localization analysis. The equations for the dynamic elastic-plastic problems are derived in terms of the parametric variational principle, which is valid for associated, non-associated and strain softening plastic constitutive models in the finite element analysis. The precise integration method, which has been widely used for discretization in time domain of the linear problems, is introduced for the solution of dynamic nonlinear equations. The new algorithm proposed is based on the combination of the parametric quadratic programming method and the precise integration method and has all the advantages in both of the algorithms. Results of numerical examples demonstrate not only the validity, but also the advantages of the algorithm proposed for the numerical solution of nonlinear dynamic problems.展开更多
The Voronoi cell finite element method (VCFEM) is adopted to overcome the limitations of the classic displacement based finite element method in the numerical simulation of heterogeneous materials. The parametric va...The Voronoi cell finite element method (VCFEM) is adopted to overcome the limitations of the classic displacement based finite element method in the numerical simulation of heterogeneous materials. The parametric variational principle and quadratic programming method are developed for elastic-plastic Voronoi finite element analysis of two-dimensional problems. Finite element formulations are derived and a standard quadratic programming model is deduced from the elastic-plastic equations. Influence of microscopic heterogeneities on the overall mechanical response of heterogeneous materials is studied in detail. The overall properties of heterogeneous materials depend mostly on the size, shape and distribution of the material phases of the microstructure. Numerical examples are presented to demonstrate the validity and effectiveness of the method developed.展开更多
Based on the idea of serendipity element,we construct and analyze the first quadratic serendipity finite volume element method for arbitrary convex polygonalmeshes in this article.The explicit construction of quadrati...Based on the idea of serendipity element,we construct and analyze the first quadratic serendipity finite volume element method for arbitrary convex polygonalmeshes in this article.The explicit construction of quadratic serendipity element shape function is introduced from the linear generalized barycentric coordinates,and the quadratic serendipity element function space based on Wachspress coordinate is selected as the trial function space.Moreover,we construct a family of unified dual partitions for arbitrary convex polygonal meshes,which is crucial to finite volume element scheme,and propose a quadratic serendipity polygonal finite volume element method with fewer degrees of freedom.Finally,under certain geometric assumption conditions,the optimal H1 error estimate for the quadratic serendipity polygonal finite volume element scheme is obtained,and verified by numerical experiments.展开更多
Dynamic positioning capability(DPCap) analysis is essential in the selection of thrusters, in their configuration, and during preliminary investigation of the positioning ability of a newly designed vessel dynamic pos...Dynamic positioning capability(DPCap) analysis is essential in the selection of thrusters, in their configuration, and during preliminary investigation of the positioning ability of a newly designed vessel dynamic positioning system.DPCap analysis can help determine the maximum environmental forces, in which the DP system can counteract in given headings. The accuracy of the DPCap analysis is determined by the precise estimation of the environmental forces as well as the effectiveness of the thrust allocation logic. This paper is dedicated to developing an effective and efficient software program for the DPCap analysis for marine vessels. Estimation of the environmental forces can be obtained by model tests, hydrodynamic computation and empirical formulas. A quadratic programming method is adopted to allocate the total thrust on every thruster of the vessel. A detailed description of the thrust allocation logic of the software program is given. The effectiveness of the new program DPCap Polar Plot(DPCPP)was validated by a DPCap analysis for a supply vessel. The present study indicates that the developed program can be used in the DPCap analysis for marine vessels. Moreover, DPCap analysis considering the thruster failure mode might give guidance to the designers of vessels whose thrusters need to be safer.展开更多
A new algorithm is developed based on the parametric variational principle for elastic-plastic analysis of Cosserat continuum. The governing equations of the classic elastic-plastic problem are regularized by adding r...A new algorithm is developed based on the parametric variational principle for elastic-plastic analysis of Cosserat continuum. The governing equations of the classic elastic-plastic problem are regularized by adding rotational degrees of freedom to the conventional translational degrees of freedom in conventional continuum mechanics. The parametric potential energy princi- ple of the Cosserat theory is developed, from which the finite element formulation of the Cosserat theory and the corresponding parametric quadratic programming model are constructed. Strain localization problems are computed and the mesh independent results are obtained.展开更多
A new algorithm for solving the three-dimensional elastic contact problem with friction is presented. The algorithm is a non-interior smoothing algorithm based on an NCP-function. The parametric variational principle ...A new algorithm for solving the three-dimensional elastic contact problem with friction is presented. The algorithm is a non-interior smoothing algorithm based on an NCP-function. The parametric variational principle and parametric quadratic programming method were applied to the analysis of three-dimensional frictional contact problem. The solution of the contact problem was finally reduced to a linear complementarity problem, which was reformulated as a system of nonsmooth equations via an NCP-function. A smoothing approximation to the nonsmooth equations was given by the aggregate function. A Newton method was used to solve the resulting smoothing nonlinear equations. The algorithm presented is easy to understand and implement. The reliability and efficiency of this algorithm are demonstrated both by the numerical experiments of LCP in mathematical way and the examples of contact problems in mechanics.展开更多
The time delay estimation (TDE) of two different received signals from the same source has attracted many interests of researchers in the field of signal processing. A method described for precisely estimating time de...The time delay estimation (TDE) of two different received signals from the same source has attracted many interests of researchers in the field of signal processing. A method described for precisely estimating time delay in this paper is based on the assumption that little priori knowledge on statistical characteristics is available for the received signals. The variance of the estimate is derived . The basic architecture of this method is to use the adaptive noise canceller, in the steady state , and to interpolate the weight coefficients by using a generalized quadratic interpolation matrix. The formula of the time delay estimation is presented . The method proposed by F.A. Reed is a special case of this method . The hardware implementation is much easier than that of the conventional time delay estimation method . The results of the system simulation and the experimental results at sea show a good agreement with the theoretical analysis.展开更多
In this paper,we mainly focus on proving the existence of lump solutions to a generalized(3+1)-dimensional nonlinear differential equation.Hirota’s bilinear method and a quadratic function method are employed to deri...In this paper,we mainly focus on proving the existence of lump solutions to a generalized(3+1)-dimensional nonlinear differential equation.Hirota’s bilinear method and a quadratic function method are employed to derive the lump solutions localized in the whole plane for a(3+1)-dimensional nonlinear differential equation.Three examples of such a nonlinear equation are presented to investigate the exact expressions of the lump solutions.Moreover,the 3d plots and corresponding density plots of the solutions are given to show the space structures of the lump waves.In addition,the breath-wave solutions and several interaction solutions of the(3+1)-dimensional nonlinear differential equation are obtained and their dynamics are analyzed.展开更多
This paper indicates the possible difficulties for applying the interior point method to NPcomplete problems,transforms an NP-complete problem into a nonconvex quadratic program and then develops some convexity theori...This paper indicates the possible difficulties for applying the interior point method to NPcomplete problems,transforms an NP-complete problem into a nonconvex quadratic program and then develops some convexity theories for it. Lastly it proposes an algorithm which uses Karmarkar's algorithm as a subroutine. The finite convergence of this algorithm is also proved.展开更多
Quantitative data are essential to an appro-priate characterization of vegetation.In the past few years,considerable attention has been paid to vegetation sampling techniques.A number of methods have been developed fo...Quantitative data are essential to an appro-priate characterization of vegetation.In the past few years,considerable attention has been paid to vegetation sampling techniques.A number of methods have been developed for plant density estimations that utilize spacing distances instead of fixed-area quadrats.In this paper,we review the main distance methods for estimating density and propose a new distance method denominated the quartered neighbor method.In this method,the sampling point is considered the center,and the area around it is divided into four quadrants.The distance from the closest individual in each quadrant to its closest neighbor in the same quadrant is measured,and the average of them is the distance we need.It is actually an integration of two old distance methods,the nearest neighbor method,and the point-centered quarter method.With our new method and an old distance method(the point-centered quarter method),we calculated the average spacing distances of the Larix principis-rupprechtii population in the larch forests of the Donglingshan Mountain.Comparing the two methods with the quadrat method,we found they were almost the same in accuracy,but the precision of the new one was better.Meanwhile,it is adequate in sampling intensity and adaptable for general use in rapid ecological survey work.展开更多
In this paper,we construct an H1-conforming quadratic finite element on convex polygonal meshes using the generalized barycentric coordinates.The element has optimal approximation rates.Using this quadratic element,tw...In this paper,we construct an H1-conforming quadratic finite element on convex polygonal meshes using the generalized barycentric coordinates.The element has optimal approximation rates.Using this quadratic element,two stable discretizations for the Stokes equations are developed,which can be viewed as the extensions of the P2-P0 and the Q2-(discontinuous)P1 elements,respectively,to polygonal meshes.Numerical results are presented,which support our theoretical claims.展开更多
In the paper,we propose a novel linearly implicit structure-preserving algorithm,which is derived by combing the invariant energy quadratization approach with the exponential time differencing method,to construct effi...In the paper,we propose a novel linearly implicit structure-preserving algorithm,which is derived by combing the invariant energy quadratization approach with the exponential time differencing method,to construct efficient and accurate time discretization scheme for a large class of Hamiltonian partial differential equations(PDEs).The proposed scheme is a linear system,and can be solved more efficient than the original energy-preserving ex-ponential integrator scheme which usually needs nonlinear iterations.Various experiments are performed to verify the conservation,efficiency and good performance at relatively large time step in long time computations.展开更多
基金supported by the National Natural Science Foundation of China(12171106)the Natural Science Foundation of Guangxi Province(2020GXNSFDA238017 and 2018GXNSFFA281007)the Shanghai Sailing Program(21YF1430300)。
文摘This paper discusses the two-block large-scale nonconvex optimization problem with general linear constraints.Based on the ideas of splitting and sequential quadratic optimization(SQO),a new feasible descent method for the discussed problem is proposed.First,we consider the problem of quadratic optimal(QO)approximation associated with the current feasible iteration point,and we split the QO into two small-scale QOs which can be solved in parallel.Second,a feasible descent direction for the problem is obtained and a new SQO-type method is proposed,namely,splitting feasible SQO(SF-SQO)method.Moreover,under suitable conditions,we analyse the global convergence,strong convergence and rate of superlinear convergence of the SF-SQO method.Finally,preliminary numerical experiments regarding the economic dispatch of a power system are carried out,and these show that the SF-SQO method is promising.
文摘In this paper,we explore bound preserving and high-order accurate local discontinuous Galerkin(LDG)schemes to solve a class of chemotaxis models,including the classical Keller-Segel(KS)model and two other density-dependent problems.We use the convex splitting method,the variant energy quadratization method,and the scalar auxiliary variable method coupled with the LDG method to construct first-order temporal accurate schemes based on the gradient flow structure of the models.These semi-implicit schemes are decoupled,energy stable,and can be extended to high accuracy schemes using the semi-implicit spectral deferred correction method.Many bound preserving DG discretizations are only worked on explicit time integration methods and are difficult to get high-order accuracy.To overcome these difficulties,we use the Lagrange multipliers to enforce the implicit or semi-implicit LDG schemes to satisfy the bound constraints at each time step.This bound preserving limiter results in the Karush-Kuhn-Tucker condition,which can be solved by an efficient active set semi-smooth Newton method.Various numerical experiments illustrate the high-order accuracy and the effect of bound preserving.
文摘The objective of the paper is to develop a new algorithm for numerical solution of dynamic elastic-plastic strain hardening/softening problems. The gradient dependent model is adopted in the numerical model to overcome the result mesh-sensitivity problem in the dynamic strain softening or strain localization analysis. The equations for the dynamic elastic-plastic problems are derived in terms of the parametric variational principle, which is valid for associated, non-associated and strain softening plastic constitutive models in the finite element analysis. The precise integration method, which has been widely used for discretization in time domain of the linear problems, is introduced for the solution of dynamic nonlinear equations. The new algorithm proposed is based on the combination of the parametric quadratic programming method and the precise integration method and has all the advantages in both of the algorithms. Results of numerical examples demonstrate not only the validity, but also the advantages of the algorithm proposed for the numerical solution of nonlinear dynamic problems.
基金Project supported by the National Natural Science Foundation of China(Nos.10225212, 10421002 and 10332010)the NCET Program provided by the Ministry of Education and the National Key Basic Research Special Foundation of China (No.2005CB321704)
文摘The Voronoi cell finite element method (VCFEM) is adopted to overcome the limitations of the classic displacement based finite element method in the numerical simulation of heterogeneous materials. The parametric variational principle and quadratic programming method are developed for elastic-plastic Voronoi finite element analysis of two-dimensional problems. Finite element formulations are derived and a standard quadratic programming model is deduced from the elastic-plastic equations. Influence of microscopic heterogeneities on the overall mechanical response of heterogeneous materials is studied in detail. The overall properties of heterogeneous materials depend mostly on the size, shape and distribution of the material phases of the microstructure. Numerical examples are presented to demonstrate the validity and effectiveness of the method developed.
基金supported by the National Natural Science Foundation of China(Nos.11871009,12271055)the Foundation of LCP and the Foundation of CAEP(CX20210044).
文摘Based on the idea of serendipity element,we construct and analyze the first quadratic serendipity finite volume element method for arbitrary convex polygonalmeshes in this article.The explicit construction of quadratic serendipity element shape function is introduced from the linear generalized barycentric coordinates,and the quadratic serendipity element function space based on Wachspress coordinate is selected as the trial function space.Moreover,we construct a family of unified dual partitions for arbitrary convex polygonal meshes,which is crucial to finite volume element scheme,and propose a quadratic serendipity polygonal finite volume element method with fewer degrees of freedom.Finally,under certain geometric assumption conditions,the optimal H1 error estimate for the quadratic serendipity polygonal finite volume element scheme is obtained,and verified by numerical experiments.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.51179103 and51709170)the 7th Generation Ultra Deep Water Drilling Unit Innovation Project and the Shanghai Sailing Program(Grant No.17YF1409700)
文摘Dynamic positioning capability(DPCap) analysis is essential in the selection of thrusters, in their configuration, and during preliminary investigation of the positioning ability of a newly designed vessel dynamic positioning system.DPCap analysis can help determine the maximum environmental forces, in which the DP system can counteract in given headings. The accuracy of the DPCap analysis is determined by the precise estimation of the environmental forces as well as the effectiveness of the thrust allocation logic. This paper is dedicated to developing an effective and efficient software program for the DPCap analysis for marine vessels. Estimation of the environmental forces can be obtained by model tests, hydrodynamic computation and empirical formulas. A quadratic programming method is adopted to allocate the total thrust on every thruster of the vessel. A detailed description of the thrust allocation logic of the software program is given. The effectiveness of the new program DPCap Polar Plot(DPCPP)was validated by a DPCap analysis for a supply vessel. The present study indicates that the developed program can be used in the DPCap analysis for marine vessels. Moreover, DPCap analysis considering the thruster failure mode might give guidance to the designers of vessels whose thrusters need to be safer.
基金Project supported by the National Natural Sciences Foundation (Nos. 50679013,10421202 and 10225212)the Program for Changjiang Scholars and Innovative Research Team in Universities of China (PCSIRT)the National Key Basic Research Special Foundation of China (No. 2005CB321704)
文摘A new algorithm is developed based on the parametric variational principle for elastic-plastic analysis of Cosserat continuum. The governing equations of the classic elastic-plastic problem are regularized by adding rotational degrees of freedom to the conventional translational degrees of freedom in conventional continuum mechanics. The parametric potential energy princi- ple of the Cosserat theory is developed, from which the finite element formulation of the Cosserat theory and the corresponding parametric quadratic programming model are constructed. Strain localization problems are computed and the mesh independent results are obtained.
文摘A new algorithm for solving the three-dimensional elastic contact problem with friction is presented. The algorithm is a non-interior smoothing algorithm based on an NCP-function. The parametric variational principle and parametric quadratic programming method were applied to the analysis of three-dimensional frictional contact problem. The solution of the contact problem was finally reduced to a linear complementarity problem, which was reformulated as a system of nonsmooth equations via an NCP-function. A smoothing approximation to the nonsmooth equations was given by the aggregate function. A Newton method was used to solve the resulting smoothing nonlinear equations. The algorithm presented is easy to understand and implement. The reliability and efficiency of this algorithm are demonstrated both by the numerical experiments of LCP in mathematical way and the examples of contact problems in mechanics.
文摘The time delay estimation (TDE) of two different received signals from the same source has attracted many interests of researchers in the field of signal processing. A method described for precisely estimating time delay in this paper is based on the assumption that little priori knowledge on statistical characteristics is available for the received signals. The variance of the estimate is derived . The basic architecture of this method is to use the adaptive noise canceller, in the steady state , and to interpolate the weight coefficients by using a generalized quadratic interpolation matrix. The formula of the time delay estimation is presented . The method proposed by F.A. Reed is a special case of this method . The hardware implementation is much easier than that of the conventional time delay estimation method . The results of the system simulation and the experimental results at sea show a good agreement with the theoretical analysis.
基金supported by the National Natural Science Foundation of China(Nos.12101572,12371256)2023 Shanxi Province Graduate Innovation Project(No.2023KY614)the 19th Graduate Science and Technology Project of North University of China(No.20231943)。
文摘In this paper,we mainly focus on proving the existence of lump solutions to a generalized(3+1)-dimensional nonlinear differential equation.Hirota’s bilinear method and a quadratic function method are employed to derive the lump solutions localized in the whole plane for a(3+1)-dimensional nonlinear differential equation.Three examples of such a nonlinear equation are presented to investigate the exact expressions of the lump solutions.Moreover,the 3d plots and corresponding density plots of the solutions are given to show the space structures of the lump waves.In addition,the breath-wave solutions and several interaction solutions of the(3+1)-dimensional nonlinear differential equation are obtained and their dynamics are analyzed.
文摘This paper indicates the possible difficulties for applying the interior point method to NPcomplete problems,transforms an NP-complete problem into a nonconvex quadratic program and then develops some convexity theories for it. Lastly it proposes an algorithm which uses Karmarkar's algorithm as a subroutine. The finite convergence of this algorithm is also proved.
基金The study was supported by the National Natural Science Foundation of China(Grant No.30870399)the Teachers’Foundation of Education Ministry of China.
文摘Quantitative data are essential to an appro-priate characterization of vegetation.In the past few years,considerable attention has been paid to vegetation sampling techniques.A number of methods have been developed for plant density estimations that utilize spacing distances instead of fixed-area quadrats.In this paper,we review the main distance methods for estimating density and propose a new distance method denominated the quartered neighbor method.In this method,the sampling point is considered the center,and the area around it is divided into four quadrants.The distance from the closest individual in each quadrant to its closest neighbor in the same quadrant is measured,and the average of them is the distance we need.It is actually an integration of two old distance methods,the nearest neighbor method,and the point-centered quarter method.With our new method and an old distance method(the point-centered quarter method),we calculated the average spacing distances of the Larix principis-rupprechtii population in the larch forests of the Donglingshan Mountain.Comparing the two methods with the quadrat method,we found they were almost the same in accuracy,but the precision of the new one was better.Meanwhile,it is adequate in sampling intensity and adaptable for general use in rapid ecological survey work.
基金supported by the NSFC grant 11671210 and 12171244.
文摘In this paper,we construct an H1-conforming quadratic finite element on convex polygonal meshes using the generalized barycentric coordinates.The element has optimal approximation rates.Using this quadratic element,two stable discretizations for the Stokes equations are developed,which can be viewed as the extensions of the P2-P0 and the Q2-(discontinuous)P1 elements,respectively,to polygonal meshes.Numerical results are presented,which support our theoretical claims.
基金supported by the National Natural Science Foundation of China(Grant Nos.12171245,11971416,11971242,12301508)by the Natural Science Foundation of Henan Province(Grant No.222300420280)+1 种基金by the Natural Science Foundation of Hunan Province(Grant No.2023JJ40656)by the Scientific Research Fund of Xuchang University(Grant No.2024ZD010).
文摘In the paper,we propose a novel linearly implicit structure-preserving algorithm,which is derived by combing the invariant energy quadratization approach with the exponential time differencing method,to construct efficient and accurate time discretization scheme for a large class of Hamiltonian partial differential equations(PDEs).The proposed scheme is a linear system,and can be solved more efficient than the original energy-preserving ex-ponential integrator scheme which usually needs nonlinear iterations.Various experiments are performed to verify the conservation,efficiency and good performance at relatively large time step in long time computations.