In this paper, a modified version of the Classical Lagrange Multiplier method is developed for convex quadratic optimization problems. The method, which is evolved from the first order derivative test for optimality o...In this paper, a modified version of the Classical Lagrange Multiplier method is developed for convex quadratic optimization problems. The method, which is evolved from the first order derivative test for optimality of the Lagrangian function with respect to the primary variables of the problem, decomposes the solution process into two independent ones, in which the primary variables are solved for independently, and then the secondary variables, which are the Lagrange multipliers, are solved for, afterward. This is an innovation that leads to solving independently two simpler systems of equations involving the primary variables only, on one hand, and the secondary ones on the other. Solutions obtained for small sized problems (as preliminary test of the method) demonstrate that the new method is generally effective in producing the required solutions.展开更多
The Lagrange multiplier method plays an important role in establishing generalized variational principles notonly in tluid mechallics. but also in elasticity. Sometimes, however, one may come across variational crisi...The Lagrange multiplier method plays an important role in establishing generalized variational principles notonly in tluid mechallics. but also in elasticity. Sometimes, however, one may come across variational crisis(somemultipliers vanish identically). failing to achieve his aim. The crisis is caused by the fact that the Inultipliers are treatedas independent variables in the process of variatioll. but after identification they become functions of the originalindependent variables. To overcome it, a Inodified Lagrange multiplier method or semi-inverse method has beenproposed to deduce generalized varistional principles. Some e-camples are given to illustrate its convenience andeffectiveness of the novel method.展开更多
The distributed Lagrange multiplier/fictitious domain(DLM/FD)-mixed finite element method is developed and analyzed in this paper for a transient Stokes interface problem with jump coefficients.The semi-and fully disc...The distributed Lagrange multiplier/fictitious domain(DLM/FD)-mixed finite element method is developed and analyzed in this paper for a transient Stokes interface problem with jump coefficients.The semi-and fully discrete DLM/FD-mixed finite element scheme are developed for the first time for this problem with a moving interface,where the arbitrary Lagrangian-Eulerian(ALE)technique is employed to deal with the moving and immersed subdomain.Stability and optimal convergence properties are obtained for both schemes.Numerical experiments are carried out for different scenarios of jump coefficients,and all theoretical results are validated.展开更多
Chemical process optimization can be described as large-scale nonlinear constrained minimization. The modified augmented Lagrange multiplier methods (MALMM) for large-scale nonlinear constrained minimization are studi...Chemical process optimization can be described as large-scale nonlinear constrained minimization. The modified augmented Lagrange multiplier methods (MALMM) for large-scale nonlinear constrained minimization are studied in this paper. The Lagrange function contains the penalty terms on equality and inequality constraints and the methods can be applied to solve a series of bound constrained sub-problems instead of a series of unconstrained sub-problems. The steps of the methods are examined in full detail. Numerical experiments are made for a variety of problems, from small to very large-scale, which show the stability and effectiveness of the methods in large-scale
problems.展开更多
A parameter-free approach is proposed to determine the Lagrange multiplier for the constraint of material volume in the level set method.It is inspired by the procedure of determining the threshold of sensitivity numb...A parameter-free approach is proposed to determine the Lagrange multiplier for the constraint of material volume in the level set method.It is inspired by the procedure of determining the threshold of sensitivity number in the BESO method.It first computes the difference between the volume of current design and the upper bound of volume.Then,the Lagrange multiplier is regarded as the threshold of sensitivity number to remove the redundant material.Numerical examples proved that this approach is effective to constrain the volume.More importantly,there is no parameter in the proposed approach,which makes it convenient to use.In addition,the convergence is stable,and there is no big oscillation.展开更多
In this paper, a unified matrix recovery model was proposed for diverse corrupted matrices. Resulting from the separable structure of the proposed model, the convex optimization problem can be solved efficiently by ad...In this paper, a unified matrix recovery model was proposed for diverse corrupted matrices. Resulting from the separable structure of the proposed model, the convex optimization problem can be solved efficiently by adopting an inexact augmented Lagrange multiplier (IALM) method. Additionally, a random projection accelerated technique (IALM+RP) was adopted to improve the success rate. From the preliminary numerical comparisons, it was indicated that for the standard robust principal component analysis (PCA) problem, IALM+RP was at least two to six times faster than IALM with an insignificant reduction in accuracy; and for the outlier pursuit (OP) problem, IALM+RP was at least 6.9 times faster, even up to 8.3 times faster when the size of matrix was 2 000×2 000.展开更多
This paper, with a finite element method, studies the interaction of a coupled incompressible fluid-rigid structure system with a free surface subjected to external wave excitations. With this fully coupled model, the...This paper, with a finite element method, studies the interaction of a coupled incompressible fluid-rigid structure system with a free surface subjected to external wave excitations. With this fully coupled model, the rigid structure is taken as "fictitious" fluid with zero strain rate. Both fluid and structure are described by velocity and pressure. The whole domain, including fluid region and structure region, is modeled by the incompressible Navier-Stokes equations which are discretized with fixed Eulerian mesh. However, to keep the structure' s rigid body shape and behavior, a rigid body constraint is enforced on the "fictitious" fluid domain by use of the Distributed Lagrange Multipher/Fictitious Domain (DLM/ FD) method which is originally introduced to solve particulate flow problems by Glowinski et al. For the verification of the model presented herein, a 2D numerical wave tank is established to simulate small amplitude wave propagations, and then numerical results are compared with analytical solutions. Finally, a 2D example of fluid-structure interaction under wave dynamic forces provides convincing evidences for the method excellent solution quality and fidelity.展开更多
In two-phase sampling, or double sampling, from a population with size N we take one, relatively large, sample size n. From this relatively large sample we take a small sub-sample size m, which usually costs more per ...In two-phase sampling, or double sampling, from a population with size N we take one, relatively large, sample size n. From this relatively large sample we take a small sub-sample size m, which usually costs more per sample unit than the first one. In double sampling with regression estimators, the sample of the first phase n is used for the estimation of the average of an auxiliary variable X, which should be strongly related to the main variable Y (which is estimated from the sub-sample m). Sampling optimization can be achieved by minimizing cost C with fixed var Y, or by finding a minimum var Y for fixed C. In this paper we optimize sampling with use of Lagrange multipliers, either by minimizing variance of Y and having predetermined cost, or by minimizing cost and having predetermined variance of Y.展开更多
In this paper,we propose a mortar element method with Lagrange multiplier for incompressible Stokes problem,i.e.,the matching constraints of velocity on mortar edges are expressed in terms of Lagrange multipliers.We a...In this paper,we propose a mortar element method with Lagrange multiplier for incompressible Stokes problem,i.e.,the matching constraints of velocity on mortar edges are expressed in terms of Lagrange multipliers.We also present P_1 noncon- forming element attached to the subdomains.By proving inf-sup condition,we derive optimal error estimates for velocity and pressure.Moreover,we obtain satisfactory approximation for normal derivatives of the velocity across the interfaces.展开更多
In this paper, a new augmented Lagrangian function with 4-piecewise linear NCP function is introduced for solving nonlinear programming problems with equality constrained and inequality constrained. It is proved that ...In this paper, a new augmented Lagrangian function with 4-piecewise linear NCP function is introduced for solving nonlinear programming problems with equality constrained and inequality constrained. It is proved that a solution of the original constrained problem and corresponding values of Lagrange multipliers can be found by solving an unconstrained minimization of the augmented Lagrange function. Meanwhile, a new Lagrangian multiplier method corresponding with new augmented Lagrangian function is proposed. And this method is implementable and convergent.展开更多
The supposedly missing dark energy of the cosmos is found quantitatively in a direct analysis without involving ordinary energy. The analysis relies on five dimensional Kaluza-Klein spacetime and a Lagrangian constrai...The supposedly missing dark energy of the cosmos is found quantitatively in a direct analysis without involving ordinary energy. The analysis relies on five dimensional Kaluza-Klein spacetime and a Lagrangian constrained by an auxiliary condition. Employing the Lagrangian multiplier method, it is found that this multiplier is equal to the dark energy of the cosmos and is given by where E is energy, m is mass, c is the speed of light, and λ is the Lagrangian multiplier. The result is in full agreement with cosmic measurements which were awarded the 2011 Nobel Prize in Physics as well as with the interpretation that dark energy is the energy of the quantum wave while ordinary energy is the energy of the quantum particle. Consequently dark energy could not be found directly using our current measurement methods because measurement leads to wave collapse leaving only the quantum particle and its ordinary energy intact.展开更多
This paper is devoted to characterizing the Riemann-Stieltjes operators and pointwise multipliers on F(p, q, s) spaces in the unit ball of C^n which contain many classical function spaces, such as the Bloch space, B...This paper is devoted to characterizing the Riemann-Stieltjes operators and pointwise multipliers on F(p, q, s) spaces in the unit ball of C^n which contain many classical function spaces, such as the Bloch space, BMOA and Q8 spaces. The boundedness and compactness of these operators on F(p, q, s) spaces are characterized by means of an embedding theorem, i.e., F(p,q, s) spaces boundedly embedded into the tent-type spaces Tp,s^∞(μ)展开更多
Recently many research works have been conducted and published regarding fractional order differential equations. There are several approaches available for numerical approximations of the solution of fractional order...Recently many research works have been conducted and published regarding fractional order differential equations. There are several approaches available for numerical approximations of the solution of fractional order diffusion equations. Spectral collocation method based on Lagrange’s basis polynomials to approximate numerical solutions of one-dimensional (1D) space fractional diffusion equations are introduced in this research paper. The proposed form of approximate solution satisfies non-zero Dirichlet’s boundary conditions on both boundaries. Collocation scheme produce a system of first order Ordinary Differential Equations (ODE) from the fractional diffusion equation. We applied this method with four different sets of collocation points to compare their performance.展开更多
The convergence analysis of a nonlinear Lagrange algorithm for solving nonlinear constrained optimization problems with both inequality and equality constraints is explored in detail. The estimates for the derivatives...The convergence analysis of a nonlinear Lagrange algorithm for solving nonlinear constrained optimization problems with both inequality and equality constraints is explored in detail. The estimates for the derivatives of the multiplier mapping and the solution mapping of the proposed algorithm are discussed via the technique of the singular value decomposition of matrix. Based on the estimates, the local convergence results and the rate of convergence of the algorithm are presented when the penalty parameter is less than a threshold under a set of suitable conditions on problem functions. Furthermore, the condition number of the Hessian of the nonlinear Lagrange function with respect to the decision variables is analyzed, which is closely related to efficiency of the algorithm. Finally, the preliminary numericM results for several typical test problems are reported.展开更多
文摘In this paper, a modified version of the Classical Lagrange Multiplier method is developed for convex quadratic optimization problems. The method, which is evolved from the first order derivative test for optimality of the Lagrangian function with respect to the primary variables of the problem, decomposes the solution process into two independent ones, in which the primary variables are solved for independently, and then the secondary variables, which are the Lagrange multipliers, are solved for, afterward. This is an innovation that leads to solving independently two simpler systems of equations involving the primary variables only, on one hand, and the secondary ones on the other. Solutions obtained for small sized problems (as preliminary test of the method) demonstrate that the new method is generally effective in producing the required solutions.
文摘The Lagrange multiplier method plays an important role in establishing generalized variational principles notonly in tluid mechallics. but also in elasticity. Sometimes, however, one may come across variational crisis(somemultipliers vanish identically). failing to achieve his aim. The crisis is caused by the fact that the Inultipliers are treatedas independent variables in the process of variatioll. but after identification they become functions of the originalindependent variables. To overcome it, a Inodified Lagrange multiplier method or semi-inverse method has beenproposed to deduce generalized varistional principles. Some e-camples are given to illustrate its convenience andeffectiveness of the novel method.
基金P.Sun was supported by NSF Grant DMS-1418806C.S.Zhang was partially supported by the National Key Research and Development Program of China(Grant No.2016YFB0201304)+1 种基金the Major Research Plan of National Natural Science Foundation of China(Grant Nos.91430215,91530323)the Key Research Program of Frontier Sciences of CAS.
文摘The distributed Lagrange multiplier/fictitious domain(DLM/FD)-mixed finite element method is developed and analyzed in this paper for a transient Stokes interface problem with jump coefficients.The semi-and fully discrete DLM/FD-mixed finite element scheme are developed for the first time for this problem with a moving interface,where the arbitrary Lagrangian-Eulerian(ALE)technique is employed to deal with the moving and immersed subdomain.Stability and optimal convergence properties are obtained for both schemes.Numerical experiments are carried out for different scenarios of jump coefficients,and all theoretical results are validated.
文摘Chemical process optimization can be described as large-scale nonlinear constrained minimization. The modified augmented Lagrange multiplier methods (MALMM) for large-scale nonlinear constrained minimization are studied in this paper. The Lagrange function contains the penalty terms on equality and inequality constraints and the methods can be applied to solve a series of bound constrained sub-problems instead of a series of unconstrained sub-problems. The steps of the methods are examined in full detail. Numerical experiments are made for a variety of problems, from small to very large-scale, which show the stability and effectiveness of the methods in large-scale
problems.
基金This research work is supported by the National Natural Science Foundation of China(Grant No.51975227).
文摘A parameter-free approach is proposed to determine the Lagrange multiplier for the constraint of material volume in the level set method.It is inspired by the procedure of determining the threshold of sensitivity number in the BESO method.It first computes the difference between the volume of current design and the upper bound of volume.Then,the Lagrange multiplier is regarded as the threshold of sensitivity number to remove the redundant material.Numerical examples proved that this approach is effective to constrain the volume.More importantly,there is no parameter in the proposed approach,which makes it convenient to use.In addition,the convergence is stable,and there is no big oscillation.
基金Supported by National Natural Science Foundation of China (No.51275348)College Students Innovation and Entrepreneurship Training Program of Tianjin University (No.201210056339)
文摘In this paper, a unified matrix recovery model was proposed for diverse corrupted matrices. Resulting from the separable structure of the proposed model, the convex optimization problem can be solved efficiently by adopting an inexact augmented Lagrange multiplier (IALM) method. Additionally, a random projection accelerated technique (IALM+RP) was adopted to improve the success rate. From the preliminary numerical comparisons, it was indicated that for the standard robust principal component analysis (PCA) problem, IALM+RP was at least two to six times faster than IALM with an insignificant reduction in accuracy; and for the outlier pursuit (OP) problem, IALM+RP was at least 6.9 times faster, even up to 8.3 times faster when the size of matrix was 2 000×2 000.
基金This study is supported by the National Natural Science Foundation of China (Grant No50579046) the Science Foundation of Tianjin Municipal Commission of Science and Technology (Grant No043114711)
文摘This paper, with a finite element method, studies the interaction of a coupled incompressible fluid-rigid structure system with a free surface subjected to external wave excitations. With this fully coupled model, the rigid structure is taken as "fictitious" fluid with zero strain rate. Both fluid and structure are described by velocity and pressure. The whole domain, including fluid region and structure region, is modeled by the incompressible Navier-Stokes equations which are discretized with fixed Eulerian mesh. However, to keep the structure' s rigid body shape and behavior, a rigid body constraint is enforced on the "fictitious" fluid domain by use of the Distributed Lagrange Multipher/Fictitious Domain (DLM/ FD) method which is originally introduced to solve particulate flow problems by Glowinski et al. For the verification of the model presented herein, a 2D numerical wave tank is established to simulate small amplitude wave propagations, and then numerical results are compared with analytical solutions. Finally, a 2D example of fluid-structure interaction under wave dynamic forces provides convincing evidences for the method excellent solution quality and fidelity.
文摘In two-phase sampling, or double sampling, from a population with size N we take one, relatively large, sample size n. From this relatively large sample we take a small sub-sample size m, which usually costs more per sample unit than the first one. In double sampling with regression estimators, the sample of the first phase n is used for the estimation of the average of an auxiliary variable X, which should be strongly related to the main variable Y (which is estimated from the sub-sample m). Sampling optimization can be achieved by minimizing cost C with fixed var Y, or by finding a minimum var Y for fixed C. In this paper we optimize sampling with use of Lagrange multipliers, either by minimizing variance of Y and having predetermined cost, or by minimizing cost and having predetermined variance of Y.
基金This work is supported by the National Natural Science Foundation of China under grants (10471067) the Scientific Research Foundation of University under grants (NY207096).
文摘In this paper,we propose a mortar element method with Lagrange multiplier for incompressible Stokes problem,i.e.,the matching constraints of velocity on mortar edges are expressed in terms of Lagrange multipliers.We also present P_1 noncon- forming element attached to the subdomains.By proving inf-sup condition,we derive optimal error estimates for velocity and pressure.Moreover,we obtain satisfactory approximation for normal derivatives of the velocity across the interfaces.
文摘In this paper, a new augmented Lagrangian function with 4-piecewise linear NCP function is introduced for solving nonlinear programming problems with equality constrained and inequality constrained. It is proved that a solution of the original constrained problem and corresponding values of Lagrange multipliers can be found by solving an unconstrained minimization of the augmented Lagrange function. Meanwhile, a new Lagrangian multiplier method corresponding with new augmented Lagrangian function is proposed. And this method is implementable and convergent.
文摘The supposedly missing dark energy of the cosmos is found quantitatively in a direct analysis without involving ordinary energy. The analysis relies on five dimensional Kaluza-Klein spacetime and a Lagrangian constrained by an auxiliary condition. Employing the Lagrangian multiplier method, it is found that this multiplier is equal to the dark energy of the cosmos and is given by where E is energy, m is mass, c is the speed of light, and λ is the Lagrangian multiplier. The result is in full agreement with cosmic measurements which were awarded the 2011 Nobel Prize in Physics as well as with the interpretation that dark energy is the energy of the quantum wave while ordinary energy is the energy of the quantum particle. Consequently dark energy could not be found directly using our current measurement methods because measurement leads to wave collapse leaving only the quantum particle and its ordinary energy intact.
基金Supported in part by the National Natural Science Foundation of China(11271359)the Fundamental Research Funds for the Central Universities(2014-Ia-037and 2015-IVA-069)
文摘This paper is devoted to characterizing the Riemann-Stieltjes operators and pointwise multipliers on F(p, q, s) spaces in the unit ball of C^n which contain many classical function spaces, such as the Bloch space, BMOA and Q8 spaces. The boundedness and compactness of these operators on F(p, q, s) spaces are characterized by means of an embedding theorem, i.e., F(p,q, s) spaces boundedly embedded into the tent-type spaces Tp,s^∞(μ)
文摘Recently many research works have been conducted and published regarding fractional order differential equations. There are several approaches available for numerical approximations of the solution of fractional order diffusion equations. Spectral collocation method based on Lagrange’s basis polynomials to approximate numerical solutions of one-dimensional (1D) space fractional diffusion equations are introduced in this research paper. The proposed form of approximate solution satisfies non-zero Dirichlet’s boundary conditions on both boundaries. Collocation scheme produce a system of first order Ordinary Differential Equations (ODE) from the fractional diffusion equation. We applied this method with four different sets of collocation points to compare their performance.
基金Supported by the National Natural Science Foundation of China(11201357,81271513 and 91324201)the Fundamental Research Funds for the Central Universities under project(2014-Ia-001)
文摘The convergence analysis of a nonlinear Lagrange algorithm for solving nonlinear constrained optimization problems with both inequality and equality constraints is explored in detail. The estimates for the derivatives of the multiplier mapping and the solution mapping of the proposed algorithm are discussed via the technique of the singular value decomposition of matrix. Based on the estimates, the local convergence results and the rate of convergence of the algorithm are presented when the penalty parameter is less than a threshold under a set of suitable conditions on problem functions. Furthermore, the condition number of the Hessian of the nonlinear Lagrange function with respect to the decision variables is analyzed, which is closely related to efficiency of the algorithm. Finally, the preliminary numericM results for several typical test problems are reported.