Earth to Moon Transfer: Direct vs Via Libration Points (L1, L2)

NASA Technical Reports Server (NTRS)

Condon, Gerald L.; Wilson, Samuel W.

2004-01-01

For some three decades, the Apollo-style mission has served as a proven baseline technique for transporting flight crews to the Moon and back with expendable hardware. This approach provides an optimal design for expeditionary missions, emphasizing operational flexibility in terms of safely returning the crew in the event of a hardware failure. However, its application is limited essentially to low-latitude lunar sites, and it leaves much to be desired as a model for exploratory and evolutionary programs that employ reusable space-based hardware. This study compares the performance requirements for a lunar orbit rendezvous mission type with one using the cislunar libration point (L1) as a stopover and staging point for access to arbitrary sites on the lunar surface. For selected constraints and mission objectives, it contrasts the relative uniformity of performance cost when the L1 staging point is used with the wide variation of cost for the Apollo-style lunar orbit rendezvous.

Communication: A simplified coupled-cluster Lagrangian for polarizable embedding.

PubMed

Krause, Katharina; Klopper, Wim

2016-01-28

A simplified coupled-cluster Lagrangian, which is linear in the Lagrangian multipliers, is proposed for the coupled-cluster treatment of a quantum mechanical system in a polarizable environment. In the simplified approach, the amplitude equations are decoupled from the Lagrangian multipliers and the energy obtained from the projected coupled-cluster equation corresponds to a stationary point of the Lagrangian.

Communication: A simplified coupled-cluster Lagrangian for polarizable embedding

DOE Office of Scientific and Technical Information (OSTI.GOV)

Krause, Katharina; Klopper, Wim, E-mail: klopper@kit.edu

A simplified coupled-cluster Lagrangian, which is linear in the Lagrangian multipliers, is proposed for the coupled-cluster treatment of a quantum mechanical system in a polarizable environment. In the simplified approach, the amplitude equations are decoupled from the Lagrangian multipliers and the energy obtained from the projected coupled-cluster equation corresponds to a stationary point of the Lagrangian.

An extended Lagrangian method

NASA Technical Reports Server (NTRS)

Liou, Meng-Sing

1993-01-01

A unique formulation of describing fluid motion is presented. The method, referred to as 'extended Lagrangian method', is interesting from both theoretical and numerical points of view. The formulation offers accuracy in numerical solution by avoiding numerical diffusion resulting from mixing of fluxes in the Eulerian description. Meanwhile, it also avoids the inaccuracy incurred due to geometry and variable interpolations used by the previous Lagrangian methods. The present method is general and capable of treating subsonic flows as well as supersonic flows. The method proposed in this paper is robust and stable. It automatically adapts to flow features without resorting to clustering, thereby maintaining rather uniform grid spacing throughout and large time step. Moreover, the method is shown to resolve multidimensional discontinuities with a high level of accuracy, similar to that found in 1D problems.

The motion near L{sub 4} equilibrium point under non-point mass primaries

DOE Office of Scientific and Technical Information (OSTI.GOV)

Huda, I. N., E-mail: ibnu.nurul@students.itb.ac.id; Utama, J. A.; Madley, D.

2015-09-30

The Circular Restricted Three-Body Problem (CRTBP) possesses five equilibrium points, that comprise three collinear (L{sub 1}, L{sub 2}, and L{sub 3}) and two triangular points (L{sub 4} and L{sub 5}). The classical study (with the primaries are point mass) suggests that the equilibrium points may cause the velocity of infinitesimal object relatively becomes zero and reveals the zero velocity curve. We study the motion of infinitesimal object near triangular equilibrium point (L{sub 4}) and determine its zero velocity curve. We extend the study by taking into account the effects of radiation of the bigger primary (q{sub 1} ≠ 1, q{submore » 2} = 1) and oblateness of the smaller primary (A{sub 1} = 0, A{sub 2} ≠ 0). The location of L{sub 4} is analytically derived then the stability of L{sub 4} and its zero velocity curves are studied numerically. Our study suggests that the oblateness and the radiation of primaries may affect the stability and zero velocity curve around L{sub 4}.« less

Lagrangian postprocessing of computational hemodynamics.

PubMed

Shadden, Shawn C; Arzani, Amirhossein

2015-01-01

Recent advances in imaging, modeling, and computing have rapidly expanded our capabilities to model hemodynamics in the large vessels (heart, arteries, and veins). This data encodes a wealth of information that is often under-utilized. Modeling (and measuring) blood flow in the large vessels typically amounts to solving for the time-varying velocity field in a region of interest. Flow in the heart and larger arteries is often complex, and velocity field data provides a starting point for investigating the hemodynamics. This data can be used to perform Lagrangian particle tracking, and other Lagrangian-based postprocessing. As described herein, Lagrangian methods are necessary to understand inherently transient hemodynamic conditions from the fluid mechanics perspective, and to properly understand the biomechanical factors that lead to acute and gradual changes of vascular function and health. The goal of the present paper is to review Lagrangian methods that have been used in post-processing velocity data of cardiovascular flows.

Lagrangian postprocessing of computational hemodynamics

PubMed Central

Shadden, Shawn C.; Arzani, Amirhossein

2014-01-01

Recent advances in imaging, modeling and computing have rapidly expanded our capabilities to model hemodynamics in the large vessels (heart, arteries and veins). This data encodes a wealth of information that is often under-utilized. Modeling (and measuring) blood flow in the large vessels typically amounts to solving for the time-varying velocity field in a region of interest. Flow in the heart and larger arteries is often complex, and velocity field data provides a starting point for investigating the hemodynamics. This data can be used to perform Lagrangian particle tracking, and other Lagrangian-based postprocessing. As described herein, Lagrangian methods are necessary to understand inherently transient hemodynamic conditions from the fluid mechanics perspective, and to properly understand the biomechanical factors that lead to acute and gradual changes of vascular function and health. The goal of the present paper is to review Lagrangian methods that have been used in post-processing velocity data of cardiovascular flows. PMID:25059889

Lagrangian space consistency relation for large scale structure

DOE Office of Scientific and Technical Information (OSTI.GOV)

Horn, Bart; Hui, Lam; Xiao, Xiao

Consistency relations, which relate the squeezed limit of an (N+1)-point correlation function to an N-point function, are non-perturbative symmetry statements that hold even if the associated high momentum modes are deep in the nonlinear regime and astrophysically complex. Recently, Kehagias & Riotto and Peloso & Pietroni discovered a consistency relation applicable to large scale structure. We show that this can be recast into a simple physical statement in Lagrangian space: that the squeezed correlation function (suitably normalized) vanishes. This holds regardless of whether the correlation observables are at the same time or not, and regardless of whether multiple-streaming is present.more » Furthermore, the simplicity of this statement suggests that an analytic understanding of large scale structure in the nonlinear regime may be particularly promising in Lagrangian space.« less

Lagrangian space consistency relation for large scale structure

DOE PAGES

Horn, Bart; Hui, Lam; Xiao, Xiao

2015-09-29

Consistency relations, which relate the squeezed limit of an (N+1)-point correlation function to an N-point function, are non-perturbative symmetry statements that hold even if the associated high momentum modes are deep in the nonlinear regime and astrophysically complex. Recently, Kehagias & Riotto and Peloso & Pietroni discovered a consistency relation applicable to large scale structure. We show that this can be recast into a simple physical statement in Lagrangian space: that the squeezed correlation function (suitably normalized) vanishes. This holds regardless of whether the correlation observables are at the same time or not, and regardless of whether multiple-streaming is present.more » Furthermore, the simplicity of this statement suggests that an analytic understanding of large scale structure in the nonlinear regime may be particularly promising in Lagrangian space.« less

An extended Lagrangian method

NASA Technical Reports Server (NTRS)

Liou, Meng-Sing

1992-01-01

A unique formulation of describing fluid motion is presented. The method, referred to as 'extended Lagrangian method', is interesting from both theoretical and numerical points of view. The formulation offers accuracy in numerical solution by avoiding numerical diffusion resulting from mixing of fluxes in the Eulerian description. Meanwhile, it also avoids the inaccuracy incurred due to geometry and variable interpolations used by the previous Lagrangian methods. Unlike the Lagrangian method previously imposed which is valid only for supersonic flows, the present method is general and capable of treating subsonic flows as well as supersonic flows. The method proposed in this paper is robust and stable. It automatically adapts to flow features without resorting to clustering, thereby maintaining rather uniform grid spacing throughout and large time step. Moreover, the method is shown to resolve multi-dimensional discontinuities with a high level of accuracy, similar to that found in one-dimensional problems.

The Lagrangian Points

ERIC Educational Resources Information Center

Linton, J. Oliver

2017-01-01

There are five unique points in a star/planet system where a satellite can be placed whose orbital period is equal to that of the planet. Simple methods for calculating the positions of these points, or at least justifying their existence, are developed.

Extended hamiltonian formalism and Lorentz-violating lagrangians

NASA Astrophysics Data System (ADS)

Colladay, Don

2017-09-01

A new perspective on the classical mechanical formulation of particle trajectories in Lorentz-violating theories is presented. Using the extended hamiltonian formalism, a Legendre Transformation between the associated covariant lagrangian and hamiltonian varieties is constructed. This approach enables calculation of trajectories using Hamilton's equations in momentum space and the Euler-Lagrange equations in velocity space away from certain singular points that arise in the theory. Singular points are naturally de-singularized by requiring the trajectories to be smooth functions of both velocity and momentum variables. In addition, it is possible to identify specific sheets of the dispersion relations that correspond to specific solutions for the lagrangian. Examples corresponding to bipartite Finsler functions are computed in detail. A direct connection between the lagrangians and the field-theoretic solutions to the Dirac equation is also established for a special case.

Uncertainty quantification in Eulerian-Lagrangian models for particle-laden flows

NASA Astrophysics Data System (ADS)

Fountoulakis, Vasileios; Jacobs, Gustaaf; Udaykumar, Hs

2017-11-01

A common approach to ameliorate the computational burden in simulations of particle-laden flows is to use a point-particle based Eulerian-Lagrangian model, which traces individual particles in their Lagrangian frame and models particles as mathematical points. The particle motion is determined by Stokes drag law, which is empirically corrected for Reynolds number, Mach number and other parameters. The empirical corrections are subject to uncertainty. Treating them as random variables renders the coupled system of PDEs and ODEs stochastic. An approach to quantify the propagation of this parametric uncertainty to the particle solution variables is proposed. The approach is based on averaging of the governing equations and allows for estimation of the first moments of the quantities of interest. We demonstrate the feasibility of our proposed methodology of uncertainty quantification of particle-laden flows on one-dimensional linear and nonlinear Eulerian-Lagrangian systems. This research is supported by AFOSR under Grant FA9550-16-1-0008.

Lagrangian space consistency relation for large scale structure

DOE Office of Scientific and Technical Information (OSTI.GOV)

Horn, Bart; Hui, Lam; Xiao, Xiao, E-mail: bh2478@columbia.edu, E-mail: lh399@columbia.edu, E-mail: xx2146@columbia.edu

Consistency relations, which relate the squeezed limit of an (N+1)-point correlation function to an N-point function, are non-perturbative symmetry statements that hold even if the associated high momentum modes are deep in the nonlinear regime and astrophysically complex. Recently, Kehagias and Riotto and Peloso and Pietroni discovered a consistency relation applicable to large scale structure. We show that this can be recast into a simple physical statement in Lagrangian space: that the squeezed correlation function (suitably normalized) vanishes. This holds regardless of whether the correlation observables are at the same time or not, and regardless of whether multiple-streaming is present.more » The simplicity of this statement suggests that an analytic understanding of large scale structure in the nonlinear regime may be particularly promising in Lagrangian space.« less

Mean-Lagrangian formalism and covariance of fluid turbulence.

PubMed

Ariki, Taketo

2017-05-01

Mean-field-based Lagrangian framework is developed for the fluid turbulence theory, which enables physically objective discussions, especially, of the history effect. Mean flow serves as a purely geometrical object of Lie group theory, providing useful operations to measure the objective rate and history integration of the general tensor field. The proposed framework is applied, on the one hand, to one-point closure model, yielding an objective expression of the turbulence viscoelastic effect. Application to two-point closure, on the other hand, is also discussed, where natural extension of known Lagrangian correlation is discovered on the basis of an extended covariance group.

Stability of triangular lagrangian points in elliptical restricted three body problem under the radiating binary systems

NASA Astrophysics Data System (ADS)

Narayan, A.; Singh, Nutan

2014-10-01

This paper studies the stability of Triangular Lagrangian points in the model of elliptical restricted three body problem, under the assumption that both the primaries are radiating. The model proposed is applicable to the well known binary systems Achird, Luyten, αCen AB, Kruger-60, Xi-Bootis. Conditional stability of the motion around the triangular points exists for 0≤ μ≤ μ ∗, where μ is the mass ratio. The method of averaging due to Grebenikov has been exploited throughout the analysis of stability of the system. The critical mass ratio depends on the combined effects of radiation of both the primaries and eccentricity of this orbit. It is found by adopting the simulation technique that the range of stability decreases as the radiation pressure parameter increases.

Fingerprints of heavy scales in electroweak effective Lagrangians

NASA Astrophysics Data System (ADS)

Pich, Antonio; Rosell, Ignasi; Santos, Joaquín; Sanz-Cillero, Juan José

2017-04-01

The couplings of the electroweak effective theory contain information on the heavy-mass scales which are no-longer present in the low-energy Lagrangian. We build a general effective Lagrangian, implementing the electroweak chiral symmetry breaking SU(2) L ⊗ SU(2) R → SU(2) L+ R , which couples the known particle fields to heavier states with bosonic quantum numbers J P = 0± and 1±. We consider colour-singlet heavy fields that are in singlet or triplet representations of the electroweak group. Integrating out these heavy scales, we analyze the pattern of low-energy couplings among the light fields which are generated by the massive states. We adopt a generic non-linear realization of the electroweak symmetry breaking with a singlet Higgs, without making any assumption about its possible doublet structure. Special attention is given to the different possible descriptions of massive spin-1 fields and the differences arising from naive implementations of these formalisms, showing their full equivalence once a proper short-distance behaviour is required.

Terrestrial glint seen from deep space: Oriented ice crystals detected from the Lagrangian point

NASA Astrophysics Data System (ADS)

Marshak, Alexander; Várnai, Tamás.; Kostinski, Alexander

2017-05-01

The Deep Space Climate Observatory (DSCOVR) spacecraft resides at the first Lagrangian point about one million miles from Earth. A polychromatic imaging camera onboard delivers nearly hourly observations of the entire sunlit face of the Earth. Many images contain unexpected bright flashes of light over both ocean and land. We construct a yearlong time series of flash latitudes, scattering angles, and oxygen absorption to demonstrate conclusively that the flashes over land are specular reflections off tiny ice platelets floating in the air nearly horizontally. Such deep space detection of tropospheric ice can be used to constrain the likelihood of oriented crystals and their contribution to Earth albedo. These glint observations also support proposals for detecting starlight glints off faint companions in our search for habitable exoplanets.

Planetary Defense From Space: Part 1-Keplerian Theory

NASA Astrophysics Data System (ADS)

Maccone, Claudio

A system of two space bases housing missiles is proposed to achieve the Planetary Defense of the Earth against dangerous asteroids and comets. We show that the layout of the Earth-Moon system with the five relevant Lagrangian (or libration) points in space leads naturally to only one, unmistakable location of these two space bases within the sphere of influence of the Earth. These locations are at the two Lagrangian points L1 (in between the Earth and the Moon) and L3 (in the direction opposite to the Moon from the Earth). We show that placing bases of missiles at L1 and L3 would cause those missiles to deflect the trajectory of asteroids by hitting them orthogonally to their impact trajectory toward the Earth, so as to maximize their deflection. We show that the confocal conics are the best class of trajectories fulfilling this orthogonal deflection requirement. An additional remark is that the theory developed in this paper is just a beginning of a larger set of future research work. In fact, while in this paper we only develop the Keplerian analytical theory of the Optimal Planetary Defense achievable from the Earth-Moon Lagrangian points L1 and L3, much more sophisticated analytical refinements would be needed to: Take into account many perturbation forces of all kinds acting on both the asteroids and missiles shot from L1 and L3; add more (non-optimal) trajectories of missiles shot from either the Lagrangian points L4 and L5 of the Earth-Moon system or from the surface of the Moon itself; encompass the full range of missiles currently available to the US (and possibly other countries) so as to really see "which asteroids could be diverted by which missiles", even in the very simplified scheme outlined here. Outlined for the first time in February 2002, our Confocal Planetary Defense concept is a Keplerian Theory that proved simple enough to catch the attention of scholars, representatives of the US Military and popular writers. These developments could

A Theoretical Framework for Lagrangian Descriptors

NASA Astrophysics Data System (ADS)

Lopesino, C.; Balibrea-Iniesta, F.; García-Garrido, V. J.; Wiggins, S.; Mancho, A. M.

This paper provides a theoretical background for Lagrangian Descriptors (LDs). The goal of achieving rigorous proofs that justify the ability of LDs to detect invariant manifolds is simplified by introducing an alternative definition for LDs. The definition is stated for n-dimensional systems with general time dependence, however we rigorously prove that this method reveals the stable and unstable manifolds of hyperbolic points in four particular 2D cases: a hyperbolic saddle point for linear autonomous systems, a hyperbolic saddle point for nonlinear autonomous systems, a hyperbolic saddle point for linear nonautonomous systems and a hyperbolic saddle point for nonlinear nonautonomous systems. We also discuss further rigorous results which show the ability of LDs to highlight additional invariants sets, such as n-tori. These results are just a simple extension of the ergodic partition theory which we illustrate by applying this methodology to well-known examples, such as the planar field of the harmonic oscillator and the 3D ABC flow. Finally, we provide a thorough discussion on the requirement of the objectivity (frame-invariance) property for tools designed to reveal phase space structures and their implications for Lagrangian descriptors.

N =1 Lagrangians for generalized Argyres-Douglas theories

NASA Astrophysics Data System (ADS)

Agarwal, Prarit; Sciarappa, Antonio; Song, Jaewon

2017-10-01

We find N = 1 Lagrangian gauge theories that flow to generalized ArgyresDouglas theories with N = 2 supersymmetry. We find that certain SU quiver gauge theories flow to generalized Argyres-Douglas theories of type ( A k-1 , A mk-1) and ( I m,km , S). We also find quiver gauge theories of SO/Sp gauge groups flowing to the ( A 2 m-1 , D 2 mk+1), ( A 2 m , D 2 m( k-1)+ k ) and D m(2 k + 2) m(2 k + 2) [ m] theories.

"Lagrangian" for a Non-Lagrangian Field Theory with N=2 Supersymmetry.

PubMed

Gadde, Abhijit; Razamat, Shlomo S; Willett, Brian

2015-10-23

We suggest that at least some of the strongly coupled N=2 quantum field theories in 4D can have a nonconformal N=1 Lagrangian description flowing to them at low energies. In particular, we construct such a description for the N=2 rank one superconformal field theory with E(6) flavor symmetry, for which a Lagrangian description was previously unavailable. We utilize this description to compute several supersymmetric partition functions.

Science objectives of the magnetic field experiment onboard Aditya-L1 spacecraft

NASA Astrophysics Data System (ADS)

Yadav, Vipin K.; Srivastava, Nandita; Ghosh, S. S.; Srikar, P. T.; Subhalakshmi, Krishnamoorthy

2018-01-01

The Aditya-L1 is first Indian solar mission scheduled to be placed in a halo orbit around the first Lagrangian point (L1) of Sun-Earth system in the year 2018-19. The approved scientific payloads onboard Aditya-L1 spacecraft includes a Fluxgate Digital Magnetometer (FGM) to measure the local magnetic field which is necessary to supplement the outcome of other scientific experiments onboard. The in-situ vector magnetic field data at L1 is essential for better understanding of the data provided by the particle and plasma analysis experiments, onboard Aditya-L1 mission. Also, the dynamics of Coronal Mass Ejections (CMEs) can be better understood with the help of in-situ magnetic field data at the L1 point region. This data will also serve as crucial input for the short lead-time space weather forecasting models. The proposed FGM is a dual range magnetic sensor on a 6 m long boom mounted on the Sun viewing panel deck and configured to deploy along the negative roll direction of the spacecraft. Two sets of sensors (tri-axial each) are proposed to be mounted, one at the tip of boom (6 m from the spacecraft) and other, midway (3 m from the spacecraft). The main science objective of this experiment is to measure the magnitude and nature of the interplanetary magnetic field (IMF) locally and to study the disturbed magnetic conditions and extreme solar events by detecting the CME from Sun as a transient event. The proposed secondary science objectives are to study the impact of interplanetary structures and shock solar wind interaction on geo-space environment and to detect low frequency plasma waves emanating from the solar corona at L1 point. This will provide a better understanding on how the Sun affects interplanetary space. In this paper, we shall give the main scientific objectives of the magnetic field experiment and brief technical details of the FGM onboard Aditya-1 spacecraft.

Leading-order classical Lagrangians for the nonminimal standard-model extension

NASA Astrophysics Data System (ADS)

Reis, J. A. A. S.; Schreck, M.

2018-03-01

In this paper, we derive the general leading-order classical Lagrangian covering all fermion operators of the nonminimal standard-model extension (SME). Such a Lagrangian is considered to be the point-particle analog of the effective field theory description of Lorentz violation that is provided by the SME. At leading order in Lorentz violation, the Lagrangian obtained satisfies the set of five nonlinear equations that govern the map from the field theory to the classical description. This result can be of use for phenomenological studies of classical bodies in gravitational fields.

Lagrangian statistics in compressible isotropic homogeneous turbulence

NASA Astrophysics Data System (ADS)

Yang, Yantao; Wang, Jianchun; Shi, Yipeng; Chen, Shiyi

2011-11-01

