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 55 records · Page 3

An IR Sounding-Based Analysis of the Saharan Air Layer in North Africa

Intense daytime surface heating over barren-to-sparsely vegetated surfaces results in dry convective mixing. In the absence of external forcing such as mountain waves, the dry convection can produce a deep, well-mixed, nearly isentropic boundary layer that becomes a well-mixed residual layer in the evening. These well-mixed layers (WML) retain their unique mid-tropospheric thermal and humidity structure for several days. To detect the SAL and characterize its properties, AIRS Level 2 Ver. 6 temperature and humidity products (2003-Present) are evaluated against rawinsondes and compared to model analysis at each of the 55 rawinsonde stations in northern Africa. To distinguish WML from Saharan air layers (WMLs of Saharan origin), the detection involved a two-step process: 1) algorithm-based detection of WMLs in dry environments (less than 7 g per kilogram mixing ratio) 2) identification of Sahara air layers (SAL) by applying Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) back trajectories to determine the history of each WML. WML occurrence rates from AIRS closely resemble that from rawinsondes, yet rates from model analysis were up to 30% higher than observations in the Sahara due to model errors. Despite the overly frequent occurrence of WMLs from model analysis, HYSPLIT trajectory analysis showed that SAL occurrence rates (given a WML exists) from rawinsondes, AIRS, and model analysis were nearly identical. Although the number of WMLs varied among the data sources, the proportion of WMLs which were classified as SAL was nearly the same. The analysis of SAL bulk properties showed that AIRS and model analysis exhibited a slight warm and moist bias relative to rawinsondes in non-Saharan locations, but model analysis was notably warmer than rawinsondes and AIRS within the Sahara. The latter result is likely associated with the dearth of available data assimilated by model analysis in the Sahara. The variability of SAL thicknesses was reasonably captured by both AIRS and model analysis, but the former favor layers than are thinner than observations. Finally, further analysis of HYSPLIT trajectories revealed that fewer than 10% and 33% of all SAL back trajectories passed through regions with notable precipitation (>100 mm accumulated along the trajectory path) or Aerosol Optical Depth (AOD greater than 0.4, 75th percentile of AOD) on average, respectively. Trajectory analysis indicated that only 57% of Saharan and 24% of non-Saharan WMLs are definitively of Saharan origin (Saharan requirement: Two consecutive days in Sahara and 24 or more of those hours within 72 hours of detection). Non-SAL WMLs either originate from local-to-regionally generated residual layers or from mid-latitude air streams that do not linger over the Sahara for a sufficient time period. Initial analysis shows these non-SAL WMLs tend to be both notably cooler and slightly moister than their SAL counter parts. Continuing analysis will address what role Saharan and non-Saharan air masses characteristics may play on local and regional environmental conditions.

Nicholls, Stephen D.↗

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↗

Characterizing Differences in the Aerosol Plume and Cloud Structure over Ascension Island During the 2016 and 2017 Biomass Burning Seasons

Marine boundary layer clouds, including the transition from stratocumulus to cumulus, are poorly represented in numerical weather prediction and general circulation models. In many cases, the complex physical relationships between marine boundary cloud morphology and the environmental conditions in which the clouds exist are not well understood. Such uncertainties arise in the presence of biomass burning carbonaceous aerosol, as is the case over the southeast Atlantic Ocean. It is likely that the absorbing and heating properties of these aerosols influence the microphysical composition and macrophysical arrangement of marine stratocumulus and trade cumulus in this region; however, this has yet to be quantified. The deployment of the Atmospheric Radiation Measurement Mobile Facility #1 (AMF1) in support of LASIC (Layered Atlantic Smoke Interactions with Clouds) provided a unique opportunity to collect observations of cloud and aerosol properties during two consecutive biomass burning seasons during July through October of 2016 and 2017 over Ascension Island (7.96 S, 14.35 W). Through the use of AMF1 observations, the Modern Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2), and back trajectories from the Hybrid Single-Particle Lagrangian Integrated Trajectory model (HYSPLIT), it will be demonstrated that differences in the atmospheric circulation during the two years result in varying aerosol conditions over Ascension Island. When the aerosol plume is overhead, the aerosol loading is higher during the 2016 season as a result of a weaker subtropical high-pressure system. Furthermore, the aerosol plume originates from central Africa in 2016, but further south in 2017. Contrasts in the season-to-season and day-to-day aerosol loading are used to categorize boundary layer cloud and sub-cloud turbulence measurements above Ascension Island using the AMF1 Doppler lidar and cloud radar.

