Numerical study about natural escape and capture routes by the Moon via Lagrangian points L1 and L2

The lunar sphere of influence, whose radius is some 66,300 km, has regions of stable orbits around the Moon and also regions that contain trajectories which, after spending some time around the Moon, escape and are later recaptured by lunar gravity. Both the escape and the capture occur along the Lagrangian equilibrium points L1 and L2. In this study, we mapped out the region of lunar influence considering the restricted three-body Earth-Moon-particle problem and the four-body Sun-Earth-Moon-particle (probe) problem. We identified the stable trajectories, and the escape and capture trajectories through the L I and L2 in plots of the eccentricity versus the semi-major axis as a function of the time that the energy of the osculating lunar trajectory in the two-body Moon-particle problem remains negative. We also investigated the properties of these routes, giving special attention to the fact that they supply a natural mechanism for performing low-energy transfers between the Earth and the Moon, and can thus be useful on a great number of future missions. (C) 2007 Published by Elsevier Ltd on behalf of COSPAR.

Partial hyperbolicity for symplectic diffeomorphisms

We prove that every robustly transitive and every stably ergodic symplectic diffeomorphism on a compact manifold admits a dominated splitting. In fact, these diffeomorphisms are partially hyperbolic. (c) 2005 Elsevier SAS. All rights reserved.

Symplectic actions on coadjoint orbits

We present a compact expression for the field theoretical actions based on the symplectic analysis of coadjoint orbits of Lie groups. The final formula for the action density α c becomes a bilinear form 〈(S, 1/λ), (y, m y)〉, where S is a 1-cocycle of the Lie group (a schwarzian type of derivative in conformai case), λ is a coefficient of the central element of the algebra and script Y sign ≡ (y, m y) is the generalized Maurer-Cartan form. In this way the action is fully determined in terms of the basic group theoretical objects. This result is illustrated on a number of examples, including the superconformal model with N = 2. In this case the method is applied to derive the N = 2 superspace generalization of the D=2 Polyakov (super-) gravity action in a manifest (2, 0) supersymmetric form. As a byproduct we also find a natural (2, 0) superspace generalization of the Beltrami equations for the (2, 0) supersymmetric world-sheet metric describing the transition from the conformal to the chiral gauge.

Deconvolução de Euler: passado, presente e futuro – um tutorial

Neste tutorial apresentamos uma revisão da deconvolução de Euler que consiste de três partes. Na primeira parte, recordamos o papel da clássica formulação da deconvolução de Euler 2D e 3D como um método para localizar automaticamente fontes de campos potenciais anômalas e apontamos as dificuldades desta formulação: a presença de uma indesejável nuvem de soluções, o critério empírico usado para determinar o índice estrutural (um parâmetro relacionado com a natureza da fonte anômala), a exeqüibilidade da aplicação da deconvolução de Euler a levantamentos magnéticos terrestres, e a determinação do mergulho e do contraste de susceptibilidade magnética de contatos geológicos (ou o produto do contraste de susceptibilidade e a espessura quando aplicado a dique fino). Na segunda parte, apresentamos as recentes melhorias objetivando minimizar algumas dificuldades apresentadas na primeira parte deste tutorial. Entre estas melhorias incluem-se: i) a seleção das soluções essencialmente associadas com observações apresentando alta razão sinal-ruído; ii) o uso da correlação entre a estimativa do nível de base da anomalia e a própria anomalia observada ou a combinação da deconvolução de Euler com o sinal analítico para determinação do índice estrutural; iii) a combinação dos resultados de (i) e (ii), permitindo estimar o índice estrutural independentemente do número de soluções; desta forma, um menor número de observações (tal como em levantamentos terrestres) pode ser usado; iv) a introdução de equações adicionais independentes da equação de Euler que permitem estimar o mergulho e o contraste de susceptibilidade das fontes magnéticas 2D. Na terceira parte apresentaremos um prognóstico sobre futuros desenvolvimentos a curto e médio prazo envolvendo a deconvolução de Euler. As principais perspectivas são: i) novos ataques aos problemas selecionados na segunda parte deste tutorial; ii) desenvolvimento de métodos que permitam considerar interferências de fontes localizadas ao lado ou acima da fonte principal, e iii) uso das estimativas de localização da fonte anômala produzidas pela deconvolução de Euler como vínculos em métodos de inversão para obter a delineação das fontes em um ambiente computacional amigável.

