Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Lagrangian analysis”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Development of an Extended Reality (XR) Tool for Earth Science Visualization

In this presentation, we will discuss our work in adapting the NASA open source XR software, the Mixed Reality Exploration Toolkit (MRET), to an earth science domain. MRET is a NASA open source XR software for rapidly building extended reality (XR) environments for NASA domain problems, e.g., pulling in CAD models of thermal vac chamber and Roman Space Telescope to do fit checks. Primarily used for hardware integration & test, we have been adapting and extending MRET for science problems. Traditionally, scientists view and analyze the result of calculated or measured observables with static 1-D, 2-D or 3-D plots. It can be difficult to identify, track and understand the evolution of key features due to poor viewing angles and the nature of flat computer screens. Additionally, numerical models, such as the NASA GEOS climate model, are almost exclusively formulated and analyzed on Eulerian grids with points fixed in space and time. However, atmospheric phenomena such as convective clouds, hurricanes and wildfire smoke plumes move with the 3-D flow field, and it is often difficult and unnatural to understand these phenomena in an Eulerian reference frame as opposed to the Lagrangian reference frame in which nature operates. As part of an Earth Science Technology Office (ESTO) proposal, we have been adapting MRET to be a scientific exploration and analysis XR tool with integrated Lagrangian Dynamics (LD) for the Goddard Earth Observing System (GEOS) numerical weather prediction model. We believe this will help scientists identify, track, and understand the evolution of Earth Science phenomena. This presentation will discuss current results in our work in developing this XR tool for Earth Science.

Thomas Grubb↗
NASA NTRS1991

Efficient high-resolution rotor wake calculations using flow field reconstruction

Recent computational studies of helicopter flowfields for acoustic applications has led to the development of exceptionally efficient analysis methods for carrying out Lagrangian simulations of the vortex wake. The capabilities of these methods are demonstrated here through calculations involving a variety of wake/rotor interactions associated with both main rotor and tail rotor flowfields. These methods incorporate a flowfield reconstruction approach efficiently to compute vortex wake interactions within arbitrary regions located near the main rotor. This paper summarizes the development of the reconstruction approach and outlines some of its potential uses, while focusing primarily on its application to rotor noise calculations. The development of the nearfield velocity corrections used in the analysis of close vortex interactions is described, and the overall structure of the current reconstruction approach is discussed.

Quackenbush, Todd R.↗
DOE OSTI2022

Cosmological analysis of three-dimensional BOSS galaxy clustering and Planck CMB lensing cross correlations via Lagrangian perturbation theory

In this work, we present a formalism for jointly fitting pre- and post-reconstruction redshift-space clustering (RSD) and baryon acoustic oscillations (BAO) plus gravitational lensing (of the CMB) that works directly with the observed 2-point statistics. The formalism is based upon (effective) Lagrangian perturbation theory and a Lagrangian bias expansion, which models RSD, BAO and galaxy-lensing cross correlations within a consistent dynamical framework. As an example we present an analysis of clustering measured by the Baryon Oscillation Spectroscopic Survey in combination with CMB lensing measured by Planck. The post-reconstruction BAO strongly constrains the distance-redshift relation, the full-shape redshift-space clustering constrains the matter density and growth rate, and CMB lensing constrains the clustering amplitude. Using only the redshift space data we obtain Ω m = 0.303 ± 0.008, $H_0$ = 69.21 ± 0.78 and $σ_8$ = 0.743 ± 0.043. The addition of lensing information, even when restricted to the Northern Galactic Cap, improves constraints to Ω m = 0.303 ± 0.008, $H_0$ = 69.21 ± 0.77 and $σ_8$ = 0.707 ± 0.035, in tension with CMB and cosmic shear constraints. The combination of Ω m and $H_0$ are consistent with Planck, though their constraints derive mostly from redshift-space clustering. The low $σ_8$ value are driven by cross correlations with CMB lensing in the low redshift bin ($\textit{z}$ ≃ 0.38) and at large angular scales, which show a 20% deficit compared to expectations from galaxy clustering alone. We conduct several systematics tests on the data and find none that could fully explain these tensions.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1989

Nonlinear finite element formulation for the large displacement analysis in multibody system dynamics

