Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Lagrangian modeling”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

Feasibility of Surfactant-Free Supported Emulsion Liquid Membrane Extraction

Supported emulsion liquid membrane (SELM) is an effective means to conduct liquid-liquid extraction. SELM extraction is particularly attractive for separation tasks in the microgravity environment where density difference between the solvent and the internal phase of the emulsion is inconsequential and a stable dispersion can be maintained without surfactant. In this research, dispersed two-phase flow in SELM extraction is modeled using the Lagrangian method. The results show that SELM extraction process in the microgravity environment can be simulated on earth by matching the density of the solvent and the stripping phase. Feasibility of surfactant-free SELM (SFSELM) extraction is assessed by studying the coalescence behavior of the internal phase in the absence of the surfactant. Although the contacting area between the solvent and the internal phase in SFSELM extraction is significantly less than the area provided by regular emulsion due to drop coalescence, it is comparable to the area provided by a typical hollow-fiber membrane. Thus, the stripping process is highly unlikely to become the rate-limiting step in SFSELM extraction. SFSELM remains an effective way to achieve simultaneous extraction and stripping and is able to eliminate the equilibrium limitation in the typical solvent extraction processes. The SFSELM design is similar to the supported liquid membrane design in some aspects.

Hu, Shih-Yao B.↗
NASA NTRS2002

Variability of ozone loss during Arctic winter (1991 to 2000) estimated from UARS Microwave Limb Sounder measurement

A comprehensive analysis of version 5 Upper Atmosphere Research Satellite (UARS) Microwave Limb Sounder (MLS) ozone data using a Lagrangian Transport (LT) model provides estimates of chemical ozone depletion for the 1991-1992 through 1997-1998 Arctic winters. These new estimates give a consistent, three-dimensional picture of ozone loss during seven Arctic winters; previous Arctic ozone loss estimates from MLS were based on various earlier data versions and were done only for late winter and only for a subset of the years observed by MLS. We find large interannual variability in the amount, timing, and patterns of ozone depletion and in the degree to which chemical loss is masked by dynamical processes.

Arctic ozone loss↗
NASA NTRS2004

Investigation of Polar Stratospheric Cloud Solid Particle Formation Mechanisms Using ILAS and AVHRR Observations in the Arctic

Satellite observations of denitrification and ice clouds in the Arctic lower stratosphere in February 1997 are used with Lagrangian microphysical box model calculations to evaluate nucleation mechanisms of solid polar stratospheric cloud (PSC) particles. The occurrences of ice clouds are not correlated in time and space with the locations of back trajectories of denitrified air masses, indicating that ice particle surfaces are not always a prerequisite for the formation of solid PSCs that lead to denitrification. In contrast, the model calculations incorporating a pseudoheterogeneous freezing process occurring at the vapor-liquid interface can quantitatively explain most of the observed denitrification when the nucleation activation free energy for nitric acid dihydrate formation is raised by only approx.10% relative to the current published values. Once nucleated, the conversion of nitric acid dihydrate to the stable trihydrate phase brings the computed levels of denitrification closer to the measurements. INDEX TERMS: 0305 Atmospheric Composition and Structure: Aerosols and particles (0345, 4801); 0320 Atmospheric Composition and Sb~lctureC: loud physics and chemistry; 0340 Atmospheric Composition and Structure: Middle atmosphere-composition and chemistry

Irie, H.↗
NASA NTRSDate not supplied

Exploring the Elevated Water Vapor Signal Associated With the Free Tropospheric Biomass Burning Plume Over the Southeast Atlantic Ocean

