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 19 records

Semi-Lagrangian integration of a grid-point shallow water model on the sphere

This paper describes a semi-Lagrangian technique for integrating the equations of motion on the global domain. The technique uses an auxiliary spherical coordinate system at each near-polar gridpoint of the latitude-longitude grid; the auxiliary system is obtained by a rotation such that the new equator passes through the gridpoint in question and the new coordinate directions coincide with those of the original system at that point. The technique was applied to the shallow water equations, incorporating a semiimplicit treatment of the adjustment terms on a C-grid, with two-time levels. A five day integration was successfully carried out for a situation involving strong cross-polar flow. No filtering or diffusion was required to maintain stability over a five day period.

Mcdonald, A.↗

Semi-Lagrangian integration of a gridpoint shallow water model on the sphere

A stable, semi-Lagrangian, semi-implicit, two-time-level, gridpoint integration scheme for the shallow water equations on the sphere is presented. A rotated spherical coordinate system is used to integrate the equations of motion at each gridpoint poleward of a certain latitude, thus overcoming problems associated with the polar singularity. The results of medium term integrations of large scale test patterns using a long time step are presented.

Mcdonald, A.↗

Southern California Megacity CO2, CH4, and CO Flux Estimates Using Ground- and Space-Based Remote Sensing and a Lagrangian Model

We estimate the overall CO2, CH4, and CO flux from the South Coast Air Basin using an inversion that couples Total Carbon Column Observing Network (TCCON) and Orbiting Carbon Observatory-2 (OCO-2) observations, with the Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) model and the Open-source Data Inventory for Anthropogenic CO2 (ODIAC). Using TCCON data we estimate the direct net CO2 flux from the So-CAB to be 104±26 Tg CO2 yr(exp -1) for the study period of July 2013–August 2016. We obtain a slightly higher estimate of 120±30 Tg CO2 yr(exp -1) using OCO-2 data. These CO2 emission estimates are on the low end of previous work. Our net CH4 (360±90 Gg CH4 y(exp -1)) flux estimate is in agreement with central values from previous top-down studies going back to 2010 (342–440 Gg CH4 yr(exp -1)). CO emissions are estimated at 487±122 Gg CO yr(exp -1), much lower than previous top-down estimates (1440 Gg CO yr(exp -1)). Given the decreasing emissions of CO, this finding is not unexpected. We perform sensitivity tests to estimate how much errors in the prior, errors in the covariance, different inversion schemes, or a coarser dynamical model influence the emission estimates. Overall, the uncertainty is estimated to be 25%, with the largest contribution from the dynamical model. Lessons learned here may help in future inversions of satellite data over urban areas.

Total Carbon Column Observing Network (TCCON)↗

Air Pollution Scenario over Pakistan: Characterization and Ranking of Extremely Polluted Cities using Long-Term Concentrations of Aerosols and Trace Gases

