Search NASA⌕ Search

SEARCH · Search NASA

Results for “Lagrangian function”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Approximate solutions for half-dark solitons in spinor non-equilibrium Polariton condensates

In this work I generalize and apply an analytical approximation to analyze 1D states of non-equilibrium spinor polariton Bose–Einstein condensates (BEC). Solutions for the condensate wave functions carrying black solitons and half-dark solitons are presented. The derivation is based on the non-conservative Lagrangian formalism for complex Ginzburg–Landau type equations (cGLE), which provides ordinary differential equations for the parameters of the dark soliton solutions in their dynamic environment. Explicit expressions for the stationary dark soliton solution are stated. Subsequently the method is extended to spin sensitive polariton condensates, which yields ordinary differential equations for the parameters of half-dark solitons. Finally a stationary case with explicit expressions for half-dark solitons is presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Variational integrators for reduced magnetohydrodynamics

Reduced magnetohydrodynamics is a simplified set of magnetohydrodynamics equations with applications to both fusion and astrophysical plasmas, possessing a noncanonical Hamiltonian structure and consequently a number of conserved functionals. We propose a new discretisation strategy for these equations based on a discrete variational principle applied to a formal Lagrangian. The resulting integrator preserves important quantities like the total energy, magnetic helicity and cross helicity exactly (up to machine precision). As the integrator is free of numerical resistivity, spurious reconnection along current sheets is absent in the ideal case. If effects of electron inertia are added, reconnection of magnetic field lines is allowed, although the resulting model still possesses a noncanonical Hamiltonian structure. After reviewing the conservation laws of the model equations, the adopted variational principle with the related conservation laws is described both at the continuous and discrete level. We verify the favourable properties of the variational integrator in particular with respect to the preservation of the invariants of the models under consideration and compare with results from the literature and those of a pseudo-spectral code.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2012

On whole Abelian model dynamics

Physics challenge is to determine the objects dynamics. However, there are two ways for deciphering the part. The first one is to search for the ultimate constituents; the second one is to understand its behaviour in whole terms. Therefore, the parts can be defined either from elementary constituents or as whole functions. Historically, science has been moving through the first aspect, however, quarks confinement and complexity are interrupting this usual approach. These relevant facts are supporting for a systemic vision be introduced. Our effort here is to study on the whole meaning through gauge theory. Consider a systemic dynamics oriented through the U(1) - systemic gauge parameter which function is to collect a fields set {l_brace}A{sub {mu}I}{r_brace}. Derive the corresponding whole gauge invariant Lagrangian, equations of motion, Bianchi identities, Noether relationships, charges and Ward-Takahashi equations. Whole Lorentz force and BRST symmetry are also studied. These expressions bring new interpretations further than the usual abelian model. They are generating a systemic system governed by 2N+ 10 classical equations plus Ward-Takahashi identities. A whole dynamics based on the notions of directive and circumstance is producing a set determinism where the parts dynamics are inserted in the whole evolution. A dynamics based on state, collective and individual equations with a systemic interdependence.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

An Immersed Finite Element Lagrangian-Eulerian Code-Coupling Framework

This report presents an assessment of immersed Eulerian-Lagrangian code-coupling techniques suitable for use in a broad range of mechanics applications. The coupling algorithm is based on an immersed finite element method that considers the Lagrangian and Eulerian overlap regions in the overall variational formulation. In this report the basic formulation details are presented followed by various aspects of the code-coupling algorithm using OpenIFEM as the Lagrangian/coupling framework. A series of representative test cases that illustrate the code-coupling algorithm are discussed. The current work provides an in-depth investigation into the immersed finite element method for the purposes of providing a rigorous coupling technique that is minimally invasive in the respective Eulerian and Lagrangian codes. A number of extensions to the base immersed finite element method have been examined. These extension include nodal and quadrature-based indicator functions, a Lagrangian volume-fraction calculation in regions of overlap, and the use of penalty constraints between the Lagrangian and Eulerian domains. A unique MPI-based coupling strategy that retains the independent MPI structure of each code has been demonstrated.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

Trispectrum from co-dimension 2(n) Galileons