In southern Africa, widespread agricultural fires produce substantial biomass burning (BB) emissions over the region. The seasonal smoke plumes associated with these emissions are then advected westward over the persistent stratocumulus cloud deck in the southeast Atlantic (SEA) Ocean, resulting in aerosol effects which vary with time and location. Much work has focused on the effects of these aerosol plumes, but previous studies have also described an elevated free tropospheric water vapor signal over the SEA. Water vapor influences climate in its own right, and it is especially important to consider atmospheric water vapor when quantifying aerosol–cloud interactions and aerosol radiative effects. Here we present airborne observations made during the NASA ORACLES (ObseRvations of Aerosols above CLouds and their intEractionS) campaign over the SEA Ocean. In observations collected from multiple independent instruments on the NASA P-3 aircraft (from near-surface to 6–7 km), we observe a strongly linear correlation between pollution indicators (carbon monoxide (CO) and aerosol loading) and atmospheric water vapor content, seen at all altitudes above the boundary layer. The focus of the current study is on the especially strong correlation observed during the ORACLES-2016 deployment (out of Walvis Bay, Namibia), but a similar relationship is also observed in the August 2017 and October 2018 ORACLES deployments. Using reanalyses from the European Centre for Medium-Range Weather Forecasts (ECMWF) and Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2), and specialized WRF-Chem simulations, we trace the plume–vapor relationship to an initial humid, smoky continental source region, where it mixes with clean, dry upper tropospheric air and then is subjected to conditions of strong westward advection, namely the southern African easterly jet (AEJ-S). Our analysis indicates that air masses likely left the continent with the same relationship between water vapor and carbon monoxide as was observed by aircraft. This linear relationship developed over the continent due to daytime convection within a deep continental boundary layer (up to ∼5–6 km) and mixing with higher-altitude air, which resulted in fairly consistent vertical gradients in CO and water vapor, decreasing with altitude and varying in time, but this water vapor does not originate as a product of the BB combustion itself. Due to a combination of conditions and mixing between the smoky, moist continental boundary layer and the dry and fairly clean upper-troposphere air above (∼6 km), the smoky, humid air is transported by strong zonal winds and then advected over the SEA (to the ORACLES flight region) following largely isentropic trajectories. Hybrid Single-Particle Lagrangian Integrated Trajectory model (HYSPLIT) back trajectories support this interpretation. This work thus gives insights into the conditions and processes which cause water vapor to covary with plume strength. Better understanding of this relationship, including how it varies spatially and temporally, is important to accurately quantify direct, semi-direct, and indirect aerosol effects over this region.

Kristina Pistone↗
NASA NTRS1978

Investigation of turbulent diffusion in the extreme lower atmosphere

Turbulent diffusion in the extreme lower layer of the atmosphere (up to 5 m) has been investigated. Turbulent flow was simulated under dry, stable and calm conditions by means of a 3.04 m diameter fan installed at a field site situated on flat grassland. The ambient wind was continuously monitored by means of a cup anemometer placed outside the wake, and the temperature distribution was measured by four thermometers placed on an 18 m tower, also outside the wake. Balloons and red smoke were used to visualize the wake flow and investigate the predominant sizes of turbulent eddies and their streamwise behavior. The mean and turbulent velocities along the turbulence line were measured using an array of hot-wire anomometers. Results provide substantial verification of a recently proposed model (Koper and Sadeh, 1975; Koper et al., 1978) relating the Lagrangian to the Eulerian turbulent velocity autocorrelation. In this model the Lagrangian autocorrelation is given by a domain integral over a set of ordinary Eulerian autocorrelations acquired simultaneously at all points within the flow field in question, which is viewed as a turbulence 'box'.

Koper, C. A., Jr.↗
NASA NTRS1988

Modeling of combustion processes of stick propellants via combined Eulerian-Lagrangian approach

This research is motivated by the improved ballistic performance of large-caliber guns using stick propellant charges. A comprehensive theoretical model for predicting the flame spreading, combustion, and grain deformation phenomena of long, unslotted stick propellants is presented. The formulation is based upon a combined Eulerian-Lagrangian approach to simulate special characteristics of the two phase combustion process in a cartridge loaded with a bundle of sticks. The model considers five separate regions consisting of the internal perforation, the solid phase, the external interstitial gas phase, and two lumped parameter regions at either end of the stick bundle. For the external gas phase region, a set of transient one-dimensional fluid-dynamic equations using the Eulerian approach is obtained; governing equations for the stick propellants are formulated using the Lagrangian approach. The motion of a representative stick is derived by considering the forces acting on the entire propellant stick. The instantaneous temperature and stress fields in the stick propellant are modeled by considering the transient axisymmetric heat conduction equation and dynamic structural analysis.

