A PRECONDITIONER FOR DOMAIN DECOMPOSITION METHODS WITH MORTAR LAGRANGE MULTIPLIERS

In this paper we consider the domain decomposition methods with mortar element Lagrange multipliers to two-dimensional elliptic problems.We shall construct a kind of simple preconditioners for the corresponding interface equation.It will be shown that condition number of the preconditioned interface matrix is almost optimal.

A PRECONDITIONER FOR THREE-DIMENSIONAL DOMAIN DECOMPOSITION METHODS WITH LAGRANGE MULTIPLIERS

In this paper we consider domain decomposition methods for three-dimensional elliptic problems with Lagrange multipliers, and construct a kind of simple preconditioner for the corresponding interface equation. It will be shown that condition number of the resulting preconditioned interface matrix is almost optimal.

SUPERCONVERGENCE OF THE NON-CONFORMING DOMAIN DECOMPOSED FINITE ELEMENT METHOD WITH LAGRANGE MULTIPLIERS FOR PARABOLIC PROBLEMS

Abstract. In this paper which is motivated by computation on parallel machine, we showthat the superconvergence results of the finite element method(FEM) with Lagrange mul-tipliers based on domain decomposition method (DDM) with nonmatching grids can becarried over to parabolic problems. The main idea of this paper is to achieve the combina-tion of parallel computational method with the higher accuracy technique by interpolationfinite element postprocessing.

Chebyshev-Lagrange Multipliers Technique for Vibration Analysis of Functionally Graded Material Beams using Various Beam Theories

In this paper,a new technique for analysing functionally graded material(FGM)beams using the Chebyshev polynomials and Lagrange multipliers with various beam theories is presented.By utilizing the inner products and the Chebyshev polynomials’orthogonality properties incorporated with Lagrange multipliers,we can combine the governing equation and boundary conditions to yield the matrix equations with explicit weighting coefficients.Numerical examples are provided for vibration analysis of various beam theories and assumptions.Based on numerical evaluations,it is revealed that the proposed technique can efficiently achieve good agreement with those of the references.

A Modified Lagrange Method for Solving Convex Quadratic Optimization Problems

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.

Modified Augmented Lagrange Multiplier Methods for Large-Scale Chemical Process Optimization

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.

Accelerated Matrix Recovery via Random Projection Based on Inexact Augmented Lagrange Multiplier Method

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.

A Parameter-Free Approach to Determine the Lagrange Multiplier in the Level Set Method by Using the BESO

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.

Modified Lagrange Multiplier Method and Generalized Variational Principle in Fluid Mechanics

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.

Fully Coupled Fluid-Structure Interaction Model Based on Distributed Lagrange Multiplier/Fictitious Domain Method

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.

Study on the tradeoff between interpretability and precision in fuzzy modeling

An approach to identifying fuzzy models considering both interpretability and precision was proposed. Firstly, interpretability issues about fuzzy models were analyzed. Then, a heuristic strategy was used to select input variables by increasing the number of input variables, and the Gustafson-Kessel fuzzy clustering algorithm, combined with the least square method, was used to identify the fuzzy model. Subsequently, an interpretability measure was described by the product of the number of input variables and the number of rules, while precision was weighted by root mean square error, and the selection objective function concerning interpretability and precision was defined. Given the maximum and minimum number of input variables and rules, a set of fuzzy models was constructed. Finally, the optimal fuzzy model was selected by the objective function, and was optimized by a genetic algorithm to achieve a good tradeoff between interpretability and precision. The performance of the proposed method was illustrated by the well-known Box-Jenkins gas furnace benchmark; the results demonstrate its validity.

Multi-phase signal setting and capacity of signalized intersection

A capacity model of multi-phase signalized intersections is derived by a stopping-line method. It is simplified with two normal situations： one situation involves one straight lane and one left-turn lane; the other situation involves two straight lanes and one left-turn lane. The results show that the capacity is mainly relative to signal cycle length, phase length, intersection layout and following time. With regard to the vehicles arrival rates, the optimal model is derived based on each phase＇s remaining time balance, and it is solved by Lagrange multipliers. Therefore, the calculation models of the optimal signal cycle length and phase lengths are derived and simplified. Compared to the existing models, the proposed model is more convenient and practical. Finally, a practical intersection is chosen and its signal cycles and phase lengths are calculated by the proposed model.

