Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Tracer equation”

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 37 records · Page 2

On the maintenance of the Marine boundary layer-free tropospheric gradient of transient tracers

A simple model of the marine boundary layer is used to arrive at a steady-state budget for trace gases in this layer. Five processes maintain the tracer balance in the layer. The first is a large-scale horizontal convergence of the tracer; the second is the pumping by clouds of the tracer out of the boundary layer into the free troposphere; the third is penetration by the tracer from the free troposphere into the marine boundary from compensating subsidence in the cloud-free air surrounding cumulus clouds; the fourth is sea surface emission or deposition of the tracer; and the fifth is in situ photochemical production or loss of the tracer within the marine boundary layer. The budget equation is used in deriving an expression for the mixing ratio difference between the two regions. A scale analysis suggests that the mixing ratio difference between the marine boundary layer and the free troposphere can be an order of magnitude smaller in disturbed weather conditions than in suppressed conditions, showing that the difference in chemical concentrations between the free troposphere and the marine boundary layer cannot be rationalized solely on the basis of chemical measurements.

Gidel, L. T.↗

A perspective of Middle-Atmosphere Dynamics (MAD) studies at the New International Equatorial Observatory (NIEO)

The equatorial region has attracted many MAD studies mainly based on data of limited locations and resolutions. Established at NIEO are: (1) Climatology of the equatorial middle atmosphere (all of the mean zonal flow, the meridional and/or east-west circulations and the planetary/gravity waves are described based on massive, reliable data statistics); (2) Troposphere-stratosphere coupling at the equator (the candidate location of NIEO is just at the stratospheric fountain area where the tracers and waves are pumped up into the middle atmosphere); and (3) Mesosphere-thermosphere coupling at the equator; thermospheric superrotation, which may be caused either by ion drag or by tidal breaking, is examined in detail by observations covering a wide altitude range from the mesosphere through the thermosphere.

Yamanaka, M. D.↗

Transport of atmospheric tracers by planetary waves during a winter stratospheric warming event - A three-dimensional model simulation

In a simple three-dimensional primitive equation model, a wave number 1 major stratospheric warming is simulated. With the aid of two idealized tracers it is shown that the transport during a major warming event is characterized by a small, well-organized tongue of subtropical air flowing around the displaced winter vortex into the polar cap and by a wide area with strong quasi-horizontal mixing (surf zone). The description of these dynamical processes requires a full three-dimensional space resolution.

Kouker, W.↗

Hydrologic Connectivity and Patch‐To‐Hillslope Scale Relations in Dryland Ecosystems

Abstract In drylands, runoff during storms redistributes water and nutrients from bare soil areas to vegetated patches, subsidizing vegetation with additional resources. The extent of this redistribution depends on the interplay between surface roughness and permeability; greater permeability in vegetated patches promotes run‐on to vegetation, but greater surface roughness diverts runoff, producing tortuous flow paths that bypass vegetation. Here, this interplay is examined in virtual experiments using the 2D Saint Venant Equations to measure runoff connectivity. Flowpaths are delineated using tracers advected by the flow. Distances between tracer sources and sinks along flowpaths measure hydrologic connectivity at two lengthscales: connectivity to the hillslope outlet and within‐slope source‐sink connectivity. Differences between these connectivity lengthscales indicate how flow may “by‐pass” vegetated patches within hillslopes. At the hillslope scale, a derived power‐law relation between the runoff coefficient and outlet connectivity describes hillslope water losses, providing a foundation for identifying landscapes likely to shed water.

54 ENVIRONMENTAL SCIENCES↗

Investigation of Tropical Transport with UARS Data

Measurements of trace constituents obtained by instruments aboard the Upper Atmosphere Research Satellite (UARS) have been used to study transport processes associated with the quasi-biennial oscillation, laterally propagating Rossby waves, and upward propagating Kelvin waves in the tropical and subtropical upper troposphere and stratosphere. Mean vertical motions, vertical diffusivities and in-mixing rates were inferred from observations of the 'tape recorder' signal in near-equatorial stratospheric water vapor. The effect of the quasi-biennial oscillation (QBO) on tracer distributions in the upper half of the stratosphere was seen in a spectacular 'staircase' pattern, predominantly in the winter hemisphere, revealing the latitudinally asymmetric nature of QBO transport due to induced mean meridional circulations and modulation of lateral mixing associated with planetary Rossby waves. The propagation of Rossby waves across the equator in the westerly phase of the QBO was seen in tracer fields and corroborating United Kingdom Meteorological Office (UKMO) analyses; a modeling study of the effect of these waves on typical QBO wind profiles was performed. Water vapor in the upper troposphere and lower stratosphere was found to exhibit signatures of the tropical intraseasonal oscillation (TIO) and faster Kelvin waves in the two regions, respectively.