Kuo, K. K.↗
NASA NTRSDate not supplied

Gas-Particle Interaction Model Development in Plume Surface Interaction Erosion and Cratering

The development of a predictive simulation capability for extra-terrestrial Plume Surface Interaction (PSI) environments is undertaken by the Fluid Dynamics Branch at NASA/MSFC under the Game Changing Development (GCD) Program funded by NASA’s Space Technology Mission Directorate (STMD). Predictive simulation capability for propulsive Martian and Lunar landing PSI requires accurate modeling of the complex high-speed plume impingement and resulting gas/particle-cloud and ejecta stream formation. Physics modeling gap analysis during formulation of the PSI project scope identified two particular phenomena of first order importance in gas-particle interactions that lacked existing modeling formulations. The first arises from the lack of models capturing the unsteady drag forces imparted on particles by the rapidly expanding supersonic compressible flow. The second addresses the fluctuating forces and dispersions imparted on both the particle and the fluid resulting from the interference and wake turbulence generated by close proximity particle clouds, dubbed the Particle Turbulent Kinetic Energy(PTKE). Their first order significance has been identified in experiments, but simulation models currently do not exist for either effect. The development of models and the eventual implementation into the Eulerian Gas-Granular Flow Solver (GGFS) simulation tools was constructed as one of four major tasks of the PSI project. In this process, data on particle kinetics and gas-particle interactions are collected from carefully designed experiments of particles embedded in jets. The effects observed in the experiments are then replicated in high-fidelity particle-resolved CFD simulations to inform the formulation of improvements to particle phase drag models for implementations in the more efficient Eulerian-Lagrangian CFD simulations. The resulting models are ultimately ported to the Eulerian-Eulerian models applied for most efficient simulations in PSI production application tools. This paper describes the activities and accomplishments of the past year in the gas-particle interaction modeling task of the PSI project.

Jeff West↗
NASA NTRS1979

Mathematical model investigation of long-term transport of ocean-dumped sewage sludge related to remote sensing

An existing, three-dimensional, Eulerian-Lagrangian finite-difference model was modified and used to examine the transport processes of dumped sewage sludge in the New York Bight. Both in situ and laboratory data were utilized in an attempt to approximate model inputs such as mean current speed, horizontal diffusion coefficients, particle size distributions, and specific gravities. The results presented are a quantitative description of the fate of a negatively buoyant sewage sludge plume resulting from continuous and instantaneous barge releases. Concentrations of the sludge near the surface were compared qualitatively with those remotely sensed. Laboratory study was performed to investigate the behavior of sewage sludge dumping in various ambient density conditions.

Kuo, C. Y.↗
DOE OSTI2024

DNS of Flame Stabilization Dynamics of a Swirl-Stabilized Spray Burner Using Sustainable Aviation Fuels

Concern for emission reduction has motivated the development of new cost-effective alternative sustainable aviation fuels (SAFs). Drop-in SAFs with blends of the certified and the alternative fuels are beneficial because they do not require engine modifications for use in current aviation engines. Combustion characteristics are of particular concern when comparing a new fuel with Jet-A. Lean blow-off and flame stabilization have been identified to be governed by the fuel propensity to autoignition, i.e. derived cetane number. However, fundamental investigation is required to identify the impact of derived cetane number on the combustion mode of flame stabilization in a realistic combustor. The combustion mode has also a direct impact on turbulent combustion modeling closures. DNS simulations are performed in the low-Mach solver of the Pele Suite called PeleLMeX. Lagrangian multi-phase modeling is used to capture the liquid spray injection of Jet-A (reference fuel) and C1 as a representative of a low cetane number SAF. Adaptive Mesh Refinement (AMR) is used to enable a more efficient simulation of a more realistic domain size and embedded boundary treatment is used to model a bluff-body geometry. Local extinction and edge flame propagation were observed for both fuels. The edge flame propagation mode was quantified in terms of a Damkohler number defined as the ratio between progress variable reaction rate and its diffusive flux. Initial analysis suggests that a mixed mode combustion occurs for the edge flame propagation, with flame propagation assisted by ignition and autoignition co-existing for both Jet-A and C1 flames. The analysis of the extinction region shows a larger progress variable for Jet-A which can further increase the local displacement speed and present a faster reignition of the stoichiometric mixture fraction.