In this work we conducted the Direct Numerical Simulation (DNS) of a forced compressible isotropic homogeneous turbulence and investigated the flow statistics from the Lagrangian point of view, namely the statistics is computed following the passive tracers trajectories. The numerical method combined the Eulerian field solver which was developed by Wang et al. (2010, J. Comp. Phys., 229, 5257-5279), and a Lagrangian module for tracking the tracers and recording the data. The Lagrangian probability density functions (p.d.f.'s) have then been calculated for both kinetic and thermodynamic quantities. In order to isolate the shearing part from the compressing part of the flow, we employed the Helmholtz decomposition to decompose the flow field (mainly the velocity field) into the solenoidal and compressive parts. The solenoidal part was compared with the incompressible case, while the compressibility effect showed up in the compressive part. The Lagrangian structure functions and cross-correlation between various quantities will also be discussed. This work was supported in part by the China's Turbulence Program under Grant No.2009CB724101.

Multi-Lagrangians for integrable systems

NASA Astrophysics Data System (ADS)

Nutku, Y.; Pavlov, M. V.

2002-03-01

We propose a general scheme to construct multiple Lagrangians for completely integrable nonlinear evolution equations that admit multi-Hamiltonian structure. The recursion operator plays a fundamental role in this construction. We use a conserved quantity higher/lower than the Hamiltonian in the potential part of the new Lagrangian and determine the corresponding kinetic terms by generating the appropriate momentum map. This leads to some remarkable new developments. We show that nonlinear evolutionary systems that admit N-fold first order local Hamiltonian structure can be cast into variational form with 2N-1 Lagrangians which will be local functionals of Clebsch potentials. This number increases to 3N-2 when the Miura transformation is invertible. Furthermore we construct a new Lagrangian for polytropic gas dynamics in 1+1 dimensions which is a free, local functional of the physical field variables, namely density and velocity, thus dispensing with the necessity of introducing Clebsch potentials entirely. This is a consequence of bi-Hamiltonian structure with a compatible pair of first and third order Hamiltonian operators derived from Sheftel's recursion operator.

Some Lagrangians for systems without a Lagrangian

NASA Astrophysics Data System (ADS)

Nucci, M. C.; Leach, P. G. L.

2011-03-01

We demonstrate how to construct many different Lagrangians for two famous examples that were deemed by Douglas (1941 Trans. Am. Math. Soc. 50 71-128) not to have a Lagrangian. Following Bateman's dictum (1931 Phys. Rev. 38 815-9), we determine different sets of equations that are compatible with those of Douglas and derivable from a variational principle.

Power corrections to the HTL effective Lagrangian of QED

NASA Astrophysics Data System (ADS)

Carignano, Stefano; Manuel, Cristina; Soto, Joan

2018-05-01

We present compact expressions for the power corrections to the hard thermal loop (HTL) Lagrangian of QED in d space dimensions. These are corrections of order (L / T) 2, valid for momenta L ≪ T, where T is the temperature. In the limit d → 3 we achieve a consistent regularization of both infrared and ultraviolet divergences, which respects the gauge symmetry of the theory. Dimensional regularization also allows us to witness subtle cancellations of infrared divergences. We also discuss how to generalize our results in the presence of a chemical potential, so as to obtain the power corrections to the hard dense loop (HDL) Lagrangian.

Forecasting Future Sea Ice Conditions: A Lagrangian Approach

DTIC Science & Technology

2015-09-30

1 DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Forecasting Future Sea Ice Conditions: A Lagrangian ...GCMs participating in IPCC AR5 agree with observed source region patterns from the satellite- derived dataset. 4- Compare Lagrangian ice... Lagrangian sea-ice back trajectories to estimate thermodynamic and dynamic (advection) ice loss. APPROACH We use a Lagrangian trajectory model to

Lagrangian-Hamiltonian unified formalism for autonomous higher order dynamical systems

NASA Astrophysics Data System (ADS)

Prieto-Martínez, Pedro Daniel; Román-Roy, Narciso

2011-09-01

The Lagrangian-Hamiltonian unified formalism of Skinner and Rusk was originally stated for autonomous dynamical systems in classical mechanics. It has been generalized for non-autonomous first-order mechanical systems, as well as for first-order and higher order field theories. However, a complete generalization to higher order mechanical systems is yet to be described. In this work, after reviewing the natural geometrical setting and the Lagrangian and Hamiltonian formalisms for higher order autonomous mechanical systems, we develop a complete generalization of the Lagrangian-Hamiltonian unified formalism for these kinds of systems, and we use it to analyze some physical models from this new point of view.

Point-spread function reconstruction in ground-based astronomy by l(1)-l(p) model.

PubMed

Chan, Raymond H; Yuan, Xiaoming; Zhang, Wenxing

2012-11-01

In ground-based astronomy, images of objects in outer space are acquired via ground-based telescopes. However, the imaging system is generally interfered by atmospheric turbulence, and hence images so acquired are blurred with unknown point-spread function (PSF). To restore the observed images, the wavefront of light at the telescope's aperture is utilized to derive the PSF. A model with the Tikhonov regularization has been proposed to find the high-resolution phase gradients by solving a least-squares system. Here we propose the l(1)-l(p) (p=1, 2) model for reconstructing the phase gradients. This model can provide sharper edges in the gradients while removing noise. The minimization models can easily be solved by the Douglas-Rachford alternating direction method of a multiplier, and the convergence rate is readily established. Numerical results are given to illustrate that the model can give better phase gradients and hence a more accurate PSF. As a result, the restored images are much more accurate when compared to the traditional Tikhonov regularization model.

Lagrangian averages, averaged Lagrangians, and the mean effects of fluctuations in fluid dynamics.

PubMed

Holm, Darryl D.

2002-06-01

We begin by placing the generalized Lagrangian mean (GLM) equations for a compressible adiabatic fluid into the Euler-Poincare (EP) variational framework of fluid dynamics, for an averaged Lagrangian. This is the Lagrangian averaged Euler-Poincare (LAEP) theorem. Next, we derive a set of approximate small amplitude GLM equations (glm equations) at second order in the fluctuating displacement of a Lagrangian trajectory from its mean position. These equations express the linear and nonlinear back-reaction effects on the Eulerian mean fluid quantities by the fluctuating displacements of the Lagrangian trajectories in terms of their Eulerian second moments. The derivation of the glm equations uses the linearized relations between Eulerian and Lagrangian fluctuations, in the tradition of Lagrangian stability analysis for fluids. The glm derivation also uses the method of averaged Lagrangians, in the tradition of wave, mean flow interaction. Next, the new glm EP motion equations for incompressible ideal fluids are compared with the Euler-alpha turbulence closure equations. An alpha model is a GLM (or glm) fluid theory with a Taylor hypothesis closure. Such closures are based on the linearized fluctuation relations that determine the dynamics of the Lagrangian statistical quantities in the Euler-alpha equations. Thus, by using the LAEP theorem, we bridge between the GLM equations and the Euler-alpha closure equations, through the small-amplitude glm approximation in the EP variational framework. We conclude by highlighting a new application of the GLM, glm, and alpha-model results for Lagrangian averaged ideal magnetohydrodynamics. (c) 2002 American Institute of Physics.

Long Term Missions at the Sun-Earth Libration Point L1: ACE, SOHO, and WIND

NASA Technical Reports Server (NTRS)

Roberts, Craig E.

2011-01-01

Three heliophysics missions - the Solar Heliospheric Observatory (SOHO), the Advanced Composition Explorer (ACE), and the Global Geoscience WIND - have been orbiting the Sun-Earth interior libration point L1 continuously since 1996, 1997, and 2004, respectively. ACE and WIND (both NASA missions) and SOHO (an ESA-NASA joint mission) are all operated from the NASA Goddard Space Flight Center Flight Dynamics Facility. While ACE and SOHO have been dedicated libration point orbiters since their launches, WIND prior to 2004 flew a remarkable 10-year deep-space trajectory that featured 38 targeted lunar flybys. The L1 orbits and the mission histories of the three spacecraft are briefly reviewed, and the station-keeping techniques and orbit maneuver experience are discussed.

Learn the Lagrangian: A Vector-Valued RKHS Approach to Identifying Lagrangian Systems.

PubMed

Cheng, Ching-An; Huang, Han-Pang

2016-12-01

We study the modeling of Lagrangian systems with multiple degrees of freedom. Based on system dynamics, canonical parametric models require ad hoc derivations and sometimes simplification for a computable solution; on the other hand, due to the lack of prior knowledge in the system's structure, modern nonparametric models in machine learning face the curse of dimensionality, especially in learning large systems. In this paper, we bridge this gap by unifying the theories of Lagrangian systems and vector-valued reproducing kernel Hilbert space. We reformulate Lagrangian systems with kernels that embed the governing Euler-Lagrange equation-the Lagrangian kernels-and show that these kernels span a subspace capturing the Lagrangian's projection as inverse dynamics. By such property, our model uses only inputs and outputs as in machine learning and inherits the structured form as in system dynamics, thereby removing the need for the mundane derivations for new systems as well as the generalization problem in learning from scratches. In effect, it learns the system's Lagrangian, a simpler task than directly learning the dynamics. To demonstrate, we applied the proposed kernel to identify the robot inverse dynamics in simulations and experiments. Our results present a competitive novel approach to identifying Lagrangian systems, despite using only inputs and outputs.

Lagrangian statistics in weakly forced two-dimensional turbulence.

PubMed

Rivera, Michael K; Ecke, Robert E

2016-01-01

Measurements of Lagrangian single-point and multiple-point statistics in a quasi-two-dimensional stratified layer system are reported. The system consists of a layer of salt water over an immiscible layer of Fluorinert and is forced electromagnetically so that mean-squared vorticity is injected at a well-defined spatial scale ri. Simultaneous cascades develop in which enstrophy flows predominately to small scales whereas energy cascades, on average, to larger scales. Lagrangian correlations and one- and two-point displacements are measured for random initial conditions and for initial positions within topological centers and saddles. Some of the behavior of these quantities can be understood in terms of the trapping characteristics of long-lived centers, the slow motion near strong saddles, and the rapid fluctuations outside of either centers or saddles. We also present statistics of Lagrangian velocity fluctuations using energy spectra in frequency space and structure functions in real space. We compare with complementary Eulerian velocity statistics. We find that simultaneous inverse energy and enstrophy ranges present in spectra are not directly echoed in real-space moments of velocity difference. Nevertheless, the spectral ranges line up well with features of moment ratios, indicating that although the moments are not exhibiting unambiguous scaling, the behavior of the probability distribution functions is changing over short ranges of length scales. Implications for understanding weakly forced 2D turbulence with simultaneous inverse and direct cascades are discussed.

Lagrangian statistics in weakly forced two-dimensional turbulence

DOE Office of Scientific and Technical Information (OSTI.GOV)

Rivera, Michael K.; Ecke, Robert E.

Measurements of Lagrangian single-point and multiple-point statistics in a quasi-two-dimensional stratified layer system are reported. The system consists of a layer of salt water over an immiscible layer of Fluorinert and is forced electromagnetically so that mean-squared vorticity is injected at a well-defined spatial scale r i. Simultaneous cascades develop in which enstrophy flows predominately to small scales whereas energy cascades, on average, to larger scales. Lagrangian correlations and one- and two-point displacements are measured for random initial conditions and for initial positions within topological centers and saddles. Some of the behavior of these quantities can be understood in termsmore » of the trapping characteristics of long-lived centers, the slow motion near strong saddles, and the rapid fluctuations outside of either centers or saddles. We also present statistics of Lagrangian velocity fluctuations using energy spectra in frequency space and structure functions in real space. We compare with complementary Eulerian velocity statistics. We find that simultaneous inverse energy and enstrophy ranges present in spectra are not directly echoed in real-space moments of velocity difference. Nevertheless, the spectral ranges line up well with features of moment ratios, indicating that although the moments are not exhibiting unambiguous scaling, the behavior of the probability distribution functions is changing over short ranges of length scales. Furthermore, implications for understanding weakly forced 2D turbulence with simultaneous inverse and direct cascades are discussed.« less

Lagrangian statistics in weakly forced two-dimensional turbulence

DOE PAGES

Rivera, Michael K.; Ecke, Robert E.

2016-01-14

Measurements of Lagrangian single-point and multiple-point statistics in a quasi-two-dimensional stratified layer system are reported. The system consists of a layer of salt water over an immiscible layer of Fluorinert and is forced electromagnetically so that mean-squared vorticity is injected at a well-defined spatial scale r i. Simultaneous cascades develop in which enstrophy flows predominately to small scales whereas energy cascades, on average, to larger scales. Lagrangian correlations and one- and two-point displacements are measured for random initial conditions and for initial positions within topological centers and saddles. Some of the behavior of these quantities can be understood in termsmore » of the trapping characteristics of long-lived centers, the slow motion near strong saddles, and the rapid fluctuations outside of either centers or saddles. We also present statistics of Lagrangian velocity fluctuations using energy spectra in frequency space and structure functions in real space. We compare with complementary Eulerian velocity statistics. We find that simultaneous inverse energy and enstrophy ranges present in spectra are not directly echoed in real-space moments of velocity difference. Nevertheless, the spectral ranges line up well with features of moment ratios, indicating that although the moments are not exhibiting unambiguous scaling, the behavior of the probability distribution functions is changing over short ranges of length scales. Furthermore, implications for understanding weakly forced 2D turbulence with simultaneous inverse and direct cascades are discussed.« less

A Chiang-type lagrangian in CP^2

NASA Astrophysics Data System (ADS)

Cannas da Silva, Ana

2018-03-01

We analyse a monotone lagrangian in CP^2 that is hamiltonian isotopic to the standard lagrangian RP^2, yet exhibits a distinguishing behaviour under reduction by one of the toric circle actions, namely it intersects transversally the reduction level set and it projects one-to-one onto a great circle in CP^1. This lagrangian thus provides an example of embedded composition fitting work of Wehrheim-Woodward and Weinstein.

Lagrangian description of warm plasmas

NASA Technical Reports Server (NTRS)

Kim, H.

1970-01-01

Efforts are described to extend the averaged Lagrangian method of describing small signal wave propagation and nonlinear wave interaction, developed by earlier workers for cold plasmas, to the more general conditions of warm collisionless plasmas, and to demonstrate particularly the effectiveness of the method in analyzing wave-wave interactions. The theory is developed for both the microscopic description and the hydrodynamic approximation to plasma behavior. First, a microscopic Lagrangian is formulated rigorously, and expanded in terms of perturbations about equilibrium. Two methods are then described for deriving a hydrodynamic Lagrangian. In the first of these, the Lagrangian is obtained by velocity integration of the exact microscopic Lagrangian. In the second, the expanded hydrodynamic Lagrangian is obtained directly from the expanded microscopic Lagrangian. As applications of the microscopic Lagrangian, the small-signal dispersion relations and the coupled mode equations are derived for all possible waves in a warm infinite, weakly inhomogeneous magnetoplasma, and their interactions are examined.

An Extended Lagrangian Method

NASA Technical Reports Server (NTRS)

Liou, Meng-Sing

1995-01-01

A unique formulation of describing fluid motion is presented. The method, referred to as 'extended Lagrangian method,' is interesting from both theoretical and numerical points of view. The formulation offers accuracy in numerical solution by avoiding numerical diffusion resulting from mixing of fluxes in the Eulerian description. The present method and the Arbitrary Lagrangian-Eulerian (ALE) method have a similarity in spirit-eliminating the cross-streamline numerical diffusion. For this purpose, we suggest a simple grid constraint condition and utilize an accurate discretization procedure. This grid constraint is only applied to the transverse cell face parallel to the local stream velocity, and hence our method for the steady state problems naturally reduces to the streamline-curvature method, without explicitly solving the steady stream-coordinate equations formulated a priori. Unlike the Lagrangian method proposed by Loh and Hui which is valid only for steady supersonic flows, the present method is general and capable of treating subsonic flows and supersonic flows as well as unsteady flows, simply by invoking in the same code an appropriate grid constraint suggested in this paper. The approach is found to be robust and stable. It automatically adapts to flow features without resorting to clustering, thereby maintaining rather uniform grid spacing throughout and large time step. Moreover, the method is shown to resolve multi-dimensional discontinuities with a high level of accuracy, similar to that found in one-dimensional problems.

Lagrangian descriptors in dissipative systems.

PubMed

Junginger, Andrej; Hernandez, Rigoberto

2016-11-09

The reaction dynamics of time-dependent systems can be resolved through a recrossing-free dividing surface associated with the transition state trajectory-that is, the unique trajectory which is bound to the barrier region for all time in response to a given time-dependent potential. A general procedure based on the minimization of Lagrangian descriptors has recently been developed by Craven and Hernandez [Phys. Rev. Lett., 2015, 115, 148301] to construct this particular trajectory without requiring perturbative expansions relative to the naive transition state point at the top of the barrier. The extension of the method to account for dissipation in the equations of motion requires additional considerations established in this paper because the calculation of the Lagrangian descriptor involves the integration of trajectories in forward and backward time. The two contributions are in general very different because the friction term can act as a source (in backward time) or sink (in forward time) of energy, leading to the possibility that information about the phase space structure may be lost due to the dominance of only one of the terms. To compensate for this effect, we introduce a weighting scheme within the Lagrangian descriptor and demonstrate that for thermal Langevin dynamics it preserves the essential phase space structures, while they are lost in the nonweighted case.

[MPLW515L point mutation in patients with myeloproliferative disease].

PubMed

Xia, Jun; Xu, Wei; Zhang, Su-Jiang; Fan, Lei; Qiao, Chun; Li, Jian-Yong

2008-12-01

In order to investigate the frequency of MPLW515L and JAK2V617F point mutations of the patients with myeloproliferative disease (MPD) in Nanjing area, MPLW515L and JAK2V617F point mutations were simultaneously detected by alleles specific polymerase chain reaction (AS-PCR) and sequencing in 190 MPD patients. The results showed that MPLW515L point mutation was detected in 1 out of 102 essential thrombocythemia (ET) patients (1.0%) and was not detected in 32 polycythemia vera (PV) patients, 13 idiopathic myelofibrosis (IMF) patients, 43 chronic myelogenous leukemia (CML) patients. JAK2V617F point mutation was detected in 20 out of 32 PV patients (62.5%), 43 out of 102 ET patients (42.2%), 5 out of 13 IMF patients (38.5%), and was not detected in 43 CML patients. It is concluded that MPLW515L point mutation exists in ET patient, but is not found in PV, IMF and CML. JAK2V617F point mutation exists in PV, ET and IMF, but not in CML.

Construction of Lagrangians and Hamiltonians from the Equation of Motion

ERIC Educational Resources Information Center

Yan, C. C.

1978-01-01

Demonstrates that infinitely many Lagrangians and Hamiltonians can be constructed from a given equation of motion. Points out the lack of an established criterion for making a proper selection. (Author/GA)

About non standard Lagrangians in cosmology

DOE Office of Scientific and Technical Information (OSTI.GOV)

Dimitrijevic, Dragoljub D.; Milosevic, Milan

A review of non standard Lagrangians present in modern cosmological models will be considered. Well known example of non standard Lagrangian is Dirac-Born-Infeld (DBI) type Lagrangian for tachyon field. Another type of non standard Lagrangian under consideration contains scalar field which describes open p-adic string tachyon and is called p-adic string theory Lagrangian. We will investigate homogenous cases of both DBI and p-adic fields and obtain Lagrangians of the standard type which have the same equations of motions as aforementioned non standard one.

Cooperative Convex Optimization in Networked Systems: Augmented Lagrangian Algorithms With Directed Gossip Communication

NASA Astrophysics Data System (ADS)

Jakovetic, Dusan; Xavier, João; Moura, José M. F.

2011-08-01

We study distributed optimization in networked systems, where nodes cooperate to find the optimal quantity of common interest, x=x^\\star. The objective function of the corresponding optimization problem is the sum of private (known only by a node,) convex, nodes' objectives and each node imposes a private convex constraint on the allowed values of x. We solve this problem for generic connected network topologies with asymmetric random link failures with a novel distributed, decentralized algorithm. We refer to this algorithm as AL-G (augmented Lagrangian gossiping,) and to its variants as AL-MG (augmented Lagrangian multi neighbor gossiping) and AL-BG (augmented Lagrangian broadcast gossiping.) The AL-G algorithm is based on the augmented Lagrangian dual function. Dual variables are updated by the standard method of multipliers, at a slow time scale. To update the primal variables, we propose a novel, Gauss-Seidel type, randomized algorithm, at a fast time scale. AL-G uses unidirectional gossip communication, only between immediate neighbors in the network and is resilient to random link failures. For networks with reliable communication (i.e., no failures,) the simplified, AL-BG (augmented Lagrangian broadcast gossiping) algorithm reduces communication, computation and data storage cost. We prove convergence for all proposed algorithms and demonstrate by simulations the effectiveness on two applications: l_1-regularized logistic regression for classification and cooperative spectrum sensing for cognitive radio networks.

78 FR 21491 - DeltaPoint Capital IV, L.P., DeltaPoint Capital IV (New York), L.P.; Notice Seeking Exemption...

Federal Register 2010, 2011, 2012, 2013, 2014

2013-04-10

... SMALL BUSINESS ADMINISTRATION [License No. 02/02-0662, 02/02-0661] DeltaPoint Capital IV, L.P., DeltaPoint Capital IV (New York), L.P.; Notice Seeking Exemption Under Section 312 of the Small Business Investment Act, Conflicts of Interest Notice is hereby given that DeltaPoint Capital IV, L.P. and DeltaPoint...

Accuracy Progressive Calculation of Lagrangian Trajectories from Gridded Velocity Field

DTIC Science & Technology

2014-01-01

traced (Vries and Doos 2001). The two types of velocity are convertible. Routine ocean data 36 assimilation systems (Galanis et al. 2006; Lozano et...1b) 54 The position of each fluid particle, R(t) = [x(t), y(t), z(t)], is specified in the Lagrangian system . 55 The connection...coordinate 233 system at the southwest corner (Fig. 6). The x*- and y*- axes point eastward and northward, 234 respectively. Here, the superscript

Lagrangian Vortices in Developing Tropical Cyclones

DTIC Science & Technology

2015-06-25

45 −6 −4 −2 0 2 4 6 8 10 Hyperbolic trajectory Vortex core Shear sheath Lagrangian...5 10 15 20 25 30 35 40 45 −1 −0.8 −0.6 −0.4 −0.2 0 0.2 0.4 0.6 0.8 1 × 10−9 (a) P13LP16E P17L 700 hPa 700 hPa 4590804102 100 80 60 40 20 0 0 5 10...15 20 25 30 35 40 45 −8 −6 −4 −2 0 2 4 6 8 10(b) Figure 4. The (a) OW (s−2) and (b) OWLag (radians) fields for the Atlantic development region at

CELFE: Coupled Eulerian-Lagrangian Finite Element program for high velocity impact. Part 1: Theory and formulation. [hydroelasto-viscoplastic model

NASA Technical Reports Server (NTRS)

Lee, C. H.

1978-01-01

A 3-D finite element program capable of simulating the dynamic behavior in the vicinity of the impact point, together with predicting the dynamic response in the remaining part of the structural component subjected to high velocity impact is discussed. The finite algorithm is formulated in a general moving coordinate system. In the vicinity of the impact point contained by a moving failure front, the relative velocity of the coordinate system will approach the material particle velocity. The dynamic behavior inside the region is described by Eulerian formulation based on a hydroelasto-viscoplastic model. The failure front which can be regarded as the boundary of the impact zone is described by a transition layer. The layer changes the representation from the Eulerian mode to the Lagrangian mode outside the failure front by varying the relative velocity of the coordinate system to zero. The dynamic response in the remaining part of the structure described by the Lagrangian formulation is treated using advanced structural analysis. An interfacing algorithm for coupling CELFE with NASTRAN is constructed to provide computational capabilities for large structures.

A few words about resonances in the electroweak effective Lagrangian

DOE Office of Scientific and Technical Information (OSTI.GOV)

Rosell, Ignasi; Pich, Antonio; Santos, Joaquín

Contrary to a widely spread believe, we have demonstrated that strongly coupled electroweak models including both a light Higgs-like boson and massive spin-1 resonances are not in conflict with experimental constraints on the oblique S and T parameters. We use an effective Lagrangian implementing the chiral symmetry breaking SU (2){sub L} ⊗ SU (2){sub R} → SU (2){sub L+R} that contains the Standard Model gauge bosons coupled to the electroweak Goldstones, one Higgs-like scalar state h with mass m{sub h} = 126 GeV and the lightest vector and axial-vector resonance multiplets V and A. We have considered the one-loop calculationmore » of S and T in order to study the viability of these strongly-coupled scenarios, being short-distance constraints and dispersive relations the main ingredients of the calculation. Once we have constrained the resonance parameters, we do a first approach to the determination of the low energy constants of the electroweak effective theory at low energies (without resonances). We show this determination in the case of the purely Higgsless bosonic Lagrangian.« less

Implications of Lagrangian transport for coupled chemistry-climate simulations

NASA Astrophysics Data System (ADS)

Stenke, A.; Dameris, M.; Grewe, V.; Garny, H.

2008-10-01

For the first time a purely Lagrangian transport algorithm is applied in a fully coupled chemistry-climate model (CCM). We use the Lagrangian scheme ATTILA for the transport of water vapour, cloud water and chemical trace species in the ECHAM4.L39(DLR)/CHEM (E39C) CCM. The advantage of the Lagrangian approach is that it is numerically non-diffusive and therefore maintains steeper and more realistic gradients than the operational semi-Lagrangian transport scheme. In case of radiatively active species changes in the simulated distributions feed back to model dynamics which in turn affect the modelled transport. The implications of the Lagrangian transport scheme for stratospheric model dynamics and tracer distributions in the upgraded model version E39C-ATTILA (E39C-A) are evaluated by comparison with observations and results of the E39C model with the operational semi-Lagrangian advection scheme. We find that several deficiencies in stratospheric dynamics in E39C seem to originate from a pronounced modelled wet bias and an associated cold bias in the extra-tropical lowermost stratosphere. The reduction of the simulated moisture and temperature bias in E39C-A leads to a significant advancement of stratospheric dynamics in terms of the mean state as well as annual and interannual variability. As a consequence of the favourable numerical characteristics of the Lagrangian transport scheme and the improved model dynamics, E39C-A generally shows more realistic stratospheric tracer distributions: Compared to E39C high stratospheric chlorine (Cly) concentrations extend further downward and agree now well with analyses derived from observations. Therefore E39C-A realistically covers the altitude of maximum ozone depletion in the stratosphere. The location of the ozonopause, i.e. the transition from low tropospheric to high stratospheric ozone values, is also clearly improved in E39C-A. Furthermore, the simulated temporal evolution of stratospheric Cly in the past is

Form of the manifestly covariant Lagrangian

NASA Astrophysics Data System (ADS)

Johns, Oliver Davis

1985-10-01

The preferred form for the manifestly covariant Lagrangian function of a single, charged particle in a given electromagnetic field is the subject of some disagreement in the textbooks. Some authors use a ``homogeneous'' Lagrangian and others use a ``modified'' form in which the covariant Hamiltonian function is made to be nonzero. We argue in favor of the ``homogeneous'' form. We show that the covariant Lagrangian theories can be understood only if one is careful to distinguish quantities evaluated on the varied (in the sense of the calculus of variations) world lines from quantities evaluated on the unvaried world lines. By making this distinction, we are able to derive the Hamilton-Jacobi and Klein-Gordon equations from the ``homogeneous'' Lagrangian, even though the covariant Hamiltonian function is identically zero on all world lines. The derivation of the Klein-Gordon equation in particular gives Lagrangian theoretical support to the derivations found in standard quantum texts, and is also shown to be consistent with the Feynman path-integral method. We conclude that the ``homogeneous'' Lagrangian is a completely adequate basis for covariant Lagrangian theory both in classical and quantum mechanics. The article also explores the analogy with the Fermat theorem of optics, and illustrates a simple invariant notation for the Lagrangian and other four-vector equations.

Lagrangian formulation of the 2(2S+1)-component model and its connection with the skew symmetric/tensor description

DOE Office of Scientific and Technical Information (OSTI.GOV)

Dvoeglazov, V.V.

1993-12-01

In the framework of the 2(2S + 1)-component theory for massless particles, the dynamical invariants have been derived from the Lagrangian density which is considered to be a 4-vector. A la Majorana interpretation of the 6-component spinors, the field operators of S=1 particles, as the left- and right-circularly polarized radiation, leads the author to the conserved quantities which are analogous to ones obtained by Lipkin and Sudbery. The scalar Lagrangian of this model is shown to be equivalent to the Lagrangian of a free massless field, introduced by Hayashi. As a consequence of a new {open_quotes}gauge{close_quotes} invariance this skew-symmetric fieldmore » describes physical particles with the longitudinal components only.« less

Long Term Missions at the Sun-Earth Libration Point L1: ACE, SOHO, and WIND

NASA Technical Reports Server (NTRS)

Roberts, Craig E.

2011-01-01

Three heliophysics missions -- the Advanced Composition Explorer (ACE), Solar Heliospheric Observatory (SOHO), and the Global Geoscience WIND -- have been orbiting the Sun-Earth interior libration point L1 continuously since 1997, 1996, and 2004, respectively. ACE and WIND (both NASA missions) and SOHO (an ESA-NASA joint mission) are all operated from the NASA Goddard Space Flight Center (GSFC). While ACE and SOHO have been dedicated libration point orbiters since their launches, WIND has had also a remarkable 10-year career flying a deep-space, multiple lunar-flyby trajectory prior to 2004. That era featured 36 targeted lunar flybys with excursions to both L1 and L2 before its final insertion in L1 orbit. A figure depicts the orbits of the three spacecraft, showing projections of the orbits onto the orthographic planes of a solar rotating ecliptic frame of reference. The SOHO orbit is a quasi-periodic halo orbit, where the frequencies of the in-plane and out-of-plane motions are practically equal. Such an orbit is seen to repeat itself with a period of approximately 178 days. For ACE and WIND, the frequencies of the in-plane and out-of-plane motions are unequal, giving rise to the characteristic Lissajous motion. ACE's orbit is of moderately small amplitude, whereas WIND's orbit is a large-amplitude Lissajous of dimensions close to those of the SOHO halo orbit. As motion about the collinear points is inherently unstable, stationkeeping maneuvers are necessary to prevent orbital decay and eventual escape from the L1 region. Though the three spacecraft are dissimilar (SOHO is a 3-axis stabilized Sun pointer, WIND is a spin-stabilized ecliptic pole pointer, and ACE is also spin-stabilized with its spin axis maintained between 4 and 20 degrees of the Sun), the stationkeeping technique for the three is fundamentally the same. The technique consists of correcting the energy of the orbit via a delta-V directed parallel or anti-parallel to the Spacecraft-to-Sun line. SOHO

Reduction of numerical diffusion in three-dimensional vortical flows using a coupled Eulerian/Lagrangian solution procedure

NASA Technical Reports Server (NTRS)

Felici, Helene M.; Drela, Mark

1993-01-01

A new approach based on the coupling of an Eulerian and a Lagrangian solver, aimed at reducing the numerical diffusion errors of standard Eulerian time-marching finite-volume solvers, is presented. The approach is applied to the computation of the secondary flow in two bent pipes and the flow around a 3D wing. Using convective point markers the Lagrangian approach provides a correction of the basic Eulerian solution. The Eulerian flow in turn integrates in time the Lagrangian state-vector. A comparison of coarse and fine grid Eulerian solutions makes it possible to identify numerical diffusion. It is shown that the Eulerian/Lagrangian approach is an effective method for reducing numerical diffusion errors.

Extreme Lagrangian acceleration in confined turbulent flow.

PubMed

Kadoch, Benjamin; Bos, Wouter J T; Schneider, Kai

2008-05-09

A Lagrangian study of two-dimensional turbulence for two different geometries, a periodic and a confined circular geometry, is presented to investigate the influence of solid boundaries on the Lagrangian dynamics. It is found that the Lagrangian acceleration is even more intermittent in the confined domain than in the periodic domain. The flatness of the Lagrangian acceleration as a function of the radius shows that the influence of the wall on the Lagrangian dynamics becomes negligible in the center of the domain, and it also reveals that the wall is responsible for the increased intermittency. The transition in the Lagrangian statistics between this region, not directly influenced by the walls, and a critical radius which defines a Lagrangian boundary layer is shown to be very sharp with a sudden increase of the acceleration flatness from about 5 to about 20.

Orbit Determination (OD) Error Analysis Results for the Triana Sun-Earth L1 Libration Point Mission and for the Fourier Kelvin Stellar Interferometer (FKSI) Sun-Earth L2 Libration Point Mission Concept

NASA Technical Reports Server (NTRS)

Marr, Greg C.

2003-01-01

The Triana spacecraft was designed to be launched by the Space Shuttle. The nominal Triana mission orbit will be a Sun-Earth L1 libration point orbit. Using the NASA Goddard Space Flight Center's Orbit Determination Error Analysis System (ODEAS), orbit determination (OD) error analysis results are presented for all phases of the Triana mission from the first correction maneuver through approximately launch plus 6 months. Results are also presented for the science data collection phase of the Fourier Kelvin Stellar Interferometer Sun-Earth L2 libration point mission concept with momentum unloading thrust perturbations during the tracking arc. The Triana analysis includes extensive analysis of an initial short arc orbit determination solution and results using both Deep Space Network (DSN) and commercial Universal Space Network (USN) statistics. These results could be utilized in support of future Sun-Earth libration point missions.

Classical Lagrangians and Finsler structures for the nonminimal fermion sector of the standard model extension

NASA Astrophysics Data System (ADS)

Schreck, M.

2016-05-01

This article is devoted to finding classical point-particle equivalents for the fermion sector of the nonminimal standard model extension (SME). For a series of nonminimal operators, such Lagrangians are derived at first order in Lorentz violation using the algebraic concept of Gröbner bases. Subsequently, the Lagrangians serve as a basis for reanalyzing the results of certain kinematic tests of special relativity that were carried out in the past century. Thereby, a number of new constraints on coefficients of the nonminimal SME is obtained. In the last part of the paper we point out connections to Finsler geometry.

Asymptotic-preserving Lagrangian approach for modeling anisotropic transport in magnetized plasmas for arbitrary magnetic fields

NASA Astrophysics Data System (ADS)

Chacon, Luis; Del-Castillo-Negrete, Diego; Hauck, Cory

2012-10-01

Modeling electron transport in magnetized plasmas is extremely challenging due to the extreme anisotropy between parallel (to the magnetic field) and perpendicular directions (χ/χ˜10^10 in fusion plasmas). Recently, a Lagrangian Green's function approach, developed for the purely parallel transport case,footnotetextD. del-Castillo-Negrete, L. Chac'on, PRL, 106, 195004 (2011)^,footnotetextD. del-Castillo-Negrete, L. Chac'on, Phys. Plasmas, 19, 056112 (2012) has been extended to the anisotropic transport case in the tokamak-ordering limit with constant density.footnotetextL. Chac'on, D. del-Castillo-Negrete, C. Hauck, JCP, submitted (2012) An operator-split algorithm is proposed that allows one to treat Eulerian and Lagrangian components separately. The approach is shown to feature bounded numerical errors for arbitrary χ/χ ratios, which renders it asymptotic-preserving. In this poster, we will present the generalization of the Lagrangian approach to arbitrary magnetic fields. We will demonstrate the potential of the approach with various challenging configurations, including the case of transport across a magnetic island in cylindrical geometry.

A Eulerian-Lagrangian Model to Simulate Two-Phase/Particulate Flows

NASA Technical Reports Server (NTRS)

Apte, S. V.; Mahesh, K.; Lundgren, T.

2003-01-01

Figure 1 shows a snapshot of liquid fuel spray coming out of an injector nozzle in a realistic gas-turbine combustor. Here the spray atomization was simulated using a stochastic secondary breakup model (Apte et al. 2003a) with point-particle approximation for the droplets. Very close to the injector, it is observed that the spray density is large and the droplets cannot be treated as point-particles. The volume displaced by the liquid in this region is significant and can alter the gas-phase ow and spray evolution. In order to address this issue, one can compute the dense spray regime by an Eulerian-Lagrangian technique using advanced interface tracking/level-set methods (Sussman et al. 1994; Tryggvason et al. 2001; Herrmann 2003). This, however, is computationally intensive and may not be viable in realistic complex configurations. We therefore plan to develop a methodology based on Eulerian-Lagrangian technique which will allow us to capture the essential features of primary atomization using models to capture interactions between the fluid and droplets and which can be directly applied to the standard atomization models used in practice. The numerical scheme for unstructured grids developed by Mahesh et al. (2003) for incompressible flows is modified to take into account the droplet volume fraction. The numerical framework is directly applicable to realistic combustor geometries. Our main objectives in this work are: Develop a numerical formulation based on Eulerian-Lagrangian techniques with models for interaction terms between the fluid and particles to capture the Kelvin- Helmholtz type instabilities observed during primary atomization. Validate this technique for various two-phase and particulate flows. Assess its applicability to capture primary atomization of liquid jets in conjunction with secondary atomization models.

Lagrangian transport in a class of three-dimensional buoyancy-driven flows

NASA Astrophysics Data System (ADS)

Contreras, Sebastian; Speetjens, Michel; Clercx, Herman

2017-11-01

The study concerns the Lagrangian dynamics of three-dimensional (3D) buoyancy-driven cavity flows under steady and laminar conditions due to a global temperature gradient imposed via an opposite hot and cold sidewall. This serves as archetypal configuration for natural-convection flows in which gravity is perpendicular to the global temperature gradient. Limited insight into the Lagrangian properties of this class of flows motivates this study. The 3D Lagrangian dynamics are investigated in terms of the generic structure of the Lagrangian flow topology that is described in terms of the Grashof number (Gr) and the Prandtl number (Pr). Gr is the principal control parameter for the flow topology: vanishing Gr yields a state of closed streamlines (integrable state); increasing Gr causes the formation of toroidal coherent structures embedded in chaotic streamlines governed by Hamiltonian mechanisms. Fluid inertia prevails for ``smaller'' Gr. A buoyancy-induced bifurcation of the flow topology occurs for ``larger'' Gr and underlies the emergence of ``secondary rolls'' and secondary tori for ``larger'' Pr. Stagnation points and corresponding manifold interactions are key to the dynamics. S.C. acknowledges financial support from Consejo Nacional de Ciencia y Tecnología (CONACYT).

Dynamics of a satellite and normalization around lagrangian points in the Neptune-Triton system

NASA Astrophysics Data System (ADS)

Yokoyama, T.; Stuchi, T. J.; Solórzano, C. R. H.; Corria, A. A.; Prado, A. F. B.; Winter, O. C.; Winter, S. M. G.; Sanchez, D. M.

Since a mission to Pluto-Charon is in progress a similar mission to Neptune-Triton system probably is just a matter of time Besides the interesting items listed in the Hammel s proposal to explore Neptune Hammel et al 2002 there are others very interesting points Triton is a very large satellite with M T M N approx 2 09 times 10 -4 where M T and M N are the masses of Triton and Neptune respectively Its current inclination with respect to Neptune s equator is about 157 345 0 This unusual high and retrograde inclination for a very large inner satellite makes this problem unique in our solar system Hammel et al 2002 also propose that a top priority in this mission should be a Neptune-Triton orbiter not just a flyby tour Therefore assuming a massless spacecraft orbiting Neptune-Triton system basically the system can be stated in terms of the classical restricted three body problem The new ingredient is the Neptune oblateness and the retrograde motion of Triton With some slight displacements the lagrangian equilibria points still exist as well as many of the properties of the classical problem In this work we first give an extensive numerical exploration in the case when the spacecraft orbits Triton considering Sun Neptune and its oblateness as disturbers In the plane bf a times bf I where bf a is the semi major axis and bf I is the inclination of the orbiter we give a plot of the stable regions where the massless can survive for thousand of years Retrograde and direct inclinations were

Normal Modes of a Lagrangian System Constrained in a Potential Well.

DTIC Science & Technology

1983-12-01

A’ -137 948 NORMAL MODES OF A LFHbRANGIAN SYSTEM CONSTRAINED INvi P0TENTIAL WELL(U WISCONSNN UNIV-MADISON MATHEMATICS RESEARCH CENTER V EN DEC F1...Carolina 27709 DT FLE OP Y UNIVERSITY OF WISCONSIN-MADISON MATHEMATICS RESEARCH CENTER NORMAL MODES OF A LAGRANGIAN SYSTEM CONSTRAINED IN A POTENTIAL WELL...respect to the norm lYE [f i + 2 yi )dtl/ 0 Since H I(S’ 1 n’) C CO(S, fle ), then the set A 1 0 is an open set in H1 (lf’) The periodic solution of

Lagrangian averaging with geodesic mean

NASA Astrophysics Data System (ADS)

Oliver, Marcel

2017-11-01

This paper revisits the derivation of the Lagrangian averaged Euler (LAE), or Euler-α equations in the light of an intrinsic definition of the averaged flow map as the geodesic mean on the volume-preserving diffeomorphism group. Under the additional assumption that first-order fluctuations are statistically isotropic and transported by the mean flow as a vector field, averaging of the kinetic energy Lagrangian of an ideal fluid yields the LAE Lagrangian. The derivation presented here assumes a Euclidean spatial domain without boundaries.

Lagrangian averaging with geodesic mean.

PubMed

Oliver, Marcel

2017-11-01

This paper revisits the derivation of the Lagrangian averaged Euler (LAE), or Euler- α equations in the light of an intrinsic definition of the averaged flow map as the geodesic mean on the volume-preserving diffeomorphism group. Under the additional assumption that first-order fluctuations are statistically isotropic and transported by the mean flow as a vector field, averaging of the kinetic energy Lagrangian of an ideal fluid yields the LAE Lagrangian. The derivation presented here assumes a Euclidean spatial domain without boundaries.

Lagrangian Perturbation Approach to the Formation of Large-scale Structure

NASA Astrophysics Data System (ADS)

Buchert, Thomas

The present lecture notes address three columns on which the Lagrangian perturbation approach to cosmological dynamics is based: 1. the formulation of a Lagrangian theory of self-gravitating flows in which the dynamics is described in terms of a single field variable; 2. the procedure, how to obtain the dynamics of Eulerian fields from the Lagrangian picture, and 3. a precise definition of a Newtonian cosmology framework in which Lagrangian perturbation solutions can be studied. While the first is a discussion of the basic equations obtained by transforming the Eulerian evolution and field equations to the Lagrangian picture, the second exemplifies how the Lagrangian theory determines the evolution of Eulerian fields including kinematical variables like expansion, vorticity, as well as the shear and tidal tensors. The third column is based on a specification of initial and boundary conditions, and in particular on the identification of the average flow of an inhomogeneous cosmology with a `Hubble-flow'. Here, we also look at the limits of the Lagrangian perturbation approach as inferred from comparisons with N-body simulations and illustrate some striking properties of the solutions.

Determination of the Minimum Delta V Transfer Trajectory from a Low Earth Orbit to a Stable Orbit Around the Lagrangian Point L4 in a Restricted Four-Body System.

DTIC Science & Technology

1979-12-01

the Moon revolves around the Earth. T. A . Heppenheimer (Ref 5) and his colleagues B. O’Leary and D. Kaplan (Ref 7) have also examined the pro- blem of...colony location and transfer trajectories. Heppenheimer found a 2/1 resonant orbit around the Earth 2S that could be reached from L2 by a Hohmann...would not allow a catcher near either of these points. In their work together, O’Leary, Kaplan and Heppenheimer found a transfer trajectory from L2

Three dimensional Lagrangian structures in the Antarctic Polar Vortex.

NASA Astrophysics Data System (ADS)

Mancho, Ana M.; Garcia-Garrido, Victor J.; Curbelo, Jezabel; Niang, Coumba; Mechoso, Carlos R.; Wiggins, Stephen

2017-04-01

Dynamical systems theory has supported the description of transport processes in fluid dynamics. For understanding trajectory patterns in chaotic advection the geometrical approach by Poincaré seeks for spatial structures that separate regions corresponding to qualitatively different types of trajectories. These structures have been referred to as Lagrangian Coherent Structures (LCS), which typically in geophysical flows are well described under the approach of incompressible 2D flows. Different tools have been used to visualize LCS. In this presentation we use Lagrangian Descriptors [1,2,3,4] (function M) for visualizing 3D Lagrangian structures in the atmosphere, in particular in the Antarctic Polar Vortex. The function M is computed in a fully 3D incompressible flow obtained from data provided by the European Centre for Medium-Range Weather Forecast and it is represented in 2D surfaces. We discuss the findings during the final warming that took place in the spring of 1979 [5]. This research is supported by MINECO grant MTM2014-56392-R. Support is acknowledged also from CSIC grant COOPB20265, U.S. NSF grant AGS-1245069 and ONR grant No. N00014- 01-1-0769. C. Niang acknowledges Fundacion Mujeres por Africa and ICMAT Severo Ochoa project SEV-2011-0087 for financial support. [1] C. Mendoza, A. M. Mancho. The hidden geometry of ocean flows. Physical Review Letters 105 (2010), 3, 038501-1-038501-4. [2] A. M. Mancho, S. Wiggins, J. Curbelo, C. Mendoza. Lagrangian Descriptors: A Method for Revealing Phase Space Structures of General Time Dependent Dynamical Systems. Communications in Nonlinear Science and Numerical Simulation. 18 (2013) 3530-3557. [3] C. Lopesino, F. Balibrea-Iniesta, S. Wiggins and A. M. Mancho. Lagrangian descriptors for two dimensional, area preserving autonomous and nonautonomous maps. Communications in Nonlinear Science and Numerical Simulations, 27 (2015) (1-3), 40-51. [4] C. Lopesino, F. Balibrea-Iniesta, V. J. García-Garrido, S. Wiggins, and A

Lagrangians for generalized Argyres-Douglas theories

NASA Astrophysics Data System (ADS)

Benvenuti, Sergio; Giacomelli, Simone

2017-10-01

We continue the study of Lagrangian descriptions of N=2 Argyres-Douglas theories. We use our recent interpretation in terms of sequential confinement to guess the Lagrangians of all the Argyres-Douglas models with Abelian three dimensional mirror. We find classes of four dimensional N=1 quivers that flow in the infrared to generalized Argyres-Douglas theories, such as the ( A k , A kN + N -1) models. We study in detail how the N=1 chiral rings map to the Coulomb and Higgs Branches of the N=2 CFT's. The three dimensional mirror RG flows are shown to land on the N=4 complete graph quivers. We also compactify to three dimensions the gauge theory dual to ( A 1, D 4), and find the expected Abelianization duality with N=4 SQED with 3 flavors.

Lagrangian study of transport of subarctic water across the Subpolar Front in the Japan Sea

NASA Astrophysics Data System (ADS)

Prants, Sergey V.; Uleysky, Michael Yu.; Budyansky, Maxim V.

2018-06-01

The southward near-surface transport of transformed subarctic water across the Subpolar Front in the Japan Sea is simulated and analyzed based on altimeter data from January 1, 1993 to December 31, 2017. Computing Lagrangian indicators for a large number of synthetic particles, advected by the AVISO velocity field, we find preferred transport pathways across the Subpolar Front. The southward transport occurs mainly in the central part of the frontal zone due to suitable dispositions of mesoscale eddies promoting propagation of subarctic water to the south. It is documented with the help of Lagrangian origin and L-maps and verified by the tracks of available drifters. The transport of transformed subarctic water to the south is compared with the transport of transformed subtropical water to the north simulated by Prants et al. (Nonlinear Process Geophys 24(1):89-99, 2017c).

Lagrangian study of transport of subarctic water across the Subpolar Front in the Japan Sea

NASA Astrophysics Data System (ADS)

Prants, Sergey V.; Uleysky, Michael Yu.; Budyansky, Maxim V.

2018-05-01

The southward near-surface transport of transformed subarctic water across the Subpolar Front in the Japan Sea is simulated and analyzed based on altimeter data from January 1, 1993 to December 31, 2017. Computing Lagrangian indicators for a large number of synthetic particles, advected by the AVISO velocity field, we find preferred transport pathways across the Subpolar Front. The southward transport occurs mainly in the central part of the frontal zone due to suitable dispositions of mesoscale eddies promoting propagation of subarctic water to the south. It is documented with the help of Lagrangian origin and L-maps and verified by the tracks of available drifters. The transport of transformed subarctic water to the south is compared with the transport of transformed subtropical water to the north simulated by Prants et al. (Nonlinear Process Geophys 24(1):89-99, 2017c).

Vorticity and symplecticity in multi-symplectic, Lagrangian gas dynamics

NASA Astrophysics Data System (ADS)

Webb, G. M.; Anco, S. C.

2016-02-01

The Lagrangian, multi-dimensional, ideal, compressible gas dynamic equations are written in a multi-symplectic form, in which the Lagrangian fluid labels, m i (the Lagrangian mass coordinates) and time t are the independent variables, and in which the Eulerian position of the fluid element {x}={x}({m},t) and the entropy S=S({m},t) are the dependent variables. Constraints in the variational principle are incorporated by means of Lagrange multipliers. The constraints are: the entropy advection equation S t = 0, the Lagrangian map equation {{x}}t={u} where {u} is the fluid velocity, and the mass continuity equation which has the form J=τ where J={det}({x}{ij}) is the Jacobian of the Lagrangian map in which {x}{ij}=\\partial {x}i/\\partial {m}j and τ =1/ρ is the specific volume of the gas. The internal energy per unit volume of the gas \\varepsilon =\\varepsilon (ρ ,S) corresponds to a non-barotropic gas. The Lagrangian is used to define multi-momenta, and to develop de Donder-Weyl Hamiltonian equations. The de Donder-Weyl equations are cast in a multi-symplectic form. The pullback conservation laws and the symplecticity conservation laws are obtained. One class of symplecticity conservation laws give rise to vorticity and potential vorticity type conservation laws, and another class of symplecticity laws are related to derivatives of the Lagrangian energy conservation law with respect to the Lagrangian mass coordinates m i . We show that the vorticity-symplecticity laws can be derived by a Lie dragging method, and also by using Noether’s second theorem and a fluid relabelling symmetry which is a divergence symmetry of the action. We obtain the Cartan-Poincaré form describing the equations and we discuss a set of differential forms representing the equation system.

Option volatility and the acceleration Lagrangian

NASA Astrophysics Data System (ADS)

Baaquie, Belal E.; Cao, Yang

2014-01-01

This paper develops a volatility formula for option on an asset from an acceleration Lagrangian model and the formula is calibrated with market data. The Black-Scholes model is a simpler case that has a velocity dependent Lagrangian. The acceleration Lagrangian is defined, and the classical solution of the system in Euclidean time is solved by choosing proper boundary conditions. The conditional probability distribution of final position given the initial position is obtained from the transition amplitude. The volatility is the standard deviation of the conditional probability distribution. Using the conditional probability and the path integral method, the martingale condition is applied, and one of the parameters in the Lagrangian is fixed. The call option price is obtained using the conditional probability and the path integral method.

78 FR 21491 - DeltaPoint Capital IV, L.P., DeltaPoint Capital IV (New York), L.P.; Notice Seeking Exemption...

Federal Register 2010, 2011, 2012, 2013, 2014

2013-04-10

... Small Business Investment Act of 1958, as amended (``the Act''), in connection with the financing of a... SMALL BUSINESS ADMINISTRATION [License No. 02/02-0662, 02/02-0661] DeltaPoint Capital IV, L.P., DeltaPoint Capital IV (New York), L.P.; Notice Seeking Exemption Under Section 312 of the Small Business...

Stochastic modeling of Lagrangian accelerations

NASA Astrophysics Data System (ADS)

Reynolds, Andy

2002-11-01

It is shown how Sawford's second-order Lagrangian stochastic model (Phys. Fluids A 3, 1577-1586, 1991) for fluid-particle accelerations can be combined with a model for the evolution of the dissipation rate (Pope and Chen, Phys. Fluids A 2, 1437-1449, 1990) to produce a Lagrangian stochastic model that is consistent with both the measured distribution of Lagrangian accelerations (La Porta et al., Nature 409, 1017-1019, 2001) and Kolmogorov's similarity theory. The later condition is found not to be satisfied when a constant dissipation rate is employed and consistency with prescribed acceleration statistics is enforced through fulfilment of a well-mixed condition.

Lagrangian transport near perturbed periodic lines in three-dimensional unsteady flows

NASA Astrophysics Data System (ADS)

Speetjens, Michel

2015-11-01

Periodic lines formed by continuous strings of periodic points are key organizing entities in the Lagrangian flow topology of certain three-dimensional (3D) time-periodic flows. Such lines generically consist of elliptic and/or hyperbolic points and thus give rise to 3D flow topologies made up of families of concentric closed trajectories embedded in chaotic regions. Weak perturbation destroys the periodic lines and causes said trajectories to coalesce into families of concentric tubes. However, emergence of isolated periodic points near the disintegrating periodic lines and/or partitioning of the original lines into elliptic and hyperbolic segments interrupt the tube formation. This yields incomplete tubes that interact with the (chaotic) environment through their open ends, resulting in intricate and essentially 3D flow topologies These phenomena have been observed in various realistic flows yet the underlying mechanisms are to date only partially understood. This study deepens insight into the (perturbed) Lagrangian dynamics of these flows by way of a linearized representation of the equations of motion near the periodic lines. Predictions on the basis of this investigation are in full (qualitative) agreement with observed behavior in the actual flows

Transfer to the Collinear Libration Point L3 in the Sun-Earth+Moon System

NASA Technical Reports Server (NTRS)

Hou, Xi-yun; Tang, Jing-shi; Liu, Lin

2007-01-01

The collinear libration point L3 of the sun-earth+moon system is an ideal place for some space missions. Although there has been a great amount of work concerning the applications of the other two collinear libration points L1 and L2, little work has been done about the point L3. In this paper, the dynamics of the libration points was briefly introduced first. Then a way to transfer the spacecraft to the collinear libration point L3 via the invariant manifolds of the other two collinear libration points was proposed. Theoretical works under the model of circular restricted three-body problem were done. For the sun-earth+moon system, this model is a good approximation. The results obtained are useful when a transfer trajectory under the real solar system is designed.

Lagrangian particle method for compressible fluid dynamics

NASA Astrophysics Data System (ADS)

Samulyak, Roman; Wang, Xingyu; Chen, Hsin-Chiang

2018-06-01

A new Lagrangian particle method for solving Euler equations for compressible inviscid fluid or gas flows is proposed. Similar to smoothed particle hydrodynamics (SPH), the method represents fluid cells with Lagrangian particles and is suitable for the simulation of complex free surface/multiphase flows. The main contributions of our method, which is different from SPH in all other aspects, are (a) significant improvement of approximation of differential operators based on a polynomial fit via weighted least squares approximation and the convergence of prescribed order, (b) a second-order particle-based algorithm that reduces to the first-order upwind method at local extremal points, providing accuracy and long term stability, and (c) more accurate resolution of entropy discontinuities and states at free interfaces. While the method is consistent and convergent to a prescribed order, the conservation of momentum and energy is not exact and depends on the convergence order. The method is generalizable to coupled hyperbolic-elliptic systems. Numerical verification tests demonstrating the convergence order are presented as well as examples of complex multiphase flows.

Functional integral for non-Lagrangian systems

NASA Astrophysics Data System (ADS)

Kochan, Denis

2010-02-01

A functional integral formulation of quantum mechanics for non-Lagrangian systems is presented. The approach, which we call “stringy quantization,” is based solely on classical equations of motion and is free of any ambiguity arising from Lagrangian and/or Hamiltonian formulation of the theory. The functionality of the proposed method is demonstrated on several examples. Special attention is paid to the stringy quantization of systems with a general A-power friction force -κq˙A. Results for A=1 are compared with those obtained in the approaches by Caldirola-Kanai, Bateman, and Kostin. Relations to the Caldeira-Leggett model and to the Feynman-Vernon approach are discussed as well.

Coherent Lagrangian swirls among submesoscale motions.

PubMed

Beron-Vera, F J; Hadjighasem, A; Xia, Q; Olascoaga, M J; Haller, G

2018-03-05

The emergence of coherent Lagrangian swirls (CLSs) among submesoscale motions in the ocean is illustrated. This is done by applying recent nonlinear dynamics tools for Lagrangian coherence detection on a surface flow realization produced by a data-assimilative submesoscale-permitting ocean general circulation model simulation of the Gulf of Mexico. Both mesoscale and submesoscale CLSs are extracted. These extractions prove the relevance of coherent Lagrangian eddies detected in satellite-altimetry-based geostrophic flow data for the arguably more realistic ageostrophic multiscale flow.

Lagrangian fluid description with simple applications in compressible plasma and gas dynamics

NASA Astrophysics Data System (ADS)

Schamel, Hans

2004-03-01

The Lagrangian fluid description, in which the dynamics of fluids is formulated in terms of trajectories of fluid elements, not only presents an alternative to the more common Eulerian description but has its own merits and advantages. This aspect, which seems to be not fully explored yet, is getting increasing attention in fluid dynamics and related areas as Lagrangian codes and experimental techniques are developed utilizing the Lagrangian point of view with the ultimate goal of a deeper understanding of flow dynamics. In this tutorial review we report on recent progress made in the analysis of compressible, more or less perfect flows such as plasmas and dilute gases. The equations of motion are exploited to get further insight into the formation and evolution of coherent structures, which often exhibit a singular or collapse type behavior occurring in finite time. It is argued that this technique of solution has a broad applicability due to the simplicity and generality of equations used. The focus is on four different topics, the physics of which being governed by simple fluid equations subject to initial and/or boundary conditions. Whenever possible also experimental results are mentioned. In the expansion of a semi-infinite plasma into a vacuum the energetic ion peak propagating supersonically towards the vacuum-as seen in laboratory experiments-is interpreted by means of the Lagrangian fluid description as a relic of a wave breaking scenario of the corresponding inviscid ion dynamics. The inclusion of viscosity is shown numerically to stabilize the associated density collapse giving rise to a well defined fast ion peak reminiscent of adhesive matter. In purely convection driven flows the Lagrangian flow velocity is given by its initial value and hence the Lagrangian velocity gradient tensor can be evaluated accurately to find out the appearance of singularities in density and vorticity and the emergence of new structures such as wavelets in one-dimension (1D

Alternative kinetic energy metrics for Lagrangian systems

NASA Astrophysics Data System (ADS)

Sarlet, W.; Prince, G.

2010-11-01

We examine Lagrangian systems on \\ {R}^n with standard kinetic energy terms for the possibility of additional, alternative Lagrangians with kinetic energy metrics different to the Euclidean one. Using the techniques of the inverse problem in the calculus of variations we find necessary and sufficient conditions for the existence of such Lagrangians. We illustrate the problem in two and three dimensions with quadratic and cubic potentials. As an aside we show that the well-known anomalous Lagrangians for the Coulomb problem can be removed by switching on a magnetic field, providing an appealing resolution of the ambiguous quantizations of the hydrogen atom.

Eulerian and Lagrangian methods for vortex tracking in 2D and 3D flows

NASA Astrophysics Data System (ADS)

Huang, Yangzi; Green, Melissa

2014-11-01

Coherent structures are a key component of unsteady flows in shear layers. Improvement of experimental techniques has led to larger amounts of data and requires of automated procedures for vortex tracking. Many vortex criteria are Eulerian, and identify the structures by an instantaneous local swirling motion in the field, which are indicated by closed or spiral streamlines or pathlines in a reference frame. Alternatively, a Lagrangian Coherent Structures (LCS) analysis is a Lagrangian method based on the quantities calculated along fluid particle trajectories. In the current work, vortex detection is demonstrated on data from the simulation of two cases: a 2D flow with a flat plate undergoing a 45 ° pitch-up maneuver and a 3D wall-bounded turbulence channel flow. Vortices are visualized and tracked by their centers and boundaries using Γ1, the Q criterion, and LCS saddle points. In the cases of 2D flow, saddle points trace showed a rapid acceleration of the structure which indicates the shedding from the plate. For channel flow, saddle points trace shows that average structure convection speed exhibits a similar trend as a function of wall-normal distance as the mean velocity profile, and leads to statistical quantities of vortex dynamics. Dr. Jeff Eldredge and his research group at UCLA are gratefully acknowledged for sharing the database of simulation for the current research. This work was supported by the Air Force Office of Scientific Research under AFOSR Award No. FA9550-14-1-0210.

A unifying framework for ghost-free Lorentz-invariant Lagrangian field theories

NASA Astrophysics Data System (ADS)

Li, Wenliang

2018-04-01

We propose a framework for Lorentz-invariant Lagrangian field theories where Ostrogradsky's scalar ghosts could be absent. A key ingredient is the generalized Kronecker delta. The general Lagrangians are reformulated in the language of differential forms. The absence of higher order equations of motion for the scalar modes stems from the basic fact that every exact form is closed. The well-established Lagrangian theories for spin-0, spin-1, p-form, spin-2 fields have natural formulations in this framework. We also propose novel building blocks for Lagrangian field theories. Some of them are novel nonlinear derivative terms for spin-2 fields. It is nontrivial that Ostrogradsky's scalar ghosts are absent in these fully nonlinear theories.

Target Lagrangian kinematic simulation for particle-laden flows.

PubMed

Murray, S; Lightstone, M F; Tullis, S

2016-09-01

The target Lagrangian kinematic simulation method was motivated as a stochastic Lagrangian particle model that better synthesizes turbulence structure, relative to stochastic separated flow models. By this method, the trajectories of particles are constructed according to synthetic turbulent-like fields, which conform to a target Lagrangian integral timescale. In addition to recovering the expected Lagrangian properties of fluid tracers, this method is shown to reproduce the crossing trajectories and continuity effects, in agreement with an experimental benchmark.

An updated Lagrangian particle hydrodynamics (ULPH) for Newtonian fluids

NASA Astrophysics Data System (ADS)

Tu, Qingsong; Li, Shaofan

2017-11-01

In this work, we have developed an updated Lagrangian particle hydrodynamics (ULPH) for Newtonian fluid. Unlike the smoothed particle hydrodynamics, the non-local particle hydrodynamics formulation proposed here is consistent and convergence. Unlike the state-based peridynamics, the discrete particle dynamics proposed here has no internal material bond between particles, and it is not formulated with respect to initial or a fixed referential configuration. In specific, we have shown that (1) the non-local update Lagrangian particle hydrodynamics formulation converges to the conventional local fluid mechanics formulation; (2) the non-local updated Lagrangian particle hydrodynamics can capture arbitrary flow discontinuities without any changes in the formulation, and (3) the proposed non-local particle hydrodynamics is computationally efficient and robust.

Near-Surface Monsoonal Circulation of the Vietnam East Sea from Lagrangian Drifters

DTIC Science & Technology

2015-09-30

Sea from Lagrangian Drifters Luca Centurioni Scripps Institution of Oceanography 9500 Gilman Drive Mail Code 0213 La Jolla, California 92103...Contribute to the study of coastal and open ocean current systems in sparsely sampled regions such us the South China Sea (SCS), using a Lagrangian ...We intend to make new Lagrangian and Eulerian observations to measure the seasonal circulation 1) in the coastal waters of Vietnam and 2) in the SCS

Microscopic Lagrangian description of warm plasmas. I - Linear wave propagation. II - Nonlinear wave interactions

NASA Technical Reports Server (NTRS)

Kim, H.; Crawford, F. W.

1977-01-01

It is pointed out that the conventional iterative analysis of nonlinear plasma wave phenomena, which involves a direct use of Maxwell's equations and the equations describing the particle dynamics, leads to formidable theoretical and algebraic complexities, especially for warm plasmas. As an effective alternative, the Lagrangian method may be applied. It is shown how this method may be used in the microscopic description of small-signal wave propagation and in the study of nonlinear wave interactions. The linear theory is developed for an infinite, homogeneous, collisionless, warm magnetoplasma. A summary is presented of a perturbation expansion scheme described by Galloway and Kim (1971), and Lagrangians to third order in perturbation are considered. Attention is given to the averaged-Lagrangian density, the action-transfer and coupled-mode equations, and the general solution of the coupled-mode equations.

L-fuzzy fixed points theorems for L-fuzzy mappings via βℱL-admissible pair.

PubMed

Rashid, Maliha; Azam, Akbar; Mehmood, Nayyar

2014-01-01

We define the concept of βℱL-admissible for a pair of L-fuzzy mappings and establish the existence of common L-fuzzy fixed point theorem. Our result generalizes some useful results in the literature. We provide an example to support our result.

The research of the coupled orbital-attitude controlled motion of celestial body in the neighborhood of the collinear libration point L1

NASA Astrophysics Data System (ADS)

Shmyrov, A.; Shmyrov, V.; Shymanchuk, D.

2017-10-01

This article considers the motion of a celestial body within the restricted three-body problem of the Sun-Earth system. The equations of controlled coupled attitude-orbit motion in the neighborhood of collinear libration point L1 are investigated. The translational orbital motion of a celestial body is described using Hill's equations of circular restricted three-body problem of the Sun-Earth system. Rotational orbital motion is described using Euler's dynamic equations and quaternion kinematic equation. We investigate the problem of stability of celestial body rotational orbital motion in relative equilibrium positions and stabilization of celestial body rotational orbital motion with proposed control laws in the neighborhood of collinear libration point L1. To study stabilization problem, Lyapunov function is constructed in the form of the sum of the kinetic energy and special "kinematic function" of the Rodriguez-Hamiltonian parameters. Numerical modeling of the controlled rotational motion of a celestial body at libration point L1 is carried out. The numerical characteristics of the control parameters and rotational motion are given.

Comparison of PD-L1 mRNA Expression Measured with the CheckPoint Typer® Assay with PD-L1 Protein Expression Assessed with Immunohistochemistry in Non-small Cell Lung Cancer.

PubMed

Erber, Ramona; Stöhr, Robert; Herlein, Stefanie; Giedl, Claudia; Rieker, Ralf Joachim; Fuchs, Florian; Ficker, Joachim H; Hartmann, Arndt; Veltrup, Elke; Wirtz, Ralph M; Brueckl, Wolfgang M

2017-12-01

Immunohistochemical (IHC) assessment of programmed death-ligand 1 (PD-L1) in non-small cell lung cancer (NSCLC) has become important since the development of anti-PD-1/-PD-L1 directed drugs. Various PD-L1 antibodies and cut-offs have been used in different trials to predict response to these drugs, thus comparison of those studies is difficult. We compared PD-L1 mRNA expression measured by RT-qPCR with PD-L1 protein expression evaluated by IHC. Moreover, we investigated the impact of different tumour tissue acquisition methods on the reliability of PD-L1 measurement techniques. NSCLC cases (N=22), including n=9 mediastinal lymph node biopsies acquired by endobronchial ultrasound-guided transbronchial needle aspiration (EBUS-TBNA) and n=5 metastases, were evaluated prospectively for PD-L1 protein on tumor cells (TC) and immune cells (IC) using E1L3N and 28-8 antibodies and PD-L1 mRNA using the CheckPoint TYPER® assay. In primary NSCLC tissues, agreement between PD-L1 mRNA and TC staining using the 28-8 antibody was excellent (ĸ=0.85, p=0.0002). Comparing both PD-L1 antibodies against each other showed a kappa value of 0.58 (p=0.0106). In EBUS-TBNA, PD-L1 mRNA correlated perfectly with the 28-8 antibody (ĸ=1.0, p=0.0023). PD-L1 mRNA levels significantly differed when comparing 28-8 TC staining of tumours >49% with 1-49% and 0% (p=0.0040; p=0.0081, respectively). In metastatic lesions, differences between PD-L1 mRNA and IHC became apparent (ĸ=0.2, p=0.2525). Testing of PD-L1 mRNA and 28-8 IHC showed an excellent agreement in NSCLC samples including mediastinal lymph node biopsies. Since PD-L1 expression in >50% TC detected by 28-8 IHC can be reliably detected by RT-qPCR, quantitative PD-L1 mRNA determination should be considered as an alternative to IHC as there is no interobserver variability in RNA results. Copyright© 2017, International Institute of Anticancer Research (Dr. George J. Delinasios), All rights reserved.

L-Fuzzy Fixed Points Theorems for L-Fuzzy Mappings via β ℱL-Admissible Pair

PubMed Central

Rashid, Maliha; Azam, Akbar

2014-01-01

We define the concept of β ℱL-admissible for a pair of L-fuzzy mappings and establish the existence of common L-fuzzy fixed point theorem. Our result generalizes some useful results in the literature. We provide an example to support our result. PMID:24688441

Special Bohr-Sommerfeld Lagrangian submanifolds

NASA Astrophysics Data System (ADS)

Tyurin, N. A.

2016-12-01

We introduce a new notion in symplectic geometry, that of speciality for Lagrangian submanifolds satisfying the Bohr- Sommerfeld condition. We show that it enables one to construct finite-dimensional moduli spaces of special Bohr- Sommerfeld Lagrangian submanifolds with respect to any ample line bundle on an algebraic variety with a Hodge metric regarded as the symplectic form. This construction can be used to study mirror symmetry.

A high-order Lagrangian-decoupling method for the incompressible Navier-Stokes equations

NASA Technical Reports Server (NTRS)

Ho, Lee-Wing; Maday, Yvon; Patera, Anthony T.; Ronquist, Einar M.

1989-01-01

A high-order Lagrangian-decoupling method is presented for the unsteady convection-diffusion and incompressible Navier-Stokes equations. The method is based upon: (1) Lagrangian variational forms that reduce the convection-diffusion equation to a symmetric initial value problem; (2) implicit high-order backward-differentiation finite-difference schemes for integration along characteristics; (3) finite element or spectral element spatial discretizations; and (4) mesh-invariance procedures and high-order explicit time-stepping schemes for deducing function values at convected space-time points. The method improves upon previous finite element characteristic methods through the systematic and efficient extension to high order accuracy, and the introduction of a simple structure-preserving characteristic-foot calculation procedure which is readily implemented on modern architectures. The new method is significantly more efficient than explicit-convection schemes for the Navier-Stokes equations due to the decoupling of the convection and Stokes operators and the attendant increase in temporal stability. Numerous numerical examples are given for the convection-diffusion and Navier-Stokes equations for the particular case of a spectral element spatial discretization.

Computing eddy-driven effective diffusivity using Lagrangian particles

DOE PAGES

Wolfram, Phillip J.; Ringler, Todd D.

2017-08-14

A novel method to derive effective diffusivity from Lagrangian particle trajectory data sets is developed and then analyzed relative to particle-derived meridional diffusivity for eddy-driven mixing in an idealized circumpolar current. Quantitative standard dispersion- and transport-based mixing diagnostics are defined, compared and contrasted to motivate the computation and use of effective diffusivity derived from Lagrangian particles. We compute the effective diffusivity by first performing scalar transport on Lagrangian control areas using stored trajectories computed from online Lagrangian In-situ Global High-performance particle Tracking (LIGHT) using the Model for Prediction Across Scales Ocean (MPAS-O). Furthermore, the Lagrangian scalar transport scheme is comparedmore » against an Eulerian scalar transport scheme. Spatially-variable effective diffusivities are computed from resulting time-varying cumulative concentrations that vary as a function of cumulative area. The transport-based Eulerian and Lagrangian effective diffusivity diagnostics are found to be qualitatively consistent with the dispersion-based diffusivity. All diffusivity estimates show a region of increased subsurface diffusivity within the core of an idealized circumpolar current and results are within a factor of two of each other. The Eulerian and Lagrangian effective diffusivities are most similar; smaller and more spatially diffused values are obtained with the dispersion-based diffusivity computed with particle clusters.« less

Computing eddy-driven effective diffusivity using Lagrangian particles

DOE Office of Scientific and Technical Information (OSTI.GOV)

Wolfram, Phillip J.; Ringler, Todd D.

A novel method to derive effective diffusivity from Lagrangian particle trajectory data sets is developed and then analyzed relative to particle-derived meridional diffusivity for eddy-driven mixing in an idealized circumpolar current. Quantitative standard dispersion- and transport-based mixing diagnostics are defined, compared and contrasted to motivate the computation and use of effective diffusivity derived from Lagrangian particles. We compute the effective diffusivity by first performing scalar transport on Lagrangian control areas using stored trajectories computed from online Lagrangian In-situ Global High-performance particle Tracking (LIGHT) using the Model for Prediction Across Scales Ocean (MPAS-O). Furthermore, the Lagrangian scalar transport scheme is comparedmore » against an Eulerian scalar transport scheme. Spatially-variable effective diffusivities are computed from resulting time-varying cumulative concentrations that vary as a function of cumulative area. The transport-based Eulerian and Lagrangian effective diffusivity diagnostics are found to be qualitatively consistent with the dispersion-based diffusivity. All diffusivity estimates show a region of increased subsurface diffusivity within the core of an idealized circumpolar current and results are within a factor of two of each other. The Eulerian and Lagrangian effective diffusivities are most similar; smaller and more spatially diffused values are obtained with the dispersion-based diffusivity computed with particle clusters.« less

Shear and shearless Lagrangian structures in compound channels

NASA Astrophysics Data System (ADS)

Enrile, F.; Besio, G.; Stocchino, A.

2018-03-01

Transport processes in a physical model of a natural stream with a composite cross-section (compound channel) are investigated by means of a Lagrangian analysis based on nonlinear dynamical system theory. Two-dimensional free surface Eulerian experimental velocity fields of a uniform flow in a compound channel form the basis for the identification of the so-called Lagrangian Coherent Structures. Lagrangian structures are recognized as the key features that govern particle trajectories. We seek for two particular class of Lagrangian structures: Shear and shearless structures. The former are generated whenever the shear dominates the flow whereas the latter behave as jet-cores. These two type of structures are detected as ridges and trenches of the Finite-Time Lyapunov Exponents fields, respectively. Besides, shearlines computed applying the geodesic theory of transport barriers mark Shear Lagrangian Coherent Structures. So far, the detection of these structures in real experimental flows has not been deeply investigated. Indeed, the present results obtained in a wide range of the controlling parameters clearly show a different behaviour depending on the shallowness of the flow. Shear and Shearless Lagrangian Structures detected from laboratory experiments clearly appear as the flow develops in shallow conditions. The presence of these Lagrangian Structures tends to fade in deep flow conditions.

Forms of null Lagrangians in field theories of continuum mechanics

NASA Astrophysics Data System (ADS)

Kovalev, V. A.; Radaev, Yu. N.

2012-02-01

The divergence representation of a null Lagrangian that is regular in a star-shaped domain is used to obtain its general expression containing field gradients of order ≤ 1 in the case of spacetime of arbitrary dimension. It is shown that for a static three-component field in the three-dimensional space, a null Lagrangian can contain up to 15 independent elements in total. The general form of a null Lagrangian in the four-dimensional Minkowski spacetime is obtained (the number of physical field variables is assumed arbitrary). A complete theory of the null Lagrangian for the n-dimensional spacetime manifold (including the four-dimensional Minkowski spacetime as a special case) is given. Null Lagrangians are then used as a basis for solving an important variational problem of an integrating factor. This problem involves searching for factors that depend on the spacetime variables, field variables, and their gradients and, for a given system of partial differential equations, ensure the equality between the scalar product of a vector multiplier by the system vector and some divergence expression for arbitrary field variables and, hence, allow one to formulate a divergence conservation law on solutions to the system.

Asymptotic-preserving Lagrangian approach for modeling anisotropic transport in magnetized plasmas

NASA Astrophysics Data System (ADS)

Chacon, Luis; Del-Castillo-Negrete, Diego

2012-03-01

Modeling electron transport in magnetized plasmas is extremely challenging due to the extreme anisotropy between parallel (to the magnetic field) and perpendicular directions (the transport-coefficient ratio χ/χ˜10^10 in fusion plasmas). Recently, a novel Lagrangian Green's function method has been proposedfootnotetextD. del-Castillo-Negrete, L. Chac'on, PRL, 106, 195004 (2011); D. del-Castillo-Negrete, L. Chac'on, Phys. Plasmas, submitted (2011) to solve the local and non-local purely parallel transport equation in general 3D magnetic fields. The approach avoids numerical pollution, is inherently positivity-preserving, and is scalable algorithmically (i.e., work per degree-of-freedom is grid-independent). In this poster, we discuss the extension of the Lagrangian Green's function approach to include perpendicular transport terms and sources. We present an asymptotic-preserving numerical formulation, which ensures a consistent numerical discretization temporally and spatially for arbitrary χ/χ ratios. We will demonstrate the potential of the approach with various challenging configurations, including the case of transport across a magnetic island in cylindrical geometry.

Optimal Sunshade Configurations for Space-Based Geoengineering near the Sun-Earth L1 Point.

PubMed

Sánchez, Joan-Pau; McInnes, Colin R

2015-01-01

Within the context of anthropogenic climate change, but also considering the Earth's natural climate variability, this paper explores the speculative possibility of large-scale active control of the Earth's radiative forcing. In particular, the paper revisits the concept of deploying a large sunshade or occulting disk at a static position near the Sun-Earth L1 Lagrange equilibrium point. Among the solar radiation management methods that have been proposed thus far, space-based concepts are generally seen as the least timely, albeit also as one of the most efficient. Large occulting structures could potentially offset all of the global mean temperature increase due to greenhouse gas emissions. This paper investigates optimal configurations of orbiting occulting disks that not only offset a global temperature increase, but also mitigate regional differences such as latitudinal and seasonal difference of monthly mean temperature. A globally resolved energy balance model is used to provide insights into the coupling between the motion of the occulting disks and the Earth's climate. This allows us to revise previous studies, but also, for the first time, to search for families of orbits that improve the efficiency of occulting disks at offsetting climate change on both global and regional scales. Although natural orbits exist near the L1 equilibrium point, their period does not match that required for geoengineering purposes, thus forced orbits were designed that require small changes to the disk attitude in order to control its motion. Finally, configurations of two occulting disks are presented which provide the same shading area as previously published studies, but achieve reductions of residual latitudinal and seasonal temperature changes.

Optimal Sunshade Configurations for Space-Based Geoengineering near the Sun-Earth L1 Point

PubMed Central

Sánchez, Joan-Pau; McInnes, Colin R.

2015-01-01

Within the context of anthropogenic climate change, but also considering the Earth’s natural climate variability, this paper explores the speculative possibility of large-scale active control of the Earth’s radiative forcing. In particular, the paper revisits the concept of deploying a large sunshade or occulting disk at a static position near the Sun-Earth L1 Lagrange equilibrium point. Among the solar radiation management methods that have been proposed thus far, space-based concepts are generally seen as the least timely, albeit also as one of the most efficient. Large occulting structures could potentially offset all of the global mean temperature increase due to greenhouse gas emissions. This paper investigates optimal configurations of orbiting occulting disks that not only offset a global temperature increase, but also mitigate regional differences such as latitudinal and seasonal difference of monthly mean temperature. A globally resolved energy balance model is used to provide insights into the coupling between the motion of the occulting disks and the Earth’s climate. This allows us to revise previous studies, but also, for the first time, to search for families of orbits that improve the efficiency of occulting disks at offsetting climate change on both global and regional scales. Although natural orbits exist near the L1 equilibrium point, their period does not match that required for geoengineering purposes, thus forced orbits were designed that require small changes to the disk attitude in order to control its motion. Finally, configurations of two occulting disks are presented which provide the same shading area as previously published studies, but achieve reductions of residual latitudinal and seasonal temperature changes. PMID:26309047

Lagrangian particle method for compressible fluid dynamics

DOE Office of Scientific and Technical Information (OSTI.GOV)

Samulyak, Roman; Wang, Xingyu; Chen, Hsin -Chiang

A new Lagrangian particle method for solving Euler equations for compressible inviscid fluid or gas flows is proposed. Similar to smoothed particle hydrodynamics (SPH), the method represents fluid cells with Lagrangian particles and is suitable for the simulation of complex free surface / multi-phase flows. The main contributions of our method, which is different from SPH in all other aspects, are (a) significant improvement of approximation of differential operators based on a polynomial fit via weighted least squares approximation and the convergence of prescribed order, (b) a second-order particle-based algorithm that reduces to the first-order upwind method at local extremalmore » points, providing accuracy and long term stability, and (c) more accurate resolution of entropy discontinuities and states at free inter-faces. While the method is consistent and convergent to a prescribed order, the conservation of momentum and energy is not exact and depends on the convergence order . The method is generalizable to coupled hyperbolic-elliptic systems. As a result, numerical verification tests demonstrating the convergence order are presented as well as examples of complex multiphase flows.« less

Lagrangian descriptors of driven chemical reaction manifolds.

PubMed

Craven, Galen T; Junginger, Andrej; Hernandez, Rigoberto

2017-08-01

The persistence of a transition state structure in systems driven by time-dependent environments allows the application of modern reaction rate theories to solution-phase and nonequilibrium chemical reactions. However, identifying this structure is problematic in driven systems and has been limited by theories built on series expansion about a saddle point. Recently, it has been shown that to obtain formally exact rates for reactions in thermal environments, a transition state trajectory must be constructed. Here, using optimized Lagrangian descriptors [G. T. Craven and R. Hernandez, Phys. Rev. Lett. 115, 148301 (2015)PRLTAO0031-900710.1103/PhysRevLett.115.148301], we obtain this so-called distinguished trajectory and the associated moving reaction manifolds on model energy surfaces subject to various driving and dissipative conditions. In particular, we demonstrate that this is exact for harmonic barriers in one dimension and this verification gives impetus to the application of Lagrangian descriptor-based methods in diverse classes of chemical reactions. The development of these objects is paramount in the theory of reaction dynamics as the transition state structure and its underlying network of manifolds directly dictate reactivity and selectivity.

Lagrangian particle method for compressible fluid dynamics

DOE PAGES

Samulyak, Roman; Wang, Xingyu; Chen, Hsin -Chiang

2018-02-09

A new Lagrangian particle method for solving Euler equations for compressible inviscid fluid or gas flows is proposed. Similar to smoothed particle hydrodynamics (SPH), the method represents fluid cells with Lagrangian particles and is suitable for the simulation of complex free surface / multi-phase flows. The main contributions of our method, which is different from SPH in all other aspects, are (a) significant improvement of approximation of differential operators based on a polynomial fit via weighted least squares approximation and the convergence of prescribed order, (b) a second-order particle-based algorithm that reduces to the first-order upwind method at local extremalmore » points, providing accuracy and long term stability, and (c) more accurate resolution of entropy discontinuities and states at free inter-faces. While the method is consistent and convergent to a prescribed order, the conservation of momentum and energy is not exact and depends on the convergence order . The method is generalizable to coupled hyperbolic-elliptic systems. As a result, numerical verification tests demonstrating the convergence order are presented as well as examples of complex multiphase flows.« less

Dehydration and Lagrangian Cold Point in the extratropical Tropopause region

NASA Astrophysics Data System (ADS)

Hoor, P.; Wernli, H.

2012-04-01

The tropopause region of the tropics and extratropics is sensitive to modifications of the radiation budget through changes of radiatively active substances like ozone and water vapour. Both may also modify the temperature structure and the strengths of the tropopause inversion layer (TIL). Stratospheric water vapour is mainly controlled by dehydration in the tropics. Ascending air masses encounter their minimum temperature in the TTL region (tropical tropopause layer) which determines the water vapour fraction which enters the stratosphere. In the lowermost stratosphere of the extratropics however, the tropical signal might be lost due to mixing with airmasses which crossed the tropopause (TST: troposphere to stratosphere) at higher temperatures, therefore carrying more water vapour to the extratropical stratosphere. We investigate statistical 90 day backward trajectories to investigate the role of dehydration at the extratropical tropopause for the water vapour budget at the tropopause at mid and high latitudes. We use a set of 800000 trajectories for summer and winter, respectively, on the basis of ECMWF-T799L91 operational data (kinematic wind fields). We analyze the trajectories for the time and locations of their cold point and TST. Our results indicate that : 1) TST and dehydration occur at different locations 2) Dehydration occurs in general before trajectories enter the stratosphere 3) Dehydration of TST trajectories can occur in northern winter after TST in the region of high tropopauses over Siberia

A Lagrangian subgrid-scale model with dynamic estimation of Lagrangian time scale for large eddy simulation of complex flows

NASA Astrophysics Data System (ADS)

Verma, Aman; Mahesh, Krishnan

2012-08-01

The dynamic Lagrangian averaging approach for the dynamic Smagorinsky model for large eddy simulation is extended to an unstructured grid framework and applied to complex flows. The Lagrangian time scale is dynamically computed from the solution and does not need any adjustable parameter. The time scale used in the standard Lagrangian model contains an adjustable parameter θ. The dynamic time scale is computed based on a "surrogate-correlation" of the Germano-identity error (GIE). Also, a simple material derivative relation is used to approximate GIE at different events along a pathline instead of Lagrangian tracking or multi-linear interpolation. Previously, the time scale for homogeneous flows was computed by averaging along directions of homogeneity. The present work proposes modifications for inhomogeneous flows. This development allows the Lagrangian averaged dynamic model to be applied to inhomogeneous flows without any adjustable parameter. The proposed model is applied to LES of turbulent channel flow on unstructured zonal grids at various Reynolds numbers. Improvement is observed when compared to other averaging procedures for the dynamic Smagorinsky model, especially at coarse resolutions. The model is also applied to flow over a cylinder at two Reynolds numbers and good agreement with previous computations and experiments is obtained. Noticeable improvement is obtained using the proposed model over the standard Lagrangian model. The improvement is attributed to a physically consistent Lagrangian time scale. The model also shows good performance when applied to flow past a marine propeller in an off-design condition; it regularizes the eddy viscosity and adjusts locally to the dominant flow features.

Lagrangian formulation for penny-shaped and Perkins-Kern geometry models

DOE Office of Scientific and Technical Information (OSTI.GOV)

Lee, W.S.

1989-09-01

This paper discusses basic theories for vertical penny-shaped and Perkins-Kern (PK) geometry models developed with a Lagrangian formulation combined with a virtual-work analysis. The Lagrangian formulation yields a pair of nonlinear equations in R/sub f/ or L/sub f/ and b/sub f/, the fracture radius or length and half-width. By introduction of a virtual-work analysis, a simple equation is obtained that can be solved numerically. This equation is written in a form that can be used to determine fracture geometry when the fluid-loss coefficient of the fracturing fluid is known. Also, this equation, coupled with a material-balance equation after shut-in, canmore » be used to analyze pressure-decline data after shut-in to determine the effective fluid-loss coefficient and fracture geometry.« less

Coupled Eulerian-Lagrangian transport of large debris by tsunamis

NASA Astrophysics Data System (ADS)

Conde, Daniel A. S.; Ferreira, Rui M. L.; Sousa Oliveira, Carlos

2016-04-01

Tsunamis are notorious for the large disruption they can cause on coastal environments, not only due to the imparted momentum of the incoming wave but also due to its capacity to transport large quantities of solid debris, either from natural or human-made sources, over great distances. A 2DH numerical model under development at CERIS-IST (Ferreira et al., 2009; Conde, 2013) - STAV2D - capable of simulating solid transport in both Eulerian and Lagrangian paradigms will be used to assess the relevance of Lagrangian-Eulerian coupling when modelling the transport of solid debris by tsunamis. The model has been previously validated and applied to tsunami scenarios (Conde, 2013), being well-suited for overland tsunami propagation and capable of handling morphodynamic changes in estuaries and seashores. The discretization scheme is an explicit Finite Volume technique employing flux-vector splitting and a reviewed Roe-Riemann solver. Source term formulations are employed in a semi-implicit way, including the two-way coupling of the Lagrangian and Eulerian solvers by means of conservative mass and momentum transfers between fluid and solid phases. The model was applied to Sines Port, a major commercial port in Portugal, where two tsunamigenic scenarios are considered: an 8.5 Mw scenario, consistent with the Great Lisbon Earthquake and Tsunami of the 1st November 1755 (Baptista, 2009), and an hypothetical 9.5 Mw worst-case scenario based on the same historical event. Open-ocean propagation of these scenarios were simulated with GeoClaw model from ClawPack (Leveque, 2011). Following previous efforts on the modelling of debris transport by tsunamis in seaports (Conde, 2015), this work discusses the sensitivity of the obtained results with respect to the phenomenological detail of the employed Eulerian-Lagrangian formulation and the resolution of the mesh used in the Eulerian solver. The results have shown that the fluid to debris mass ratio is the key parameter regarding the

Application of Quasi-Lagrangian Diagnostics and FGGE Data in a Study of East-Coast Cyclogenesis.

DTIC Science & Technology

1981-09-01

QUASI-LAGRANGIAN UDIAGNOSTICS AND FGGE DATA IN A STUDY OF EAST-COAST CYCLOGENESIS by Donald A. Roman September 1981 Thesis Advisor: Dr. Carlyle H. Wash...REPORT DOCUMENTATION PAGE DEVFOR COMPLETING FORM I.M9PR UM19 2. GOV? AccaS "@ "mO L.VIPICUTS CA? ALOG MUUNW S TyP =F NOo’BP io OE4. TITLC (nd Sisle . ye i...n.J000 Application of Quasi-Lagrangian Master’s Thesis Diagnostics and FGGE Data in a Study 0#4 September 1981 East-Coast Cyclogenesis C PlaPranIWG

The Cerebro-oculo-facio-skeletal Syndrome Point Mutation F231L in the ERCC1 DNA Repair Protein Causes Dissociation of the ERCC1-XPF Complex*

PubMed Central

Faridounnia, Maryam; Wienk, Hans; Kovačič, Lidija; Folkers, Gert E.; Jaspers, Nicolaas G. J.; Kaptein, Robert; Hoeijmakers, Jan H. J.; Boelens, Rolf

2015-01-01

The ERCC1-XPF heterodimer, a structure-specific DNA endonuclease, is best known for its function in the nucleotide excision repair (NER) pathway. The ERCC1 point mutation F231L, located at the hydrophobic interaction interface of ERCC1 (excision repair cross-complementation group 1) and XPF (xeroderma pigmentosum complementation group F), leads to severe NER pathway deficiencies. Here, we analyze biophysical properties and report the NMR structure of the complex of the C-terminal tandem helix-hairpin-helix domains of ERCC1-XPF that contains this mutation. The structures of wild type and the F231L mutant are very similar. The F231L mutation results in only a small disturbance of the ERCC1-XPF interface, where, in contrast to Phe231, Leu231 lacks interactions stabilizing the ERCC1-XPF complex. One of the two anchor points is severely distorted, and this results in a more dynamic complex, causing reduced stability and an increased dissociation rate of the mutant complex as compared with wild type. These data provide a biophysical explanation for the severe NER deficiencies caused by this mutation. PMID:26085086

The Cerebro-oculo-facio-skeletal Syndrome Point Mutation F231L in the ERCC1 DNA Repair Protein Causes Dissociation of the ERCC1-XPF Complex.

PubMed

Faridounnia, Maryam; Wienk, Hans; Kovačič, Lidija; Folkers, Gert E; Jaspers, Nicolaas G J; Kaptein, Robert; Hoeijmakers, Jan H J; Boelens, Rolf

2015-08-14

The ERCC1-XPF heterodimer, a structure-specific DNA endonuclease, is best known for its function in the nucleotide excision repair (NER) pathway. The ERCC1 point mutation F231L, located at the hydrophobic interaction interface of ERCC1 (excision repair cross-complementation group 1) and XPF (xeroderma pigmentosum complementation group F), leads to severe NER pathway deficiencies. Here, we analyze biophysical properties and report the NMR structure of the complex of the C-terminal tandem helix-hairpin-helix domains of ERCC1-XPF that contains this mutation. The structures of wild type and the F231L mutant are very similar. The F231L mutation results in only a small disturbance of the ERCC1-XPF interface, where, in contrast to Phe(231), Leu(231) lacks interactions stabilizing the ERCC1-XPF complex. One of the two anchor points is severely distorted, and this results in a more dynamic complex, causing reduced stability and an increased dissociation rate of the mutant complex as compared with wild type. These data provide a biophysical explanation for the severe NER deficiencies caused by this mutation. © 2015 by The American Society for Biochemistry and Molecular Biology, Inc.

Parallel computing using a Lagrangian formulation

NASA Technical Reports Server (NTRS)

Liou, May-Fun; Loh, Ching Yuen

1991-01-01

A new Lagrangian formulation of the Euler equation is adopted for the calculation of 2-D supersonic steady flow. The Lagrangian formulation represents the inherent parallelism of the flow field better than the common Eulerian formulation and offers a competitive alternative on parallel computers. The implementation of the Lagrangian formulation on the Thinking Machines Corporation CM-2 Computer is described. The program uses a finite volume, first-order Godunov scheme and exhibits high accuracy in dealing with multidimensional discontinuities (slip-line and shock). By using this formulation, a better than six times speed-up was achieved on a 8192-processor CM-2 over a single processor of a CRAY-2.

Parallel computing using a Lagrangian formulation

NASA Technical Reports Server (NTRS)

Liou, May-Fun; Loh, Ching-Yuen

1992-01-01

This paper adopts a new Lagrangian formulation of the Euler equation for the calculation of two dimensional supersonic steady flow. The Lagrangian formulation represents the inherent parallelism of the flow field better than the common Eulerian formulation and offers a competitive alternative on parallel computers. The implementation of the Lagrangian formulation on the Thinking Machines Corporation CM-2 Computer is described. The program uses a finite volume, first-order Godunov scheme and exhibits high accuracy in dealing with multidimensional discontinuities (slip-line and shock). By using this formulation, we have achieved better than six times speed-up on a 8192-processor CM-2 over a single processor of a CRAY-2.

Modeling the Global Coronal Field with Simulated Synoptic Magnetograms from Earth and the Lagrange Points L3, L4, and L5

NASA Astrophysics Data System (ADS)

Petrie, Gordon; Pevtsov, Alexei; Schwarz, Andrew; DeRosa, Marc

2018-06-01

The solar photospheric magnetic flux distribution is key to structuring the global solar corona and heliosphere. Regular full-disk photospheric magnetogram data are therefore essential to our ability to model and forecast heliospheric phenomena such as space weather. However, our spatio-temporal coverage of the photospheric field is currently limited by our single vantage point at/near Earth. In particular, the polar fields play a leading role in structuring the large-scale corona and heliosphere, but each pole is unobservable for {>} 6 months per year. Here we model the possible effect of full-disk magnetogram data from the Lagrange points L4 and L5, each extending longitude coverage by 60°. Adding data also from the more distant point L3 extends the longitudinal coverage much further. The additional vantage points also improve the visibility of the globally influential polar fields. Using a flux-transport model for the solar photospheric field, we model full-disk observations from Earth/L1, L3, L4, and L5 over a solar cycle, construct synoptic maps using a novel weighting scheme adapted for merging magnetogram data from multiple viewpoints, and compute potential-field models for the global coronal field. Each additional viewpoint brings the maps and models into closer agreement with the reference field from the flux-transport simulation, with particular improvement at polar latitudes, the main source of the fast solar wind.

Lagrangian analysis. Modern tool of the dynamics of solids

NASA Astrophysics Data System (ADS)

Cagnoux, J.; Chartagnac, P.; Hereil, P.; Perez, M.; Seaman, L.

Lagrangian specificity of the required measurements is assured by the fact that a transducer enclosed within a solid material is necessarily linked in motion to the particles of the material which surround it. This Lagrangian instrumentation is described in the second chapter. The authors are concerned with the techniques considered today to be the most effective. These are, for stress : piezoresistive gauges (50 Ω and low impedance) and piezoelectric techniques (PVF2 gauges, quartz transducers) ; and for particle velocity : electromagnetic gauges, VISAR and IDL Doppler laser interferometers. In each case both the physical principles as well as techniques of use are set out in detail. For the most part, the authors use their own experience to describe the calibration of these instrumentation systems and to compare their characteristics : measurement range, response time, accuracy, useful recording time, detection area... These characteristics should be taken into account by the physicist when he has to choose the instrumentation systems best adapted to the Lagrangian analysis he intends to apply to any given material. The discussion at the end of chapter 2 should guide his choice both for plane and spherical one-dimensional motions. The third chapter examines to what extent the accuracy of Lagrangian analysis is affected by the accuracies of the numerical analysis methods and experimental techniques. By means of a discussion of different cases of analysis, the authors want to make the reader aware of the different kinds of sources of errors that may be encountered. This work brings up to date the state of studies on Lagrangian analysis methods based on a wide review of bibliographical sources together with the contribution made to research in this field by the four authors themselves in the course of the last ten years. Le formage des métaux par explosif, la consolidation dynamique des poudres, la balistique terminale, l'abattage des roches par explosif, sont autant d

Dual domain material point method for multiphase flows

NASA Astrophysics Data System (ADS)

Zhang, Duan

2017-11-01

Although the particle-in-cell method was first invented in the 60's for fluid computations, one of its later versions, the material point method, is mostly used for solid calculations. Recent development of the multi-velocity formulations for multiphase flows and fluid-structure interactions requires the Lagrangian capability of the method be combined with Eulerian calculations for fluids. Because of different numerical representations of the materials, additional numerical schemes are needed to ensure continuity of the materials. New applications of the method to compute fluid motions have revealed numerical difficulties in various versions of the method. To resolve these difficulties, the dual domain material point method is introduced and improved. Unlike other particle based methods, the material point method uses both Lagrangian particles and Eulerian mesh, therefore it avoids direct communication between particles. With this unique property and the Lagrangian capability of the method, it is shown that a multiscale numerical scheme can be efficiently built based on the dual domain material point method. In this talk, the theoretical foundation of the method will be introduced. Numerical examples will be shown. Work sponsored by the next generation code project of LANL.

Lagrangian Assimilation of Satellite Data for Climate Studies in the Arctic

NASA Technical Reports Server (NTRS)

Lindsay, Ronald W.; Zhang, Jin-Lun; Stern, Harry

2004-01-01

Under this grant we have developed and tested a new Lagrangian model of sea ice. A Lagrangian model keeps track of material parcels as they drift in the model domain. Besides providing a natural framework for the assimilation of Lagrangian data, it has other advantages: 1) a model that follows material elements is well suited for a medium such as sea ice in which an element retains its identity for a long period of time; 2) model cells can be added or dropped as needed, allowing the spatial resolution to be increased in areas of high variability or dense observations; 3) ice from particular regions, such as the marginal seas, can be marked and traced for a long time; and 4) slip lines in the ice motion are accommodated more naturally because there is no internal grid. Our work makes use of these strengths of the Lagrangian formulation.

A Vertically Lagrangian Finite-Volume Dynamical Core for Global Models

NASA Technical Reports Server (NTRS)

Lin, Shian-Jiann

2003-01-01

A finite-volume dynamical core with a terrain-following Lagrangian control-volume discretization is described. The vertically Lagrangian discretization reduces the dimensionality of the physical problem from three to two with the resulting dynamical system closely resembling that of the shallow water dynamical system. The 2D horizontal-to-Lagrangian-surface transport and dynamical processes are then discretized using the genuinely conservative flux-form semi-Lagrangian algorithm. Time marching is split- explicit, with large-time-step for scalar transport, and small fractional time step for the Lagrangian dynamics, which permits the accurate propagation of fast waves. A mass, momentum, and total energy conserving algorithm is developed for mapping the state variables periodically from the floating Lagrangian control-volume to an Eulerian terrain-following coordinate for dealing with physical parameterizations and to prevent severe distortion of the Lagrangian surfaces. Deterministic baroclinic wave growth tests and long-term integrations using the Held-Suarez forcing are presented. Impact of the monotonicity constraint is discussed.

Adjoint of the global Eulerian-Lagrangian coupled atmospheric transport model (A-GELCA v1.0): development and validation

NASA Astrophysics Data System (ADS)

Belikov, Dmitry A.; Maksyutov, Shamil; Yaremchuk, Alexey; Ganshin, Alexander; Kaminski, Thomas; Blessing, Simon; Sasakawa, Motoki; Gomez-Pelaez, Angel J.; Starchenko, Alexander

2016-02-01

We present the development of the Adjoint of the Global Eulerian-Lagrangian Coupled Atmospheric (A-GELCA) model that consists of the National Institute for Environmental Studies (NIES) model as an Eulerian three-dimensional transport model (TM), and FLEXPART (FLEXible PARTicle dispersion model) as the Lagrangian Particle Dispersion Model (LPDM). The forward tangent linear and adjoint components of the Eulerian model were constructed directly from the original NIES TM code using an automatic differentiation tool known as TAF (Transformation of Algorithms in Fortran; http://www.FastOpt.com, with additional manual pre- and post-processing aimed at improving transparency and clarity of the code and optimizing the performance of the computing, including MPI (Message Passing Interface). The Lagrangian component did not require any code modification, as LPDMs are self-adjoint and track a significant number of particles backward in time in order to calculate the sensitivity of the observations to the neighboring emission areas. The constructed Eulerian adjoint was coupled with the Lagrangian component at a time boundary in the global domain. The simulations presented in this work were performed using the A-GELCA model in forward and adjoint modes. The forward simulation shows that the coupled model improves reproduction of the seasonal cycle and short-term variability of CO2. Mean bias and standard deviation for five of the six Siberian sites considered decrease roughly by 1 ppm when using the coupled model. The adjoint of the Eulerian model was shown, through several numerical tests, to be very accurate (within machine epsilon with mismatch around to ±6 e-14) compared to direct forward sensitivity calculations. The developed adjoint of the coupled model combines the flux conservation and stability of an Eulerian discrete adjoint formulation with the flexibility, accuracy, and high resolution of a Lagrangian backward trajectory formulation. A-GELCA will be incorporated

Lagrangian motion, coherent structures, and lines of persistent material strain.

PubMed

Samelson, R M

2013-01-01

Lagrangian motion in geophysical fluids may be strongly influenced by coherent structures that support distinct regimes in a given flow. The problems of identifying and demarcating Lagrangian regime boundaries associated with dynamical coherent structures in a given velocity field can be studied using approaches originally developed in the context of the abstract geometric theory of ordinary differential equations. An essential insight is that when coherent structures exist in a flow, Lagrangian regime boundaries may often be indicated as material curves on which the Lagrangian-mean principal-axis strain is large. This insight is the foundation of many numerical techniques for identifying such features in complex observed or numerically simulated ocean flows. The basic theoretical ideas are illustrated with a simple, kinematic traveling-wave model. The corresponding numerical algorithms for identifying candidate Lagrangian regime boundaries and lines of principal Lagrangian strain (also called Lagrangian coherent structures) are divided into parcel and bundle schemes; the latter include the finite-time and finite-size Lyapunov exponent/Lagrangian strain (FTLE/FTLS and FSLE/FSLS) metrics. Some aspects and results of oceanographic studies based on these approaches are reviewed, and the results are discussed in the context of oceanographic observations of dynamical coherent structures.

Comparison of updated Lagrangian FEM with arbitrary Lagrangian Eulerian method for 3D thermo-mechanical extrusion of a tube profile

NASA Astrophysics Data System (ADS)

Kronsteiner, J.; Horwatitsch, D.; Zeman, K.

2017-10-01

Thermo-mechanical numerical modelling and simulation of extrusion processes faces several serious challenges. Large plastic deformations in combination with a strong coupling of thermal with mechanical effects leads to a high numerical demand for the solution as well as for the handling of mesh distortions. The two numerical methods presented in this paper also reflect two different ways to deal with mesh distortions. Lagrangian Finite Element Methods (FEM) tackle distorted elements by building a new mesh (called re-meshing) whereas Arbitrary Lagrangian Eulerian (ALE) methods use an "advection" step to remap the solution from the distorted to the undistorted mesh. Another difference between conventional Lagrangian and ALE methods is the separate treatment of material and mesh in ALE, allowing the definition of individual velocity fields. In theory, an ALE formulation contains the Eulerian formulation as a subset to the Lagrangian description of the material. The investigations presented in this paper were dealing with the direct extrusion of a tube profile using EN-AW 6082 aluminum alloy and a comparison of experimental with Lagrangian and ALE results. The numerical simulations cover the billet upsetting and last until one third of the billet length is extruded. A good qualitative correlation of experimental and numerical results could be found, however, major differences between Lagrangian and ALE methods concerning thermo-mechanical coupling lead to deviations in the thermal results.

Gravity, Time, and Lagrangians

NASA Astrophysics Data System (ADS)

Huggins, Elisha

2010-11-01

Feynman mentioned to us that he understood a topic in physics if he could explain it to a college freshman, a high school student, or a dinner guest. Here we will discuss two topics that took us a while to get to that level. One is the relationship between gravity and time. The other is the minus sign that appears in the Lagrangian. (Why would one subtract potential energy from kinetic energy?) In this paper we discuss a thought experiment that relates gravity and time. Then we use a Feynman thought experiment to explain the minus sign in the Lagrangian. Our surprise was that these two topics are related.

An asymptotic-preserving Lagrangian algorithm for the time-dependent anisotropic heat transport equation

DOE Office of Scientific and Technical Information (OSTI.GOV)

Chacon, Luis; del-Castillo-Negrete, Diego; Hauck, Cory D.

2014-09-01

We propose a Lagrangian numerical algorithm for a time-dependent, anisotropic temperature transport equation in magnetized plasmas in the large guide field regime. The approach is based on an analytical integral formal solution of the parallel (i.e., along the magnetic field) transport equation with sources, and it is able to accommodate both local and non-local parallel heat flux closures. The numerical implementation is based on an operator-split formulation, with two straightforward steps: a perpendicular transport step (including sources), and a Lagrangian (field-line integral) parallel transport step. Algorithmically, the first step is amenable to the use of modern iterative methods, while themore » second step has a fixed cost per degree of freedom (and is therefore scalable). Accuracy-wise, the approach is free from the numerical pollution introduced by the discrete parallel transport term when the perpendicular to parallel transport coefficient ratio X ⊥ /X ∥ becomes arbitrarily small, and is shown to capture the correct limiting solution when ε = X⊥L 2 ∥/X1L 2 ⊥ → 0 (with L∥∙ L⊥ , the parallel and perpendicular diffusion length scales, respectively). Therefore, the approach is asymptotic-preserving. We demonstrate the capabilities of the scheme with several numerical experiments with varying magnetic field complexity in two dimensions, including the case of transport across a magnetic island.« less

Time-resolved large-scale volumetric pressure fields of an impinging jet from dense Lagrangian particle tracking

NASA Astrophysics Data System (ADS)

Huhn, F.; Schanz, D.; Manovski, P.; Gesemann, S.; Schröder, A.

2018-05-01

Time-resolved volumetric pressure fields are reconstructed from Lagrangian particle tracking with high seeding concentration using the Shake-The-Box algorithm in a perpendicular impinging jet flow with exit velocity U=4 m/s (Re˜ 36,000) and nozzle-plate spacing H/D=5. Helium-filled soap bubbles are used as tracer particles which are illuminated with pulsed LED arrays. A large measurement volume has been covered (cloud of tracked particles in a volume of 54 L, ˜ 180,000 particles). The reconstructed pressure field has been validated against microphone recordings at the wall with high correlation coefficients up to 0.88. In a reduced measurement volume (13 L), dense Lagrangian particle tracking is shown to be feasable up to the maximal possible jet velocity of U=16 m/s.

Water-Channel Estimation of Eulerian and Lagrangian Time Scales of the Turbulence in Idealized Two-Dimensional Urban Canopies

NASA Astrophysics Data System (ADS)

Di Bernardino, Annalisa; Monti, Paolo; Leuzzi, Giovanni; Querzoli, Giorgio

2017-11-01

Lagrangian and Eulerian statistics are obtained from a water-channel experiment of an idealized two-dimensional urban canopy flow in neutral conditions. The objective is to quantify the Eulerian (TE) and Lagrangian (TL) time scales of the turbulence above the canopy layer as well as to investigate their dependence on the aspect ratio of the canopy, AR, as the latter is the ratio of the width ( W) to the height ( H) of the canyon. Experiments are also conducted for the case of flat terrain, which can be thought of as equivalent to a classical one-directional shear flow. The values found for the Eulerian time scales on flat terrain are in agreement with previous numerical results found in the literature. It is found that both the streamwise and vertical components of the Lagrangian time scale, T_u^L and T_w^L , follow Raupach's linear law within the constant-flux layer. The same holds true for T_w^L in both the canopies analyzed (AR= 1 and AR= 2) and also for T_u^L when AR = 1. In contrast, for AR = 2, T_u^L follows Raupach's law only above z=2H. Below that level, T_u^L is nearly constant with height, showing at z=H a value approximately one order of magnitude greater than that found for AR = 1. It is shown that the assumption usually adopted for flat terrain, that β =TL/TE is proportional to the inverse of the turbulence intensity, also holds true even for the canopy flow in the constant-flux layer. In particular, γ /i_u fits well β _u =T_u^L /T_u^E in both the configurations by choosing γ to be 0.35 (here, i_u =σ _u / \\bar{u} , where \\bar{u} and σ _u are the mean and the root-mean-square of the streamwise velocity component, respectively). On the other hand, β _w =T_w^L /T_w^E follows approximately γ /i_w =0.65/( {σ _w /\\bar{u} } ) for z > 2H, irrespective of the AR value. The second main objective is to estimate other parameters of interest in dispersion studies, such as the eddy diffusivity of momentum (KT) and the Kolmogorov constant (C_0). It is

Lagrangian ocean analysis: Fundamentals and practices

DOE PAGES

van Sebille, Erik; Griffies, Stephen M.; Abernathey, Ryan; ...

2017-11-24

Lagrangian analysis is a powerful way to analyse the output of ocean circulation models and other ocean velocity data such as from altimetry. In the Lagrangian approach, large sets of virtual particles are integrated within the three-dimensional, time-evolving velocity fields. A variety of tools and methods for this purpose have emerged, over several decades. Here, we review the state of the art in the field of Lagrangian analysis of ocean velocity data, starting from a fundamental kinematic framework and with a focus on large-scale open ocean applications. Beyond the use of explicit velocity fields, we consider the influence of unresolvedmore » physics and dynamics on particle trajectories. We comprehensively list and discuss the tools currently available for tracking virtual particles. We then showcase some of the innovative applications of trajectory data, and conclude with some open questions and an outlook. Our overall goal of this review paper is to reconcile some of the different techniques and methods in Lagrangian ocean analysis, while recognising the rich diversity of codes that have and continue to emerge, and the challenges of the coming age of petascale computing.« less

Lagrangian ocean analysis: Fundamentals and practices

NASA Astrophysics Data System (ADS)

van Sebille, Erik; Griffies, Stephen M.; Abernathey, Ryan; Adams, Thomas P.; Berloff, Pavel; Biastoch, Arne; Blanke, Bruno; Chassignet, Eric P.; Cheng, Yu; Cotter, Colin J.; Deleersnijder, Eric; Döös, Kristofer; Drake, Henri F.; Drijfhout, Sybren; Gary, Stefan F.; Heemink, Arnold W.; Kjellsson, Joakim; Koszalka, Inga Monika; Lange, Michael; Lique, Camille; MacGilchrist, Graeme A.; Marsh, Robert; Mayorga Adame, C. Gabriela; McAdam, Ronan; Nencioli, Francesco; Paris, Claire B.; Piggott, Matthew D.; Polton, Jeff A.; Rühs, Siren; Shah, Syed H. A. M.; Thomas, Matthew D.; Wang, Jinbo; Wolfram, Phillip J.; Zanna, Laure; Zika, Jan D.

2018-01-01

Lagrangian analysis is a powerful way to analyse the output of ocean circulation models and other ocean velocity data such as from altimetry. In the Lagrangian approach, large sets of virtual particles are integrated within the three-dimensional, time-evolving velocity fields. Over several decades, a variety of tools and methods for this purpose have emerged. Here, we review the state of the art in the field of Lagrangian analysis of ocean velocity data, starting from a fundamental kinematic framework and with a focus on large-scale open ocean applications. Beyond the use of explicit velocity fields, we consider the influence of unresolved physics and dynamics on particle trajectories. We comprehensively list and discuss the tools currently available for tracking virtual particles. We then showcase some of the innovative applications of trajectory data, and conclude with some open questions and an outlook. The overall goal of this review paper is to reconcile some of the different techniques and methods in Lagrangian ocean analysis, while recognising the rich diversity of codes that have and continue to emerge, and the challenges of the coming age of petascale computing.

Lagrangian ocean analysis: Fundamentals and practices

DOE Office of Scientific and Technical Information (OSTI.GOV)

van Sebille, Erik; Griffies, Stephen M.; Abernathey, Ryan

Lagrangian analysis is a powerful way to analyse the output of ocean circulation models and other ocean velocity data such as from altimetry. In the Lagrangian approach, large sets of virtual particles are integrated within the three-dimensional, time-evolving velocity fields. A variety of tools and methods for this purpose have emerged, over several decades. Here, we review the state of the art in the field of Lagrangian analysis of ocean velocity data, starting from a fundamental kinematic framework and with a focus on large-scale open ocean applications. Beyond the use of explicit velocity fields, we consider the influence of unresolvedmore » physics and dynamics on particle trajectories. We comprehensively list and discuss the tools currently available for tracking virtual particles. We then showcase some of the innovative applications of trajectory data, and conclude with some open questions and an outlook. Our overall goal of this review paper is to reconcile some of the different techniques and methods in Lagrangian ocean analysis, while recognising the rich diversity of codes that have and continue to emerge, and the challenges of the coming age of petascale computing.« less

Preliminary Design Considerations for Access and Operations in Earth-Moon L1/L2 Orbits

NASA Technical Reports Server (NTRS)

Folta, David C.; Pavlak, Thomas A.; Haapala, Amanda F.; Howell, Kathleen C.

2013-01-01

Within the context of manned spaceflight activities, Earth-Moon libration point orbits could support lunar surface operations and serve as staging areas for future missions to near-Earth asteroids and Mars. This investigation examines preliminary design considerations including Earth-Moon L1/L2 libration point orbit selection, transfers, and stationkeeping costs associated with maintaining a spacecraft in the vicinity of L1 or L2 for a specified duration. Existing tools in multi-body trajectory design, dynamical systems theory, and orbit maintenance are leveraged in this analysis to explore end-to-end concepts for manned missions to Earth-Moon libration points.

Lagrangian continuum dynamics in ALEGRA.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Wong, Michael K. W.; Love, Edward

Alegra is an ALE (Arbitrary Lagrangian-Eulerian) multi-material finite element code that emphasizes large deformations and strong shock physics. The Lagrangian continuum dynamics package in Alegra uses a Galerkin finite element spatial discretization and an explicit central-difference stepping method in time. The goal of this report is to describe in detail the characteristics of this algorithm, including the conservation and stability properties. The details provided should help both researchers and analysts understand the underlying theory and numerical implementation of the Alegra continuum hydrodynamics algorithm.

Carbon Dioxide Compensation Points of Flowering Plants 1

PubMed Central

Krenzer, Eugene G.; Moss, Dale N.; Crookston, R. Kent

1975-01-01

Carbon dioxide compensation points of several hundred species of monocotyledons and dicotyledons have been measured during the course of various experiments in our laboratory over a period of several years. These have been classified into two groups: high, compensation points of 40 μl/l or greater; and low, compensation points of 10 μl/l or less. They are listed alphabetically both by families and species for monocotyledons and dicotyledons. Only two species did not unequivocally fit into the above established groups. These were Moricandia arvensis (L.) DC., which had an average compensation point of 26 μl/l and Panicum milioides Nees ex Trin., which was variable, but most often equilibrated between 12 to 20 μl/l CO2. PMID:16659272

Predictability of the Lagrangian Motion in the Upper Ocean

NASA Astrophysics Data System (ADS)

Piterbarg, L. I.; Griffa, A.; Griffa, A.; Mariano, A. J.; Ozgokmen, T. M.; Ryan, E. H.

2001-12-01

The complex non-linear dynamics of the upper ocean leads to chaotic behavior of drifter trajectories in the ocean. Our study is focused on estimating the predictability limit for the position of an individual Lagrangian particle or a particle cluster based on the knowledge of mean currents and observations of nearby particles (predictors). The Lagrangian prediction problem, besides being a fundamental scientific problem, is also of great importance for practical applications such as search and rescue operations and for modeling the spread of fish larvae. A stochastic multi-particle model for the Lagrangian motion has been rigorously formulated and is a generalization of the well known "random flight" model for a single particle. Our model is mathematically consistent and includes a few easily interpreted parameters, such as the Lagrangian velocity decorrelation time scale, the turbulent velocity variance, and the velocity decorrelation radius, that can be estimated from data. The top Lyapunov exponent for an isotropic version of the model is explicitly expressed as a function of these parameters enabling us to approximate the predictability limit to first order. Lagrangian prediction errors for two new prediction algorithms are evaluated against simple algorithms and each other and are used to test the predictability limits of the stochastic model for isotropic turbulence. The first algorithm is based on a Kalman filter and uses the developed stochastic model. Its implementation for drifter clusters in both the Tropical Pacific and Adriatic Sea, showed good prediction skill over a period of 1-2 weeks. The prediction error is primarily a function of the data density, defined as the number of predictors within a velocity decorrelation spatial scale from the particle to be predicted. The second algorithm is model independent and is based on spatial regression considerations. Preliminary results, based on simulated, as well as, real data, indicate that it performs

Lagrangian Timescales of Southern Ocean Upwelling in a Hierarchy of Model Resolutions

NASA Astrophysics Data System (ADS)

Drake, Henri F.; Morrison, Adele K.; Griffies, Stephen M.; Sarmiento, Jorge L.; Weijer, Wilbert; Gray, Alison R.

2018-01-01

In this paper we study upwelling pathways and timescales of Circumpolar Deep Water (CDW) in a hierarchy of models using a Lagrangian particle tracking method. Lagrangian timescales of CDW upwelling decrease from 87 years to 31 years to 17 years as the ocean resolution is refined from 1° to 0.25° to 0.1°. We attribute some of the differences in timescale to the strength of the eddy fields, as demonstrated by temporally degrading high-resolution model velocity fields. Consistent with the timescale dependence, we find that an average Lagrangian particle completes 3.2 circumpolar loops in the 1° model in comparison to 0.9 loops in the 0.1° model. These differences suggest that advective timescales and thus interbasin merging of upwelling CDW may be overestimated by coarse-resolution models, potentially affecting the skill of centennial scale climate change projections.

Analytical solution of the problem of a shock wave in the collapsing gas in Lagrangian coordinates

NASA Astrophysics Data System (ADS)

Kuropatenko, V. F.; Shestakovskaya, E. S.

2016-10-01

It is proposed the exact solution of the problem of a convergent shock wave and gas dynamic compression in a spherical vessel with an impermeable wall in Lagrangian coordinates. At the initial time the speed of cold ideal gas is equal to zero, and a negative velocity is set on boundary of the sphere. When t > t0 the shock wave spreads from this point into the gas. The boundary of the sphere will move under the certain law correlated with the motion of the shock wave. The trajectories of the gas particles in Lagrangian coordinates are straight lines. The equations determining the structure of the gas flow between the shock front and gas border have been found as a function of time and Lagrangian coordinate. The dependence of the entropy on the velocity of the shock wave has been found too. For Lagrangian coordinates the problem is first solved. It is fundamentally different from previously known formulations of the problem of the self-convergence of the self-similar shock wave to the center of symmetry and its reflection from the center, which was built up for the infinite area in Euler coordinates.

Two-dimensional Lagrangian simulation of suspended sediment

USGS Publications Warehouse

Schoellhamer, David H.

1988-01-01

A two-dimensional laterally averaged model for suspended sediment transport in steady gradually varied flow that is based on the Lagrangian reference frame is presented. The layered Lagrangian transport model (LLTM) for suspended sediment performs laterally averaged concentration. The elevations of nearly horizontal streamlines and the simulation time step are selected to optimize model stability and efficiency. The computational elements are parcels of water that are moved along the streamlines in the Lagrangian sense and are mixed with neighboring parcels. Three applications show that the LLTM can accurately simulate theoretical and empirical nonequilibrium suspended sediment distributions and slug injections of suspended sediment in a laboratory flume.

78 FR 32294 - DeltaPoint Capital IV, L.P., DeltaPoint Capital IV (New York), L.P., License No. 02/02-0662,02/02...

Federal Register 2010, 2011, 2012, 2013, 2014

2013-05-29

... Small Business Investment Act of 1958, as amended (``the Act''), in connection with the financing of a... SMALL BUSINESS ADMINISTRATION DeltaPoint Capital IV, L.P., DeltaPoint Capital IV (New York), L.P., License No. 02/02-0662,02/02-0661; Notice Seeking Exemption Under Section 312 of the Small Business...

An online-coupled NWP/ACT model with conserved Lagrangian levels

NASA Astrophysics Data System (ADS)

Sørensen, B.; Kaas, E.; Lauritzen, P. H.

2012-04-01

Numerical weather and climate modelling is under constant development. Semi-implicit semi-Lagrangian (SISL) models have proven to be numerically efficient in both short-range weather forecasts and climate models, due to the ability to use long time steps. Chemical/aerosol feedback mechanism are becoming more and more relevant in NWP as well as climate models, since the biogenic and anthropogenic emissions can have a direct effect on the dynamics and radiative properties of the atmosphere. To include chemical feedback mechanisms in the NWP models, on-line coupling is crucial. In 3D semi-Lagrangian schemes with quasi-Lagrangian vertical coordinates the Lagrangian levels are remapped to Eulerian model levels each time step. This remapping introduces an undesirable tendency to smooth sharp gradients and creates unphysical numerical diffusion in the vertical distribution. A semi-Lagrangian advection method is introduced, it combines an inherently mass conserving 2D semi-Lagrangian scheme, with a SISL scheme employing both hybrid vertical coordinates and a fully Lagrangian vertical coordinate. This minimizes the vertical diffusion and thus potentially improves the simulation of the vertical profiles of moisture, clouds, and chemical constituents. Since the Lagrangian levels suffer from traditional Lagrangian limitations caused by the convergence and divergence of the flow, remappings to the Eulerian model levels are generally still required - but this need only be applied after a number of time steps - unless dynamic remapping methods are used. For this several different remapping methods has been implemented. The combined scheme is mass conserving, consistent, and multi-tracer efficient.

Spontaneous T-cell responses against the immune check point programmed-death-ligand 1 (PD-L1) in patients with chronic myeloproliferative neoplasms correlate with disease stage and clinical response.

PubMed

Holmström, Morten Orebo; Riley, Caroline Hasselbalch; Skov, Vibe; Svane, Inge Marie; Hasselbalch, Hans Carl; Andersen, Mads Hald

2018-01-01

The Chronic Myeloproliferative Neoplasms (MPN) are cancers characterized by hyperinflammation and immune deregulation. Concurrently, the expression of the immune check point programmed death ligand 1 (PD-L1) is induced by inflammation. In this study we report on the occurrence of spontaneous T cell responses against a PD-L1 derived epitope in patients with MPN. We show that 71% of patients display a significant immune response against PD-L1, and patients with advanced MPN have significantly fewer and weaker PD-L1 specific immune responses compared to patients with non-advanced MPN. The PD-L1 specific T cell responses are CD4 + T cell responses, and by gene expression analysis we show that expression of PD-L1 is enhanced in patients with MPN. This could imply that the tumor specific immune response in MPN could be enhanced by vaccination with PD-L1 derived epitopes by boosting the anti-regulatory immune response hereby allowing tumor specific T cell to exert anti-tumor immunity.

Lagrangian modeling of global atmospheric methane (1990-2012)

NASA Astrophysics Data System (ADS)

Arfeuille, Florian; Henne, Stephan; Brunner, Dominik

2016-04-01

In the MAIOLICA-II project, the lagrangian particle model FLEXPART is used to simulate the global atmospheric methane over the 1990-2012 period. In this lagrangian framework, 3 million particles are permanently transported based on winds from ERA-interim. The history of individual particles can be followed allowing for a comprehensive analysis of transport pathways and timescales. The link between sources (emissions) and receptors (measurement stations) is then established in a straightforward manner, a prerequisite for source inversion problems. FLEXPART was extended to incorporate the methane loss by reaction with OH, soil uptake and stratospheric loss reactions with prescribed Cl and O(1d) radicals. Sources are separated into 245 different tracers, depending on source origin (anthropogenic, wetlands, rice, biomass burning, termites, wild animals, oceans, volcanoes), region of emission, and time since emission (5 age classes). The inversion method applied is a fixed-lag Kalman smoother similar to that described in Bruhwiler et al. [2005]. Results from the FLEXPART global methane simulation and from the subsequent inversion will be presented. Results notably suggest: - A reduction in methane growth rates due to diminished wetland emissions and anthropogenic European emission in 1990-1993. - A second decrease in 1995-1996 is also mainly attributed to these two emission categories. - A reduced increase in Chinese anthropogenic emissions after 2003 compared to EDGAR inventories. - Large South American wetlands emissions during the entire period. Bruhwiler, L. M. P., Michalak, A. M., Peters, W., Baker, D. F. & Tans, P. 2005: An improved Kalman smoother fore atmospheric inversions, Atmos Chem Phys, 5, 2691-2702.

Lagrangian chaos in three- dimensional steady buoyancy-driven flows

NASA Astrophysics Data System (ADS)

Contreras, Sebastian; Speetjens, Michel; Clercx, Herman

2016-11-01

Natural convection plays a key role in fluid dynamics owing to its ubiquitous presence in nature and industry. Buoyancy-driven flows are prototypical systems in the study of thermal instabilities and pattern formation. The differentially heated cavity problem has been widely studied for the investigation of buoyancy-induced oscillatory flow. However, far less attention has been devoted to the three-dimensional Lagrangian transport properties in such flows. This study seeks to address this by investigating Lagrangian transport in the steady flow inside a cubic cavity differentially-heated from the side. The theoretical and numerical analysis expands on previously reported similarities between the current flow and lid-driven flows. The Lagrangian dynamics are controlled by the Péclet number (Pe) and the Prandtl number (Pr). Pe controls the behaviour qualitatively in that growing Pe progressively perturbs the integable state (Pe =0), thus paving the way to chaotic dynamics. Pr plays an entirely quantitative role in that Pr<1 and Pr>1 amplifies and diminishes, respectively, the perturbative effect of non-zero Pe. S.C. acknowledges financial support from Consejo Nacional de Ciencia y Tecnología (CONACYT).

Bayesian Lagrangian Data Assimilation and Drifter Deployment Strategies

NASA Astrophysics Data System (ADS)

Dutt, A.; Lermusiaux, P. F. J.

2017-12-01

Ocean currents transport a variety of natural (e.g. water masses, phytoplankton, zooplankton, sediments, etc.) and man-made materials and other objects (e.g. pollutants, floating debris, search and rescue, etc.). Lagrangian Coherent Structures (LCSs) or the most influential/persistent material lines in a flow, provide a robust approach to characterize such Lagrangian transports and organize classic trajectories. Using the flow-map stochastic advection and a dynamically-orthogonal decomposition, we develop uncertainty prediction schemes for both Eulerian and Lagrangian variables. We then extend our Bayesian Gaussian Mixture Model (GMM)-DO filter to a joint Eulerian-Lagrangian Bayesian data assimilation scheme. The resulting nonlinear filter allows the simultaneous non-Gaussian estimation of Eulerian variables (e.g. velocity, temperature, salinity, etc.) and Lagrangian variables (e.g. drifter/float positions, trajectories, LCSs, etc.). Its results are showcased using a double-gyre flow with a random frequency, a stochastic flow past a cylinder, and realistic ocean examples. We further show how our Bayesian mutual information and adaptive sampling equations provide a rigorous efficient methodology to plan optimal drifter deployment strategies and predict the optimal times, locations, and types of measurements to be collected.

Scale-by-scale contributions to Lagrangian particle acceleration

NASA Astrophysics Data System (ADS)

Lalescu, Cristian C.; Wilczek, Michael

2017-11-01

Fluctuations on a wide range of scales in both space and time are characteristic of turbulence. Lagrangian particles, advected by the flow, probe these fluctuations along their trajectories. In an effort to isolate the influence of the different scales on Lagrangian statistics, we employ direct numerical simulations (DNS) combined with a filtering approach. Specifically, we study the acceleration statistics of tracers advected in filtered fields to characterize the smallest temporal scales of the flow. Emphasis is put on the acceleration variance as a function of filter scale, along with the scaling properties of the relevant terms of the Navier-Stokes equations. We furthermore discuss scaling ranges for higher-order moments of the tracer acceleration, as well as the influence of the choice of filter on the results. Starting from the Lagrangian tracer acceleration as the short time limit of the Lagrangian velocity increment, we also quantify the influence of filtering on Lagrangian intermittency. Our work complements existing experimental results on intermittency and accelerations of finite-sized, neutrally-buoyant particles: for the passive tracers used in our DNS, feedback effects are neglected such that the spatial averaging effect is cleanly isolated.

Chaotic Lagrangian models for turbulent relative dispersion.

PubMed

Lacorata, Guglielmo; Vulpiani, Angelo

2017-04-01

A deterministic multiscale dynamical system is introduced and discussed as a prototype model for relative dispersion in stationary, homogeneous, and isotropic turbulence. Unlike stochastic diffusion models, here trajectory transport and mixing properties are entirely controlled by Lagrangian chaos. The anomalous "sweeping effect," a known drawback common to kinematic simulations, is removed through the use of quasi-Lagrangian coordinates. Lagrangian dispersion statistics of the model are accurately analyzed by computing the finite-scale Lyapunov exponent (FSLE), which is the optimal measure of the scaling properties of dispersion. FSLE scaling exponents provide a severe test to decide whether model simulations are in agreement with theoretical expectations and/or observation. The results of our numerical experiments cover a wide range of "Reynolds numbers" and show that chaotic deterministic flows can be very efficient, and numerically low-cost, models of turbulent trajectories in stationary, homogeneous, and isotropic conditions. The mathematics of the model is relatively simple, and, in a geophysical context, potential applications may regard small-scale parametrization issues in general circulation models, mixed layer, and/or boundary layer turbulence models as well as Lagrangian predictability studies.

Chaotic Lagrangian models for turbulent relative dispersion

NASA Astrophysics Data System (ADS)

Lacorata, Guglielmo; Vulpiani, Angelo

2017-04-01

A deterministic multiscale dynamical system is introduced and discussed as a prototype model for relative dispersion in stationary, homogeneous, and isotropic turbulence. Unlike stochastic diffusion models, here trajectory transport and mixing properties are entirely controlled by Lagrangian chaos. The anomalous "sweeping effect," a known drawback common to kinematic simulations, is removed through the use of quasi-Lagrangian coordinates. Lagrangian dispersion statistics of the model are accurately analyzed by computing the finite-scale Lyapunov exponent (FSLE), which is the optimal measure of the scaling properties of dispersion. FSLE scaling exponents provide a severe test to decide whether model simulations are in agreement with theoretical expectations and/or observation. The results of our numerical experiments cover a wide range of "Reynolds numbers" and show that chaotic deterministic flows can be very efficient, and numerically low-cost, models of turbulent trajectories in stationary, homogeneous, and isotropic conditions. The mathematics of the model is relatively simple, and, in a geophysical context, potential applications may regard small-scale parametrization issues in general circulation models, mixed layer, and/or boundary layer turbulence models as well as Lagrangian predictability studies.

Lagrangian theory of structure formation in relativistic cosmology. IV. Lagrangian approach to gravitational waves

NASA Astrophysics Data System (ADS)

Al Roumi, Fosca; Buchert, Thomas; Wiegand, Alexander

2017-12-01

The relativistic generalization of the Newtonian Lagrangian perturbation theory is investigated. In previous works, the perturbation and solution schemes that are generated by the spatially projected gravitoelectric part of the Weyl tensor were given to any order of the perturbations, together with extensions and applications for accessing the nonperturbative regime. We here discuss more in detail the general first-order scheme within the Cartan formalism including and concentrating on the gravitational wave propagation in matter. We provide master equations for all parts of Lagrangian-linearized perturbations propagating in the perturbed spacetime, and we outline the solution procedure that allows one to find general solutions. Particular emphasis is given to global properties of the Lagrangian perturbation fields by employing results of Hodge-de Rham theory. We here discuss how the Hodge decomposition relates to the standard scalar-vector-tensor decomposition. Finally, we demonstrate that we obtain the known linear perturbation solutions of the standard relativistic perturbation scheme by performing two steps: first, by restricting our solutions to perturbations that propagate on a flat unperturbed background spacetime and, second, by transforming to Eulerian background coordinates with truncation of nonlinear terms.

Investigating Trojan Asteroids at the L4/L5 Sun-Earth Lagrange Points

NASA Technical Reports Server (NTRS)

John, K. K.; Graham, L. D.; Abell, P. A.

2015-01-01

Investigations of Earth's Trojan asteroids will have benefits for science, exploration, and resource utilization. By sending a small spacecraft to the Sun-Earth L4 or L5 Lagrange points to investigate near-Earth objects, Earth's Trojan population can be better understood. This could lead to future missions for larger precursor spacecraft as well as human missions. The presence of objects in the Sun-Earth L4 and L5 Lagrange points has long been suspected, and in 2010 NASA's Wide-field Infrared Survey Explorer (WISE) detected a 300 m object. To investigate these Earth Trojan asteroid objects, it is both essential and feasible to send spacecraft to these regions. By exploring a wide field area, a small spacecraft equipped with an IR camera could hunt for Trojan asteroids and other Earth co-orbiting objects at the L4 or L5 Lagrange points in the near-term. By surveying the region, a zeroth-order approximation of the number of objects could be obtained with some rough constraints on their diameters, which may lead to the identification of potential candidates for further study. This would serve as a precursor for additional future robotic and human exploration targets. Depending on the inclination of these potential objects, they could be used as proving areas for future missions in the sense that the delta-V's to get to these targets are relatively low as compared to other rendezvous missions. They can serve as platforms for extended operations in deep space while interacting with a natural object in microgravity. Theoretically, such low inclination Earth Trojan asteroids exist. By sending a spacecraft to L4 or L5, these likely and potentially accessible targets could be identified.

Imposing a Lagrangian Particle Framework on an Eulerian Hydrodynamics Infrastructure in Flash

NASA Technical Reports Server (NTRS)

Dubey, A.; Daley, C.; ZuHone, J.; Ricker, P. M.; Weide, K.; Graziani, C.

2012-01-01

In many astrophysical simulations, both Eulerian and Lagrangian quantities are of interest. For example, in a galaxy cluster merger simulation, the intracluster gas can have Eulerian discretization, while dark matter can be modeled using particles. FLASH, a component-based scientific simulation code, superimposes a Lagrangian framework atop an adaptive mesh refinement Eulerian framework to enable such simulations. The discretization of the field variables is Eulerian, while the Lagrangian entities occur in many different forms including tracer particles, massive particles, charged particles in particle-in-cell mode, and Lagrangian markers to model fluid structure interactions. These widely varying roles for Lagrangian entities are possible because of the highly modular, flexible, and extensible architecture of the Lagrangian framework. In this paper, we describe the Lagrangian framework in FLASH in the context of two very different applications, Type Ia supernovae and galaxy cluster mergers, which use the Lagrangian entities in fundamentally different ways.

Imposing a Lagrangian Particle Framework on an Eulerian Hydrodynamics Infrastructure in FLASH

NASA Astrophysics Data System (ADS)

Dubey, A.; Daley, C.; ZuHone, J.; Ricker, P. M.; Weide, K.; Graziani, C.

2012-08-01

In many astrophysical simulations, both Eulerian and Lagrangian quantities are of interest. For example, in a galaxy cluster merger simulation, the intracluster gas can have Eulerian discretization, while dark matter can be modeled using particles. FLASH, a component-based scientific simulation code, superimposes a Lagrangian framework atop an adaptive mesh refinement Eulerian framework to enable such simulations. The discretization of the field variables is Eulerian, while the Lagrangian entities occur in many different forms including tracer particles, massive particles, charged particles in particle-in-cell mode, and Lagrangian markers to model fluid-structure interactions. These widely varying roles for Lagrangian entities are possible because of the highly modular, flexible, and extensible architecture of the Lagrangian framework. In this paper, we describe the Lagrangian framework in FLASH in the context of two very different applications, Type Ia supernovae and galaxy cluster mergers, which use the Lagrangian entities in fundamentally different ways.

Infrared and visible image fusion based on total variation and augmented Lagrangian.

PubMed

Guo, Hanqi; Ma, Yong; Mei, Xiaoguang; Ma, Jiayi

2017-11-01

This paper proposes a new algorithm for infrared and visible image fusion based on gradient transfer that achieves fusion by preserving the intensity of the infrared image and then transferring gradients in the corresponding visible one to the result. The gradient transfer suffers from the problems of low dynamic range and detail loss because it ignores the intensity from the visible image. The new algorithm solves these problems by providing additive intensity from the visible image to balance the intensity between the infrared image and the visible one. It formulates the fusion task as an l 1 -l 1 -TV minimization problem and then employs variable splitting and augmented Lagrangian to convert the unconstrained problem to a constrained one that can be solved in the framework of alternating the multiplier direction method. Experiments demonstrate that the new algorithm achieves better fusion results with a high computation efficiency in both qualitative and quantitative tests than gradient transfer and most state-of-the-art methods.

COLAcode: COmoving Lagrangian Acceleration code

NASA Astrophysics Data System (ADS)

Tassev, Svetlin V.

2016-02-01

COLAcode is a serial particle mesh-based N-body code illustrating the COLA (COmoving Lagrangian Acceleration) method; it solves for Large Scale Structure (LSS) in a frame that is comoving with observers following trajectories calculated in Lagrangian Perturbation Theory (LPT). It differs from standard N-body code by trading accuracy at small-scales to gain computational speed without sacrificing accuracy at large scales. This is useful for generating large ensembles of accurate mock halo catalogs required to study galaxy clustering and weak lensing; such catalogs are needed to perform detailed error analysis for ongoing and future surveys of LSS.

Source apportion of atmospheric particulate matter: a joint Eulerian/Lagrangian approach.

PubMed

Riccio, A; Chianese, E; Agrillo, G; Esposito, C; Ferrara, L; Tirimberio, G

2014-12-01

PM2.5 samples were collected during an annual monitoring campaign (January 2012-January 2013) in the urban area of Naples, one of the major cities in Southern Italy. Samples were collected by means of a standard gravimetric sampler (Tecora Echo model) and characterized from a chemical point of view by ion chromatography. As a result, 143 samples together with their ionic composition have been collected. We extend traditional source apportionment techniques, usually based on multivariate factor analysis, interpreting the chemical analysis results within a Lagrangian framework. The Hybrid Single-Particle Lagrangian Integrated Trajectory Model (HYSPLIT) model was used, providing linkages to the source regions in the upwind areas. Results were analyzed in order to quantify the relative weight of different source types/areas. Model results suggested that PM concentrations are strongly affected not only by local emissions but also by transboundary emissions, especially from the Eastern and Northern European countries and African Saharan dust episodes.

Lagrangian geometrical optics of nonadiabatic vector waves and spin particles

DOE PAGES

Ruiz, D. E.; Dodin, I. Y.

2015-07-29

Linear vector waves, both quantum and classical, experience polarization-driven bending of ray trajectories and polarization dynamics that can be interpreted as the precession of the "wave spin". Here, both phenomena are governed by an effective gauge Hamiltonian vanishing in leading-order geometrical optics. This gauge Hamiltonian can be recognized as a generalization of the Stern-Gerlach Hamiltonian that is commonly known for spin-1/2 quantum particles. The corresponding reduced Lagrangians for continuous nondissipative waves and their geometrical-optics rays are derived from the fundamental wave Lagrangian. The resulting Euler-Lagrange equations can describe simultaneous interactions of N resonant modes, where N is arbitrary, and leadmore » to equations for the wave spin, which happens to be an (N 2 - 1)-dimensional spin vector. As a special case, classical equations for a Dirac particle (N = 2) are deduced formally, without introducing additional postulates or interpretations, from the Dirac quantum Lagrangian with the Pauli term. The model reproduces the Bargmann-Michel-Telegdi equations with added Stern-Gerlach force.« less

Nonunitary Lagrangians and Unitary Non-Lagrangian Conformal Field Theories.

PubMed

Buican, Matthew; Laczko, Zoltan

2018-02-23

In various dimensions, we can sometimes compute observables of interacting conformal field theories (CFTs) that are connected to free theories via the renormalization group (RG) flow by computing protected quantities in the free theories. On the other hand, in two dimensions, it is often possible to algebraically construct observables of interacting CFTs using free fields without the need to explicitly construct an underlying RG flow. In this Letter, we begin to extend this idea to higher dimensions by showing that one can compute certain observables of an infinite set of unitary strongly interacting four-dimensional N=2 superconformal field theories (SCFTs) by performing simple calculations involving sets of nonunitary free four-dimensional hypermultiplets. These free fields are distant cousins of the Majorana fermion underlying the two-dimensional Ising model and are not obviously connected to our interacting theories via an RG flow. Rather surprisingly, this construction gives us Lagrangians for particular observables in certain subsectors of many "non-Lagrangian" SCFTs by sacrificing unitarity while preserving the full N=2 superconformal algebra. As a by-product, we find relations between characters in unitary and nonunitary affine Kac-Moody algebras. We conclude by commenting on possible generalizations of our construction.

Nonunitary Lagrangians and Unitary Non-Lagrangian Conformal Field Theories

NASA Astrophysics Data System (ADS)

Buican, Matthew; Laczko, Zoltan

2018-02-01

In various dimensions, we can sometimes compute observables of interacting conformal field theories (CFTs) that are connected to free theories via the renormalization group (RG) flow by computing protected quantities in the free theories. On the other hand, in two dimensions, it is often possible to algebraically construct observables of interacting CFTs using free fields without the need to explicitly construct an underlying RG flow. In this Letter, we begin to extend this idea to higher dimensions by showing that one can compute certain observables of an infinite set of unitary strongly interacting four-dimensional N =2 superconformal field theories (SCFTs) by performing simple calculations involving sets of nonunitary free four-dimensional hypermultiplets. These free fields are distant cousins of the Majorana fermion underlying the two-dimensional Ising model and are not obviously connected to our interacting theories via an RG flow. Rather surprisingly, this construction gives us Lagrangians for particular observables in certain subsectors of many "non-Lagrangian" SCFTs by sacrificing unitarity while preserving the full N =2 superconformal algebra. As a by-product, we find relations between characters in unitary and nonunitary affine Kac-Moody algebras. We conclude by commenting on possible generalizations of our construction.

An adaptive reconstruction for Lagrangian, direct-forcing, immersed-boundary methods

NASA Astrophysics Data System (ADS)

Posa, Antonio; Vanella, Marcos; Balaras, Elias

2017-12-01

Lagrangian, direct-forcing, immersed boundary (IB) methods have been receiving increased attention due to their robustness in complex fluid-structure interaction problems. They are very sensitive, however, on the selection of the Lagrangian grid, which is typically used to define a solid or flexible body immersed in a fluid flow. In the present work we propose a cost-efficient solution to this problem without compromising accuracy. Central to our approach is the use of isoparametric mapping to bridge the relative resolution requirements of Lagrangian IB, and Eulerian grids. With this approach, the density of surface Lagrangian markers, which is essential to properly enforce boundary conditions, is adapted dynamically based on the characteristics of the underlying Eulerian grid. The markers are not stored and the Lagrangian data-structure is not modified. The proposed scheme is implemented in the framework of a moving least squares reconstruction formulation, but it can be adapted to any Lagrangian, direct-forcing formulation. The accuracy and robustness of the approach is demonstrated in a variety of test cases of increasing complexity.

Interaction of the stream from L 1 with the outer edge of the accretion disk in a cataclysmic variable

NASA Astrophysics Data System (ADS)

Kaigorodov, P. V.; Bisikalo, D. V.; Kurbatov, E. P.

2017-08-01

Vertical oscillations of the gas at the outer edge of the accretion disk in a semi-detached binary due to interaction with the stream of matter from the inner Lagrangian point L 1 are considered. Mixing of the matter from the stream from L 1 with matter of the disk halo results in the formation of a system of two diverging shocks and a contact discontinuity, or so-called "hot line". The passage of matter through the region of the hot line leads to an increase in its vertical velocity and a thickening of the disk at phases 0.7-0.8. Subsequently, the matter moving along the outer edge of the disk also experiences vertical oscillations, forming secondary maxima at phases 0.2-0.4. It is shown that, for systems with component mass ratios of 0.6, these oscillations will be amplified with each passage of the matter through the hotline zone, while the observations will be quenched in systems with component mass ratios 0.07 and 7. The most favorable conditions for the flow of matter from the stream through the edge of the disk arise for component mass ratios 0.62. A theoretical relation between the phases of disk thickenings and the component mass ratio of the system is derived.

Scalar curvature of Lagrangian Riemannian submersions and their harmonicity

NASA Astrophysics Data System (ADS)

Eken Meri˙ç, Şemsi; Kiliç, Erol; Sağiroğlu, Yasemi˙n

In this paper, we consider a Lagrangian Riemannian submersion from a Hermitian manifold to a Riemannian manifold and establish some basic inequalities to obtain relationships between the intrinsic and extrinsic invariants for such a submersion. Indeed, using these inequalities, we provide necessary and sufficient conditions for which a Lagrangian Riemannian submersion π has totally geodesic or totally umbilical fibers. Moreover, we study the harmonicity of Lagrangian Riemannian submersions and obtain a characterization for such submersions to be harmonic.

A Lagrangian effective field theory

DOE Office of Scientific and Technical Information (OSTI.GOV)

Vlah, Zvonimir; White, Martin; Aviles, Alejandro

We have continued the development of Lagrangian, cosmological perturbation theory for the low-order correlators of the matter density field. We provide a new route to understanding how the effective field theory (EFT) of large-scale structure can be formulated in the Lagrandian framework and a new resummation scheme, comparing our results to earlier work and to a series of high-resolution N-body simulations in both Fourier and configuration space. The `new' terms arising from EFT serve to tame the dependence of perturbation theory on small-scale physics and improve agreement with simulations (though with an additional free parameter). We find that all ofmore » our models fare well on scales larger than about two to three times the non-linear scale, but fail as the non-linear scale is approached. This is slightly less reach than has been seen previously. At low redshift the Lagrangian model fares as well as EFT in its Eulerian formulation, but at higher z the Eulerian EFT fits the data to smaller scales than resummed, Lagrangian EFT. Furthermore, all the perturbative models fare better than linear theory.« less

A Lagrangian effective field theory

DOE Office of Scientific and Technical Information (OSTI.GOV)

Vlah, Zvonimir; White, Martin; Aviles, Alejandro, E-mail: zvlah@stanford.edu, E-mail: mwhite@berkeley.edu, E-mail: aviles@berkeley.edu

We have continued the development of Lagrangian, cosmological perturbation theory for the low-order correlators of the matter density field. We provide a new route to understanding how the effective field theory (EFT) of large-scale structure can be formulated in the Lagrandian framework and a new resummation scheme, comparing our results to earlier work and to a series of high-resolution N-body simulations in both Fourier and configuration space. The 'new' terms arising from EFT serve to tame the dependence of perturbation theory on small-scale physics and improve agreement with simulations (though with an additional free parameter). We find that all ofmore » our models fare well on scales larger than about two to three times the non-linear scale, but fail as the non-linear scale is approached. This is slightly less reach than has been seen previously. At low redshift the Lagrangian model fares as well as EFT in its Eulerian formulation, but at higher z the Eulerian EFT fits the data to smaller scales than resummed, Lagrangian EFT. All the perturbative models fare better than linear theory.« less

A Lagrangian effective field theory

DOE PAGES

Vlah, Zvonimir; White, Martin; Aviles, Alejandro

2015-09-02

We have continued the development of Lagrangian, cosmological perturbation theory for the low-order correlators of the matter density field. We provide a new route to understanding how the effective field theory (EFT) of large-scale structure can be formulated in the Lagrandian framework and a new resummation scheme, comparing our results to earlier work and to a series of high-resolution N-body simulations in both Fourier and configuration space. The `new' terms arising from EFT serve to tame the dependence of perturbation theory on small-scale physics and improve agreement with simulations (though with an additional free parameter). We find that all ofmore » our models fare well on scales larger than about two to three times the non-linear scale, but fail as the non-linear scale is approached. This is slightly less reach than has been seen previously. At low redshift the Lagrangian model fares as well as EFT in its Eulerian formulation, but at higher z the Eulerian EFT fits the data to smaller scales than resummed, Lagrangian EFT. Furthermore, all the perturbative models fare better than linear theory.« less

Early Mission Maneuver Operations for the Deep Space Climate Observatory Sun-Earth L1 Libration Point Mission

NASA Technical Reports Server (NTRS)

Roberts, Craig; Case, Sara; Reagoso, John; Webster, Cassandra

2015-01-01

The Deep Space Climate Observatory mission launched on February 11, 2015, and inserted onto a transfer trajectory toward a Lissajous orbit around the Sun-Earth L1 libration point. This paper presents an overview of the baseline transfer orbit and early mission maneuver operations leading up to the start of nominal science orbit operations. In particular, the analysis and performance of the spacecraft insertion, mid-course correction maneuvers, and the deep-space Lissajous orbit insertion maneuvers are discussed, com-paring the baseline orbit with actual mission results and highlighting mission and operations constraints..

Lagrangian methods of cosmic web classification

NASA Astrophysics Data System (ADS)

Fisher, J. D.; Faltenbacher, A.; Johnson, M. S. T.

2016-05-01

The cosmic web defines the large-scale distribution of matter we see in the Universe today. Classifying the cosmic web into voids, sheets, filaments and nodes allows one to explore structure formation and the role environmental factors have on halo and galaxy properties. While existing studies of cosmic web classification concentrate on grid-based methods, this work explores a Lagrangian approach where the V-web algorithm proposed by Hoffman et al. is implemented with techniques borrowed from smoothed particle hydrodynamics. The Lagrangian approach allows one to classify individual objects (e.g. particles or haloes) based on properties of their nearest neighbours in an adaptive manner. It can be applied directly to a halo sample which dramatically reduces computational cost and potentially allows an application of this classification scheme to observed galaxy samples. Finally, the Lagrangian nature admits a straightforward inclusion of the Hubble flow negating the necessity of a visually defined threshold value which is commonly employed by grid-based classification methods.

Evaluation of wastewater contaminant transport in surface waters using verified Lagrangian sampling.

PubMed

Antweiler, Ronald C; Writer, Jeffrey H; Murphy, Sheila F

2014-02-01

Contaminants released from wastewater treatment plants can persist in surface waters for substantial distances. Much research has gone into evaluating the fate and transport of these contaminants, but this work has often assumed constant flow from wastewater treatment plants. However, effluent discharge commonly varies widely over a 24-hour period, and this variation controls contaminant loading and can profoundly influence interpretations of environmental data. We show that methodologies relying on the normalization of downstream data to conservative elements can give spurious results, and should not be used unless it can be verified that the same parcel of water was sampled. Lagrangian sampling, which in theory samples the same water parcel as it moves downstream (the Lagrangian parcel), links hydrologic and chemical transformation processes so that the in-stream fate of wastewater contaminants can be quantitatively evaluated. However, precise Lagrangian sampling is difficult, and small deviations - such as missing the Lagrangian parcel by less than 1h - can cause large differences in measured concentrations of all dissolved compounds at downstream sites, leading to erroneous conclusions regarding in-stream processes controlling the fate and transport of wastewater contaminants. Therefore, we have developed a method termed "verified Lagrangian" sampling, which can be used to determine if the Lagrangian parcel was actually sampled, and if it was not, a means for correcting the data to reflect the concentrations which would have been obtained had the Lagrangian parcel been sampled. To apply the method, it is necessary to have concentration data for a number of conservative constituents from the upstream, effluent, and downstream sites, along with upstream and effluent concentrations that are constant over the short-term (typically 2-4h). These corrections can subsequently be applied to all data, including non-conservative constituents. Finally, we show how data

Bayesian Nonlinear Assimilation of Eulerian and Lagrangian Coastal Flow Data

DTIC Science & Technology

2015-09-30

Lagrangian Coastal Flow Data Dr. Pierre F.J. Lermusiaux Department of Mechanical Engineering Center for Ocean Science and Engineering Massachusetts...Develop and apply theory, schemes and computational systems for rigorous Bayesian nonlinear assimilation of Eulerian and Lagrangian coastal flow data...coastal ocean fields, both in Eulerian and Lagrangian forms. - Further develop and implement our GMM-DO schemes for robust Bayesian nonlinear estimation

Euler-Lagrangian computation for estuarine hydrodynamics

USGS Publications Warehouse

Cheng, Ralph T.

1983-01-01

The transport of conservative and suspended matter in fluid flows is a phenomenon of Lagrangian nature because the process is usually convection dominant. Nearly all numerical investigations of such problems use an Eulerian formulation for the convenience that the computational grids are fixed in space and because the vast majority of field data are collected in an Eulerian reference frame. Several examples are given in this paper to illustrate a modeling approach which combines the advantages of both the Eulerian and Lagrangian computational techniques.

On the lagrangian 1-form structure of the hyperbolic calogero-moser system

NASA Astrophysics Data System (ADS)

Jairuk, Umpon; Tanasittikosol, Monsit; Yoo-Kong, Sikarin

2017-06-01

In this work, we present the Lagrangian 1-form structure of the hyperbolic Calogero-Moser system in both discrete-time level and continuous-time level. The discrete-time hyperbolic Calogero-Moser system is obtained by considering pole reduction of the semi-discrete Kadomtsev-Petviashvili (KP) equation. Furthermore, it is shown that the hyperbolic Calogero-Moser system possesses the key relation, known as the discrete-time closure relation. This relation is a consequence of the compatibility property of the temporal Lax matrices. The continuous-time hierarchy of the hyperbolic Calogero-Moser system is obtained by taking two successive continuum limits, namely, the skewed limit and full limit. With these successive limits, the continuous-time closure relation is also obtained and is shown to hold at the continuous level.

Sea Fog Forecasting with Lagrangian Models

NASA Astrophysics Data System (ADS)

Lewis, J. M.

2014-12-01

In 1913, G. I. Taylor introduced us to a Lagrangian view of sea fog formation. He conducted his study off the coast of Newfoundland in the aftermath of the Titanic disaster. We briefly review Taylor's classic work and then apply these same principles to a case of sea fog formation and dissipation off the coast of California. The resources used in this study consist of: 1) land-based surface and upper-air observations, 2) NDBC (National Data Buoy Center) observations from moored buoys equipped to measure dew point temperature as well as the standard surface observations at sea (wind, sea surface temperature, pressure, and air temperature), 3) satellite observations of cloud, and 4) a one-dimensional (vertically directed) boundary layer model that tracks with the surface air motion and makes use of sophisticated turbulence-radiation parameterizations. Results of the investigation indicate that delicate interplay and interaction between the radiation and turbulence processes makes accurate forecasts of sea fog onset unlikely in the near future. This pessimistic attitude stems from inadequacy of the existing network of observations and uncertainties in modeling dynamical processes within the boundary layer.

An Eulerian/Lagrangian coupling procedure for three-dimensional vortical flows

NASA Technical Reports Server (NTRS)

Felici, Helene M.; Drela, Mark

1993-01-01

A coupled Eulerian/Lagrangian method is presented for the reduction of numerical diffusion observed in solutions of 3D vortical flows using standard Eulerian finite-volume time-marching procedures. A Lagrangian particle tracking method, added to the Eulerian time-marching procedure, provides a correction of the Eulerian solution. In turn, the Eulerian solution is used to integrate the Lagrangian state-vector along the particles trajectories. While the Eulerian solution ensures the conservation of mass and sets the pressure field, the particle markers describe accurately the convection properties and enhance the vorticity and entropy capturing capabilities of the Eulerian solver. The Eulerian/Lagrangian coupling strategies are discussed and the combined scheme is tested on a constant stagnation pressure flow in a 90 deg bend and on a swirling pipe flow. As the numerical diffusion is reduced when using the Lagrangian correction, a vorticity gradient augmentation is identified as a basic problem of this inviscid calculation.

Thermostating extended Lagrangian Born-Oppenheimer molecular dynamics.

PubMed

Martínez, Enrique; Cawkwell, Marc J; Voter, Arthur F; Niklasson, Anders M N

2015-04-21

Extended Lagrangian Born-Oppenheimer molecular dynamics is developed and analyzed for applications in canonical (NVT) simulations. Three different approaches are considered: the Nosé and Andersen thermostats and Langevin dynamics. We have tested the temperature distribution under different conditions of self-consistent field (SCF) convergence and time step and compared the results to analytical predictions. We find that the simulations based on the extended Lagrangian Born-Oppenheimer framework provide accurate canonical distributions even under approximate SCF convergence, often requiring only a single diagonalization per time step, whereas regular Born-Oppenheimer formulations exhibit unphysical fluctuations unless a sufficiently high degree of convergence is reached at each time step. The thermostated extended Lagrangian framework thus offers an accurate approach to sample processes in the canonical ensemble at a fraction of the computational cost of regular Born-Oppenheimer molecular dynamics simulations.

Adaptive reconnection-based arbitrary Lagrangian Eulerian method

DOE Office of Scientific and Technical Information (OSTI.GOV)

Bo, Wurigen; Shashkov, Mikhail

We present a new adaptive Arbitrary Lagrangian Eulerian (ALE) method. This method is based on the reconnection-based ALE (ReALE) methodology of Refs. [35], [34] and [6]. The main elements in a standard ReALE method are: an explicit Lagrangian phase on an arbitrary polygonal (in 2D) mesh in which the solution and positions of grid nodes are updated; a rezoning phase in which a new grid is defined by changing the connectivity (using Voronoi tessellation) but not the number of cells; and a remapping phase in which the Lagrangian solution is transferred onto the new grid. Furthermore, in the standard ReALEmore » method, the rezoned mesh is smoothed by using one or several steps toward centroidal Voronoi tessellation, but it is not adapted to the solution in any way.« less

Adaptive reconnection-based arbitrary Lagrangian Eulerian method

DOE PAGES

Bo, Wurigen; Shashkov, Mikhail

2015-07-21

We present a new adaptive Arbitrary Lagrangian Eulerian (ALE) method. This method is based on the reconnection-based ALE (ReALE) methodology of Refs. [35], [34] and [6]. The main elements in a standard ReALE method are: an explicit Lagrangian phase on an arbitrary polygonal (in 2D) mesh in which the solution and positions of grid nodes are updated; a rezoning phase in which a new grid is defined by changing the connectivity (using Voronoi tessellation) but not the number of cells; and a remapping phase in which the Lagrangian solution is transferred onto the new grid. Furthermore, in the standard ReALEmore » method, the rezoned mesh is smoothed by using one or several steps toward centroidal Voronoi tessellation, but it is not adapted to the solution in any way.« less

Vortex dynamics and Lagrangian statistics in a model for active turbulence.

PubMed

James, Martin; Wilczek, Michael

2018-02-14

Cellular suspensions such as dense bacterial flows exhibit a turbulence-like phase under certain conditions. We study this phenomenon of "active turbulence" statistically by using numerical tools. Following Wensink et al. (Proc. Natl. Acad. Sci. U.S.A. 109, 14308 (2012)), we model active turbulence by means of a generalized Navier-Stokes equation. Two-point velocity statistics of active turbulence, both in the Eulerian and the Lagrangian frame, is explored. We characterize the scale-dependent features of two-point statistics in this system. Furthermore, we extend this statistical study with measurements of vortex dynamics in this system. Our observations suggest that the large-scale statistics of active turbulence is close to Gaussian with sub-Gaussian tails.

On Lagrangian residual currents with applications in south San Francisco Bay, California

USGS Publications Warehouse

Cheng, Ralph T.; Casulli, Vincenzo

1982-01-01

The Lagrangian residual circulation has often been introduced as the sum of the Eulerian residual circulation and the Stokes' drift. Unfortunately, this definition of the Lagrangian residual circulation is conceptually incorrect because both the Eulerian residual circulation and the Stokes' drift are Eulerian variables. In this paper a classification of various residual variables are reviewed and properly defined. The Lagrangian residual circulation is then studied by means of a two-stage formulation of a computer model. The tidal circulation is first computed in a conventional Eulerian way, and then the Lagrangian residual circulation is determined by a method patterned after the method of markers and cells. To demonstrate properties of the Lagrangian residual circulation, application of this approach in South San Francisco Bay, California, is considered. With the aid of the model results, properties of the Eulerian and Lagrangian residual circulation are examined. It can be concluded that estimation of the Lagrangian residual circulation from Eulerian data may lead to unacceptable error, particularly in a tidal estuary where the tidal excursion is of the same order of magnitude as the length scale of the basin. A direction calculation of the Lagrangian residual circulation must be made and has been shown to be feasible.

General Lagrangian formulation for higher spin fields with arbitrary index symmetry. 2. Fermionic fields

NASA Astrophysics Data System (ADS)

Reshetnyak, A.

2013-04-01

We continue the construction of a Lagrangian description of irreducible half-integer higher-spin representations of the Poincare group with an arbitrary Young tableaux having k rows, on a basis of the BRST-BFV approach suggested for bosonic fields in our first article [I.L. Buchbinder, A. Reshetnyak, Nucl. Phys. B 862 (2012) 270, arXiv:1110.5044 [hep-th

Parent formulation at the Lagrangian level

NASA Astrophysics Data System (ADS)

Grigoriev, Maxim

2011-07-01

The recently proposed first-order parent formalism at the level of equations of motion is specialized to the case of Lagrangian systems. It is shown that for diffeomorphism-invariant theories the parent formulation takes the form of an AKSZ-type sigma model. The proposed formulation can be also seen as a Lagrangian version of the BV-BRST extension of the Vasiliev unfolded approach. We also discuss its possible interpretation as a multidimensional generalization of the Hamiltonian BFV-BRST formalism. The general construction is illustrated by examples of (parametrized) mechanics, relativistic particle, Yang-Mills theory, and gravity.

Thermostating extended Lagrangian Born-Oppenheimer molecular dynamics

DOE Office of Scientific and Technical Information (OSTI.GOV)

Martínez, Enrique; Cawkwell, Marc J.; Voter, Arthur F.

Here, Extended Lagrangian Born-Oppenheimer molecular dynamics is developed and analyzed for applications in canonical (NVT) simulations. Three different approaches are considered: the Nosé and Andersen thermostats and Langevin dynamics. We have tested the temperature distribution under different conditions of self-consistent field (SCF) convergence and time step and compared the results to analytical predictions. We find that the simulations based on the extended Lagrangian Born-Oppenheimer framework provide accurate canonical distributions even under approximate SCF convergence, often requiring only a single diagonalization per time step, whereas regular Born-Oppenheimer formulations exhibit unphysical fluctuations unless a sufficiently high degree of convergence is reached atmore » each time step. Lastly, the thermostated extended Lagrangian framework thus offers an accurate approach to sample processes in the canonical ensemble at a fraction of the computational cost of regular Born-Oppenheimer molecular dynamics simulations.« less

Thermostating extended Lagrangian Born-Oppenheimer molecular dynamics

DOE PAGES

Martínez, Enrique; Cawkwell, Marc J.; Voter, Arthur F.; ...

2015-04-21

Here, Extended Lagrangian Born-Oppenheimer molecular dynamics is developed and analyzed for applications in canonical (NVT) simulations. Three different approaches are considered: the Nosé and Andersen thermostats and Langevin dynamics. We have tested the temperature distribution under different conditions of self-consistent field (SCF) convergence and time step and compared the results to analytical predictions. We find that the simulations based on the extended Lagrangian Born-Oppenheimer framework provide accurate canonical distributions even under approximate SCF convergence, often requiring only a single diagonalization per time step, whereas regular Born-Oppenheimer formulations exhibit unphysical fluctuations unless a sufficiently high degree of convergence is reached atmore » each time step. Lastly, the thermostated extended Lagrangian framework thus offers an accurate approach to sample processes in the canonical ensemble at a fraction of the computational cost of regular Born-Oppenheimer molecular dynamics simulations.« less

The spatial gradients in the solar wind and IMF in the vicinity of the first Lagrangian point

NASA Astrophysics Data System (ADS)

Lai, H.; Russell, C. T.; Riley, P.

2017-12-01

To verify the accuracy of predicted solar wind conditions at L1, we need to know how accurate our measurements are as well as the spatial gradients of solar wind properties since the data are not obtained precisely at the L1 point. With ACE, Wind, and DSCOVR currently taking measurements in the vicinity of L1, we first need to test whether their responses to the solar wind are the same and if not, to determine which data are most accurate. Secondly, we need to study the coherency scales of the solar wind properties, which determine the scale over which the measurements can be accurately extrapolated. By comparing the measurements during large solar wind structures (e.g. CMEs), we find that the magnetic fields from all spacecraft are measured accurately, but the plasma parameters can be significantly different from one spacecraft to another. By examining the sum of magnetic and plasma thermal pressure across tangential discontinuities, we find that the density and temperature measurements from Wind and DSCOVR do show pressure continuity as expected while ACE does not. Since plasma data from DSCOVR have a greater variability about the mean and have many data gaps, we believe that data from Wind should be used whenever available. We find that strength of the magnetic field and zero levels of the various magnetometers are consistent, but the direction of the magnetic field can change significantly in the cross-flow direction. Thus, over the separation distance of spacecraft near L1, large changes in the IMF direction can appear between spacecraft even though the IMF is accurately measured. In contrast, the plasma parameters, when measured accurately, are spatially uniform over about 100Re and may be extrapolated well. Our results can also be applied to improving future space weather mission design. A constellation of cubesats with magnetometers would be needed to determine the IMF impinging on the magnetosphere. Fewer plasma instruments are needed to determine the

Comparing L1 and L2 Texts and Writers in First-Year Composition

ERIC Educational Resources Information Center

Eckstein, Grant; Ferris, Dana

2018-01-01

Scholars have at various points discussed the needs of second language (L2) writers enrolled in "mainstream" composition courses where they are mixed with native (L1) English speakers. Other researchers have investigated the experiences of L2 writers in mainstream classes and the perceptions of their instructors about their abilities and…

Contractions of AdS brane algebra and superGalileon Lagrangians

DOE Office of Scientific and Technical Information (OSTI.GOV)

Kamimura, Kiyoshi; Onda, Seiji

2013-06-15

We examine AdS Galileon Lagrangians using the method of nonlinear realization. By contractions (1) flat curvature limit, (2) non-relativistic brane algebra limit, and (3) (1) + (2) limits we obtain DBI, Newton-Hoock, and Galilean Galileons, respectively. We make clear how these Lagrangians appear as invariant 4-forms and/or pseudo-invariant Wess-Zumino (WZ) terms using Maurer-Cartan (MC) equations on the coset G/SO(3, 1). We show the equations of motion are written in terms of the MC forms only and explain why the inverse Higgs condition is obtained as the equation of motion for all cases. The supersymmetric extension is also examined using amore » supercoset SU(2, 2 Double-Vertical-Line 1)/(SO(3, 1) Multiplication-Sign U(1)) and five WZ forms are constructed. They are reduced to the corresponding five Galileon WZ forms in the bosonic limit and are candidates for supersymmetric Galileon action.« less

Assimilating Eulerian and Lagrangian data in traffic-flow models

NASA Astrophysics Data System (ADS)

Xia, Chao; Cochrane, Courtney; DeGuire, Joseph; Fan, Gaoyang; Holmes, Emma; McGuirl, Melissa; Murphy, Patrick; Palmer, Jenna; Carter, Paul; Slivinski, Laura; Sandstede, Björn

2017-05-01

Data assimilation of traffic flow remains a challenging problem. One difficulty is that data come from different sources ranging from stationary sensors and camera data to GPS and cell phone data from moving cars. Sensors and cameras give information about traffic density, while GPS data provide information about the positions and velocities of individual cars. Previous methods for assimilating Lagrangian data collected from individual cars relied on specific properties of the underlying computational model or its reformulation in Lagrangian coordinates. These approaches make it hard to assimilate both Eulerian density and Lagrangian positional data simultaneously. In this paper, we propose an alternative approach that allows us to assimilate both Eulerian and Lagrangian data. We show that the proposed algorithm is accurate and works well in different traffic scenarios and regardless of whether ensemble Kalman or particle filters are used. We also show that the algorithm is capable of estimating parameters and assimilating real traffic observations and synthetic observations obtained from microscopic models.

Effective Lagrangians and Current Algebra in Three Dimensions

NASA Astrophysics Data System (ADS)

Ferretti, Gabriele

In this thesis we study three dimensional field theories that arise as effective Lagrangians of quantum chromodynamics in Minkowski space with signature (2,1) (QCD3). In the first chapter, we explain the method of effective Langrangians and the relevance of current algebra techniques to field theory. We also provide the physical motivations for the study of QCD3 as a toy model for confinement and as a theory of quantum antiferromagnets (QAF). In chapter two, we derive the relevant effective Lagrangian by studying the low energy behavior of QCD3, paying particular attention to how the global symmetries are realized at the quantum level. In chapter three, we show how baryons arise as topological solitons of the effective Lagrangian and also show that their statistics depends on the number of colors as predicted by the quark model. We calculate mass splitting and magnetic moments of the soliton and find logarithmic corrections to the naive quark model predictions. In chapter four, we drive the current algebra of the theory. We find that the current algebra is a co -homologically non-trivial generalization of Kac-Moody algebras to three dimensions. This fact may provide a new, non -perturbative way to quantize the theory. In chapter five, we discuss the renormalizability of the model in the large-N expansion. We prove the validity of the non-renormalization theorem and compute the critical exponents in a specific limiting case, the CP^ {N-1} model with a Chern-Simons term. Finally, chapter six contains some brief concluding remarks.

A macroscopic plasma Lagrangian and its application to wave interactions and resonances

NASA Technical Reports Server (NTRS)

Peng, Y. K. M.

1974-01-01

The derivation of a macroscopic plasma Lagrangian is considered, along with its application to the description of nonlinear three-wave interaction in a homogeneous plasma and linear resonance oscillations in a inhomogeneous plasma. One approach to obtain the Lagrangian is via the inverse problem of the calculus of variations for arbitrary first and second order quasilinear partial differential systems. Necessary and sufficient conditions for the given equations to be Euler-Lagrange equations of a Lagrangian are obtained. These conditions are then used to determine the transformations that convert some classes of non-Euler-Lagrange equations to Euler-Lagrange equation form. The Lagrangians for a linear resistive transmission line and a linear warm collisional plasma are derived as examples. Using energy considerations, the correct macroscopic plasma Lagrangian is shown to differ from the velocity-integrated low Lagrangian by a macroscopic potential energy that equals twice the particle thermal kinetic energy plus the energy lost by heat conduction.

Sensitivity Analysis of a Lagrangian Sea Ice Model

NASA Astrophysics Data System (ADS)

Rabatel, Matthias; Rampal, Pierre; Bertino, Laurent; Carrassi, Alberto; Jones, Christopher K. R. T.

2017-04-01

Large changes in the Arctic sea ice have been observed in the last decades in terms of the ice thickness, extension and drift. Understanding the mechanisms behind these changes is of paramount importance to enhance our modeling and forecasting capabilities. For 40 years, models have been developed to describe the non-linear dynamical response of the sea ice to a number of external and internal factors. Nevertheless, there still exists large deviations between predictions and observations. There are related to incorrect descriptions of the sea ice response and/or to the uncertainties about the different sources of information: parameters, initial and boundary conditions and external forcing. Data assimilation (DA) methods are used to combine observations with models, and there is nowadays an increasing interest of DA for sea-ice models and observations. We consider here the state-of-the art sea-ice model, neXtSIM te{Rampal2016a}, which is based on a time-varying Lagrangian mesh and makes use of the Elasto-Brittle rheology. Our ultimate goal is designing appropriate DA scheme for such a modelling facility. This contribution reports about the first milestone along this line: a sensitivity analysis in order to quantify forecast error to guide model development and to set basis for further Lagrangian DA methods. Specific features of the sea-ice dynamics in relation to the wind are thus analysed. Virtual buoys are deployed across the Arctic domain and their trajectories of motion are analysed. The simulated trajectories are also compared to real buoys trajectories observed. The model response is also compared with that one from a model version not including internal forcing to highlight the role of the rheology. Conclusions and perspectives for the general DA implementation are also discussed. \\bibitem{Rampal2016a} P. Rampal, S. Bouillon, E. Ólason, and M. Morlighem. ne{X}t{SIM}: a new {L}agrangian sea ice model. The Cryosphere, 10 (3): 1055-1073, 2016.

Spacecraft Maneuvering at the Sun/Earth-Moon L2 Libration Point

NASA Astrophysics Data System (ADS)

Shahid, Kamran

of the proposed control techniques for spacecraft maneuvering using solar radiation pressure at the L2 libration point. Stationkeeping accuracies of 50m and formation maintenance accuracies of less than 1m are possible using solar radiation pressure at a sub-L2 libration point. The benefits of these control techniques include increasing libration point mission lifetimes and doubling payload mass fractions as compared to conventional propulsion methods.

Next generation extended Lagrangian first principles molecular dynamics

NASA Astrophysics Data System (ADS)

Niklasson, Anders M. N.

2017-08-01

Extended Lagrangian Born-Oppenheimer molecular dynamics [A. M. N. Niklasson, Phys. Rev. Lett. 100, 123004 (2008)] is formulated for general Hohenberg-Kohn density-functional theory and compared with the extended Lagrangian framework of first principles molecular dynamics by Car and Parrinello [Phys. Rev. Lett. 55, 2471 (1985)]. It is shown how extended Lagrangian Born-Oppenheimer molecular dynamics overcomes several shortcomings of regular, direct Born-Oppenheimer molecular dynamics, while improving or maintaining important features of Car-Parrinello simulations. The accuracy of the electronic degrees of freedom in extended Lagrangian Born-Oppenheimer molecular dynamics, with respect to the exact Born-Oppenheimer solution, is of second-order in the size of the integration time step and of fourth order in the potential energy surface. Improved stability over recent formulations of extended Lagrangian Born-Oppenheimer molecular dynamics is achieved by generalizing the theory to finite temperature ensembles, using fractional occupation numbers in the calculation of the inner-product kernel of the extended harmonic oscillator that appears as a preconditioner in the electronic equations of motion. Material systems that normally exhibit slow self-consistent field convergence can be simulated using integration time steps of the same order as in direct Born-Oppenheimer molecular dynamics, but without the requirement of an iterative, non-linear electronic ground-state optimization prior to the force evaluations and without a systematic drift in the total energy. In combination with proposed low-rank and on the fly updates of the kernel, this formulation provides an efficient and general framework for quantum-based Born-Oppenheimer molecular dynamics simulations.

Next generation extended Lagrangian first principles molecular dynamics.

PubMed

Niklasson, Anders M N

2017-08-07

Extended Lagrangian Born-Oppenheimer molecular dynamics [A. M. N. Niklasson, Phys. Rev. Lett. 100, 123004 (2008)] is formulated for general Hohenberg-Kohn density-functional theory and compared with the extended Lagrangian framework of first principles molecular dynamics by Car and Parrinello [Phys. Rev. Lett. 55, 2471 (1985)]. It is shown how extended Lagrangian Born-Oppenheimer molecular dynamics overcomes several shortcomings of regular, direct Born-Oppenheimer molecular dynamics, while improving or maintaining important features of Car-Parrinello simulations. The accuracy of the electronic degrees of freedom in extended Lagrangian Born-Oppenheimer molecular dynamics, with respect to the exact Born-Oppenheimer solution, is of second-order in the size of the integration time step and of fourth order in the potential energy surface. Improved stability over recent formulations of extended Lagrangian Born-Oppenheimer molecular dynamics is achieved by generalizing the theory to finite temperature ensembles, using fractional occupation numbers in the calculation of the inner-product kernel of the extended harmonic oscillator that appears as a preconditioner in the electronic equations of motion. Material systems that normally exhibit slow self-consistent field convergence can be simulated using integration time steps of the same order as in direct Born-Oppenheimer molecular dynamics, but without the requirement of an iterative, non-linear electronic ground-state optimization prior to the force evaluations and without a systematic drift in the total energy. In combination with proposed low-rank and on the fly updates of the kernel, this formulation provides an efficient and general framework for quantum-based Born-Oppenheimer molecular dynamics simulations.

A purely Lagrangian method for computing linearly-perturbed flows in spherical geometry

NASA Astrophysics Data System (ADS)

Jaouen, Stéphane

2007-07-01

In many physical applications, one wishes to control the development of multi-dimensional instabilities around a one-dimensional (1D) complex flow. For predicting the growth rates of these perturbations, a general numerical approach is viable which consists in solving simultaneously the one-dimensional equations and their linearized form for three-dimensional perturbations. In Clarisse et al. [J.-M. Clarisse, S. Jaouen, P.-A. Raviart, A Godunov-type method in Lagrangian coordinates for computing linearly-perturbed planar-symmetric flows of gas dynamics, J. Comp. Phys. 198 (2004) 80-105], a class of Godunov-type schemes for planar-symmetric flows of gas dynamics has been proposed. Pursuing this effort, we extend these results to spherically symmetric flows. A new method to derive the Lagrangian perturbation equations, based on the canonical form of systems of conservation laws with zero entropy flux [B. Després, Lagrangian systems of conservation laws. Invariance properties of Lagrangian systems of conservation laws, approximate Riemann solvers and the entropy condition, Numer. Math. 89 (2001) 99-134; B. Després, C. Mazeran, Lagrangian gas dynamics in two dimensions and Lagrangian systems, Arch. Rational Mech. Anal. 178 (2005) 327-372] is also described. It leads to many advantages. First of all, many physical problems we are interested in enter this formalism (gas dynamics, two-temperature plasma equations, ideal magnetohydrodynamics, etc.) whatever is the geometry. Secondly, a class of numerical entropic schemes is available for the basic flow [11]. Last, linearizing and devising numerical schemes for the perturbed flow is straightforward. The numerical capabilities of these methods are illustrated on three test cases of increasing difficulties and we show that - due to its simplicity and its low computational cost - the Linear Perturbations Code (LPC) is a powerful tool to understand and predict the development of hydrodynamic instabilities in the linear regime.

Evaluation of wastewater contaminant transport in surface waters using verified Lagrangian sampling

USGS Publications Warehouse

Antweiler, Ronald C.; Writer, Jeffrey H.; Murphy, Sheila F.

2014-01-01

Contaminants released from wastewater treatment plants can persist in surface waters for substantial distances. Much research has gone into evaluating the fate and transport of these contaminants, but this work has often assumed constant flow from wastewater treatment plants. However, effluent discharge commonly varies widely over a 24-hour period, and this variation controls contaminant loading and can profoundly influence interpretations of environmental data. We show that methodologies relying on the normalization of downstream data to conservative elements can give spurious results, and should not be used unless it can be verified that the same parcel of water was sampled. Lagrangian sampling, which in theory samples the same water parcel as it moves downstream (the Lagrangian parcel), links hydrologic and chemical transformation processes so that the in-stream fate of wastewater contaminants can be quantitatively evaluated. However, precise Lagrangian sampling is difficult, and small deviations – such as missing the Lagrangian parcel by less than 1 h – can cause large differences in measured concentrations of all dissolved compounds at downstream sites, leading to erroneous conclusions regarding in-stream processes controlling the fate and transport of wastewater contaminants. Therefore, we have developed a method termed “verified Lagrangian” sampling, which can be used to determine if the Lagrangian parcel was actually sampled, and if it was not, a means for correcting the data to reflect the concentrations which would have been obtained had the Lagrangian parcel been sampled. To apply the method, it is necessary to have concentration data for a number of conservative constituents from the upstream, effluent, and downstream sites, along with upstream and effluent concentrations that are constant over the short-term (typically 2–4 h). These corrections can subsequently be applied to all data, including non-conservative constituents. Finally, we

Microscopic Lagrangian description of warm plasmas. IV - Macroscopic approximation

NASA Technical Reports Server (NTRS)

Kim, H.; Crawford, F. W.

1983-01-01

The averaged-Lagrangian method is applied to linear wave propagation and nonlinear three-wave interaction in a warm magnetoplasma, in the macroscopic approximation. The microscopic Lagrangian treated by Kim and Crawford (1977) and by Galloway and Crawford (1977) is first expanded to third order in perturbation. Velocity integration is then carried out, before applying Hamilton's principle to obtain a general description of wave propagation and coupling. The results are specialized to the case of interaction between two electron plasma waves and an Alfven wave. The method is shown to be more powerful than the alternative possibility of working from the beginning with a macroscopic Lagrangian density.

Lagrangian acceleration statistics in a turbulent channel flow

NASA Astrophysics Data System (ADS)

Stelzenmuller, Nickolas; Polanco, Juan Ignacio; Vignal, Laure; Vinkovic, Ivana; Mordant, Nicolas

2017-05-01

Lagrangian acceleration statistics in a fully developed turbulent channel flow at Reτ=1440 are investigated, based on tracer particle tracking in experiments and direct numerical simulations. The evolution with wall distance of the Lagrangian velocity and acceleration time scales is analyzed. Dependency between acceleration components in the near-wall region is described using cross-correlations and joint probability density functions. The strong streamwise coherent vortices typical of wall-bounded turbulent flows are shown to have a significant impact on the dynamics. This results in a strong anisotropy at small scales in the near-wall region that remains present in most of the channel. Such statistical properties may be used as constraints in building advanced Lagrangian stochastic models to predict the dispersion and mixing of chemical components for combustion or environmental studies.

Effective Lagrangian in de Sitter spacetime

NASA Astrophysics Data System (ADS)

Kitamoto, Hiroyuki; Kitazawa, Yoshihisa

2017-01-01

Scale invariant fluctuations of metric are a universal feature of quantum gravity in de Sitter spacetime. We construct an effective Lagrangian which summarizes their implications on local physics by integrating superhorizon metric fluctuations. It shows infrared quantum effects are local and render fundamental couplings time dependent. We impose Lorenz invariance on the effective Lagrangian as it is required by the principle of general covariance. We show that such a requirement leads to unique physical predictions by fixing the quantization ambiguities. We explain how the gauge parameter dependence of observables is canceled. In particular the relative evolution speed of the couplings are shown to be gauge invariant.

Quantization of Non-Lagrangian Systems

NASA Astrophysics Data System (ADS)

Kochan, Denis

A novel method for quantization of non-Lagrangian (open) systems is proposed. It is argued that the essential object, which provides both classical and quantum evolution, is a certain canonical two-form defined in extended velocity space. In this setting classical dynamics is recovered from the stringy-type variational principle, which employs umbilical surfaces instead of histories of the system. Quantization is then accomplished in accordance with the introduced variational principle. The path integral for the transition probability amplitude (propagator) is rearranged to a surface functional integral. In the standard case of closed (Lagrangian) systems the presented method reduces to the standard Feynman's approach. The inverse problem of the calculus of variation, the problem of quantization ambiguity and the quantum mechanics in the presence of friction are analyzed in detail.

The S-Lagrangian and a theory of homeostasis in living systems

NASA Astrophysics Data System (ADS)

Sandler, U.; Tsitolovsky, L.

2017-04-01

A major paradox of living things is their ability to actively counteract degradation in a continuously changing environment or being injured through homeostatic protection. In this study, we propose a dynamic theory of homeostasis based on a generalized Lagrangian approach (S-Lagrangian), which can be equally applied to physical and nonphysical systems. Following discoverer of homeostasis Cannon (1935), we assume that homeostasis results from tendency of the organisms to decrease of the stress and avoid of death. We show that the universality of homeostasis is a consequence of analytical properties of the S-Lagrangian, while peculiarities of the biochemical and physiological mechanisms of homeostasis determine phenomenological parameters of the S-Lagrangian. Additionally, we reveal that plausible assumptions about S-Lagrangian features lead to good agreement between theoretical descriptions and observed homeostatic behavior. Here, we have focused on homeostasis of living systems, however, the proposed theory is also capable of being extended to social systems.

A radiative seesaw model with higher order terms under an alternative U(1)B-L

NASA Astrophysics Data System (ADS)

Nomura, Takaaki; Okada, Hiroshi

2018-06-01

We propose a model based on an alternative U(1) B - L gauge symmetry with 5 dimensional operators in the Lagrangian, and we construct the neutrino masses at one-loop level, and discuss lepton flavor violations, dark matter, and the effective number of neutrino species due to two massless particles in our model. Then we search allowed region to satisfy the current experimental data of neutrino oscillation and lepton flavor violations without conflict of several constraints such as stability of dark matter and the effective number of neutrino species, depending on normal hierarchy and inverted one.

Sigma decomposition: the CP-odd Lagrangian

NASA Astrophysics Data System (ADS)

Hierro, I. M.; Merlo, L.; Rigolin, S.

2016-04-01

In Alonso et al., JHEP 12 (2014) 034, the CP-even sector of the effective chiral Lagrangian for a generic composite Higgs model with a symmetric coset has been constructed, up to four momenta. In this paper, the CP-odd couplings are studied within the same context. If only the Standard Model bosonic sources of custodial symmetry breaking are considered, then at most six independent operators form a basis. One of them is the weak- θ term linked to non-perturbative sources of CP violation, while the others describe CP-odd perturbative couplings between the Standard Model gauge bosons and an Higgs-like scalar belonging to the Goldstone boson sector. The procedure is then applied to three distinct exemplifying frameworks: the original SU(5)/SO(5) Georgi-Kaplan model, the minimal custodial-preserving SO(5)/SO(4) model and the minimal SU(3)/(SU(2) × U(1)) model, which intrinsically breaks custodial symmetry. Moreover, the projection of the high-energy electroweak effective theory to the low-energy chiral effective Lagrangian for a dynamical Higgs is performed, uncovering strong relations between the operator coefficients and pinpointing the differences with the elementary Higgs scenario.

Point mutations in acetylcholinesterase 1 associated with chlorpyrifos resistance in the brown planthopper, Nilaparvata lugens Stål.

PubMed

Zhang, Y; Yang, B; Li, J; Liu, M; Liu, Z

2017-08-01

Insecticide resistance frequently results from target-site insensitivity, such as point mutations in acetylcholinesterases (AChEs) for resistance to organophosphates and carbamates. From a field-originated population of Nilaparvata lugens, a major rice pest, a resistant population (R9) was obtained by nine-generation continuous selection with chlorpyrifos. From the same field population, a relatively susceptible population (S9) was also constructed through rearing without any insecticides. Compared to the susceptible strain, Sus [medium lethal dose (LC 50 ) = 0.012 mg/l], R9 had a resistance ratio (RR) of 253.08-fold, whereas the RR of S9 was only 2.25-fold. Piperonyl butoxide and triphenyl phosphate synergized chlorpyrifos in R9 less than three-fold, indicating other important mechanisms for high resistance. The target-site insensitivity was supported by the key property differences of crude AChEs between R9 and S9. Compared to S9, three mutations (G119S, F331C and I332L) were detected in NlAChE1 from individuals of the R9 and field populations, but no mutation was detected in NlAChE2. G119S and F331C could decreased insecticide sensitivities in recombinant NlAChE1, whereas I332L took effect through increasing the influence of F331C on target insensitivity. F331C might be deleterious because of its influence on the catalytic efficiency of NlAChE1, whereas I332L would decrease these adverse effects and maintain the normal functions of AChEs. © 2017 The Royal Entomological Society.

Stress-Induced Acoustic Variation in L2 and L1 Spanish Vowels.

PubMed

Romanelli, Sofía; Menegotto, Andrea; Smyth, Ron

2018-05-28

We assessed the effect of lexical stress on the duration and quality of Spanish word-final vowels /a, e, o/ produced by American English late intermediate learners of L2 Spanish, as compared to those of native L1 Argentine Spanish speakers. Participants read 54 real words ending in /a, e, o/, with either final or penultimate lexical stress, embedded in a text and a word list. We measured vowel duration and both F1 and F2 frequencies at 3 temporal points. stressed vowels were longer than unstressed vowels, in Spanish L1 and L2. L1 and L2 Spanish stressed /a/ and /e/ had higher F1 values than their unstressed counterparts. Only the L2 speakers showed evidence of rising offglides for /e/ and /o/. The L2 and L1 Spanish vowel space was compressed in the absence of stress. Lexical stress affected the vowel quality of L1 and L2 Spanish vowels. We provide an up-to-date account of the formant trajectories of Argentine River Plate Spanish word-final /a, e, o/ and offer experimental support to the claim that stress affects the quality of Spanish vowels in word-final contexts. © 2018 S. Karger AG, Basel.

Lagrangian methods in nonlinear plasma wave interaction

NASA Technical Reports Server (NTRS)

Crawford, F. W.

1980-01-01

Analysis of nonlinear plasma wave interactions is usually very complicated, and simplifying mathematical approaches are highly desirable. The application of averaged-Lagrangian methods offers a considerable reduction in effort, with improved insight into synchronism and conservation (Manley-Rowe) relations. This chapter indicates how suitable Lagrangian densities have been defined, expanded, and manipulated to describe nonlinear wave-wave and wave-particle interactions in the microscopic, macroscopic and cold plasma models. Recently, further simplifications have been introduced by the use of techniques derived from Lie algebra. These and likely future developments are reviewed briefly.

Using Lagrangian Coherent Structures to understand coastal water quality

NASA Astrophysics Data System (ADS)

Fiorentino, L. A.; Olascoaga, M. J.; Reniers, A.; Feng, Z.; Beron-Vera, F. J.; MacMahan, J. H.

2012-09-01

The accumulation of pollutants near the shoreline can result in low quality coastal water with negative effects on human health. To understand the role of mixing by tidal flows in coastal water quality we study the nearshore Lagrangian circulation. Specifically, we reveal Lagrangian Coherent Structures (LCSs), i.e., distinguished material curves which shape global mixing patterns and thus act as skeletons of the Lagrangian circulation. This is done using the recently developed geodesic theory of transport barriers. Particular focus is placed on Hobie Beach, a recreational subtropical marine beach located in Virginia Key, Miami, Florida. According to studies of water quality, Hobie Beach is characterized by high microbial levels. Possible sources of pollution in Hobie Beach include human bather shedding, dog fecal matter, runoff, and sand efflux at high tides. Consistent with the patterns formed by satellite-tracked drifter trajectories, the LCSs extracted from simulated currents reveal a Lagrangian circulation favoring the retention near the shoreline of pollutants released along the shoreline, which can help explain the low quality water registered at Hobie Beach.

Integration over families of Lagrangian submanifolds in BV formalism

NASA Astrophysics Data System (ADS)

Mikhailov, Andrei

2018-03-01

Gauge fixing is interpreted in BV formalism as a choice of Lagrangian submanifold in an odd symplectic manifold (the BV phase space). A natural construction defines an integration procedure on families of Lagrangian submanifolds. In string perturbation theory, the moduli space integrals of higher genus amplitudes can be interpreted in this way. We discuss the role of gauge symmetries in this construction. We derive the conditions which should be imposed on gauge symmetries for the consistency of our integration procedure. We explain how these conditions behave under the deformations of the worldsheet theory. In particular, we show that integrated vertex operator is actually an inhomogeneous differential form on the space of Lagrangian submanifolds.

Symmetries in Lagrangian Dynamics

ERIC Educational Resources Information Center

Ferrario, Carlo; Passerini, Arianna

2007-01-01

In the framework of Noether's theorem, a distinction between Lagrangian and dynamical symmetries is made, in order to clarify some aspects neglected by textbooks. An intuitive setting of the concept of invariance of differential equations is presented. The analysis is completed by deriving the symmetry properties in the motion of a charged…

Testing higher-order Lagrangian perturbation theory against numerical simulation. 1: Pancake models

NASA Technical Reports Server (NTRS)

Buchert, T.; Melott, A. L.; Weiss, A. G.

1993-01-01

We present results showing an improvement of the accuracy of perturbation theory as applied to cosmological structure formation for a useful range of quasi-linear scales. The Lagrangian theory of gravitational instability of an Einstein-de Sitter dust cosmogony investigated and solved up to the third order is compared with numerical simulations. In this paper we study the dynamics of pancake models as a first step. In previous work the accuracy of several analytical approximations for the modeling of large-scale structure in the mildly non-linear regime was analyzed in the same way, allowing for direct comparison of the accuracy of various approximations. In particular, the Zel'dovich approximation (hereafter ZA) as a subclass of the first-order Lagrangian perturbation solutions was found to provide an excellent approximation to the density field in the mildly non-linear regime (i.e. up to a linear r.m.s. density contrast of sigma is approximately 2). The performance of ZA in hierarchical clustering models can be greatly improved by truncating the initial power spectrum (smoothing the initial data). We here explore whether this approximation can be further improved with higher-order corrections in the displacement mapping from homogeneity. We study a single pancake model (truncated power-spectrum with power-spectrum with power-index n = -1) using cross-correlation statistics employed in previous work. We found that for all statistical methods used the higher-order corrections improve the results obtained for the first-order solution up to the stage when sigma (linear theory) is approximately 1. While this improvement can be seen for all spatial scales, later stages retain this feature only above a certain scale which is increasing with time. However, third-order is not much improvement over second-order at any stage. The total breakdown of the perturbation approach is observed at the stage, where sigma (linear theory) is approximately 2, which corresponds to the

The Monotonic Lagrangian Grid for Rapid Air-Traffic Evaluation

NASA Technical Reports Server (NTRS)

Kaplan, Carolyn; Dahm, Johann; Oran, Elaine; Alexandrov, Natalia; Boris, Jay

2010-01-01

The Air Traffic Monotonic Lagrangian Grid (ATMLG) is presented as a tool to evaluate new air traffic system concepts. The model, based on an algorithm called the Monotonic Lagrangian Grid (MLG), can quickly sort, track, and update positions of many aircraft, both on the ground (at airports) and in the air. The underlying data structure is based on the MLG, which is used for sorting and ordering positions and other data needed to describe N moving bodies and their interactions. Aircraft that are close to each other in physical space are always near neighbors in the MLG data arrays, resulting in a fast nearest-neighbor interaction algorithm that scales as N. Recent upgrades to ATMLG include adding blank place-holders within the MLG data structure, which makes it possible to dynamically change the MLG size and also improves the quality of the MLG grid. Additional upgrades include adding FAA flight plan data, such as way-points and arrival and departure times from the Enhanced Traffic Management System (ETMS), and combining the MLG with the state-of-the-art strategic and tactical conflict detection and resolution algorithms from the NASA-developed Stratway software. In this paper, we present results from our early efforts to couple ATMLG with the Stratway software, and we demonstrate that it can be used to quickly simulate air traffic flow for a very large ETMS dataset.

Acoustic streaming: an arbitrary Lagrangian-Eulerian perspective.

PubMed

Nama, Nitesh; Huang, Tony Jun; Costanzo, Francesco

2017-08-25

We analyse acoustic streaming flows using an arbitrary Lagrangian Eulerian (ALE) perspective. The formulation stems from an explicit separation of time scales resulting in two subproblems: a first-order problem, formulated in terms of the fluid displacement at the fast scale, and a second-order problem, formulated in terms of the Lagrangian flow velocity at the slow time scale. Following a rigorous time-averaging procedure, the second-order problem is shown to be intrinsically steady, and with exact boundary conditions at the oscillating walls. Also, as the second-order problem is solved directly for the Lagrangian velocity, the formulation does not need to employ the notion of Stokes drift, or any associated post-processing, thus facilitating a direct comparison with experiments. Because the first-order problem is formulated in terms of the displacement field, our formulation is directly applicable to more complex fluid-structure interaction problems in microacoustofluidic devices. After the formulation's exposition, we present numerical results that illustrate the advantages of the formulation with respect to current approaches.

Generalized Lagrangian coherent structures

NASA Astrophysics Data System (ADS)

Balasuriya, Sanjeeva; Ouellette, Nicholas T.; Rypina, Irina I.

2018-06-01

The notion of a Lagrangian Coherent Structure (LCS) is by now well established as a way to capture transient coherent transport dynamics in unsteady and aperiodic fluid flows that are known over finite time. We show that the concept of an LCS can be generalized to capture coherence in other quantities of interest that are transported by, but not fully locked to, the fluid. Such quantities include those with dynamic, biological, chemical, or thermodynamic relevance, such as temperature, pollutant concentration, vorticity, kinetic energy, plankton density, and so on. We provide a conceptual framework for identifying the Generalized Lagrangian Coherent Structures (GLCSs) associated with such evolving quantities. We show how LCSs can be seen as a special case within this framework, and provide an overarching discussion of various methods for identifying LCSs. The utility of this more general viewpoint is highlighted through a variety of examples. We also show that although LCSs approximate GLCSs in certain limiting situations under restrictive assumptions on how the velocity field affects the additional quantities of interest, LCSs are not in general sufficient to describe their coherent transport.

Conceptual Design of a Synoptic Interplanetary Monitor Platform at L sub 1 (SIMPL).

DTIC Science & Technology

1985-11-01

solar events. -159- . . . .. . 105 II1II" -I .5 year mission at Earth-Sun- libration point plus transfer orbit eDashed line is approximate true dose as...Design .. ...................................... 27 4.1 The L Libration Point .......................... 27 4.2 L Orbit Options...34) to provide power, attitude control, communications, and other support to maintain the instruments in a halo orbit around the L libration point ; 4. a

Lagrangian Statistics and Intermittency in Gulf of Mexico.

PubMed

Lin, Liru; Zhuang, Wei; Huang, Yongxiang

2017-12-12

Due to the nonlinear interaction between different flow patterns, for instance, ocean current, meso-scale eddies, waves, etc, the movement of ocean is extremely complex, where a multiscale statistics is then relevant. In this work, a high time-resolution velocity with a time step 15 minutes obtained by the Lagrangian drifter deployed in the Gulf of Mexico (GoM) from July 2012 to October 2012 is considered. The measured Lagrangian velocity correlation function shows a strong daily cycle due to the diurnal tidal cycle. The estimated Fourier power spectrum E(f) implies a dual-power-law behavior which is separated by the daily cycle. The corresponding scaling exponents are close to -1.75 and -2.75 respectively for the time scale larger (resp. 0.1 ≤ f ≤ 0.4 day -1 ) and smaller (resp. 2 ≤ f ≤ 8 day -1 ) than 1 day. A Hilbert-based approach is then applied to this data set to identify the possible multifractal property of the cascade process. The results show an intermittent dynamics for the time scale larger than 1 day, while a less intermittent dynamics for the time scale smaller than 1 day. It is speculated that the energy is partially injected via the diurnal tidal movement and then transferred to larger and small scales through a complex cascade process, which needs more studies in the near future.

Extending geometrical optics: A Lagrangian theory for vector waves

NASA Astrophysics Data System (ADS)

Ruiz, D. E.

2016-10-01

Even diffraction aside, the commonly known equations of geometrical optics (GO) are not entirely accurate. GO considers wave rays as classical particles, which are completely described by their coordinates and momenta, but rays have another degree of freedom, namely, polarization. As a result, wave rays can behave as particles with spin. A well-known example of polarization dynamics is wave-mode conversion, which can be interpreted as rotation of the (classical) ``wave spin.'' However, there are other less-known manifestations of the wave spin, such as polarization precession and polarization-driven bending of ray trajectories. This talk presents recent advances in extending and reformulating GO as a first-principle Lagrangian theory, whose effective-gauge Hamiltonian governs both mentioned polarization phenomena simultaneously. Examples and numerical results are presented. When applied to classical waves, the theory correctly predicts the polarization-driven divergence of left- and right- polarized electromagnetic waves in isotropic media, such as dielectrics and nonmagnetized plasmas. In the case of particles with spin, the formalism also yields a point-particle Lagrangian model for the Dirac electron, i.e. the relativistic spin-1/2 electron, which includes both the Stern-Gerlach spin potential and the Bargmann-Michel-Telegdi spin precession. Additionally, the same theory contributes, perhaps unexpectedly, to the understanding of ponderomotive effects in both wave and particle dynamics; e.g., the formalism allows to obtain the ponderomotive Hamiltonian for a Dirac electron interacting with an arbitrarily large electromagnetic laser field with spin effects included. Supported by the NNSA SSAA Program through DOE Research Grant No. DE-NA0002948, by the U.S. DOE through Contract No. DE-AC02-09CH11466, and by the U.S. DOD NDSEG Fellowship through Contract No. 32-CFR-168a.

78 FR 69516 - DeltaPoint Capital IV, L.P.; License No. 02/02-0662; Notice Seeking Exemption Under the Small...

Federal Register 2010, 2011, 2012, 2013, 2014

2013-11-19

... 107.730). DeltaPoint Capital IV, L.P. provided financing to BioMaxx, Inc., 1 Fishers Road, Suite 160...., an Associate of DeltaPoint Capital IV, L.P., owns more than ten percent of BioMaxx, Inc. Therefore...

Influence of compressibility on the Lagrangian statistics of vorticity-strain-rate interactions.

PubMed

Danish, Mohammad; Sinha, Sawan Suman; Srinivasan, Balaji

2016-07-01

The objective of this study is to investigate the influence of compressibility on Lagrangian statistics of vorticity and strain-rate interactions. The Lagrangian statistics are extracted from "almost" time-continuous data sets of direct numerical simulations of compressible decaying isotropic turbulence by employing a cubic spline-based Lagrangian particle tracker. We study the influence of compressibility on Lagrangian statistics of alignment in terms of compressibility parameters-turbulent Mach number, normalized dilatation-rate, and flow topology. In comparison to incompressible turbulence, we observe that the presence of compressibility in a flow field weakens the alignment tendency of vorticity toward the largest strain-rate eigenvector. Based on the Lagrangian statistics of alignment conditioned on dilatation and topology, we find that the weakened tendency of alignment observed in compressible turbulence is because of a special group of fluid particles that have an initially negligible dilatation-rate and are associated with stable-focus-stretching topology.

A two-field modified Lagrangian formulation for robust simulations of extrinsic cohesive zone models

NASA Astrophysics Data System (ADS)

Cazes, F.; Coret, M.; Combescure, A.

2013-06-01

This paper presents the robust implementation of a cohesive zone model based on extrinsic cohesive laws (i.e. laws involving an infinite initial stiffness). To this end, a two-field Lagrangian weak formulation in which cohesive tractions are chosen as the field variables along the crack's path is presented. Unfortunately, this formulation cannot model the infinite compliance of the broken elements accurately, and no simple criterion can be defined to determine the loading-unloading change of state at the integration points of the cohesive elements. Therefore, a modified Lagrangian formulation using a fictitious cohesive traction instead of the classical cohesive traction as the field variable is proposed. Thanks to this change of variable, the cohesive law becomes an increasing function of the equivalent displacement jump, which eliminates the problems mentioned previously. The ability of the proposed formulations to simulate fracture accurately and without field oscillations is investigated through three numerical test examples.

Lagrangian transport properties of pulmonary interfacial flows

PubMed Central

Smith, Bradford J.; Lukens, Sarah; Yamaguchi, Eiichiro; Gaver, Donald P.

2012-01-01

Disease states characterized by airway fluid occlusion and pulmonary surfactant insufficiency, such as respiratory distress syndrome, have a high mortality rate. Understanding the mechanics of airway reopening, particularly involving surfactant transport, may provide an avenue to increase patient survival via optimized mechanical ventilation waveforms. We model the occluded airway as a liquid-filled rigid tube with the fluid phase displaced by a finger of air that propagates with both mean and sinusoidal velocity components. Finite-time Lyapunov exponent (FTLE) fields are employed to analyse the convective transport characteristics, taking note of Lagrangian coherent structures (LCSs) and their effects on transport. The Lagrangian perspective of these techniques reveals flow characteristics that are not readily apparent by observing Eulerian measures. These analysis techniques are applied to surfactant-free velocity fields determined computationally, with the boundary element method, and measured experimentally with micro particle image velocimetry (μ-PIV). We find that the LCS divides the fluid into two regimes, one advected upstream (into the thin residual film) and the other downstream ahead of the advancing bubble. At higher oscillatory frequencies particles originating immediately inside the LCS experience long residence times at the air–liquid interface, which may be conducive to surfactant transport. At high frequencies a well-mixed attractor region is identified; this volume of fluid cyclically travels along the interface and into the bulk fluid. The Lagrangian analysis is applied to velocity data measured with 0.01 mg ml−1 of the clinical pulmonary surfactant Infasurf in the bulk fluid, demonstrating flow field modifications with respect to the surfactant-free system that were not visible in the Eulerian frame. PMID:23049141

Effective Lagrangian in nonlinear electrodynamics and its properties of causality and unitarity

DOE Office of Scientific and Technical Information (OSTI.GOV)

Shabad, Anatoly E.; Usov, Vladimir V.

2011-05-15

In nonlinear electrodynamics, by implementing the causality principle as the requirement that the group velocity of elementary excitations over a background field should not exceed the speed of light in the vacuum c=1, and the unitarity principle as the requirement that the residue of the propagator should be nonnegative, we establish the positive convexity of the effective Lagrangian on the class of constant fields, also the positivity of all characteristic dielectric and magnetic permittivity constants that are derivatives of the effective Lagrangian with respect to the field invariants. Violation of the general principles by the one-loop approximation in QED atmore » exponentially large magnetic field is analyzed, resulting in complex energy ghosts that signal the instability of the magnetized vacuum. Superluminal excitations (tachyons) appear, too, but for the magnetic field exceeding its instability threshold. Also other popular Lagrangians are tested to establish that the ones leading to spontaneous vacuum magnetization possess wrong convexity.« less

Development and evaluation of GRAL-C dispersion model, a hybrid Eulerian-Lagrangian approach capturing NO-NO 2-O 3 chemistry

NASA Astrophysics Data System (ADS)

Oettl, Dietmar; Uhrner, Ulrich

2011-02-01

Based on two recent publications using Lagrangian dispersion models to simulate NO-NO 2-O 3 chemistry for industrial plumes, a similar modified approach was implemented using GRAL-C ( Graz Lagrangian Model with Chemistry) and tested on two urban applications. In the hybrid dispersion model GRAL-C, the transport and turbulent diffusion of primary species such as NO and NO 2 are treated in a Lagrangian framework while those of O 3 are treated in an Eulerian framework. GRAL-C was employed on a one year street canyon simulation in Berlin and on a four-day simulation during a winter season in Graz, the second biggest city in Austria. In contrast to Middleton D.R., Jones A.R., Redington A.L., Thomson D.J., Sokhi R.S., Luhana L., Fisher B.E.A. (2008. Lagrangian modelling of plume chemistry for secondary pollutants in large industrial plumes. Atmospheric Environment 42, 415-427) and Alessandrini S., Ferrero E. (2008. A Lagrangian model with chemical reactions: application in real atmosphere. Proceedings of the 12th Int. Conf. on Harmonization within atmospheric dispersion modelling for regulatory purposes. Croatian Meteorological Journal, 43, ISSN: 1330-0083, 235-239) the treatment of ozone was modified in order to facilitate urban scale simulations encompassing dense road networks. For the street canyon application, modelled daily mean NO x/NO 2 concentrations deviated by +0.4%/-15% from observations, while the correlations for NO x and NO 2 were 0.67 and 0.76 respectively. NO 2 concentrations were underestimated in summer, but were captured well for other seasons. In Graz a fair agreement for NO x and NO 2 was obtained between observed and modelled values for NO x and NO 2. Simulated diurnal cycles of NO 2 and O 3 matched observations reasonably well, although O 3 was underestimated during the day. A possible explanation here might lie in the non-consideration of volatile organic compounds (VOCs) chemistry.

Generalized extended Lagrangian Born-Oppenheimer molecular dynamics

DOE PAGES

Niklasson, Anders M. N.; Cawkwell, Marc J.

2014-10-29

Extended Lagrangian Born-Oppenheimer molecular dynamics based on Kohn-Sham density functional theory is generalized in the limit of vanishing self-consistent field optimization prior to the force evaluations. The equations of motion are derived directly from the extended Lagrangian under the condition of an adiabatic separation between the nuclear and the electronic degrees of freedom. We show how this separation is automatically fulfilled and system independent. The generalized equations of motion require only one diagonalization per time step and are applicable to a broader range of materials with improved accuracy and stability compared to previous formulations.

Lagrangian solution of supersonic real gas flows

NASA Technical Reports Server (NTRS)

Loh, Ching-Yuen; Liou, Meng-Sing

1993-01-01

The present extention of a Lagrangian approach of the Riemann solution procedure, which was originally proposed for perfect gases, to real gases, is nontrivial and requires the development of an exact real-gas Riemann solver for the Lagrangian form of the conservation laws. Calculations including complex wave interactions of various types were conducted to test the accuracy and robustness of the approach. Attention is given to the case of 2D oblique waves' capture, where a slip line is clearly in evidence; the real gas effect is demonstrated in the case of a generic engine nozzle.

A Lagrangian Analysis of a Developing and Non-Developing Disturbance Observed During the PREDICT Experiment

DTIC Science & Technology

2012-12-03

paper provides an introduction of Lagrangian techniques for locating flow boundaries that encompass regions of recirculation in time- dependent flows...the low- to mid- level embryonic vortex from adverse conditions, while the 1The glossary on NOAA’s Hurricane Research Division’s web - site uses...wave or disturbance. This paper provides an introduction of Lagrangian techniques for locating flow boundaries that encompass regions of recirculation

Intermittent Lagrangian velocities and accelerations in three-dimensional porous medium flow.

PubMed

Holzner, M; Morales, V L; Willmann, M; Dentz, M

2015-07-01

Intermittency of Lagrangian velocity and acceleration is a key to understanding transport in complex systems ranging from fluid turbulence to flow in porous media. High-resolution optical particle tracking in a three-dimensional (3D) porous medium provides detailed 3D information on Lagrangian velocities and accelerations. We find sharp transitions close to pore throats, and low flow variability in the pore bodies, which gives rise to stretched exponential Lagrangian velocity and acceleration distributions characterized by a sharp peak at low velocity, superlinear evolution of particle dispersion, and double-peak behavior in the propagators. The velocity distribution is quantified in terms of pore geometry and flow connectivity, which forms the basis for a continuous-time random-walk model that sheds light on the observed Lagrangian flow and transport behaviors.

Elimination of artificial grid distortion and hourglass-type motions by means of Lagrangian subzonal masses and pressures

DOE Office of Scientific and Technical Information (OSTI.GOV)

Caramana, E.J.; Shashkov, M.J.

1997-12-31

The bane of Lagrangian hydrodynamics calculations is premature breakdown of the grid topology that results in severe degradation of accuracy and run termination often long before the assumption of Lagrangian zonal mass ceased to be valid. At short spatial grid scales this is usually referred to by the terms hourglass mode or keystone motion associated in particular with underconstrained grids such as quadrilaterals and hexahedrons in two and three dimensions, respectively. At longer spatial scales relative to the grid spacing there is what is referred to ubiquitously as spurious vorticity, or the long-thin zone problem. In both cases the resultmore » is anomalous grid distortion and tangling that has nothing to do with the actual solution, as would be the case for turbulent flow. In this work the authors show how such motions can be eliminated by the proper use of subzonal Lagrangian masses, and associated densities and pressures. These subzonal masses arise in a natural way from the fact that they require the mass associated with the nodal grid point to be constant in time. This is addition to the usual assumption of constant, Lagrangian zonal mass in staggered grid hydrodynamics scheme. The authors show that with proper discretization of subzonal forces resulting from subzonal pressures, hourglass motion and spurious vorticity can be eliminated for a very large range of problems. Finally the authors are presenting results of calculations of many test problems.« less

Lagrangian Approach to Study Catalytic Fluidized Bed Reactors

NASA Astrophysics Data System (ADS)

Madi, Hossein; Hossein Madi Team; Marcelo Kaufman Rechulski Collaboration; Christian Ludwig Collaboration; Tilman Schildhauer Collaboration

2013-03-01

Lagrangian approach of fluidized bed reactors is a method, which simulates the movement of catalyst particles (caused by the fluidization) by changing the gas composition around them. Application of such an investigation is in the analysis of the state of catalysts and surface reactions under quasi-operando conditions. The hydrodynamics of catalyst particles within a fluidized bed reactor was studied to improve a Lagrangian approach. A fluidized bed methanation employed in the production of Synthetic Natural Gas from wood was chosen as the case study. The Lagrangian perspective was modified and improved to include different particle circulation patterns, which were investigated through this study. Experiments were designed to evaluate the concepts of the model. The results indicate that the setup is able to perform the designed experiments and a good agreement between the simulation and the experimental results were observed. It has been shown that fluidized bed reactors, as opposed to fixed beds, can be used to avoid the deactivation of the methanation catalyst due to carbon deposits. Carbon deposition on the catalysts tested with the Lagrangian approach was investigated by temperature programmed oxidation (TPO) analysis of ex-situ catalyst samples. This investigation was done to identify the effects of particles velocity and their circulation patterns on the amount and type of deposited carbon on the catalyst surface. Ecole Polytechnique Federale de Lausanne(EPFL), Paul Scherrer Institute (PSI)

Lagrangian statistics and flow topology in forced two-dimensional turbulence.

PubMed

Kadoch, B; Del-Castillo-Negrete, D; Bos, W J T; Schneider, K

2011-03-01

A study of the relationship between Lagrangian statistics and flow topology in fluid turbulence is presented. The topology is characterized using the Weiss criterion, which provides a conceptually simple tool to partition the flow into topologically different regions: elliptic (vortex dominated), hyperbolic (deformation dominated), and intermediate (turbulent background). The flow corresponds to forced two-dimensional Navier-Stokes turbulence in doubly periodic and circular bounded domains, the latter with no-slip boundary conditions. In the double periodic domain, the probability density function (pdf) of the Weiss field exhibits a negative skewness consistent with the fact that in periodic domains the flow is dominated by coherent vortex structures. On the other hand, in the circular domain, the elliptic and hyperbolic regions seem to be statistically similar. We follow a Lagrangian approach and obtain the statistics by tracking large ensembles of passively advected tracers. The pdfs of residence time in the topologically different regions are computed introducing the Lagrangian Weiss field, i.e., the Weiss field computed along the particles' trajectories. In elliptic and hyperbolic regions, the pdfs of the residence time have self-similar algebraic decaying tails. In contrast, in the intermediate regions the pdf has exponential decaying tails. The conditional pdfs (with respect to the flow topology) of the Lagrangian velocity exhibit Gaussian-like behavior in the periodic and in the bounded domains. In contrast to the freely decaying turbulence case, the conditional pdfs of the Lagrangian acceleration in forced turbulence show a comparable level of intermittency in both the periodic and the bounded domains. The conditional pdfs of the Lagrangian curvature are characterized, in all cases, by self-similar power-law behavior with a decay exponent of order -2.

The Lagrangian-Hamiltonian formalism for higher order field theories

NASA Astrophysics Data System (ADS)

Vitagliano, Luca

2010-06-01

We generalize the Lagrangian-Hamiltonian formalism of Skinner and Rusk to higher order field theories on fiber bundles. As a byproduct we solve the long standing problem of defining, in a coordinate free manner, a Hamiltonian formalism for higher order Lagrangian field theories. Namely, our formalism does only depend on the action functional and, therefore, unlike previously proposed ones, is free from any relevant ambiguity.

Lagrangian and Hamiltonian constraints for guiding-center Hamiltonian theories

DOE Office of Scientific and Technical Information (OSTI.GOV)

Tronko, Natalia; Brizard, Alain J.

A consistent guiding-center Hamiltonian theory is derived by Lie-transform perturbation method, with terms up to second order in magnetic-field nonuniformity. Consistency is demonstrated by showing that the guiding-center transformation presented here satisfies separate Jacobian and Lagrangian constraints that have not been explored before. A new first-order term appearing in the guiding-center phase-space Lagrangian is identified through a calculation of the guiding-center polarization. It is shown that this new polarization term also yields a simpler expression of the guiding-center toroidal canonical momentum, which satisfies an exact conservation law in axisymmetric magnetic geometries. Finally, an application of the guiding-center Lagrangian constraint onmore » the guiding-center Hamiltonian yields a natural interpretation for its higher-order corrections.« less

Dissipative inertial transport patterns near coherent Lagrangian eddies in the ocean.

PubMed

Beron-Vera, Francisco J; Olascoaga, María J; Haller, George; Farazmand, Mohammad; Triñanes, Joaquín; Wang, Yan

2015-08-01

Recent developments in dynamical systems theory have revealed long-lived and coherent Lagrangian (i.e., material) eddies in incompressible, satellite-derived surface ocean velocity fields. Paradoxically, observed drifting buoys and floating matter tend to create dissipative-looking patterns near oceanic eddies, which appear to be inconsistent with the conservative fluid particle patterns created by coherent Lagrangian eddies. Here, we show that inclusion of inertial effects (i.e., those produced by the buoyancy and size finiteness of an object) in a rotating two-dimensional incompressible flow context resolves this paradox. Specifically, we obtain that anticyclonic coherent Lagrangian eddies attract (repel) negatively (positively) buoyant finite-size particles, while cyclonic coherent Lagrangian eddies attract (repel) positively (negatively) buoyant finite-size particles. We show how these results explain dissipative-looking satellite-tracked surface drifter and subsurface float trajectories, as well as satellite-derived Sargassum distributions.

Direct observations of L-I-H and H-I-L transitions with the X-point reciprocating probe in ASDEX Upgrade

DOE Office of Scientific and Technical Information (OSTI.GOV)

Müller, S. H.; Conway, G. D.; Birkenmeier, G.

A reciprocating Langmuir probe was used to directly measure the behavior of turbulence and flows in the X-point region during transitions between low-(L) and high-confinement (H) mode in ASDEX Upgrade. The probe traverses the divertor horizontally in 140 ms, typically 2–5 cm below the X-point. Toroidal Mach number, density, floating potential (ϕ{sub f}), and electron temperature (T{sub e}) are measured. In the regime accessible to the probe (P{sub inj}<1.5 MW, line-integrated core density <4×10{sup 19} m{sup −2}), the L-H transition features an intermediate phase (I-phase), characterized by limit-cycle oscillations at 0.5–3 kHz [Conway et al., Phys. Rev. Lett. 106, 065001 (2011)]. The probe measurements revealmore » that this pulsing affects both the density and the toroidal Mach number. It is present in both the low-(LFS) and high-field sides (HFS) of the scrape-off layer, while high-amplitude broadband turbulence usually dominates the private-flux region. Profile comparisons between L-mode and I-phase show lower density in pulsing regions and small shifts in T{sub e}, directed oppositely on LFS and HFS, which are compensated by shifts in ϕ{sub f} to yield a surprisingly unchanged plasma potential profile. Directly observed L-I-phase transitions reveal that the onset of the pulsing is preceded by a fast 50% density drop in the HFS X-point region. Back transitions to L-mode occur essentially symmetrically, with the pulsing stopping first, followed by a fast recovery to L-mode density levels in the divertor.« less

Deep Space Transportation System Using the Sun-Earth L2 Point

NASA Technical Reports Server (NTRS)

Matsumoto, Michihiro

2007-01-01

Recently, various kinds of planetary explorations have become more feasible, taking the advantage of low thrust propulsion means such as ion engines that have come into practical use. The field of space activity has now been expanded even to the rim of the outer solar system. In this context, the Japan Aerospace Exploration Agency (JAXA) has started investigating a Deep Space Port built at the L2 Lagrange point in the Sun-Earth system. For the purpose of making the deep space port practically useful, there is a need to establish a method to making spaceship depart and return from/to the port. This paper first discusses the escape maneuvers originating from the L2 point under the restricted three-body problem. Impulsive maneuvers from the L2 point are extensively studied here, and using the results, optimal low-thrust escape strategies are synthesized. Furthermore, this paper proposes the optimal escape and acceleration maneuvers schemes using Electric Delta-V Earth Gravity Assist (EDVEGA) technique.

Lagrangian Observations and Modeling of Marine Larvae

NASA Astrophysics Data System (ADS)

Paris, Claire B.; Irisson, Jean-Olivier

2017-04-01

Just within the past two decades, studies on the early-life history stages of marine organisms have led to new paradigms in population dynamics. Unlike passive plant seeds that are transported by the wind or by animals, marine larvae have motor and sensory capabilities. As a result, marine larvae have a tremendous capacity to actively influence their dispersal. This is continuously revealed as we develop new techniques to observe larvae in their natural environment and begin to understand their ability to detect cues throughout ontogeny, process the information, and use it to ride ocean currents and navigate their way back home, or to a place like home. We present innovative in situ and numerical modeling approaches developed to understand the underlying mechanisms of larval transport in the ocean. We describe a novel concept of a Lagrangian platform, the Drifting In Situ Chamber (DISC), designed to observe and quantify complex larval behaviors and their interactions with the pelagic environment. We give a brief history of larval ecology research with the DISC, showing that swimming is directional in most species, guided by cues as diverse as the position of the sun or the underwater soundscape, and even that (unlike humans!) larvae orient better and swim faster when moving as a group. The observed Lagrangian behavior of individual larvae are directly implemented in the Connectivity Modeling System (CMS), an open source Lagrangian tracking application. Simulations help demonstrate the impact that larval behavior has compared to passive Lagrangian trajectories. These methodologies are already the base of exciting findings and are promising tools for documenting and simulating the behavior of other small pelagic organisms, forecasting their migration in a changing ocean.

Markov Chain Monte Carlo from Lagrangian Dynamics.

PubMed

Lan, Shiwei; Stathopoulos, Vasileios; Shahbaba, Babak; Girolami, Mark

2015-04-01

Hamiltonian Monte Carlo (HMC) improves the computational e ciency of the Metropolis-Hastings algorithm by reducing its random walk behavior. Riemannian HMC (RHMC) further improves the performance of HMC by exploiting the geometric properties of the parameter space. However, the geometric integrator used for RHMC involves implicit equations that require fixed-point iterations. In some cases, the computational overhead for solving implicit equations undermines RHMC's benefits. In an attempt to circumvent this problem, we propose an explicit integrator that replaces the momentum variable in RHMC by velocity. We show that the resulting transformation is equivalent to transforming Riemannian Hamiltonian dynamics to Lagrangian dynamics. Experimental results suggests that our method improves RHMC's overall computational e ciency in the cases considered. All computer programs and data sets are available online (http://www.ics.uci.edu/~babaks/Site/Codes.html) in order to allow replication of the results reported in this paper.

Toroidal regularization of the guiding center Lagrangian

DOE PAGES

Burby, J. W.; Ellison, C. L.

2017-11-22

In the Lagrangian theory of guiding center motion, an effective magnetic field B* = B+ (m/e)v ∥∇ x b appears prominently in the equations of motion. Because the parallel component of this field can vanish, there is a range of parallel velocities where the Lagrangian guiding center equations of motion are either ill-defined or very badly behaved. Moreover, the velocity dependence of B* greatly complicates the identification of canonical variables and therefore the formulation of symplectic integrators for guiding center dynamics. Here, this letter introduces a simple coordinate transformation that alleviates both these problems simultaneously. In the new coordinates, themore » Liouville volume element is equal to the toroidal contravariant component of the magnetic field. Consequently, the large-velocity singularity is completely eliminated. Moreover, passing from the new coordinate system to canonical coordinates is extremely simple, even if the magnetic field is devoid of flux surfaces. We demonstrate the utility of this approach in regularizing the guiding center Lagrangian by presenting a new and stable one-step variational integrator for guiding centers moving in arbitrary time-dependent electromagnetic fields.« less

Toroidal regularization of the guiding center Lagrangian

DOE Office of Scientific and Technical Information (OSTI.GOV)

Burby, J. W.; Ellison, C. L.

In the Lagrangian theory of guiding center motion, an effective magnetic field B* = B+ (m/e)v ∥∇ x b appears prominently in the equations of motion. Because the parallel component of this field can vanish, there is a range of parallel velocities where the Lagrangian guiding center equations of motion are either ill-defined or very badly behaved. Moreover, the velocity dependence of B* greatly complicates the identification of canonical variables and therefore the formulation of symplectic integrators for guiding center dynamics. Here, this letter introduces a simple coordinate transformation that alleviates both these problems simultaneously. In the new coordinates, themore » Liouville volume element is equal to the toroidal contravariant component of the magnetic field. Consequently, the large-velocity singularity is completely eliminated. Moreover, passing from the new coordinate system to canonical coordinates is extremely simple, even if the magnetic field is devoid of flux surfaces. We demonstrate the utility of this approach in regularizing the guiding center Lagrangian by presenting a new and stable one-step variational integrator for guiding centers moving in arbitrary time-dependent electromagnetic fields.« less

Spectral-clustering approach to Lagrangian vortex detection.

PubMed

Hadjighasem, Alireza; Karrasch, Daniel; Teramoto, Hiroshi; Haller, George

2016-06-01

One of the ubiquitous features of real-life turbulent flows is the existence and persistence of coherent vortices. Here we show that such coherent vortices can be extracted as clusters of Lagrangian trajectories. We carry out the clustering on a weighted graph, with the weights measuring pairwise distances of fluid trajectories in the extended phase space of positions and time. We then extract coherent vortices from the graph using tools from spectral graph theory. Our method locates all coherent vortices in the flow simultaneously, thereby showing high potential for automated vortex tracking. We illustrate the performance of this technique by identifying coherent Lagrangian vortices in several two- and three-dimensional flows.

A view on coupled cluster perturbation theory using a bivariational Lagrangian formulation.

PubMed

Kristensen, Kasper; Eriksen, Janus J; Matthews, Devin A; Olsen, Jeppe; Jørgensen, Poul

2016-02-14

We consider two distinct coupled cluster (CC) perturbation series that both expand the difference between the energies of the CCSD (CC with single and double excitations) and CCSDT (CC with single, double, and triple excitations) models in orders of the Møller-Plesset fluctuation potential. We initially introduce the E-CCSD(T-n) series, in which the CCSD amplitude equations are satisfied at the expansion point, and compare it to the recently developed CCSD(T-n) series [J. J. Eriksen et al., J. Chem. Phys. 140, 064108 (2014)], in which not only the CCSD amplitude, but also the CCSD multiplier equations are satisfied at the expansion point. The computational scaling is similar for the two series, and both are term-wise size extensive with a formal convergence towards the CCSDT target energy. However, the two series are different, and the CCSD(T-n) series is found to exhibit a more rapid convergence up through the series, which we trace back to the fact that more information at the expansion point is utilized than for the E-CCSD(T-n) series. The present analysis can be generalized to any perturbation expansion representing the difference between a parent CC model and a higher-level target CC model. In general, we demonstrate that, whenever the parent parameters depend upon the perturbation operator, a perturbation expansion of the CC energy (where only parent amplitudes are used) differs from a perturbation expansion of the CC Lagrangian (where both parent amplitudes and parent multipliers are used). For the latter case, the bivariational Lagrangian formulation becomes more than a convenient mathematical tool, since it facilitates a different and faster convergent perturbation series than the simpler energy-based expansion.

Modeling of Multi-Kilovolt X-ray Driven Ablation and Closure of Pinholes during Point-Projection Backlit Imaging

NASA Astrophysics Data System (ADS)

Bullock, A. B.; Landen, O. L.; Bradley, D. K.

2000-10-01

Pinhole-assisted point-projection backlighting of large samples with few µm pinholes can result in pinhole closure due to x-ray driven ablation of the high Z pinhole substrate, thereby potentially limiting the usefulness of this imaging method. The results of a previous study[1] using streaked 1-D backlit imaging of 25 mm W wires at the OMEGA laser are compared to simulations produced by HYADES, a 1-D Lagrangian hydrodynamics code. Interestingly, the observed image resolution stays fixed while the pinhole transmission drops within 1-2 ns, suggesting rapid filling by a long scale-length of low density substrate material. These results will be compared to time-dependent HYADES predictions of pinhole closure timescales and resolution. 1 A.B. bullock, D.K. Bradley, and O.L. Landen, to be published in RSI (2001). *This work was performed under the auspices of the U.S. Department of Energy by University of California/Lawrence Livermore National Laboratory under Contract No. W-7405-Eng-48.

Gravity, Time, and Lagrangians

ERIC Educational Resources Information Center

Huggins, Elisha

2010-01-01

Feynman mentioned to us that he understood a topic in physics if he could explain it to a college freshman, a high school student, or a dinner guest. Here we will discuss two topics that took us a while to get to that level. One is the relationship between gravity and time. The other is the minus sign that appears in the Lagrangian. (Why would one…

Insights into the three-dimensional Lagrangian geometry of the Antarctic polar vortex

NASA Astrophysics Data System (ADS)

Curbelo, Jezabel; José García-Garrido, Víctor; Mechoso, Carlos Roberto; Mancho, Ana Maria; Wiggins, Stephen; Niang, Coumba

2017-07-01

In this paper we study the three-dimensional (3-D) Lagrangian structures in the stratospheric polar vortex (SPV) above Antarctica. We analyse and visualize these structures using Lagrangian descriptor function M. The procedure for calculation with reanalysis data is explained. Benchmarks are computed and analysed that allow us to compare 2-D and 3-D aspects of Lagrangian transport. Dynamical systems concepts appropriate to 3-D, such as normally hyperbolic invariant curves, are discussed and applied. In order to illustrate our approach we select an interval of time in which the SPV is relatively undisturbed (August 1979) and an interval of rapid SPV changes (October 1979). Our results provide new insights into the Lagrangian structure of the vertical extension of the stratospheric polar vortex and its evolution. Our results also show complex Lagrangian patterns indicative of strong mixing processes in the upper troposphere and lower stratosphere. Finally, during the transition to summer in the late spring, we illustrate the vertical structure of two counterrotating vortices, one the polar and the other an emerging one, and the invariant separatrix that divides them.

Lagrangian-averaged model for magnetohydrodynamic turbulence and the absence of bottlenecks.

PubMed

Pietarila Graham, Jonathan; Mininni, Pablo D; Pouquet, Annick

2009-07-01

We demonstrate that, for the case of quasiequipartition between the velocity and the magnetic field, the Lagrangian-averaged magnetohydrodynamics (LAMHD) alpha model reproduces well both the large-scale and the small-scale properties of turbulent flows; in particular, it displays no increased (superfilter) bottleneck effect with its ensuing enhanced energy spectrum at the onset of the subfilter scales. This is in contrast to the case of the neutral fluid in which the Lagrangian-averaged Navier-Stokes alpha model is somewhat limited in its applications because of the formation of spatial regions with no internal degrees of freedom and subsequent contamination of superfilter-scale spectral properties. We argue that, as the Lorentz force breaks the conservation of circulation and enables spectrally nonlocal energy transfer (associated with Alfvén waves), it is responsible for the absence of a viscous bottleneck in magnetohydrodynamics (MHD), as compared to the fluid case. As LAMHD preserves Alfvén waves and the circulation properties of MHD, there is also no (superfilter) bottleneck found in LAMHD, making this method capable of large reductions in required numerical degrees of freedom; specifically, we find a reduction factor of approximately 200 when compared to a direct numerical simulation on a large grid of 1536;{3} points at the same Reynolds number.

[L-THIA-based management design for controlling urban non-point source pollution].

PubMed

Guo, Qing-Hai; Yang, Liu; Ke-Ming, Ma

2007-11-01

L-THIA Model was used to simulate the amounts of NPS pollutants in 2 catchments of Sanjiao watershed (Sj1, Sj2) in Hanyang district, and the total simulated amount of NPS loads in Sj1 and Sj2 were 1.82 x 10(4) kg, 1.38 x 10(5) kg, respectively. Based on the theory of resource-sink" and interaction of pattern with process, a series of BMPs, including green roof, grassland, porous pavement, infiltration trench, vegetative filter strip and wet pond, were optimized, and effects of BMPs were simulated along the surface runoff pathway. The results show that total pollutants outputs entering Sj1 and Sj2 account for 14.65% and 6.57%, respectively. Combining L-THIA model and BMPs in series is a proper measure for non-point source pollution control and urban development planning at watershed or region scale.

Stochastic Lagrangian dynamics for charged flows in the E-F regions of ionosphere

DOE Office of Scientific and Technical Information (OSTI.GOV)

Tang Wenbo; Mahalov, Alex

2013-03-15

We develop a three-dimensional numerical model for the E-F region ionosphere and study the Lagrangian dynamics for plasma flows in this region. Our interest rests on the charge-neutral interactions and the statistics associated with stochastic Lagrangian motion. In particular, we examine the organizing mixing patterns for plasma flows due to polarized gravity wave excitations in the neutral field, using Lagrangian coherent structures (LCS). LCS objectively depict the flow topology-the extracted attractors indicate generation of ionospheric density gradients, due to accumulation of plasma. Using Lagrangian measures such as the finite-time Lyapunov exponents, we locate the Lagrangian skeletons for mixing in plasma,more » hence where charged fronts are expected to appear. With polarized neutral wind, we find that the corresponding plasma velocity is also polarized. Moreover, the polarized velocity alone, coupled with stochastic Lagrangian motion, may give rise to polarized density fronts in plasma. Statistics of these trajectories indicate high level of non-Gaussianity. This includes clear signatures of variance, skewness, and kurtosis of displacements taking polarized structures aligned with the gravity waves, and being anisotropic.« less

Observing Coronal Mass Ejections from the Sun-Earth L5 Point

NASA Astrophysics Data System (ADS)

Gopalswamy, N.; Davila, J. M.; St Cyr, O. C.

2013-12-01

Coronal mass ejections (CMEs) are the most energetic phenomenon in the heliosphere and are known to be responsible for severe space weather. Most of the current knowledge on CMEs accumulated over the past few decades has been derived from observations made from the Sun-Earth line, which is not the ideal vantage point to observe Earth-affecting CMEs (Gopalswamy et al., 2011a,b). The STEREO mission viewed CMEs from points away from the Sun-Earth line and demonstrated the importance of such observations in understanding the three-dimensional structure of CMEs and their true kinematics. In this paper, we show that it is advantageous to observe CMEs from the Sun-Earth L5 point in studying CMEs that affect Earth. In particular, these observations are important in identifying that part of the CME that is likely to arrive at Earth. L5 observations are critical for several aspects of CME studies such as: (i) they can also provide near-Sun space speed of CMEs, which is an important input for modeling Earth-arriving CMEs, (ii) backside and frontside CMEs can be readily distinguished even without inner coronal imagers, and (iii) preceding CMEs in the path of Earth-affecting CMEs can be identified for a better estimate of the travel time, which may not be possible from the Sun-Earth line. We also discuss how the L5 vantage point compares with the Sun-Earth L4 point for observing Earth-affecting CMEs. References Gopalswamy, N., Davila, J. M., St. Cyr, O. C., Sittler, E. C., Auchère, F., Duvall, T. L., Hoeksema, J. T., Maksimovic, M., MacDowall, R. J., Szabo, A., Collier, M. R. (2011a), Earth-Affecting Solar Causes Observatory (EASCO): A potential International Living with a Star Mission from Sun-Earth L5 JASTP 73, 658-663, DOI: 10.1016/j.jastp.2011.01.013 Gopalswamy, N., Davila, J. M., Auchère, F., Schou, J., Korendyke, C. M. Shih, A., Johnston, J. C., MacDowall, R. J., Maksimovic, M., Sittler, E., et al. (2011b), Earth-Affecting Solar Causes Observatory (EASCO): a mission at

On the Lagrangian description of unsteady boundary layer separation. Part 1: General theory

NASA Technical Reports Server (NTRS)

Vandommelen, Leon L.; Cowley, Stephen J.

1989-01-01

Although unsteady, high-Reynolds number, laminar boundary layers have conventionally been studied in terms of Eulerian coordinates, a Lagrangian approach may have significant analytical and computational advantages. In Lagrangian coordinates the classical boundary layer equations decouple into a momentum equation for the motion parallel to the boundary, and a hyperbolic continuity equation (essentially a conserved Jacobian) for the motion normal to the boundary. The momentum equations, plus the energy equation if the flow is compressible, can be solved independently of the continuity equation. Unsteady separation occurs when the continuity equation becomes singular as a result of touching characteristics, the condition for which can be expressed in terms of the solution of the momentum equations. The solutions to the momentum and energy equations remain regular. Asymptotic structures for a number of unsteady 3-D separating flows follow and depend on the symmetry properties of the flow. In the absence of any symmetry, the singularity structure just prior to separation is found to be quasi 2-D with a displacement thickness in the form of a crescent shaped ridge. Physically the singularities can be understood in terms of the behavior of a fluid element inside the boundary layer which contracts in a direction parallel to the boundary and expands normal to it, thus forcing the fluid above it to be ejected from the boundary layer.

Bifurcations and degenerate periodic points in a three dimensional chaotic fluid flow

DOE Office of Scientific and Technical Information (OSTI.GOV)

Smith, L. D., E-mail: lachlan.smith@monash.edu; CSIRO Mineral Resources, Clayton, Victoria 3800; Rudman, M.

2016-05-15

Analysis of the periodic points of a conservative periodic dynamical system uncovers the basic kinematic structure of the transport dynamics and identifies regions of local stability or chaos. While elliptic and hyperbolic points typically govern such behaviour in 3D systems, degenerate (parabolic) points also play an important role. These points represent a bifurcation in local stability and Lagrangian topology. In this study, we consider the ramifications of the two types of degenerate periodic points that occur in a model 3D fluid flow. (1) Period-tripling bifurcations occur when the local rotation angle associated with elliptic points is reversed, creating a reversalmore » in the orientation of associated Lagrangian structures. Even though a single unstable point is created, the bifurcation in local stability has a large influence on local transport and the global arrangement of manifolds as the unstable degenerate point has three stable and three unstable directions, similar to hyperbolic points, and occurs at the intersection of three hyperbolic periodic lines. The presence of period-tripling bifurcation points indicates regions of both chaos and confinement, with the extent of each depending on the nature of the associated manifold intersections. (2) The second type of bifurcation occurs when periodic lines become tangent to local or global invariant surfaces. This bifurcation creates both saddle–centre bifurcations which can create both chaotic and stable regions, and period-doubling bifurcations which are a common route to chaos in 2D systems. We provide conditions for the occurrence of these tangent bifurcations in 3D conservative systems, as well as constraints on the possible types of tangent bifurcation that can occur based on topological considerations.« less

Lagrangian trajectories, residual currents and rectification process in the Northern Gulf of California

NASA Astrophysics Data System (ADS)

Rodríguez, Pablo Alonso; Carbajal, Noel; Rodríguez, Juan Heberto Gaviño

2017-07-01

Considering a semi-implicit approximation of the Coriolis terms, a numerical solution of the vertically integrated equations of motion is proposed. To test the two-dimensional numerical model, several experiments for the calculation of Euler, Stokes and Lagrange residual currents in the Gulf of California were carried out. To estimate the Lagrangian residual current, trajectories of particles were also simulated. The applied tidal constituents were M2, S2, K2, N2, K1, P1 and O1. At spring tides, strong tidal velocities occur in the northern half of the gulf. In this region of complex geometry, depths change from a few meter in the northern shelf zone to more than 3000 m in the southern part. In the archipelago region, the presence of islands alters amplitude and direction of tidal currents producing a rectification process which is reflected in a clockwise circulation around Tiburón Island in the Lagrangian residual current. The rectification process is explained by the superposition of the Euler and Stokes residual currents. Residual current patterns show several cyclonic and anticyclonic gyres in the Northern Gulf of California. Numerical experiments for individual and combinations of several tidal constituents revealed a large variability of Lagrangian trajectories.

Lagrangians and Euler morphisms from connections on the frame bundle

NASA Astrophysics Data System (ADS)

Kurek, J.; Mikulski, W. M.

2011-07-01

We classify all natural operators transforming torsion free classical linear connections ∇ on m-dimensional manifolds M into r-th order Lagrangians λ(∇) and Euler morphisms E(∇) on the linear frame bundle P1M. We also briefly write how this classification result can be generalized on higher order frame bundles PkM instead of P1M.

Asymptotic-preserving Lagrangian approach for modeling anisotropic transport in magnetized plasmas

NASA Astrophysics Data System (ADS)

Chacon, Luis; Del-Castillo-Negrete, Diego

2011-10-01

Modeling electron transport in magnetized plasmas is extremely challenging due to the extreme anisotropy introduced by the presence of the magnetic field (χ∥ /χ⊥ ~1010 in fusion plasmas). Recently, a novel Lagrangian method has been proposed to solve the local and non-local purely parallel transport equation in general 3D magnetic fields. The approach avoids numerical pollution (in fact, it respects transport barriers -flux surfaces- exactly by construction), is inherently positivity-preserving, and is scalable algorithmically (i.e., work per degree-of-freedom is grid-independent). In this poster, we discuss the extension of the Lagrangian approach to include perpendicular transport and sources. We present an asymptotic-preserving numerical formulation that ensures a consistent numerical discretization temporally and spatially for arbitrary χ∥ /χ⊥ ratios. This is of importance because parallel and perpendicular transport terms in the transport equation may become comparable in regions of the plasma (e.g., at incipient islands), while remaining disparate elsewhere. We will demonstrate the potential of the approach with various challenging configurations, including the case of transport across a magnetic island in cylindrical geometry. D. del-Castillo-Negrete, L. Chacón, PRL, 106, 195004 (2011); DPP11 invited talk by del-Castillo-Negrete.

Residual sweeping errors in turbulent particle pair diffusion in a Lagrangian diffusion model.

PubMed

Malik, Nadeem A

2017-01-01

Thomson, D. J. & Devenish, B. J. [J. Fluid Mech. 526, 277 (2005)] and others have suggested that sweeping effects make Lagrangian properties in Kinematic Simulations (KS), Fung et al [Fung J. C. H., Hunt J. C. R., Malik N. A. & Perkins R. J. J. Fluid Mech. 236, 281 (1992)], unreliable. However, such a conclusion can only be drawn under the assumption of locality. The major aim here is to quantify the sweeping errors in KS without assuming locality. Through a novel analysis based upon analysing pairs of particle trajectories in a frame of reference moving with the large energy containing scales of motion it is shown that the normalized integrated error [Formula: see text] in the turbulent pair diffusivity (K) due to the sweeping effect decreases with increasing pair separation (σl), such that [Formula: see text] as σl/η → ∞; and [Formula: see text] as σl/η → 0. η is the Kolmogorov turbulence microscale. There is an intermediate range of separations 1 < σl/η < ∞ in which the error [Formula: see text] remains negligible. Simulations using KS shows that in the swept frame of reference, this intermediate range is large covering almost the entire inertial subrange simulated, 1 < σl/η < 105, implying that the deviation from locality observed in KS cannot be atributed to sweeping errors. This is important for pair diffusion theory and modeling. PACS numbers: 47.27.E?, 47.27.Gs, 47.27.jv, 47.27.Ak, 47.27.tb, 47.27.eb, 47.11.-j.

Two-point correlation function for Dirichlet L-functions

NASA Astrophysics Data System (ADS)

Bogomolny, E.; Keating, J. P.

2013-03-01

The two-point correlation function for the zeros of Dirichlet L-functions at a height E on the critical line is calculated heuristically using a generalization of the Hardy-Littlewood conjecture for pairs of primes in arithmetic progression. The result matches the conjectured random-matrix form in the limit as E → ∞ and, importantly, includes finite-E corrections. These finite-E corrections differ from those in the case of the Riemann zeta-function, obtained in Bogomolny and Keating (1996 Phys. Rev. Lett. 77 1472), by certain finite products of primes which divide the modulus of the primitive character used to construct the L-function in question.

Users manual for a one-dimensional Lagrangian transport model

USGS Publications Warehouse

Schoellhamer, D.H.; Jobson, H.E.

1986-01-01

A Users Manual for the Lagrangian Transport Model (LTM) is presented. The LTM uses Lagrangian calculations that are based on a reference frame moving with the river flow. The Lagrangian reference frame eliminates the need to numerically solve the convective term of the convection-diffusion equation and provides significant numerical advantages over the more commonly used Eulerian reference frame. When properly applied, the LTM can simulate riverine transport and decay processes within the accuracy required by most water quality studies. The LTM is applicable to steady or unsteady one-dimensional unidirectional flows in fixed channels with tributary and lateral inflows. Application of the LTM is relatively simple and optional capabilities improve the model 's convenience. Appendices give file formats and three example LTM applications that include the incorporation of the QUAL II water quality model 's reaction kinetics into the LTM. (Author 's abstract)

Uncovering the Geometry of Barrierless Reactions Using Lagrangian Descriptors.

PubMed

Junginger, Andrej; Hernandez, Rigoberto

2016-03-03

Transition-state theories describing barrierless chemical reactions, or more general activated problems, are often hampered by the lack of a saddle around which the dividing surface can be constructed. For example, the time-dependent transition-state trajectory uncovering the nonrecrossing dividing surface in thermal reactions in the framework of the Langevin equation has relied on perturbative approaches in the vicinity of the saddle. We recently obtained an alternative approach using Lagrangian descriptors to construct time-dependent and recrossing-free dividing surfaces. This is a nonperturbative approach making no reference to a putative saddle. Here we show how the Lagrangian descriptor can be used to obtain the transition-state geometry of a dissipated and thermalized reaction across barrierless potentials. We illustrate the method in the case of a 1D Brownian motion for both barrierless and step potentials; however, the method is not restricted and can be directly applied to different kinds of potentials and higher dimensional systems.

Blow-up solutions for L 2 supercritical gKdV equations with exactly k blow-up points

NASA Astrophysics Data System (ADS)

Lan, Yang

2017-08-01

In this paper we consider the slightly L 2-supercritical gKdV equations \\partialt u+(uxx+u\\vert u\\vert p-1)_x=0 , with the nonlinearity 5 and 0<\\varepsilon\\ll 1 . In the previous work of the author, we know that there exists a stable self-similar blow-up dynamics for slightly L 2-supercritical gKdV equations. Such solutions can be viewed as solutions with a single blow-up point. In this paper we will prove the existence of solutions with multiple blow-up points, and give a description of the formation of the singularity near the blow-up time.

Design of the stabilizing control of the orbital motion in the vicinity of the collinear libration point L1 using the analytical representation of the invariant manifold

NASA Astrophysics Data System (ADS)

Maliavkin, G. P.; Shmyrov, A. S.; Shmyrov, V. A.

2018-05-01

Vicinities of collinear libration points of the Sun-Earth system are currently quite attractive for the space navigation. Today, various projects on placing of spacecrafts observing the Sun in the L1 libration point and telescopes in L2 have been implemented (e.g. spacecrafts "WIND", "SOHO", "Herschel", "Planck"). Collinear libration points being unstable leads to the problem of stabilization of a spacecraft's motion. Laws of stabilizing motion control in vicinity of L1 point can be constructed using the analytical representation of a stable invariant manifold. Efficiency of these control laws depends on the precision of the representation. Within the model of Hill's approximation of the circular restricted three-body problem in the rotating geocentric coordinate system one can obtain the analytical representation of an invariant manifold filled with bounded trajectories in a form of series in terms of powers of the phase variables. Approximate representations of the orders from the first to the fourth inclusive can be used to construct four laws of stabilizing feedback motion control under which trajectories approach the manifold. By virtue of numerical simulation the comparison can be made: how the precision of the representation of the invariant manifold influences the efficiency of the control, expressed by energy consumptions (characteristic velocity). It shows that using approximations of higher orders in constructing the control laws can significantly reduce the energy consumptions on implementing the control compared to the linear approximation.

Immune Checkpoint PD-1/PD-L1: Is There Life Beyond Antibodies?

PubMed

Konstantinidou, Markella; Zarganes-Tzitzikas, Tryfon; Magiera-Mularz, Katarzyna; Holak, Tad A; Dömling, Alexander

2018-04-23

The PD-1/PD-L1 interaction has emerged as a significant target in cancer immunotherapy. Current medications include monoclonal antibodies, which have shown impressive clinical results in the treatment of several types of tumors. The cocrystal structure of human PD-1 and PD-L1 is expected to be a valuable starting point for the design of novel inhibitors, along with the recent crystal structures with monoclonal antibodies, small molecules, and macrocycles. © 2018 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.

Lagrangian and Eulerian statistics obtained from direct numerical simulations of homogeneous turbulence

NASA Technical Reports Server (NTRS)

Squires, Kyle D.; Eaton, John K.

1991-01-01

Direct numerical simulation is used to study dispersion in decaying isotropic turbulence and homogeneous shear flow. Both Lagrangian and Eulerian data are presented allowing direct comparison, but at fairly low Reynolds number. The quantities presented include properties of the dispersion tensor, isoprobability contours of particle displacement, Lagrangian and Eulerian velocity autocorrelations and time scale ratios, and the eddy diffusivity tensor. The Lagrangian time microscale is found to be consistently larger than the Eulerian microscale, presumably due to the advection of the small scales by the large scales in the Eulerian reference frame.

Dirac and Pauli form factors of nucleons using nonlocal chiral effective Lagrangian

NASA Astrophysics Data System (ADS)

He, Fangcheng; Wang, Ping

2017-11-01

Dirac and Pauli form factors are investigated in the relativistic chiral effective Lagrangian. The octet and decuplet intermediate states are included in the one-loop calculation. The 4-dimensional regulator is introduced to deal with the divergence. Different from the non-relativistic case, this 4-dimensional regulator is generated from the nonlocal Lagrangian with the gauge link, which guarantees local gauge invariance. As a result, additional diagrams appear which ensure electric charge 1 and 0 for proton and neutron respectively. The obtained Dirac and Pauli form factors of the nucleons are all reasonable up to relatively large Q 2. Supported by National Natural Science Foundation of China (11475186) and Sino-German CRC 110 (NSFC 11621131001)

A contemporary look at Hermann Hankel's 1861 pioneering work on Lagrangian fluid dynamics

NASA Astrophysics Data System (ADS)

Frisch, Uriel; Grimberg, Gérard; Villone, Barbara

2017-12-01

The present paper is a companion to the paper by Villone and Rampf (2017), titled "Hermann Hankel's On the general theory of motion of fluids, an essay including an English translation of the complete Preisschrift from 1861" together with connected documents [Eur. Phys. J. H 42, 557-609 (2017)]. Here we give a critical assessment of Hankel's work, which covers many important aspects of fluid dynamics considered from a Lagrangian-coordinates point of view: variational formulation in the spirit of Hamilton for elastic (barotropic) fluids, transport (we would now say Lie transport) of vorticity, the Lagrangian significance of Clebsch variables, etc. Hankel's work is also put in the perspective of previous and future work. Hence, the action spans about two centuries: from Lagrange's 1760-1761 Turin paper on variational approaches to mechanics and fluid mechanics problems to Arnold's 1966 founding paper on the geometrical/variational formulation of incompressible flow. The 22-year-old Hankel - who was to die 12 years later — emerges as a highly innovative master of mathematical fluid dynamics, fully deserving Riemann's assessment that his Preisschrift contains "all manner of good things."

Incomplete augmented Lagrangian preconditioner for steady incompressible Navier-Stokes equations.

PubMed

Tan, Ning-Bo; Huang, Ting-Zhu; Hu, Ze-Jun

2013-01-01

An incomplete augmented Lagrangian preconditioner, for the steady incompressible Navier-Stokes equations discretized by stable finite elements, is proposed. The eigenvalues of the preconditioned matrix are analyzed. Numerical experiments show that the incomplete augmented Lagrangian-based preconditioner proposed is very robust and performs quite well by the Picard linearization or the Newton linearization over a wide range of values of the viscosity on both uniform and stretched grids.

Incomplete Augmented Lagrangian Preconditioner for Steady Incompressible Navier-Stokes Equations

PubMed Central

Tan, Ning-Bo; Huang, Ting-Zhu; Hu, Ze-Jun

2013-01-01

An incomplete augmented Lagrangian preconditioner, for the steady incompressible Navier-Stokes equations discretized by stable finite elements, is proposed. The eigenvalues of the preconditioned matrix are analyzed. Numerical experiments show that the incomplete augmented Lagrangian-based preconditioner proposed is very robust and performs quite well by the Picard linearization or the Newton linearization over a wide range of values of the viscosity on both uniform and stretched grids. PMID:24235888

Programmers manual for a one-dimensional Lagrangian transport model

USGS Publications Warehouse

Schoellhamer, D.H.; Jobson, H.E.

1986-01-01

A one-dimensional Lagrangian transport model for simulating water-quality constituents such as temperature, dissolved oxygen , and suspended sediment in rivers is presented in this Programmers Manual. Lagrangian transport modeling techniques, the model 's subroutines, and the user-written decay-coefficient subroutine are discussed in detail. Appendices list the program codes. The Programmers Manual is intended for the model user who needs to modify code either to adapt the model to a particular need or to use reaction kinetics not provided with the model. (Author 's abstract)

Lorentz Invariance of Gravitational Lagrangians in the Space of Reference Frames

NASA Astrophysics Data System (ADS)

Cognola, G.

1980-06-01

The recently proposed theories of gravitation in the space of reference frames S are based on a Lagrangian invariant with respect to the homogeneous Lorentz group. However, in theories of this kind, the Lorentz invariance is not a necessary consequence of some physical principles, as in the theories formulated in space-time, but rather a purely esthetic request. In the present paper, we give a systematic method for the construction of gravitational theories in the space S, without assuming a priori the Lorentz invariance of the Lagrangian. The Einstein-Cartan equations of gravitation are obtained requiring only that the Lagrangian is invariant under proper rotations and has particular transformation properties under space reflections and space-time dilatations

Critical N = (1, 1) general massive supergravity

NASA Astrophysics Data System (ADS)

Deger, Nihat Sadik; Moutsopoulos, George; Rosseel, Jan

2018-04-01

In this paper we study the supermultiplet structure of N = (1, 1) General Massive Supergravity at non-critical and critical points of its parameter space. To do this, we first linearize the theory around its maximally supersymmetric AdS3 vacuum and obtain the full linearized Lagrangian including fermionic terms. At generic values, linearized modes can be organized as two massless and 2 massive multiplets where supersymmetry relates them in the standard way. At critical points logarithmic modes appear and we find that in three of such points some of the supersymmetry transformations are non-invertible in logarithmic multiplets. However, in the fourth critical point, there is a massive logarithmic multiplet with invertible supersymmetry transformations.

PD-L1 gene polymorphisms and low serum level of PD-L1 protein are associated to type 1 diabetes in Chile.

PubMed

Pizarro, Carolina; García-Díaz, Diego F; Codner, Ethel; Salas-Pérez, Francisca; Carrasco, Elena; Pérez-Bravo, Francisco

2014-11-01

Type 1 diabetes (T1D) has a complex etiology in which genetic and environmental factors are involved, whose interactions have not yet been completely clarified. In this context, the role in PD-1 pathway and its ligands 1 and 2 (PD-L1 and PD-L2) have been proposed as candidates in several autoimmune diseases. The aim of this work was to determine the allele and haplotype frequency of six gene polymorphisms of PD-ligands (PD-L1 and PD-L2) in Chilean T1D patients and their effect on serum levels of PD-L1 and autoantibody profile (GAD65 and IA2). This study cohort comprised 205 T1D patients and 205 normal children. We performed genotypic analysis of PD-L1 and PD-L2 genes by TaqMan method. Determination of anti-GAD65 and anti-IA-2 autoantibodies was performed by ELISA. The PD-L1 serum levels were measured. The allelic distribution of PD-L1 variants (rs2297137 and rs4143815) showed differences between T1D patients and controls (p = 0.035 and p = 0.022, respectively). No differences were detected among the PD-L2 polymorphisms, and only the rs16923189 showed genetic variation. T1D patients showed decreased serum levels of PD-L1 compared to controls: 1.42 [0.23-7.45] ng/mL versus 3.35 [0.49-5.89] ng/mL (p < 0.025). In addition, the CGG haplotype in PD-L1 associated with T1D (constructed from rs822342, rs2297137 and rs4143815 polymorphisms) showed an OR = 1.44 [1.08 to 1.93]. Finally, no association of these genetic variants was observed with serum concentrations of PD ligands or auto-antibody profile, although a correlation between PD-L1 ligand serum concentration and the age at disease onset was detected. Two polymorphism of PD-L1 are presented in different allelic variants between T1D and healthy subjects, also PDL-1 serum levels are significantly lowered in diabetics patients. Moreover, the age of onset of the disease determine differences between serum ligand levels in diabetics, being lower in younger. These results points to a possible establishment of

Deconstructing field-induced ketene isomerization through Lagrangian descriptors.

PubMed

Craven, Galen T; Hernandez, Rigoberto

2016-02-07

The time-dependent geometrical separatrices governing state transitions in field-induced ketene isomerization are constructed using the method of Lagrangian descriptors. We obtain the stable and unstable manifolds of time-varying transition states as dynamic phase space objects governing configurational changes when the ketene molecule is subjected to an oscillating electric field. The dynamics of the isomerization reaction are modeled through classical trajectory studies on the Gezelter-Miller potential energy surface and an approximate dipole moment model which is coupled to a time-dependent electric field. We obtain a representation of the reaction geometry, over varying field strengths and oscillation frequencies, by partitioning an initial phase space into basins labeled according to which product state is reached at a given time. The borders between these basins are in agreement with those obtained using Lagrangian descriptors, even in regimes exhibiting chaotic dynamics. Major outcomes of this work are: validation and extension of a transition state theory framework built from Lagrangian descriptors, elaboration of the applicability for this theory to periodically- and aperiodically-driven molecular systems, and prediction of regimes in which isomerization of ketene and its derivatives may be controlled using an external field.

3 Lectures: "Lagrangian Models", "Numerical Transport Schemes", and "Chemical and Transport Models"

NASA Technical Reports Server (NTRS)

Douglass, A.

2005-01-01

The topics for the three lectures for the Canadian Summer School are Lagrangian Models, numerical transport schemes, and chemical and transport models. In the first lecture I will explain the basic components of the Lagrangian model (a trajectory code and a photochemical code), the difficulties in using such a model (initialization) and show some applications in interpretation of aircraft and satellite data. If time permits I will show some results concerning inverse modeling which is being used to evaluate sources of tropospheric pollutants. In the second lecture I will discuss one of the core components of any grid point model, the numerical transport scheme. I will explain the basics of shock capturing schemes, and performance criteria. I will include an example of the importance of horizontal resolution to polar processes. We have learned from NASA's global modeling initiative that horizontal resolution matters for predictions of the future evolution of the ozone hole. The numerical scheme will be evaluated using performance metrics based on satellite observations of long-lived tracers. The final lecture will discuss the evolution of chemical transport models over the last decade. Some of the problems with assimilated winds will be demonstrated, using satellite data to evaluate the simulations.

Vorticity-divergence semi-Lagrangian global atmospheric model SL-AV20: dynamical core

NASA Astrophysics Data System (ADS)

Tolstykh, Mikhail; Shashkin, Vladimir; Fadeev, Rostislav; Goyman, Gordey

2017-05-01

SL-AV (semi-Lagrangian, based on the absolute vorticity equation) is a global hydrostatic atmospheric model. Its latest version, SL-AV20, provides global operational medium-range weather forecast with 20 km resolution over Russia. The lower-resolution configurations of SL-AV20 are being tested for seasonal prediction and climate modeling. The article presents the model dynamical core. Its main features are a vorticity-divergence formulation at the unstaggered grid, high-order finite-difference approximations, semi-Lagrangian semi-implicit discretization and the reduced latitude-longitude grid with variable resolution in latitude. The accuracy of SL-AV20 numerical solutions using a reduced lat-lon grid and the variable resolution in latitude is tested with two idealized test cases. Accuracy and stability of SL-AV20 in the presence of the orography forcing are tested using the mountain-induced Rossby wave test case. The results of all three tests are in good agreement with other published model solutions. It is shown that the use of the reduced grid does not significantly affect the accuracy up to the 25 % reduction in the number of grid points with respect to the regular grid. Variable resolution in latitude allows us to improve the accuracy of a solution in the region of interest.

Identifying Lagrangian fronts with favourable fishery conditions

NASA Astrophysics Data System (ADS)

Prants, S. V.; Budyansky, M. V.; Uleysky, M. Yu.

2014-08-01

Lagrangian fronts (LFs) in the ocean are defined as boundaries between surface waters with strongly different Lagrangian properties. They can be accurately detected in a given velocity field by computing synoptic maps for displacements of synthetic tracers and other Lagrangian indicators. We use Pacific saury catch and location data for a number of commercial fishery seasons in the region of the northwest Pacific with one of the richest fishery in the world. It is shown statistically that the saury fishing grounds with maximal catches are not randomly distributed over the region but located mainly along the sharp LFs where productive cold waters of the Oyashio Current, warmer waters of the southern branch of the Soya Current, and waters of warm-core Kuroshio rings converge. Computation of those fronts in altimetric geostrophic velocity fields both in the years with the First and Second Oyashio Intrusions shows that in spite of different oceanographic conditions LF locations may serve as good indicators of potential fishing grounds. Possible biophysical reasons for saury aggregation near sharp LFs are discussed. We propose a mechanism for effective export of nutrient rich waters based on stretching of material lines in the vicinity of hyperbolic objects in the ocean. The developed method, based on identifying LFs in any velocity fields, is quite general and may be applied to find potential fishing grounds for the other pelagic fish.

Eulerian Lagrangian Adaptive Fup Collocation Method for solving the conservative solute transport in heterogeneous porous media

NASA Astrophysics Data System (ADS)

Gotovac, Hrvoje; Srzic, Veljko

2014-05-01

Contaminant transport in natural aquifers is a complex, multiscale process that is frequently studied using different Eulerian, Lagrangian and hybrid numerical methods. Conservative solute transport is typically modeled using the advection-dispersion equation (ADE). Despite the large number of available numerical methods that have been developed to solve it, the accurate numerical solution of the ADE still presents formidable challenges. In particular, current numerical solutions of multidimensional advection-dominated transport in non-uniform velocity fields are affected by one or all of the following problems: numerical dispersion that introduces artificial mixing and dilution, grid orientation effects, unresolved spatial and temporal scales and unphysical numerical oscillations (e.g., Herrera et al, 2009; Bosso et al., 2012). In this work we will present Eulerian Lagrangian Adaptive Fup Collocation Method (ELAFCM) based on Fup basis functions and collocation approach for spatial approximation and explicit stabilized Runge-Kutta-Chebyshev temporal integration (public domain routine SERK2) which is especially well suited for stiff parabolic problems. Spatial adaptive strategy is based on Fup basis functions which are closely related to the wavelets and splines so that they are also compactly supported basis functions; they exactly describe algebraic polynomials and enable a multiresolution adaptive analysis (MRA). MRA is here performed via Fup Collocation Transform (FCT) so that at each time step concentration solution is decomposed using only a few significant Fup basis functions on adaptive collocation grid with appropriate scales (frequencies) and locations, a desired level of accuracy and a near minimum computational cost. FCT adds more collocations points and higher resolution levels only in sensitive zones with sharp concentration gradients, fronts and/or narrow transition zones. According to the our recent achievements there is no need for solving the large

A coupled Eulerian/Lagrangian method for the solution of three-dimensional vortical flows

NASA Technical Reports Server (NTRS)

Felici, Helene Marie

1992-01-01

A coupled Eulerian/Lagrangian method is presented for the reduction of numerical diffusion observed in solutions of three-dimensional rotational flows using standard Eulerian finite-volume time-marching procedures. A Lagrangian particle tracking method using particle markers is added to the Eulerian time-marching procedure and provides a correction of the Eulerian solution. In turn, the Eulerian solutions is used to integrate the Lagrangian state-vector along the particles trajectories. The Lagrangian correction technique does not require any a-priori information on the structure or position of the vortical regions. While the Eulerian solution ensures the conservation of mass and sets the pressure field, the particle markers, used as 'accuracy boosters,' take advantage of the accurate convection description of the Lagrangian solution and enhance the vorticity and entropy capturing capabilities of standard Eulerian finite-volume methods. The combined solution procedures is tested in several applications. The convection of a Lamb vortex in a straight channel is used as an unsteady compressible flow preservation test case. The other test cases concern steady incompressible flow calculations and include the preservation of turbulent inlet velocity profile, the swirling flow in a pipe, and the constant stagnation pressure flow and secondary flow calculations in bends. The last application deals with the external flow past a wing with emphasis on the trailing vortex solution. The improvement due to the addition of the Lagrangian correction technique is measured by comparison with analytical solutions when available or with Eulerian solutions on finer grids. The use of the combined Eulerian/Lagrangian scheme results in substantially lower grid resolution requirements than the standard Eulerian scheme for a given solution accuracy.

Chiral Lagrangian with Heavy Quark-Diquark Symmetry

DOE Office of Scientific and Technical Information (OSTI.GOV)

Jie Hu; Thomas Mehen

2005-11-29

We construct a chiral Lagrangian for doubly heavy baryons and heavy mesons that is invariant under heavy quark-diquark symmetry at leading order and includes the leading O(1/m{sub Q}) symmetry violating operators. The theory is used to predict the electromagnetic decay width of the J=3/2 member of the ground state doubly heavy baryon doublet. Numerical estimates are provided for doubly charm baryons. We also calculate chiral corrections to doubly heavy baryon masses and strong decay widths of low lying excited doubly heavy baryons.

Vocabulary Development in Norwegian L1 and L2 Learners in the Kindergarten-School Transition

ERIC Educational Resources Information Center

Karlsen, Jannicke; Lyster, Solveig-Alma Halaas; Lervåg, Arne

2017-01-01

This study examined the vocabulary development of Norwegian second language (L2) learners with Urdu/Punjabi as their first language (L1) at two time-points from kindergarten to primary school, and compared it to the vocabulary development of monolingual Norwegian children. Using path models, the associations between number of picture books in the…

Symmetries of SU(2) Skyrmion in Hamiltonian and Lagrangian Approaches

NASA Astrophysics Data System (ADS)

Hong, Soon-Tae; Kim, Yong-Wan; Park, Young-Jai

We apply the Batalin-Fradkin-Tyutin (BFT) method to the SU(2) Skyrmion to study the full symmetry structure of the model at the first-class Hamiltonian level. On the other hand, we also analyze the symmetry structure of the action having the WZ term, which corresponds to this Hamiltonian, in the framework of the Lagrangian approach. Furthermore, following the BFV formalism we derive the BRST invariant gauge fixed Lagrangian from the above extended action.

Seakeeping with the semi-Lagrangian particle finite element method

NASA Astrophysics Data System (ADS)

Nadukandi, Prashanth; Servan-Camas, Borja; Becker, Pablo Agustín; Garcia-Espinosa, Julio

2017-07-01

The application of the semi-Lagrangian particle finite element method (SL-PFEM) for the seakeeping simulation of the wave adaptive modular vehicle under spray generating conditions is presented. The time integration of the Lagrangian advection is done using the explicit integration of the velocity and acceleration along the streamlines (X-IVAS). Despite the suitability of the SL-PFEM for the considered seakeeping application, small time steps were needed in the X-IVAS scheme to control the solution accuracy. A preliminary proposal to overcome this limitation of the X-IVAS scheme for seakeeping simulations is presented.

A new method to calibrate Lagrangian model with ASAR images for oil slick trajectory.

PubMed

Tian, Siyu; Huang, Xiaoxia; Li, Hongga

2017-03-15

Since Lagrangian model coefficients vary with different conditions, it is necessary to calibrate the model to obtain optimal coefficient combination for special oil spill accident. This paper focuses on proposing a new method to calibrate Lagrangian model with time series of Envisat ASAR images. Oil slicks extracted from time series images form a detected trajectory of special oil slick. Lagrangian model is calibrated by minimizing the difference between simulated trajectory and detected trajectory. mean center position distance difference (MCPD) and rotation difference (RD) of Oil slicks' or particles' standard deviational ellipses (SDEs) are calculated as two evaluations. The two parameters are taken to evaluate the performance of Lagrangian transport model with different coefficient combinations. This method is applied to Penglai 19-3 oil spill accident. The simulation result with calibrated model agrees well with related satellite observations. It is suggested the new method is effective to calibrate Lagrangian model. Copyright © 2016 Elsevier Ltd. All rights reserved.

Chaotic Dynamics in a Low-Energy Transfer Strategy to the Equilateral Equilibrium Points in the Earth-Moon System

NASA Astrophysics Data System (ADS)

Salazar, F. J. T.; Macau, E. E. N.; Winter, O. C.

In the frame of the equilateral equilibrium points exploration, numerous future space missions will require maximization of payload mass, simultaneously achieving reasonable transfer times. To fulfill this request, low-energy non-Keplerian orbits could be used to reach L4 and L5 in the Earth-Moon system instead of high energetic transfers. Previous studies have shown that chaos in physical systems like the restricted three-body Earth-Moon-particle problem can be used to direct a chaotic trajectory to a target that has been previously considered. In this work, we propose to transfer a spacecraft from a circular Earth Orbit in the chaotic region to the equilateral equilibrium points L4 and L5 in the Earth-Moon system, exploiting the chaotic region that connects the Earth with the Moon and changing the trajectory of the spacecraft (relative to the Earth) by using a gravity assist maneuver with the Moon. Choosing a sequence of small perturbations, the time of flight is reduced and the spacecraft is guided to a proper trajectory so that it uses the Moon's gravitational force to finally arrive at a desired target. In this study, the desired target will be an orbit about the Lagrangian equilibrium points L4 or L5. This strategy is not only more efficient with respect to thrust requirement, but also its time transfer is comparable to other known transfer techniques based on time optimization.

Extended Lagrangian Excited State Molecular Dynamics.

PubMed

Bjorgaard, J A; Sheppard, D; Tretiak, S; Niklasson, A M N

2018-02-13

An extended Lagrangian framework for excited state molecular dynamics (XL-ESMD) using time-dependent self-consistent field theory is proposed. The formulation is a generalization of the extended Lagrangian formulations for ground state Born-Oppenheimer molecular dynamics [Phys. Rev. Lett. 2008 100, 123004]. The theory is implemented, demonstrated, and evaluated using a time-dependent semiempirical model, though it should be generally applicable to ab initio theory. The simulations show enhanced energy stability and a significantly reduced computational cost associated with the iterative solutions of both the ground state and the electronically excited states. Relaxed convergence criteria can therefore be used both for the self-consistent ground state optimization and for the iterative subspace diagonalization of the random phase approximation matrix used to calculate the excited state transitions. The XL-ESMD approach is expected to enable numerically efficient excited state molecular dynamics for such methods as time-dependent Hartree-Fock (TD-HF), Configuration Interactions Singles (CIS), and time-dependent density functional theory (TD-DFT).

Differential geometry based solvation model II: Lagrangian formulation.

PubMed

Chen, Zhan; Baker, Nathan A; Wei, G W

2011-12-01

Solvation is an elementary process in nature and is of paramount importance to more sophisticated chemical, biological and biomolecular processes. The understanding of solvation is an essential prerequisite for the quantitative description and analysis of biomolecular systems. This work presents a Lagrangian formulation of our differential geometry based solvation models. The Lagrangian representation of biomolecular surfaces has a few utilities/advantages. First, it provides an essential basis for biomolecular visualization, surface electrostatic potential map and visual perception of biomolecules. Additionally, it is consistent with the conventional setting of implicit solvent theories and thus, many existing theoretical algorithms and computational software packages can be directly employed. Finally, the Lagrangian representation does not need to resort to artificially enlarged van der Waals radii as often required by the Eulerian representation in solvation analysis. The main goal of the present work is to analyze the connection, similarity and difference between the Eulerian and Lagrangian formalisms of the solvation model. Such analysis is important to the understanding of the differential geometry based solvation model. The present model extends the scaled particle theory of nonpolar solvation model with a solvent-solute interaction potential. The nonpolar solvation model is completed with a Poisson-Boltzmann (PB) theory based polar solvation model. The differential geometry theory of surfaces is employed to provide a natural description of solvent-solute interfaces. The optimization of the total free energy functional, which encompasses the polar and nonpolar contributions, leads to coupled potential driven geometric flow and PB equations. Due to the development of singularities and nonsmooth manifolds in the Lagrangian representation, the resulting potential-driven geometric flow equation is embedded into the Eulerian representation for the purpose of

Differential geometry based solvation model II: Lagrangian formulation

PubMed Central

Chen, Zhan; Baker, Nathan A.; Wei, G. W.

2010-01-01

Solvation is an elementary process in nature and is of paramount importance to more sophisticated chemical, biological and biomolecular processes. The understanding of solvation is an essential prerequisite for the quantitative description and analysis of biomolecular systems. This work presents a Lagrangian formulation of our differential geometry based solvation model. The Lagrangian representation of biomolecular surfaces has a few utilities/advantages. First, it provides an essential basis for biomolecular visualization, surface electrostatic potential map and visual perception of biomolecules. Additionally, it is consistent with the conventional setting of implicit solvent theories and thus, many existing theoretical algorithms and computational software packages can be directly employed. Finally, the Lagrangian representation does not need to resort to artificially enlarged van der Waals radii as often required by the Eulerian representation in solvation analysis. The main goal of the present work is to analyze the connection, similarity and difference between the Eulerian and Lagrangian formalisms of the solvation model. Such analysis is important to the understanding of the differential geometry based solvation model. The present model extends the scaled particle theory (SPT) of nonpolar solvation model with a solvent-solute interaction potential. The nonpolar solvation model is completed with a Poisson-Boltzmann (PB) theory based polar solvation model. The differential geometry theory of surfaces is employed to provide a natural description of solvent-solute interfaces. The minimization of the total free energy functional, which encompasses the polar and nonpolar contributions, leads to coupled potential driven geometric flow and Poisson-Boltzmann equations. Due to the development of singularities and nonsmooth manifolds in the Lagrangian representation, the resulting potential-driven geometric flow equation is embedded into the Eulerian representation for

Lagrangian formulation of the general relativistic Poynting-Robertson effect

NASA Astrophysics Data System (ADS)

De Falco, Vittorio; Battista, Emmanuele; Falanga, Maurizio

2018-04-01

We propose the Lagrangian formulation for describing the motion of a test particle in a general relativistic, stationary, and axially symmetric spacetime. The test particle is also affected by a radiation field, modeled as a coherent flux of photons traveling along the null geodesics of the background spacetime, including the general relativistic Poynting-Robertson effect. The innovative part of this work is to prove the existence of the potential linked to the dissipative action caused by the Poynting-Robertson effect in general relativity through the help of an integrating factor, depending on the energy of the system. Generally, such kinds of inverse problems involving dissipative effects might not admit a Lagrangian formulation; especially, in general relativity, there are no examples of such attempts in the literature so far. We reduce this general relativistic Lagrangian formulation to the classic case in the weak-field limit. This approach facilitates further studies in improving the treatment of the radiation field, and it contains, for example, some implications for a deeper comprehension of the gravitational waves.

Identification and Lagrangian analysis of oceanographic structures favorable for fishery of neon flying squid ( Ommastrephes bartramii) in the South Kuril area

NASA Astrophysics Data System (ADS)

Budyansky, M. V.; Prants, S. V.; Samko, E. V.; Uleysky, M. Yu.

2017-09-01

Based on the AVISO velocity field, we compute daily synoptic Lagrangian maps in the South Kuril area for the fishery seasons of 1998, 1999, and 2001-2005 from available catching data on neon flying squid (NFS). With the help of drift maps for artificial particles, we found that the majority of NFS fishing grounds featuring maximum catches are situated near large-scale Lagrangian intrusions: tongues of water penetrating the surrounding water of other Lagrangian properties. It is shown that the NFS catch locations tend to accumulate at places where waters with different magnitudes of certain Lagrangian indicators converge, mix, and produce filaments, swirls, and tendrils typical of chaotic advection. Potential NFS fishing grounds are mainly located near (1) Lagrangian intrusions of the Subarctic front, (2) intrusions of Okhotsk Sea and Oyashio waters around mesoscale anticyclones east of Hokkaido with subsequent penetration of catch locations inside eddies and (3) intrusions of subtropical waters into the central part of the South Kuril area due to interaction with eddies of different size and polarity. Possible reasons for increased biological production and fishery in the vicinity of Lagrangian intrusions are discussed.

Perturbed Equations of Motion for Formation Flight Near the Sun-Earth L2 Point

NASA Technical Reports Server (NTRS)

Segerman, Alan M.; Zedd, Michael F.

2005-01-01

This Memorandum Report consists of a compilation of three individual reports, of increasing complexity, describing investigations of formation flight of spacecraft in the vicinity of the L2 Sun-Earth 1ibration point. The individual reports form the following parts of this compilation: - Introduction to the relative motion of spacecraft about the Sun-Earth L2 Point - Linear and quadratic modelling and solution of the relative motion - Modelling the Perturbations - Elliptical Earth Orbit, Lunar Gravity, Solar Radiation Pressure, Thrusters. The three parts are self-contained, with somewhat, varying notation and terminology. After fair1y significant literature searches: this new work (of Parts 2 and 3) is deemed to be unique because it describes the primary perturbations to the description of relative motion between nearby spacecraft. The effect of the elliptical motion of the Earth about the Sun was verified to be the dominant perturbation to the circular restricted three body problem. Contributions due to lunar gravity and solar radiation pressure are seen to have much smaller effect.

A Constructive Approach to Regularity of Lagrangian Trajectories for Incompressible Euler Flow in a Bounded Domain

NASA Astrophysics Data System (ADS)

Besse, Nicolas; Frisch, Uriel

2017-04-01

The 3D incompressible Euler equations are an important research topic in the mathematical study of fluid dynamics. Not only is the global regularity for smooth initial data an open issue, but the behaviour may also depend on the presence or absence of boundaries. For a good understanding, it is crucial to carry out, besides mathematical studies, high-accuracy and well-resolved numerical exploration. Such studies can be very demanding in computational resources, but recently it has been shown that very substantial gains can be achieved first, by using Cauchy's Lagrangian formulation of the Euler equations and second, by taking advantage of analyticity results of the Lagrangian trajectories for flows whose initial vorticity is Hölder-continuous. The latter has been known for about 20 years (Serfati in J Math Pures Appl 74:95-104, 1995), but the combination of the two, which makes use of recursion relations among time-Taylor coefficients to obtain constructively the time-Taylor series of the Lagrangian map, has been achieved only recently (Frisch and Zheligovsky in Commun Math Phys 326:499-505, 2014; Podvigina et al. in J Comput Phys 306:320-342, 2016 and references therein). Here we extend this methodology to incompressible Euler flow in an impermeable bounded domain whose boundary may be either analytic or have a regularity between indefinite differentiability and analyticity. Non-constructive regularity results for these cases have already been obtained by Glass et al. (Ann Sci Éc Norm Sup 45:1-51, 2012). Using the invariance of the boundary under the Lagrangian flow, we establish novel recursion relations that include contributions from the boundary. This leads to a constructive proof of time-analyticity of the Lagrangian trajectories with analytic boundaries, which can then be used subsequently for the design of a very high-order Cauchy-Lagrangian method.

Variational Lagrangian data assimilation in open channel networks

NASA Astrophysics Data System (ADS)

Wu, Qingfang; Tinka, Andrew; Weekly, Kevin; Beard, Jonathan; Bayen, Alexandre M.

2015-04-01

This article presents a data assimilation method in a tidal system, where data from both Lagrangian drifters and Eulerian flow sensors were fused to estimate water velocity. The system is modeled by first-order, hyperbolic partial differential equations subject to periodic forcing. The estimation problem can then be formulated as the minimization of the difference between the observed variables and model outputs, and eventually provide the velocity and water stage of the hydrodynamic system. The governing equations are linearized and discretized using an implicit discretization scheme, resulting in linear equality constraints in the optimization program. Thus, the flow estimation can be formed as an optimization problem and efficiently solved. The effectiveness of the proposed method was substantiated by a large-scale field experiment in the Sacramento-San Joaquin River Delta in California. A fleet of 100 sensors developed at the University of California, Berkeley, were deployed in Walnut Grove, CA, to collect a set of Lagrangian data, a time series of positions as the sensors moved through the water. Measurements were also taken from Eulerian sensors in the region, provided by the United States Geological Survey. It is shown that the proposed method can effectively integrate Lagrangian and Eulerian measurement data, resulting in a suited estimation of the flow variables within the hydraulic system.

Lagrangian predictability characteristics of an Ocean Model

NASA Astrophysics Data System (ADS)

Lacorata, Guglielmo; Palatella, Luigi; Santoleri, Rosalia

2014-11-01

The Mediterranean Forecasting System (MFS) Ocean Model, provided by INGV, has been chosen as case study to analyze Lagrangian trajectory predictability by means of a dynamical systems approach. To this regard, numerical trajectories are tested against a large amount of Mediterranean drifter data, used as sample of the actual tracer dynamics across the sea. The separation rate of a trajectory pair is measured by computing the Finite-Scale Lyapunov Exponent (FSLE) of first and second kind. An additional kinematic Lagrangian model (KLM), suitably treated to avoid "sweeping"-related problems, has been nested into the MFS in order to recover, in a statistical sense, the velocity field contributions to pair particle dispersion, at mesoscale level, smoothed out by finite resolution effects. Some of the results emerging from this work are: (a) drifter pair dispersion displays Richardson's turbulent diffusion inside the [10-100] km range, while numerical simulations of MFS alone (i.e., without subgrid model) indicate exponential separation; (b) adding the subgrid model, model pair dispersion gets very close to observed data, indicating that KLM is effective in filling the energy "mesoscale gap" present in MFS velocity fields; (c) there exists a threshold size beyond which pair dispersion becomes weakly sensitive to the difference between model and "real" dynamics; (d) the whole methodology here presented can be used to quantify model errors and validate numerical current fields, as far as forecasts of Lagrangian dispersion are concerned.

Eulerian and Lagrangian approaches to multidimensional condensation and collection

NASA Astrophysics Data System (ADS)

Li, Xiang-Yu; Brandenburg, A.; Haugen, N. E. L.; Svensson, G.

2017-06-01

Turbulence is argued to play a crucial role in cloud droplet growth. The combined problem of turbulence and cloud droplet growth is numerically challenging. Here an Eulerian scheme based on the Smoluchowski equation is compared with two Lagrangian superparticle (or superdroplet) schemes in the presence of condensation and collection. The growth processes are studied either separately or in combination using either two-dimensional turbulence, a steady flow or just gravitational acceleration without gas flow. Good agreement between the different schemes for the time evolution of the size spectra is observed in the presence of gravity or turbulence. The Lagrangian superparticle schemes are found to be superior over the Eulerian one in terms of computational performance. However, it is shown that the use of interpolation schemes such as the cloud-in-cell algorithm is detrimental in connection with superparticle or superdroplet approaches. Furthermore, the use of symmetric over asymmetric collection schemes is shown to reduce the amount of scatter in the results. For the Eulerian scheme, gravitational collection is rather sensitive to the mass bin resolution, but not so in the case with turbulence.