A generalized theory of multi-field Galileons has been recently put forward. This model stems from the ongoing effort to embed generic Galileon theories within brane constructions. Such an approach has proved very useful in connecting interesting and essential features of these theories with geometric properties of the branes embedding. We investigate the cosmological implications of a very restrictive multi-field Galileon theory whose leading interaction is solely quartic in the scalar field π and lends itself nicely to an interesting cosmology. The bispectrum is characterized by a naturally small amplitude (f{sub NL}∼<1) and an equilateral shape-function. The trispectrum of curvature fluctuations has features which are quite distinctive with respect to their P(X,φ) counterpart. We also show that, despite an absent cubic Lagrangian in the full theory, non-Gaussianities in this model cannot produce the combination of a small bispectrum alongside with a large trispectrum. We further expand on this point to draw a lesson on what having a symmetry in the full background independent theory entails at the level of fluctuations and vice-versa.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Born-Infeld and charged black holes with non-linear source in f(T) gravity

We investigate f(T) theory coupled with a nonlinear source of electrodynamics, for a spherically symmetric and static spacetime in 4D. We re-obtain the Born-Infeld and Reissner-Nordstrom-AdS solutions. We generalize the no-go theorem for any content that obeys the relationship T{sup 0}{sub 0}=T{sup 1}{sub 1} for the energy-momentum tensor and a given set of tetrads. Our results show new classes of solutions where the metrics are related through b(r)=−Na(r). We do the introductory analysis showing that solutions are that of asymptotically flat black holes, with a singularity at the origin of the radial coordinate, covered by a single event horizon. We also reconstruct the action for this class of solutions and obtain the functional form f(T)=f{sub 0}(−T){sup (N+3)/[2(N+1)]} and L{sub NED}=L{sub 0}(−F){sup (N+3)/[2(N+1)]}. Using the Lagrangian density of Born-Infeld, we obtain a new class of charged black holes where the action reads f(T)=−16β{sub BI}[1−√1+(T/4β{sub BI})].

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Scale separation for multi-scale modeling of free-surface and two-phase flows with the conservative sharp interface method

In this paper we present a scale separation approach for multi-scale modeling of free-surface and two-phase flows with complex interface evolution. By performing a stimulus-response operation on the level-set function representing the interface, separation of resolvable and non-resolvable interface scales is achieved efficiently. Uniform positive and negative shifts of the level-set function are used to determine non-resolvable interface structures. Non-resolved interface structures are separated from the resolved ones and can be treated by a mixing model or a Lagrangian-particle model in order to preserve mass. Resolved interface structures are treated by the conservative sharp-interface model. Since the proposed scale separation approach does not rely on topological information, unlike in previous work, it can be implemented in a straightforward fashion into a given level set based interface model. A number of two- and three-dimensional numerical tests demonstrate that the proposed method is able to cope with complex interface variations accurately and significantly increases robustness against underresolved interface structures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2011

Prototyping method for Bragg-type atom interferometers

We present a method for rapid modeling of new Bragg ultracold atom-interferometer (AI) designs useful for assessing the performance of such interferometers. The method simulates the overall effect on the condensate wave function in a given AI design using two separate elements. These are (1) modeling the effect of a Bragg pulse on the wave function and (2) approximating the evolution of the wave function during the intervals between the pulses. The actual sequence of these pulses and intervals is then followed to determine the approximate final wave function from which the interference pattern can be calculated. The exact evolution between pulses is assumed to be governed by the Gross-Pitaevskii (GP) equation whose solution is approximated using a Lagrangian variational method to facilitate rapid estimation of performance. The method presented here is an extension of an earlier one that was used to analyze the results of an experiment [J. E. Simsarian et al., Phys. Rev. Lett. 85, 2040 (2000)], where the phase of a Bose-Einstein condensate was measured using a Mach-Zehnder-type Bragg AI. We have developed both 1D and 3D versions of this method and we have determined their validity by comparing their predicted interference patterns with those obtained by numerical integration of the 1D GP equation and with the results of the above experiment. We find excellent agreement between the 1D interference patterns predicted by this method and those found by the GP equation. We show that we can reproduce all of the results of that experiment without recourse to an ad hoc velocity-kick correction needed by the earlier method, including some experimental results that the earlier model did not predict. We also found that this method provides estimates of 1D interference patterns at least four orders-of-magnitude faster than direct numerical solution of the 1D GP equation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2012