adaptive mesh refinement↗
DOE OSTI2023

Lagrangian and Eulerian Supersaturation Statistics in Turbulent Cloudy Rayleigh–Bénard Convection: Applications for LES Subgrid Modeling

Abstract Turbulent fluctuations of scalar and velocity fields are critical for cloud microphysical processes, e.g., droplet activation and size distribution evolution, and can therefore influence cloud radiative forcing and precipitation formation. Lagrangian and Eulerian water vapor, temperature, and supersaturation statistics are investigated in direct numerical simulations (DNS) of turbulent Rayleigh–Bénard convection in the Pi Convection Cloud Chamber to provide a foundation for parameterizing subgrid-scale fluctuations in atmospheric models. A subgrid model for water vapor and temperature variances and covariance and supersaturation variance is proposed, valid for both clear and cloudy conditions. Evaluation of phase change contributions through an a priori test using DNS data shows good performance of the model. Supersaturation is a nonlinear function of temperature and water vapor, and relative external fluxes of water vapor and heat (e.g., during entrainment-mixing and phase change) influence turbulent supersaturation fluctuations. Although supersaturation has autocorrelation and structure functions similar to the independent scalars (temperature and water vapor), the autocorrelation time scale of supersaturation differs. Relative scalar fluxes in DNS without cloud make supersaturation PDFs less skewed than the adiabatic case, where they are highly negatively skewed. However, droplet condensation changes the PDF shape response: it becomes positively skewed for the adiabatic case and negatively skewed when the sidewall relative fluxes are large. Condensation also increases correlations between water vapor and temperature in the presence of relative scalar fluxes but decreases correlations for the adiabatic case. These changes in correlation suppress supersaturation variability for the nonadiabatic cases and increase it for the adiabatic case. Implications of this work for subgrid microphysics modeling using a Lagrangian stochastic scheme are also discussed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Seasonal variation of the surface cross-shelf exchange in the northern South China Sea: a Lagrangian perspective

Previous studies on cross-shelf exchange, predominantly adopted an Eulerian perspective, struggled to identify water sources and pathways. Using a high-resolution regional ocean modeling system (ROMS) and Lagrangian particle tracking, this study systematically investigates the seasonal variation and dynamics of surface cross-shelf exchange in the northern South China Sea (NSCS) from a Lagrangian perspective. Based on daily released 30-day drifter trajectories we identify the key pathways, source regions for surface cross-shelf exchange, revealing pronounced seasonal variability. Results show the surface cross-shelf exchange generally following monsoon-driven Ekman transport. However, local dynamics, especially topographic modulation, can outweigh the expected Ekman-driven transport, producing surface exchange patterns opposite to that predicted from the prevailing winds. Topographic effects vary across different scales. In the coastal waters of western Guangdong during winter, despite downwelling-favorable winds, the modulation of alongshore currents by island topography induces an eastward pressure gradient. This gradient, through geostrophic balance, drives offshore flows opposite to wind-driven Ekman onshore transport. Furthermore, the eastern widened shelf exhibits a distinct seasonal variation of cross-shelf exchange, with strong offshore transport (opposite to the direction of Ekman transport) in winter and exceedingly weak exchange in summer. Analysis of the underlying mechanisms reveals that this winter offshore transport is primarily attributed to geostrophic flows driven by surface pressure gradient that is jointly modulated by the Kuroshio intrusion and local widened shelf topography, and enhanced by cumulative submesoscale processes. In summer, a persistent and strong along-isobath jet acts as a dynamic barrier, effectively suppressing the exchange. These findings highlight the important role of topography in regulating surface material transport, and have important implications for predicting the advection and dispersion of plankton or oil spills over the continental shelf influenced by monsoon.

Hao, Dongliang↗
NASA NTRS1997