A total Lagrangian finite element formulation for the deformable bodies in multibody mechanical systems that undergo finite relative rotations is developed. The deformable bodies are discretized using finite element methods. The shape functions that are used to describe the displacement field are required to include the rigid body modes that describe only large translational displacements. This does not impose any limitations on the technique because most commonly used shape functions satisfy this requirement. The configuration of an element is defined using four sets of coordinate systems: Body, Element, Intermediate element, Global. The body coordinate system serves as a unique standard for the assembly of the elements forming the deformable body. The element coordinate system is rigidly attached to the element and therefore it translates and rotates with the element. The intermediate element coordinate system, whose axes are initially parallel to the element axes, has an origin which is rigidly attached to the origin of the body coordinate system and is used to conveniently describe the configuration of the element in undeformed state with respect to the body coordinate system.

Rismantab-Sany, J.↗
DOE OSTI2019

Lagrangian Evolution of the Northeast Pacific Marine Boundary Layer Structure and Cloud during CSET

Flight data from the Cloud System Evolution over the Trades (CSET) campaign over the Pacific stratocumulus-to-cumulus transition are organized into 18 Lagrangian cases suitable for study and future modeling, made possible by the use of a track-and-resample flight strategy. Analysis of these cases shows that 2-day Lagrangian coherence of long-lived species (CO and O 3 ) is high ($r$ = 0.93 and 0.73, respectively), but that of subcloud aerosol, MBL depth, and cloud properties is limited. Although they span a wide range in meteorological conditions, most sampled air masses show a clear transition when considering 2-day changes in cloudiness (−31% averaged over all cases), MBL depth (+560 m), estimated inversion strength (EIS; −2.2 K), and decoupling, agreeing with previous satellite studies and theory. Changes in precipitation and droplet number were less consistent. The aircraft-based analysis is augmented by geostationary satellite retrievals and reanalysis data along Lagrangian trajectories between aircraft sampling times, documenting the evolution of cloud fraction, cloud droplet number concentration, EIS, and MBL depth. An expanded trajectory set spanning the summer of 2015 is used to show that the CSET-sampled air masses were representative of the season, with respect to EIS and cloud fraction. Two Lagrangian case studies attractive for future modeling are presented with aircraft and satellite data. The first features a clear Sc–Cu transition involving MBL deepening and decoupling with decreasing cloud fraction, and the second undergoes a much slower cloud evolution despite a greater initial depth and decoupling state. Finally, potential causes for the differences in evolution are explored, including free-tropospheric humidity, subsidence, surface fluxes, and microphysics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Chiral algebra, localization, modularity, surface defects, and all that

We study the 2D vertex operator algebra (VOA) construction in 4D N = 2 superconformal field theories on S 3 × S 1 , focusing on both old puzzles and new observations. The VOA lives on a two-torus T 2 ⊂ S 3 × S 1 , it is 1 2 Z -graded, and this torus is equipped with the natural choice of spin structure (1,0) for the Z + 1 2 -graded operators, corresponding to the NS sector vacuum character. By analyzing the possible refinements of the Schur index that preserves the VOA, we find that it admits discrete deformations, which allows access to the remaining spin structures (1,1), (0,1), and (0,0), of which the latter two involve the inclusion of a particular surface defect. For Lagrangian theories, we perform the detailed analysis: we describe the natural supersymmetric background, perform localization, and derive the gauged symplectic boson action on a torus in any spin structure. In the absence of flavor fugacities, the 2D and 4D path integrals precisely match, including the Casimir factors. We further analyze the 2D theory: we identify its integration cycle and the two-point functions and interpret flavor holonomies as screening charges in the VOA. Next, we make some observations about modularity; the T-transformation acts on our four partition functions and lifts to a large diffeomorphism on S 3 × S 1 . More interestingly, we generalize the four partition functions on the torus to an infinite family labeled by both the spin structure and the integration cycle inside the complexified maximal torus of the gauge group. Members of this family transform into one another under the full modular group, and we confirm the recent observation that the S-transform of the Schur index in Lagrangian theories exhibits logarithmic behavior. Finally, we comment on how locally our background reproduces the Ω-background.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Energy and momentum conservation in the Euler–Poincaré formulation of local Vlasov–Maxwell-type systems