Pakistan ranks third in the world in terms of mortality attributable to air pollution, with aerosol mass concentrations (PM2.5) consistently well above WHO (World Health Organization) air quality guidelines (AQG). However, regulation is dependent on a sparse network of air quality monitoring stations and insufficient ground data. This study utilizes long-term observations of aerosols and trace gases to characterize and rank the air pollution scenarios and pollution characteristics of 80 selected cities in Pakistan. Datasets used include (1) the Aqua and Terra (AquaTerra) MODIS (Moderate Resolution Imaging Spectroradiometer) Level 2 Collection 6.1 merged Dark Target and Deep Blue (DTB) aerosol optical depth (AOD) retrieval products; (2) the CAMS (Copernicus Atmosphere Monitoring Service) reanalysis PM1, PM2.5, and PM10 data; (3) the MERRA-2 (Modern-Era Retrospective analysis for Research and Applications, Version 2) reanalysis PM2.5 data, (4) the OMI (Ozone Monitoring Instrument) tropospheric vertical column density (TVCD) of nitrogen dioxide (NO2), and VCD of sulfur dioxide (SO2) in the Planetary Boundary Layer (PBL), (5) the VIIRS (Visible Infrared Imaging Radiometer Suite) Nighttime Lights data, (6) MODIS Collection 6 Version 2 global monthly fire location data (MCD14ML), (7) population density, (8) MODIS Level 3 Collection 6 land cover types, (9) AERONET (AErosol RObotic NETwork) Version 3 Level 2.0 data, and (10) ground-based PM2.5 concentrations from air quality monitoring stations. Potential Source Contribution Function (PSCF) analyses were performed by integrating with ground-based PM2.5 concentrations and the NOAA (National Oceanic and Atmospheric Administration) HYSPLIT (Hybrid Single-Particle Lagrangian Integrated Trajectory) air parcel back trajectories to identify potential pollution source areas which are responsible for extreme air pollution in Pakistan. Results show that the ranking of the top polluted cities depends on the type of pollutant considered and the metric used. For example, Jhang, Multan, and Vehari were characterized as the top three polluted cities in Pakistan when considering AquaTerra DTB AOD products; for PM1, PM2.5, and PM10 Lahore, Gujranwala, and Okara were the top three; for tropospheric NO2 VCD Lahore, Rawalpindi, and Islamabad and for PBL SO2 VCD Lahore, Mirpur, and Gujranwala. The results demonstrate that Pakistan’s entire population has been exposed to high PM2.5 concentrations for many years, with a mean annual value of 54.7 μg/cu. m, over all Pakistan from 2003 to 2020. This value exceeds Pakistan’s National Environmental Quality Standards (Pak-NEQS, i.e., <15 μg/cu. m annual mean) for ambient air defined by the Pakistan Environmental Protection Agency (Pak-EPA) as well as the WHO Interim Target-1 (i.e., mean annual PM2.5 <35 μg/cu. m). The spatial analyses of the concentrations of aerosols and trace gases in terms of population density, nighttime lights, land cover types, and fire location data, and the PSCF analysis indicate that Pakistan’s air quality is strongly affected by anthropogenic sources inside of Pakistan, with contributions from surrounding countries. Statistically significant positive (increasing) trends in PM1, PM2.5, PM10, tropospheric NO2 VCD, and SO2 VCD were observed in ~89%, ~67%, ~48%, 91%, and ~88% of the Pakistani cities (80 cities), respectively. This comprehensive analysis of aerosol and trace gas levels, their characteristics in spatio-temporal domains, and their trends over Pakistan, is the first of its kind. Results will be helpful to the Ministry of Climate Change (Government of Pakistan), Pak-EPA, SUPARCO (Pakistan Space and Upper Atmosphere Research Commission), policymakers, and the local research community to mitigate air pollution and its effects on human health.

Muhammad Bilal↗

Simulation of cryogenic liquid flows with vapor bubbles

Liquid flows in rocket engine components (such as bearings, seals, and pumps) often involve the formation of vapor bubbles due to local superheating of the fluid (either boiling or cavitation). Under the present effort, an analysis has been developed for liquid flows with vapor bubbles, based on a combined Eulerian-Lagrangian technique, in which the continuous (liquid) phase is treated by solving a system of Eulerian conservation equations, while the discrete (vapor bubble) phase is dealt with by integrating Lagrangian equations of motion in computational coordinates. Vapor bubbles of changing size can be accommodated easily by this analysis, and models for the simulation of bubble formation, growth, and motion have been included. The effect of bubble motion and other bubble processes on the continuous (liquid) phase has been accounted for by appropriate bubble mass, momentum, and energy interchange source terms in the Eulerian conservation equations. To demonstrate the viability of the resulting procedure, the cavitating flow of liquid oxygen through a simplified model of a labyrinth seal has been successfully calculated.

De Jong, Frederik J.↗

Turbulent transport measurements with a laser Doppler velocimeter

The power spectrum of phototube current from a laser Doppler velocimeter operating in the heterodyne mode has been computed. The spectrum is obtained in terms of the space time correlation function of the fluid. The spectral width and shape predicted by the theory are in agreement with experiment. For normal operating parameters the time average spectrum contains information only for times shorter than the Lagrangian integral time scale of the turbulence. To examine the long time behavior, one must use either extremely small scattering angles, much longer wavelength radiation or a different mode of signal analysis, e.g., FM detection.

Edwards, R. V.↗

Turbulent transport measurements with a laser Doppler velocimeter.

The power spectrum of phototube current from a laser Doppler velocimeter operating in the heterodyne mode has been computed. The spectral width and shape predicted by the theory are in agreement with experiment. For normal operating parameters the time-average spectrum contains information only for times shorter than the Lagrangian-integral time scale of the turbulence. To examine the long-time behavior, one must use either extremely small scattering angles, much-longer-wavelength radiation, or a different mode of signal analysis, e.g., FM detection.

Edwards, R. V.↗

Repeated-cascade theory of strong turbulence in a magnetized plasma