Functions and vector fields on C(CPn)-singular manifolds

Let M2n+1 be a C(CPn) -singular manifold. We study functions and vector fields with isolated singularities on M2n+1. A C(CPn) -singular manifold is obtained from a smooth manifold M2n+1 with boundary in the form of a disjoint union of complex projective spaces CPn boolean OR CPn boolean OR ... boolean OR CPn with subsequent capture of a cone over each component of the boundary. Let M2n+1 be a compact C(CPn) -singular manifold with k singular points. The Euler characteristic of M2n+1 is equal to chi(M2n+1) = k(1 - n)/2. Let M2n+1 be a C(CPn)-singular manifold with singular points m(1), ..., m(k). Suppose that, on M2n+1, there exists an almost smooth vector field V (x) with finite number of zeros m(1), ..., m(k), x(1), ..., x(1). Then chi(M2n+1) = Sigma(l)(i=1) ind(x(i)) + Sigma(k)(i=1) ind(m(i)).

Optimal reactive power flow via the modified barrier Lagrangian function approach

A new approach called the Modified Barrier Lagrangian Function (MBLF) to solve the Optimal Reactive Power Flow problem is presented. In this approach, the inequality constraints are treated by the Modified Barrier Function (MBF) method, which has a finite convergence property: i.e. the optimal solution in the MBF method can actually be in the bound of the feasible set. Hence, the inequality constraints can be precisely equal to zero. Another property of the MBF method is that the barrier parameter does not need to be driven to zero to attain the solution. Therefore, the conditioning of the involved Hessian matrix is greatly enhanced. In order to show this, a comparative analysis of the numeric conditioning of the Hessian matrix of the MBLF approach, by the decomposition in singular values, is carried out. The feasibility of the proposed approach is also demonstrated with comparative tests to Interior Point Method (IPM) using various IEEE test systems and two networks derived from Brazilian generation/transmission system. The results show that the MBLF method is computationally more attractive than the IPM in terms of speed, number of iterations and numerical conditioning. (C) 2011 Elsevier B.V. All rights reserved.

Nontwist symplectic maps in tokamaks

We review symplectic nontwist maps that we have introduced to describe Lagrangian transport properties in magnetically confined plasmas in tokamaks. These nontwist maps are suitable to describe the formation and destruction of transport barriers in the shearless region (i.e., near the curve where the twist condition does not hold). The maps can be used to investigate two kinds of problems in plasmas with non-monotonic field profiles: the first is the chaotic magnetic field line transport in plasmas with external resonant perturbations. The second problem is the chaotic particle drift motion caused by electrostatic drift waves. The presented analytical maps, derived from plasma models with equilibrium field profiles and control parameters that are commonly measured in plasma discharges, can be used to investigate long-term transport properties. (C) 2011 Elsevier B.V. All rights reserved.

Lê–Greuel type formula for the Euler obstruction and applications

The Euler obstruction of a function f can be viewed as a generalization of the Milnor number for functions defined on singular spaces. In this work, using the Euler obstruction of a function, we establish several Lê–Greuel type formulas for germs f:(X,0)→(C,0) and g:(X,0)→(C,0). We give applications when g is a generic linear form and when f and g have isolated singularities.

Moduli consensuali nei rapporti tra Stato e autonomie territoriali e ruolo delle Conferenze