The action principle by Low for the classic Vlasov–Maxwell system contains a mix of Eulerian and Lagrangian variables. This renders the Noether analysis of reparametrization symmetries inconvenient, especially since the well-known energy- and momentum-conservation laws for the system are expressed in terms of Eulerian variables only. While an Euler–Poincaré formulation of Vlasov–Maxwell-type systems, effectively starting with Low's action and using constrained variations for the Eulerian description of particle motion, has been known for a while Cendra et al, it is hard to come by a documented derivation of the related energy- and momentum-conservation laws in the spirit of the Euler–Poincaré machinery. To our knowledge only one such derivation exists in the literature so far, dealing with the so-called guiding-center Vlasov–Darwin system Sugama et al. The present exposition discusses a generic class of local Vlasov–Maxwell-type systems, with a conscious choice of adopting the language of differential geometry to exploit the Euler–Poincaré framework to its full extent. After reviewing the transition from a Lagrangian picture to an Eulerian one, we demonstrate how symmetries generated by isometries in space lead to conservation laws for linear- and angular-momentum density and how symmetry by time translation produces a conservation law for energy density. Here, we also discuss what happens if no symmetries exist. Finally, two explicit examples will be given—the classic Vlasov–Maxwell and the drift-kinetic Vlasov–Maxwell—and the results expressed in the language of regular vector calculus for familiarity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Energy-Casimir, dynamically accessible, and Lagrangian stability of extended magnetohydrodynamic equilibria

The formal stability analysis of Eulerian extended magnetohydrodynamics (XMHD) equilibria is considered within the noncanonical Hamiltonian framework by means of the energy-Casimir variational principle and the dynamically accessible stability method. Specifically, we find explicit sufficient stability conditions for axisymmetric XMHD and Hall MHD (HMHD) equilibria with toroidal flow and for equilibria with arbitrary flow under constrained perturbations. The dynamically accessible, second-order variation of the Hamiltonian, which can potentially provide explicit stability criteria for generic equilibria, is also obtained. Moreover, we examine the Lagrangian stability of the general quasineutral two-fluid model written in terms of MHD-like variables, by finding the action and the Hamiltonian functionals of the linearized dynamics, working within a mixed Lagrangian-Eulerian framework. Upon neglecting electron mass, we derive a HMHD energy principle, and in addition, the perturbed induction equation arises from Hamilton’s equations of motion in view of a consistency condition for the relation between the perturbed magnetic potential and the canonical variables.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Multi-parametric analysis for mixed integer linear programming: An application to transmission upgrade and congestion management

Upgrading the capacity of existing transmission lines is essential for meeting the growing energy demands, facilitating the integration of renewable energy, and ensuring the security of the transmission system. This study focuses on the selection of lines whose capacities and by how much should be expanded from the perspective of the Independent System Operators (ISOs) to minimize the total system cost. We employ advanced multi-parametric programming and an enhanced branch-and-bound algorithm to address complex mixed-integer linear programming (MILP) problems, considering multi-period time constraints and physical limitations of generators and transmission lines. To characterize the various decisions in transmission expansion, we model the increased capacity of existing lines as parameters within a specified range. This study first relaxes the binary variables to continuous variables and applies the Lagrange method and Karush-Kuhn-Tucker (KKT) conditions to obtain optimal solutions and identify critical regions associated with active and inactive constraints. Moreover, we extend the traditional branch-and-bound (B&B) method by determining the problem’s upper and lower bounds at each node of the B&B decision tree, helping to manage computational challenges in large-scale MILP problems. Here, we compare the difference between the upper and lower bounds to obtain an approximate optimal solution within the decision-makers’ tolerable error range. In addition, the first derivative of the objective function on the parameters of each line is used to inform the selection of lines for easing congestion and maximizing social welfare. Finally, the capacity upgrades are selected by weighing the reductions in system costs against the expense of upgrading line capacities. The findings are supported by numerical simulations and provide transmission-line planners with decision-making guidance.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Eulerian and Lagrangian second-order statistics of superfluid 4 He grid turbulence