The cosmology of the Fab-Four

We have recently proposed a novel self tuning mechanism to alleviate the famous cosmological constant problem, based on the general scalar tensor theory proposed by Horndeski. The self-tuning model ends up consisting of four geometric terms in the action, with each term containing a free potential function of the scalar field; the four together being labeled as the Fab-Four. In this paper we begin the important task of deriving the cosmology associated with the Fab-Four Lagrangian. Performing a phase plane analysis of the system we are able to obtain a number of fixed points for the system, with some remarkable new solutions emerging from the trade-off between the various potentials. As well as obtaining inflationary solutions we also find conventional radiation/matter-like solutions, but in regimes where the energy density is dominated by a cosmological constant, and where we do not have any explicit forms of radiation or matter. Stability conditions for matter solutions are obtained and we show how it is possible for there to exist an extended period of 'matter domination' opening up the possibility that we can generate cosmological structures, and recover a consistent cosmology even in the presence of a large cosmological constant.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

A class of customized proximal point algorithms for linearly constrained convex optimization

In this paper, we propose a class of customized proximal point algorithms for linearly constrained convex optimization problems. The algorithms are implementable, provided that the proximal operator of the objective function is easy to evaluate. We show that, with special setting of the algorithmic scalar, our algorithms contain the customized proximal point algorithm (He et al., Optim Appl 56:559–572, 2013), the linearized augmented Lagrangian method (Yang and Yuan, Math Comput 82:301–329, 2013), the Bregman Operator Splitting algorithm (Zhang et al., SIAM J Imaging Sci 3:253–276, 2010) as special cases. The global convergence and worst-case convergence rate measured by the iteration complexity are established for the proposed algorithms. Numerical results demonstrate that the algorithms work well for a wide range of the scalar.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Dynamic Programming and Error Estimates for Stochastic Control Problems with Maximum Cost

This work is concerned with stochastic optimal control for a running maximum cost. A direct approach based on dynamic programming techniques is studied leading to the characterization of the value function as the unique viscosity solution of a second order Hamilton–Jacobi–Bellman (HJB) equation with an oblique derivative boundary condition. A general numerical scheme is proposed and a convergence result is provided. Error estimates are obtained for the semi-Lagrangian scheme. These results can apply to the case of lookback options in finance. Moreover, optimal control problems with maximum cost arise in the characterization of the reachable sets for a system of controlled stochastic differential equations. Some numerical simulations on examples of reachable analysis are included to illustrate our approach.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Ergodic Lagrangian dynamics in a superhero universe

We present a fictional scenario that, while undeniably whimsical, provides the foundation for a unique exercise in extended problem solving, physics analysis, and quantitative model development. Starting with the foundational premise of the Wild Cards shared-world superhero universe, we demonstrate how a variety of concepts appropriate to the advanced undergraduate level—ergodicity, functional analysis, Lagrangian mechanics, and the ever-important simplifying approximation—can be combined into a rich, coherent mathematical model. The goal of this case study is to develop a useful pedagogical exercise in exploring an open-ended research question that presents, at first glance, no clear path forward. Being both eclectic and lengthy, this exercise offers a unique way for students to apply their core physics and mathematics education. It is perhaps best used within a senior honors seminar or within a brief (e.g., January term) elective class.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Some aspects of the five-dimensional Lovelock black hole spacetime: Strong homotopy retract, perihelion precession and quasistationary levels