Dunkerton, Timothy J.↗

A comparison of two- and three-dimensional tracer transport within a stratospheric circulation model

Use of the residual circulation for stratospheric tracer transport has been compared to a fully three-dimensional calculation. The wind fields used in this study were obtained from a global, semispectral, primitive equation model, extending from 10 to 100 km in altitude. Comparisons were done with a passive tracer and an ozone-like substance over a two-month period corresponding to a Northern Hemisphere winter. It was found that the use of the residual circulation can lead to errors in the tracer concentrations of about a factor of 2. The error is made up of two components. One is fluctuating with a period of approximately one month and reflects directly the wave transience that occurs on that time-scale. The second part is increasing steadily over the integration period and results from an overestimate of the vertical transport by the residual circulation. Furthermore, the equatorward and upward mixing that occurs with transport by the three-dimensional circulation at low latitudes is not well reproduced when the residual circulation is used.

Schneider, H.-R.↗

A Wavelet based Suboptimal Kalman Filter for Assimilation of Stratospheric Chemical Tracer Observations

A suboptimal Kalman filter system which evolves error covariances in terms of a truncated set of wavelet coefficients has been developed for the assimilation of chemical tracer observations of CH4. This scheme projects the discretized covariance propagation equations and covariance matrix onto an orthogonal set of compactly supported wavelets. Wavelet representation is localized in both location and scale, which allows for efficient representation of the inherently anisotropic structure of the error covariances. The truncation is carried out in such a way that the resolution of the error covariance is reduced only in the zonal direction, where gradients are smaller. Assimilation experiments which last 24 days, and used different degrees of truncation were carried out. These reduced the covariance size by 90, 97 and 99 % and the computational cost of covariance propagation by 80, 93 and 96 % respectively. The difference in both error covariance and the tracer field between the truncated and full systems over this period were found to be not growing in the first case, and growing relatively slowly in the later two cases. The largest errors in the tracer fields were found to occur in regions of largest zonal gradients in the constituent field. This results indicate that propagation of error covariances for a global two-dimensional data assimilation system are currently feasible. Recommendations for further reduction in computational cost are made with the goal of extending this technique to three-dimensional global assimilation systems.

Tangborn, Andrew↗

Vertical motion as inferred by the VAS 6.7 micron water vapor channel

The change of topography of the upper-most water vapor layer is used as a tracer to directly estimate vertical motion in subtropical atmospheres. An equation for vertical motion is derived using the relationship between 6.7-micron brightness temperature and the altitude of the upper layer of water vapor. This equation is used to estimate the vertical motion over the southeastern U.S. for July, 1984. The obtained values are compared to the vertical motion derived kinematically from radiosondes.

Stout, John↗

Isentropic mass exchange between the Tropics and extratropics in the stratosphere

The isentropic mass exchange between the Tropics and extratropics in the stratosphere is investigated with a semi-Lagrangian transport model for the periods from 1 June to 31 October 1992 and from 1 December 1992 to 30 April 1993 using winds from the U.K. Meteorological Office data assimilation system. Calculations with an idealized, initially zonally symmetric tracer show that in the middle and upper stratosphere the bulk of tropical air is transported into the midlatitudes of the winter hemisphere although there exist quasi-permeable barriers in the subtropics. The transport takes place in the form of planetary-scale 'tongues' of material that are drawn poleward in association with the episodic amplification of planetary-scale waves in high latitudes of the winter hemisphere. Once air of tropical origin is transported to the midlatitudes it is irreversibly mixed with the midlatitude air in the 'surf zone.' Air of tropical origin can, however, hardly penetrate into the interior of the winter polar vortex until the breakdown of the vortex. Transport of tropical air into the midlatitudes of the summer hemisphere is strongly inhibited. In the lower stratosphere, tropical air is transported into the northern and southern midlatitudes. During the period from 1 June to 31 October 1992, the amount of tropical air transported into the Northern Hemisphere is, however, much smaller than that transported into the Southern Hemisphere, and there exist strong gradients in the tracer field in the equatorial region, indicating that there is a quasi-permeable barrier to cross-equator mass exchange. During the period from 1 December 1992 to 30 April 1993, on the other hand, roughly the same amounts of tropical air are transported into the Southern Hemisphere, and there exist strong gradients in the tracer field in the equatorial region, indicating that there is a quasi-permeable barrier to cross-equator mass exchange. During the period from 1 December 1992 to 30 April 1993, on the other hand, roughly the same amounts of tropical air are transported into the northern midlatitudes and into the southern midlatitudes, and there exist no significant transport barriers in the equatorial area.

