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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS2019

Status and Examples for the Version 06 IMERG Multi-Satellite Products

After five years of development following the launch of the Global Precipitation Measurement (GPM) missionCore Observatory, the GPM data products are now being extended across the joint Tropical Rainfall MeasuringMission (TRMM) and GPM eras. Version 06 of the U.S. GPM team's Integrated Multi-satellitE Retrievals forGPM (IMERG) merged precipitation product provides a consistent intercalibration for all precipitation productscomputed from individual satellites with the TRMM and GPM Core Observatory sensors as the TRMM- andGPM-era calibrators, respectively, and incorporates monthly surface gauge data. One major change in the basicIMERG algorithm for V06 is that precipitation motion vectors (used to drive the quasi-Lagrangian interpolation,or "morphing") are computed by tracking vertically integrated vapor (TQV) fields analyzed in MERRA2 andGEOS5. This innovation provides globally complete coverage, expanding IMERG's coverage beyond the 60N-Slatitude band previously provided by IR-based vectors, although precipitation over snowy/icy surfaces is stillmasked out as unreliable. A second innovation is that the Quality Index (QI) data field computed for the half-hourlydatasets has been refined to include estimates of correlation at microwave overpass times.We will summarize the processing status for V06 IMERG, for which the retrospective processing shouldbe actively advancing at meeting time. We will show early examples of performance. For example, the TQVmotion vectors are typically slightly better than the IR-based vectors at all latitudes. The transition across theTRMM/GPM data boundary will be discussed, including the necessity of filling in the TRMM-based calibrationsover the latitude band 35-65 in each hemisphere. The notional schedule for the eventual retirement of thepredecessor TRMM Multi-satellite Precipitation Analysis (TMPA) multi-satellite dataset will be updated as well.

Huffman, George J.↗
DOE OSTI2021

Slow manifolds of classical Pauli particle enable structure-preserving geometric algorithms for guiding center dynamics

Since variational symplectic integrators for the guiding center was proposed, structure-preserving geometric algorithms have become an active research field in plasma physics. We found that the slow manifolds of the classical Pauli particle enable a family of structure-preserving geometric algorithms for guiding center dynamics with long-term stability and accuracy. This discovery overcomes the difficulty associated with the unstable parasitic modes for variational symplectic integrators when applied to the degenerate guiding center Lagrangian. It is a pleasant surprise that Pauli's Hamiltonian for electrons, which predated the Dirac equation and marks the beginning of particle physics, reappears in classical physics as an effective algorithm for solving an important plasma physics problem. This technique is applicable to other degenerate Lagrangians reduced from regular Lagrangians.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1974

Numerical calculations of velocity and pressure distribution around oscillating airfoils

An analytical procedure based on the Navier-Stokes equations was developed for analyzing and representing properties of unsteady viscous flow around oscillating obstacles. A variational formulation of the vorticity transport equation was discretized in finite element form and integrated numerically. At each time step of the numerical integration, the velocity field around the obstacle was determined for the instantaneous vorticity distribution from the finite element solution of Poisson's equation. The time-dependent boundary conditions around the oscillating obstacle were introduced as external constraints, using the Lagrangian Multiplier Technique, at each time step of the numerical integration. The procedure was then applied for determining pressures around obstacles oscillating in unsteady flow. The obtained results for a cylinder and an airfoil were illustrated in the form of streamlines and vorticity and pressure distributions.

Bratanow, T.↗
NASA NTRS1983

Ignition of confined gaseous mixtures by hot surfaces and hot wires

Ignition times and spatial and temporal variations of temperature and concentration in gaseous mixtures confined between two infinite parallel walls or two infinite cylinders have been obtained by numerical integration of the appropriate conservation equations written in Lagrangian coordinates. Ignition times and ignition energies are presented for the case of an isothermal wall in terms of the initial mixture pressure and equivalence ratio for both one and two-step chemical reaction mechanisms. The numerical results indicate that there is a critical mixture pressure for which the ignition time is minimum. The values of this critical pressure are larger (smaller) than 1 atm for the one- (two-) step reaction mechanism. The critical pressure for the ignition time is not equal to the critical pressure for the ignition energy. The ignition time and energy decrease with the equivalence ratio within a certain range and then remain constant.

