A Lagrangian Identification of the Main Sources of Moisture Affecting Northeastern Brazil during Its Pre-Rainy and Rainy Seasons

This work examines the sources of moisture affecting the semi-arid Brazilian Northeast (NEB) during its pre-rainy and rainy season (JFMAM) through a Lagrangian diagnosis method. The FLEXPART model identifies the humidity contributions to the moisture budget over a region through the continuous computation of changes in the specific humidity along back or forward trajectories up to 10 days period. The numerical experiments were done for the period that spans between 2000 and 2004 and results were aggregated on a monthly basis. Results show that besides a minor local recycling component, the vast majority of moisture reaching NEB area is originated in the south Atlantic basin and that the nearby wet Amazon basin bears almost no impact. Moreover, although the maximum precipitation in the ""Poligono das Secas'' region (PS) occurs in March and the maximum precipitation associated with air parcels emanating from the South Atlantic towards PS is observed along January to March, the highest moisture contribution from this oceanic region occurs slightly later (April). A dynamical analysis suggests that the maximum precipitation observed in the PS sector does not coincide with the maximum moisture supply probably due to the combined effect of the Walker and Hadley cells in inhibiting the rising motions over the region in the months following April.

A Lagrangian identification of major sources of moisture over Central Brazil and La Plata Basin

[1] This work examines the main sources of moisture over Central Brazil and La Plata Basin during the year through a new Lagrangian diagnosis method which identifies the humidity contributions to the moisture budget over a region. This methodology computes budgets of evaporation minus precipitation by calculating changes in the specific humidity along back-trajectories for the previous 10 d. The origin of all air masses residing over each region was tracked during a period of 5 years (2000-2004). These regions were selected because they coincide with two centers of action of a known dipole precipitation variability mode observed in different temporal scales (from intra seasonal up to inter decadal timescales) and are related to the climatic variability of the South American Monsoon System. The results suggested the importance of the tropical south Atlantic as a moisture source for Central Brazil, and of recycling for La Plata basin. It seems that the Tropical South Atlantic plays an important role as a moisture source for Central Brazil and La Plata basin along the year, particularly during the austral summer. The north Atlantic is also an additional source for both regions during the austral summer.

Reduced gradient method combined with augmented Lagrangian and barrier for the optimal power flow problem

A new approach for solving the optimal power flow (OPF) problem is established by combining the reduced gradient method and the augmented Lagrangian method with barriers and exploring specific characteristics of the relations between the variables of the OPF problem. Computer simulations on IEEE 14-bus and IEEE 30-bus test systems illustrate the method. (c) 2007 Elsevier Inc. All rights reserved.

A Lagrangian relaxation approach to a coupled lot-sizing and cutting stock problem

Industrial production processes involving both lot-sizing and cutting stock problems are common in many industrial settings. However, they are usually treated in a separate way, which could lead to costly production plans. In this paper, a coupled mathematical model is formulated and a heuristic method based on Lagrangian relaxation is proposed. Computational results prove its effectiveness. (C) 2009 Elsevier B.V. All rights reserved.

Global minimization using an Augmented Lagrangian method with variable lower-level constraints

A novel global optimization method based on an Augmented Lagrangian framework is introduced for continuous constrained nonlinear optimization problems. At each outer iteration k the method requires the epsilon(k)-global minimization of the Augmented Lagrangian with simple constraints, where epsilon(k) -> epsilon. Global convergence to an epsilon-global minimizer of the original problem is proved. The subproblems are solved using the alpha BB method. Numerical experiments are presented.

Structured minimal-memory inexact quasi-Newton method and secant preconditioners for augmented Lagrangian optimization

Augmented Lagrangian methods for large-scale optimization usually require efficient algorithms for minimization with box constraints. On the other hand, active-set box-constraint methods employ unconstrained optimization algorithms for minimization inside the faces of the box. Several approaches may be employed for computing internal search directions in the large-scale case. In this paper a minimal-memory quasi-Newton approach with secant preconditioners is proposed, taking into account the structure of Augmented Lagrangians that come from the popular Powell-Hestenes-Rockafellar scheme. A combined algorithm, that uses the quasi-Newton formula or a truncated-Newton procedure, depending on the presence of active constraints in the penalty-Lagrangian function, is also suggested. Numerical experiments using the Cute collection are presented.

Improving ultimate convergence of an augmented Lagrangian method

Optimization methods that employ the classical Powell-Hestenes-Rockafellar augmented Lagrangian are useful tools for solving nonlinear programming problems. Their reputation decreased in the last 10 years due to the comparative success of interior-point Newtonian algorithms, which are asymptotically faster. In this research, a combination of both approaches is evaluated. The idea is to produce a competitive method, being more robust and efficient than its `pure` counterparts for critical problems. Moreover, an additional hybrid algorithm is defined, in which the interior-point method is replaced by the Newtonian resolution of a Karush-Kuhn-Tucker (KKT) system identified by the augmented Lagrangian algorithm. The software used in this work is freely available through the Tango Project web page:http://www.ime.usp.br/similar to egbirgin/tango/.

Augmented Lagrangian methods under the constant positive linear dependence constraint qualification

