Search NASA⌕ Search

SEARCH · Search NASA

Results for “Lagrangian integration”

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 127 records · Page 7

Large-Deformation Displacement Transfer Functions for Shape Predictions of Highly Flexible Slender Aerospace Structures

Large deformation displacement transfer functions were formulated for deformed shape predictions of highly flexible slender structures like aircraft wings. In the formulation, the embedded beam (depth wise cross section of structure along the surface strain sensing line) was first evenly discretized into multiple small domains, with surface strain sensing stations located at the domain junctures. Thus, the surface strain (bending strains) variation within each domain could be expressed with linear of nonlinear function. Such piecewise approach enabled piecewise integrations of the embedded beam curvature equations [classical (Eulerian), physical (Lagrangian), and shifted curvature equations] to yield closed form slope and deflection equations in recursive forms.

Deflection equations↗

Self-Consistent Relativistic Electron Scattering using the Sherlock Scattering Model for X-ray Diagnostics

We present on a new, self-consistent, arbitrary-temperature Romberg integration scheme for modeling electron scattering in materials in a LANL Lagrangian Shock Hydro (LSH) code. Electron beam-target interactions are fundamental to a wide range of scientific and technological applications. When high-energy electron beams hit their target, they may scatter, deposit energy, or ionize the source. These processes govern the behavior and outcomes in nanotechnology manufacturing, electron microscopy, and modern X-ray diagnostics. Simulating these interactions is essential for interpreting experimental results, predicting material responses, and designing efficient tools and experiments. At Los Alamos, this is done using a LSH code, which is a multi-dimension, multi-material, massively parallel, multi-physics code used to simulate applications from asteroid impacts to electron beam interactions. By effectively and efficiently modeling the way that electrons scatter from the beam we can bolster these simulations and more accurately predict experimental outcomes. The model currently implemented in the LSH of interest is based on work by Papp and does not self-consistently preserve momentum in the slightly relativistic regime; here we adopt a model proposed by Braams and Karney and implement a Romberg integration scheme to compute the diffusion tensor. In this paper we will provide background on the Braams-Karney diffusion tensor as well as the Romberg integration scheme we employed to numerically solve for it. We will show that our integration scheme is accurate in solving for the set of scalar potentials used to re-express the diffusion tensor in differential form, and in solving for the diffusion coefficients in the larger LSH code. By using this diffusion tensor rather than the existing Papp one, and numerically integrating it with a Romberg method, we produce much more accurate, self-consistent results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Monitoring estuarine circulation and ocean waste dispersion using an integrated satellite-aircraft-drogue approach

The author has identified the following significant results. An inexpensive, integrated drogue-aircraft-satellite approach was developed which is based on the Lagrangian technique and employs remotely tracked drogues and dyes together with satellite observation of natural tracers, such as suspended sediment. Results include current circulation studies in Delaware Bay in support of an oil slick movement model; investigations of the dispersion and movement of acid wastes dumped 40 miles off the Delaware coast; and coastal current circulation. In each case, the integrated drogue-aircraft-satellite approach compares favorably with other techniques on the basis of accuracy, cost effectiveness, and performance under severe weather conditions.

Klemas, V.↗

Vortex simulation of three-dimensional, spotlike disturbances in a laminar boundary layer

The growth of a turbulent spot in a laminar boundary layer as the spot evolves from a localized disturbance in the layer, is simulated numerically using a three-dimensional vortex filament description of the vorticity field. The filaments are marked with a sequence of node points which are tracked in a Lagrangian reference frame. Velocity computation is done by Biot-Savart integration. Although some discrepancies with experiment appear to exist in the near wall region, the gross properties of the spot, including the velocities of the leading and trailing edges and the velocity perturbations away from the wall, are in good agreement with experiment.

Leonard, A.↗

Vortex simulation of three-dimensional, spotlike disturbances in a laminar boundary layer

The growth of a turbulent spot in a laminar boundary layer, as the spot evolves from a localized disturbance in the layer, is simulated numerically using a three-dimensional vortex filament description of the vorticity field. The filaments are marked with a sequence of mode points which are tracked in a Lagrangian reference frame. Velocity computation is done by Biot-Savart integration. Although some discrepancies with experiment appear to exist in the near wall region, the gross properties of the spot, including the velocities of the leading and trailing edges and the velocity perturbations away from the wall, are in good agreement with experiment.