Ramos, J. I.↗
NASA NTRS1992

Finite element Euler calculations of unsteady transonic cascade flows

A Galerkin finite element procedure incorporating an explicit Runge-Kutta time-stepping scheme has been developed in this work to solve unsteady transonic flow in cascades. The computational domain is discretized by a globally unstructured but locally structured blade-fitted deformable mesh. The Galerkin approximation is applied to the unsteady Euler equations based on a mixed Eulerian-Lagrangian description. The semi-discretized equations are integrated forward in time using a multistage Runge-Kutta scheme. An artificial dissipation operator of the type proposed by Jameson is adapted in the current scheme to capture shocks and suppress nonphysical oscillations. Phase-shifted boundary conditions are used to reduce the computational domain to a single reference passage. Results for both steady and unsteady transonic flows through cascades are presented and compared to existing finite volume solutions.

Hsiao, Chingteng↗
NASA NTRS1994

A Lagrangian dynamic subgrid-scale model turbulence

A new formulation of the dynamic subgrid-scale model is tested in which the error associated with the Germano identity is minimized over flow pathlines rather than over directions of statistical homogeneity. This procedure allows the application of the dynamic model with averaging to flows in complex geometries that do not possess homogeneous directions. The characteristic Lagrangian time scale over which the averaging is performed is chosen such that the model is purely dissipative, guaranteeing numerical stability when coupled with the Smagorinsky model. The formulation is tested successfully in forced and decaying isotropic turbulence and in fully developed and transitional channel flow. In homogeneous flows, the results are similar to those of the volume-averaged dynamic model, while in channel flow, the predictions are superior to those of the plane-averaged dynamic model. The relationship between the averaged terms in the model and vortical structures (worms) that appear in the LES is investigated. Computational overhead is kept small (about 10 percent above the CPU requirements of the volume or plane-averaged dynamic model) by using an approximate scheme to advance the Lagrangian tracking through first-order Euler time integration and linear interpolation in space.

Meneveau, C.↗
NASA NTRS2003

Nonlinear Shell Modeling of Thin Membranes with Emphasis on Structural Wrinkling

Thin solar sail membranes of very large span are being envisioned for near-term space missions. One major design issue that is inherent to these very flexible structures is the formation of wrinkling patterns. Structural wrinkles may deteriorate a solar sail's performance and, in certain cases, structural integrity. In this paper, a geometrically nonlinear, updated Lagrangian shell formulation is employed using the ABAQUS finite element code to simulate the formation of wrinkled deformations in thin-film membranes. The restrictive assumptions of true membranes, i.e. Tension Field theory (TF), are not invoked. Two effective modeling strategies are introduced to facilitate convergent solutions of wrinkled equilibrium states. Several numerical studies are carried out, and the results are compared with recent experimental data. Good agreement is observed between the numerical simulations and experimental data.

Tessler, Alexander↗
DOE OSTI2023

QES-Plume v1.0: a Lagrangian dispersion model

Low-cost simulations providing accurate predictions of transport of airborne material in urban areas, vegetative canopies, and complex terrain are demanding because of the small-scale heterogeneity of the features influencing the mean flow and turbulence fields. Common models used to predict turbulent transport of passive scalars are based on the Lagrangian stochastic dispersion model. The Quick Environmental Simulation (QES) tool is a low-computational-cost framework developed to provide high-resolution wind and concentration fields in a variety of complex atmospheric-boundary-layer environments. Part of the framework, QES-Plume, is a Lagrangian dispersion code that uses a time-implicit integration scheme to solve the generalized Langevin equations which require mean flow and turbulence fields. Here, QES-Plume is driven by QES-Winds, a 3D fast-response model that computes mass-consistent wind fields around buildings, vegetation, and hills using empirical parameterizations, and QES-Turb, a local-mixing-length turbulence model. In this paper, the particle dispersion model is presented and validated against analytical solutions to examine QES-Plume’s performance under idealized conditions. In particular, QES-Plume is evaluated against a classical Gaussian plume model for an elevated continuous point-source release in uniform flow, the Lagrangian scaling of dispersion in isotropic turbulence, and a non-Gaussian plume model for an elevated continuous point-source release in a power-law boundary-layer flow. In these cases, QES-Plume yields a maximum relative error below 6 % when compared with analytical solutions. In addition, the model is tested against wind-tunnel data for a uniform array of cubical buildings. QES-Plume exhibits good agreement with the experiment with 99 % of matched zeros and 59 % of the predicted concentrations falling within a factor of 2 of the experimental concentrations. Furthermore, results also emphasize the importance of using high-quality turbulence models for particle dispersion in complex environments. Finally, QES-Plume demonstrates excellent computational performance.

