Search NASA⌕ Search

SEARCH · Search NASA

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Pressure Gain, Stability, and Operability of Methane/Syngas Based RDEs Under Steady and Transient Conditions (Final Project Report)

The scope of this work addresses key issues associated with losses associated with the detonation wave and other processes internal to the RDE operation, as well as it develops modeling tools for the evaluation of these losses and exhaust emissions in RDEs. The main challenge in studying RDEs is that RDE performance is highly reliant on the specifics of the design so much so that simple/canonical systems alone cannot provide useful engineering information, but practical RDE designs are sufficiently complex and involve extreme operational environments that detailed access either experimentally (laser diagnostics, for instance) or computationally (direct numerical simulations) are as yet to become practical. To overcome this challenge, we have conducted a combined experimental/simulation/analytical study investigating key phenomena that control the characteristics of operation of RDEs. As a result, the study has developed tools and methods that can be used to evaluate performance and design approaches using reduced-physics models, with the assumptions validated using detailed simulations, and the model prediction tested using experimental observations. The specific objectives of the research were: (1) Develop and demonstrate a low-loss fully axial injection concept, taking advantage of stratification effects to alter the detonation structure and position the wave favorably within the combustor; (2) Obtain stability and operability characteristics of an RDE across operating conditions to aid in the development of operability and performance rules for the operations of other systems; and (3) Develop quantitative metrics for performance gain as well as quantitative description of the loss mechanisms through a combination of diagnostics development, reduced-order modeling, and detailed simulations. The work conducted here has made contribution on design of low-loss inlets that has broad application within the power generation industry for use with pressure gain combustion. The operability and stability of different designs, while focusing on axial air inlet designs, has been analyzed. The effect of nozzle and injection conditions was studied. Models and simulations of exhaust emissions, focusing on NOx emission has been developed and used to investigate how operation of the RDE affect NOx production using Lagrangian analysis of RDE simulations. This work has built on previous programs, with the goal of further understanding operation of RDEs and elevate the readiness of design consideration. In addition, a suite of diagnostic and modeling tools have been developed to obtain quantitative metrics on performance based on measurements, which can be readily transferred to other experimental configurations.

08 HYDROGEN↗
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 OSTI2025

Stability analysis of the Eulerian–Lagrangian finite volume methods for nonlinear hyperbolic equations in one space dimension

In this paper, we construct a novel Eulerian–Lagrangian finite volume (ELFV) method for nonlinear scalar hyperbolic equations in one space dimension. It is well known that the exact solutions to such problems may contain shocks though the initial conditions are smooth, and direct numerical methods may suffer from restricted time step sizes. To relieve the restriction, we propose an ELFV method, where the space-time domain was separated by the partition lines originated from the cell interfaces whose slopes are obtained following the Rakine–Hugoniot junmp condition. Unfortunately, to avoid the intersection of the partition lines, the time step sizes are still limited. To fix this gap, we detect effective troubled cells (ETCs) and carefully design the influence region of each ETC, within which the partitioned space-time regions are merged together to form a new one. Then with the new partition of the space-time domain, we theoretically prove that the proposed first-order scheme with Euler forward time discretization is total-variation-diminishing and maximum-principle-preserving with at least twice larger time step constraints than the classical first order Eulerian method for Burgers’ equation. Numerical experiments verify the optimality of the designed time step sizes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Developing a Lagrangian Frame Transformation on Satellite Data to Study Cloud Microphysical Transitions in Arctic Marine Cold Air Outbreaks

Abstract Arctic marine cold air outbreaks (CAOs) generate distinct and dynamic cloud regimes due to intense air‐sea interactions. To understand the temporal evolution of CAO cloud properties and compare different CAO events, a Lagrangian perspective is particularly useful. We developed a novel technique that enables the conversion of inherently Eulerian satellite data into a Lagrangian framework, combining the broad spatiotemporal coverage of satellite observations with the advantages of Lagrangian tracking. This technique was applied to eight CAO cases associated with a recent field campaign. Our results reveal a striking contrast among the cases in terms of cloud‐top phase transitions, providing new insights into the evolution of CAO cloud properties.

Lagrangian analysis↗
DOE OSTI2024

Analysis of DESI × DES using the Lagrangian effective theory of LSS

