Search NASA⌕ Search

SEARCH · Search NASA

Results for “Eulerian vertical advection”

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.

12 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2021

An Evaluation of Size-Resolved Cloud Microphysics Scheme Numerics for use with Radar Observations Part II: Condensation and Evaporation

Accurate numerical modeling of clouds and precipitation is essential for weather forecasting and climate change research. While size-resolved (bin) cloud microphysics models predict particle size distributions without imposing shapes, results are subject to artificial size distribution broadening owing to numerical diffusion associated with various processes. Whereas Part 1 addressed collision-coalescence, here we investigate numerical diffusion that occurs in solving condensation and evaporation. In a parcel model framework, all of the numerical schemes examined converge to one solution of condensation and evaporation as the mass grid is refined, and the advection-based schemes are recommended over the reassigning schemes. Including Eulerian vertical advection in a column limits the convergence to some extent, but that limitation occurs at a sufficiently fine mass grid, and the number of iterations in solving vertical advection should be minimized to reduce numerical diffusion. Insubstantial numerical diffusion in solving condensation can be amplified if collision-coalescence is also active, which in turn can be substantially diminished if turbulence effects on collision are incorporated. Large-eddy simulations of a drizzling stratocumulus field reveal that changes in moments of Doppler spectra obtained using different mass grids are consistent with those obtained from the simplified framework, and that spectral moments obtained using a mass grid designed to effectively reduce numerical diffusion are generally closer to observations. Notable differences between the simulations and observations still exist, and our results suggest a need to consider whether factors other than numerical diffusion in the fundamental process schemes employed can cause such differences.

Cloud microphysics↗
NASA NTRS2016

Hemispheric Differences in Tropical Lower Stratospheric Transport and Tracers Annual Cycle

Transport of long-lived tracers (such as O, CO, and N O) in the lower stratosphere largely determines the composition of the entire stratosphere. Stratospheric transport includes the mean residual circulation (with air rising in the tropics and sinking in the polar and middle latitudes), plus two-way isentropic (quasi-horizontal) mixing by eddies. However, the relative importance of two transport components remains uncertain. Previous studies quantified the relative role of these processes based on tropics-wide average characteristics under common assumption of well-mixed tropics. However, multiple instruments provide us with evidence that show significant differences in the seasonal cycle of ozone between the Northern (0-20N) and Southern (0-20S) tropical (NT and ST respectively) lower stratosphere. In this study we investigate these differences in tracer seasonality and quantify transport processes affecting tracers annual cycle amplitude using simulations from Goddard Earth Observing System Chemistry Climate Model (GEOSCCM) and Whole Atmosphere Community Climate Model (WACCM) and compare them to observations from the Microwave Limb Sounder (MLS) on the Aura satellite. We detect the observed contrast between the ST and NT in GEOSCCM and WACCM: annual cycle in ozone and other chemical tracers is larger in the NT than in the ST but opposite is true for the annual cycle in vertical advection. Ozone budgets in the models, analyzed based on the Transformed Eulerian Mean (TEM) framework, demonstrate a major role of quasi-horizontal mixing vertical advection in determining the NTST ozone distribution and behavior. Analysis of zonal variations in the NT and ST ozone annual cycles further suggests important role of North American and Asian Summer Monsoons (associated with strong isentropic mixing) on the lower stratospheric ozone in the NT. Furthermore, multi model comparison shows that most CCMs reproduce the observed characteristic of ozone annual cycle quite well. Thus, latitudinal variations within the tropics have to be considered in order to understand the balance between upwelling and quasi- horizontal mixing in the tropical lower stratosphere and the paradigm of well mixed tropics has to be reconsidered.

Chemistry-Climate Models↗
NASA NTRS1996

An Atmospheric Tape Recorder: The Imprint of Tropical Tropopause Temperatures on Stratospheric Water Vapor