Collow, Allison B. Marquardt↗

The Impact of MISR-Derived Injection Height Initialization on Wildfire and Volcanic Plume Dispersion in the HYSPLIT Model

The dispersion of particles from wildfires, volcanic eruptions, dust storms, and other aerosol sources can affect air quality and other environmental factors downwind. Aerosol injection height is one source attribute that mediates downwind dispersion, as wind speed and direction can vary dramatically with elevation. Using plume heights derived from space-based, multi-angle imaging, we examine the impact of initializing plumes with satellite-measured vs. nominal (model-calculated or VAAC observations) injection height on the simulated dispersion of six large aerosol plumes. When there are significant differences in nominal vs. satellite-derived particle injection heights, or if one injection height is within the planetary boundary layer (PBL) and the other is above the PBL, differences in simulation results can arise. In the cases studied with significant nominal vs. satellite-derived injection height differences, the NOAA Air Resources Laboratory's Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model tends to represent plume evolution better if the injection height in the model is constrained by hyper-stereo satellite retrievals.

Plume heights↗

Radiative Heating Rate Profiles over the Southeast Atlantic Ocean during the 2016 and 2017 Biomass Burning Seasons

Marine boundary layer clouds, including the transition from stratocumulus to cumulus, are poorly represented in numerical weather prediction and general circulation models. Further uncertainties in the cloud structure arise in the presence of biomass burning carbonaceous aerosol, as is the case over the southeast Atlantic Ocean where biomass burning aerosol is transported from the African continent. As the aerosol plume progresses across the southeast Atlantic Ocean, radiative heating within the aerosol layer has the potential to alter the thermodynamic environment and therefore the cloud structure; however, limited work has been done to quantify this along the trajectory of the aerosol plume in the region. The deployment of the First Atmospheric Radiation Measurement Mobile Facility in support of the Layered Atlantic Smoke Interactions with Clouds field campaign provided a unique opportunity to collect observations of cloud and aerosol properties during two consecutive biomass burning seasons during July through October of 2016 and 2017 over Ascension Island (7.96 °S, 14.35 °W). Using observed profiles of temperature, humidity, and clouds from the field campaign, alongside aerosol optical properties from MERRA-2 as input for the Rapid Radiation Transfer Model, profiles of the radiative heating rate due to aerosols and clouds were computed. Radiative heating is also assessed across the southeast Atlantic Ocean using an ensemble of back trajectories from the Hybrid Single Particle Lagrangian Integrated Trajectory Model. Idealized experiments using the Rapid Radiation Transfer Model with and without aerosols and a range of values for the single scattering albedo demonstrate that shortwave heating within the aerosol layer above Ascension Island can locally range between 2 and 8 K per day depending on the aerosol optical properties, though impacts of the aerosol can be felt elsewhere in the atmospheric column. When considered under clear conditions, the aerosol has a cooling effect at the TOA and based on the observed cloud properties at Ascension Island, the cloud albedo is not large enough to overcome this. Shortwave radiative heating due to biomass burning aerosol is not balanced by additional longwave cooling, and the net radiative impact results in a stabilization of the lower troposphere. However, these results are extremely sensitive to the single scattering albedo assumptions in models.

Allison B. Marquardt Collow↗

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↗

Assessment of aerosol burden over Ghana

