Atmospheric chemistry and transport modeling in the outer solar system

NASA Astrophysics Data System (ADS)

Lee, Yuan-Tai (Anthony)

2001-11-01

This thesis consists of 1-D and 2-D photochemical- dynamical modeling in the upper atmospheres of outer planets. For 1-D modeling, a unified hydrocarbon photochemical model has been studied in Jupiter, Saturn, Uranus, Neptune, and Titan, by comparing with the Voyager observations, and the recent measurements of methyl radicals by ISO in Saturn and Neptune. The CH3 observation implies a kinetically sensitive test to the measured and estimated hydrocarbon rate constants at low temperatures. We identify the key reactions that control the concentrations of CH3 in the model, such as the three-body recombination reaction, CH3 + CH3 + M --> C 2H6 + M, and the recycling reaction H + CH3 + M --> CH4 + M. The results show reasonable agreement with ISO values. In Chapter 4, the detection of PH3 in the lower stratosphere and upper troposphere of Jupiter has provided a photochemical- dynamical coupling model to derive the eddy diffusion coefficient in the upper troposphere of Jupiter. Using a two-layers photochemical model with updated photodissociation cross-sections and chemical rate constants for NH3 and PH 3, we find that the upper tropospheric eddy diffusion coefficient <10 5 cm2 sec-1, and the deeper tropospheric value >106 cm2 sec-1, are required to match the derived PH3 vertical profile by the observation. The best-fit functional form derivation of eddy diffusion coefficient in the upper troposphere of Jupiter above 400 mbar is K = 2.0 × 104 (n/2.2 × 1019)-0.5 cm 2 sec-1. On the other hand, Chapter 5 demonstrates a dynamical-only 2-D model of C2H6 providing a complete test for the current 2-D transport models in Jovian lower stratosphere and upper troposphere (270 to 0.1 mbar pressure levels). Different combinations of residual advection, horizontal eddy dispersion, and vertical eddy mixing are examined at different latitudes.

Energy balance constraints on gravity wave induced eddy diffusion in the mesosphere and lower thermosphere

NASA Technical Reports Server (NTRS)

Strobel, D. F.; Apruzese, J. P.; Schoeberl, M. R.

1985-01-01

The constraints on turbulence improved by the mesospheric heat budget are reexamined, and the sufficiency of the theoretical evidence to support the hypothesis that the eddy Prandtl number is greater than one in the mesosphere is considered. The mesopause thermal structure is calculated with turbulent diffusion coefficients commonly used in chemical models and deduced from mean zonal wind deceleration. It is shown that extreme mesopause temperatures of less than 100 K are produced by the large net cooling. The results demonstrate the importance of the Prandtl number for mesospheric turbulence.

VERTICAL DIFFUSION IN SMALL STRATIFIED LAKES: DATA AND ERROR ANALYSIS

EPA Science Inventory

Water temperature profiles were measured at 2-min intervals in a stratified temperate lake with a surface area of 0.06 km2 and a aximum depth of 10 m from May 7 to August 9, 1989. he data were used to calculate the vertical eddy diffusion coefficient K2 in the hypolimnion. he dep...

A Semi-Analytical Model for Dispersion Modelling Studies in the Atmospheric Boundary Layer

NASA Astrophysics Data System (ADS)

Gupta, A.; Sharan, M.

2017-12-01

The severe impact of harmful air pollutants has always been a cause of concern for a wide variety of air quality analysis. The analytical models based on the solution of the advection-diffusion equation have been the first and remain the convenient way for modeling air pollutant dispersion as it is easy to handle the dispersion parameters and related physics in it. A mathematical model describing the crosswind integrated concentration is presented. The analytical solution to the resulting advection-diffusion equation is limited to a constant and simple profiles of eddy diffusivity and wind speed. In practice, the wind speed depends on the vertical height above the ground and eddy diffusivity profiles on the downwind distance from the source as well as the vertical height. In the present model, a method of eigen-function expansion is used to solve the resulting partial differential equation with the appropriate boundary conditions. This leads to a system of first order ordinary differential equations with a coefficient matrix depending on the downwind distance. The solution of this system, in general, can be expressed in terms of Peano-baker series which is not easy to compute, particularly when the coefficient matrix becomes non-commutative (Martin et al., 1967). An approach based on Taylor's series expansion is introduced to find the numerical solution of first order system. The method is applied to various profiles of wind speed and eddy diffusivities. The solution computed from the proposed methodology is found to be efficient and accurate in comparison to those available in the literature. The performance of the model is evaluated with the diffusion datasets from Copenhagen (Gryning et al., 1987) and Hanford (Doran et al., 1985). In addition, the proposed method is used to deduce three dimensional concentrations by considering the Gaussian distribution in crosswind direction, which is also evaluated with diffusion data corresponding to a continuous point source.

Evaluation of diffusion models in breast cancer.

PubMed

Panek, Rafal; Borri, Marco; Orton, Matthew; O'Flynn, Elizabeth; Morgan, Veronica; Giles, Sharon L; deSouza, Nandita; Leach, Martin O; Schmidt, Maria A

2015-08-01

The purpose of this study is to investigate whether the microvascular pseudodiffusion effects resulting with non-monoexponential behavior are present in breast cancer, taking into account tumor spatial heterogeneity. Additionally, methodological factors affecting the signal in low and high diffusion-sensitizing gradient ranges were explored in phantom studies. The effect of eddy currents and accuracy of b-value determination using a multiple b-value diffusion-weighted MR imaging sequence were investigated in test objects. Diffusion model selection and noise were then investigated in volunteers (n = 5) and breast tumor patients (n = 21) using the Bayesian information criterion. 54.3% of lesion voxels were best fitted by a monoexponential, 26.2% by a stretched-exponential, and 19.5% by a biexponential intravoxel incoherent motion (IVIM) model. High correlation (0.92) was observed between diffusion coefficients calculated using mono- and stretched-exponential models and moderate (0.59) between monoexponential and IVIM (medians: 0.96/0.84/0.72 × 10(-3) mm(2)/s, respectively). Distortion due to eddy currents depended on the direction of the diffusion gradient and displacement varied between 1 and 6 mm for high b-value images. Shift in the apparent diffusion coefficient due to intrinsic field gradients was compensated for by averaging diffusion data obtained from opposite directions. Pseudodiffusion and intravoxel heterogeneity effects were not observed in approximately half of breast cancer and normal tissue voxels. This result indicates that stretched and IVIM models should be utilized in regional analysis rather than global tumor assessment. Cross terms between diffusion-sensitization gradients and other imaging or susceptibility-related gradients are relevant in clinical protocols, supporting the use of geometric averaging of diffusion-weighted images acquired with diffusion-sensitization gradients in opposite directions.