We describe observations of tropical stratospheric water vapor q that show clear evidence of large-scale upward advection of the signal from annual fluctuations in the effective 'entry mixing ratio' q(sub E) of air entering the tropical stratosphere. In other words, air is 'marked,' on emergence above the highest cloud tops, like a signal recorded on an upward moving magnetic tape. We define q(sub E) as the mean water vapor mixing ratio, at the tropical tropopause, of air that will subsequently rise and enter the stratospheric 'overworld' at about 400 K. The observations show a systematic phase lag, increasing with altitude, between the annual cycle in q(sub E) and the annual cycle in q at higher altitudes. The observed phase lag agrees with the phase lag calculated assuming advection by the transformed Eulerian-mean vertical velocity of a q(sub E) crudely estimated from 100-hPa temperatures, which we use as a convenient proxy for tropopause temperatures. The phase agreement confirms the overall robustness of the calculation and strongly supports the tape recorder hypothesis. Establishing a quantitative link between q(sub E) and observed tropopause temperatures, however, proves difficult because the process of marking the tape depends subtly on both small- and large-scale processes. The tape speed, or large-scale upward advection speed, has a substantial annual variation and a smaller variation due to the quasi-biennial oscillation, which delays or accelerates the arrival of the signal by a month or two in the middle stratosphere. As the tape moves upward, the signal is attenuated with an e-folding time of about 7 to 9 months between 100 and 50 hPa and about 15 to 18 months between 50 and 20 hPa, constraining possible orders of magnitude both of vertical diffusion K(sub z) and of rates of mixing in from the extratropics. For instance, if there were no mixing in, then K(sub z) would be in the range 0.03-0.09 m(exp 2)/s; this is an upper bound on K(sub z).

Mote, Philip W.↗
NASA NTRS2017

Hemispheric Differences in the Annual Cycle of Tropical Lower Stratosphere Transport and Tracers

Transport in the tropical lower stratosphere plays a major role in determining the composition of the entire stratosphere. Previous studies that quantified the relative role of transport processes have generally assumed well-mixed tropics and focused on tropical-wide average characteristics. However, it has recently been shown that there is a hemispheric difference in the annual cycle of tropical lower stratosphere ozone and other tracers, with a larger amplitude in the northern tropics (NT) than in the southern tropics (ST). In this study, we examined the ability of chemistry climate models (CCMs) to reproduce the hemispheric differences in ozone (O3) and other tracers (i.e., hydrochloric acid, or HCl and nitrous oxide, or N2O), and then use the CCMs to examine the cause of these differences. Examination of CCM simulations from the CCMVal-2 project shows that the majority of the CCMs produce the observed feature of a larger annual cycle in the NT than ST O3 and other tracers. However, only around a third of the models produce an ozone annual cycle similar to that observed. Transformed Eulerian Mean analysis of two of the CCMs shows that seasonality in vertical advection drives the seasonality in ST O3 and N2O while seasonality of horizontal mixing drives the seasonality in NT O3 and N2O, with a large increase in horizontal mixing during northern summer (associated with the Asian monsoon). Thus, latitudinal and longitudinal variations within the tropics have to be considered to fully understand the balance between transport processes in tropical lower stratosphere.

transport↗
NASA NTRS2019

Antarctic Ice Shelf Thickness Change from Multimission Lidar Mapping

We calculate rates of ice thickness change and bottom melt for ice shelves in West Antarctica and the Antarctic Peninsula from a combination of elevation measurements from NASA–CECS Antarctic ice mapping campaigns and NASA Operation IceBridge corrected for oceanic processes from measurements and models, surface velocity measurements from synthetic aperture radar, and high-resolution outputs from regional climate models. The ice thickness change rates are calculated in a Lagrangian reference frame to reduce the effects from advection of sharp vertical features, such as cracks and crevasses, that can saturate Eulerian-derived estimates. We use our method over different ice shelves in Antarctica, which vary in terms of size, repeat coverage from airborne altimetry, and dominant processes governing their recent changes. We find that the Larsen-C Ice Shelf is close to steady state over our observation period with spatial variations in ice thickness largely due to the flux divergence of the shelf. Firn and surface processes are responsible for some short-term variability in ice thickness of the Larsen-C Ice Shelf over the time period. The Wilkins Ice Shelf is sensitive to short-timescale coastal and upper-ocean processes, and basal melt is the dominant contributor to the ice thickness change over the period. At the Pine Island Ice Shelf in the critical region near the grounding zone, we find that ice shelf thickness change rates exceed 40m/yr, with the change dominated by strong submarine melting. Regions near the grounding zones of the Dotson and Crosson ice shelves are decreasing in thickness at rates greater than 40m/yr, also due to intense basal melt. NASA–CECS Antarctic ice mapping and NASA Operation IceBridge campaigns provide validation datasets for floating ice shelves at moderately high resolution when coregistered using Lagrangian methods.