Although air pollution in Ghana is ranked number one in environmental health threats to public health and sixth to cause of deaths, routine monitoring is rare. This paper presents fourteen years (2005-2018) assessment of aerosol optical depth (AOD) at 3 km resolution from MODIS Aqua and Terra satellites to ascertain the Spatio-temporal and seasonal distribution of aerosols over Ghana and its major cities. The MODIS AOD at 3 km were validated against ground-based Aerosol Robotic Network (AERONET) AODs to ascertain the suitability of the MODIS 3 km data for air quality application in the region. The contribution of distant aerosols to city aerosol loadings was also assessed with Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) backscatter model. A moderate-high aerosol burden (AODs ~ 0.50) was observed over Ghana with a significant contribution from the pre-monsoon season. City centres of Takoradi and Kumasi showed higher aerosol loads (AODs ~ 0.80) than Accra and Tamale. The HYSPLIT model showed that distant or transported aerosol sources to the city centres were of both marine and land generated origins. Linear regression analysis between MODIS AOD and AERONET AOD showed a reasonably good correlation of ~ 0.60 for Aqua and Terra. From the validation analysis, both Aqua and Terra satellites can be used for air quality monitoring over Ghana; however, more ground research must be conducted to ascertain better aerosol model assumptions for the region.

Aerosols↗

Enhanced Monitoring and Forecasting of South Asian Air Quality Episodes with Multi-Sensor Satellite Products and Dispersion Modeling

Air pollution and particulates represent a serious environmental, public health, and overarching societal concern in the Hindu Kush Himalaya (HKH) region of south-central Asia, especially during the boreal winter and spring months. Frequent contributors to poor air quality of the region include dust and particulate matter transport from the Middle East region and western India, persistent nocturnal fog and smog during the stable dry monsoon months, and biomass burning during the pre-monsoon months of boreal spring. Nocturnal and multi-day persistent stable fog and smog events can occur during the stable conditions of the dry monsoon months that lead to hazardous visibility reductions (e.g., at major airports such as Delhi, India)and poor air quality, particularly in urban corridors.Our project goal is therefore to design and implement a robust air quality and chemistry observation and modeling system using the Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) and the Weather Research and Forecasting coupled with Chemistry (WRF-Chem) models. This system effectively assimilates aerosol and trace gas retrievals from geostationary and polar-orbiting satellites to advance the monitoring and prediction capabilities of air quality and visibility reductions throughout the HKH region. For this presentation, we highlight several[multi-spectral] geostationary satellite products and dispersion modeling capabilities being developed and implemented for improved situational awareness and forecasting of air quality episodes.Launched in December 2018, the Advanced Meteorological Imager (AMI) aboard the Geostationary-Korean Multi-Purpose Satellite-2A (GEO-KOMPSAT-2A, or GK2A) provides 16 channels of multi-spectral information similar to the U.S. GOES-16/17 satellites. Proven community recipes for generating multi-spectral Red-Green-Blue (RGB) composite products to detect specific meteorological phenomena are applied to the AMI data for the HKH region. We implemented four such RGB products consisting of truecolor with Rayleigh correction to aid in smoke detection; dust RGB for monitoring larger particulate matter; nighttime microphysics for fog/smog detection;and natural color fire RGB and the shortwave 3.8 micron channel for fire hot-spot detection.We additionally configured HYSPLIT for dust/sand dispersion and concentration forecasts over HKH. The HYSPLIT dispersion simulations are initially being conducted with two different methodologies: (1) dust/sand lofting based on initializing plume releases from visual inspection of the dust RGB imagery; and (2) dust/sand release and dispersion using the internal HYSPLIT algorithm that determines lofting from input wind speeds and prescribed land use. The goal is to develop a near real-time HYSPLIT modeling solution for routinely simulating the transport and dispersion of dust plumes across the HKH region. This presentation illustrates these products and capabilities for various use cases from the 2019 to 2021 period, highlighting a particularly unhealthy episode during late March 2021 across Nepal.