A study of atmospheric diffusion from the LANDSAT imagery. [pollution transport over the ocean

NASA Technical Reports Server (NTRS)

Dejesusparada, N. (Principal Investigator); Viswanadham, Y.; Torsani, J. A.

1981-01-01

LANDSAT multispectral scanner data of the smoke plumes which originated in eastern Cabo Frio, Brazil and crossed over into the Atlantic Ocean, are analyzed to illustrate how high resolution LANDSAT imagery can aid meteorologists in evaluating specific air pollution events. The eleven LANDSAT images selected are for different months and years. The results show that diffusion is governed primarily by water and air temperature differences. With colder water, low level air is very stable and the vertical diffusion is minimal; but water warmer than the air induces vigorous diffusion. The applicability of three empirical methods for determining the horizontal eddy diffusivity coefficient in the Gaussian plume formula was evaluated with the estimated standard deviation of the crosswind distribution of material in the plume from the LANDSAT imagery. The vertical diffusion coefficient in stable conditions is estimated using Weinstock's formulation. These results form a data base for use in the development and validation of meso scale atmospheric diffusion models.

Steepest Ascent Low/Non-Low-Frequency Ratio in Empirical Mode Decomposition to Separate Deterministic and Stochastic Velocities From a Single Lagrangian Drifter

NASA Astrophysics Data System (ADS)

Chu, Peter C.

2018-03-01

SOund Fixing And Ranging (RAFOS) floats deployed by the Naval Postgraduate School (NPS) in the California Current system from 1992 to 2001 at depth between 150 and 600 m (http://www.oc.nps.edu/npsRAFOS/) are used to study 2-D turbulent characteristics. Each drifter trajectory is adaptively decomposed using the empirical mode decomposition (EMD) into a series of intrinsic mode functions (IMFs) with corresponding specific scale for each IMF. A new steepest ascent low/non-low-frequency ratio is proposed in this paper to separate a Lagrangian trajectory into low-frequency (nondiffusive, i.e., deterministic) and high-frequency (diffusive, i.e., stochastic) components. The 2-D turbulent (or called eddy) diffusion coefficients are calculated on the base of the classical turbulent diffusion with mixing length theory from stochastic component of a single drifter. Statistical characteristics of the calculated 2-D turbulence length scale, strength, and diffusion coefficients from the NPS RAFOS data are presented with the mean values (over the whole drifters) of the 2-D diffusion coefficients comparable to the commonly used diffusivity tensor method.

The altitude distribution of the Venus ultraviolet nightglow and implications on vertical transport

NASA Technical Reports Server (NTRS)

Gerard, J. C.; Stewart, A. I. F.; Bougher, S. W.

1981-01-01

The altitude distribution of the nitric oxide nightglow was measured with an ultraviolet spectrometer on board Pioneer Venus, in order to study the effects of the distribution on the Venus nightside lower thermosphere transport properties. Limb profiles were obtained with an 8 ms integration period on several orbits near periapsis. The observations were made between P minus 2 min and P plus 4 min, where altitude ranges between 150 and 350 km, and latitude varies from 24 degrees N to 9 degrees S. A method independent of the spacecraft attitude data was used to fit the observed limb profiles, and to find the altitude of the maximum of the layer (115 plus or minus 2 km), and the topside scale height (about 3 km). It is shown that downward transport by diffusion alone is not sufficient, and if vertical motion is parameterized by eddy diffusion, an eddy diffusion coefficient is deduced from the altitude of the layer.

A model for calculating the vertical distribution of the atmospheric electric potential in the exchange layer in a maritime clean atmosphere

NASA Astrophysics Data System (ADS)

Kulkarni, M. N.; Kamra, A. K.

2012-11-01

A theoretical model is developed for calculating the vertical distribution of atmospheric electric potential in exchange layer of maritime clean atmosphere. The transport of space charge in electrode layer acts as a convective generator in this model and plays a major role in determining potential distribution in vertical. Eddy diffusion is the main mechanism responsible for the distribution of space charge in vertical. Our results show that potential at a particular level increases with increase in the strength of eddy diffusion under similar conditions. A method is suggested to estimate columnar resistance, the ionospheric potential and the vertical atmospheric electric potential distribution in exchange layer from measurements of total air-earth current density and surface electric field made over oceans. The results are validated and found to be in very good agreement with the previous aircraft measurements. Different parameters involved in the proposed methodology can be determined either theoretically, as in the present work, or experimentally using the near surface atmospheric electrical measurements or using some other surface-based measurement technique such as LIDAR. A graphical relationship between the atmospheric eddy diffusion coefficient and height of exchange layer obtained from atmospheric electrical approach, is reported.

Observational Inferences of Lateral Eddy Diffusivity in the Halocline of the Beaufort Gyre

NASA Astrophysics Data System (ADS)

Meneghello, Gianluca; Marshall, John; Cole, Sylvia T.; Timmermans, Mary-Louise

2017-12-01

Using Ekman pumping rates mediated by sea ice in the Arctic Ocean's Beaufort Gyre (BG), the magnitude of lateral eddy diffusivities required to balance downward pumping is inferred. In this limit—that of vanishing residual-mean circulation—eddy-induced upwelling exactly balances downward pumping. The implied eddy diffusivity varies spatially and decays with depth, with values of 50-400 m2/s. Eddy diffusivity estimated using mixing length theory applied to BG mooring data exhibits a similar decay with depth and range of values from 100 m2/s to more than 600 m2/s. We conclude that eddy diffusivities in the BG are likely large enough to balance downward Ekman pumping, arresting the deepening of the gyre and suggesting that eddies play a zero-order role in buoyancy and freshwater budgets of the BG.

Anisotropic Mesoscale Eddy Transport in Ocean General Circulation Models

NASA Astrophysics Data System (ADS)

Reckinger, S. J.; Fox-Kemper, B.; Bachman, S.; Bryan, F.; Dennis, J.; Danabasoglu, G.

2014-12-01

Modern climate models are limited to coarse-resolution representations of large-scale ocean circulation that rely on parameterizations for mesoscale eddies. The effects of eddies are typically introduced by relating subgrid eddy fluxes to the resolved gradients of buoyancy or other tracers, where the proportionality is, in general, governed by an eddy transport tensor. The symmetric part of the tensor, which represents the diffusive effects of mesoscale eddies, is universally treated isotropically in general circulation models. Thus, only a single parameter, namely the eddy diffusivity, is used at each spatial and temporal location to impart the influence of mesoscale eddies on the resolved flow. However, the diffusive processes that the parameterization approximates, such as shear dispersion, potential vorticity barriers, oceanic turbulence, and instabilities, typically have strongly anisotropic characteristics. Generalizing the eddy diffusivity tensor for anisotropy extends the number of parameters to three: a major diffusivity, a minor diffusivity, and the principal axis of alignment. The Community Earth System Model (CESM) with the anisotropic eddy parameterization is used to test various choices for the newly introduced parameters, which are motivated by observations and the eddy transport tensor diagnosed from high resolution simulations. Simply setting the ratio of major to minor diffusivities to a value of five globally, while aligning the major axis along the flow direction, improves biogeochemical tracer ventilation and reduces global temperature and salinity biases. These effects can be improved even further by parameterizing the anisotropic transport mechanisms in the ocean.

Anisotropic mesoscale eddy transport in ocean general circulation models

NASA Astrophysics Data System (ADS)

Reckinger, Scott; Fox-Kemper, Baylor; Bachman, Scott; Bryan, Frank; Dennis, John; Danabasoglu, Gokhan

2014-11-01

In modern climate models, the effects of oceanic mesoscale eddies are introduced by relating subgrid eddy fluxes to the resolved gradients of buoyancy or other tracers, where the proportionality is, in general, governed by an eddy transport tensor. The symmetric part of the tensor, which represents the diffusive effects of mesoscale eddies, is universally treated isotropically. However, the diffusive processes that the parameterization approximates, such as shear dispersion and potential vorticity barriers, typically have strongly anisotropic characteristics. Generalizing the eddy diffusivity tensor for anisotropy extends the number of parameters from one to three: major diffusivity, minor diffusivity, and alignment. The Community Earth System Model (CESM) with the anisotropic eddy parameterization is used to test various choices for the parameters, which are motivated by observations and the eddy transport tensor diagnosed from high resolution simulations. Simply setting the ratio of major to minor diffusivities to a value of five globally, while aligning the major axis along the flow direction, improves biogeochemical tracer ventilation and reduces temperature and salinity biases. These effects can be improved by parameterizing the oceanic anisotropic transport mechanisms.

Isopycnal mixing by mesoscale eddies significantly impacts oceanic anthropogenic carbon uptake

NASA Astrophysics Data System (ADS)

Gnanadesikan, Anand; Pradal, Marie-Aude; Abernathey, Ryan

2015-06-01

Anthropogenic carbon dioxide uptake varies across Earth System Models for reasons that have remained obscure. When varied within a single model, the lateral eddy mixing coefficient ARedi produces a range of uptake similar to the modeled range. The highest uptake, resulting from a simulation with a constant ARedi of 2400 m2/s, simulates 15% more historical carbon uptake than a model with ARedi = 400 m2/s. A sudden doubling in carbon dioxide produces a 21% range in carbon uptake across the models. Two spatially dependent representations of ARedi produce uptake that lies in the middle of the range of constant values despite predicting very large values in the subtropical gyres. One-dimensional diffusive models of the type used for integrated assessments can be fit to the simulations, with ARedi accounting for a substantial fraction of the effective vertical diffusion. Such models, however, mask significant regional changes in stratification and biological carbon storage.

Observational Inferences of Lateral Eddy Diffusivity in the Halocline of the Beaufort Gyre

NASA Astrophysics Data System (ADS)

Meneghello, G.; Marshall, J.; Cole, S. T.; Timmermans, M. L.

2017-12-01

Using Ekman pumping rates mediated by sea-ice in the Arctic Ocean's Beaufort Gyre (BG), the magnitude of lateral eddy diffusivities required to balance downward pumping is inferred. In this limit — that of vanishing residual-mean circulation — eddy-induced upwelling exactly balances downward pumping. The implied eddy diffusivity varies spatially with values of 50-400 m2/s, and decays with depth. Eddy diffusivity estimated using mixing length theory applied to BG mooring data exhibits a similar range of values from 100 m2/s to more than 600 m2/s, and also decays with depth. We conclude that eddy diffusivities in the BG are likely large enough to balance downward Ekman pumping, arresting the deepening of the gyre and suggesting that eddies play a zero-order role in buoyancy and freshwater budgets of the BG.

Simulations of eddy kinetic energy transport in barotropic turbulence

NASA Astrophysics Data System (ADS)

Grooms, Ian

2017-11-01

Eddy energy transport in rotating two-dimensional turbulence is investigated using numerical simulation. Stochastic forcing is used to generate an inhomogeneous field of turbulence and the time-mean energy profile is diagnosed. An advective-diffusive model for the transport is fit to the simulation data by requiring the model to accurately predict the observed time-mean energy distribution. Isotropic harmonic diffusion of energy is found to be an accurate model in the case of uniform, solid-body background rotation (the f plane), with a diffusivity that scales reasonably well with a mixing-length law κ ∝V ℓ , where V and ℓ are characteristic eddy velocity and length scales. Passive tracer dynamics are added and it is found that the energy diffusivity is 75 % of the tracer diffusivity. The addition of a differential background rotation with constant vorticity gradient β leads to significant changes to the energy transport. The eddies generate and interact with a mean flow that advects the eddy energy. Mean advection plus anisotropic diffusion (with reduced diffusivity in the direction of the background vorticity gradient) is moderately accurate for flows with scale separation between the eddies and mean flow, but anisotropic diffusion becomes a much less accurate model of the transport when scale separation breaks down. Finally, it is observed that the time-mean eddy energy does not look like the actual eddy energy distribution at any instant of time. In the future, stochastic models of the eddy energy transport may prove more useful than models of the mean transport for predicting realistic eddy energy distributions.

Efficacy of distortion correction on diffusion imaging: comparison of FSL eddy and eddy_correct using 30 and 60 directions diffusion encoding.

PubMed

Yamada, Haruyasu; Abe, Osamu; Shizukuishi, Takashi; Kikuta, Junko; Shinozaki, Takahiro; Dezawa, Ko; Nagano, Akira; Matsuda, Masayuki; Haradome, Hiroki; Imamura, Yoshiki

2014-01-01

Diffusion imaging is a unique noninvasive tool to detect brain white matter trajectory and integrity in vivo. However, this technique suffers from spatial distortion and signal pileup or dropout originating from local susceptibility gradients and eddy currents. Although there are several methods to mitigate these problems, most techniques can be applicable either to susceptibility or eddy-current induced distortion alone with a few exceptions. The present study compared the correction efficiency of FSL tools, "eddy_correct" and the combination of "eddy" and "topup" in terms of diffusion-derived fractional anisotropy (FA). The brain diffusion images were acquired from 10 healthy subjects using 30 and 60 directions encoding schemes based on the electrostatic repulsive forces. For the 30 directions encoding, 2 sets of diffusion images were acquired with the same parameters, except for the phase-encode blips which had opposing polarities along the anteroposterior direction. For the 60 directions encoding, non-diffusion-weighted and diffusion-weighted images were obtained with forward phase-encoding blips and non-diffusion-weighted images with the same parameter, except for the phase-encode blips, which had opposing polarities. FA images without and with distortion correction were compared in a voxel-wise manner with tract-based spatial statistics. We showed that images corrected with eddy and topup possessed higher FA values than images uncorrected and corrected with eddy_correct with trilinear (FSL default setting) or spline interpolation in most white matter skeletons, using both encoding schemes. Furthermore, the 60 directions encoding scheme was superior as measured by increased FA values to the 30 directions encoding scheme, despite comparable acquisition time. This study supports the combination of eddy and topup as a superior correction tool in diffusion imaging rather than the eddy_correct tool, especially with trilinear interpolation, using 60 directions encoding scheme.

PSF mapping-based correction of eddy-current-induced distortions in diffusion-weighted echo-planar imaging.

PubMed

In, Myung-Ho; Posnansky, Oleg; Speck, Oliver

2016-05-01

To accurately correct diffusion-encoding direction-dependent eddy-current-induced geometric distortions in diffusion-weighted echo-planar imaging (DW-EPI) and to minimize the calibration time at 7 Tesla (T). A point spread function (PSF) mapping based eddy-current calibration method is newly presented to determine eddy-current-induced geometric distortions even including nonlinear eddy-current effects within the readout acquisition window. To evaluate the temporal stability of eddy-current maps, calibration was performed four times within 3 months. Furthermore, spatial variations of measured eddy-current maps versus their linear superposition were investigated to enable correction in DW-EPIs with arbitrary diffusion directions without direct calibration. For comparison, an image-based eddy-current correction method was additionally applied. Finally, this method was combined with a PSF-based susceptibility-induced distortion correction approach proposed previously to correct both susceptibility and eddy-current-induced distortions in DW-EPIs. Very fast eddy-current calibration in a three-dimensional volume is possible with the proposed method. The measured eddy-current maps are very stable over time and very similar maps can be obtained by linear superposition of principal-axes eddy-current maps. High resolution in vivo brain results demonstrate that the proposed method allows more efficient eddy-current correction than the image-based method. The combination of both PSF-based approaches allows distortion-free images, which permit reliable analysis in diffusion tensor imaging applications at 7T. © 2015 Wiley Periodicals, Inc.

3-D residual eddy current field characterisation: applied to diffusion weighted magnetic resonance imaging.

PubMed

O'Brien, Kieran; Daducci, Alessandro; Kickler, Nils; Lazeyras, Francois; Gruetter, Rolf; Feiweier, Thorsten; Krueger, Gunnar

2013-08-01

Clinical use of the Stejskal-Tanner diffusion weighted images is hampered by the geometric distortions that result from the large residual 3-D eddy current field induced. In this work, we aimed to predict, using linear response theory, the residual 3-D eddy current field required for geometric distortion correction based on phantom eddy current field measurements. The predicted 3-D eddy current field induced by the diffusion-weighting gradients was able to reduce the root mean square error of the residual eddy current field to ~1 Hz. The model's performance was tested on diffusion weighted images of four normal volunteers, following distortion correction, the quality of the Stejskal-Tanner diffusion-weighted images was found to have comparable quality to image registration based corrections (FSL) at low b-values. Unlike registration techniques the correction was not hindered by low SNR at high b-values, and results in improved image quality relative to FSL. Characterization of the 3-D eddy current field with linear response theory enables the prediction of the 3-D eddy current field required to correct eddy current induced geometric distortions for a wide range of clinical and high b-value protocols.

Residual fields from extinct dynamos

NASA Astrophysics Data System (ADS)

Parker, E. N.

The generation of magnetic fields in convective zones of declining vigor and/or thickness is considered, the goal being to explain the magnetic fields observed in A-stars. The investigation is restricted to kinematical dynamos in order to show some of the many possibilities, which depend on the assumed conditions of decline of the convection. The examples illustrate the quantitative detail required to describe the convection in order to extract any firm conclusions concerning specific stars. The first example treats the basic problem of diffusion from a layer of declining thickness. The second has a buoyant rise added to the field in the layer. The third deals with plane dynamo waves in a region with declining eddy diffusivity, dynamo coefficient, and large-scale shear. It is noted that the dynamo number may increase or decrease with declining convection, with an increase expected if the large-scale shear does not decline as rapidly as the eddy diffusivity. It is shown that one of the components of the field may increase without bound even when the dynamo number declines to zero.

The effects of intermittency on statistical characteristics of turbulence and scale similarity of breakdown coefficients

NASA Astrophysics Data System (ADS)

Novikov, E. A.

1990-05-01

The influence of intermittency on turbulent diffusion is expressed in terms of the statistics of the dissipation field. The high-order moments of relative diffusion are obtained by using the concept of scale similarity of the breakdown coefficients (bdc). The method of bdc is useful for obtaining new models and general results, which then can be expressed in terms of multifractals. In particular, the concavity and other properties of spectral codimension are proved. Special attention is paid to the logarithmically periodic modulations. The parametrization of small-scale intermittent turbulence, which can be used for large-eddy simulation, is presented. The effect of molecular viscosity is taken into account in the spirit of the renorm group, but without spectral series, ɛ expansion, and fictitious random forces.

Towards a more consistent picture of isopycnal mixing in climate models

NASA Astrophysics Data System (ADS)

Gnanadesikan, A.; Pradal, M. A. S.; Koszalka, I.; Abernathey, R. P.

2014-12-01

The stirring of tracers by mesoscale eddies along isopycnal surfaces is often represented in coarse-resolution models by the Redi diffusion parameter ARedi. Theoretical treatments of ARedi often assume it should scale as the eddy energy or the growth rate of mesoscale eddies,. producing a picture where it is high in boundary currents and low )of order a few hundred m2/s) in the gyre interiors. However, observational estimates suggest that ARedi should be very large (of order thousands of m2/s) in the gyre interior. We present results of recent simulations comparing a range of spatially constant values ARedi (with values of 400, 800, 1200 and 2400 m2/s) to a spatially resolved estimate based on altimetry and a zonally averaged version of the same estimate. In general, increasing the ARedi coefficient destratifies and warms the high latitudes. Relative to our control simulation, the spatially dependent coefficient is lower in the Southern Ocean, but high in the North Pacific, and so the temperature changes mirror this. We also examine the response of ocean hypoxia to these changes. In general, the zonally averaged version of the altimetry-based estimate of ARedi does not capture the full 2d representation.