Coupling of Sph and Finite Element Codes for Multi-Layer Orbital Debris Shield Design

Particle-based hydrodynamics models offer distinct advantages over Eulerian and Lagrangian hydrocodes in particular shock physics applications. Particle models are designed to avoid the mesh distortion and state variable diffusion problems which can hinder the effective use of Lagrangian and Eulerian codes respectively. However conventional particle-in-cell and smooth particle hydrodynamics methods employ particles which are actually moving interpolation points. A new particle-based modeling methodology, termed Hamiltonian particle hydrodynamics, was developed by Fahrenthold and Koo (1997) to provide an alternative, fully Lagrangian, energy-based approach to shock physics simulations. This alternative formulation avoids the tensile and boundary instabilities associated with standard smooth particle hydrodynamics formulations and the diffusive grid- to-particle mapping schemes characteristic of particle-in-cell methods. In the work described herein, the method of Fahrenthold and Koo has been extended, by coupling the aforementioned hydrodynamic particle model to a hexahedral finite element based description of the continuum dynamics. The resulting continuum model retains all of the features (including general contact-impact effects) of Hamiltonian particle hydrodynamics, while in addition accounting for tensile strength, plasticity, and damage effects important in the simulation of hypervelocity impact on orbital debris shielding. A three dimensional, vectorized, and autotasked implementation of the extended particle method described here has been coded for application to orbital debris shielding design. Source code for the pre-processor (PREP), analysis code (EXOS), post-processor (POST), and rezoner (ZONE), have been delivered separately, along with a User's Guide describing installation and application of the software.

Fahrenthold, Eric P.↗
DOE OSTI2021

Scalable In Situ Computation of Lagrangian Representations via Local Flow Maps

In situ computation of Lagrangian flow maps to enable post hoc time-varying vector field analysis has recently become an active area of research. However, the current literature is largely limited to theoretical settings and lacks a solution to address scalability of the technique in distributed memory. To improve scalability, we propose and evaluate the benefits and limitations of a simple, yet novel, performance optimization. Our proposed optimization is a communication-free model resulting in local Lagrangian flow maps, requiring no message passing or synchronization between processes, intrinsically improving scalability, and thereby reducing overall execution time and alleviating the encumbrance placed on simulation codes from communication overheads. To evaluate our approach, we computed Lagrangian flow maps for four time-varying simulation vector fields and investigated how execution time and reconstruction accuracy are impacted by the number of GPUs per compute node, the total number of compute nodes, particles per rank, and storage intervals. Our study consisted of experiments computing Lagrangian flow maps with up to 67M particle trajectories over 500 cycles and used as many as 2048 GPUs across 512 compute nodes. In all, our study contributes an evaluation of a communication-free model as well as a scalability study of computing distributed Lagrangian flow maps at scale using in situ infrastructure on a modern supercomputer.

Sane, Sudhanshu↗
DOE OSTI2025

Emergence of the polydeterminant in QCD

A generalization of the determinant appears in particle physics in effective Lagrangian interaction terms that model the chiral anomaly in quantum chromodynamics (Giacosa et al. in Phys Rev D 97(9):091901, 2018, Phys Rev D 109(7):L071502, 2024), in particular in connection to mesons. This polydeterminant function, known in the mathematical literature as a mixed discriminant, associates N distinct N x N complex matrices into a complex number and reduces to the usual determinant when all matrices are taken as equal. Here, we explore the main properties of the polydeterminant applied to (quantum) fields by using a formalism and a language close to high-energy physics approaches. We discuss its use as a tool to write down novel chiral anomalous Lagrangian terms and present an explicit illustrative model for mesons. Finally, the extension of the polydeterminant as a function of tensors is shown.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Connecting the Radiative Influences of Aerosol upon the Mass Flux Profiles of Shallow Cumuli across the Southeast Atlantic Ocean Basin and its Boundaries (Final Report)