Leonard, A.↗

A Vertically Lagrangian Finite-Volume Dynamical Core for Global Models

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.

Lin, Shian-Jiann↗

Effect of Space Radiation on TES Detectors Performance

The Athena mission and its X-IFU instrument (X-ray Integral Field Unit) will be positioned at the Sun-Earth Lagrangian point L1, where it will be subject to solar wind (low flux plasma of 95% protons and 5% alpha particles) with energy below 0.1 MeV, and to galactic cosmic rays and solar flares (energetic protons, alpha particles, and electrons) with energies up to hundreds of MeV for protons and the GeV for heavier ions. Some of these particles will go through the satellite and hit the focal plane assembly and hence the detectors. These detectors will be TES (Transition-edge sensor) microcalorimeters, flown for the first time in such an environment. In order to ensure the performance of this type of detectors throughout the duration of such mission, it is critical to study the impact of the radiation on their behavior. Indeed, although a lot of reference material exist for semiconductor detectors such as CCDs, little is currently known about the impact of radiation on TES detectors. These energetic events could cause local heating or damage to the detectors and affect their performance. In this work, we describe how we designed a test campaign to assess the impact of L1 radiation on TES detectors for Athena/X-IFU-like missions and present the results of the tests. Analyses includes assessing changes in the pulse shapes and energy resolution of the detectors measured at 55 mK after several radiation dose steps performed at 4 K.

Sophie Beaumont↗

Tropical amplitudes for colored Lagrangians

Recently a new formulation for scattering amplitudes in Tr(Φ 3 ) theory has been given based on simple combinatorial ideas in the space of kinematic data. This allows all-loop integrated amplitudes to be expressed as “curve integrals” defined using tropical building blocks — the “headlight functions”. This paper shows how the formulation extends to the amplitudes of more general Lagrangians. We will present a number of different ways of introducing tropical “numerator functions” that allow us to describe general Lagrangian interactions. The simplest family of these “tropical numerators” computes the amplitudes of interesting Lagrangians with infinitely many interactions. We also describe methods for tropically formulating the amplitudes for general Lagrangians. One uses a variant of “Wick contraction” to glue together numerator factors for general interaction vertices. Another uses a natural characterization of polygons on surfaces to give a novel combinatorial description of all possible diagrams associated with arbitrary valence interactions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

CIRCLES: Congestion Impacts Reduction via CAV-in-the-loop Lagrangian Energy Smoothing

The energy efficiency of today’s vehicular mobility relies on the un-integrated combination of i) control via static assets (traffic lights, metering, variable speed limits, etc.); and ii) onboard vehicle automation (adaptive cruise control (ACC), ecodriving, etc.). These two families of control were not co-designed and are not engineered to work in coordination. Recent studies have shown i) limitations of controls, and even ii) negative impacts of ACC. This project focused on the technology development, implementation and prototyping, and validation of Mobile Traffic Control (MTC). MTC can be viewed as an extension of classical traffic control (in which static infrastructure actuates traffic flow). In the MTC paradigm, automated vehicles actuate the entire flow via their behavior, offering enhanced possibilities to optimize the energy footprint of traffic, if designed correctly. We set out to demonstrate for the first time that considerably reduced fuel consumption of all vehicles in traffic can be achieved via distributed control of a small proportion of CAVs. Compared to baseline vehicular technologies, our work offers a significant design departure: control algorithms for the CAVs consider the impact one vehicle can have on overall traffic, improving resulting overall fuel consumption. We focus on using a few vehicles (as traffic controllers via CAV technology) to improve the energy efficiency of traffic flow to further optimize energy efficiency. A live-traffic demonstration in November 2022 featured the deployment of 100 specially-equipped CAVs, constituting an approximate local penetration rate upwards of 2%.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Fast Multigrid Reduction-in-Time for Advection via Modified Semi-Lagrangian Coarse-Grid Operators