Here, we use particle-tracking velocimetry to study Eulerian and Lagrangian second-order statistics of superfluid 4 He grid turbulence. The Eulerian energy spectra at scales larger than the mean distance between quantum vortex lines behave classically with close to Kolmogorov-1941 scaling and are almost isotropic. The Lagrangian second-order structure functions and frequency power spectra, measured at scales comparable with the intervortex distance, demonstrate a sharp transition from nearly classical behavior to a regime dominated by the motion of quantum vortex lines. Employing the homogeneity of the flow, we verify a set of relations that connect various second-order statistical objects that stress different aspects of turbulent behavior, allowing a multifaceted analysis. We use the two-way bridge relations between Eulerian energy spectra and second-order structure functions to reconstruct the energy spectrum from the known second-order velocity structure function and vice versa. The Lagrangian frequency spectrum reconstructed from the measured Eulerian spectrum using the Eulerian-Lagrangian bridge differs from the measured Lagrangian spectrum in the quasiclassical range, which calls for further investigation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Sierra/SolidMechanics 4.54 Goodyear Specific

Sierra/SolidMechanics (Sierra/SM) is a Lagrangian, three-dimensional finite element analysis code for solids and structures subjected to extensive contact and large deformations, encompassing explicit and implicit dynamic as well as quasistatic loading regimes. This document supplements the primary Sierra/SM 4.54 User's Guide, describing capabilities specific to Goodyear analysis use cases, including additional implicit solver options, material models, finite element formulations, and contact settings.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Sierra/SolidMechanics 4.52 Goodyear Specific

Sierra/SolidMechanics (Sierra/SM) is a Lagrangian, three-dimensional finite element analysis code for solids and structures subjected to extensive contact and large deformations, encompassing explicit and implicit dynamic as well as quasistatic loading regimes. This document supplements the primary Sierra/SM 4.52 User's Guide, describing capabilities specific to Goodyear analysis use cases, including additional implicit solver options, material models, finite element formulations, and contact settings.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Sierra/SolidMechanics 5.0 Goodyear Specific

Sierra/SolidMechanics (Sierra/SM) is a Lagrangian, three-dimensional finite element analysis code for solids and structures subjected to extensive contact and large deformations, encompassing explicit and implicit dynamic as well as quasistatic loading regimes. This document supplements the primary Sierra/SM 5.0 User’s Guide, describing capabilities specific to Goodyear analysis use cases, including additional implicit solver options, material models, finite element formulations, and contact settings.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Sierra/SolidMechanics 4.56 Goodyear User's Guide

Sierra/SolidMechanics (Sierra/SM) is a Lagrangian, three-dimensional finite element analysis code for solids and structures subjected to extensive contact and large deformations, encompassing explicit and implicit dynamic as well as quasistatic loading regimes. This document supplements the primary Sierra/SM 4.56 User’s Guide, describing capabilities specific to Goodyear analysis use cases, including additional implicit solver options, material models, finite element formulations, and contact settings.

42 ENGINEERING↗
NASA NTRS1981

Mixed models and reduced/selective integration displacement models for nonlinear analysis of curved beams

Simple mixed models are developed for use in the geometrically nonlinear analysis of deep arches. A total Lagrangian description of the arch deformation is used, the analytical formulation being based on a form of the nonlinear deep arch theory with the effects of transverse shear deformation included. The fundamental unknowns comprise the six internal forces and generalized displacements of the arch, and the element characteristic arrays are obtained by using Hellinger-Reissner mixed variational principle. The polynomial interpolation functions employed in approximating the forces are one degree lower than those used in approximating the displacements, and the forces are discontinuous at the interelement boundaries. Attention is given to the equivalence between the mixed models developed herein and displacement models based on reduced integration of both the transverse shear and extensional energy terms. The advantages of mixed models over equivalent displacement models are summarized. Numerical results are presented to demonstrate the high accuracy and effectiveness of the mixed models developed and to permit a comparison of their performance with that of other mixed models reported in the literature.

Noor, A. K.↗
NASA NTRS1992

A two-phase restricted equilibrium model for combustion of metalized solid propellants