58 GEOSCIENCES↗
NASA NTRS1990

Behaviour of plate elements based on the first-order shear deformation theory

A new analytical technique to assess the performance of shear deformable elements is presented, which makes it possible to determine a priori whether a given element will lock when used to model thin structures. The role that shear constraints play in determining the behavior of thin elements was established by comparing the results of key numerical tests with the predictions of element behavior made by studying the form of the shear constraints. Conclusions regarding locking behavior and the effects of reduced integration in thin shear deformable elements are presented, including the findings (1) that singularity of the shear stiffness matrix is not sufficient to avoid locking; (2) that the effect of mesh refinement on an element that contains spurious constraints is two-fold; and (3) that reduced integration does not remove spurious constraints but rather relaxes them. The results of the study are in agreement with previous studies of Mindlin plate elements in regarding Lagrangian elements as superior to serendipity elements when either full or reduced integration is employed.

Averill, R. C.↗
DOE OSTI2021

Shock-ramp analysis test problem

Quasi-isentropic (ramp) compression is now a well-established experimental method and so are the analysis techniques to give Lagrangian sound speed, pressure, and density along the sample material's isentrope. A shock followed by ramp compression is a natural extension to investigate, for example, shock melt and refreeze on compression, or isentropes of states off the Hugoniot or principal isentrope. In practice, graded-density impactors produce initial shocks, compression by shaped laser pulses may be unable to produce a smooth pressure increase from zero, and incidental perturbations on the drive pulse may also give rise to shocks, so robust shock-ramp analysis methods will be needed. Appropriate analysis methods are needed for shock-ramp experiments, based on those for quasi-isentropic compression, and these require validation. This paper describes three different analyses of a shock-ramp test problem, including an assessment of their estimated errors. The methods tested were based on hydrodynamic characteristics or integration backward in space. All methods gave the known Lagrangian sound speed to within ~1%, and pressure and volume to within less than 2% and 1%, demonstrating that the analysis methods of isentropic compression experiments can be confidently extended to the analysis of shock and ramp compression.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Integration and validation of some modules for modelling of high-speed chemically reactive flows in two-phase gas-droplet mixtures

Three modules are integrated into the built-in OpenFOAM rhoCentralFoam solver towards accurate and efficient modelling of high-speed chemically reactive flows in two-phase gas-droplet mixtures within the OpenFOAM 10.0 framework. The first module is the mixture-averaged diffusion model. The second module is the built-in OpenFOAM Lagrangian solver coupled with optimised droplet drag coefficient and convective heat transfer coefficient sub-models. The last module is a sparse stiff chemistry solver based on dynamic adaptive hybrid integration (AHI-S). The optimised droplet sub-models are first verified in correct implementation for subsequent simulations in this work. Further, they show good accuracy against experimental and analytical data in the modelling of ammonia droplet acceleration and cooling in the flowing and/or low-temperature air. The accuracy and efficiency gains related to the mixture-averaged diffusion model and the AHI-S chemistry solver are examined by simulating 1-D detonation propagation in ammonia droplet-free/laden ammoniaoxygen mixtures. Numerical results of detonation propagation speed, gaseous temperature, density, and species distributions around the induction zone show good agreement with experimental data and analytical solutions. Compared to the built-in OpenFOAM diffusion model, the mixture-averaged diffusion model provides different numerical predictions of pulsating instabilities in detonation propagation. It shows better accuracy in depicting the detonation structure within the droplet-free section attributed to improved multi-component diffusion modelling. Compared to the built-in OpenFOAM solver EulerImplicit (backward Euler), the AHI-S chemistry solver reduces the computational cost by around 50%. It achieves satisfactory accuracy in calculating detonation propagation speed within the droplet-free section with the optimal efficiency when the safety factor, β, equals 0.5.