A two-dimensional Navier-Stokes equation of vorticity in fluid turbulence is used to model drift turbulence in a plasma with a strong constant magnetic field and a constant mean density gradient. The nonlinear eddy diffusivity is described by a time-integrated Lagrangian correlation of velocities, and the repeated-cascade method is employed to choose the rank accounting for nearest-neighbor interactions, to calculate the Lagrangian correlation, and to close the correlation hierarchy. As a result, the diffusivity becomes dependent on the plasma's induced diffusion and is represented by a memory chain that is cut off by similarity and inertial randomization. Spectral laws relating the kinetic-energy spectrum to the -5, -5/2, -3, and -11 powers of wavenumber are derived for the velocity subranges of production, approach to inertia, inertia, and dissipation, respectively. It is found that the diffusivity is proportional to some inverse power of the magnetic field, that power being 1, 2/3, 5/6, and 2, respectively, for the four velocity subranges.

Tchen, C. M.↗

Rain estimation from satellites - Areal rainfall-rain area relation

The relationship between the areal rainfall rate and the rain area was investigated using rainfall rate data collected in the so-called B-scale during the first two phases of GATE. The data for the first two phases, each lasting about three weeks, are based mainly on the C-band radar observations carried out on board a research vessel at the center of the B-scale area 200 km in diameter, which is located at roughly 8.5 N, 23.5 W. It is demonstrated that the areal average rainfall rate correlates significantly with the rain area if the area under consideration is large compared to the Lagrangian integral scale.

Chiu, Long S.↗

Configuring the HYSPLIT Model for National Weather Service Forecast Office and Spaceflight Meteorology Group Applications

For expedience in delivering dispersion guidance in the diversity of operational situations, National Weather Service Melbourne (MLB) and Spaceflight Meteorology Group (SMG) are becoming increasingly reliant on the PC-based version of the HYSPLIT model run through a graphical user interface (GUI). While the GUI offers unique advantages when compared to traditional methods, it is difficult for forecasters to run and manage in an operational environment. To alleviate the difficulty in providing scheduled real-time trajectory and concentration guidance, the Applied Meteorology Unit (AMU) configured a Linux version of the Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) (HYSPLIT) model that ingests the National Centers for Environmental Prediction (NCEP) guidance, such as the North American Mesoscale (NAM) and the Rapid Update Cycle (RUC) models. The AMU configured the HYSPLIT system to automatically download the NCEP model products, convert the meteorological grids into HYSPLIT binary format, run the model from several pre-selected latitude/longitude sites, and post-process the data to create output graphics. In addition, the AMU configured several software programs to convert local Weather Research and Forecast (WRF) model output into HYSPLIT format.

Dreher, Joseph G.↗

Applied Meteorology Unit (AMU)

This report summarizes the Applied Meteorology Unit (AMU) activities for the first quarter of Fiscal Year 2010 (October - December 2009). A detailed project schedule is included in the Appendix. Included tasks are: (1) Peak Wind Tool for User Launch Commit Criteria (LCC), (2) Objective Lightning Probability Tool, Phase III, (3) Peak Wind Tool for General Forecasting, Phase II, (4) Upgrade Summer Severe Weather Tool in Meteorological Interactive Data Display System (MIDDS), (5) Advanced Regional Prediction System (ARPS) Data Analysis System (ADAS) Update and Maintainability, (5) Verify 12-km resolution North American Model (MesoNAM) Performance, and (5) Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) Graphical User Interface.

Bauman, William↗

Applied Meteorology Unit (AMU) Quarterly Report - Fourth Quarter FY-09

This report summarizes the Applied Meteorology Unit (AMU) activities for the fourth quarter of Fiscal Year 2009 (July - September 2009). Tasks reports include: (1) Peak Wind Tool for User Launch Commit Criteria (LCC), (2) Objective Lightning Probability Tool. Phase III, (3) Peak Wind Tool for General Forecasting. Phase II, (4) Update and Maintain Advanced Regional Prediction System (ARPS) Data Analysis System (ADAS), (5) Verify MesoNAM Performance (6) develop a Graphical User Interface to update selected parameters for the Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLlT)

Bauman, William↗

Configuring the HYSPLIT Model for National Weather Service Forecast Office and Spaceflight Meteorology Group Applications