Uncertainty for calculating transport on Titan: A probabilistic description of bimolecular diffusion parameters

NASA Astrophysics Data System (ADS)

Plessis, S.; McDougall, D.; Mandt, K.; Greathouse, T.; Luspay-Kuti, A.

2015-11-01

Bimolecular diffusion coefficients are important parameters used by atmospheric models to calculate altitude profiles of minor constituents in an atmosphere. Unfortunately, laboratory measurements of these coefficients were never conducted at temperature conditions relevant to the atmosphere of Titan. Here we conduct a detailed uncertainty analysis of the bimolecular diffusion coefficient parameters as applied to Titan's upper atmosphere to provide a better understanding of the impact of uncertainty for this parameter on models. Because temperature and pressure conditions are much lower than the laboratory conditions in which bimolecular diffusion parameters were measured, we apply a Bayesian framework, a problem-agnostic framework, to determine parameter estimates and associated uncertainties. We solve the Bayesian calibration problem using the open-source QUESO library which also performs a propagation of uncertainties in the calibrated parameters to temperature and pressure conditions observed in Titan's upper atmosphere. Our results show that, after propagating uncertainty through the Massman model, the uncertainty in molecular diffusion is highly correlated to temperature and we observe no noticeable correlation with pressure. We propagate the calibrated molecular diffusion estimate and associated uncertainty to obtain an estimate with uncertainty due to bimolecular diffusion for the methane molar fraction as a function of altitude. Results show that the uncertainty in methane abundance due to molecular diffusion is in general small compared to eddy diffusion and the chemical kinetics description. However, methane abundance is most sensitive to uncertainty in molecular diffusion above 1200 km where the errors are nontrivial and could have important implications for scientific research based on diffusion models in this altitude range.

Efficacy of Distortion Correction on Diffusion Imaging: Comparison of FSL Eddy and Eddy_Correct Using 30 and 60 Directions Diffusion Encoding

PubMed Central

Yamada, Haruyasu; Abe, Osamu; Shizukuishi, Takashi; Kikuta, Junko; Shinozaki, Takahiro; Dezawa, Ko; Nagano, Akira; Matsuda, Masayuki; Haradome, Hiroki; Imamura, Yoshiki

2014-01-01

Diffusion imaging is a unique noninvasive tool to detect brain white matter trajectory and integrity in vivo. However, this technique suffers from spatial distortion and signal pileup or dropout originating from local susceptibility gradients and eddy currents. Although there are several methods to mitigate these problems, most techniques can be applicable either to susceptibility or eddy-current induced distortion alone with a few exceptions. The present study compared the correction efficiency of FSL tools, “eddy_correct” and the combination of “eddy” and “topup” in terms of diffusion-derived fractional anisotropy (FA). The brain diffusion images were acquired from 10 healthy subjects using 30 and 60 directions encoding schemes based on the electrostatic repulsive forces. For the 30 directions encoding, 2 sets of diffusion images were acquired with the same parameters, except for the phase-encode blips which had opposing polarities along the anteroposterior direction. For the 60 directions encoding, non–diffusion-weighted and diffusion-weighted images were obtained with forward phase-encoding blips and non–diffusion-weighted images with the same parameter, except for the phase-encode blips, which had opposing polarities. FA images without and with distortion correction were compared in a voxel-wise manner with tract-based spatial statistics. We showed that images corrected with eddy and topup possessed higher FA values than images uncorrected and corrected with eddy_correct with trilinear (FSL default setting) or spline interpolation in most white matter skeletons, using both encoding schemes. Furthermore, the 60 directions encoding scheme was superior as measured by increased FA values to the 30 directions encoding scheme, despite comparable acquisition time. This study supports the combination of eddy and topup as a superior correction tool in diffusion imaging rather than the eddy_correct tool, especially with trilinear interpolation, using 60 directions encoding scheme. PMID:25405472

Scalar rate correlation at a turbulent liquid free surface - A two-regime correlation for high Schmidt numbers

NASA Technical Reports Server (NTRS)

Khoo, Boo-Cheong; Sonin, Ain A.

1992-01-01

An experimental correlation is derived for gas absorption at a turbulent, shear-free liquid interface. The correlation is expressed in terms of the liquid-side turbulence intensity, liquid-side macroscale, and the properties of the diffusing gas and solvent. The transfer coefficient increases linearly with rms velocity up to a point where the eddy Reynolds number reaches a critical (Schmidt number dependent) value. At higher velocities, there is a more rapid linear rise. The slope of the lower Reynolds number region is proportional to the square root of the diffusivity; at Reynolds numbers much higher than that of the break point, the slope becomes independent of diffusivity.

Double-spin-echo diffusion weighting with a modified eddy current adjustment.

PubMed

Finsterbusch, Jürgen

2010-04-01

Magnetic field inhomogeneities like eddy current-related gradient fields cause geometric distortions in echo-planar imaging (EPI). This in particular affects diffusion-weighted imaging where these distortions vary with the direction of the diffusion weighting and hamper the accurate determination of diffusion parameters. The double-spin-echo preparation often used aims to reduce the cumulative eddy current effect by adjusting the diffusion-weighting gradient pulse durations to the time constant of the dominant eddy current contribution. However, eddy currents with a variety of time constants may be present and cause residual distortions. Here, a modification is proposed where the two bipolar gradient pairs of the preparation are adjusted independently to different time constants. At the expense of a slightly prolonged echo time, residual geometric distortions and correspondingly increased values of the diffusion anisotropy can be reduced as is demonstrated in phantoms and the human brain. Thus, it may help to improve the reliability of diffusion-weighted EPI. Copyright 2010 Elsevier Inc. All rights reserved.

Vertical eddy diffusivity as a control parameter in the tropical Pacific

NASA Astrophysics Data System (ADS)

Martinez Avellaneda, N.; Cornuelle, B.

2011-12-01

Ocean models suffer from errors in the treatment of turbulent sub-grid-scale motions responsible for mixing and energy dissipation. Unrealistic small-scale physics in models can have large-scale consequences, such as biases in the upper ocean temperature, a symptom of poorly-simulated upwelling, currents and air-sea interactions. This is of special importance in the tropical Pacific Ocean (TP), which is home to energetic air-sea interactions that affect global climate. It has been shown in a number of studies that the simulated ENSO variability is highly dependent on the state of the ocean (e.g.: background mixing). Moreover, the magnitude of the vertical numerical diffusion is of primary importance in properly reproducing the Pacific equatorial thermocline. This work is part of a NASA-funded project to estimate the space- and time-varying ocean mixing coefficients in an eddy-permitting (1/3dgr) model of the TP to obtain an improved estimate of its time-varying circulation and its underlying dynamics. While an estimation procedure for the TP (26dgr S - 30dgr N) in underway using the MIT general circulation model, complementary adjoint-based sensitivity studies have been carried out for the starting ocean state from Forget (2010). This analysis aids the interpretation of the estimated mixing coefficients and possible error compensation. The focus of the sensitivity tests is the Equatorial Undercurrent and sub-thermocline jets (i.e., Tsuchiya Jets), which have been thought to have strong dependence on vertical diffusivity and should provide checks on the estimated mixing parameters. In order to build intuition for the vertical diffusivity adjoint results in the TP, adjoint and forward perturbed simulations were carried out for an idealized sharp thermocline in a rectangular domain.

Evaluation of a scalar eddy transport coefficient based on geometric constraints

NASA Astrophysics Data System (ADS)

Bachman, S. D.; Marshall, D. P.; Maddison, J. R.; Mak, J.

2017-01-01