Many iterative parallel-in-time algorithms have been shown to be highly efficient for diffusion-dominated partial differential equations (PDEs) but are inefficient or even divergent when applied to advection-dominated PDEs. We consider the application of the multigrid reduction-in-time (MGRIT) algorithm to linear advection PDEs. Here, the key to efficient time integration with this method is using a coarse-grid operator that provides a sufficiently accurate approximation to the so-called ideal coarse-grid operator. For certain classes of semi-Lagrangian discretizations, we present a novel semi-Lagrangian-based coarse-grid operator that leads to fast and scalable multilevel time integration of linear advection PDEs. The coarse-grid operator is composed of a semi-Lagrangian discretization followed by a correction term, with the correction designed so that the leading-order truncation error of the composite operator is approximately equal to that of the ideal coarse-grid operator. Parallel results show substantial speed-ups over sequential time integration for variable-wave-speed advection problems in one and two spatial dimensions, and using high-order discretizations up to order five. The proposed approach establishes the first practical method that provides small and scalable MGRIT iteration counts for advection problems.

97 MATHEMATICS AND COMPUTING↗

Safe Deep Reinforcement Learning for Active Distribution System Model Predictive Control with EVs and DERs

The temporal and spatial mismatch between PV generation and electric vehicle (EV) charging and discharging may cause voltage violations in active distribution networks. Despite the widespread use of deep reinforcement learning (DRL) in power system optimization and control, it lacks guarantees on constraint satisfaction during both training and deployment. This paper proposes a Lagrangian-based safe DRL approach for model predictive control (MPC) of active distribution systems with large-scale integration of PVs, EVs, and energy storage systems (ESSs). A Transformer-LSTM time-series model is proposed to forecast EV charging demand, which is then formulated as a constraint to ensure charging requirements are met. Using this prediction, a Lagrangian-based safe soft actor-critic (SAC) framework is developed for real-time control in a three-phase unbalanced distribution system, enforcing voltage safety constraints while optimizing the cumulative net reward. By integrating the forecasting model with multi-period constraints, the proposed framework jointly coordinates PV systems, EV charging and discharging, and ESS scheduling within the MPC horizon. Numerical experiments on a modified IEEE 123-bus system with real-world data show that, under a high PV penetration scenario, the proposed method increases the net reward by 30.74% and reduces average voltage violations from 0.0011 p.u. to 0.0002 p.u. compared with standard SAC. Compared with the optimal power flow (OPF) approach, it achieves similar voltage security while yielding lower line losses. It also maintains real-time control capability, reducing operation latency to 53.21 ms per 15-minute control interval. The proposed method remains effective under varying PV/EV penetrations and load conditions.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Variational energy principle for compressible, baroclinic flow. 1: First and second variations of total kinetic action

The case of a cold gas in the absence of external force fields is considered. Since the only energy involved is kinetic energy, the total kinetic action (i.e., the space-time integral of the kinetic energy density) should serve as the total free-energy functional in this case, and as such should be a local minimum for all possible fluctuations about stable flow. This conjecture is tested by calculating explicit, manifestly covariant expressions for the first and second variations of the total kinetic action in the context of Lagrangian kinematics. The general question of the correlation between physical stability and the convexity of any action integral that can be interpreted as the total free-energy functional of the flow is discussed and illustrated for the cases of rectillinear and rotating shearing flows.

Schmid, L. A.↗

The fermionic universal one-loop effective action

Recent development of path integral matching techniques based on the covariant derivative expansion has made manifest a universal structure of one-loop effective Lagrangians. The universal terms can be computed once and for all to serve as a reference for one-loop matching calculations and to ease their automation. Here we present the fermionic universal one-loop effective action (UOLEA), resulting from integrating out heavy fermions (Dirac or Majorana) with scalar, pseudo-scalar, vector and axial-vector couplings. We also clarify the relation of the new terms computed here to terms previously computed in the literature and those that remain to complete the UOLEA. Our results can be readily used to efficiently obtain analytical expressions for effective operators arising from heavy fermion loops.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Effective interactions of the open bosonic string via field theory

Abstract We describe a method to extract an effective Lagrangian description for open bosonic strings, at zero transcendentality. The method relies on a particular formulation of its scattering amplitudes derived from color-kinematics duality. More precisely, starting from a (DF) 2 + YM quantum field theory, we integrate out all the massive degrees of freedom to generate an expansion in the inverse string tensionα′. We explicitly compute the Lagrangian terms through$$ \mathcal{O} $$ O (α′ 4 ), and target the sector of operators proportional toF 4 to all orders inα′.