Sutterley, Tyler C.↗
NASA NTRS2018

The NASA Eulerian Snow on Sea Ice Model (NESOSIM) v1.0: initial model development and analysis

The NASA Eulerian Snow On Sea Ice Model (NESOSIM) is a new, open-source snow budget model that is currentlyconfigured to produce daily estimates of the depth and density of snow on sea ice across the Arctic Ocean throughthe accumulation season. NESOSIM has been developed in a three-dimensional Eulerian framework and includes two(vertical) snow layers and several simple parameterizations (accumulation, wind packing, advection_divergence, blowingsnow lost to leads) to represent key sources and sinks of snow on sea ice. The model is forced with daily inputs of snowfall and near-surface winds (from reanalyses), sea ice concentration (from satellite passive microwave data) and sea ice drift (from satellite feature tracking) during the accumulation season (August through April). In this study, we present the NESOSIM formulation, calibration efforts, sensitivity studies and validation efforts across an Arctic Ocean domain (100 km horizontal resolution). The simulated snow depth and density are calibrated with in situ data collected on drifting ice stations during the 1980s. NESOSIM shows strong agreement with the in situ seasonal cycles of snow depth and density, and shows good (moderate) agreement with the regional snow depth (density) distributions. NESOSIM is run for a contemporary period (2000 to 2015), with the results showing strong sensitivity to the reanalysis derived snowfall forcing data, with the Modern-Era Retrospective analysis for Research and Applications (MERRA) and the Japanese Meteorological Agency 55-year reanalysis (JRA-55) forced snow depths generally higher than ERA-Interim,and the Arctic System Reanalysis (ASR) generally lower. We also generate and force NESOSIM with a consensus "median" daily snowfall dataset from these reanalyses.The results are compared against snow depth estimates derived from NASA's Operation IceBridge (OIB) snow radar data from 2009 to 2015, showing moderate_strong correlations and root mean squared errors of 10 cm depending on the OIB snow depth product analyzed, similar to the comparisons between OIB snow depths and the commonly used modified Warren snow depth climatology. Potential improvements to this initial NESOSIM formulation are discussed in the hopes of improving the accuracy and reliability of these simulated snow depths and densities.

Petty, Alek A.↗
NASA NTRS1981

An objective method for determining the generalized transport tensor for two-dimensional Eulerian models

An objective method for deriving the components of a generalized transport tensor for a two-dimensional model is presented. Representative meridional and vertical velocities and thermodynamic scalars at a uniform grid are used to reduce the problem to the solution of two flux equations for two unknowns. One unknown is the stream-function, coefficient of an antisymmetric tensor, which corrects the Eulerian mean motions for Stokes drift; the other is a time constant, which converts the deviatory velocity tensor to a symmetric transport tensor. The complete asymmetric tensor, called a transport tensor, has a divergence which yields both advection and diffusion by the deviatory velocities. Advantages and disadvantages of Lagrangian and Eulerian averages are discussed, and meridional-vertical velocity correlations are provided.

Danielsen, E. F.↗
NASA NTRS2000

Sensitivity of Age-of-Air Calculations to the Choice of Advection Scheme