Ad oltre un decennio dalla riforma costituzionale del 2001, la non compiuta attuazione della stessa sembra condurre ad una nuova ed assai prossima revisione costituzionale. Nel mutato quadro costituzionale di quest’ultimo periodo l’attenzione e l’interesse sono stati richiamati da una grande quantità di accordi, di intese e di altri moduli consensuali, introdotti ed esitati nei rapporti tra Stato ed autonomie territoriali. Dall’esame del sistema delle Conferenze si è evidenziata la loro indispensabilità ai fini del coordinamento delle azioni politico –amministrative delle autonomie territoriali, ma è anche venuta fuori l’esigenza di ricorrere frequentemente ad altri tipi di moduli consensuali, spesso non tipizzati in ambito legislativo. I principi di sussidiarietà e di leale collaborazione inducono non di rado ad assumere, nella concretezza, forme eterogenee poco chiare e confuse (accordo, intesa, concerto, parere) che, pur considerati i contributi offerti dalla dottrina e dalla giurisprudenza, richiederebbero un intervento da parte del legislatore. L’approfondimento tematico ha rilevato, nel contesto di rapporti tra Stato ed autonomie territoriali, un movimento per così dire ondulatorio tra accese spinte autonomistiche e rimarcate esigenze centralistiche, lasciando ancora nella prospettiva una operativa ed efficiente armonia istituzionale.

Long-Range Pollution Transport: Trans-Atlantic Mechanisms and Lagrangian Modeling Methods

Over the past several decades, it has become apparent that anthropogenic activities have resulted in the large-scale enhancement of the levels of many trace gases throughout the troposphere. More recently, attention has been given to the transport pathway taken by these emissions as they are dispersed throughout the atmosphere. The transport pathway determines the physical characteristics of emissions plumes and therefore plays an important role in the chemical transformations that can occur downwind of source regions. For example, the production of ozone (O3) is strongly dependent upon the transport its precursors undergo. O3 can initially be formed within air masses while still over polluted source regions. These polluted air masses can experience continued O3 production or O3 destruction downwind, depending on the air mass's chemical and transport characteristics. At present, however, there are a number of uncertainties in the relationships between transport and O3 production in the North Atlantic lower free troposphere. The first phase of the study presented here used measurements made at the Pico Mountain observatory and model simulations to determine transport pathways for US emissions to the observatory. The Pico Mountain observatory was established in the summer of 2001 in order to address the need to understand the relationships between transport and O3 production. Measurements from the observatory were analyzed in conjunction with model simulations from the Lagrangian particle dispersion model (LPDM), FLEX-PART, in order to determine the transport pathway for events observed at the Pico Mountain observatory during July 2003. A total of 16 events were observed, 4 of which were analyzed in detail. The transport time for these 16 events varied from 4.5 to 7 days, while the transport altitudes over the ocean ranged from 2-8 km, but were typically less than 3 km. In three of the case studies, eastward advection and transport in a weak warm conveyor belt (WCB) airflow was responsible for the export of North American emissions into the FT, while transport in the FT was governed by easterly winds driven by the Azores/Bermuda High (ABH) and transient northerly lows. In the fourth case study, North American emissions were lofted to 6-8 km in a WCB before being entrained in the same cyclone's dry airstream and transported down to the observatory. The results of this study show that the lower marine FT may provide an important transport environment where O3 production may continue, in contrast to transport in the marine boundary layer, where O3 destruction is believed to dominate. The second phase of the study presented here focused on improving the analysis methods that are available with LPDMs. While LPDMs are popular and useful for the analysis of atmospheric trace gas measurements, identifying the transport pathway of emissions from their source to a receptor (the Pico Mountain observatory in our case) using the standard gridded model output, particularly during complex meteorological scenarios can be difficult can be difficult or impossible. The transport study in phase 1 was limited to only 1 month out of more than 3 years of available data and included only 4 case studies out of the 16 events specifically due to this confounding factor. The second phase of this study addressed this difficulty by presenting a method to clearly and easily identify the pathway taken by only those emissions that arrive at a receptor at a particular time, by combining the standard gridded output from forward (i.e., concentrations) and backward (i.e., residence time) LPDM simulations, greatly simplifying similar analyses. The ability of the method to successfully determine the source-to-receptor pathway, restoring this Lagrangian information that is lost when the data are gridded, is proven by comparing the pathway determined from this method with the particle trajectories from both the forward and backward models. A sample analysis is also presented, demonstrating that this method is more accurate and easier to use than existing methods using standard LPDM products. Finally, we discuss potential future work that would be possible by combining the backward LPDM simulation with gridded data from other sources (e.g., chemical transport models) to obtain a Lagrangian sampling of the air that will eventually arrive at a receptor.