Remote Sensing↗

Oklahoma Health and Air Quality: Mapping Air Quality Using NASA Earth Observations to Investigate Recent Increases in Ozone Concentrations

Tropospheric ozone (O3) is formed by anthropogenic pollutants interacting with sunlight and is considered harmful to human health in high concentrations. In the summer of 2018, the Oklahoma Department of Environmental Quality (DEQ) measured unexpected spikes in O3 in Seiling, Oklahoma, with concentrations exceeding those measured in bustling Oklahoma City and Tulsa. The DEQ tracks air quality using ground monitors and does not utilize Earth observation data in its monitoring or analysis. This project used remotely sensed data to investigate these2018 air quality anomalies, identifying possible causes. We analyzed atmospheric data from Terra and Aqua Moderate Resolution Imaging Spectroradiometer (MODIS), and Sentinel-5P Tropospheric Ozone Monitoring Instrument (TROPOMI) in conjunction with ground-based measurements of tropospheric ozone (O3), nitrogen dioxide (NO2), methane (CH4), carbon monoxide (CO), formaldehyde (HCHO) and aerosol optical depth (AOD). We compared Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model simulations and Earth observation visualizations to pinpoint ozone spike causes. We also generated models to identify contributing factors to variations in ground ozone concentrations in our study area. The results point to a variety of ozone spike causes, primarily from outside of the state, and support the placement of additionalNO2, O3, and CO monitors to the southeast of Seiling. These analyses can help guide the placement of future monitors in the ground monitoring network and inform air quality regulations in Oklahoma.

Carolina Rosales↗

Analysis of Atmospheric Conditions Responsible for an Ozone Exceedance Event in Southeast Virginia on June 15, 2022

On June 15, 2022, the Virginia Department of Environmental Quality monitoring site at Suffolk/Holland, a rural site in southeast Virginia, recorded its first ozone exceedance since 2016. An ozone exceedance day occurs when the daily maximum 8-hour average surface ozone concentration is greater than 70 ppb. On this day, it was 75 ppb. This event is also noteworthy due to the rapid increase in surface ozone between 6 am and 9 am EDT as well as the hourly maximum ozone concentration of 82 ppb that was measured at 1 pm EDT. In this analysis, we utilize various observational and model data to determine the atmospheric conditions responsible for this ozone exceedance event. The analysis is conducted in two parts: (1) an evaluation of the accuracy of the GEOS-CF and WRF-Chem model forecasts and (2) an investigation of the mechanisms responsible for the high surface ozone at the Virginia DEQ Suffolk/Holland monitoring site. Comparisons of model forecasted ozone time-height cross sections with measured ozone by lidars from the Tropospheric Ozone Lidar Network (TOLNet) at NASA Langley and NASA Goddard reveal that both models forecasted ozone reasonably well near the surface. Model and observational data reveal that an upper-level ridge was present over much of the eastern United States, with a broad anticyclonic circulation near the surface, resulting in north-northeasterly flow over southeast Virginia. Back trajectory calculations using the Hybrid Single-Particle Lagrangian Integrated Trajectory model (HYSPLIT) and model data suggest that a low-level plume rich in ozone and precursor species originating from the New York City and New Jersey region was transported along the east coast between June 14 and June 15, arriving in southeast Virginia around 6 am EDT on June 15. In the early morning hours, this plume mixed down to the surface, elevating the background ozone mixing ratio as well as the mixing ratios of several precursor species. Other potential contributions are also explored and discussed. Lastly, the authors note that this work represents the analysis of the authors and not the Virginia Department of Environmental Quality.

Daniel B Phoenix↗

Expanding SPoRT RGBs and Machine Learning Techniques to Enhance Air Quality Monitoring in Southern Asia