In this work we use Lagrangian perturbation theory to analyze the harmonic space galaxy clustering signal of the Bright Galaxy Survey (BGS) and luminous red galaxies (LRGs) targeted by the dark energy spectroscopic instrument (DESI), combined with the galaxy-galaxy lensing signal measured around these galaxies using Dark Energy Survey Year 3 source galaxies. The BGS and LRG galaxies are extremely well characterized by DESI spectroscopy and, as a result, lens galaxy redshift uncertainty and photometric systematics contribute negligibly to the error budget of our “2 × 2-point” analysis. On the modeling side, this work represents the first application of the SPINOSAURUS code, implementing an effective field theory model for galaxy intrinsic alignments, and we additionally introduce a new scheme ( MAIAR ) for marginalizing over the large uncertainties in the redshift evolution of the intrinsic alignment signal. Furthermore, this is the first application of a hybrid effective field theory model for galaxy bias based on the A emulus 𝜈 simulations. Our main result is a measurement of the amplitude of the lensing signal, 𝑆 8 = 𝜎 8 ⁢(Ω 𝑚 /0.3) 0.5 = 0.85⁢0$^{+0.042}_{−0.050}$, consistent with values of this parameter derived from the primary cosmic microwave background. This constraint is artificially improved by a factor of 51% if we assume a more standard, but restrictive parametrization for the redshift evolution and sample dependence of the intrinsic alignment signal, and 63% if we additionally assume the nonlinear alignment model. Furthermore, we show that when fixing the cosmological model to the best-fit values from Planck PR4 there is > 5⁢𝜎 evidence for a deviation of the evolution of the intrinsic alignment signal from the functional form that is usually assumed in cosmic shear and galaxy-galaxy lensing studies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Forecasting for ESCAPE: A Multi-Institution Hybrid Forecasting and Nowcasting Operation for Sea-Breeze Convection Supporting a Ground-Based and Airborne Field Campaign

The Experiment of Sea-Breeze Convection, Aerosols, Precipitation and Environment (ESCAPE) field project deployed two aircraft and ground-based assets in the vicinity of Houston, Texas, between 27 May and 2 July 2022, examining how meteorological conditions, dynamics, and aerosols control the initiation, early growth stage, and evolution of coastal convective clouds. To ensure that airborne- and ground-based assets were deployed appropriately, a forecasting and nowcasting team was formed. Daily forecasts guided real-time decision-making by assessing synoptic weather conditions, environmental aerosol, and a variety of atmospheric modeling data to assign a probability for meeting specific ESCAPE campaign objectives. During the research flights, a small team of forecasters provided “nowcasting” support by analyzing radar, satellite, and new model data in real time. The nowcasting team proved invaluable to the campaign operation, as sometimes changing environmental conditions affected, for example, the timing of convective initiation. In addition to the success of the forecasting and nowcasting teams, the ESCAPE campaign offered a unique “testbed” opportunity where in-person and virtual support both contributed to campaign objectives. The forecasting and nowcasting teams were each composed of new and experienced forecasters alike, where new forecasters were given invaluable experience that would otherwise be difficult to attain. Both teams received training on forecast models, map analysis, Hybrid Single-Particle Lagrangian Integrated Trajectory model (HYSPLIT), and thermodynamic sounding analysis before the beginning of the campaign. In this article, the ESCAPE forecasting and nowcasting teams reflect on these experiences, providing potentially useful advice for future field campaigns requiring forecasting and nowcasting support in a hybrid virtual/in-person framework.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Most general EFTs from spurion analysis Hilbert series and minimal lepton flavor violation

We derive a saturation theorem for general Effective Field Theories (EFTs) constructed using spurion analysis. Let S be a set of spurion fields introduced to organize the breaking of a global symmetry G f , and H S be the subgroup of G f that remains unbroken under a generic vacuum expectation value ⟨S⟩; we show that the EFT Lagrangian constructed from the spurion analysis saturates the EFT Lagrangian without the spurions but restricted to H S invariance, provided that arbitrary powers of the spurion fields are allowed. As examples, we study several implementations of the Minimal Lepton Flavor Violation (MLFV) principle, corresponding to various origins of the neutrino masses. In each scenario, we compute the Hilbert series to obtain the numbers of independent lepton flavor covariants that appear in the corresponding EFT at mass dimension 6. These numbers agree with the number of H S invariants in the EFT without the spurions, demonstrating the saturation theorem.