The age of air has recently emerged as a diagnostic of atmospheric transport unaffected by chemical parameterizations, and the features in the age distributions computed in models have been interpreted in terms of the models' large-scale circulation field. This study shows, however, that in addition to the simulated large-scale circulation, three-dimensional age calculations can also be affected by the choice of advection scheme employed in solving the tracer continuity equation, Specifically, using the 3.0deg latitude X 3.6deg longitude and 40 vertical level version of the Geophysical Fluid Dynamics Laboratory SKYHI GCM and six online transport schemes ranging from Eulerian through semi-Lagrangian to fully Lagrangian, it will be demonstrated that the oldest ages are obtained using the nondiffusive centered-difference schemes while the youngest ages are computed with a semi-Lagrangian transport (SLT) scheme. The centered- difference schemes are capable of producing ages older than 10 years in the mesosphere, thus eliminating the "young bias" found in previous age-of-air calculations. At this stage, only limited intuitive explanations can be advanced for this sensitivity of age-of-air calculations to the choice of advection scheme, In particular, age distributions computed online with the National Center for Atmospheric Research Community Climate Model (MACCM3) using different varieties of the SLT scheme are substantially older than the SKYHI SLT distribution. The different varieties, including a noninterpolating-in-the-vertical version (which is essentially centered-difference in the vertical), also produce a narrower range of age distributions than the suite of advection schemes employed in the SKYHI model. While additional MACCM3 experiments with a wider range of schemes would be necessary to provide more definitive insights, the older and less variable MACCM3 age distributions can plausibly be interpreted as being due to the semi-implicit semi-Lagrangian dynamics employed in the MACCM3. This type of dynamical core (employed with a 60-min time step) is likely to reduce SLT's interpolation errors that are compounded by the short-term variability characteristic of the explicit centered-difference dynamics employed in the SKYHI model (time step of 3 min). In the extreme case of a very slowly varying circulation, the choice of advection scheme has no effect on two-dimensional (latitude-height) age-of-air calculations, owing to the smooth nature of the transport circulation in 2D models. These results suggest that nondiffusive schemes may be the preferred choice for multiyear simulations of tracers not overly sensitive to the requirement of monotonicity (this category includes many greenhouse gases). At the same time, age-of-air calculations offer a simple quantitative diagnostic of a scheme's long-term diffusive properties and may help in the evaluation of dynamical cores in multiyear integrations. On the other hand, the sensitivity of the computed ages to the model numerics calls for caution in using age of air as a diagnostic of a GCM's large-scale circulation field.

Eluszkiewicz, Janusz↗
NASA NTRS1993

Solution of mixed convection heat transfer from isothermal in-line fins

Transient and steady state combined natural and forced convective flows over two in-line finite thickness fins (louvers) in a vertical channel are numerically solved using two methods. The first method of solution is based on the 'Simple Arbitrary Lagrangian Eulerian' (SALE) technique which incorporates mainly two computational phases: (1) a Lagrangian phase in which the velocity field is updated by the effects of all forces, and (2) an Eulerian phase that executes all advective fluxes of mass, momentum and energy. The second method of solution uses the finite element code entitled FIDAP. In the first part, comparison of the results by FIDAP, SALE, and available experimental work were done and discussed for steady state forced convection over louvered fins. Good agreements were deduced between the three sets of results especially for the flow over a single fin. In the second part and in the absence of experimental literature, the numerical predictions were extended to the transient transports and to the opposing flow where pressure drop is reversed. Results are presented and discussed for heat transfer and pressure drop in assisting and opposing mixed convection flows.

Khalilollahi, Amir↗
NASA NTRS1993

Eddy mixing coefficient values derived from simulations with the NASA Ames Mars GCM