Air pollution poses significant environmental, public health, and societal concerns in the Hindu Kush Himalaya (HKH) region of south-central Asia, notably during the dry monsoon months (~November to May). Key contributors to poor air quality include dust from the Middle East and western India, persistent nocturnal fog/smog, and biomass burning. To address this issue, we established a robust air quality and chemistry observation and modeling product suite utilizing multi-spectral red-green-blue (RGB) composite satellite products from Korea’s GEO-KOMPSAT-2A satellite, the Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model for dust transport forecasts, and the Weather Research and Forecasting coupled with Chemistry (WRF-Chem) model to predict aerosols and chemical species concentrations. Our team employed similar RGB recipes transitioned by the NASA Short-term Prediction Research and Transition (SPoRT) Center for the GOES-R era products over the Western Hemisphere, with significant success in depicting dust and nocturnal fog / low clouds, and to a lesser extent smoke and fire hot spots. We will extend these capabilities for the HKH region by applying an artificial intelligence (AI) model that objectively identifies dust from multi-spectral satellite data over the Southwestern United States. The AI model will be calibrated for automated dust detection over HKH from GEO-KOMPSAT-2A satellite data, with a goal of developing a similar AI model for objectively identifying smoke as well. The ultimate goal of this effort is to enhance dust and smoke predictions in the region by establishing improved emission initializations in HYSPLIT and/or WRF-Chem through the automated AI detection of dust and smoke.

Jonathan L. Case↗

Twistor Actions for Integrable Systems

Many integrable systems can be reformulated as holomorphic vector bundles on twistor space. This is a powerful organizing principle in the theory of integrable systems. One shortcoming is that it is formulated at the level of the equations of motion. From this perspective, it is mysterious that integrable systems have Lagrangians. In this paper, we study a Chern-Simons action on twistor space and use it to derive the Lagrangians of some integrable sigma models. Our focus is on examples that come from dimensionally reduced gravity and supergravity. The dimensional reduction of general relativity to two spacetime dimensions is an integrable coset sigma model coupled to a dilaton and 2d gravity. The dimensional reduction of supergravity to two spacetime dimensions is an integrable coset sigma model coupled to matter fermions, a dilaton, and 2d supergravity. We derive Lax operators and Lagrangians for these 2d integrable systems using the Chern-Simons theory on twistor space. In the supergravity example, we use an extended setup in which twistor Chern-Simons theory is coupled to a pair of matter fermions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A macroscopic plasma Lagrangian and its application to wave interactions and resonances

The derivation of a macroscopic plasma Lagrangian is considered, along with its application to the description of nonlinear three-wave interaction in a homogeneous plasma and linear resonance oscillations in a inhomogeneous plasma. One approach to obtain the Lagrangian is via the inverse problem of the calculus of variations for arbitrary first and second order quasilinear partial differential systems. Necessary and sufficient conditions for the given equations to be Euler-Lagrange equations of a Lagrangian are obtained. These conditions are then used to determine the transformations that convert some classes of non-Euler-Lagrange equations to Euler-Lagrange equation form. The Lagrangians for a linear resistive transmission line and a linear warm collisional plasma are derived as examples. Using energy considerations, the correct macroscopic plasma Lagrangian is shown to differ from the velocity-integrated low Lagrangian by a macroscopic potential energy that equals twice the particle thermal kinetic energy plus the energy lost by heat conduction.

Peng, Y. K. M.↗

An Eulerian/Lagrangian coupling procedure for three-dimensional vortical flows

A coupled Eulerian/Lagrangian method is presented for the reduction of numerical diffusion observed in solutions of 3D vortical flows using standard Eulerian finite-volume time-marching procedures. A Lagrangian particle tracking method, added to the Eulerian time-marching procedure, provides a correction of the Eulerian solution. In turn, the Eulerian solution is used to integrate the Lagrangian state-vector along the particles trajectories. While the Eulerian solution ensures the conservation of mass and sets the pressure field, the particle markers describe accurately the convection properties and enhance the vorticity and entropy capturing capabilities of the Eulerian solver. The Eulerian/Lagrangian coupling strategies are discussed and the combined scheme is tested on a constant stagnation pressure flow in a 90 deg bend and on a swirling pipe flow. As the numerical diffusion is reduced when using the Lagrangian correction, a vorticity gradient augmentation is identified as a basic problem of this inviscid calculation.