42 ENGINEERING↗
DOE OSTI2023

A sharp interface Lagrangian-Eulerian method for flexible-body fluid-structure interaction

This paper introduces a sharp-interface approach to simulating fluid-structure interaction (FSI) involving flexible bodies described by general nonlinear material models and across a broad range of mass density ratios. This new flexible-body immersed Lagrangian-Eulerian (ILE) scheme extends our prior work on integrating partitioned and immersed approaches to rigid-body FSI. Our numerical approach incorporates the geometrical and domain solution flexibility of the immersed boundary (IB) method with an accuracy comparable to body-fitted approaches that sharply resolve flows and stresses up to the fluid-structure interface. Unlike many IB methods, our ILE formulation uses distinct momentum equations for the fluid and solid subregions with a Dirichlet-Neumann coupling strategy that connects fluid and solid subproblems through simple interface conditions. As in earlier work, we use approximate Lagrange multiplier forces to treat the kinematic interface conditions along the fluid-structure interface. This penalty approach simplifies the linear solvers needed by our formulation by introducing two representations of the fluid-structure interface, one that moves with the fluid and another that moves with the structure, that are connected by stiff springs. This approach also enables the use of multi-rate time stepping, which allows us to use different time step sizes for the fluid and structure subproblems. Our fluid solver relies on an immersed interface method (IIM) for discrete surfaces to impose stress jump conditions along complex interfaces while enabling the use of fast structured-grid solvers for the incompressible Navier-Stokes equations. The dynamics of the volumetric structural mesh are determined using a standard finite element approach to large-deformation nonlinear elasticity via a nearly incompressible solid mechanics formulation. This formulation also readily accommodates compressible structures with a constant total volume, and it can handle fully compressible solid structures for cases in which at least part of the solid boundary does not contact the incompressible fluid. Selected grid convergence studies demonstrate second-order convergence in volume conservation and in the pointwise discrepancies between corresponding positions of the two interface representations as well as between first and second-order convergence in the structural displacements. The time stepping scheme is also demonstrated to yield second-order convergence. To assess and validate the robustness and accuracy of the new algorithm, comparisons are made with computational and experimental FSI benchmarks. Test cases include both smooth and sharp geometries in various flow conditions. Furthermore, we also demonstrate the capabilities of this methodology by applying it to model the transport and capture of a geometrically realistic, deformable blood clot in an inferior vena cava filter.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A Lagrangian dual method for two-stage robust optimization with binary uncertainties

This report presents a new exact method to calculate worst-case parameter realizations in two-stage robust optimization problems with categorical or binary-valued uncertain data. Traditional exact algorithms for these problems, notably Benders decomposition and column-and-constraint generation, compute worst-case parameter realizations by solving mixed-integer bilinear optimization subproblems. However, their numerical solution can be computationally expensive not only due to their resulting large size after reformulating the bilinear terms, but also because decision-independent bounds on their variables are typically unknown. We propose an alternative Lagrangian dual method that circumvents these difficulties and is readily integrated in either algorithm. We specialize the method to problems where the binary parameters switch on or off constraints as these are commonly encountered in applications, and discuss extensions to problems that lack relatively complete recourse and to those with integer recourse. Numerical experiments provide evidence of significant computational improvements over existing methods.

42 ENGINEERING↗
NASA NTRS2018

Extratropical Cyclone Precipitation Life Cycles: A Satellite-Based Analysis