Chen, Ping↗

Approximation methods for the solution of inverse problems in lake and sea sediment core analysis

A theoretical model employing one-dimensional (depth) transport equations to describe vertical redistribution of ocean-floor and lake-floor sediment (particulates, volcanic ash, microtektites, or radioactive tracers) by episodic and nonepisodic events including bioturbation is developed analytically and demonstrated. The principles underlying the model are explained; the model equations are derived; the inverse problem of identifying the depth-dependent bioturbation coefficient is addressed; two approximation theorems are presented; and numerical results for two sample problems are presented graphically. It is suggested that compatification, porosity effects, and depth-dependent sedimentation be taken into account when formulating future models.

Banks, H. T.↗

Islet: interpolation semi-Lagrangian element-based transport

Abstract. Advection of trace species, or tracers, also called tracer transport, in models of the atmosphere and other physical domains is an important and potentially computationally expensive part of a model's dynamical core. Semi-Lagrangian (SL) advection methods are efficient because they permit a time step much larger than the advective stability limit for explicit Eulerian methods without requiring the solution of a globally coupled system of equations as implicit Eulerian methods do. Thus, to reduce the computational expense of tracer transport, dynamical cores often use SL methods to advect tracers. The class of interpolation semi-Lagrangian (ISL) methods contains potentially extremely efficient SL methods. We describe a finite-element ISL transport method that we call the interpolation semi-Lagrangian element-based transport (Islet) method, such as for use with atmosphere models discretized using the spectral element method. The Islet method uses three grids that share an element grid: a dynamics grid supporting, for example, the Gauss–Legendre–Lobatto basis of degree three; a physics parameterizations grid with a configurable number of finite-volume subcells per element; and a tracer grid supporting use of Islet bases with particular basis again configurable. This method provides extremely accurate tracer transport and excellent diagnostic values in a number of verification problems.

58 GEOSCIENCES↗

Parameterization of linear wave chemical transport in planetary atmospheres by eddy diffusion

Two-dimensional transport by linear waves is parameterized by wave velocity and temperature correlations and nonconservative properties of the atmosphere in a residual Eulerian framework. Explicit expressions derived for the eddy tracer flux vector are written in eddy diffusion tensor form, but it is demonstrated that strong coupling of chemistry and temperature through the rate constants, as is the case for ozone, requires additional terms proportional to the mixing ratio and not to its gradient. A variety of limiting cases are discussed. In the residual Eulerian framework the K(yz), K(zy) components are nonzero for a nonconservative atmosphere and/or a reactive tracer. Tracers have common transport coefficients in the limit where their chemistry is negligible in comparison to transience. It is shown that the usual one-dimensional vertical eddy diffusion continuity equation is rigorously valid only when the temperature dependence of the rate constants is negligible and the tracer has strong vertical stratification relative to its horizontal variability.

Strobel, D. F.↗

Modeling Insights into Deuterium Excess as an Indicator of Water Vapor Source Conditions