Felici, Helene M.↗

A coupled Eulerian/Lagrangian method for the solution of three-dimensional vortical flows

A coupled Eulerian/Lagrangian method is presented for the reduction of numerical diffusion observed in solutions of three-dimensional rotational flows using standard Eulerian finite-volume time-marching procedures. A Lagrangian particle tracking method using particle markers is added to the Eulerian time-marching procedure and provides a correction of the Eulerian solution. In turn, the Eulerian solutions is used to integrate the Lagrangian state-vector along the particles trajectories. The Lagrangian correction technique does not require any a-priori information on the structure or position of the vortical regions. While the Eulerian solution ensures the conservation of mass and sets the pressure field, the particle markers, used as 'accuracy boosters,' take advantage of the accurate convection description of the Lagrangian solution and enhance the vorticity and entropy capturing capabilities of standard Eulerian finite-volume methods. The combined solution procedures is tested in several applications. The convection of a Lamb vortex in a straight channel is used as an unsteady compressible flow preservation test case. The other test cases concern steady incompressible flow calculations and include the preservation of turbulent inlet velocity profile, the swirling flow in a pipe, and the constant stagnation pressure flow and secondary flow calculations in bends. The last application deals with the external flow past a wing with emphasis on the trailing vortex solution. The improvement due to the addition of the Lagrangian correction technique is measured by comparison with analytical solutions when available or with Eulerian solutions on finer grids. The use of the combined Eulerian/Lagrangian scheme results in substantially lower grid resolution requirements than the standard Eulerian scheme for a given solution accuracy.

Felici, Helene Marie↗

Using XR for Improving Scientific Discovery With Numerical Weather Models

Our work explores the use of extended reality (XR) to im- prove scientific discovery with numerical weather/climate models that inform Earth science digital twins, specifically the NASA Goddard Earth Observing System (GEOS) global atmospheric model. The overall project is named the Vi- sualization And Lagrangian dynamics Immersive eXtended Reality Toolkit (VALIXR), which has two main areas of focus: (1) enhancing the understanding of and interaction with model output data through advanced visualizations in the XR environment, and (2) the integration of Lagrangian dynamics into the GEOS model, which allows a natural, feature-specific analysis of Earth science phenomena as op- posed to traditional, fixed-point Eulerian dynamics. Here, we report initial work on these focus areas.

Thomas Grubb↗

A direct transcription-based multiple shooting formulation for dynamic optimization

The growing need for fast and efficient solution techniques for solving dynamic optimization problems is driven by a broad spectrum of applications in scheduling and control. We suggest a novel framework for dynamic optimization that utilizes a multiple shooting “backbone” with discrete rather than continuous subproblems, thereby eliminating need for repeated time-integration. A Lagrangian relaxation (LR)-based decomposition scheme is proposed, which dualizes the state continuity requirements between subproblems and enables parallel solution of the problem. We demonstrate the applicability of our method on two case studies: the Van der Pol oscillator and a batch reactor.

42 ENGINEERING↗

Application of Lagrangian techniques for calculating the on-axis rotational transform

The Floquet exponents of periodic field lines are studied through the variations of the magnetic action on the magnetic axis, which is assumed to be elliptical. The near-axis formalism developed by Mercier, Solov'ev and Shafranov is combined with a Lagrangian approach. The on-axis Floquet exponent is shown to coincide with the on-axis rotational transform. A discrete solution suitable for numerical implementation is introduced, which gives the Floquet exponents as solutions to an eigenvalue problem. This discrete formalism expresses the exponents as the eigenvalues of a $6$ X $6$ matrix.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