An Eulerian-Lagrangian two-phase approach was adopted to model the multi-phase reacting internal flow in a solid rocket with a metalized propellant. An Eulerian description was used to analyze the motion of the continuous phase which includes the gas as well as the small (micron-sized) particulates, while a Lagrangian description is used for the analysis of the discrete phase which consists of the larger particulates in the motor chamber. The particulates consist of Al and Al2O3 such that the particulate composition is 100 percent Al at injection from the propellant surface with Al2O3 fraction increasing due to combustion along the particle trajectory. An empirical model is used to compute the combustion rate for agglomerates while the continuous phase chemistry is treated using chemical equilibrium. The computer code was used to simulate the reacting flow in a solid rocket motor with an AP/HTPB/Al propellant. The computed results show the existence of an extended combustion zone in the chamber rather than a thin reaction region. The presence of the extended combustion zone results in the chamber flow field and chemical being far from isothermal (as would be predicted by a surface combustion assumption). The temperature in the chamber increases from about 2600 K at the propellant surface to about 3350 K in the core. Similarly the chemical composition and the density of the propellant gas also show spatially non-uniform distribution in the chamber. The analysis developed under the present effort provides a more sophisticated tool for solid rocket internal flow predictions than is presently available, and can be useful in studying apparent anomalies and improving the simple correlations currently in use. The code can be used in the analysis of combustion efficiency, thermal load in the internal insulation, plume radiation, etc.

Sabnis, J. S.↗
NASA NTRS2018

Improved Rice Residue Burning Emissions Estimates: Accounting for Practice-Specific Emission Factors in Air Pollution Assessments of Vietnam

In Southeast Asia and Vietnam, rice residues are routinely burned after the harvest to prepare fields for the next season. Specific to Vietnam, the two prevalent burning practices include: a). piling the residues after hand harvesting; b). burning the residues without piling, after machine harvesting. In this study, we synthesized field and laboratory studies from the literature on rice residue burning emission factors for Particulate Matter less than 2.5 microns (PM2.5). We found significant differences in the resulting burning-practice specific emission factors, with 16.9 grams per square kilogram (plus or minus 6.9) for pile burning and 8.8 grams per square kilogram (plus or minus 3.5) for non-pile burning. We calculated burning practice specific emissions based on rice area data, region-specific fuel-loading factors, combined emission factors, and estimates of burning from the literature. Our results for year 2015 estimate 180 gigagrams of PM2.5 result from the pile burning method and 130 gigagrams result from non-pile burning method, with the most-likely current emission scenario of 150 gigagrams PM2.5 emissions for Vietnam. For comparison purposes, we calculated emissions using generalized agricultural emission factors employed in global biomass burning studies. These results estimate 80 gigagrams PM2.5, which is only 44 percent of the pile burning-based estimates, suggesting underestimation in previous studies. We compare our emissions to an existing all-combustion sources inventory, results show emissions account for 14-18 percent of Vietnam's total PM2.5 depending on burning practice. Within the highly-urbanized and cloud-covered Hanoi Capital region (HCR), we use rice area from Sentinel-1A to derive spatially-explicit emissions and indirectly estimate residue burning dates. Results from HYSPLIT (Hybrid Single-Particle Lagrangian Integrated Trajectory) back-trajectory analysis stratified by season show autumn has most emission trajectories originating in the North, while spring has most originating in the South, suggesting the latter may have bigger impact on air quality. From these results, we highlight locations where emission mitigation efforts could be focused and suggest measures for pollutant mitigation. Our study demonstrates the need to account for emissions variation due to different burning practices.

Air Pollution↗
NASA NTRS2012

Ice Accretion Modeling using an Eulerian Approach for Droplet Impingement

A three-dimensional Eulerian analysis has been developed for modeling droplet impingement on lifting bodes. The Eulerian model solves the conservation equations of mass and momentum to obtain the droplet flow field properties on the same mesh used in CFD simulations. For complex configurations such as a full rotorcraft, the Eulerian approach is more efficient because the Lagrangian approach would require a significant amount of seeding for accurate estimates of collection efficiency. Simulations are done for various benchmark cases such as NACA0012 airfoil, MS317 airfoil and oscillating SC2110 airfoil to illustrate its use. The present results are compared with results from the Lagrangian approach used in an industry standard analysis called LEWICE.

Kim, Joe Woong↗