Effective Field Theories↗
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 OSTI2025

Analysis of Σ* via isospin selective reaction 𝐾 𝐿 ⁢𝑝 →𝜋 + ⁢Σ 0

The isospin-selective reaction 𝐾 𝐿 ⁢𝑝 →𝜋 + ⁢Σ 0 provides a clean probe for investigating 𝐼 =1 Σ* resonances. In this work, we perform an analysis of this reaction using an effective Lagrangian approach for the first time, incorporating the well-established Σ⁡(1189)⁢1/2 + , Σ⁡(1385)⁢3/2 + , Σ⁡(1670)⁢3/2 − , and Σ⁡(1775)⁢5/2 − states, while also exploring contributions from other unestablished states. By fitting the available differential cross-section and recoil polarization data, adhering to partial-wave phase conventions same as the Particle Data Group, we find that besides the established resonances, contributions from Σ⁡(1660)⁢1/2 + , Σ⁡(1580)⁢3/2 − , and a Σ*⁢(1/2 − ) improve the description. Notably, a Σ*⁢(1/2 − ) resonance with mass around 1.54 GeV, consistent with Σ⁡(1620)⁢1/2 − , is found to be essential for describing the data in this channel, a stronger indication than found in previous analyses focusing on 𝜋⁢Λ final states. While providing complementary support for Σ⁡(1660)⁢1/2 + and Σ⁡(1580)⁢3/2 − , our results highlight the importance of the Σ⁡(1620)⁢1/2 − region in 𝐾 𝐿 ⁢𝑝 →𝜋 + ⁢Σ 0 . Future high-precision measurements are needed to solidify these findings and further constrain the Σ* spectrum.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Hyporheic-zone Processes and Stream Oxygen Dynamics: Insights from a Multiscale Reactive Transport Model: Modeling Archive