Deuterium excess (d) is interpreted in conventional paleoclimate reconstructions as a tracer of oceanic source region conditions, such as temperature, where precipitation originates. Previous studies have adopted co-isotopic approaches to estimate past changes in both site and oceanic source temperatures for ice core sites using empirical relationships derived from conceptual distillation models, particularly Mixed Cloud Isotopic Models (MCIMs). However, the relationship between d and oceanic surface conditions remains unclear in past contexts. We investigate this climate-isotope relationship for sites in Greenland and Antarctica using multiple simulations of the water isotope-enabled Goddard Institute for Space Studies (GISS) ModelE-R general circulation model and apply a novel suite of model vapor source distribution (VSD) tracers to assess d as a proxy for source temperature variability under a range of climatic conditions. Simulated average source temperatures determined by the VSDs are compared to synthetic source temperature estimates calculated using MCIM equations linking d to source region conditions. We show that although deuterium excess is generally a faithful tracer of source temperatures as estimated by the MCIM approach, large discrepancies in the isotope-climate relationship occur around Greenland during the Last Glacial Maximum simulation, when precipitation seasonality and moisture source regions were notably different from present. This identified sensitivity in d as a source temperature proxy suggests that quantitative climate reconstructions from deuterium excess should be treated with caution for some sites when boundary conditions are significantly different from the present day. Also, the exclusion of the influence of humidity and other evaporative source changes in MCIM regressions may be a limitation of quantifying source temperature fluctuations from deuterium excess in some instances.

simulation↗

In defense of Ertel's potential vorticity and its general applicability as a meteorological tracer

A potential vorticity theorem and its two summary statements published by Haynes and McIntyre are challenged conceptually by equations, discussions and examples. The apparent simplification proposed by the authors to convert from a mass to volume integral, i.e., by cancelling density against the specific volume in the potential vorticity, changes the physical significance of the integrand. It no longer is the potential vorticity. The resulting mean for either a bulk Eulerian or Lagrangian system is then not analogous to a mixing ratio and therefore not independent of the broad spectrum of internal waves, the independence that makes Ertel's potential vorticity so valuable either as a stratospheric tracer or as a predictive or diagnostic, large scale, meteorological variable.

Danielsen, Edwin F.↗

Cosmological constraints from the DESI DR1 joint power spectrum and bispectrum analysis

We derive cosmological parameter constraints from the Dark Energy Spectroscopic Instrument (DESI) Data Release 1 (DR1) galaxy clustering data, based on a joint full-shape analysis of the power spectrum multipoles and the bispectrum monopole using the ShapeFit framework. This is the follow-up of our previous work, in which we obtained for the first time constraints on the ShapeFit parameters using the bispectrum of DESI DR1. Here we present the first ShapeFit cosmological inference results using the bispectrum of DESI DR1. We recover values for the matter density parameter and Hubble constant of respectively $Ω_m=0.310\pm0.012$ and $H_0=[68.92\pm0.97]\,\mathrm{km\, s^{-1} Mpc^{-1}}$, consistent with previous results from the full DESI DR1 dataset that did not use the bispectrum signal. The inclusion of the bispectrum significantly tightens the constraints on the amplitude of fluctuations, reducing the error-bars in $\ln(A_s\times10^{10})$ by approximately 20%, compared to using the power spectrum alone. We also explore extended cosmological models by performing fits for the evolving dark energy equation of state $w_0w_a$, and the sum of neutrino masses $\sum m_ν$. In these cases, we obtain constraints slightly larger than the ones from previous works from the DESI collaboration, due to not combining the full-shape results with other probes in all tracers. We find no strong evidence of deviations from standard $Λ$CDM, with the dark energy equation-of-state remaining within 2$σ$ from a cosmological constant $Λ$, and the neutrino mass being consistent with the normal hierarchy, $\sum m_ν<0.1\,[eV]$ at 95% confidence limit. These constraints are broadly consistent with other DESI DR1 analyses, thus validating the robustness of the ShapeFit compression approach and the inclusion of the bispectrum for cosmological inference.