Modeling and analysis of rigid multibody systems with driving constraints and frictional translation joints

An approach is proposed for modeling and anal- yses of rigid multibody systems with frictional translation joints and driving constraints. The geometric constraints of translational joints with small clearance are treated as bilat- eral constraints by neglecting the impact between sliders and guides. Firstly, the normal forces acting on sliders, the driv- ing constraint forces （or moments） and the constraint forces of smooth revolute joints are all described by complementary conditions. The frictional contacts are characterized by a set- valued force law of Coulomb＇s dry friction. Combined with the theory of the horizontal linear complementarity problem （HLCP）, an event-driven scheme is used to detect the transi- tions of the contact situation between sliders and guides, and the stick-slip transitions of sliders, respectively. And then, all constraint forces in the system can be computed easily. Secondly, the dynamic equations of multibody systems are written at the acceleration-force level by the Lagrange multiplier technique, and the Baumgarte stabilization method is used to reduce the constraint drift. Finally, a numerical example is given to show some non-smooth dynamical behaviors of the studied system. The obtained results validate the feasibility of algorithm and the effect of constraint stabilization.

Airy, Beltrami, Maxwell, Einstein and Lanczos Potentials Revisited

The purpose of this paper is to revisit the well known potentials, also called stress functions, needed in order to study the parametrizations of the stress equations, respectively provided by G.B. Airy (1863) for 2-dimensional elasticity, then by E. Beltrami (1892), J.C. Maxwell (1870) for 3-dimensional elasticity, finally by A. Einstein (1915) for 4-dimensional elasticity, both with a variational procedure introduced by C. Lanczos (1949, 1962) in order to relate potentials to Lagrange multipliers. Using the methods of Algebraic Analysis, namely mixing differential geometry with homological algebra and combining the double duality test involved with the Spencer cohomology, we shall be able to extend these results to an arbitrary situation with an arbitrary dimension n. We shall also explain why double duality is perfectly adapted to variational calculus with differential constraints as a way to eliminate the corresponding Lagrange multipliers. For example, the canonical parametrization of the stress equations is just described by the formal adjoint of the components of the linearized Riemann tensor considered as a linear second order differential operator but the minimum number of potentials needed is equal to for any minimal parametrization, the Einstein parametrization being “in between” with potentials. We provide all the above results without even using indices for writing down explicit formulas in the way it is done in any textbook today, but it could be strictly impossible to obtain them without using the above methods. We also revisit the possibility (Maxwell equations of electromagnetism) or the impossibility (Einstein equations of gravitation) to obtain canonical or minimal parametrizations for various equations of physics. It is nevertheless important to notice that, when n and the algorithms presented are known, most of the calculations can be achieved by using computers for the corresponding symbolic computations. Finally, though the paper is mathematically oriented as it aims providing new insights towards the mathematical foundations of general relativity, it is written in a rather self-contained way.

A FICTITIOUS DOMAIN METHOD FOR DIRICHLET PROBLEM AND ITS APPLICATIONS TO GENER ALIZED STOKES PROBLEM

This paper discusses a fictitious domain method for the linear Dirichlet problem and its applications to the generalized Stokes problem. This method treats Dirichlet boundary condit ion via a Lagrange multiplier technique and is well suited to the no-slip bound ary condition in viscous flow problems. In order to improve the accuracy of solu tions, meshes are refined according to the a posteriori error estimate. The mini -element discretization is applied to solve the generalized Stokes problem. Fin ally, some numerical results to validate this method are presented for partial d ifferential equations with Dirichlet boundary condition.

A Local Deviation Constraint Based Non-Rigid Structure From Motion Approach