Physics↗

Optimal control/structure integrated design of a flexible space platform with articulated appendages

A complete set of closed form symbolic equations are developed for a 2-D model of a system comprised of a flexible platform with articulated appendages. The system considered has a flexible two link manipulator, with joint compliance, at one end of the flexible platform and an articulated antenna on the other end. The model is formulated using a Lagrangian approach which incorporates homogeneous transformation matrices, similar to those typically seen in the robotics literature, for the kinematic representation. The integrated control/structure design strategy considered emphasizes rapid system response and minimal structural weight. Constraints are included to keep antenna pointing error to a minimum, and at the same time achieve accurate end point positioning of the manipulator with minimal oscillations.

Kelkar, Atul G.↗

An asymptotic-preserving semi-Lagrangian algorithm for the anisotropic heat transport equation with arbitrary magnetic fields

Here, we extend the recently proposed semi-Lagrangian algorithm for the extremely anisotropic heat transport equation [Chacón et al., J. Comput. Phys ., 272 (2014)] to deal with arbitrary magnetic field topologies. The original scheme (which showed remarkable numerical properties) was valid for the so-called tokamak-ordering regime, in which the magnetic field magnitude was not allowed to vary much along field lines. The proposed extension maintains the attractive features of the original scheme (including the analytical Green's function, which is critical for tractability) with minor modifications, while allowing for completely general magnetic fields. The accuracy and generality of the approach are demonstrated by numerical experiment with an analytical manufactured solution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Lagrangian mixed layer modeling of the western equatorial Pacific

Processes that control the upper ocean thermohaline structure in the western equatorial Pacific are examined using a Lagrangian mixed layer model. The one-dimensional bulk mixed layer model of Garwood (1977) is integrated along the trajectories derived from a nonlinear 1 1/2 layer reduced gravity model forced with actual wind fields. The Global Precipitation Climatology Project (GPCP) data are used to estimate surface freshwater fluxes for the mixed layer model. The wind stress data which forced the 1 1/2 layer model are used for the mixed layer model. The model was run for the period 1987-1988. This simple model is able to simulate the isothermal layer below the mixed layer in the western Pacific warm pool and its variation. The subduction mechanism hypothesized by Lukas and Lindstrom (1991) is evident in the model results. During periods of strong South Equatorial Current, the warm and salty mixed layer waters in the central Pacific are subducted below the fresh shallow mixed layer in the western Pacific. However, this subduction mechanism is not evident when upwelling Rossby waves reach the western equatorial Pacific or when a prominent deepening of the mixed layer occurs in the western equatorial Pacific or when a prominent deepening of the mixed layer occurs in the western equatorial Pacific due to episodes of strong wind and light precipitation associated with the El Nino-Southern Oscillation. Comparison of the results between the Lagrangian mixed layer model and a locally forced Eulerian mixed layer model indicated that horizontal advection of salty waters from the central Pacific strongly affects the upper ocean salinity variation in the western Pacific, and that this advection is necessary to maintain the upper ocean thermohaline structure in this region.

Shinoda, Toshiaki↗

Vortex methods for two- and three-dimensional flow simulations

The point vortex and vortex blob methods for two dimensional flows are presented. Several results are discussed concerning the numerical analysis of the latter scheme, e.g., the preservation of globally conserved quantities and the analysis of the spatial discretization error resulting from the convection of fixed blobs of vorticity. An application to the two dimensional mixing layer is briefly described. The contour dynamics method is also discussed. The simulation of three dimensional flows with vortex methods is discussed. A natural way to represent the vorticity is in the form of closed tubes of filaments of vorticity, although other schemes are examined. Applications to aircraft trailing vortices and to a turbulent spot in a laminar boundary layer are presented. Hybrid schemes that use an Eulerian mesh to solve the Poisson equation for the velocity field are discussed. The goal of these schemes is to avoid the high cost of the Biot-Savart integration if many vortex elements are used while enjoying most of the advantages of pure Lagrangian schemes.

Leonard, A.↗