The Distribution of the Irreducibles in an Algebraic Number Field

The objective of this thesis is to study the distribution of the number of principal ideals generated by an irreducible element in an algebraic number field, namely in the non-unique factorization ring of integers of such a field. In particular we are investigating the size of M(x), defined as M ( x ) =∑ (α) α irred.|N (α)|≤≠ 1, where x is any positive real number and N (α) is the norm of α. We finally obtain asymptotic results for hl(x).

Velocities and elastic moduli of basalts from Pacific DSDP sites

Compressional (Vp) and shear (Vs) wave velocities have been measured to 10 kb in 32 cores of basalt from 14 Pacific sites of the Deep Sea Drilling Project. Both Vp and V s show wide ranges (3.70 to 6.38 km/sec for Vp and 1.77 to 3.40 km/sec for V s at 0.5 kb) which are linearly related to density and sea floor age, confirming earlier findings by Christensen and Salisbury of decreasing velocity with progressive submarine weathering based on studies of basalts from five sites in the Atlantic. Combined Pacific and Atlantic data give rates of decreasing velocity of -1.89 and -1.35 km/sec per 100 my for Vp and Vs respectively. New analyses of oceanic seismic refraction data indicate a decrease in layer 2 velocities with age similar to that observed in the laboratory, suggesting that weathering penetrates to several hundred meters in many regions and is largely responsible for the extreme range and variability of layer 2 refraction velocities.

Reduced order models for industrial thermo-fluid problems