Values of eddy mixing coefficients, especially that for vertical mixing, are of particular importance for various photochemical models of the Mars atmosphere, which are either globally or zonally averaged. These models represent atmospheric transport processes in terms of eddy diffusion. While this is not appropriate for the advective transport by the mean meridional circulation (in general, a part of the eddy diffusion tensor actually does represent advection--a 'correction' to the Eulerian-mean circulation), it can be applied to simulate both the small- and larger-scale eddy transports of an atmospheric constituent. An effort is underway to estimate values of the eddy mixing coefficients for the Mars atmosphere using circulation data generated with the NASA Ames Mars GCM. This model simulates the three-dimensional winds in the atmosphere, and these are then used as inputs for a tracer transport model. The latter model has previously been used in interactive dust transport experiments with the Ames GCM. Idealized tracer transport experiments, with a conservative or nearly conservative hypothetical tracer, are performed to generate data from which the eddy diffusion coefficients can be estimated. Carrying out matched pairs of tracer experiments, using two very different initial states, permits all four components of the diffusion tensor to be determined. A large number of GCM experiments have been conducted, spanning a range of seasons and atmospheric dust loading, allowing the eddy mixing coefficients to be estimated for a variety of atmospheric conditions.

Barnes, J. R.↗
NASA NTRS1977

Development and application of a three dimensional numerical model for predicting pollutant and sediment transport using an Eulerian-Lagrangian marker particle technique

A computer coded Lagrangian marker particle in Eulerian finite difference cell solution to the three dimensional incompressible mass transport equation, Water Advective Particle in Cell Technique, WAPIC, was developed, verified against analytic solutions, and subsequently applied in the prediction of long term transport of a suspended sediment cloud resulting from an instantaneous dredge spoil release. Numerical results from WAPIC were verified against analytic solutions to the three dimensional incompressible mass transport equation for turbulent diffusion and advection of Gaussian dye releases in unbounded uniform and uniformly sheared uni-directional flow, and for steady-uniform plug channel flow. WAPIC was utilized to simulate an analytic solution for non-equilibrium sediment dropout from an initially vertically uniform particle distribution in one dimensional turbulent channel flow.

Pavish, D. L.↗
NASA NTRS2023

Assessment of Oceanographic Conditions during the North Atlantic EXport Processes in the Ocean from RemoTe Sensing (EXPORTS) Field Campaign

This manuscript presents an overview of NASA’s EXport Processes in the Ocean from Remote Sensing 2021 Field Campaign in the North Atlantic (EXPORTS NA) and provides quantitative and dynamical descriptions of the physical processes modulating water transformations during the study. A major programmatic goal was to conduct the sampling in a Lagrangian mode so that ocean ecological and biogeochemical changes can be observed independent from physical advective processes. To accomplish this goal, EXPORTS NA conducted a multi-ship, multi-asset field sampling program within a retentive, anticyclonic mode water eddy. Beneath depths of ∼100 m, Lagrangian sampling assets remained within the eddy core waters (ECWs) throughout the experiment, demonstrating that the ECWs within the mode water eddy were retentive. However, strong westerly winds from four storm events deepened the mixed layer (ML) of the surface core waters (SCWs) above the eddy’s mode water core by 25-40 m and exchanged some of the SCWs with surface waters outside of the eddy via Ekman transport. Estimates of flushing times ranged from 5-8 days, with surface exchange fractions ranging from 20-75% and were consistent with particle tracking advected by combined geostrophic and Ekman velocities. The relative contributions of horizontal and vertical advection on changes in ECW tracers depended on the horizontal and vertical gradients of that tracer. For example, horizontal advection played a large role in ECW salinity fluxes, while vertical entrainment played a larger role in the fluxes of nutrients into SCW ML. Each storm injected nutrients and low oxygen waters into the ML, after which the surface ocean ecosystem responded by reducing nutrient concentrations and increasing %O2 saturation levels. Overall, ECW values of chlorophyll and POC were the largest at the onset of the field program and decreased throughout the campaign. The analysis presented provides a physical oceanographic context for the many measurements made during the EXPORTS NA field campaign while illustrating the many challenges of conducting a production-flux experiment, even in a Lagrangian frame, and the inherent uncertainties of interpreting biological carbon pump observations that were collected in a Eulerian frame of reference.

oceanography↗