Two Augmented Lagrangian algorithms for solving KKT systems are introduced. The algorithms differ in the way in which penalty parameters are updated. Possibly infeasible accumulation points are characterized. It is proved that feasible limit points that satisfy the Constant Positive Linear Dependence constraint qualification are KKT solutions. Boundedness of the penalty parameters is proved under suitable assumptions. Numerical experiments are presented.

Instability of equilibrium points of some Lagrangian systems

In this work we show that, if L is a natural Lagrangian system such that the k-jet of the potential energy ensures it does not have a minimum at the equilibrium and such that its Hessian has rank at least n - 2, then there is an asymptotic trajectory to the associated equilibrium point and so the equilibrium is unstable. This applies, in particular, to analytic potentials with a saddle point and a Hessian with at most 2 null eigenvalues. The result is proven for Lagrangians in a specific form, and we show that the class of Lagrangians we are interested can be taken into this specific form by a subtle change of spatial coordinates. We also consider the extension of this results to systems subjected to gyroscopic forces. (C) 2008 Elsevier Inc. All rights reserved.

Minimal Lagrangian surfaces in the tangent bundle of a Riemannian surface

Given an oriented Riemannian surface (Sigma, g), its tangent bundle T Sigma enjoys a natural pseudo-Kahler structure, that is the combination of a complex structure 2, a pseudo-metric G with neutral signature and a symplectic structure Omega. We give a local classification of those surfaces of T Sigma which are both Lagrangian with respect to Omega and minimal with respect to G. We first show that if g is non-flat, the only such surfaces are affine normal bundles over geodesics. In the flat case there is, in contrast, a large set of Lagrangian minimal surfaces, which is described explicitly. As an application, we show that motions of surfaces in R(3) or R(1)(3) induce Hamiltonian motions of their normal congruences, which are Lagrangian surfaces in TS(2) or TH(2) respectively. We relate the area of the congruence to a second-order functional F = f root H(2) - K dA on the original surface. (C) 2010 Elsevier B.V. All rights reserved.

Construction of Hamiltonian-Minimal Lagrangian submanifolds in Complex Euclidean Space

We describe several families of Lagrangian submanifolds in complex Euclidean space which are H-minimal, i.e. critical points of the volume functional restricted to Hamiltonian variations. We make use of various constructions involving planar, spherical and hyperbolic curves, as well as Legendrian submanifolds of the odd-dimensional unit sphere.

An algebraic theory for generalized Jordan chains and partial signatures in the Lagrangian Grassmannian

We discuss an algebraic theory for generalized Jordan chains and partial signatures, that are invariants associated to sequences of symmetric bilinear forms on a vector space. We introduce an intrinsic notion of partial signatures in the Lagrangian Grassmannian of a symplectic space that does not use local coordinates, and we give a formula for the Maslov index of arbitrary real analytic paths in terms of partial signatures.

Unconstrained Finite Element for Geometrical Nonlinear Dynamics of Shells

This paper presents a positional FEM formulation to deal with geometrical nonlinear dynamics of shells. The main objective is to develop a new FEM methodology based on the minimum potential energy theorem written regarding nodal positions and generalized unconstrained vectors not displacements and rotations. These characteristics are the novelty of the present work and avoid the use of large rotation approximations. A nondimensional auxiliary coordinate system is created, and the change of configuration function is written following two independent mappings from which the strain energy function is derived. This methodology is called positional and, as far as the authors' knowledge goes, is a new procedure to approximated geometrical nonlinear structures. In this paper a proof for the linear and angular momentum conservation property of the Newmark beta algorithm is provided for total Lagrangian description. The proposed shell element is locking free for elastic stress-strain relations due to the presence of linear strain variation along the shell thickness. The curved, high-order element together with an implicit procedure to solve nonlinear equations guarantees precision in calculations. The momentum conserving, the locking free behavior, and the frame invariance of the adopted mapping are numerically confirmed by examples. Copyright (C) 2009 H. B. Coda and R. R. Paccola.

Nonquadratic gauge fixing and global gauge invariance in the effective action

The possibility of having a gauge fixing term in the effective Lagrangian that is not a quadratic expression has been explored in spin-two theories so as to have a propagator that is both traceless and transverse. We first show how this same approach can be used in spontaneously broken gauge theories as an alternate to the 't Hooft gauge fixing which avoids terms quadratic in the scalar fields. This ""nonquadratic"" gauge fixing in the effective action results in two complex fermionic and one real bosonic ghost field. A global gauge invariance involving a fermionic gauge parameter, analogous to the usual Becchi-Rouet-Stora-Tyutin invariance, is present in this effective action.

Position-dependent noncommutativity in quantum mechanics

The model of the position-dependent noncommutativity in quantum mechanics is proposed. We start with given commutation relations between the operators of coordinates [(x) over cap (i), (x) over cap (j)] = omega(ij) ((x) over cap), and construct the complete algebra of commutation relations, including the operators of momenta. The constructed algebra is a deformation of a standard Heisenberg algebra and obeys the Jacobi identity. The key point of our construction is a proposed first-order Lagrangian, which after quantization reproduces the desired commutation relations. Also we study the possibility to localize the noncommutativity.