The Atlantic Ocean covers approximately 25% of Earth’s surface and the atmosphere above it is home to a complex array of clouds and aerosols that have important influences on regional and global weather and climate. These influences must be accurately depicted in short range, medium range, seasonal, and climate forecast models. Conditions over the tropical Atlantic are particularly complex due to continental scale plumes of dust from the Sahara Desert and smoke from agricultural burning in Africa that drift across the Atlantic Ocean basin toward the Americas. These plumes are often found meandering in the lower atmosphere above shallow tropical clouds that form above the ocean surface, presumably mingling with these clouds on occasion due to convective mixing processes. Elevated dust and smoke particles absorb incoming sunlight and substantially warm the marine atmosphere in the layer in which they are present. This warming may alter the thermal stability of the marine atmosphere and may throttle or enhance the development of clouds, change their internal structure and the rate at which they precipitate. Alternatively, it may isolate the lower atmosphere from drier layers above enabling water vapor to accumulate near the ocean surface potentially leading to the development of deeper convection. Our study was organized around the principal concept of determining the impact of African smoke and dust plumes upon cloud development at ASI and understanding how a popular shallow convection parameterization used in models responds to the presence of this aerosol. We analyzed observations collected during the US Department of Energy (DOE) Layered Atlantic Smoke Interactions with Clouds (LASIC) campaign using the Atmospheric Radiation Measurement Program’s Mobile Facility #1 (AMF-1), which was deployed to Ascension Island (ASI) in Southeastern Atlantic for a one-year period. Ascension Island is immediately downwind from an African source of these plumes, but far enough removed to enable the lower atmosphere to have reacted to their presence. A main goal of our study was to compute radiative heating rate profiles over Ascension Island and along the trajectory from the biomass burning regions along coastal Africa to Ascension Island. To compute the radiative heating rate due to aerosols and clouds, we employed observed profiles of temperature, humidity, and clouds from LASIC alongside aerosol optical properties from the Modern-era Retrospective analysis for Research and Applications Version 2 (MERRA-2), as input for the Rapid Radiation Transfer Model (RRTM). Radiative heating was also assessed across the southeast Atlantic Ocean using an ensemble of back trajectories from the Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) model. We were successful in this effort and the resulting publication is already being cited despite its relatively short lifetime in the literature. The second part of our study involved the process-level representation of the clouds observed over the Southeastern Atlantic in models. Our initial task, upon which the balance of this portion of the study depended, was to evaluate the representativeness of the clouds observed by AMF-1. What orographic influence did Ascension Island have on the measured cloud structure? We set out to answer this question by trying to separate observations that were clearly indicative of orographic forcing from those that were representative of open ocean. The siting of AMF-1, which was 340-m above sea-level on the slope of a steep escarpment, proved demonstrably problematic for cloud process measurements. Despite considerable effort using multiple approaches including artificial intelligence (AI), we were unable to successfully compensate for the orographic forcing present at the AMF-1 site on ASI and, hence, produce process-level summaries of the convective mass flux representative of open ocean over the Southeastern Atlantic. Even though the project has officially ended, we are making a final attempt using two new types of clustering algorithm (i.e., AI) on the recommendation of a former student, who is an expert in this area. Since we have only recently begun to test these new AI schemes, the recommendations from the second part of our project outlined below are based on our experience at the time of this report.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

What Controls Crystal Diversity and Microphysical Variability in Cirrus Clouds?

Abstract Variability of ice microphysical properties like crystal size and density in cirrus clouds is important for climate through its impact on radiative forcing, but challenging to represent in models. For the first time, recent laboratory experiments of particle growth (tied to crystal morphology via deposition density) are combined with a state‐of‐the‐art Lagrangian particle‐based microphysics model in large‐eddy simulations to examine sources of microphysical variability in cirrus. Simulated particle size distributions compare well against balloon‐borne observations. Overall, microphysical variability is dominated by variability in the particles' thermodynamic histories. However, diversity in crystal morphology notably increases spatial variability of mean particle size and density, especially at mid‐levels in the cloud. Little correlation between instantaneous crystal properties and supersaturation occurs even though the modeled particle morphology is directly tied to supersaturation based on laboratory measurements. Thus, the individual thermodynamic paths of each particle, not the instantaneous conditions, control the evolution of particle properties.

54 ENVIRONMENTAL SCIENCES↗