A new method is presented to generate reduced order models (ROMs) in Fluid Dynamics problems of industrial interest. The method is based on the expansion of the flow variables in a Proper Orthogonal Decomposition (POD) basis, calculated from a limited number of snapshots, which are obtained via Computational Fluid Dynamics (CFD). Then, the POD-mode amplitudes are calculated as minimizers of a properly defined overall residual of the equations and boundary conditions. The method includes various ingredients that are new in this field. The residual can be calculated using only a limited number of points in the flow field, which can be scattered either all over the whole computational domain or over a smaller projection window. The resulting ROM is both computationally efficient(reconstructed flow fields require, in cases that do not present shock waves, less than 1 % of the time needed to compute a full CFD solution) and flexible(the projection window can avoid regions of large localized CFD errors).Also, for problems related with aerodynamics, POD modes are obtained from a set of snapshots calculated by a CFD method based on the compressible Navier Stokes equations and a turbulence model (which further more includes some unphysical stabilizing terms that are included for purely numerical reasons), but projection onto the POD manifold is made using the inviscid Euler equations, which makes the method independent of the CFD scheme. In addition, shock waves are treated specifically in the POD description, to avoid the need of using a too large number of snapshots. Various definitions of the residual are also discussed, along with the number and distribution of snapshots, the number of retained modes, and the effect of CFD errors. The method is checked and discussed on several test problems that describe (i) heat transfer in the recirculation region downstream of a backwards facing step, (ii) the flow past a two-dimensional airfoil in both the subsonic and transonic regimes, and (iii) the flow past a three-dimensional horizontal tail plane. The method is both efficient and numerically robust in the sense that the computational effort is quite small compared to CFD and results are both reasonably accurate and largely insensitive to the definition of the residual, to CFD errors, and to the CFD method itself, which may contain artificial stabilizing terms. Thus, the method is amenable for practical engineering applications. Resumen Se presenta un nuevo método para generar modelos de orden reducido (ROMs) aplicado a problemas fluidodinámicos de interés industrial. El nuevo método se basa en la expansión de las variables fluidas en una base POD, calculada a partir de un cierto número de snapshots, los cuales se han obtenido gracias a simulaciones numéricas (CFD). A continuación, las amplitudes de los modos POD se calculan minimizando un residual global adecuadamente definido que combina las ecuaciones y las condiciones de contorno. El método incluye varios ingredientes que son nuevos en este campo de estudio. El residual puede calcularse utilizando únicamente un número limitado de puntos del campo fluido. Estos puntos puede encontrarse dispersos a lo largo del dominio computacional completo o sobre una ventana de proyección. El modelo ROM obtenido es tanto computacionalmente eficiente (en aquellos casos que no presentan ondas de choque reconstruir los campos fluidos requiere menos del 1% del tiempo necesario para calcular una solución CFD) como flexible (la ventana de proyección puede escogerse de forma que evite contener regiones con errores en la solución CFD localizados y grandes). Además, en problemas aerodinámicos, los modos POD se obtienen de un conjunto de snapshots calculados utilizando un código CFD basado en la versión compresible de las ecuaciones de Navier Stokes y un modelo de turbulencia (el cual puede incluir algunos términos estabilizadores sin sentido físico que se añaden por razones puramente numéricas), aunque la proyección en la variedad POD se hace utilizando las ecuaciones de Euler, lo que hace al método independiente del esquema utilizado en el código CFD. Además, las ondas de choque se tratan específicamente en la descripción POD para evitar la necesidad de utilizar un número demasiado grande de snapshots. Varias definiciones del residual se discuten, así como el número y distribución de los snapshots,el número de modos retenidos y el efecto de los errores debidos al CFD. El método se comprueba y discute para varios problemas de evaluación que describen (i) la transferencia de calor en la región de recirculación aguas abajo de un escalón, (ii) el flujo alrededor de un perfil bidimensional en regímenes subsónico y transónico y (iii) el flujo alrededor de un estabilizador horizontal tridimensional. El método es tanto eficiente como numéricamente robusto en el sentido de que el esfuerzo computacional es muy pequeño comparado con el requerido por el CFD y los resultados son razonablemente precisos y muy insensibles a la definición del residual, los errores debidos al CFD y al método CFD en sí mismo, el cual puede contener términos estabilizadores artificiales. Por lo tanto, el método puede utilizarse en aplicaciones prácticas de ingeniería.

Examples of PDEs all whose points are characteristic

Let π : FM ! M be the bundle of linear frames of a manifold M. A basis Lijk , j < k, of diffeomorphism invariant Lagrangians on J1 (FM) was determined in [J. Muñoz Masqué, M. E. Rosado, Invariant variational problems on linear frame bundles, J. Phys. A35 (2002) 2013-2036]. The notion of a characteristic hypersurface for an arbitrary first-order PDE system on an ar- bitrary bred manifold π : P → M, is introduced and for the systems dened by the Euler-Lagrange equations of Lijk every hypersurface is shown to be characteristic. The Euler-Lagrange equations of the natural basis of Lagrangian densities Lijk on the bundle of linear frames of a manifold M which are invariant under diffeomorphisms, are shown to be an underdetermined PDEs systems such that every hypersurface of M is characteristic for such equations. This explains why these systems cannot be written in the Cauchy-Kowaleska form, although they are known to be formally integrable by using the tools of geometric theory of partial differential equations, see [J. Muñoz Masqué, M. E. Rosado, Integrability of the eld equations of invariant variational problems on linear frame bundles, J. Geom. Phys. 49 (2004), 119-155]