A suite of idealized models is used to evaluate and compare several previously proposed scalings for the eddy transport coefficient in downgradient mesoscale eddy closures. Of special interest in this comparison is a scaling introduced as part of the eddy parameterization framework of Marshall et al. (2012), which is derived using the inherent geometry of the Eliassen-Palm eddy flux tensor. The primary advantage of using this coefficient in a downgradient closure is that all dimensional terms are explicitly specified and the only uncertainty is a nondimensional parameter, α, which is bounded by one in magnitude. In each model a set of passive tracers is initialized, whose flux statistics are used to invert for the eddy-induced tracer transport. Unlike previous work, where this technique has been employed to diagnose the tensor coefficient of a linear flux-gradient relationship, the idealization of these models allows the lateral eddy transport to be described by a scalar coefficient. The skill of the extant scalings is then measured by comparing their predicted values against the coefficients diagnosed using this method. The Marshall et al. (2012), scaling is shown to scale most closely with the diagnosed coefficients across all simulations. It is shown that the skill of this scaling is due to its functional dependence on the total eddy energy, and that this scaling provides an excellent match to the diagnosed fluxes even in the limit of constant α. Possible extensions to this work, including how to incorporate the resultant transport coefficient into the Gent and McWilliams parameterization, are discussed.

Correction of eddy current distortions in high angular resolution diffusion imaging.

PubMed

Zhuang, Jiancheng; Lu, Zhong-Lin; Vidal, Christine Bouteiller; Damasio, Hanna

2013-06-01

To correct distortions caused by eddy currents induced by large diffusion gradients during high angular resolution diffusion imaging without any auxiliary reference scans. Image distortion parameters were obtained by image coregistration, performed only between diffusion-weighted images with close diffusion gradient orientations. A linear model that describes distortion parameters (translation, scale, and shear) as a function of diffusion gradient directions was numerically computed to allow individualized distortion correction for every diffusion-weighted image. The assumptions of the algorithm were successfully verified in a series of experiments on phantom and human scans. Application of the proposed algorithm in high angular resolution diffusion images markedly reduced eddy current distortions when compared to results obtained with previously published methods. The method can correct eddy current artifacts in the high angular resolution diffusion images, and it avoids the problematic procedure of cross-correlating images with significantly different contrasts resulting from very different gradient orientations or strengths. Copyright © 2012 Wiley Periodicals, Inc.

Results of rocket measurements of D-region ionization over Thumba in MAP

NASA Technical Reports Server (NTRS)

Chakrabarty, D. K.; Beig, G.; Garg, S. C.; Subrahmanyam, P.; Zalpuri, K. S.; Somayajulu, Y. V.; Rao, M. N. M.; Tandel, C. B.; Murlikrishna, T. R.

1989-01-01

Under MAP, two rockets were launched from Thumba (8.5 N, 76.8 E) around 1030 hrs Lt with identical payloads on 7 and 10 March 1986 for D region studies. Positive ion densities were measured by spherical probe and Gerdien condenser and electron densities were measured by Langmuir probe and propagation experiments. In both flights a valley in ionization height profile was noticed around 83 km. The density of ionization at this altitude was about 4 x 10(2) cu cm. A detailed positive ion-chemical scheme was used to reproduce the measured ionization height profiles. The density of NO needed to reproduce the valley in ionization at 83 km came around 5 x 10(5) cu cm. A photochemical treatment without diffusion process was found inadequate to explain this value of NO. Calculations showed that the value of vertical eddy diffusion needed to reproduce the value of NO was around 10(6)sq cm/s. Interestingly, the same value of eddy diffusion coefficient was obtained when derived in the manner described by Thrane and his coworkers using only the positive ion current data of spherical probes.

Lateral eddy diffusivity estimates from simulated and observed drifter trajectories: a case study for the Agulhas Current system

NASA Astrophysics Data System (ADS)

Rühs, Siren; Zhurbas, Victor; Durgadoo, Jonathan V.; Biastoch, Arne

2017-04-01

The Lagrangian description of fluid motion by sets of individual particle trajectories is extensively used to characterize connectivity between distinct oceanic locations. One important factor influencing the connectivity is the average rate of particle dispersal, generally quantified as Lagrangian diffusivity. In addition to Lagrangian observing programs, Lagrangian analyses are performed by advecting particles with the simulated flow field of ocean general circulation models (OGCMs). However, depending on the spatio-temporal model resolution, not all scale-dependent processes are explicitly resolved in the simulated velocity fields. Consequently, the dispersal of advective Lagrangian trajectories has been assumed not to be sufficiently diffusive compared to observed particle spreading. In this study we present a detailed analysis of the spatially variable lateral eddy diffusivity characteristics of advective drifter trajectories simulated with realistically forced OGCMs and compare them with estimates based on observed drifter trajectories. The extended Agulhas Current system around South Africa, known for its intricate mesoscale dynamics, serves as a test case. We show that a state-of-the-art eddy-resolving OGCM indeed features theoretically derived dispersion characteristics for diffusive regimes and realistically represents Lagrangian eddy diffusivity characteristics obtained from observed surface drifter trajectories. The estimates for the maximum and asymptotic lateral single-particle eddy diffusivities obtained from the observed and simulated drifter trajectories show a good agreement in their spatial pattern and magnitude. We further assess the sensitivity of the simulated lateral eddy diffusivity estimates to the temporal and lateral OGCM output resolution and examine the impact of the different eddy diffusivity characteristics on the Lagrangian connectivity between the Indian Ocean and the South Atlantic.

Observed and Simulated Eddy Diffusivity Upstream of the Drake Passage

NASA Astrophysics Data System (ADS)

Tulloch, R.; Ferrari, R. M.; Marshall, J.

2012-12-01

Estimates of eddy diffusivity in the Southern Ocean are poorly constrained due to lack of observations. We compare the first direct estimate of isopycnal eddy diffusivity upstream of the Drake Passage (from Ledwell et al. 2011) with a numerical simulation. The estimate is computed from a point tracer release as part of the Diapycnal and Isopycnal Mixing Experiment in the Southern Ocean (DIMES). We find that the observational diffusivity estimate of about 500m^2/s at 1500m depth is close to that computed in a data-constrained, 1/20th of a degree simulation of the Drake Passage region. This tracer estimate also agrees with Lagrangian float calculations in the model. The role of mean flow suppression of eddy diffusivity at shallower depths will also be discussed.

Application of a planetary wave breaking parameterization to stratospheric circulation statistics

NASA Technical Reports Server (NTRS)

Randel, William J.; Garcia, Rolando R.

1994-01-01

The planetary wave parameterization scheme developed recently by Garcia is applied to statospheric circulation statistics derived from 12 years of National Meteorological Center operational stratospheric analyses. From the data a planetary wave breaking criterion (based on the ratio of the eddy to zonal mean meridional potential vorticity (PV) gradients), a wave damping rate, and a meridional diffusion coefficient are calculated. The equatorward flank of the polar night jet during winter is identified as a wave breaking region from the observed PV gradients; the region moves poleward with season, covering all high latitudes in spring. Derived damping rates maximize in the subtropical upper stratosphere (the 'surf zone'), with damping time scales of 3-4 days. Maximum diffusion coefficients follow the spatial patterns of the wave breaking criterion, with magnitudes comparable to prior published estimates. Overall, the observed results agree well with the parameterized calculations of Garcia.

ANALYTICAL SOLUTIONS OF THE ATMOSPHERIC DIFFUSION EQUATION WITH MULTIPLE SOURCES AND HEIGHT-DEPENDENT WIND SPEED AND EDDY DIFFUSIVITIES. (R825689C072)

EPA Science Inventory

Abstract