Although extratropical cyclones are the most common midlatitude storms, the relationship between their precipitation life cycle and dynamical strength life cycle has not been thoroughly analyzed. Given that thermodynamic heat exchanges associated with precipitation impact cyclone circulation, there is a need to understand the precipitation/dynamics relationship. Based on Integrated MultisatellitE Retrievals for Global Precipitation MeasurementGPM precipitation and Lagrangian cyclone tracks, the precipitation maximum occurs prior to the dynamical strength maximum 70% of the time. The lag in timing is consistent with the difference in cyclone precipitable water vapor at the two peaks. Conditional subsetting of the cyclone composites shows that if the precipitable water vapor distribution is constrained to be equal throughout the composite life cycle, the precipitation peak occurs very near the time of the peak in cyclone dynamical strength. Thus, the boost in dynamical strength caused by latent heat associated with precipitation manifest itself with little to no time lag.

Extratropical cyclone↗
NASA NTRSDate not supplied

IMERG Multi-Satellite Products Across Two Decades

The Version 06 Global Precipitation Measurement (GPM) mission products were completed over the last year, capping five years of development since the launch of the GPM Core Observatory, and covering the joint Tropical Rainfall Measuring Mission (TRMM) and GPM eras with consistently processed algorithms. The U.S. GPM team’s Integrated Multi-satellitE Retrievals for GPM (IMERG) merged precipitation product enforces a consistent intercalibration for all precipitation products computed from individual satellites with the TRMM and GPM Core Observatory sensors as the TRMM- and GPM-era calibrators, respectively, and incorporates monthly surface gauge data in the Final (research) product. Mid-latitude calibrations during the TRMM era necessarily are more approximate because TRMM only covered the latitude band 35°N-S, while GPM covers 65°N-S. Starting in V06, IMERG employs precipitation motion vectors (used to drive the quasi-Lagrangian interpolation, or “morphing”) that are computed by tracking the vertically integrated vapor as analyzed in MERRA2 and GEOS FP. This approach covers the entire globe, expanding coverage beyond the 60°N-S latitude band provided by IR-based vectors in previous versions, although we choose to mask out microwave-based precipitation over snowy/icy surfaces as unreliable. We will provide examples of performance for the V06 IMERG products, including comparison with the long-term record of GPCP and TMPA, showing higher values by about 8% in the latitude band 50°N-S over oceans; diurnal cycle, demonstrating improvement over previous versions; and daily precipitation PDFs for the entire record, showing a shift at the TRMM/GPM boundary, as well as interannual variations. These analyses have important implications for the utility of V06 IMERG data for long-record calculations. Finally, we will review the retirement of the predecessor TMPA multi-satellite dataset.

George J. Huffman↗
DOE OSTI2024

The maintenance of coherent vortex topology by Lagrangian chaos in drift-Rossby wave turbulence

Here, this work introduces the “potential vorticity bucket brigade,” a mechanism for explaining the resilience of vortex structures in magnetically confined fusion plasmas and geophysical flows. Drawing parallels with zonal jet formation, we show how inhomogeneous patterns of mixing can reinforce, rather than destroy non-zonal flow structure. We accomplish this through an exact stochastic Lagrangian representation of vorticity transport, together with a near-integrability property, which relates coherent flow topology to fluid relabeling symmetries. We demonstrate these ideas in the context of gradient-driven magnetized plasma turbulence, though the tools we develop here are model-agnostic and applicable beyond the system studied here.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Nonlinear Lifshitz photon theory in condensed matter systems

We present an interacting theory of a U(1) gauge boson with a quadratic dispersion relation, which we call the "nonlinear Lifshitz photon theory.'' The Lifshitz photon is a three-dimensional generalization of the Tkachenko mode in rotating superfluids. Starting from the Wigner crystal of charged particles coupled to a dynamical U(1) gauge field, after integrating out gapped degrees of freedom, we arrive at the Lagrangian for the nonlinear Lifshitz photon. The symmetries of the theory include a global U(1) 1-form symmetry and nonlinearly realized "magnetic" translation and rotation symmetries. The interaction terms in the theory lead to the decay of the Lifshitz photon, the rate of which we estimate. We show that the Wilson loop, which plays the role of the order parameter of the spontaneous breaking of the 1-form global symmetry, deviates from the perimeter law by an additional logarithmic factor. Here, we explore potential connections to other condensed matter systems, with a particular focus on quantum spin ice and ferromagnets. Finally, we generalize our theory to higher dimensions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
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↗