Highlights: • Some topological aspects are explored. • The perihelion precession is calculated. • The general Heun functions are applied to obtain analytical solutions. • The quasibound states and the radial wave eigenfunctions are obtained. In this work we explore some mathematical physics aspects of the spherically symmetric Lovelock black hole in high dimensions. Intended for this aim, we thoroughly consider the metric corresponding to the five-dimensional Lovelock black hole spacetime. We construct the strong retractions by the geodesic equations on the background under consideration. As a result, from the topological point of view, we construct the theory of strong homotopy retract, which will allow us, in principle, to better understand some of its suitable applications on astrophysics and cosmology, in particular, in the analysis of the spacetime singularities. We find the solutions of the equation of motion for both radial and angular coordinates, and then we describe the outer (“exterior”) and lower (“interior”) apparent horizons. Indeed, the outer apparent horizon is the last surface from which the light waves could still escape from the black hole. Thus, it is meaningful to analyze some physical phenomena related to quantum particles propagating outside the exterior apparent horizon, in particular, we discuss the quasistationary levels of scalar fields and their radial wave functions, which are given in terms of the general Heun functions. We also calculate the perihelion precession in this background.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2012

Parallel heat transport in integrable and chaotic magnetic fields

The study of transport in magnetized plasmas is a problem of fundamental interest in controlled fusion, space plasmas, and astrophysics research. Three issues make this problem particularly challenging: (i) The extreme anisotropy between the parallel (i.e., along the magnetic field), {chi}{sub ||} , and the perpendicular, {chi}{sub Up-Tack }, conductivities ({chi}{sub ||} /{chi}{sub Up-Tack} may exceed 10{sup 10} in fusion plasmas); (ii) Nonlocal parallel transport in the limit of small collisionality; and (iii) Magnetic field lines chaos which in general complicates (and may preclude) the construction of magnetic field line coordinates. Motivated by these issues, we present a Lagrangian Green's function method to solve the local and non-local parallel transport equation applicable to integrable and chaotic magnetic fields in arbitrary geometry. The method avoids by construction the numerical pollution issues of grid-based algorithms. The potential of the approach is demonstrated with nontrivial applications to integrable (magnetic island), weakly chaotic (Devil's staircase), and fully chaotic magnetic field configurations. For the latter, numerical solutions of the parallel heat transport equation show that the effective radial transport, with local and non-local parallel closures, is non-diffusive, thus casting doubts on the applicability of quasilinear diffusion descriptions. General conditions for the existence of non-diffusive, multivalued flux-gradient relations in the temperature evolution are derived.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Using multi-dimensional Smolyak interpolation to make a sum-of-products potential

We propose a new method for obtaining potential energy surfaces in sum-of-products (SOP) form. If the number of terms is small enough, a SOP potential surface significantly reduces the cost of quantum dynamics calculations by obviating the need to do multidimensional integrals by quadrature. The method is based on a Smolyak interpolation technique and uses polynomial-like or spectral basis functions and 1D Lagrange-type functions. When written in terms of the basis functions from which the Lagrange-type functions are built, the Smolyak interpolant has only a modest number of terms. The ideas are tested for HONO (nitrous acid)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Roche-lobe Overflow in Eccentric Planet–Star Systems

Many giant exoplanets are found near their Roche limit and in mildly eccentric orbits. In this study, we examine the fate of such planets through Roche-lobe overflow as a function of the physical properties of the binary components, including the eccentricity and the asynchronicity of the rotating planet. We use a direct three-body integrator to compute the trajectories of the lost mass in the ballistic limit and investigate the possible outcomes. We find three different outcomes for the mass transferred through the Lagrangian point L {sub 1}: (1) self-accretion by the planet, (2) direct impact on the stellar surface, and (3) disk formation around the star. We explore the parameter space of the three different regimes and find that at low eccentricities, e≲0.2, mass overflow leads to disk formation for most systems, while, for higher eccentricities or retrograde orbits, self-accretion is the only possible outcome. We conclude that the assumption often made in previous work that when a planet overflows its Roche lobe it is quickly disrupted and accreted by the star is not always valid.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2011

New Developments in MadGraph/MadEvent

We here present some recent developments of MadGraph/MadEvent since the latest published version, 4.0. These developments include: Jet matching with Pythia parton showers for both Standard Model and Beyond the Standard Model processes, decay chain functionality, decay width calculation and decay simulation, process generation for the Grid, a package for calculation of quarkonium amplitudes, calculation of Matrix Element weights for experimental events, automatic dipole subtraction for next-to-leading order calculations, and an interface to FeynRules, a package for automatic calculation of Feynman rules and model files from the Lagrangian of any New Physics model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