Three-dimensional analytical solutions of the atmospheric diffusion equation with multiple sources and height-dependent wind speed and eddy diffusivities are derived in a systematic fashion. For homogeneous Neumann (total reflection), Dirichlet (total adsorpti...

ANALYTICAL SOLUTIONS OF THE ATMOSPHERIC DIFFUSION EQUATION WITH MULTIPLE SOURCES AND HEIGHT-DEPENDENT WIND SPEED AND EDDY DIFFUSIVITIES. (R825689C048)

EPA Science Inventory

Abstract

Three-dimensional analytical solutions of the atmospheric diffusion equation with multiple sources and height-dependent wind speed and eddy diffusivities are derived in a systematic fashion. For homogeneous Neumann (total reflection), Dirichlet (total adsorpti...

Diffusion-prepared stimulated-echo turbo spin echo (DPsti-TSE): An eddy current-insensitive sequence for three-dimensional high-resolution and undistorted diffusion-weighted imaging.

PubMed

Zhang, Qinwei; Coolen, Bram F; Versluis, Maarten J; Strijkers, Gustav J; Nederveen, Aart J

2017-07-01

In this study, we present a new three-dimensional (3D), diffusion-prepared turbo spin echo sequence based on a stimulated-echo read-out (DPsti-TSE) enabling high-resolution and undistorted diffusion-weighted imaging (DWI). A dephasing gradient in the diffusion preparation module and rephasing gradients in the turbo spin echo module create stimulated echoes, which prevent signal loss caused by eddy currents. Near to perfect agreement of apparent diffusion coefficient (ADC) values between DPsti-TSE and diffusion-weighted echo planar imaging (DW-EPI) was demonstrated in both phantom transient signal experiments and phantom imaging experiments. High-resolution and undistorted DPsti-TSE was demonstrated in vivo in prostate and carotid vessel wall. 3D whole-prostate DWI was achieved with four b values in only 6 min. Undistorted ADC maps of the prostate peripheral zone were obtained at low and high imaging resolutions with no change in mean ADC values [(1.60 ± 0.10) × 10 -3 versus (1.60 ± 0.02) × 10 -3  mm 2 /s]. High-resolution 3D DWI of the carotid vessel wall was achieved in 12 min, with consistent ADC values [(1.40 ± 0.23) × 10 -3  mm 2 /s] across different subjects, as well as slice locations through the imaging volume. This study shows that DPsti-TSE can serve as a robust 3D diffusion-weighted sequence and is an attractive alternative to the traditional two-dimensional DW-EPI approaches. Copyright © 2017 John Wiley & Sons, Ltd.

An investigation of the solar zenith angle variation of D-region ionization

NASA Technical Reports Server (NTRS)

Ratnasiri, P. A. J.; Sechrist, C. F., Jr.

1975-01-01

Model calculations are carried out with a view to interpreting the solar zenith angle variation of D-region ionization. A model is developed for the neutral chemistry including the transport terms relating to molecular and eddy diffusion. The diurnal behavior is described of the minor neutral constituents formed in an oxygen-hydrogen-nitrogen atmosphere, in the height interval between 30 and 120 km. Computations carried out for two cases of the eddy diffusion coefficients models indicate that the constituents which are important for the D-region positive-ion chemistry do not show a significant variation with zenith angle for values up to 75 deg over the D-region heights. In the ion chemistry model, ion-pair production rates are calculated for solar X-rays between 1 A and 100 A, EUV radiations from 100 A up to the Lyman-alpha line, precipitating electrons, and galactic cosmic rays. The solar zenith angle variation of the positive-ion composition, negative-ion composition, and the electron densities are described up to 75 deg zenith angle, in the height interval between 60 and 100 km.

Eddy current-nulled convex optimized diffusion encoding (EN-CODE) for distortion-free diffusion tensor imaging with short echo times.

PubMed

Aliotta, Eric; Moulin, Kévin; Ennis, Daniel B

2018-02-01

To design and evaluate eddy current-nulled convex optimized diffusion encoding (EN-CODE) gradient waveforms for efficient diffusion tensor imaging (DTI) that is free of eddy current-induced image distortions. The EN-CODE framework was used to generate diffusion-encoding waveforms that are eddy current-compensated. The EN-CODE DTI waveform was compared with the existing eddy current-nulled twice refocused spin echo (TRSE) sequence as well as monopolar (MONO) and non-eddy current-compensated CODE in terms of echo time (TE) and image distortions. Comparisons were made in simulations, phantom experiments, and neuro imaging in 10 healthy volunteers. The EN-CODE sequence achieved eddy current compensation with a significantly shorter TE than TRSE (78 versus 96 ms) and a slightly shorter TE than MONO (78 versus 80 ms). Intravoxel signal variance was lower in phantoms with EN-CODE than with MONO (13.6 ± 11.6 versus 37.4 ± 25.8) and not different from TRSE (15.1 ± 11.6), indicating good robustness to eddy current-induced image distortions. Mean fractional anisotropy values in brain edges were also significantly lower with EN-CODE than with MONO (0.16 ± 0.01 versus 0.24 ± 0.02, P < 1 x 10 -5 ) and not different from TRSE (0.16 ± 0.01 versus 0.16 ± 0.01, P = nonsignificant). The EN-CODE sequence eliminated eddy current-induced image distortions in DTI with a TE comparable to MONO and substantially shorter than TRSE. Magn Reson Med 79:663-672, 2018. © 2017 International Society for Magnetic Resonance in Medicine. © 2017 International Society for Magnetic Resonance in Medicine.

Computation of turbulent boundary layer flows with an algebraic stress turbulence model

NASA Technical Reports Server (NTRS)

Kim, Sang-Wook; Chen, Yen-Sen

1986-01-01

An algebraic stress turbulence model is presented, characterized by the following: (1) the eddy viscosity expression is derived from the Reynolds stress turbulence model; (2) the turbulent kinetic energy dissipation rate equation is improved by including a production range time scale; and (3) the diffusion coefficients for turbulence equations are adjusted so that the kinetic energy profile extends further into the free stream region found in most experimental data. The turbulent flow equations were solved using a finite element method. Examples include: fully developed channel flow, fully developed pipe flow, flat plate boundary layer flow, plane jet exhausting into a moving stream, circular jet exhausting into a moving stream, and wall jet flow. Computational results compare favorably with experimental data for most of the examples considered. Significantly improved results were obtained for the plane jet flow, the circular jet flow, and the wall jet flow; whereas the remainder are comparable to those obtained by finite difference methods using the standard kappa-epsilon turbulence model. The latter seems to be promising with further improvement of the expression for the eddy viscosity coefficient.

A B-B-G-K-Y framework for fluid turbulence

NASA Technical Reports Server (NTRS)

Montgomery, D.

1975-01-01

A kinetic theory for fluid turbulence is developed from the Liouville equation and the associated BBGKY hierarchy. Real and imaginary parts of Fourier coefficients of fluid variables play the roles of particles. Closure is achieved by the assumption of negligible five-coefficient correlation functions and probability distributions of Fourier coefficients are the basic variables of the theory. An additional approximation leads to a closed-moment description similar to the so-called eddy-damped Markovian approximation. A kinetic equation is derived for which conservation laws and an H-theorem can be rigorously established, the H-theorem implying relaxation of the absolute equilibrium of Kraichnan. The equation can be cast in the Fokker-Planck form, and relaxation times estimated from its friction and diffusion coefficients. An undetermined parameter in the theory is the free decay time for triplet correlations. Some attention is given to the inclusion of viscous damping and external driving forces.

Quantifying residual, eddy, and mean flow effects on mixing in an idealized circumpolar current

DOE PAGES

Wolfram, Phillip J.; Ringler, Todd D.

2017-07-13

Meridional diffusivity is assessed in this paper for a baroclinically unstable jet in a high-latitudeIdealized Circumpolar Current (ICC) using the Model for Prediction Across Scales-Ocean (MPAS-O) and the online Lagrangian In-situ Global High-performance particle Tracking (LIGHT) diagnostic via space-time dispersion of particle clusters over 120 monthly realizations of O(10 6) particles on 11 potential density surfaces. Diffusivity in the jet reaches values of O(6000 m 2 s -1) and is largest near the critical layer supporting mixing suppression and critical layer theory. Values in the vicinity of the shelf break are suppressed to O(100 m 2 s -1) due tomore » the presence of westward slope front currents. Diffusivity attenuates less rapidly with depth in the jet than both eddy velocity and kinetic energy scalings would suggest. Removal of the mean flow via high-pass filtering shifts the nonlinear parameter (ratio of the eddy velocity to eddy phase speed) into the linear wave regime by increasing the eddy phase speed via the depth-mean flow. Low-pass filtering, in contrast, quantifies the effect of mean shear. Diffusivity is decomposed into mean flow shear, linear waves, and the residual nonhomogeneous turbulence components, where turbulence dominates and eddy-produced filamentation strained by background mean shear enhances mixing, accounting for ≥ 80% of the total diffusivity relative to mean shear [O(100 m 2 s -1)], linear waves [O(1000 m 2 s -1)], and undecomposed full diffusivity [O(6000 m 2 s -1)]. Finally, diffusivity parameterizations accounting for both the nonhomogeneous turbulence residual and depth variability are needed.« less

Improved convection compensating pulsed field gradient spin-echo and stimulated-echo methods.

PubMed

Sørland, G H; Seland, J G; Krane, J; Anthonsen, H W

2000-02-01

The need for convection compensating methods in NMR has been manifested through an increasing number of publications related to the subject over the past few years (J. Magn. Reson. 125, 372 (1997); 132, 13 (1998); 131, 126 (1998); 118, 50 (1996); 133, 379 (1998)). When performing measurements at elevated temperature, small convection currents may give rise to erroneous values of the diffusion coefficient. In work with high resolution NMR spectroscopy, the application of magnetic field gradients also introduces an eddy-current magnetic field which may result in errors in phase and baseline in the FFT-spectra. The eddy current field has been greatly suppressed by the application of bipolar magnetic field gradients. However, when introducing bipolar magnetic field gradients, the pulse sequence is lengthened significantly. This has recently been pointed out as a major drawback because of the loss of coherence and of NMR-signal due to transverse relaxation processes. Here we present modified convection compensating pulsed field gradient double spin echo and double stimulated echo sequences which suppress the eddy-current magnetic field without increasing the duration of the pulse sequences. Copyright 2000 Academic Press.

Characteristic eddy decomposition of turbulence in a channel

NASA Technical Reports Server (NTRS)

Moin, Parviz; Moser, Robert D.

1991-01-01

The proper orthogonal decomposition technique (Lumley's decomposition) is applied to the turbulent flow in a channel to extract coherent structures by decomposing the velocity field into characteristic eddies with random coefficients. In the homogeneous spatial directions, a generaliztion of the shot-noise expansion is used to determine the characteristic eddies. In this expansion, the Fourier coefficients of the characteristic eddy cannot be obtained from the second-order statistics. Three different techniques are used to determine the phases of these coefficients. They are based on: (1) the bispectrum, (2) a spatial compactness requirement, and (3) a functional continuity argument. Results from these three techniques are found to be similar in most respects. The implications of these techniques and the shot-noise expansion are discussed. The dominant eddy is found to contribute as much as 76 percent to the turbulent kinetic energy. In both 2D and 3D, the characteristic eddies consist of an ejection region straddled by streamwise vortices that leave the wall in the very short streamwise distance of about 100 wall units.

Artificial fluid properties for large-eddy simulation of compressible turbulent mixing

NASA Astrophysics Data System (ADS)

Cook, Andrew W.

2007-05-01

An alternative methodology is described for large-eddy simulation (LES) of flows involving shocks, turbulence, and mixing. In lieu of filtering the governing equations, it is postulated that the large-scale behavior of a LES fluid, i.e., a fluid with artificial properties, will be similar to that of a real fluid, provided the artificial properties obey certain constraints. The artificial properties consist of modifications to the shear viscosity, bulk viscosity, thermal conductivity, and species diffusivity of a fluid. The modified transport coefficients are designed to damp out high wavenumber modes, close to the resolution limit, without corrupting lower modes. Requisite behavior of the artificial properties is discussed and results are shown for a variety of test problems, each designed to exercise different aspects of the models. When combined with a tenth-order compact scheme, the overall method exhibits excellent resolution characteristics for turbulent mixing, while capturing shocks and material interfaces in a crisp fashion.

A Physical Mechanism for the Asymmetry in Top-Down and Bottom-Up Diffusion.

NASA Astrophysics Data System (ADS)

Wyngaard, J. C.

1987-04-01

Recent large-eddy simulations of the vertical diffusion of a passive, conservative scalar through the convective boundary layer (CBL) show strikingly different eddy diffusivity profiles in the `top-down' and `bottom-up' cases. These results indicate that for a given turbulent velocity field and associated scalar flux, the mean change in scalar mixing ratio across the CBL is several times larger if the flux originates at the top of the boundary layer (i.e., in top-down diffusion) rather than at the bottom. The large-eddy simulation (LES) data show that this asymmetry is due to a breakdown of the eddy-diffusion concept.A simple updraft-downdraft model of the CBL reveals a physical mechanism that could cause this unexpected behavior. The large, positive skewness of the convectively driven vertical velocity gives an appreciably higher probability of downdrafts than updrafts; this excess probability of downdrafts, interacting with the time changes of the mean mixing ratio caused by the nonstationarity of the bottom-up and top-down diffusion processes, decreases the equilibrium value of mean mixing-ratio jump across the mixed layer in the bottom-up case and increases it in the top-down case. The resulting diffusion asymmetry agrees qualitatively with that found through LES.

A Realizable Reynolds Stress Algebraic Equation Model

NASA Technical Reports Server (NTRS)

Shih, Tsan-Hsing; Zhu, Jiang; Lumley, John L.

1993-01-01

The invariance theory in continuum mechanics is applied to analyze Reynolds stresses in high Reynolds number turbulent flows. The analysis leads to a turbulent constitutive relation that relates the Reynolds stresses to the mean velocity gradients in a more general form in which the classical isotropic eddy viscosity model is just the linear approximation of the general form. On the basis of realizability analysis, a set of model coefficients are obtained which are functions of the time scale ratios of the turbulence to the mean strain rate and the mean rotation rate. The coefficients will ensure the positivity of each component of the mean rotation rate. These coefficients will ensure the positivity of each component of the turbulent kinetic energy - realizability that most existing turbulence models fail to satisfy. Separated flows over backward-facing step configurations are taken as applications. The calculations are performed with a conservative finite-volume method. Grid-independent and numerical diffusion-free solutions are obtained by using differencing schemes of second-order accuracy on sufficiently fine grids. The calculated results are compared in detail with the experimental data for both mean and turbulent quantities. The comparison shows that the present proposal significantly improves the predictive capability of K-epsilon based two equation models. In addition, the proposed model is able to simulate rotational homogeneous shear flows with large rotation rates which all conventional eddy viscosity models fail to simulate.

The dynamics of oceanic fronts. I - The Gulf Stream

NASA Technical Reports Server (NTRS)

Kao, T. W.

1980-01-01

The establishment and maintenance of the mean hydrographic properties of large-scale density fronts in the upper ocean is considered. The dynamics is studied by posing an initial value problem starting with a near-surface discharge of buoyant water with a prescribed density deficit into an ambient stationary fluid of uniform density; full time dependent diffusion and Navier-Stokes equations are then used with constant eddy diffusion and viscosity coefficients, together with a constant Coriolis parameter. Scaling analysis reveals three independent scales of the problem including the radius of deformation of the inertial length, buoyancy length, and diffusive length scales. The governing equations are then suitably scaled and the resulting normalized equations are shown to depend on the Ekman number alone for problems of oceanic interest. It is concluded that the mean Gulf Stream dynamics can be interpreted in terms of a solution of the Navier-Stokes and diffusion equations, with the cross-stream circulation responsible for the maintenance of the front; this mechanism is suggested for the maintenance of the Gulf Stream dynamics.

Modeling and calculation of turbulent lifted diffusion flames

DOE Office of Scientific and Technical Information (OSTI.GOV)

Sanders, J.P.H.; Lamers, A.P.G.G.

1994-01-01

Liftoff heights of turbulent diffusion flames have been modeled using the laminar diffusion flamelet concept of Peters and Williams. The strain rate of the smallest eddies is used as the stretch describing parameter, instead of the more common scalar dissipation rate. The h(U) curve, which is the mean liftoff height as a function of fuel exit velocity can be accurately predicted, while this was impossible with the scalar dissipation rate. Liftoff calculations performed in the flames as well as in the equivalent isothermal jets, using a standard k-[epsilon] turbulence model yield approximately the same correct slope for the h(U) curvemore » while the offset has to be reproduced by choosing an appropriate coefficient in the strain rate model. For the flame calculations a model for the pdf of the fluctuating flame base is proposed. The results are insensitive to its width. The temperature field is qualitatively different from the field calculated by Bradley et al. who used a premixed flamelet model for diffusion flames.« less

Horizontal mixing in the Southern Ocean from Argo float trajectories

NASA Astrophysics Data System (ADS)

Roach, Christopher J.; Balwada, Dhruv; Speer, Kevin

2016-08-01

We provide the first observational estimate of the circumpolar distribution of cross-stream eddy diffusivity at 1000 m in the Southern Ocean using Argo float trajectories. We show that Argo float trajectories, from the float surfacing positions, can be used to estimate lateral eddy diffusivities in the ocean and that these estimates are comparable to those obtained from RAFOS floats, where they overlap. Using the Southern Ocean State Estimate (SOSE) velocity fields to advect synthetic particles with imposed behavior that is "Argo-like" and "RAFOS-like" diffusivity estimates from both sets of synthetic particles agreed closely at the three dynamically very different test sites, the Kerguelen Island region, the Southeast Pacific Ocean, and the Scotia Sea, and support our approach. Observed cross-stream diffusivities at 1000 m, calculated from Argo float trajectories, ranged between 300 and 2500 m2 s-1, with peaks corresponding to topographic features associated with the Scotia Sea, the Kerguelen Plateau, the Campbell Plateau, and the Southeast Pacific Ridge. These observational estimates agree with previous regional estimates from the Diapycnal and Isopycnal Mixing Experiment in the Southern Ocean (DIMES) near the Drake Passage, and other estimates from natural tracers (helium), inverse modeling studies, and current meter measurements. These estimates are also compared to the suppressed eddy diffusivity in the presence of mean flows. The comparison suggests that away from regions of strong topographic steering suppression explains both the structure and magnitude of eddy diffusivity but that eddy diffusivities in the regions of topographic steering are greater than what would be theoretically expected and the ACC experiences localized enhanced cross-stream mixing in these regions.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Wolfram, Phillip J.; Ringler, Todd D.

Meridional diffusivity is assessed in this paper for a baroclinically unstable jet in a high-latitudeIdealized Circumpolar Current (ICC) using the Model for Prediction Across Scales-Ocean (MPAS-O) and the online Lagrangian In-situ Global High-performance particle Tracking (LIGHT) diagnostic via space-time dispersion of particle clusters over 120 monthly realizations of O(10 6) particles on 11 potential density surfaces. Diffusivity in the jet reaches values of O(6000 m 2 s -1) and is largest near the critical layer supporting mixing suppression and critical layer theory. Values in the vicinity of the shelf break are suppressed to O(100 m 2 s -1) due tomore » the presence of westward slope front currents. Diffusivity attenuates less rapidly with depth in the jet than both eddy velocity and kinetic energy scalings would suggest. Removal of the mean flow via high-pass filtering shifts the nonlinear parameter (ratio of the eddy velocity to eddy phase speed) into the linear wave regime by increasing the eddy phase speed via the depth-mean flow. Low-pass filtering, in contrast, quantifies the effect of mean shear. Diffusivity is decomposed into mean flow shear, linear waves, and the residual nonhomogeneous turbulence components, where turbulence dominates and eddy-produced filamentation strained by background mean shear enhances mixing, accounting for ≥ 80% of the total diffusivity relative to mean shear [O(100 m 2 s -1)], linear waves [O(1000 m 2 s -1)], and undecomposed full diffusivity [O(6000 m 2 s -1)]. Finally, diffusivity parameterizations accounting for both the nonhomogeneous turbulence residual and depth variability are needed.« less

Eddy current compensated double diffusion encoded (DDE) MRI.

PubMed

Mueller, Lars; Wetscherek, Andreas; Kuder, Tristan Anselm; Laun, Frederik Bernd

2017-01-01

Eddy currents might lead to image distortions in diffusion-weighted echo planar imaging. A method is proposed to reduce their effects on double diffusion encoding (DDE) MRI experiments and the thereby derived microscopic fractional anisotropy (μFA). The twice-refocused spin echo scheme was adapted for DDE measurements. To assess the effect of individual diffusion encodings on the image distortions, measurements of a grid of plastic rods in water were performed. The effect of eddy current compensation on μFA measurements was evaluated in the brains of six healthy volunteers. The use of an eddy current compensation reduced the signal variation. As expected, the distortions caused by the second encoding were larger than those of the first encoding, entailing a stronger need to compensate for them. For an optimal result, however, both encodings had to be compensated. The artifact reduction strongly improved the measurement of the μFA in ventricles and gray matter by reducing the overestimation. An effect of the compensation on absolute μFA values in white matter was not observed. It is advisable to compensate both encodings in DDE measurements for eddy currents. Magn Reson Med 77:328-335, 2017. © 2015 Wiley Periodicals, Inc. © 2015 Wiley Periodicals, Inc.

Characterization and correction of eddy-current artifacts in unipolar and bipolar diffusion sequences using magnetic field monitoring.

PubMed

Chan, Rachel W; von Deuster, Constantin; Giese, Daniel; Stoeck, Christian T; Harmer, Jack; Aitken, Andrew P; Atkinson, David; Kozerke, Sebastian

2014-07-01

Diffusion tensor imaging (DTI) of moving organs is gaining increasing attention but robust performance requires sequence modifications and dedicated correction methods to account for system imperfections. In this study, eddy currents in the "unipolar" Stejskal-Tanner and the velocity-compensated "bipolar" spin-echo diffusion sequences were investigated and corrected for using a magnetic field monitoring approach in combination with higher-order image reconstruction. From the field-camera measurements, increased levels of second-order eddy currents were quantified in the unipolar sequence relative to the bipolar diffusion sequence while zeroth and linear orders were found to be similar between both sequences. Second-order image reconstruction based on field-monitoring data resulted in reduced spatial misalignment artifacts and residual displacements of less than 0.43 mm and 0.29 mm (in the unipolar and bipolar sequences, respectively) after second-order eddy-current correction. Results demonstrate the need for second-order correction in unipolar encoding schemes but also show that bipolar sequences benefit from second-order reconstruction to correct for incomplete intrinsic cancellation of eddy-currents. Copyright © 2014 The Authors. Published by Elsevier Inc. All rights reserved.

Simulation of Deep Convective Clouds with the Dynamic Reconstruction Turbulence Closure

NASA Astrophysics Data System (ADS)

Shi, X.; Chow, F. K.; Street, R. L.; Bryan, G. H.

2017-12-01

The terra incognita (TI), or gray zone, in simulations is a range of grid spacing comparable to the most energetic eddy diameter. Spacing in mesoscale and simulations is much larger than the eddies, and turbulence is parameterized with one-dimensional vertical-mixing. Large eddy simulations (LES) have grid spacing much smaller than the energetic eddies, and use three-dimensional models of turbulence. Studies of convective weather use convection-permitting resolutions, which are in the TI. Neither mesoscale-turbulence nor LES models are designed for the TI, so TI turbulence parameterization needs to be discussed. Here, the effects of sub-filter scale (SFS) closure schemes on the simulation of deep tropical convection are evaluated by comparing three closures, i.e. Smagorinsky model, Deardorff-type TKE model and the dynamic reconstruction model (DRM), which partitions SFS turbulence into resolvable sub-filter scales (RSFS) and unresolved sub-grid scales (SGS). The RSFS are reconstructed, and the SGS are modeled with a dynamic eddy viscosity/diffusivity model. The RSFS stresses/fluxes allow backscatter of energy/variance via counter-gradient stresses/fluxes. In high-resolution (100m) simulations of tropical convection use of these turbulence models did not lead to significant differences in cloud water/ice distribution, precipitation flux, or vertical fluxes of momentum and heat. When model resolutions are coarsened, the Smagorinsky and TKE models overestimate cloud ice and produces large-amplitude downward heat flux in the middle troposphere (not found in the high-resolution simulations). This error is a result of unrealistically large eddy diffusivities, i.e., the eddy diffusivity of the DRM is on the order of 1 for the coarse resolution simulations, the eddy diffusivity of the Smagorinsky and TKE model is on the order of 100. Splitting the eddy viscosity/diffusivity scalars into vertical and horizontal components by using different length scales and strain rate components helps to reduce the errors, but does not completely remedy the problem. In contrast, the coarse resolution simulations using the DRM produce results that are more consistent with the high-resolution results, suggesting that the DRM is a more appropriate turbulence model for simulating convection in the TI.

A solar dynamo surface wave at the interface between convection and nonuniform rotation

NASA Technical Reports Server (NTRS)

Parker, E. N.

1993-01-01

A simple dynamo surface wave is presented to illustrate the basic principles of a dynamo operating in the thin layer of shear and suppressed eddy diffusion beneath the cyclonic convection in the convection zone of the sun. It is shown that the restriction of the shear delta(Omega)/delta(r) to a region below the convective zone provides the basic mode with a greatly reduced turbulent diffusion coefficient in the region of strong azimuthal field. The dynamo takes on the character of a surface wave tied to the lower surface z = 0 of the convective zone. There is a substantial body of evidence suggesting a fibril state for the principal flux bundles beneath the surface of the sun, with fundamental implications for the solar dynamo.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Na, Ji Sung; Koo, Eunmo; Munoz-Esparza, Domingo

High-resolution large-eddy simulation of the flow over a large wind farm (64 wind turbines) is performed using the HIGRAD/FIRETEC-WindBlade model, which is a high-performance computing wind turbine–atmosphere interaction model that uses the Lagrangian actuator line method to represent rotating turbine blades. These high-resolution large-eddy simulation results are used to parameterize the thrust and power coefficients that contain information about turbine interference effects within the wind farm. Those coefficients are then incorporated into the WRF (Weather Research and Forecasting) model in order to evaluate interference effects in larger-scale models. In the high-resolution WindBlade wind farm simulation, insufficient distance between turbines createsmore » the interference between turbines, including significant vertical variations in momentum and turbulent intensity. The characteristics of the wake are further investigated by analyzing the distribution of the vorticity and turbulent intensity. Quadrant analysis in the turbine and post-turbine areas reveals that the ejection motion induced by the presence of the wind turbines is dominant compared to that in the other quadrants, indicating that the sweep motion is increased at the location where strong wake recovery occurs. Regional-scale WRF simulations reveal that although the turbulent mixing induced by the wind farm is partly diffused to the upper region, there is no significant change in the boundary layer depth. The velocity deficit does not appear to be very sensitive to the local distribution of turbine coefficients. However, differences of about 5% on parameterized turbulent kinetic energy were found depending on the turbine coefficient distribution. Furthermore, turbine coefficients that consider interference in the wind farm should be used in wind farm parameterization for larger-scale models to better describe sub-grid scale turbulent processes.« less

Turbulent kinetics of a large wind farm and their impact in the neutral boundary layer

DOE PAGES

Na, Ji Sung; Koo, Eunmo; Munoz-Esparza, Domingo; ...

2015-12-28

High-resolution large-eddy simulation of the flow over a large wind farm (64 wind turbines) is performed using the HIGRAD/FIRETEC-WindBlade model, which is a high-performance computing wind turbine–atmosphere interaction model that uses the Lagrangian actuator line method to represent rotating turbine blades. These high-resolution large-eddy simulation results are used to parameterize the thrust and power coefficients that contain information about turbine interference effects within the wind farm. Those coefficients are then incorporated into the WRF (Weather Research and Forecasting) model in order to evaluate interference effects in larger-scale models. In the high-resolution WindBlade wind farm simulation, insufficient distance between turbines createsmore » the interference between turbines, including significant vertical variations in momentum and turbulent intensity. The characteristics of the wake are further investigated by analyzing the distribution of the vorticity and turbulent intensity. Quadrant analysis in the turbine and post-turbine areas reveals that the ejection motion induced by the presence of the wind turbines is dominant compared to that in the other quadrants, indicating that the sweep motion is increased at the location where strong wake recovery occurs. Regional-scale WRF simulations reveal that although the turbulent mixing induced by the wind farm is partly diffused to the upper region, there is no significant change in the boundary layer depth. The velocity deficit does not appear to be very sensitive to the local distribution of turbine coefficients. However, differences of about 5% on parameterized turbulent kinetic energy were found depending on the turbine coefficient distribution. Furthermore, turbine coefficients that consider interference in the wind farm should be used in wind farm parameterization for larger-scale models to better describe sub-grid scale turbulent processes.« less

Effects of Southern Hemisphere Wind Changes on the Meridional Overturning Circulation in Ocean Models.

PubMed

Gent, Peter R

2016-01-01

Observations show that the Southern Hemisphere zonal wind stress maximum has increased significantly over the past 30 years. Eddy-resolving ocean models show that the resulting increase in the Southern Ocean mean flow meridional overturning circulation (MOC) is partially compensated by an increase in the eddy MOC. This effect can be reproduced in the non-eddy-resolving ocean component of a climate model, providing the eddy parameterization coefficient is variable and not a constant. If the coefficient is a constant, then the Southern Ocean mean MOC change is balanced by an unrealistically large change in the Atlantic Ocean MOC. Southern Ocean eddy compensation means that Southern Hemisphere winds cannot be the dominant mechanism driving midlatitude North Atlantic MOC variability.

On the use of water phantom images to calibrate and correct eddy current induced artefacts in MR diffusion tensor imaging.

PubMed

Bastin, M E; Armitage, P A

2000-07-01

The accurate determination of absolute measures of diffusion anisotropy in vivo using single-shot, echo-planar imaging techniques requires the acquisition of a set of high signal-to-noise ratio, diffusion-weighted images that are free from eddy current induced image distortions. Such geometric distortions can be characterized and corrected in brain imaging data using magnification (M), translation (T), and shear (S) distortion parameters derived from separate water phantom calibration experiments. Here we examine the practicalities of using separate phantom calibration data to correct high b-value diffusion tensor imaging data by investigating the stability of these distortion parameters, and hence the eddy currents, with time. It is found that M, T, and S vary only slowly with time (i.e., on the order of weeks), so that calibration scans need not be performed after every patient examination. This not only minimises the scan time required to collect the calibration data, but also the computational time needed to characterize these eddy current induced distortions. Examples of how measurements of diffusion anisotropy are improved using this post-processing scheme are also presented.

North Atlantic Coast Comprehensive Study (NACCS) Coastal Storm Model Simulations: Waves and Water Levels

DTIC Science & Technology

2015-08-01

published in the NGA’s DNCs, with distinct values assigned to areas of sand, gravel, clay , etc. ERDC/CHL TR-15-14 94 6.5.2 Lateral eddy viscosity As with...6.5.1 Manning’s n bottom friction coefficient ................................................................... 93 6.5.2 Lateral eddy viscosity ...this study include (1) Manning’s n bottom friction coefficient, (2) lateral eddy viscosity , (3) land cover effects on winds (also referred to as

The Martian climate and energy balance models with CO2/H2O atmospheres

NASA Technical Reports Server (NTRS)

Hoffert, M. I.

1986-01-01

The analysis begins with a seasonal energy balance model (EBM) for Mars. This is used to compute surface temperature versus x = sin(latitude) and time over the seasonal cycle. The core model also computes the evolving boundaries of the CO2 icecaps, net sublimational/condensation rates, and the resulting seasonal pressure wave. Model results are compared with surface temperature and pressure history data at Viking lander sites, indicating fairly good agreement when meridional heat transport is represented by a thermal diffusion coefficient D approx. 0.015 W/sq. m/K. Condensational wind distributions are also computed. An analytic model of Martian wind circulation is then proposed, as an extension of the EMB, which incorporates vertical wind profiles containing an x-dependent function evaluated by substitution in the equation defining the diffusion coefficient. This leads to a parameterization of D(x) and of the meridional circulation which recovers the high surface winds predicted by dynamic Mars atmosphere models (approx. 10 m/sec). Peak diffusion coefficients, D approx. 0.6 w/sq m/K, are found over strong Hadley zones - some 40 times larger than those of high-latitude baroclinic eddies. When the wind parameterization is used to find streamline patterns over Martian seasons, the resulting picture shows overturning hemispheric Hadley cells crossing the equator during solstices, and attaining peak intensities during the south summer dust storm season, while condensational winds are most important near the polar caps.

Chemical transport models: the combined non-local diffusion and mixing schemes, and calculation of in-canopy resistance for dry deposition fluxes.

PubMed

Mihailovic, Dragutin T; Alapaty, Kiran; Podrascanin, Zorica

2009-03-01

Improving the parameterization of processes in the atmospheric boundary layer (ABL) and surface layer, in air quality and chemical transport models. To do so, an asymmetrical, convective, non-local scheme, with varying upward mixing rates is combined with the non-local, turbulent, kinetic energy scheme for vertical diffusion (COM). For designing it, a function depending on the dimensionless height to the power four in the ABL is suggested, which is empirically derived. Also, we suggested a new method for calculating the in-canopy resistance for dry deposition over a vegetated surface. The upward mixing rate forming the surface layer is parameterized using the sensible heat flux and the friction and convective velocities. Upward mixing rates varying with height are scaled with an amount of turbulent kinetic energy in layer, while the downward mixing rates are derived from mass conservation. The vertical eddy diffusivity is parameterized using the mean turbulent velocity scale that is obtained by the vertical integration within the ABL. In-canopy resistance is calculated by integration of inverse turbulent transfer coefficient inside the canopy from the effective ground roughness length to the canopy source height and, further, from its the canopy height. This combination of schemes provides a less rapid mass transport out of surface layer into other layers, during convective and non-convective periods, than other local and non-local schemes parameterizing mixing processes in the ABL. The suggested method for calculating the in-canopy resistance for calculating the dry deposition over a vegetated surface differs remarkably from the commonly used one, particularly over forest vegetation. In this paper, we studied the performance of a non-local, turbulent, kinetic energy scheme for vertical diffusion combined with a non-local, convective mixing scheme with varying upward mixing in the atmospheric boundary layer (COM) and its impact on the concentration of pollutants calculated with chemical and air-quality models. In addition, this scheme was also compared with a commonly used, local, eddy-diffusivity scheme. Simulated concentrations of NO2 by the COM scheme and new parameterization of the in-canopy resistance are closer to the observations when compared to those obtained from using the local eddy-diffusivity scheme. Concentrations calculated with the COM scheme and new parameterization of in-canopy resistance, are in general higher and closer to the observations than those obtained by the local, eddy-diffusivity scheme (on the order of 15-22%). To examine the performance of the scheme, simulated and measured concentrations of a pollutant (NO2) were compared for the years 1999 and 2002. The comparison was made for the entire domain used in simulations performed by the chemical European Monitoring and Evaluation Program Unified model (version UNI-ACID, rv2.0) where schemes were incorporated.

Distribution of Upper Circumpolar Deep Water on the warming continental shelf of the West Antarctic Peninsula

NASA Astrophysics Data System (ADS)

Couto, Nicole; Martinson, Douglas G.; Kohut, Josh; Schofield, Oscar

2017-07-01

We use autonomous underwater vehicles to characterize the spatial distribution of Upper Circumpolar Deep Water (UCDW) on the continental shelf of the West Antarctic Peninsula (WAP) and present the first near-synoptic measurements of mesoscale features (eddies) containing UCDW on the WAP. Thirty-three subsurface eddies with widths on the order of 10 km were detected during four glider deployments. Each eddy contributed an average of 5.8 × 1016 J to the subpycnocline waters, where a cross-shelf heat flux of 1.37 × 1019 J yr-1 is required to balance the diffusive loss of heat to overlying winter water and to the near-coastal waters. Approximately two-thirds of the heat coming onto the shelf diffuses across the pycnocline and one-third diffuses to the coastal waters; long-term warming of the subpycnocline waters is a small residual of this balance. Sixty percent of the profiles that contained UCDW were part of a coherent eddy. Between 20% and 53% of the lateral onshore heat flux to the WAP can be attributed to eddies entering Marguerite Trough, a feature in the southern part of the shelf which is known to be an important conduit for UCDW. A northern trough is identified as additional important location for eddy intrusion.

Vertical Eddy Diffusivity as a Control Parameter in the Tropical Pacific Ocean

NASA Astrophysics Data System (ADS)

Martinez Avellaneda, N.; Cornuelle, B.; Mazloff, M. R.; Stammer, D.

2012-12-01

Ocean models suffer from errors in the treatment of turbulent sub-grid scale motions causing mixing and energy dissipation. Unrealistic small-scale features in models can have large-scale consequences, such as biases in the upper ocean temperature, a symptom of poorly-simulated upwelling, currents and air-sea interactions. This is of special importance in the tropical Pacific Ocean, which is home to energetic air-sea interactions that affect global climate. It has been shown in a number of studies that the simulated ENSO variability is highly dependent on the state of the ocean (e.g.: background mixing). Moreover, the magnitude of the vertical numerical diffusion is of primary importance in properly reproducing the Pacific equatorial thermocline. Yet, it is a common practice to use spatially uniform mixing parameters in ocean simulations. This work is part of a NASA-funded project to estimate the space-varying ocean mixing coefficients in an eddy-permitting model of the tropical Pacific. The usefulness of assimilation techniques in estimating mixing parameters has been previously explored (e.g.: Stammer, 2005, Ferreira et al., 2005). The authors also demonstrated that the spatial structure of the Equatorial Undercurrent (EUC) could be improved by adjusting wind-stress and surface buoyancy flux within their error bounds. In our work, we address the important question of whether adjusting mixing parameterizations can bring about similar improvements. To that end, an eddy-permitting state estimate for the tropical Pacific is developed using the MIT general circulation model and its adjoint where the vertical diffusivity is set as a control parameter. Complementary adjoint-based sensitivity results show strong sensitivities of the Tropical Pacific thermocline (thickness and location) and the EUC transport to the vertical diffusivity in the tropics. Argo, CTD, XBT and mooring in-situ data, as well as TMI SST and altimetry observations are assimilated in order to reduce the misfit between the model simulations and the ocean observations. Model domain topography of 1/3dgr of spatial resolution interpolated from ETOPO 2. The first and the last color levels represent regions shallower than 100m and deeper than 5000m, respectively

Turbulent flow separation in three-dimensional asymmetric diffusers

NASA Astrophysics Data System (ADS)

Jeyapaul, Elbert

2011-12-01

Turbulent three-dimensional flow separation is more complicated than 2-D. The physics of the flow is not well understood. Turbulent flow separation is nearly independent of the Reynolds number, and separation in 3-D occurs at singular points and along convergence lines emanating from these points. Most of the engineering turbulence research is driven by the need to gain knowledge of the flow field that can be used to improve modeling predictions. This work is motivated by the need for a detailed study of 3-D separation in asymmetric diffusers, to understand the separation phenomena using eddy-resolving simulation methods, assess the predictability of existing RANS turbulence models and propose modeling improvements. The Cherry diffuser has been used as a benchmark. All existing linear eddy-viscosity RANS models k--o SST,k--epsilon and v2- f fail in predicting such flows, predicting separation on the wrong side. The geometry has a doubly-sloped wall, with the other two walls orthogonal to each other and aligned with the diffuser inlet giving the diffuser an asymmetry. The top and side flare angles are different and this gives rise to different pressure gradient in each transverse direction. Eddyresolving simulations using the Scale adaptive simulation (SAS) and Large Eddy Simulation (LES) method have been used to predict separation in benchmark diffuser and validated. A series of diffusers with the same configuration have been generated, each having the same streamwise pressure gradient and parametrized only by the inlet aspect ratio. The RANS models were put to test and the flow physics explored using SAS-generated flow field. The RANS model indicate a transition in separation surface from top sloped wall to the side sloped wall at an inlet aspect ratio much lower than observed in LES results. This over-sensitivity of RANS models to transverse pressure gradients is due to lack of anisotropy in the linear Reynolds stress formulation. The complexity of the flow separation is due to effects of lateral straining, streamline curvature, secondary flow of second kind, transverse pressure gradient on turbulence. Resolving these effects is possible with anisotropy turbulence models as the Explicit Algebraic Reynolds stress model (EARSM). This model has provided accurate prediction of streamwise and transverse velocity, however the wall pressure is under predicted. An improved EARSM model is developed by correcting the coefficients, which predicts a more accurate wall pressure. There exists scope for improvement of this model, by including convective effects and dynamics of velocity gradient invariants.

Ice haze, snow, and the Mars water cycle

NASA Technical Reports Server (NTRS)

Kahn, Ralph

1990-01-01

Light curves and extinction profiles derived from Martian limb observations are used to constrain the atmospheric temperature structure in regions of the atmosphere with thin haze and to analyze the haze particle properties and atmospheric eddy mixing. Temperature between 170 and 190 K are obtained for three cases at levels in the atmosphere ranging from 20 to 50 km. Eddy diffusion coefficients around 100,000 sq cm/s, typical of a nonconvecting atmosphere, are derived in the haze regions at times when the atmosphere is relatively clear of dust. This parameter apparently changes by more than three orders of magnitude with season and local conditions. The derived particle size parameter varies systematically by more than an order of magnitude with condensation level, in such a way that the characteristic fall time is always about one Martian day. Ice hazes provide a mechanism for scavenging water vapor in the thin Mars atmosphere and may play a key role in the seasonal cycle of water on Mars.

Impact of straight, unconnected, radially-oriented, and tapered mesopores on column efficiency: A theoretical investigation.

PubMed

Gritti, Fabrice

2017-02-17

Superficially porous particles (SPPs) can be prepared from a pseudomorphic transformation (PMT) which produces straight, unconnected, and radially-oriented mesopores (ROMs). ROMs can be either both ends open in fully porous particles (FPPs) or one-end-closed in SPPs. The impact of ROMs on the longitudinal diffusion (B/u), solid-liquid mass transfer resistance (C s u), and on the eddy dispersion (A(u)) height equivalent to a theoretical plate (HETP) of 3D randomly packed columns was investigated based on theoretical viewpoints. Torquato's theory of effective diffusion in packed beds (B term), Giddings' coupling theory of eddy dispersion (A term), and Giddings' generalized nonequilibrium theory (C s term) are applied to make predictions. First, it is found that the A term is nearly independent on the internal structure of the particle. Secondly, in the absence of flow, infinitely narrow and both ends open (no constriction effect) ROMs induce an internal hindrance factor of 23 regarding diffusion along the axial direction. Experimental data reveal that one-end-closed and 80Å wide ROMs in SPPs lead to a measurable internal hindrance factor of 27 regarding diffusion in the porous shell. Thirdly, above the optimum speed, the C s coefficient is dependent on the geometry (cylinders, cones, etc.) of the ROMs: when ROMs are conical in SPPs, C s is expected to decrease by 80% with respect to cylindrical ROMs. From an application perspective, PMT-SPPs prepared with narrow ROMs are well suited for the analysis of small molecules at or below optimum speed (lowest B term) while PMT-SPPs made of wide and conical ROMs are ideal for the analysis of large molecules above optimum speed (smallest C s term). Copyright © 2017. Published by Elsevier B.V.

Enhancement of diffusive transport in oscillatory flows

NASA Technical Reports Server (NTRS)

Knobloch, E.; Merryfield, W. J.

1992-01-01

The theory of transport of passive scalars in oscillatory flows is reexamined. The differences between transport in standing and traveling waves are emphasized. Both Lagrangian and Eulerian diffusivities are calculated, and the conditions for their applicability are discussed. Numerical simulations are conducted to understand the expulsion of gradients from time-dependent eddies and the resulting transport. The results indicate that it is the Eulerian diffusivity that is of primary relevance for describing enhanced transport on spatial scales larger than that of the eddies.

Characteristic-eddy decomposition of turbulence in a channel

NASA Technical Reports Server (NTRS)

Moin, Parviz; Moser, Robert D.

1989-01-01

Lumley's proper orthogonal decomposition technique is applied to the turbulent flow in a channel. Coherent structures are extracted by decomposing the velocity field into characteristic eddies with random coefficients. A generalization of the shot-noise expansion is used to determine the characteristic eddies in homogeneous spatial directions. Three different techniques are used to determine the phases of the Fourier coefficients in the expansion: (1) one based on the bispectrum, (2) a spatial compactness requirement, and (3) a functional continuity argument. Similar results are found from each of these techniques.

Electric discharge synthesis of HCN in simulated Jovian atmospheres

NASA Technical Reports Server (NTRS)

Stribling, Roscoe; Miller, Stanley L.

1987-01-01

Corona discharge is presently considered as a possible source of the HCN detected in the Jovian atmosphere at 2.2 x 10 to the -7th moles/sq cm column density, for the cases of gas mixtures containing H2, CH4, and NH3, with H2/CH4 ratios from 4.4 to 1585. A 3:1 ratio of corona discharge to lightning energy similar to that of the earth is applied to Jupiter. Depending on the lightning energy available on Jupiter and the eddy diffusion coefficients in the synthesis region, HCN column densities generated by corona discharge could account for about 10 percent of the HCN observed.

Turbulent mixing within the Kuroshio in the Tokara Strait

NASA Astrophysics Data System (ADS)

Tsutsumi, Eisuke; Matsuno, Takeshi; Lien, Ren-Chieh; Nakamura, Hirohiko; Senjyu, Tomoharu; Guo, Xinyu

2017-09-01

Turbulent mixing and background current were observed using a microstructure profiler and acoustic Doppler current profilers in the Tokara Strait, where many seamounts and small islands exist within the route of the Kuroshio in the East China Sea. Vertical structure and water properties of the Kuroshio were greatly modified downstream from shallow seamounts. In the lee of a seamount crest at 200 m depth, the modification made the flow tend to shear instability, and the vertical eddy diffusivity is enhanced by nearly 100 times that of the upstream site, to Kρ ˜ O(10-3)-O(10-2) m2 s-1. A one-dimensional diffusion model using the observed eddy diffusivity reproduced the observed downstream evolution of the temperature-salinity profile. However, the estimated diffusion time-scale is at least 10 times longer than the observed advection time-scale. This suggests that the eddy diffusivity reaches to O(10-1) m2 s-1 in the vicinity of the seamount. At a site away from the abrupt topography, eddy diffusivity was also elevated to O(10-3) m2 s-1, and was associated with shear instability presumably induced by the Kuroshio shear and near-inertial internal-wave shear. Our study suggests that a better prediction of current, water-mass properties, and nutrients within the Kuroshio requires accurate understanding and parameterization of flow-topography interaction such as internal hydraulics, the associated internal-wave processes, and turbulent mixing processes.

A submesoscale coherent vortex in the Ligurian Sea: From dynamical barriers to biological implications

NASA Astrophysics Data System (ADS)

Bosse, Anthony; Testor, Pierre; Mayot, Nicolas; Prieur, Louis; D'Ortenzio, Fabrizio; Mortier, Laurent; Le Goff, Hervé; Gourcuff, Claire; Coppola, Laurent; Lavigne, Héloïse; Raimbault, Patrick

2017-08-01

In June 2013, a glider equipped with oxygen and fluorescence sensors has been used to extensively sample an anticyclonic Submesoscale Coherent Vortex (SCV) in the Ligurian Sea (NW Mediterranean Sea). Those measurements are complemented by full-depth CTD casts (T, S, and oxygen) and water samples documenting nutrients and phytoplankton pigments within the SCV and outside. The SCV has a very homogeneous core of oxygenated waters between 300 and 1200 m formed 4.5 months earlier during the winter deep convection event. It has a strong dynamical signature with peak velocities at 700 m depth of 13.9 cm s-1 in cyclogeostrophic balance. The eddy has a small radius of 6.2 km corresponding to high Rossby number of -0.45. The vorticity at the eddy center reaches -0.8f. Cross-stream isopycnic diffusion of tracers between the eddy core and the surroundings is found to be very limited due to dynamical barriers set by the SCV associated with a diffusivity coefficient of about 0.2 m2 s-1. The deep core is nutrients-depleted with concentrations of nitrate, phosphate, and silicate, 13-18% lower than the rich surrounding waters. However, the nutriclines are shifted of about 20-50 m toward the surface thus increasing the nutrients availability for phytoplankton. Chlorophyll-a concentrations at the deep chlorophyll maximum are subsequently about twice bigger as compared to outside. Pigments further reveal the predominance of nanophytoplankton inside the eddy and an enhancement of the primary productivity. This study demonstrates the important impact of postconvective SCVs on nutrients distribution and phytoplankton community, as well as on the subsequent primary production and carbon sequestration.