The National Weather Service Forecast Office in Melbourne, FL (NWS MLB) is responsible for providing meteorological support to state and county emergency management agencies across East Central Florida in the event of incidents involving the significant release of harmful chemicals, radiation, and smoke from fires and/or toxic plumes into the atmosphere. NWS MLB uses the National Oceanic and Atmospheric Administration Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model to provide trajectory, concentration, and deposition guidance during such events. Accurate and timely guidance is critical for decision makers charged with protecting the health and well-being of populations at risk. Information that can describe the geographic extent of areas possibly affected by a hazardous release, as well as to indicate locations of primary concern, offer better opportunity for prompt and decisive action. In addition, forecasters at the NWS Spaceflight Meteorology Group (SMG) have expressed interest in using the HYSPLIT model to assist with Weather Flight Rules during Space Shuttle landing operations. In particular, SMG would provide low and mid-level HYSPLIT trajectory forecasts for cumulus clouds associated with smoke plumes, and high-level trajectory forecasts for thunderstorm anvils. Another potential benefit for both NWS MLB and SMG is using the HYSPLIT model concentration and deposition guidance in fog situations.

Dreher, Joseph↗

Configuring a Graphical User Interface for Managing Local HYSPLIT Model Runs Through AWIPS

Responding to incidents involving the release of harmful airborne pollutants is a continual challenge for Weather Forecast Offices in the National Weather Service. When such incidents occur, current protocol recommends forecaster-initiated requests of NOAA's Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model output through the National Centers of Environmental Prediction to obtain critical dispersion guidance. Individual requests are submitted manually through a secured web site, with desired multiple requests submitted in sequence, for the purpose of obtaining useful trajectory and concentration forecasts associated with the significant release of harmful chemical gases, radiation, wildfire smoke, etc., into local the atmosphere. To help manage the local HYSPLIT for both routine and emergency use, a graphical user interface was designed for operational efficiency. The interface allows forecasters to quickly determine the current HYSPLIT configuration for the list of predefined sites (e.g., fixed sites and floating sites), and to make any necessary adjustments to key parameters such as Input Model. Number of Forecast Hours, etc. When using the interface, forecasters will obtain desired output more confidently and without the danger of corrupting essential configuration files.

Wheeler, mark M.↗

Correlation of DIAL Ozone Observations with Lightning

The purpose of this project is to see whether ozone maxima measured by the DIfferential Absorption Lidar (DIAL) instrument in Huntsville, AL may be traced back to lightning events occurring 24‐ 48 hours beforehand. The methodology is to start with lidar measurements of ozone from DIAL as well as ozonesonde measurements. The HYbrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) model is then used to determine the origin of these ozone maxima 24‐48 hours prior. Data from the National Lightning Detection Network (NLDN) are used to examine the presence/absence of lightning along the trajectory. This type of analysis suggests that lightning‐produced NOx may be responsible for some of the ozone maxima over Huntsville.

Peterson, Harold↗

Correlation of DIAL Ozone Observations with Lightning

The purpose of this project is to see whether ozone maxima measured by the DIfferential Absorption Lidar (DIAL) instrument in Huntsville, AL may be traced back to lightning events occurring 24‐ 48 hours beforehand. The methodology is to start with lidar measurements of ozone from DIAL as well as ozonesonde measurements. The HYbrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) model is then used to determine the origin of these ozone maxima 24‐48 hours prior. Data from the National Lightning Detection Network (NLDN) are used to examine the presence/absence of lightning along the trajectory. This type of analysis suggests that lightning‐produced NOx may be responsible for some of the ozone maxima over Huntsville.

Peterson, Harold S.↗

Correlation of DIAL Ozone Observations with Lightning

The purpose of this project is to see whether ozone maxima measured by the DIfferential Absorption Lidar (DIAL) instrument in Huntsville, AL may be traced back to lightning events occurring 24‐48 hours beforehand. The methodology is to start with lidar measurements of ozone from DIAL. The HYbrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) model is then used to determine the origin of these ozone maxima 24‐48 hours prior. Data from the National Lightning Detection Network (NLDN) are used to examine the presence/absence of lightning along the trajectory. This type of analysis suggests that lightning‐produced NOx may be responsible for some of the ozone maxima over Huntsville.

Peterson, Harold↗

Correlation of DIAL Ozone Observations with Lightning

The purpose of this project is to see whether ozone maxima measured by the DIfferential Absorption Lidar (DIAL) instrument in Huntsville, AL may be traced back to lightning events occurring 24‐48 hours beforehand. The methodology is to start with lidar measurements of ozone from DIAL. The HYbrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) model is then used to determine the origin of these ozone maxima 24‐48 hours prior. Data from the National Lightning Detection Network (NLDN) are used to examine the presence/absence of lightning along the trajectory. This type of analysis suggests that lightning‐produced NOx may be responsible for some of the ozone maxima over Huntsville.

Peterson, Harold S.↗