Novell-Masot, S. [ICC, Barcelona U.; Geneva U., De↗

Cosmological constraints from the DESI DR1 joint power spectrum and bispectrum analysis

We derive cosmological parameter constraints from the Dark Energy Spectroscopic Instrument (DESI) Data Release 1 (DR1) galaxy clustering data, based on a joint full-shape analysis of the power spectrum multipoles and the bispectrum monopole using the ShapeFit framework. This is the follow-up of our previous work, in which we obtained for the first time constraints on the ShapeFit parameters using the bispectrum of DESI DR1. Here we present the first ShapeFit cosmological inference results using the bispectrum of DESI DR1. We recover values for the matter density parameter and Hubble constant of respectively $Ω_m=0.310\pm0.012$ and $H_0=[68.92\pm0.97]\,\mathrm{km\, s^{-1} Mpc^{-1}}$, consistent with previous results from the full DESI DR1 dataset that did not use the bispectrum signal. The inclusion of the bispectrum significantly tightens the constraints on the amplitude of fluctuations, reducing the error-bars in $\ln(A_s\times10^{10})$ by approximately 20%, compared to using the power spectrum alone. We also explore extended cosmological models by performing fits for the evolving dark energy equation of state $w_0w_a$, and the sum of neutrino masses $\sum m_ν$. In these cases, we obtain constraints slightly larger than the ones from previous works from the DESI collaboration, due to not combining the full-shape results with other probes in all tracers. We find no strong evidence of deviations from standard $Λ$CDM, with the dark energy equation-of-state remaining within 2$σ$ from a cosmological constant $Λ$, and the neutrino mass being consistent with the normal hierarchy, $\sum m_ν<0.1\,[eV]$ at 95% confidence limit. These constraints are broadly consistent with other DESI DR1 analyses, thus validating the robustness of the ShapeFit compression approach and the inclusion of the bispectrum for cosmological inference.

Novell-Masot, S. [ICC, Barcelona U.; Geneva U., De↗

EMPDF : inferring the Milky Way mass with data-driven distribution function in phase space

We introduce the emPDF (empirical distribution function), a novel dynamical modelling method that infers the gravitational potential from kinematic tracers with optimal statistical efficiency under the minimal assumption of steady state. emPDF determines the best-fitting potential by maximizing the similarity between instantaneous kinematics and the time-averaged phase-space distribution function (DF), which is empirically constructed from observation upon the theoretical foundation of oPDF (Han et al. 2016). This approach eliminates the need for presumed functional forms of DFs or orbit libraries required by conventional DF- or orbit-based methods. emPDF stands out for its flexibility, efficiency, and capability in handling observational effects, making it preferable to the popular Jeans equation or other minimal assumption methods, especially for the Milky Way (MW) outer halo where tracers often have limited sample size and poor data quality. We apply emPDF to infer the MW mass profile using Gaia DR3 data of satellite galaxies and globular clusters, obtaining enclosed masses of M (,r) = 26±8, 46±8, 90±13⁠, and 149±40 x 10 10 M ⊙ at r = 30, 50, 100⁠, and 200 kpc, respectively. These are consistent with the updated constraints from simulation-informed DF fitting (Li et al. 2020). While the simulation-informed DF offers superior precision owing to the additional information extracted from simulations, emPDF is independent of such supplementary knowledge and applicable to general tracer populations. emPDF is currently implemented for tracers with complete 6D kinematics within spherical potentials, but it can potentially be extended to address more general problems.

Astrophysics of Galaxies (astro-ph.GA)↗

Investigation of the Influence of Microgravity on Transport Mechanism in a Virtual Spaceflight Chamber: A Flight Definition Program

A need exists for understanding precisely how particles move and interact in a fluid in the absence of gravity. Such understanding is required, for example, for modeling and predicting crystal growth in space where crystals grow from solution around nucleation sites as well as for any study of particles or bubbles in liquids or in experiments where particles are used as tracers for mapping microconvection. We have produced an exact solution to the general equation of motion of particles at extremely low Reynolds number in microgravity that covers a wide range of interesting conditions. We have also developed diagnostic tools and experimental techniques to test the validity of the general equation . This program, which started in May, 1998, will produce the flight definition for an experiment in a microgravity environment of space to validate the theoretical model. We will design an experiment with the help of the theoretical model that is optimized for testing the model, measuring g, g-jitter, and other microgravity phenomena. This paper describes the goals, rational, and approach for the flight definition program. The first objective of this research is to understand the physics of particle interactions with fluids and other particles in low Reynolds number flows in microgravity. Secondary objectives are to (1) observe and quantify g-jitter effects and microconvection on particles in fluids, (2) validate an exact solution to the general equation of motion of a particle in a fluid, and (3) to characterize the ability of isolation tables to isolate experiments containing particle in liquids. The objectives will be achieved by recording a large number of holograms of particle fields in microgravity under controlled conditions, extracting the precise three-dimensional position of all of the particles as a function of time and examining the effects of all parameters on the motion of the particles. The feasibility for achieving these results has already been established in the ongoing ground-based NRA, which led to the "virtual spaceflight chamber" concept.

Trolinger, James D.↗