In many traditional non-rigid structure from motion(NRSFM)approaches,the estimation results of part feature points may significantly deviate from their true values because only the overall estimation error is considered in their models.Aimed at solving this issue,a local deviation-constrained-based column-space-fitting approach is proposed in this paper to alleviate estimation deviation.In our work,an effective model is first constructed with two terms:the overall estimation error,which is computed by a linear subspace representation,and a constraint term,which is based on the variance of the reconstruction error for each frame.Furthermore,an augmented Lagrange multipliers(ALM)iterative algorithm is presented to optimize the proposed model.Moreover,a convergence analysis is performed with three steps for the optimization process.As both the overall estimation error and the local deviation are utilized,the proposed method can achieve a good estimation performance and a relatively uniform estimation error distribution for different feature points.Experimental results on several widely used synthetic sequences and real sequences demonstrate the effectiveness and feasibility of the proposed algorithm.

Direct Assessment of Interlaminar Stresses in Composite Multilayered Plates Using a Layer-Wise Mixed Finite Element Model

An accurate determination of intedaminar transversal stresses in composite multilayered plates, especially near free-edge, is of great importance in the study of inter-ply damage modes, mainly in the initiation and growth of delamination. In this paper, interlaminar stresses are determined by layer-wise mixed finite element model. Each layer is analyzed as an isolated one where the displacement continuity is ensured by means of Lagrange multipliers （which represent the statics variables）. This procedure allows the authors to work with any single plate model, obtaining the interlaminar stresses directly without loss of precision. The FSDT （first shear deformation theory） with transverse normal strain effects included is assumed in each layer, but Lagrange polynomials are used to describe the kinematic instead of Taylor＇s polynomial functions of the thickness coordinates, as is common. This expansion allows the authors to pose the interlaminar displacements compatibility simpler than the second one. The in-plane domain of the plate is discretized by four-node quadrilateral elements, both to the field of displacement and to the Lagrange multipliers. The mixed interpolation of tensorial components technique is applied to avoid the shear-locking in the finite element model. Several examples were carried out and the results have been satisfactorily compared with those available in the literature.

Determination of Natural Frequencies of Multilayered Plates by Mixed Finite Elements

Natural frequencies for multilayer plates are calculated by mixed finite element method. The main object of this paper is to use the mixed model for multilayer plates, analyzing each layer as an isolated plate, where the continuity of displacements is achieved by Lagrange multipliers （representing static variables）. This procedure allows us to work with any model for single plate （so as to ensure the proper behavior of each layer）, and the complexity of the multilayer system is avoided by ensuring the condition of displacements by the Lagrange multipliers （static variables）. The plate is discretized by finite element modeling based on a primary hybrid model, where the domain is divided by quadrilateral, both for the displacement field and static variables. This mixed element for plates was implemented and several examples of vibrations have been verified successfully by the results obtained by other methods in the literature.

Lagrangian Relaxation Method for Multiobjective Optimization Methods: Solution Approaches

This paper introduces the Lagrangian relaxation method to solve multiobjective optimization problems. It is often required to use the appropriate technique to determine the Lagrangian multipliers in the relaxation method that leads to finding the optimal solution to the problem. Our analysis aims to find a suitable technique to generate Lagrangian multipliers, and later these multipliers are used in the relaxation method to solve Multiobjective optimization problems. We propose a search-based technique to generate Lagrange multipliers. In our paper, we choose a suitable and well-known scalarization method that transforms the original multiobjective into a scalar objective optimization problem. Later, we solve this scalar objective problem using Lagrangian relaxation techniques. We use Brute force techniques to sort optimum solutions. Finally, we analyze the results, and efficient methods are recommended.

Practical constrained least-square algorithm for moving source location using TDOA and FDOA measurements

By utilizing the time difference of arrival （TDOA） and frequency difference of arrival （FDOA） measurements of signals received at a number of receivers, a constrained least-square （CLS） algorithm for estimating the position and velocity of a moving source is proposed. By utilizing the Lagrange multipliers technique, the known relation between the intermediate variables and the source location coordinates could be exploited to constrain the solution. And without requiring apriori knowledge of TDOA and FDOA measurement noises, the proposed algorithm can satisfy the demand of practical applications. Additionally, on basis of con- volute and polynomial rooting operations, the Lagrange multipliers can be obtained efficiently and robustly allowing real-time imple- mentation and global convergence. Simulation results show that the proposed estimator achieves remarkably better performance than the two-step weighted least square （WLS） approach especially for higher measurement noise level.