This archive contains the data and Python scripts required to reproduce the analyses and figures in the study: Gomez-Velez, J. D., Rathore, S. S., Cohen, M. J., & Painter, S. L. (2025). Hyporheic-zone Processes and Stream Oxygen Dynamics: Insights from a Multiscale Reactive Transport Model. Submitted to Water Resources Research. The analysis utilizes the subgrid model Advection Dispersion Equation with Lagrangian Subgrids (ADELS) implemented in the Advanced Terrestrial Simulator (ATS; https://amanzi.github.io/ats/stable/). In this case, the ATS and Amanzi versions are (1) ATS version 1.5.1_f5ba18f8 and (2) Amanzi version 1.6-dev_53444cca4. The repository includes a Jupyter Notebook and the necessary data (Pandas DataFrames stored as pickle files) to generate the figures for the manuscript. Additionally, it contains Python scripts to create ATS input files, run the ATS simulations, and post-process the results. Finally, it provides routines for parameter estimation using the Single-Station Metabolism (SSM) model with the Differential Evolution Adaptive Metropolis (DREAM) Markov Chain Monte Carlo (MCMC) algorithm with ZS enhancements (DREAM-ZS).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

An Empirical Quantile Estimation Approach for Chance-Constrained Nonlinear Optimization Problems

We investigate an empirical quantile estimation approach to solve chance-constrained nonlinear optimization problems. Our approach is based on the reformulation of the chance constraint as an equivalent quantile constraint to provide stronger signals on the gradient. In this approach, the value of the quantile function is estimated empirically from samples drawn from the random parameters, and the gradient of the quantile function is estimated via a finite-difference approximation on top of the quantile-function-value estimation. We establish a convergence theory of this approach within the framework of an augmented Lagrangian method for solving general nonlinear constrained optimization problems. The foundation of the convergence analysis is a concentration property of the empirical quantile process, and the analysis is divided based on whether or not the quantile function is differentiable. In contrast to the sampling-and-smoothing approach used in the literature, the method developed in this paper does not involve any smoothing function and hence the quantile-function gradient approximation is easier to implement and there are less accuracy-control parameters to tune. Furthermore, we demonstrate the effectiveness of this approach and compare it with a smoothing method for the quantile-gradient estimation. Numerical investigation shows that the two approaches are competitive for certain problem instances.

Applied Probability↗
DOE OSTI2026

LeWRON: Agentic Analysis of Electroweak Phase Transitions

The electroweak phase transition (EWPT) is a central topic in particle physics and cosmology, connecting collider phenomenology, baryogenesis, and gravitational-wave observatories. Its analysis requires a technically demanding, convention-sensitive, and model-dependent pipeline, from constructing the finite-temperature effective potential to tracking thermal histories, computing bubble nucleation rates, and predicting gravitational-wave spectra. We present LeWRON (Learning ElectroWeak phase tRansitiON), an agentic framework that orchestrates this pipeline starting from an input Lagrangian. LeWRON combines audited toolbox construction with an Explorer module that uses the generated model-specific code for further analysis, including scans and plots. Intermediate analytic outputs are checked by auditor agents and stored as structured artifacts, enabling reproducible human inspection and downstream use through both a command-line interface and a public Python API. The framework supports a reproduction mode, which infers conventions from the literature and reproduces published results, and a discovery mode, which guides users through structured checkpoints for new models. We demonstrate LeWRON across representative beyond-the-Standard-Model scenarios and release the code on GitHub.

Wang, Isaac R. [Fermilab] (ORCID:000000030789218X)↗
DOE OSTI2024

Optimization of the artificial viscosity in Lagrangian staggered discretization codes. Modeling 1D stand-alone shock - case study

We have developed new measures of errors for numerical shock. The new approach is based on analysis of the structure function, and separation of the errors related to oscillations and shock width, which also include error in the position of the ”center” of the numerical shock. We have demonstrated that those measures correctly characterize the numerical solution. We introduced an objective function in, which both types of errors are weighted, and presented optimal values of the coefficients of the linear and quadratic viscosity for different weights and different Mach numbers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Investigating Characteristic Droplet Size Distributions in Large Eddy Simulations of Stratocumulus Clouds

Cloud processes relevant to radiative and precipitation properties depend on the shape of the cloud droplet size distribution. Recent holographic observations revealed that cloud droplet populations do not have the same size distribution shapes throughout but form regions of characteristic distributions with similar microphysical properties. We investigate the existence and properties of these characteristic distributions within Large‐Eddy Simulations of stratocumulus clouds using Lagrangian and bin microphysics schemes. Distribution types are identified, revealing localized characteristic distributions that vary on the scale of the largest convective cell for simulations with bin microphysics. The results from the Lagrangian microphysics scheme hint at similar behavior. Compared to observations, the simulated clouds are much more uniform. Analysis of the LES results suggests a connection to the local entrainment rate, so the poorly resolved entrainment interface in LES may be a cause of the uniformity. The uniformity of the large‐scale forcing could also be a factor.

cloud droplet size distributions↗
DOE OSTI2025

Going Off Grid: A Comparative Study of the Lagrangian and Eulerian Perspectives of New Particle Formation Events

New particle formation and growth (NPF&G) is the process by which ultrafine particles are formed from gas-phase precursors. NPF&G is the dominant source of global aerosol number with important influences on climate. Most observations of NPF&G events are conducted at stationary sites; however, NPF&G observed from stationary sites is influenced by gradual or rapid changes in the air masses passing over the site, complicating NPF&G analysis. In this work, we use observations and a 3D aerosol model to compare aerosol size distributions at a stationary site (Southern Great Plains [SGP] observatory, Oklahoma, USA) and along Lagrangian trajectories crossing the site. The model simulates the NPF&G events reasonably well at SGP. Using the model to compare the Lagrangian and stationary perspectives, we can explain previously unanalyzable days with some evidence of NPF&G as either non-event or analyzable NPF&G days. We find most of the unanalyzable NPF&G days are due to isolated and inhomogeneous NPF&G occurring upwind of the stationary site, often in the outflow of urban regions. Finally, we compare formation rates of 3 nm particles, growth rates, and the survival probability of 3 nm particles growing to 25 nm between the stationary and Lagrangian perspectives. Because of the much larger number of analyzable days along the Lagrangian trajectories, this perspective potentially provides more robust statistics and better characterization of NPF&G event extremes. Our method for extracting chemical/physical properties along Lagrangian trajectories from 3D models can be applied to a wide range of science questions.

O’Donnell, Samuel E. [Colorado State Univ., Fort C↗
DOE OSTI2024

Development of the ARM Lagrangian Large-Scale Forcing Data (ARMLAGTRAJ) Value-Added Product Based on the lagtraj Framework

The Atmospheric Radiation Measurement (ARM) large-scale forcing data developed based on the constrained variational analysis (VARANAL) value-added product (VAP) (Zhang and Lin 1997, Zhang et al. 2001, Xie et al. 2004, Tang et al. 2019) has been widely used for single-column models (SCMs), cloud-resolving models (CRMs), and large-eddy simulation models (LESs) to understand and improve physical processes in models. Recently, the U.S. Department of Energy (DOE) ARM user facility conducted several major field campaigns using ship-based moving observational platforms. For example, the Marine ARM GPCI Investigation of Clouds (MAGIC) field campaign focused on the role of subtropical marine-boundary layer (MBL) clouds, and the Multidisciplinary Drifting Observatory for the Study of Arctic Climate (MOSAiC) field campaign aimed to improve understanding of the coupled climate systems in the Arctic. Observations from moving platforms are critical to provide a comprehensive characterization of coupled-system processes associated with all stages of the cloud and/or sea-ice life cycle. Traditional ARM large-scale forcing data have been developed at fixed locations. They need to be extended to include these moving platforms to address data needs for ship-based field campaigns or to support LES modeling in a Lagrangian framework. With these considerations in mind, we develop ARM-type Lagrangian large-scale forcing data sets based on the lagtraj framework (Boeing et al. 2020) with notable enhancements in generating forcings that are more suitable for ARM field campaigns. The lagtraj is a novel tool that generates forcings for LES and SCM simulation in both Lagrangian and Eulerian perspective. This technical report focuses on the major changes we performed on the lagtraj algorithm and provides an overview of the ARM Lagrangian Large-Scale Forcing Data (ARMLAGTRAJ) value-added products.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

CoCoMET v1.0: a unified open-source toolkit for atmospheric object tracking and analysis

Advances in performance and analysis capabilities have accelerated the development of object tracking algorithms for atmospheric research. This has resulted in a growing number of studies using Lagrangian tracking techniques to analyze the evolution of atmospheric phenomena and the underlying processes. However, the increasing complexity and variety of tracking algorithms present a steep learning curve for new users and make it difficult for existing users to compare algorithm performance. We introduce CoCoMET (Community Cloud Model Evaluation Toolkit), an open-source toolkit that addresses these issues. CoCoMET simplifies the process of running multiple tracking algorithms simultaneously and analyzing objects in both model and observational datasets by specifying parameters in a single configuration file. It standardizes input data from different sources into a consistent format and unifies the tracking output across algorithms. CoCoMET enhances the functionality of existing tracking methods by calculating additional properties such as cell growth and dissipation rates, perimeter, surface area, convexity, and irregularity. In addition, CoCoMET includes a novel method for identifying mergers and splits in 2D and 3D tracks and supports the integration of Eulerian/stationary datasets external to the tracking data for process studies. Its potential utility is demonstrated through examples of model intercomparison, model evaluation against observations, and comparisons between tracking algorithms. Designed for open-source environments, CoCoMET will continue to expand with future releases, incorporating more input data types and tracking algorithms.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Source Levels of In‐Cloud Air in Shallow Cumulus: Consistency Between Paluch Diagram and Lagrangian Particle Tracking

Abstract The Paluch diagram is a widely used tool for interpreting aircraft measurements of shallow cumulus clouds. A prior study conducted by Heus et al. (2008,https://doi.org/10.1175/2008jas2572.1) concluded that the source levels of in‐cloud air inferred from the Paluch diagram exhibit biases, sometimes of several hundred meters, in comparison to those derived from Lagrangian particle tracking. In this short study we revisit this comparison. The results indicate that the upper source levels of in‐cloud air determined from the Lagrangian Particle Tracking and the Paluch diagram are consistent, and the choice of statistical methods is crucial. The significance of this research lies in confirming the reliability of the Paluch analysis, enabling its confident application to aircraft data.

Meteorology & Atmospheric Sciences↗