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.

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.

The quantitative impact of the mesopore size on the mass transfer mechanism of the new 1.9μm fully porous Titan-C18 particles. I: analysis of small molecules.

PubMed

Gritti, Fabrice; Guiochon, Georges

2015-03-06

Previous data have shown that could deliver a minimum reduced plate height as small as 1.7. Additionally, the reduction of the mesopore size after C18 derivatization and the subsequent restriction for sample diffusivity across the Titan-C18 particles were found responsible for the unusually small value of the experimental optimum reduced velocity (5 versus 10 for conventional particles) and for the large values of the average reduced solid-liquid mass transfer resistance coefficients (0.032 versus 0.016) measured for a series of seven n-alkanophenones. The improvements in column efficiency made by increasing the average mesopore size of the Titan silica from 80 to 120Å are investigated from a quantitative viewpoint based on the accurate measurements of the reduced coefficients (longitudinal diffusion, trans-particle mass transfer resistance, and eddy diffusion) and of the intra-particle diffusivity, pore, and surface diffusion for the same series of n-alkanophenone compounds. The experimental results reveal an increase (from 0% to 30%) of the longitudinal diffusion coefficients for the same sample concentration distribution (from 0.25 to 4) between the particle volume and the external volume of the column, a 40% increase of the intra-particle diffusivity for the same sample distribution (from 1 to 7) between the particle skeleton volume and the bulk phase, and a 15-30% decrease of the solid-liquid mass transfer coefficient for the n-alkanophenone compounds. Pore and surface diffusion are increased by 60% and 20%, respectively. The eddy dispersion term and the maximum column efficiency (295000plates/m) remain virtually unchanged. The rate of increase of the total plate height with increasing the chromatographic speed is reduced by 20% and it is mostly controlled (75% and 70% for 80 and 120Å pore size) by the flow rate dependence of the eddy dispersion term. Copyright © 2015 Elsevier B.V. All rights reserved.

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.

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.

Computational investigation of longitudinal diffusion, eddy dispersion, and trans-particle mass transfer in bulk, random packings of core-shell particles with varied shell thickness and shell diffusion coefficient.

PubMed

Daneyko, Anton; Hlushkou, Dzmitry; Baranau, Vasili; Khirevich, Siarhei; Seidel-Morgenstern, Andreas; Tallarek, Ulrich

2015-08-14

In recent years, chromatographic columns packed with core-shell particles have been widely used for efficient and fast separations at comparatively low operating pressure. However, the influence of the porous shell properties on the mass transfer kinetics in core-shell packings is still not fully understood. We report on results obtained with a modeling approach to simulate three-dimensional advective-diffusive transport in bulk random packings of monosized core-shell particles, covering a range of reduced mobile phase flow velocities from 0.5 up to 1000. The impact of the effective diffusivity of analyte molecules in the porous shell and the shell thickness on the resulting plate height was investigated. An extension of Giddings' theory of coupled eddy dispersion to account for retention of analyte molecules due to stagnant regions in porous shells with zero mobile phase flow velocity is presented. The plate height equation involving a modified eddy dispersion term excellently describes simulated data obtained for particle-packings with varied shell thickness and shell diffusion coefficient. It is confirmed that the model of trans-particle mass transfer resistance of core-shell particles by Kaczmarski and Guiochon [42] is applicable up to a constant factor. We analyze individual contributions to the plate height from different mass transfer mechanisms in dependence of the shell parameters. The simulations demonstrate that a reduction of plate height in packings of core-shell relative to fully porous particles arises mainly due to reduced trans-particle mass transfer resistance and transchannel eddy dispersion. Copyright © 2015 Elsevier B.V. All rights reserved.

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.

Dynamics, Stability, and Evolutionary Patterns of Mesoscale Intrathermocline Vortices

DTIC Science & Technology

2016-12-01

physical oceanography, namely, the link between the basin-scale forcing of the ocean by air-sea fluxes and the dissipation of energy and thermal variance...at the microscale. 14. SUBJECT TERMS Meddy, intrathermocline, double diffusion, energy cascade, eddy, MITgcm, numerical simulation, interleaving...lateral intrusions, lateral diffusivity, heat flux 15. NUMBER OF PAGES 69 16. PRICE CODE 17. SECURITY CLASSIFICATION OF REPORT Unclassified 18

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...

Reynolds-Stress Budgets in an Impinging Shock Wave/Boundary-Layer Interaction

NASA Technical Reports Server (NTRS)

Vyas, Manan A.; Yoder, Dennis A.; Gaitonde, Datta V.

2018-01-01

Implicit large-eddy simulation (ILES) of a shock wave/boundary-layer interaction (SBLI) was performed. Comparisons with experimental data showed a sensitivity of the current prediction to the modeling of the sidewalls. This was found to be common among various computational studies in the literature where periodic boundary conditions were used in the spanwise direction, as was the case in the present work. Thus, although the experiment was quasi-two-dimensional, the present simulation was determined to be two-dimensional. Quantities present in the exact equation of the Reynolds-stress transport, i.e., production, molecular diffusion, turbulent transport, pressure diffusion, pressure strain, dissipation, and turbulent mass flux were calculated. Reynolds-stress budgets were compared with past large-eddy simulation and direct numerical simulation datasets in the undisturbed portion of the turbulent boundary layer to validate the current approach. The budgets in SBLI showed the growth in the production term for the primary normal stress and energy transfer mechanism was led by the pressure strain term in the secondary normal stresses. The pressure diffusion term, commonly assumed as negligible by turbulence model developers, was shown to be small but non-zero in the normal stress budgets, however it played a key role in the primary shear stress budget.

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.

Turbulent eddy diffusion models in exposure assessment - Determination of the eddy diffusion coefficient.

PubMed

Shao, Yuan; Ramachandran, Sandhya; Arnold, Susan; Ramachandran, Gurumurthy

2017-03-01

The use of the turbulent eddy diffusion model and its variants in exposure assessment is limited due to the lack of knowledge regarding the isotropic eddy diffusion coefficient, D T . But some studies have suggested a possible relationship between D T and the air changes per hour (ACH) through a room. The main goal of this study was to accurately estimate D T for a range of ACH values by minimizing the difference between the concentrations measured and predicted by eddy diffusion model. We constructed an experimental chamber with a spatial concentration gradient away from the contaminant source, and conducted 27 3-hr long experiments using toluene and acetone under different air flow conditions (0.43-2.89 ACHs). An eddy diffusion model accounting for chamber boundary, general ventilation, and advection was developed. A mathematical expression for the slope based on the geometrical parameters of the ventilation system was also derived. There is a strong linear relationship between D T and ACH, providing a surrogate parameter for estimating D T in real-life settings. For the first time, a mathematical expression for the relationship between D T and ACH has been derived that also corrects for non-ideal conditions, and the calculated value of the slope between these two parameters is very close to the experimentally determined value. The values of D T obtained from the experiments are generally consistent with values reported in the literature. They are also independent of averaging time of measurements, allowing for comparison of values obtained from different measurement settings. These findings make the use of turbulent eddy diffusion models for exposure assessment in workplace/indoor environments more practical.

Large eddy simulation on buoyant gas diffusion near building

NASA Astrophysics Data System (ADS)

Tominaga, Yoshihide; Murakami, Shuzo; Mochida, Akashi

1992-12-01

Large eddy simulations on turbulent diffusion of buoyant gases near a building model are carried out for three cases in which the densimetric Froude Number (Frd) was specified at - 8.6, zero and 8.6 respectively. The accuracy of these simulations is examined by comparing the numerically predicted results with wind tunnel experiments conducted. Two types of sub-grid scale models, the standard Smagorinsky model (type 1) and the modified Smagorinsky model (type 2) are compared. The former does not take account of the production of subgrid energy by buoyancy force but the latter incorporates this effect. The latter model (type 2) gives more accurate results than those given by the standard Smagorinsky model (type 1) in terms of the distributions of kappa greater than sign C less than sign greater than sign C(sup - 2) less than sign.

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

Matrix exponential-based closures for the turbulent subgrid-scale stress tensor.

PubMed

Li, Yi; Chevillard, Laurent; Eyink, Gregory; Meneveau, Charles

2009-01-01

Two approaches for closing the turbulence subgrid-scale stress tensor in terms of matrix exponentials are introduced and compared. The first approach is based on a formal solution of the stress transport equation in which the production terms can be integrated exactly in terms of matrix exponentials. This formal solution of the subgrid-scale stress transport equation is shown to be useful to explore special cases, such as the response to constant velocity gradient, but neglecting pressure-strain correlations and diffusion effects. The second approach is based on an Eulerian-Lagrangian change of variables, combined with the assumption of isotropy for the conditionally averaged Lagrangian velocity gradient tensor and with the recent fluid deformation approximation. It is shown that both approaches lead to the same basic closure in which the stress tensor is expressed as the matrix exponential of the resolved velocity gradient tensor multiplied by its transpose. Short-time expansions of the matrix exponentials are shown to provide an eddy-viscosity term and particular quadratic terms, and thus allow a reinterpretation of traditional eddy-viscosity and nonlinear stress closures. The basic feasibility of the matrix-exponential closure is illustrated by implementing it successfully in large eddy simulation of forced isotropic turbulence. The matrix-exponential closure employs the drastic approximation of entirely omitting the pressure-strain correlation and other nonlinear scrambling terms. But unlike eddy-viscosity closures, the matrix exponential approach provides a simple and local closure that can be derived directly from the stress transport equation with the production term, and using physically motivated assumptions about Lagrangian decorrelation and upstream isotropy.

Effects of Mesoscale Eddies in the Active Mixed Layer: Test of the Parametrisation in Eddy Resolving Simulations

NASA Technical Reports Server (NTRS)

Luneva, M. V.; Clayson, C. A.; Dubovikov, Mikhail

2015-01-01

In eddy resolving simulations, we test a mixed layer mesoscale parametrisation, developed recently by Canuto and Dubovikov [Ocean Model., 2011, 39, 200-207]. With no adjustable parameters, the parametrisation yields the horizontal and vertical mesoscale fluxes in terms of coarse-resolution fields and eddy kinetic energy (EKE). We compare terms of the parametrisation diagnosed from coarse-grained fields with the eddy mesoscale fluxes diagnosed directly from the high resolution model. An expression for the EKE in terms of mean fields has also been found to get a closed parametrisation in terms of the mean fields only. In 40 numerical experiments we simulated two types of flows: idealised flows driven by baroclinic instabilities only, and more realistic flows, driven by wind and surface fluxes as well as by inflow-outflow. The diagnosed quasi-instantaneous horizontal and vertical mesoscale buoyancy fluxes (averaged over 1-2 degrees and 10 days) demonstrate a strong scatter typical for turbulent flows, however, the fluxes are positively correlated with the parametrisation with higher (0.5-0.74) correlations at the experiments with larger baroclinic radius Rossby. After being averaged over 3-4 months, diffusivities diagnosed from the eddy resolving simulations are consistent with the parametrisation for a broad range of parameters. Diagnosed vertical mesoscale fluxes restratify mixed layer and are in a good agreement with the parametrisation unless vertical turbulent mixing in the upper layer becomes strong enough in comparison with mesoscale advection. In the latter case, numerical simulations demonstrate that the deviation of the fluxes from the parametrisation is controlled by dimensionless parameter estimating the ratio of vertical turbulent mixing term to mesoscale advection. An analysis using a modified omega-equation reveals that the effects of the vertical mixing of vorticity is responsible for the two-three fold amplification of vertical mesoscale flux. Possible physical mechanisms, responsible for the amplification of vertical mesoscale flux are discussed.

An Eddy-Diffusivity Mass-flux (EDMF) closure for the unified representation of cloud and convective processes

NASA Astrophysics Data System (ADS)

Tan, Z.; Schneider, T.; Teixeira, J.; Lam, R.; Pressel, K. G.

2014-12-01

Sub-grid scale (SGS) closures in current climate models are usually decomposed into several largely independent parameterization schemes for different cloud and convective processes, such as boundary layer turbulence, shallow convection, and deep convection. These separate parameterizations usually do not converge as the resolution is increased or as physical limits are taken. This makes it difficult to represent the interactions and smooth transition among different cloud and convective regimes. Here we present an eddy-diffusivity mass-flux (EDMF) closure that represents all sub-grid scale turbulent, convective, and cloud processes in a unified parameterization scheme. The buoyant updrafts and precipitative downdrafts are parameterized with a prognostic multiple-plume mass-flux (MF) scheme. The prognostic term for the mass flux is kept so that the life cycles of convective plumes are better represented. The interaction between updrafts and downdrafts are parameterized with the buoyancy-sorting model. The turbulent mixing outside plumes is represented by eddy diffusion, in which eddy diffusivity (ED) is determined from a turbulent kinetic energy (TKE) calculated from a TKE balance that couples the environment with updrafts and downdrafts. Similarly, tracer variances are decomposed consistently between updrafts, downdrafts and the environment. The closure is internally coupled with a probabilistic cloud scheme and a simple precipitation scheme. We have also developed a relatively simple two-stream radiative scheme that includes the longwave (LW) and shortwave (SW) effects of clouds, and the LW effect of water vapor. We have tested this closure in a single-column model for various regimes spanning stratocumulus, shallow cumulus, and deep convection. The model is also run towards statistical equilibrium with climatologically relevant large-scale forcings. These model tests are validated against large-eddy simulation (LES) with the same forcings. The comparison of results verifies the capacity of this closure to realistically represent different cloud and convective processes. Implementation of the closure in an idealized GCM allows us to study cloud feedbacks to climate change and to study the interactions between clouds, convections, and the large-scale circulation.

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 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.

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.

Modeling short-term concentration fluctuations of semi-volatile pollutants in the soil-plant-atmosphere system.

PubMed

Bao, Zhongwen; Haberer, Christina M; Maier, Uli; Beckingham, Barbara; Amos, Richard T; Grathwohl, Peter

2016-11-01

Temperature changes can drive cycling of semi-volatile pollutants between different environmental compartments (e.g. atmosphere, soil, plants). To evaluate the impact of daily temperature changes on atmospheric concentration fluctuations we employed a physically based model coupling soil, plants and the atmosphere, which accounts for heat transport, effective gas diffusion, sorption and biodegradation in the soil as well as eddy diffusion and photochemical oxidation in the atmospheric boundary layer of varying heights. The model results suggest that temperature-driven re-volatilization and uptake in soils cannot fully explain significant diurnal concentration fluctuations of atmospheric pollutants as for example observed for polychlorinated biphenyls (PCBs). This holds even for relatively low water contents (high gas diffusivity) and high sorption capacity of the topsoil (high organic carbon content and high pollutant concentration in the topsoil). Observed concentration fluctuations, however, can be easily matched if a rapidly-exchanging environmental compartment, such as a plant layer, is introduced. At elevated temperatures, plants release organic pollutants, which are rapidly distributed in the atmosphere by eddy diffusion. For photosensitive compounds, e.g. some polycyclic aromatic hydrocarbons (PAHs), decreasing atmospheric concentrations would be expected during daytime for the bare soil scenario. This decline is buffered by a plant layer, which acts as a ground-level reservoir. The modeling results emphasize the importance of a rapidly-exchanging compartment above ground to explain short-term atmospheric concentration fluctuations. Copyright © 2016 Elsevier B.V. All rights reserved.

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.

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.

Dynamically consistent parameterization of mesoscale eddies. Part III: Deterministic approach

NASA Astrophysics Data System (ADS)

Berloff, Pavel

2018-07-01

This work continues development of dynamically consistent parameterizations for representing mesoscale eddy effects in non-eddy-resolving and eddy-permitting ocean circulation models and focuses on the classical double-gyre problem, in which the main dynamic eddy effects maintain eastward jet extension of the western boundary currents and its adjacent recirculation zones via eddy backscatter mechanism. Despite its fundamental importance, this mechanism remains poorly understood, and in this paper we, first, study it and, then, propose and test its novel parameterization. We start by decomposing the reference eddy-resolving flow solution into the large-scale and eddy components defined by spatial filtering, rather than by the Reynolds decomposition. Next, we find that the eastward jet and its recirculations are robustly present not only in the large-scale flow itself, but also in the rectified time-mean eddies, and in the transient rectified eddy component, which consists of highly anisotropic ribbons of the opposite-sign potential vorticity anomalies straddling the instantaneous eastward jet core and being responsible for its continuous amplification. The transient rectified component is separated from the flow by a novel remapping method. We hypothesize that the above three components of the eastward jet are ultimately driven by the small-scale transient eddy forcing via the eddy backscatter mechanism, rather than by the mean eddy forcing and large-scale nonlinearities. We verify this hypothesis by progressively turning down the backscatter and observing the induced flow anomalies. The backscatter analysis leads us to formulating the key eddy parameterization hypothesis: in an eddy-permitting model at least partially resolved eddy backscatter can be significantly amplified to improve the flow solution. Such amplification is a simple and novel eddy parameterization framework implemented here in terms of local, deterministic flow roughening controlled by single parameter. We test the parameterization skills in an hierarchy of non-eddy-resolving and eddy-permitting modifications of the original model and demonstrate, that indeed it can be highly efficient for restoring the eastward jet extension and its adjacent recirculation zones. The new deterministic parameterization framework not only combines remarkable simplicity with good performance but also is dynamically transparent, therefore, it provides a powerful alternative to the common eddy diffusion and emerging stochastic parameterizations.

Temperature boundary layer profiles in turbulent Rayleigh-Benard convection

NASA Astrophysics Data System (ADS)

Ching, Emily S. C.; Emran, Mohammad S.; Horn, Susanne; Shishkina, Olga

2017-11-01

Classical boundary-layer theory for steady flows cannot adequately describe the boundary layer profiles in turbulent Rayleigh-Benard convection. We have developed a thermal boundary layer equation which takes into account fluctuations in terms of an eddy thermal diffusivity. Based on Prandtl's mixing length ideas, we relate the eddy thermal diffusivity to the stream function. With this proposed relation, we can solve the thermal boundary layer equation and obtain a closed-form expression for the dimensionless mean temperature profile in terms of two independent parameters: θ(ξ) =1/b∫0b ξ [ 1 +3a3/b3(η - arctan(η)) ] - c dη , where ξ is the similarity variable and the parameters a, b, and c are related by the condition θ(∞) = 1 . With a proper choice of the parameters, our predictions of the temperature profile are in excellent agreement with the results of our direct numerical simulations for a wide range of Prandtl numbers (Pr), from Pr=0.01 to Pr=2547.9. OS, ME and SH acknowledge the financial support by the Deutsche Forschungsgemeinschaft (DFG) under Grants Sh405/4-2 (Heisenberg fellowship), Sh405/3-2 and Ho 5890/1-1, respectively.

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.

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.

Multi-scale observations of the variability of magmatic CO2 emissions, Mammoth Mountain, CA, USA

NASA Astrophysics Data System (ADS)

Lewicki, J. L.; Hilley, G. E.

2014-09-01

One of the primary indicators of volcanic unrest at Mammoth Mountain is diffuse emission of magmatic CO2, which can effectively track this unrest if its variability in space and time and relationship to near-surface meteorological and hydrologic phenomena versus those occurring at depth beneath the mountain are understood. In June-October 2013, we conducted accumulation chamber soil CO2 flux surveys and made half-hourly CO2 flux measurements with automated eddy covariance and accumulation chamber (auto-chamber) instrumentation at the largest area of diffuse CO2 degassing on Mammoth Mountain (Horseshoe Lake tree kill; HLTK). Estimated CO2 emission rates for HLTK based on 20 June, 30 July, and 24-25 October soil CO2 flux surveys were 165, 172, and 231 t d- 1, respectively. The average (June-October) CO2 emission rate estimated for this area was 123 t d- 1 based on an inversion of 4527 eddy covariance CO2 flux measurements and corresponding modeled source weight functions. Average daily eddy covariance and auto-chamber CO2 fluxes consistently declined over the four-month observation time. Wavelet analysis of auto-chamber CO2 flux and environmental parameter time series was used to evaluate the periodicity of, and local correlation between these variables in time-frequency space. Overall, CO2 emissions at HLTK were highly dynamic, displaying short-term (hourly to weekly) temporal variability related to meteorological and hydrologic changes, as well as long-term (monthly to multi-year) variations related to migration of CO2-rich magmatic fluids beneath the volcano. Accumulation chamber soil CO2 flux surveys were also conducted in the four additional areas of diffuse CO2 degassing on Mammoth Mountain in July-August 2013. Summing CO2 emission rates for all five areas yielded a total for the mountain of 311 t d- 1, which may suggest that emissions returned to 1998-2009 levels, following an increase from 2009 to 2011.

Modeling thermospheric neutral density

NASA Astrophysics Data System (ADS)

Qian, Liying

Satellite drag prediction requires determination of thermospheric neutral density. The NCAR Thermosphere-Ionosphere-Electrodynamics General Circulation Model (TIEGCM) and the global-mean Thermosphere-Ionosphere-Mesosphere-Electrodynamics General Circulation Model (TIMEGCM) were used to quantify thermospheric neutral density and its variations, focusing on annual/semiannual variation, the effect of using measured solar irradiance on model calculations of solar-cycle variation, and global change in the thermosphere. Satellite drag data and the MSIS00 empirical model were utilized to compare to the TIEGCM simulations. The TIEGCM simulations indicated that eddy diffusion and its annual/semiannual variation is a mechanism for annual/semiannual density variation in the thermosphere. It was found that eddy diffusion near the turbopause can effectively influence thermospheric neutral density. Eddy diffusion, together with annual insolation variation and large-scale circulation, generated global annual/semiannual density variation observed by satellite drag. Using measured solar irradiance as solar input for the TIEGCM improved the solar-cycle dependency of the density calculation shown in F10.7 -based thermospheric empirical models. It has been found that the empirical models overestimate density at low solar activity. The TIEGCM simulations did not show such solar-cycle dependency. Using historic measurements of CO2 and F 10.7, simulations of the global-mean TIMEGCM showed that thermospheric neutral density at 400 km had an average long-term decrease of 1.7% per decade from 1970 to 2000. A forecast of density decrease for solar cycle 24 suggested that thermospheric density will decrease at 400 km from present to the end of solar cycle 24 at a rate of 2.7% per decade. Reduction in thermospheric density causes less atmospheric drag on earth-orbiting space objects. The implication of this long-term decrease of thermospheric neutral density is that it will increase the lifetime of satellites, but also it will increase the amount of space junk.

Multi-scale observations of the variability of magmatic CO2 emissions, Mammoth Mountain, CA, USA

USGS Publications Warehouse

Lewicki, Jennifer L.; Hilley, George E.

2014-01-01

One of the primary indicators of volcanic unrest at Mammoth Mountain is diffuse emission of magmatic CO2, which can effectively track this unrest if its variability in space and time and relationship to near-surface meteorological and hydrologic phenomena versus those occurring at depth beneath the mountain are understood. In June–October 2013, we conducted accumulation chamber soil CO2 flux surveys and made half-hourly CO2 flux measurements with automated eddy covariance and accumulation chamber (auto-chamber) instrumentation at the largest area of diffuse CO2 degassing on Mammoth Mountain (Horseshoe Lake tree kill; HLTK). Estimated CO2 emission rates for HLTK based on 20 June, 30 July, and 24–25 October soil CO2 flux surveys were 165, 172, and 231 t d− 1, respectively. The average (June–October) CO2 emission rate estimated for this area was 123 t d− 1 based on an inversion of 4527 eddy covariance CO2 flux measurements and corresponding modeled source weight functions. Average daily eddy covariance and auto-chamber CO2 fluxes consistently declined over the four-month observation time. Wavelet analysis of auto-chamber CO2 flux and environmental parameter time series was used to evaluate the periodicity of, and local correlation between these variables in time–frequency space. Overall, CO2 emissions at HLTK were highly dynamic, displaying short-term (hourly to weekly) temporal variability related to meteorological and hydrologic changes, as well as long-term (monthly to multi-year) variations related to migration of CO2-rich magmatic fluids beneath the volcano. Accumulation chamber soil CO2 flux surveys were also conducted in the four additional areas of diffuse CO2 degassing on Mammoth Mountain in July–August 2013. Summing CO2 emission rates for all five areas yielded a total for the mountain of 311 t d− 1, which may suggest that emissions returned to 1998–2009 levels, following an increase from 2009 to 2011.

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.

Maximum entropy production principle for geostrophic turbulence

NASA Astrophysics Data System (ADS)

Sommeria, J.; Bouchet, F.; Chavanis, P. H.

2003-04-01

In 2D turbulence, complex stirring leads to the formation of steady organized states, once fine scale fluctuations have been filtered out. This self-organization can be explained in terms of statistical equilibrium for vorticity, as the most likely outcome of vorticity parcel rearrangements with the constraints of the conservation laws. A mixing entropy describing the vorticity rearrangements is introduced. Extension to the shallow water system has been proposed by Chavanis P.H. and Sommeria J. (2002), Phys. Rev. E. Generalization to multi-layer geostrophic flows is formally straightforward. Outside equilibrium, eddy fluxes should drive the system toward equilibrium, in the spirit of non equilibrium linear thermodynamics. This can been formalized in terms of a principle of maximum entropy production (MEP), as shown by Robert and Sommeria (1991), Phys. Rev. Lett. 69. Then a parameterization of eddy fluxes is obtained, involving an eddy diffusivity plus a drift term acting at larger scale. These two terms balance each other at equilibrium, resulting in a non trivial steady flow, which is the mean state of the statistical equilibrium. Applications of this eddy parametrization will be presented, in the context of oceanic circulation and Jupiter's Great Red Spot. Quantitative tests will be discussed, obtained by comparisons with direct numerical simulations. Kinetic models, inspired from plasma physics, provide a more precise description of the relaxation toward equilibrium, as shown by Chavanis P.H. 2000 ``Quasilinear theory of the 2D Euler equation'', Phys. Rev. Lett. 84. This approach provides relaxation equations with a form similar to the MEP, but not identical. In conclusion, the MEP provides the right trends of the system but its precise justification remains elusive.

A Multi-wavenumber Theory for Eddy Diffusivities: Applications to the DIMES Region

NASA Astrophysics Data System (ADS)

Chen, R.; Gille, S. T.; McClean, J.; Flierl, G.; Griesel, A.

2014-12-01

Climate models are sensitive to the representation of ocean mixing processes. This has motivated recent efforts to collect observations aimed at improving mixing estimates and parameterizations. The US/UK field program Diapycnal and Isopycnal Mixing Experiment in the Southern Ocean (DIMES), begun in 2009, is providing such estimates upstream of and within the Drake Passage. This region is characterized by topography, and strong zonal jets. In previous studies, mixing length theories, based on the assumption that eddies are dominated by a single wavenumber and phase speed, were formulated to represent the estimated mixing patterns in jets. However, in spite of the success of the single wavenumber theory in some other scenarios, it does not effectively predict the vertical structures of observed eddy diffusivities in the DIMES area. Considering that eddy motions encompass a wide range of wavenumbers, which all contribute to mixing, in this study we formulated a multi-wavenumber theory to predict eddy mixing rates. We test our theory for a domain encompassing the entire Southern Ocean. We estimated eddy diffusivities and mixing lengths from one million numerical floats in a global eddying model. These float-based mixing estimates were compared with the predictions from both the single-wavenumber and the multi-wavenumber theories. Our preliminary results in the DIMES area indicate that, compared to the single-wavenumber theory, the multi-wavenumber theory better predicts the vertical mixing structures in the vast areas where the mean flow is weak; however in the intense jet region, both theories have similar predictive skill.

A quantitative comparison of two methods to correct eddy current-induced distortions in DT-MRI.

PubMed

Muñoz Maniega, Susana; Bastin, Mark E; Armitage, Paul A

2007-04-01

Eddy current-induced geometric distortions of single-shot, diffusion-weighted, echo-planar (DW-EP) images are a major confounding factor to the accurate determination of water diffusion parameters in diffusion tensor MRI (DT-MRI). Previously, it has been suggested that these geometric distortions can be removed from brain DW-EP images using affine transformations determined from phantom calibration experiments using iterative cross-correlation (ICC). Since this approach was first described, a number of image-based registration methods have become available that can also correct eddy current-induced distortions in DW-EP images. However, as yet no study has investigated whether separate eddy current calibration or image-based registration provides the most accurate way of removing these artefacts from DT-MRI data. Here we compare how ICC phantom calibration and affine FLIRT (http://www.fmrib.ox.ac.uk), a popular image-based multi-modal registration method that can correct both eddy current-induced distortions and bulk subject motion, perform when registering DW-EP images acquired with different slice thicknesses (2.8 and 5 mm) and b-values (1000 and 3000 s/mm(2)). With the use of consistency testing, it was found that ICC was a more robust algorithm for correcting eddy current-induced distortions than affine FLIRT, especially at high b-value and small slice thickness. In addition, principal component analysis demonstrated that the combination of ICC phantom calibration (to remove eddy current-induced distortions) with rigid body FLIRT (to remove bulk subject motion) provided a more accurate registration of DT-MRI data than that achieved by affine FLIRT.

Plasma transport in the Io torus - The importance of microscopic diffusion

NASA Technical Reports Server (NTRS)

Mei, YI; Thorne, Richard M.

1991-01-01

This paper considers the question of whether the distribution of mass in the Io plasma torus is consistent with the concept of interchange eddy transport. Specifically, the flux tube content exhibits a gradual decrease with increasing radial distance from the source near Io without any evidence for substantial density irregularity associated with the plasma source or loss. Using a simple one-dimensional numerical model to simulate macroscopic interchange eddy transport, it is demonstrated that this smooth equilibrium distribution of mass can occur but only with the inclusion of a minimal level of small scale microscopic mixing at a rate approaching Bohm diffusion. Otherwise, the system exhibits a chaotic appearance which never approaches an equilibrium distribution. Various physical mechanisms for the microscopic diffusion process which is required to provide a sufficiently rapid mixing of material between the macroscopic eddies are discussed.

Correction for Eddy Current-Induced Echo-Shifting Effect in Partial-Fourier Diffusion Tensor Imaging.

PubMed

Truong, Trong-Kha; Song, Allen W; Chen, Nan-Kuei

2015-01-01

In most diffusion tensor imaging (DTI) studies, images are acquired with either a partial-Fourier or a parallel partial-Fourier echo-planar imaging (EPI) sequence, in order to shorten the echo time and increase the signal-to-noise ratio (SNR). However, eddy currents induced by the diffusion-sensitizing gradients can often lead to a shift of the echo in k-space, resulting in three distinct types of artifacts in partial-Fourier DTI. Here, we present an improved DTI acquisition and reconstruction scheme, capable of generating high-quality and high-SNR DTI data without eddy current-induced artifacts. This new scheme consists of three components, respectively, addressing the three distinct types of artifacts. First, a k-space energy-anchored DTI sequence is designed to recover eddy current-induced signal loss (i.e., Type 1 artifact). Second, a multischeme partial-Fourier reconstruction is used to eliminate artificial signal elevation (i.e., Type 2 artifact) associated with the conventional partial-Fourier reconstruction. Third, a signal intensity correction is applied to remove artificial signal modulations due to eddy current-induced erroneous T2(∗) -weighting (i.e., Type 3 artifact). These systematic improvements will greatly increase the consistency and accuracy of DTI measurements, expanding the utility of DTI in translational applications where quantitative robustness is much needed.

Correction for Eddy Current-Induced Echo-Shifting Effect in Partial-Fourier Diffusion Tensor Imaging

PubMed Central

Truong, Trong-Kha; Song, Allen W.; Chen, Nan-kuei

2015-01-01

In most diffusion tensor imaging (DTI) studies, images are acquired with either a partial-Fourier or a parallel partial-Fourier echo-planar imaging (EPI) sequence, in order to shorten the echo time and increase the signal-to-noise ratio (SNR). However, eddy currents induced by the diffusion-sensitizing gradients can often lead to a shift of the echo in k-space, resulting in three distinct types of artifacts in partial-Fourier DTI. Here, we present an improved DTI acquisition and reconstruction scheme, capable of generating high-quality and high-SNR DTI data without eddy current-induced artifacts. This new scheme consists of three components, respectively, addressing the three distinct types of artifacts. First, a k-space energy-anchored DTI sequence is designed to recover eddy current-induced signal loss (i.e., Type 1 artifact). Second, a multischeme partial-Fourier reconstruction is used to eliminate artificial signal elevation (i.e., Type 2 artifact) associated with the conventional partial-Fourier reconstruction. Third, a signal intensity correction is applied to remove artificial signal modulations due to eddy current-induced erroneous T 2 ∗-weighting (i.e., Type 3 artifact). These systematic improvements will greatly increase the consistency and accuracy of DTI measurements, expanding the utility of DTI in translational applications where quantitative robustness is much needed. PMID:26413505

Multi-Scale Modeling and the Eddy-Diffusivity/Mass-Flux (EDMF) Parameterization

NASA Astrophysics Data System (ADS)

Teixeira, J.

2015-12-01

Turbulence and convection play a fundamental role in many key weather and climate science topics. Unfortunately, current atmospheric models cannot explicitly resolve most turbulent and convective flow. Because of this fact, turbulence and convection in the atmosphere has to be parameterized - i.e. equations describing the dynamical evolution of the statistical properties of turbulence and convection motions have to be devised. Recently a variety of different models have been developed that attempt at simulating the atmosphere using variable resolution. A key problem however is that parameterizations are in general not explicitly aware of the resolution - the scale awareness problem. In this context, we will present and discuss a specific approach, the Eddy-Diffusivity/Mass-Flux (EDMF) parameterization, that not only is in itself a multi-scale parameterization but it is also particularly well suited to deal with the scale-awareness problems that plague current variable-resolution models. It does so by representing small-scale turbulence using a classic Eddy-Diffusivity (ED) method, and the larger-scale (boundary layer and tropospheric-scale) eddies as a variety of plumes using the Mass-Flux (MF) concept.

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.

Critical role for mesoscale eddy diffusion in supplying oxygen to hypoxic ocean waters

NASA Astrophysics Data System (ADS)

Gnanadesikan, Anand; Bianchi, Daniele; Pradal, Marie-Aude

2013-10-01

of the oceanic lateral eddy diffusion coefficient Aredi vary by more than an order of magnitude, ranging from less than a few hundred m2/s to thousands of m2/s. This uncertainty has first-order implications for the intensity of oceanic hypoxia, which is poorly simulated by the current generation of Earth System Models. Using satellite-based estimate of oxygen consumption in hypoxic waters to estimate the required diffusion coefficient for these waters gives a value of order 1000 m2/s. Varying Aredi across a suite of Earth System Models yields a broadly consistent result given a thermocline diapycnal diffusion coefficient of 1 × 10-5 m2/s.

Large Eddy Simulation of a Supercritical Turbulent Mixing Layer

NASA Astrophysics Data System (ADS)

Sheikhi, Reza; Hadi, Fatemeh; Safari, Mehdi

2017-11-01

Supercritical turbulent flows are relevant to a wide range of applications such as supercritical power cycles, gas turbine combustors, rocket propulsion and internal combustion engines. Large eddy simulation (LES) analysis of such flows involves solving mass, momentum, energy and scalar transport equations with inclusion of generalized diffusion fluxes. These equations are combined with a real gas equation of state and the corresponding thermodynamic mixture variables. Subgrid scale models are needed for not only the conventional convective terms but also the additional high pressure effects arising due to the nonlinearity associated with generalized diffusion fluxes and real gas equation of state. In this study, LES is carried out to study the high pressure turbulent mixing of methane with carbon dioxide in a temporally developing mixing layer under supercritical condition. LES results are assessed by comparing with data obtained from direct numerical simulation (DNS) of the same layer. LES predictions agree favorably with DNS data and represent several key supercritical turbulent flow features such as high density gradient regions. Supported by DOE Grant SC0017097; computational support is provided by DOE National Energy Research Scientific Computing Center.

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.

On the Effective Thermal Conductivity of Frost Considering Mass Diffusion and Eddy Convection

NASA Technical Reports Server (NTRS)

Kandula, Max

2010-01-01

A physical model for the effective thermal conductivity of water frost is proposed for application to the full range of frost density. The proposed model builds on the Zehner-Schlunder one-dimensional formulation for porous media appropriate for solid-to-fluid thermal conductivity ratios less than about 1000. By superposing the effects of mass diffusion and eddy convection on stagnant conduction in the fluid, the total effective thermal conductivity of frost is shown to be satisfactorily described. It is shown that the effects of vapor diffusion and eddy convection on the frost conductivity are of the same order. The results also point out that idealization of the frost structure by cylindrical inclusions offers a better representation of the effective conductivity of frost as compared to spherical inclusions. Satisfactory agreement between the theory and the measurements for the effective thermal conductivity of frost is demonstrated for a wide range of frost density and frost temperature.

The dependence of the oceans MOC on mesoscale eddy diffusivities: A model study

NASA Technical Reports Server (NTRS)

Marshall, John; Scott, Jeffery R.; Romanou, Anastasia; Kelley, Maxwell; Leboissetier, Anthony

2017-01-01

The dependence of the depth and strength of the ocean's global meridional overturning cells (MOC) on the specification of mesoscale eddy diffusivity (K) is explored in two ocean models. The GISS and MIT ocean models are driven by the same prescribed forcing fields, configured in similar ways, spun up to equilibrium for a range of K 's and the resulting MOCs mapped and documented. Scaling laws implicit in modern theories of the MOC are used to rationalize the results. In all calculations the K used in the computation of eddy-induced circulation and that used in the representation of eddy stirring along neutral surfaces, is set to the same value but is changed across experiments. We are able to connect changes in the strength and depth of the Atlantic MOC, the southern ocean upwelling MOC, and the deep cell emanating from Antarctica, to changes in K.

Understanding the Effects of Lower Boundary Conditions and Eddy Diffusion on the Ionosphere-Thermosphere System

NASA Astrophysics Data System (ADS)

Malhotra, G.; Ridley, A. J.; Marsh, D. R.; Wu, C.; Paxton, L. J.

2017-12-01

The exchange of energy between lower atmospheric regions with the ionosphere-thermosphere (IT) system is not well understood. A number of studies have observed day-to-day and seasonal variabilities in the difference between data and model output of various IT parameters. It is widely speculated that the forcing from the lower atmosphere, variability in weather systems and gravity waves that propagate upward from troposphere into the upper mesosphere and lower thermosphere (MLT) may be responsible for these spatial and temporal variations in the IT region, but their exact nature is unknown. These variabilities can be interpreted in two ways: variations in state (density, temperature, wind) of the upper mesosphere or spatial and temporal changes in the small-scale mixing, or Eddy diffusion that is parameterized within the model.In this study, firstly, we analyze the sensitivity of the thermospheric and ionospheric states - neutral densities, O/N2, total electron content (TEC), peak electron density, and peak electron height - to various lower boundary conditions in the Global Ionosphere Thermosphere Model (GITM). We use WACCM-X and GSWM to drive the lower atmospheric boundary in GITM at 100 km, and compare the results with the current MSIS-driven version of GITM, analyzing which of these simulations match the measurements from GOCE, GUVI, CHAMP, and GPS-derived TEC best. Secondly, we analyze the effect of eddy diffusion in the IT system. The turbulence due to eddy mixing cannot be directly measured and it is a challenge to completely characterize its linear and non-linear effects from other influences, since the eddy diffusion both influences the composition through direct mixing and the temperature structure due to turbulent conduction changes. In this study we input latitudinal and seasonal profiles of eddy diffusion into GITM and then analyze the changes in the thermospheric and ionospheric parameters. These profiles will be derived from both WACC-X simulations and direct observations of errors between the model and data such as GUVI O/N2 ratios and TEC data. In each case, the model results will be compared to data to determine the improvement.

Modeling mesoscale eddies

NASA Astrophysics Data System (ADS)

Canuto, V. M.; Dubovikov, M. S.

Mesoscale eddies are not resolved in coarse resolution ocean models and must be modeled. They affect both mean momentum and scalars. At present, no generally accepted model exists for the former; in the latter case, mesoscales are modeled with a bolus velocity u∗ to represent a sink of mean potential energy. However, comparison of u∗(model) vs. u∗ (eddy resolving code, [J. Phys. Ocean. 29 (1999) 2442]) has shown that u∗(model) is incomplete and that additional terms, "unrelated to thickness source or sinks", are required. Thus far, no form of the additional terms has been suggested. To describe mesoscale eddies, we employ the Navier-Stokes and scalar equations and a turbulence model to treat the non-linear interactions. We then show that the problem reduces to an eigenvalue problem for the mesoscale Bernoulli potential. The solution, which we derive in analytic form, is used to construct the momentum and thickness fluxes. In the latter case, the bolus velocity u∗ is found to contain two types of terms: the first type entails the gradient of the mean potential vorticity and represents a positive contribution to the production of mesoscale potential energy; the second type of terms, which is new, entails the velocity of the mean flow and represents a negative contribution to the production of mesoscale potential energy, or equivalently, a backscatter process whereby a fraction of the mesoscale potential energy is returned to the original reservoir of mean potential energy. This type of terms satisfies the physical description of the additional terms given by [J. Phys. Ocean. 29 (1999) 2442]. The mesoscale flux that enters the momentum equations is also contributed by two types of terms of the same physical nature as those entering the thickness flux. The potential vorticity flux is also shown to contain two types of terms: the first is of the gradient-type while the other terms entail the velocity of the mean flow. An expression is derived for the mesoscale diffusivity κM and for the mesoscale kinetic energy K in terms of the large-scale fields. The predicted κM( z) agrees with that of heuristic models. The complete mesoscale model in isopycnal coordinates is presented in Appendix D and can be used in coarse resolution ocean global circulation models.

Demonstration of Nonlinearity Bias in the Measurement of the Apparent Diffusion Coefficient in Multicenter Trials

PubMed Central

Malyarenko, Dariya; Newitt, David; Wilmes, Lisa; Tudorica, Alina; Helmer, Karl G.; Arlinghaus, Lori R.; Jacobs, Michael A.; Jajamovich, Guido; Taouli, Bachir; Yankeelov, Thomas E.; Huang, Wei; Chenevert, Thomas L.

2015-01-01

Purpose Characterize system-specific bias across common magnetic resonance imaging (MRI) platforms for quantitative diffusion measurements in multicenter trials. Methods Diffusion weighted imaging (DWI) was performed on an ice-water phantom along the superior-inferior (SI) and right-left (RL) orientations spanning ±150 mm. The same scanning protocol was implemented on 14 MRI systems at seven imaging centers. The bias was estimated as a deviation of measured from known apparent diffusion coefficient (ADC) along individual DWI directions. The relative contributions of gradient nonlinearity, shim errors, imaging gradients and eddy currents were assessed independently. The observed bias errors were compared to numerical models. Results The measured systematic ADC errors scaled quadratically with offset from isocenter, and ranged between −55% (SI) and 25% (RL). Nonlinearity bias was dependent on system design and diffusion gradient direction. Consistent with numerical models, minor ADC errors (±5%) due to shim, imaging and eddy currents were mitigated by double echo DWI and image co-registration of individual gradient directions. Conclusion The analysis confirms gradient nonlinearity as a major source of spatial DW bias and variability in off-center ADC measurements across MRI platforms, with minor contributions from shim, imaging gradients and eddy currents. The developed protocol enables empiric description of systematic bias in multicenter quantitative DWI studies. PMID:25940607

Demonstration of nonlinearity bias in the measurement of the apparent diffusion coefficient in multicenter trials.

PubMed

Malyarenko, Dariya I; Newitt, David; J Wilmes, Lisa; Tudorica, Alina; Helmer, Karl G; Arlinghaus, Lori R; Jacobs, Michael A; Jajamovich, Guido; Taouli, Bachir; Yankeelov, Thomas E; Huang, Wei; Chenevert, Thomas L

2016-03-01

Characterize system-specific bias across common magnetic resonance imaging (MRI) platforms for quantitative diffusion measurements in multicenter trials. Diffusion weighted imaging (DWI) was performed on an ice-water phantom along the superior-inferior (SI) and right-left (RL) orientations spanning ± 150 mm. The same scanning protocol was implemented on 14 MRI systems at seven imaging centers. The bias was estimated as a deviation of measured from known apparent diffusion coefficient (ADC) along individual DWI directions. The relative contributions of gradient nonlinearity, shim errors, imaging gradients, and eddy currents were assessed independently. The observed bias errors were compared with numerical models. The measured systematic ADC errors scaled quadratically with offset from isocenter, and ranged between -55% (SI) and 25% (RL). Nonlinearity bias was dependent on system design and diffusion gradient direction. Consistent with numerical models, minor ADC errors (± 5%) due to shim, imaging and eddy currents were mitigated by double echo DWI and image coregistration of individual gradient directions. The analysis confirms gradient nonlinearity as a major source of spatial DW bias and variability in off-center ADC measurements across MRI platforms, with minor contributions from shim, imaging gradients and eddy currents. The developed protocol enables empiric description of systematic bias in multicenter quantitative DWI studies. © 2015 Wiley Periodicals, Inc.

Long-term Trends and Variability of Eddy Activities in the South China Sea

NASA Astrophysics Data System (ADS)

Zhang, M.; von Storch, H.

2017-12-01

For constructing empirical downscaling models and projecting possible future states of eddy activities in the South China Sea (SCS), long-term statistical characteristics of the SCS eddy are needed. We use a daily global eddy-resolving model product named STORM covering the period of 1950-2010. This simulation has employed the MPI-OM model with a mean horizontal resolution of 10km and been driven by the NCEP reanalysis-1 data set. An eddy detection and tracking algorithm operating on the gridded sea surface height anomaly (SSHA) fields was developed. A set of parameters for the criteria in the SCS are determined through sensitivity tests. Our method detected more than 6000 eddy tracks in the South China Sea. For all of them, eddy diameters, track length, eddy intensity, eddy lifetime and eddy frequency were determined. The long-term trends and variability of those properties also has been derived. Most of the eddies propagate westward. Nearly 100 eddies travel longer than 1000km, and over 800 eddies have a lifespan of more than 2 months. Furthermore, for building the statistical empirical model, the relationship between the SCS eddy statistics and the large-scale atmospheric and oceanic phenomena has been investigated.

Diagnosing isopycnal diffusivity in an eddying, idealized midlatitude ocean basin via Lagrangian, in Situ, Global, High-Performance Particle Tracking (LIGHT)

DOE PAGES

Wolfram, Phillip J.; Ringler, Todd D.; Maltrud, Mathew E.; ...

2015-08-01

Isopycnal diffusivity due to stirring by mesoscale eddies in an idealized, wind-forced, eddying, midlatitude ocean basin is computed using Lagrangian, in Situ, Global, High-Performance Particle Tracking (LIGHT). Simulation is performed via LIGHT within the Model for Prediction across Scales Ocean (MPAS-O). Simulations are performed at 4-, 8-, 16-, and 32-km resolution, where the first Rossby radius of deformation (RRD) is approximately 30 km. Scalar and tensor diffusivities are estimated at each resolution based on 30 ensemble members using particle cluster statistics. Each ensemble member is composed of 303 665 particles distributed across five potential density surfaces. Diffusivity dependence upon modelmore » resolution, velocity spatial scale, and buoyancy surface is quantified and compared with mixing length theory. The spatial structure of diffusivity ranges over approximately two orders of magnitude with values of O(10 5) m 2 s –1 in the region of western boundary current separation to O(10 3) m 2 s –1 in the eastern region of the basin. Dominant mixing occurs at scales twice the size of the first RRD. Model resolution at scales finer than the RRD is necessary to obtain sufficient model fidelity at scales between one and four RRD to accurately represent mixing. Mixing length scaling with eddy kinetic energy and the Lagrangian time scale yield mixing efficiencies that typically range between 0.4 and 0.8. In conclusion, a reduced mixing length in the eastern region of the domain relative to the west suggests there are different mixing regimes outside the baroclinic jet region.« less

Stochastic Convection Parameterizations: The Eddy-Diffusivity/Mass-Flux (EDMF) Approach (Invited)

NASA Astrophysics Data System (ADS)

Teixeira, J.

2013-12-01

In this presentation it is argued that moist convection parameterizations need to be stochastic in order to be realistic - even in deterministic atmospheric prediction systems. A new unified convection and boundary layer parameterization (EDMF) that optimally combines the Eddy-Diffusivity (ED) approach for smaller-scale boundary layer mixing with the Mass-Flux (MF) approach for larger-scale plumes is discussed. It is argued that for realistic simulations stochastic methods have to be employed in this new unified EDMF. Positive results from the implementation of the EDMF approach in atmospheric models are presented.

Computing eddy-driven effective diffusivity using Lagrangian particles

DOE PAGES

Wolfram, Phillip J.; Ringler, Todd D.

2017-08-14

A novel method to derive effective diffusivity from Lagrangian particle trajectory data sets is developed and then analyzed relative to particle-derived meridional diffusivity for eddy-driven mixing in an idealized circumpolar current. Quantitative standard dispersion- and transport-based mixing diagnostics are defined, compared and contrasted to motivate the computation and use of effective diffusivity derived from Lagrangian particles. We compute the effective diffusivity by first performing scalar transport on Lagrangian control areas using stored trajectories computed from online Lagrangian In-situ Global High-performance particle Tracking (LIGHT) using the Model for Prediction Across Scales Ocean (MPAS-O). Furthermore, the Lagrangian scalar transport scheme is comparedmore » against an Eulerian scalar transport scheme. Spatially-variable effective diffusivities are computed from resulting time-varying cumulative concentrations that vary as a function of cumulative area. The transport-based Eulerian and Lagrangian effective diffusivity diagnostics are found to be qualitatively consistent with the dispersion-based diffusivity. All diffusivity estimates show a region of increased subsurface diffusivity within the core of an idealized circumpolar current and results are within a factor of two of each other. The Eulerian and Lagrangian effective diffusivities are most similar; smaller and more spatially diffused values are obtained with the dispersion-based diffusivity computed with particle clusters.« less

Computing eddy-driven effective diffusivity using Lagrangian particles

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

Wolfram, Phillip J.; Ringler, Todd D.

A novel method to derive effective diffusivity from Lagrangian particle trajectory data sets is developed and then analyzed relative to particle-derived meridional diffusivity for eddy-driven mixing in an idealized circumpolar current. Quantitative standard dispersion- and transport-based mixing diagnostics are defined, compared and contrasted to motivate the computation and use of effective diffusivity derived from Lagrangian particles. We compute the effective diffusivity by first performing scalar transport on Lagrangian control areas using stored trajectories computed from online Lagrangian In-situ Global High-performance particle Tracking (LIGHT) using the Model for Prediction Across Scales Ocean (MPAS-O). Furthermore, the Lagrangian scalar transport scheme is comparedmore » against an Eulerian scalar transport scheme. Spatially-variable effective diffusivities are computed from resulting time-varying cumulative concentrations that vary as a function of cumulative area. The transport-based Eulerian and Lagrangian effective diffusivity diagnostics are found to be qualitatively consistent with the dispersion-based diffusivity. All diffusivity estimates show a region of increased subsurface diffusivity within the core of an idealized circumpolar current and results are within a factor of two of each other. The Eulerian and Lagrangian effective diffusivities are most similar; smaller and more spatially diffused values are obtained with the dispersion-based diffusivity computed with particle clusters.« less

Eddy Effects in the General Circulation, Spanning Mean Currents, Mesoscale Eddies, and Topographic Generation, Including Submesoscale Nests

DTIC Science & Technology

2013-09-30

bottom form stress (pressure force) and bottom boundary layers – all the aspects associated with turbulent flows over steep topography in the presence of...filaments, and eddies; topographic current separation, form stress , and submesoscale vortex generation; Our work on isoneutral diffusion for tracers...Bump region, are due to the contribution of the bottom stress curl. Fig. 4 shows how the Gulf Stream path is directly linked to the Bottom Pressure

Comparison of GEOS-5 AGCM planetary boundary layer depths computed with various definitions

NASA Astrophysics Data System (ADS)

McGrath-Spangler, E. L.; Molod, A.

2014-07-01

Accurate models of planetary boundary layer (PBL) processes are important for forecasting weather and climate. The present study compares seven methods of calculating PBL depth in the GEOS-5 atmospheric general circulation model (AGCM) over land. These methods depend on the eddy diffusion coefficients, bulk and local Richardson numbers, and the turbulent kinetic energy. The computed PBL depths are aggregated to the Köppen-Geiger climate classes, and some limited comparisons are made using radiosonde profiles. Most methods produce similar midday PBL depths, although in the warm, moist climate classes the bulk Richardson number method gives midday results that are lower than those given by the eddy diffusion coefficient methods. Additional analysis revealed that methods sensitive to turbulence driven by radiative cooling produce greater PBL depths, this effect being most significant during the evening transition. Nocturnal PBLs based on Richardson number methods are generally shallower than eddy diffusion coefficient based estimates. The bulk Richardson number estimate is recommended as the PBL height to inform the choice of the turbulent length scale, based on the similarity to other methods during the day, and the improved nighttime behavior.

Biogenic mixing induced by intermediate Reynolds number swimming in stratified fluids

PubMed Central

Wang, Shiyan; Ardekani, Arezoo M.

2015-01-01

We study fully resolved motion of interacting swimmers in density stratified fluids using an archetypal swimming model called “squirmer”. The intermediate Reynolds number regime is particularly important, because the vast majority of organisms in the aphotic ocean (i.e. regions that are 200 m beneath the sea surface) are small (mm-cm) and their motion is governed by the balance of inertial and viscous forces. Our study shows that the mixing efficiency and the diapycnal eddy diffusivity, a measure of vertical mass flux, within a suspension of squirmers increases with Reynolds number. The mixing efficiency is in the range of O(0.0001–0.04) when the swimming Reynolds number is in the range of O(0.1–100). The values of diapycnal eddy diffusivity and Cox number are two orders of magnitude larger for vertically swimming cells compared to horizontally swimming cells. For a suspension of squirmers in a decaying isotropic turbulence, we find that the diapycnal eddy diffusivity enhances due to the strong viscous dissipation generated by squirmers as well as the interaction of squirmers with the background turbulence. PMID:26628288

Biogenic mixing induced by intermediate Reynolds number swimming in stratified fluids.

PubMed

Wang, Shiyan; Ardekani, Arezoo M

2015-12-02

We study fully resolved motion of interacting swimmers in density stratified fluids using an archetypal swimming model called "squirmer". The intermediate Reynolds number regime is particularly important, because the vast majority of organisms in the aphotic ocean (i.e. regions that are 200 m beneath the sea surface) are small (mm-cm) and their motion is governed by the balance of inertial and viscous forces. Our study shows that the mixing efficiency and the diapycnal eddy diffusivity, a measure of vertical mass flux, within a suspension of squirmers increases with Reynolds number. The mixing efficiency is in the range of O(0.0001-0.04) when the swimming Reynolds number is in the range of O(0.1-100). The values of diapycnal eddy diffusivity and Cox number are two orders of magnitude larger for vertically swimming cells compared to horizontally swimming cells. For a suspension of squirmers in a decaying isotropic turbulence, we find that the diapycnal eddy diffusivity enhances due to the strong viscous dissipation generated by squirmers as well as the interaction of squirmers with the background turbulence.

Comparison of GEOS-5 AGCM Planetary Boundary Layer Depths Computed with Various Definitions

NASA Technical Reports Server (NTRS)

Mcgrath-Spangler, E. L.; Molod, A.

2014-01-01

Accurate models of planetary boundary layer (PBL) processes are important for forecasting weather and climate. The present study compares seven methods of calculating PBL depth in the GEOS-5 atmospheric general circulation model (AGCM) over land. These methods depend on the eddy diffusion coefficients, bulk and local Richardson numbers, and the turbulent kinetic energy. The computed PBL depths are aggregated to the Koppen climate classes, and some limited comparisons are made using radiosonde profiles. Most methods produce similar midday PBL depths, although in the warm, moist climate classes, the bulk Richardson number method gives midday results that are lower than those given by the eddy diffusion coefficient methods. Additional analysis revealed that methods sensitive to turbulence driven by radiative cooling produce greater PBL depths, this effect being most significant during the evening transition. Nocturnal PBLs based on Richardson number are generally shallower than eddy diffusion coefficient based estimates. The bulk Richardson number estimate is recommended as the PBL height to inform the choice of the turbulent length scale, based on the similarity to other methods during the day, and the improved nighttime behavior.

Comparison of GEOS-5 AGCM planetary boundary layer depths computed with various definitions

NASA Astrophysics Data System (ADS)

McGrath-Spangler, E. L.; Molod, A.

2014-03-01

Accurate models of planetary boundary layer (PBL) processes are important for forecasting weather and climate. The present study compares seven methods of calculating PBL depth in the GEOS-5 atmospheric general circulation model (AGCM) over land. These methods depend on the eddy diffusion coefficients, bulk and local Richardson numbers, and the turbulent kinetic energy. The computed PBL depths are aggregated to the Köppen climate classes, and some limited comparisons are made using radiosonde profiles. Most methods produce similar midday PBL depths, although in the warm, moist climate classes, the bulk Richardson number method gives midday results that are lower than those given by the eddy diffusion coefficient methods. Additional analysis revealed that methods sensitive to turbulence driven by radiative cooling produce greater PBL depths, this effect being most significant during the evening transition. Nocturnal PBLs based on Richardson number are generally shallower than eddy diffusion coefficient based estimates. The bulk Richardson number estimate is recommended as the PBL height to inform the choice of the turbulent length scale, based on the similarity to other methods during the day, and the improved nighttime behavior.

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.

Large-Eddy Simulations of Radiatively Driven Convection: Sensitivities to the Representation of Small Scales.

NASA Astrophysics Data System (ADS)

Stevens, Bjorn; Moeng, Chin-Hoh; Sullivan, Peter P.

1999-12-01

Large-eddy simulations of a smoke cloud are examined with respect to their sensitivity to small scales as manifest in either the grid spacing or the subgrid-scale (SGS) model. Calculations based on a Smagorinsky SGS model are found to be more sensitive to the effective resolution of the simulation than are calculations based on the prognostic turbulent kinetic energy (TKE) SGS model. The difference between calculations based on the two SGS models is attributed to the advective transport, diffusive transport, and/or time-rate-of-change terms in the TKE equation. These terms are found to be leading order in the entrainment zone and allow the SGS TKE to behave in a way that tends to compensate for changes that result in larger or smaller resolved scale entrainment fluxes. This compensating behavior of the SGS TKE model is attributed to the fact that changes that reduce the resolved entrainment flux (viz., values of the eddy viscosity in the upper part of the PBL) simultaneously tend to increase the buoyant production of SGS TKE in the radiatively destabilized portion of the smoke cloud. Increased production of SGS TKE in this region then leads to increased amounts of transported, or fossil, SGS TKE in the entrainment zone itself, which in turn leads to compensating increases in the SGS entrainment fluxes. In the Smagorinsky model, the absence of a direct connection between SGS TKE in the entrainment and radiatively destabilized zones prevents this compensating mechanism from being active, and thus leads to calculations whose entrainment rate sensitivities as a whole reflect the sensitivities of the resolved-scale fluxes to values of upper PBL eddy viscosities.

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.

Can a numerically stable subgrid-scale model for turbulent flow computation be ideally accurate?: a preliminary theoretical study for the Gaussian filtered Navier-Stokes equations.

PubMed

Ida, Masato; Taniguchi, Nobuyuki

2003-09-01

This paper introduces a candidate for the origin of the numerical instabilities in large eddy simulation repeatedly observed in academic and practical industrial flow computations. Without resorting to any subgrid-scale modeling, but based on a simple assumption regarding the streamwise component of flow velocity, it is shown theoretically that in a channel-flow computation, the application of the Gaussian filtering to the incompressible Navier-Stokes equations yields a numerically unstable term, a cross-derivative term, which is similar to one appearing in the Gaussian filtered Vlasov equation derived by Klimas [J. Comput. Phys. 68, 202 (1987)] and also to one derived recently by Kobayashi and Shimomura [Phys. Fluids 15, L29 (2003)] from the tensor-diffusivity subgrid-scale term in a dynamic mixed model. The present result predicts that not only the numerical methods and the subgrid-scale models employed but also only the applied filtering process can be a seed of this numerical instability. An investigation concerning the relationship between the turbulent energy scattering and the unstable term shows that the instability of the term does not necessarily represent the backscatter of kinetic energy which has been considered a possible origin of numerical instabilities in large eddy simulation. The present findings raise the question whether a numerically stable subgrid-scale model can be ideally accurate.

Photochemistry of Saturn's Atmosphere. 1; Hydrocarbon Chemistry and Comparisons with ISO Observations

NASA Technical Reports Server (NTRS)

Moses, Julianne I.; Bezard, Bruno; Lellouch, Emmanuel; Gladstone, G. Randall; Feuchtgruber, Helmut; Allen, Mark

2000-01-01

To investigate the details of hydrocarbon photochemistry on Saturn, we have developed a one-dimensional diurnally averaged model that couples hydrocarbon and oxygen photochemistry, molecular and eddy diffusion, radiative transfer, and condensation. The model results are compared with observations from the Infrared Space Observatory (ISO) to place tighter constraints on molecular abundances, to better define Saturn's eddy diffusion coefficient profile, and to identify important chemical schemes that control the abundances of the observable hydrocarbons in Saturn's upper atmosphere. From the ISO observations, we determine that the column 12 densities of CH3, CH3C2H, and C4H2 above 10 mbar are 4 (sup +2) (sub -1.5) x 10 (exp 13) cm (sup -2), (1.1 plus or minus 0.3) x 10 (exp 15) cm (exp -2), and (1.2 plus or minus 0.3) x 10 (exp 14) cm (sup -2), respectively. The observed ISO emission features also indicate C2H2 mixing ratios of 1.2 (sup +0.9) (sub -0.6) x 10 (exp -6) at 0.3 mbar and (2.7 plus or minus 0.8) x 10 (exp -7) at 1.4 mbar, and a C2H6 mixing ratio of (9 plus or minus 2.5) x 10 (exp -6) at 0.5 mbar. Upper limits are provided for C2H4, CH2CCH2, C3H8, and C6H2 sensitivity of the model results to variations in the eddy diffusion coefficient profile, the solar flux, the CH4 photolysis branching ratios, the atomic hydrogen influx, and key reaction rates are discussed in detail. We find that C4H2 and CH3C2H are particularly good tracers of important chemical processes and physical conditions in Saturn's upper atmosphere, and C2H6 is a good tracer of the eddy diffusion coefficient in Saturn's lower stratosphere. The eddy diffusion coefficient must be smaller than approximately 3 x 10 (exp 4) sq cm s (sup -1) at pressures greater than 1 mbar in order to reproduce the C2H6 abundance inferred from ISO observations. The eddy diffusion coefficients in the upper stratosphere could be constrained by observations of CH3 radicals if the low-temperature chemistry of CH3 were better understood. We also discuss the implications of our modeling for aerosol formation in Saturn's lower stratosphere-diacetylene, butane, and water condense between approximately 1 and 300 mbar in our model and will dominate stratospheric haze formation at nonauroral latitudes. Our photochemical models will be useful for planning observational sequences and for analyzing data from the upcoming Cassini mission.

Multilayer integral method for simulation of eddy currents in thin volumes of arbitrary geometry produced by MRI gradient coils.

PubMed

Sanchez Lopez, Hector; Freschi, Fabio; Trakic, Adnan; Smith, Elliot; Herbert, Jeremy; Fuentes, Miguel; Wilson, Stephen; Liu, Limei; Repetto, Maurizio; Crozier, Stuart

2014-05-01

This article aims to present a fast, efficient and accurate multi-layer integral method (MIM) for the evaluation of complex spatiotemporal eddy currents in nonmagnetic and thin volumes of irregular geometries induced by arbitrary arrangements of gradient coils. The volume of interest is divided into a number of layers, wherein the thickness of each layer is assumed to be smaller than the skin depth and where one of the linear dimensions is much smaller than the remaining two dimensions. The diffusion equation of the current density is solved both in time-harmonic and transient domain. The experimentally measured magnetic fields produced by the coil and the induced eddy currents as well as the corresponding time-decay constants were in close agreement with the results produced by the MIM. Relevant parameters such as power loss and force induced by the eddy currents in a split cryostat were simulated using the MIM. The proposed method is capable of accurately simulating the current diffusion process inside thin volumes, such as the magnet cryostat. The method permits the priori-calculation of optimal pre-emphasis parameters. The MIM enables unified designs of gradient coil-magnet structures for an optimal mitigation of deleterious eddy current effects. Copyright © 2013 Wiley Periodicals, Inc.

Parameterized and resolved Southern Ocean eddy compensation

NASA Astrophysics Data System (ADS)

Poulsen, Mads B.; Jochum, Markus; Nuterman, Roman

2018-04-01

The ability to parameterize Southern Ocean eddy effects in a forced coarse resolution ocean general circulation model is assessed. The transient model response to a suite of different Southern Ocean wind stress forcing perturbations is presented and compared to identical experiments performed with the same model in 0.1° eddy-resolving resolution. With forcing of present-day wind stress magnitude and a thickness diffusivity formulated in terms of the local stratification, it is shown that the Southern Ocean residual meridional overturning circulation in the two models is different in structure and magnitude. It is found that the difference in the upper overturning cell is primarily explained by an overly strong subsurface flow in the parameterized eddy-induced circulation while the difference in the lower cell is mainly ascribed to the mean-flow overturning. With a zonally constant decrease of the zonal wind stress by 50% we show that the absolute decrease in the overturning circulation is insensitive to model resolution, and that the meridional isopycnal slope is relaxed in both models. The agreement between the models is not reproduced by a 50% wind stress increase, where the high resolution overturning decreases by 20%, but increases by 100% in the coarse resolution model. It is demonstrated that this difference is explained by changes in surface buoyancy forcing due to a reduced Antarctic sea ice cover, which strongly modulate the overturning response and ocean stratification. We conclude that the parameterized eddies are able to mimic the transient response to altered wind stress in the high resolution model, but partly misrepresent the unperturbed Southern Ocean meridional overturning circulation and associated heat transports.

Challenges in modeling spatiotemporally varying phytoplankton blooms in the Northwestern Arabian Sea and Gulf of Oman

NASA Astrophysics Data System (ADS)

Sedigh Marvasti, S.; Gnanadesikan, A.; Bidokhti, A. A.; Dunne, J. P.; Ghader, S.

2016-02-01

Recent years have shown an increase in harmful algal blooms in the Northwest Arabian Sea and Gulf of Oman, raising the question of whether climate change will accelerate this trend. This has led us to examine whether the Earth System Models used to simulate phytoplankton productivity accurately capture bloom dynamics in this region - both in terms of the annual cycle and interannual variability. Satellite data (SeaWIFS ocean color) show two climatological blooms in this region, a wintertime bloom peaking in February and a summertime bloom peaking in September. On a regional scale, interannual variability of the wintertime bloom is dominated by cyclonic eddies which vary in location from one year to another. Two coarse (1°) models with the relatively complex biogeochemistry (TOPAZ) capture the annual cycle but neither eddies nor the interannual variability. An eddy-resolving model (GFDL CM2.6) with a simpler biogeochemistry (miniBLING) displays larger interannual variability, but overestimates the wintertime bloom and captures eddy-bloom coupling in the south but not in the north. The models fail to capture both the magnitude of the wintertime bloom and its modulation by eddies in part because of their failure to capture the observed sharp thermocline and/or nutricline in this region. When CM2.6 is able to capture such features in the Southern part of the basin, eddies modulate diffusive nutrient supply to the surface (a mechanism not previously emphasized in the literature). For the model to simulate the observed wintertime blooms within cyclones, it will be necessary to represent this relatively unusual nutrient structure as well as the cyclonic eddies. This is a challenge in the Northern Arabian Sea as it requires capturing the details of the outflow from the Persian Gulf - something that is poorly done in global models.

High Efficiency, Low Distortion 3D Diffusion Tensor Imaging with Variable Density Spiral Fast Spin Echoes (3D DW VDS RARE)

PubMed Central

Frank, Lawrence R.; Jung, Youngkyoo; Inati, Souheil; Tyszka, J. Michael; Wong, Eric C.

2009-01-01

We present an acquisition and reconstruction method designed to acquire high resolution 3D fast spin echo diffusion tensor images while mitigating the major sources of artifacts in DTI - field distortions, eddy currents and motion. The resulting images, being 3D, are of high SNR, and being fast spin echoes, exhibit greatly reduced field distortions. This sequence utilizes variable density spiral acquisition gradients, which allow for the implementation of a self-navigation scheme by which both eddy current and motion artifacts are removed. The result is that high resolution 3D DTI images are produced without the need for eddy current compensating gradients or B0 field correction. In addition, a novel method for fast and accurate reconstruction of the non-Cartesian data is employed. Results are demonstrated in the brains of normal human volunteers. PMID:19778618

A two-equation model for heat transport in wall turbulent shear flows

NASA Astrophysics Data System (ADS)

Nagano, Y.; Kim, C.

1988-08-01

A new proposal for closing the energy equation is presented at the two-equation level of turbulence modeling. The eddy diffusivity concept is used in modeling. However, just as the eddy viscosity is determined from solutions of the k and epsilon equations, so the eddy diffusivity for heat is given as functions of temperature variance, and the dissipation rate of temperature fluctuations, together with k and epsilon. Thus, the proposed model does not require any questionable assumptions for the 'turbulent Prandtl number'. Modeled forms of the equations are developed to account for the physical effects of molecular Prandtl number and near-wall turbulence. The model is tested by application to a flat-plate boundary layer, the thermal entrance region of a pipe, and the turbulent heat transfer in fluids over a wide range of the Prandtl number. Agreement with the experiment is generally very satisfactory.

Vertical eddy diffusion coefficient from the LANDSAT imagery

NASA Technical Reports Server (NTRS)

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

1982-01-01

Analysis of five stable cases of the smoke plumes that originated in eastern Cabo Frio (22 deg 59'S; 42 deg 02'W), Brazil using LANDSAT imagery is presented for different months and years. From these images the lateral standard deviation (sigma sub y) and the lateral eddy diffusion coefficient (K sub y) are obtained from the formula based on Taylor's theory of diffusion by continuous moment. The rate of kinetic energy dissipation (e) is evaluated from the diffusion parameters sigma sub y and K sub y. Then, the vertical diffusion coefficient (K sub z) is estimated using Weinstock's formulation. These results agree well with the previous experimental values obtained over water surfaces by various workers. Values of e and K sub z show the weaker mixing processes in the marine stable boundary layer. The data sample is apparently to small to include representative active turbulent regions because such regions are so intermittent in time and in space. These results form a data base for use in the development and validation of mesoscale atmospheric diffusion models.

A sensitivity analysis of low salinity habitats simulated by a hydrodynamic model in the Manatee River estuary in Florida, USA

NASA Astrophysics Data System (ADS)

Chen, XinJian

2012-06-01

This paper presents a sensitivity study of simulated availability of low salinity habitats by a hydrodynamic model for the Manatee River estuary located in the southwest portion of the Florida peninsula. The purpose of the modeling study was to establish a regulatory minimum freshwater flow rate required to prevent the estuarine ecosystem from significant harm. The model used in the study was a multi-block model that dynamically couples a three-dimensional (3D) hydrodynamic model with a laterally averaged (2DV) hydrodynamic model. The model was calibrated and verified against measured real-time data of surface elevation and salinity at five stations during March 2005-July 2006. The calibrated model was then used to conduct a series of scenario runs to investigate effects of the flow reduction on salinity distributions in the Manatee River estuary. Based on simulated salinity distribution in the estuary, water volumes, bottom areas and shoreline lengths for salinity less than certain predefined values were calculated and analyzed to help establish the minimum freshwater flow rate for the estuarine system. The sensitivity analysis conducted during the modeling study for the Manatee River estuary examined effects of the bottom roughness, ambient vertical eddy viscosity/diffusivity, horizontal eddy viscosity/diffusivity, and ungauged flow on the model results and identified the relative importance of these model parameters (input data) to the outcome of the availability of low salinity habitats. It is found that the ambient vertical eddy viscosity/diffusivity is the most influential factor controlling the model outcome, while the horizontal eddy viscosity/diffusivity is the least influential one.

Improving our process understanding of methane emissions from a mid-latitude reservoir by combining eddy covariance monitoring with spatial surveys

EPA Science Inventory

Reservoirs are a globally important source of methane (CH4) to the atmosphere, but measuring CH4 emission rates from reservoirs is difficult due to the spatial and temporal variability of the various emission pathways, including ebullition and diffusion. We used the eddy covarian...

DNS/LES Simulations of Separated Flows at High Reynolds Numbers

NASA Technical Reports Server (NTRS)

Balakumar, P.

2015-01-01

Direct numerical simulations (DNS) and large-eddy simulations (LES) simulations of flow through a periodic channel with a constriction are performed using the dynamic Smagorinsky model at two Reynolds numbers of 2800 and 10595. The LES equations are solved using higher order compact schemes. DNS are performed for the lower Reynolds number case using a fine grid and the data are used to validate the LES results obtained with a coarse and a medium size grid. LES simulations are also performed for the higher Reynolds number case using a coarse and a medium size grid. The results are compared with an existing reference data set. The DNS and LES results agreed well with the reference data. Reynolds stresses, sub-grid eddy viscosity, and the budgets for the turbulent kinetic energy are also presented. It is found that the turbulent fluctuations in the normal and spanwise directions have the same magnitude. The turbulent kinetic energy budget shows that the production peaks near the separation point region and the production to dissipation ratio is very high on the order of five in this region. It is also observed that the production is balanced by the advection, diffusion, and dissipation in the shear layer region. The dominant term is the turbulent diffusion that is about two times the molecular dissipation.

Investigating the Eddy Diffusivity Concept in the Coastal Ocean

NASA Astrophysics Data System (ADS)

Rypina, I.; Kirincich, A.; Lentz, S. J.; Sundermeyer, M. A.

2016-12-01

We test the validity, utility, and limitations of the lateral eddy diffusivity concept in a coastal environment through analyzing data from coupled drifter and dye releases within the footprint of a high-resolution (800 m) high-frequency radar south of Martha's Vineyard, Massachusetts. Specifically, we investigate how well a combination of radar-based velocities and drifter-derived diffusivities can reproduce observed dye spreading over an 8-h time interval. A drifter-based estimate of an anisotropic diffusivity tensor is used to parameterize small-scale motions that are unresolved and under-resolved by the radar system. This leads to a significant improvement in the ability of the radar to reproduce the observed dye spreading. Our drifter-derived diffusivity estimates are O(10 m2/s), are consistent with the diffusivity inferred from aerial images of the dye taken using the quadcopter-mounted digital camera during the dye release, and are roughly an order of magnitude larger than diffusivity estimates of Okubo (O(1 m2/s)) for similar spatial scales ( 1 km). Despite the fact that the drifter-based diffusivity approach was successful in improving the ability of the radar to reproduce the observed dye spreading, the dispersion of drifters was, for the most part, not consistent with the diffusive spreading regime.

Transient tidal eddy motion in the western Gulf of Maine, part 1: Primary structure

NASA Astrophysics Data System (ADS)

Brown, W. S.; Marques, G. M.

2013-07-01

High frequency radar-derived surface current maps of the Great South Channel (GSC) in the western Gulf of Maine in 2005 revealed clockwise (CW) and anticlockwise (ACW) eddy motion associated with the strong regional tidal currents. To better elucidate the kinematics and dynamics of these transient tidal eddy motions, an observational and modeling study was conducted during the weakly stratified conditions of winter 2008-2009. Our moored bottom pressure and ADCP current measurements in 13m depth were augmented by historical current measurements in about 30m in documenting the dominance of highly polarized M2 semidiurnal currents in our nearshore study region. The high-resolution finite element coastal ocean model (QUODDY) - forced by the five principal tidal constituents - produced maps depicting the formation and evolution of the CW and ACW eddy motions that regularly follow maximum ebb and flood flows, respectively. Observation versus model current comparison required that the model bottom current drag coefficient be set to at an unusually high Cd=0.01 - suggesting the importance of form drag in the study region. The observations and model results were consistent in diagnosing CW or ACW eddy motions that (a) form nearshore in the coastal boundary layer (CBL) for about 3h after the respective tidal current maxima and then (b) translate southeastward across the GSC along curved 50m isobath at speeds of about 25m/s. Observation-based and model-based momentum budget estimates were consistent in showing a first order forced semidiurnal standing tidal wave dynamics (like the adjacent Gulf of Maine) which was modulated by adverse pressure gradient/bottom stress forcing to generate the eddy motions. Observation-based estimates of terms in the transport vorticity budget showed that in the shallower Inner Zone subregion (average depth=23m) that the diffusion of nearshore vorticity was dominant in feeding the growth of eddy motion vorticity; while in the somewhat deeper Outer Zone subregion (33m) bottom current lateral shear and water column stretching/squashing was significant in modulating the eddy motion. We conclude that the transient eddy motions in the GSC region are phase eddies that accompany the change of tide across the GSC and are (1) generated by bottom stress gradients in the shallower nearshore - an issue which needs to be better understood for improved future forecasting.

Flame-vortex interactions imaged in microgravity

NASA Technical Reports Server (NTRS)

Driscoll, James F.; Dahm, Werner J. A.; Sichel, Martin

1995-01-01

The scientific objective is to obtain high quality color-enhanced digital images of a vortex exerting aerodynamic strain on premixed and nonpremixed flames with the complicating effects of buoyancy removed. The images will provide universal (buoyancy free) scaling relations that are required to improve several types of models of turbulent combustion, including KIVA-3, discrete vortex, and large-eddy simulations. The images will be used to help quantify several source terms in the models, including those due to flame stretch, flame-generated vorticity, flame curvature, and preferential diffusion, for a range of vortex sizes and flame conditions. The experiment is an ideal way to study turbulence-chemistry interactions and isolate the effect of vortices of different sizes and strengths in a repeatable manner. A parallel computational effort is being conducted which considers full chemistry and preferential diffusion.

Air-Parcel Residence Times Within Forest Canopies

NASA Astrophysics Data System (ADS)

Gerken, Tobias; Chamecki, Marcelo; Fuentes, Jose D.

2017-10-01

We present a theoretical model, based on a simple model of turbulent diffusion and first-order chemical kinetics, to determine air-parcel residence times and the out-of-canopy export of reactive gases emitted within forest canopies under neutral conditions. Theoretical predictions of the air-parcel residence time are compared to values derived from large-eddy simulation for a range of canopy architectures and turbulence levels under neutral stratification. Median air-parcel residence times range from a few sec in the upper canopy to approximately 30 min near the ground and the distribution of residence times is skewed towards longer times in the lower canopy. While the predicted probability density functions from the theoretical model and large-eddy simulation are in good agreement with each other, the theoretical model requires only information on canopy height and eddy diffusivities inside the canopy. The eddy-diffusivity model developed additionally requires the friction velocity at canopy top and a parametrized profile of the standard deviation of vertical velocity. The theoretical model of air-parcel residence times is extended to include first-order chemical reactions over a range of of Damköhler numbers ( Da) characteristic of plant-emitted hydrocarbons. The resulting out-of-canopy export fractions range from near 1 for Da =10^{-3} to less than 0.3 at Da = 10. These results highlight the necessity for dense and tall forests to include the impacts of air-parcel residence times when calculating the out-of-canopy export fraction for reactive trace gases.

Estimates of the lateral eddy diffusivity in the Indian Ocean as derived from drifter data

NASA Astrophysics Data System (ADS)

Zhurbas, V. M.; Lyzhkov, D. A.; Kuzmina, N. P.

2014-05-01

The Global Drifter Program data set is applied to develop 2° × 2° bin estimates of the lateral eddy diffusivity K in the Indian Ocean (IO) by means of a modification of the Davis approach. The calculations were performed relative to the seasonal change in the mean currents, which is especially important in the case of monsoon-driven circulation in the IO. Estimates of K were found to be below 1 × 104 m2/s almost every-where in the IO. The spatial variations of K were analyzed in relation to the instabilities of the ocean circulation.

A theoretically based determination of bowen-ratio fetch requirements

USGS Publications Warehouse

Stannard, D.I.

1997-01-01

Determination of fetch requirements for accurate Bowen-ratio measurements of latent- and sensible-heat fluxes is more involved than for eddy-correlation measurements because Bowen-ratio sensors are located at two heights, rather than just one. A simple solution to the diffusion equation is used to derive an expression for Bowen-ratio fetch requirements, downwind of a step change in surface fluxes. These requirements are then compared to eddy-correlation fetch requirements based on the same diffusion equation solution. When the eddy-correlation and upper Bowen-ratio sensor heights are equal, and the available energy upwind and downwind of the step change is constant, the Bowen-ratio method requires less fetch than does eddy correlation. Differences in fetch requirements between the two methods are greatest over relatively smooth surfaces. Bowen-ratio fetch can be reduced significantly by lowering the lower sensor, as well as the upper sensor. The Bowen-ratio fetch model was tested using data from a field experiment where multiple Bowen-ratio systems were deployed simultaneously at various fetches and heights above a field of bermudagrass. Initial comparisons were poor, but improved greatly when the model was modified (and operated numerically) to account for the large roughness of the upwind cotton field.

Simulations of the stratocumulus-topped boundary layer with a third-order closure model

NASA Technical Reports Server (NTRS)

Moeng, C. H.; Randall, D. A.

1984-01-01

A third order closure model is proposed by Andre et al. (1982), in which the time rate of change terms, the relaxation and rapid effects for the pressure related terms, and the clipping approximation are included along with the quasi-normal closure, to study turbulence in a cloudy layer which is cooled radiatively from above. A spurious oscillation which is strongest near the inversion occurs. An analysis of the problem shows that the oscillation arises from the mean gradient and buoyancy terms of the triple moment equations; these terms are largest near the cloud top. The oscillation is physical, rather than computational. In nature the oscillation is effectively damped, by a mechanism which apparently is not included in our model. In the stably stratified layer just above the mixed layer top, turbulence can excite gravity waves, whose energy is radiated away. Because the closure assumption for the pressure terms does not take into account the transport of wave energy, the model generates spurious oscillations. Damping of the oscillations is possible by introducing diffusion terms into the triple moment equations. With a large enough choice for the diffusion coefficient, the oscillation is effectively eliminated. The results are quite sensitive to the ad hoc eddy coefficient.

Development of Eddy Current Techniques for Detection of Deep Fatigue Cracks in Multi-Layer Airframe Components

NASA Technical Reports Server (NTRS)

Wincheski, Russell A.

2008-01-01

Thick, multi-layer aluminum structure has been widely used in aircraft design in critical wing splice areas. The multi-layer structure generally consists of three or four aluminum layers with different geometry and varying thickness, which are held together with fasteners. The detection of cracks under fasteners with ultrasonic techniques in subsurface layers away from the skin is impeded primarily by interlayer bonds and faying sealant condition. Further, assessment of such sealant condition is extremely challenging in terms of complexity of structure, limited access, and inspection cost. Although Eddy current techniques can be applied on in-service aircraft from the exterior of the skin without knowing sealant condition, the current eddy current techniques are not able to detect defects with wanted sensitivity. In this work a series of low frequency eddy current probes have been designed, fabricated and tested for this application. A probe design incorporating a shielded magnetic field sensor concentrically located in the interior of a drive coil has been employed to enable a localized deep diffusion of the electromagnetic field into the part under test. Due to the required low frequency inspections, probes have been testing using a variety of magnetic field sensors (pickup coil, giant magneto-resistive, anisotropic magneto-resistive, and spin-dependent tunneling). The probe designs as well as capabilities based upon a target inspection for sub-layer cracking in an airframe wing spar joint is presented.

A priori study of subgrid-scale features in turbulent Rayleigh-Bénard convection

NASA Astrophysics Data System (ADS)

Dabbagh, F.; Trias, F. X.; Gorobets, A.; Oliva, A.

2017-10-01

At the crossroad between flow topology analysis and turbulence modeling, a priori studies are a reliable tool to understand the underlying physics of the subgrid-scale (SGS) motions in turbulent flows. In this paper, properties of the SGS features in the framework of a large-eddy simulation are studied for a turbulent Rayleigh-Bénard convection (RBC). To do so, data from direct numerical simulation (DNS) of a turbulent air-filled RBC in a rectangular cavity of aspect ratio unity and π spanwise open-ended distance are used at two Rayleigh numbers R a ∈{1 08,1 010 } [Dabbagh et al., "On the evolution of flow topology in turbulent Rayleigh-Bénard convection," Phys. Fluids 28, 115105 (2016)]. First, DNS at Ra = 108 is used to assess the performance of eddy-viscosity models such as QR, Wall-Adapting Local Eddy-viscosity (WALE), and the recent S3PQR-models proposed by Trias et al. ["Building proper invariants for eddy-viscosity subgrid-scale models," Phys. Fluids 27, 065103 (2015)]. The outcomes imply that the eddy-viscosity modeling smoothes the coarse-grained viscous straining and retrieves fairly well the effect of the kinetic unfiltered scales in order to reproduce the coherent large scales. However, these models fail to approach the exact evolution of the SGS heat flux and are incapable to reproduce well the further dominant rotational enstrophy pertaining to the buoyant production. Afterwards, the key ingredients of eddy-viscosity, νt, and eddy-diffusivity, κt, are calculated a priori and revealed positive prevalent values to maintain a turbulent wind essentially driven by the mean buoyant force at the sidewalls. The topological analysis suggests that the effective turbulent diffusion paradigm and the hypothesis of a constant turbulent Prandtl number are only applicable in the large-scale strain-dominated areas in the bulk. It is shown that the bulk-dominated rotational structures of vortex-stretching (and its synchronous viscous dissipative structures) hold the highest positive values of νt; however, the zones of backscatter energy and counter-gradient heat transport are related to the areas of compressed focal vorticity. More arguments have been attained through a priori investigation of the alignment trends imposed by existing parameterizations for the SGS heat flux, tested here inside RBC. It is shown that the parameterizations based linearly on the resolved thermal gradient are invalid in RBC. Alternatively, the tensor-diffusivity approach becomes a crucial choice of modeling the SGS heat flux, in particular, the tensorial diffusivity that includes the SGS stress tensor. This and other crucial scrutinies on a future modeling to the SGS heat flux in RBC are sought.

Large eddy simulation modeling of particle-laden flows in complex terrain

NASA Astrophysics Data System (ADS)

Salesky, S.; Giometto, M. G.; Chamecki, M.; Lehning, M.; Parlange, M. B.

2017-12-01

The transport, deposition, and erosion of heavy particles over complex terrain in the atmospheric boundary layer is an important process for hydrology, air quality forecasting, biology, and geomorphology. However, in situ observations can be challenging in complex terrain due to spatial heterogeneity. Furthermore, there is a need to develop numerical tools that can accurately represent the physics of these multiphase flows over complex surfaces. We present a new numerical approach to accurately model the transport and deposition of heavy particles in complex terrain using large eddy simulation (LES). Particle transport is represented through solution of the advection-diffusion equation including terms that represent gravitational settling and inertia. The particle conservation equation is discretized in a cut-cell finite volume framework in order to accurately enforce mass conservation. Simulation results will be validated with experimental data, and numerical considerations required to enforce boundary conditions at the surface will be discussed. Applications will be presented in the context of snow deposition and transport, as well as urban dispersion.

Radial distribution of the contributions to band broadening of a silica-based semi-preparative monolithic column.

PubMed

Abia, Jude A; Mriziq, Khaled S; Guiochon, Georges A

2009-04-01

Using an on-column local electrochemical microdetector operated in the amperometric mode, band elution profiles were recorded at different radial locations at the exit of a 10 mm id, 100 mm long silica-based monolithic column. HETP plots were then acquired at each of these locations, and all these results were fitted to the Knox equation. This provided a spatial distribution of the values of the eddy diffusion (A), the molecular diffusion (B), and the resistance to the kinetics of mass transfer (C) terms. Results obtained indicate that the wall region yields higher A values and smaller C values than the central core region. Significant radial fluctuations of these contributions to band broadening occur throughout the exit column cross-section. This phenomenon is due to the structural radial heterogeneity of the column.

A second-order closure analysis of turbulent diffusion flames. [combustion physics

NASA Technical Reports Server (NTRS)

Varma, A. K.; Fishburne, E. S.; Beddini, R. A.

1977-01-01

A complete second-order closure computer program for the investigation of compressible, turbulent, reacting shear layers was developed. The equations for the means and the second order correlations were derived from the time-averaged Navier-Stokes equations and contain third order and higher order correlations, which have to be modeled in terms of the lower-order correlations to close the system of equations. In addition to fluid mechanical turbulence models and parameters used in previous studies of a variety of incompressible and compressible shear flows, a number of additional scalar correlations were modeled for chemically reacting flows, and a typical eddy model developed for the joint probability density function for all the scalars. The program which is capable of handling multi-species, multistep chemical reactions, was used to calculate nonreacting and reacting flows in a hydrogen-air diffusion flame.

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

Wolfram, Phillip J.; Ringler, Todd D.; Maltrud, Mathew E.

Isopycnal diffusivity due to stirring by mesoscale eddies in an idealized, wind-forced, eddying, midlatitude ocean basin is computed using Lagrangian, in Situ, Global, High-Performance Particle Tracking (LIGHT). Simulation is performed via LIGHT within the Model for Prediction across Scales Ocean (MPAS-O). Simulations are performed at 4-, 8-, 16-, and 32-km resolution, where the first Rossby radius of deformation (RRD) is approximately 30 km. Scalar and tensor diffusivities are estimated at each resolution based on 30 ensemble members using particle cluster statistics. Each ensemble member is composed of 303 665 particles distributed across five potential density surfaces. Diffusivity dependence upon modelmore » resolution, velocity spatial scale, and buoyancy surface is quantified and compared with mixing length theory. The spatial structure of diffusivity ranges over approximately two orders of magnitude with values of O(10 5) m 2 s –1 in the region of western boundary current separation to O(10 3) m 2 s –1 in the eastern region of the basin. Dominant mixing occurs at scales twice the size of the first RRD. Model resolution at scales finer than the RRD is necessary to obtain sufficient model fidelity at scales between one and four RRD to accurately represent mixing. Mixing length scaling with eddy kinetic energy and the Lagrangian time scale yield mixing efficiencies that typically range between 0.4 and 0.8. In conclusion, a reduced mixing length in the eastern region of the domain relative to the west suggests there are different mixing regimes outside the baroclinic jet region.« less

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.

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.

QNSE Theory of Turbulence in Rotating Fluids and the Nastrom & Gage Spectrum

NASA Astrophysics Data System (ADS)

Galperin, B.

2017-12-01

An analytical theory of turbulence, the quasi-normal scale elimination (QNSE), has been developed for neutrally stratified rotating flows. The theory provides near-first principle framework for the representation of flow anisotropization under the action of rotation. The anisotropization reveals itself in the emergence of different eddy viscosities and eddy diffusivities in different directions and directional dependence of the kinetic and potential energies spectra. In addition, there are also phenomena of componentality, eddy viscosities are different for different velocity components, and the onset of the inverse energy cascade. The anisotropization increases with increasing scale. The characteristic scales for the crossover between the turbulence and inertial wave domains is the Woods scale, LΩ = [ɛ/(2Ω)3)]1/2, ɛ being the rate of the viscous dissipation, which is analogous to the Ozmidov scale in flows with stable stratification. Rapid rotation renders the horizontal eddy viscosity negative, and in order to preserve it positive, a weak rotation limit is invoked. Within that limit, an analytical theory of the transition from the Kolmogorov to a rotation-dominated turbulence regime is developed. The dispersion relation of linear inertial waves is unaffected by turbulence while all one-dimensional energy spectra undergo steepening from the Kolmogorov -5/3 to the -3 slope. The longitudinal and transverse spectra are congruent with the famous atmospheric spectra by Nastrom & Gage. Thus, for the first time, these spectra are obtained within an analytical theory. QNSE explains the latitudinal dependence of the spectra and lends itself for practical applications in simulations of atmospheric and oceanic flows as it produces closed expressions for the eddy viscosities and eddy diffusivities. The Nastrom & Gage spectra also apply to the oceanic flows.

Kinetic investigation of narrow-bore columns packed with prototype sub-2 μm superficially porous particles with various shell thickness.

PubMed

Gritti, Fabrice; Omamogho, Jesse; Guiochon, Georges

2011-10-07

The recent successful breakthrough of sub-3 μm shell particles in HPLC has triggered considerable research efforts toward the design of new brands of core-shell particles. We investigated the mass transfer mechanism of a few analytes in narrow-bore columns packed with prototype 1.7 μm shell particles, made of 1.0, 1.2, and 1.4 μm solid nonporous cores surrounded by porous shells 350, 250, and 150 nm thick, respectively. Three probe solutes, uracil, naphthalene, and insulin, were chosen to assess the kinetic performance of these columns. Inverse size exclusion chromatography, peak parking experiments, and the numerical integration of the experimental peak profiles were carried out in order to measure the external, internal, and total column porosities, the true bulk diffusion coefficients of these analytes, the height equivalent to a theoretical plate, the longitudinal diffusion term, and the trans-particle mass transfer resistance term. The residual eddy diffusion term was measured by difference. The results show the existence of important trans-column velocity biases (7%) possibly due to the presence of particle multiplets in the slurry mixture used during the packing process. Our results illustrates some of the difficulties encountered by scientists preparing and packing shell particles into narrow-bore columns. Copyright © 2011 Elsevier B.V. All rights reserved.

Anisotropic Shear Dispersion Parameterization for Mesoscale Eddy Transport

NASA Astrophysics Data System (ADS)

Reckinger, S. J.; Fox-Kemper, B.

2016-02-01

The effects of mesoscale eddies are universally treated isotropically in general circulation models. However, the processes that the parameterization approximates, such as shear dispersion, typically have strongly anisotropic characteristics. The Gent-McWilliams/Redi mesoscale eddy parameterization is extended for anisotropy and tested using 1-degree Community Earth System Model (CESM) simulations. The sensitivity of the model to anisotropy includes a reduction of temperature and salinity biases, a deepening of the southern ocean mixed-layer depth, and improved ventilation of biogeochemical tracers, particularly in oxygen minimum zones. The parameterization is further extended to include the effects of unresolved shear dispersion, which sets the strength and direction of anisotropy. The shear dispersion parameterization is similar to drifter observations in spatial distribution of diffusivity and high-resolution model diagnosis in the distribution of eddy flux orientation.

Comparison of Four Mixed Layer Mesoscale Parameterizations and the Equation for an Arbitrary Tracer

NASA Technical Reports Server (NTRS)

Canuto, V. M.; Dubovikov, M. S.

2011-01-01

In this paper we discuss two issues, the inter-comparison of four mixed layer mesoscale parameterizations and the search for the eddy induced velocity for an arbitrary tracer. It must be stressed that our analysis is limited to mixed layer mesoscales since we do not treat sub-mesoscales and small turbulent mixing. As for the first item, since three of the four parameterizations are expressed in terms of a stream function and a residual flux of the RMT formalism (residual mean theory), while the fourth is expressed in terms of vertical and horizontal fluxes, we needed a formalism to connect the two formulations. The standard RMT representation developed for the deep ocean cannot be extended to the mixed layer since its stream function does not vanish at the ocean's surface. We develop a new RMT representation that satisfies the surface boundary condition. As for the general form of the eddy induced velocity for an arbitrary tracer, thus far, it has been assumed that there is only the one that originates from the curl of the stream function. This is because it was assumed that the tracer residual flux is purely diffusive. On the other hand, we show that in the case of an arbitrary tracer, the residual flux has also a skew component that gives rise to an additional bolus velocity. Therefore, instead of only one bolus velocity, there are now two, one coming from the curl of the stream function and other from the skew part of the residual flux. In the buoyancy case, only one bolus velocity contributes to the mean buoyancy equation since the residual flux is indeed only diffusive.

Joint correction of Nyquist artifact and minuscule motion-induced aliasing artifact in interleaved diffusion weighted EPI data using a composite two-dimensional phase correction procedure

PubMed Central

Chang, Hing-Chiu; Chen, Nan-kuei

2016-01-01

Diffusion-weighted imaging (DWI) obtained with interleaved echo-planar imaging (EPI) pulse sequence has great potential of characterizing brain tissue properties at high spatial-resolution. However, interleaved EPI based DWI data may be corrupted by various types of aliasing artifacts. First, inconsistencies in k-space data obtained with opposite readout gradient polarities result in Nyquist artifact, which is usually reduced with 1D phase correction in post-processing. When there exist eddy current cross terms (e.g., in oblique-plane EPI), 2D phase correction is needed to effectively reduce Nyquist artifact. Second, minuscule motion induced phase inconsistencies in interleaved DWI scans result in image-domain aliasing artifact, which can be removed with reconstruction procedures that take shot-to-shot phase variations into consideration. In existing interleaved DWI reconstruction procedures, Nyquist artifact and minuscule motion-induced aliasing artifact are typically removed subsequently in two stages. Although the two-stage phase correction generally performs well for non-oblique plane EPI data obtained from well-calibrated system, the residual artifacts may still be pronounced in oblique-plane EPI data or when there exist eddy current cross terms. To address this challenge, here we report a new composite 2D phase correction procedure, which effective removes Nyquist artifact and minuscule motion induced aliasing artifact jointly in a single step. Our experimental results demonstrate that the new 2D phase correction method can much more effectively reduce artifacts in interleaved EPI based DWI data as compared with the existing two-stage artifact correction procedures. The new method robustly enables high-resolution DWI, and should prove highly valuable for clinical uses and research studies of DWI. PMID:27114342

Large Eddy Simulation of Gravitational Effects on Transitional and Turbulent Gas-Jet Diffusion Flames

NASA Technical Reports Server (NTRS)

Givi, Peyman; Jaberi, Farhad A.

2001-01-01

The basic objective of this work is to assess the influence of gravity on "the compositional and the spatial structures" of transitional and turbulent diffusion flames via large eddy simulation (LES), and direct numerical simulation (DNS). The DNS is conducted for appraisal of the various closures employed in LES, and to study the effect of buoyancy on the small scale flow features. The LES is based on our "filtered mass density function"' (FMDF) model. The novelty of the methodology is that it allows for reliable simulations with inclusion of "realistic physics." It also allows for detailed analysis of the unsteady large scale flow evolution and compositional flame structure which is not usually possible via Reynolds averaged simulations.

A simple predictive model for the structure of the oceanic pycnocline

PubMed

Gnanadesikan

1999-03-26

A simple theory for the large-scale oceanic circulation is developed, relating pycnocline depth, Northern Hemisphere sinking, and low-latitude upwelling to pycnocline diffusivity and Southern Ocean winds and eddies. The results show that Southern Ocean processes help maintain the global ocean structure and that pycnocline diffusion controls low-latitude upwelling.

Eddy diffusivity of quasi-neutrally-buoyant inertial particles

NASA Astrophysics Data System (ADS)

Martins Afonso, Marco; Muratore-Ginanneschi, Paolo; Gama, Sílvio M. A.; Mazzino, Andrea

2018-04-01

We investigate the large-scale transport properties of quasi-neutrally-buoyant inertial particles carried by incompressible zero-mean periodic or steady ergodic flows. We show how to compute large-scale indicators such as the inertial-particle terminal velocity and eddy diffusivity from first principles in a perturbative expansion around the limit of added-mass factor close to unity. Physically, this limit corresponds to the case where the mass density of the particles is constant and close in value to the mass density of the fluid, which is also constant. Our approach differs from the usual over-damped expansion inasmuch as we do not assume a separation of time scales between thermalization and small-scale convection effects. For a general flow in the class of incompressible zero-mean periodic velocity fields, we derive closed-form cell equations for the auxiliary quantities determining the terminal velocity and effective diffusivity. In the special case of parallel flows these equations admit explicit analytic solution. We use parallel flows to show that our approach sheds light onto the behavior of terminal velocity and effective diffusivity for Stokes numbers of the order of unity.

The Dynamics of Helium and its Impact on the Upper Thermosphere

NASA Astrophysics Data System (ADS)

Sutton, E. K.; Thayer, J. P.; Wang, W.; Solomon, S. C.; Schmidt, F.

2015-12-01

The TIE-GCM was recently augmented to include helium and argon, two approximately inert species that can be used as tracers of dynamics in the thermosphere. The former species is treated as a major species due to its large abundance near the upper boundary. The effects of exospheric transport are also included in order to simulate realistic seasonal and latitudinal helium distributions. The latter species is treated as a classical minor species, imparting absolutely no forces on the background atmosphere. In this study, we examine the interplay of the various dynamical terms - i.e. background circulation, molecular and Eddy diffusion - as they drive departures from the distributions that would be expected under assumptions of diffusive equilibrium. As this has implications on the formulation of all semi-empirical thermospheric models, we use this understanding to identify the conditions under which helium can significantly affect nowcasts and forecasts of neutral density.

A multi-scalar PDF approach for LES of turbulent spray combustion

NASA Astrophysics Data System (ADS)

Raman, Venkat; Heye, Colin

2011-11-01

A comprehensive joint-scalar probability density function (PDF) approach is proposed for large eddy simulation (LES) of turbulent spray combustion and tests are conducted to analyze the validity and modeling requirements. The PDF method has the advantage that the chemical source term appears closed but requires models for the small scale mixing process. A stable and consistent numerical algorithm for the LES/PDF approach is presented. To understand the modeling issues in the PDF method, direct numerical simulation of a spray flame at three different fuel droplet Stokes numbers and an equivalent gaseous flame are carried out. Assumptions in closing the subfilter conditional diffusion term in the filtered PDF transport equation are evaluated for various model forms. In addition, the validity of evaporation rate models in high Stokes number flows is analyzed.

Rossby wave activity in a two-dimensional model - Closure for wave driving and meridional eddy diffusivity

NASA Technical Reports Server (NTRS)

Hitchman, Matthew H.; Brasseur, Guy

1988-01-01

A parameterization of the effects of Rossby waves in the middle atmosphere is proposed for use in two-dimensional models. By adding an equation for conservation of Rossby wave activity, closure is obtained for the meridional eddy fluxes and body force due to Rossby waves. Rossby wave activity is produced in a climatological fashion at the tropopause, is advected by a group velocity which is determined solely by model zonal winds, and is absorbed where it converges. Absorption of Rossby wave activity causes both an easterly torque and an irreversible mixing of potential vorticity, represented by the meridional eddy diffusivity, K(yy). The distribution of Rossby wave driving determines the distribution of K(yy), which is applied to all of the chemical constituents. This provides a self-consistent coupling of the wave activity with the winds, tracer distributions and the radiative field. Typical winter stratospheric values for K(yy) of 2 million sq m/sec are obtained. Poleward tracer advection is enhanced and meridional tracer gradients are reduced where Rossby wave activity is absorbed in the model.

Fractional Diffusion Analysis of the Electromagnetic Field In Fractured Media Part II: 2.5-D Approach

NASA Astrophysics Data System (ADS)

Ge, J.; Everett, M. E.; Weiss, C. J.

2012-12-01

A 2.5D finite difference (FD) frequency-domain modeling algorithm based on the theory of fractional diffusion of electromagnetic (EM) fields generated by a loop source lying above a fractured geological medium is addressed in this paper. The presence of fractures in the subsurface, usually containing highly conductive pore fluids, gives rise to spatially hierarchical flow paths of induced EM eddy currents. The diffusion of EM eddy currents in such formations is anomalous, generalizing the classical Gaussian process described by the conventional Maxwell equations. Based on the continuous time random walk (CTRW) theory, the diffusion of EM eddy currents in a rough medium is governed by the fractional Maxwell equations. Here, we model the EM response of a 2D subsurface containing fractured zones, with a 3D loop source, which results the so-called 2.5D model geometry. The governing equation in the frequency domain is converted using Fourier transform into k domain along the strike direction (along which the model conductivity doesn't vary). The resulting equation system is solved by the multifrontal massively parallel solver (MUMPS). The data obtained is then converted back to spatial domain and the time domain. We find excellent agreement between the FD and analytic solutions for a rough halfspace model. Then FD solutions are calculated for a 2D fault zone model with variable conductivity and roughness. We compare the results with responses from several classical models and explore the relationship between the roughness and the spatial density of the fracture distribution.

Magnetic loss and B(H) behaviour of non-oriented electrical sheets under a trapezoidal exciting field

NASA Astrophysics Data System (ADS)

Kedous-Lebouc, A.; Errard, S.; Cornut, B.; Brissonneau, P.

1994-05-01

The excess loss and hysteresis response of electrical steel are measured and discussed in the case of trapezoidal field excitation similar to the current provided by a current commutation supply of a self-synchronous rotating machine. Three industrial non-oriented SiFe samples of different magnetic grades and thicknesses are tested using an automatic Epstein frame equipment. The losses and the unusual observed B( H) loops are analysed in terms of the rate of change of the field, the diffusion of the induction inside the sheet and by the calculation of the theoretical hysteresis cycles due to the eddy currents.

Development of a second order closure model for computation of turbulent diffusion flames

NASA Technical Reports Server (NTRS)

Varma, A. K.; Donaldson, C. D.

1974-01-01

A typical eddy box model for the second-order closure of turbulent, multispecies, reacting flows developed. The model structure was quite general and was valid for an arbitrary number of species. For the case of a reaction involving three species, the nine model parameters were determined from equations for nine independent first- and second-order correlations. The model enabled calculation of any higher-order correlation involving mass fractions, temperatures, and reaction rates in terms of first- and second-order correlations. Model predictions for the reaction rate were in very good agreement with exact solutions of the reaction rate equations for a number of assumed flow distributions.

Comparison between the intra-particle diffusivity in the hydrophilic interaction chromatography and reversed phase liquid chromatography modes. Impact on the column efficiency.

PubMed

Gritti, Fabrice; Guiochon, Georges

2013-07-05

The effective diffusion coefficients of five low molecular weigh compounds (naphthalene, uracil, uridine, adenosine, and cytosine) were measured at room temperature in a 4.6mm×100mm column packed with 3.5μm XBridge HILIC particles. The mobile phase was an acetonitrile-water mixture (92.5/7.5, v/v) containing 10mM ammonium acetate and 0.02% acetic acid. Using a physically reliable model of effective diffusion in binary composite media (Torquato's model), accurate estimates of the intra-particle diffusivities in the HILIC particles were obtained as a function of the retention of these analytes. The HILIC diffusion coefficients were compared to those previously obtained for endcapped RPLC-C18 particles (5.0μm Gemini-C18). The experimental results confirm that adsorption sites are not localized in RPLC whereas they are so in the HILIC mode. In contrast to RPLC columns, HILIC columns provide longitudinal diffusion B/u terms that increase very little with increasing retention factors. This confirms the absence of surface diffusion in HILIC. The impact of intra-particle diffusivity on the column efficiency was projected in HILIC and RPLC on the basis of the measured intra-particle diffusivities and on the well established theory of band broadening in particulate columns. Accordingly, RPLC columns generate short-range eddy dispersion and solid-liquid mass transfer resistance Cu terms that increase less than do HILIC column with increasing retention factors. The HETP contribution caused by the trans-column structure heterogeneity is smaller in the HILIC than in the RPLC modes because the transverse excursion length is smaller in HILIC. Even though the overall column efficiencies are comparable in HILIC and RPLC, this study shows that the individual mass transfer phenomena are inherently different in the HILIC and the RPLC retention modes. Copyright © 2013 Elsevier B.V. All rights reserved.

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...

Boundary Layer Model for Air Pollutant Concentrations Due to Highway Traffic

ERIC Educational Resources Information Center

Ragland, Kenneth W.; Peirce, J. Jeffrey

1975-01-01

A numerical solution of the three-dimensional steady-state diffusion equation for a finite width line source is presented. The wind speed and eddy diffusivity as a function of height above the roadway are obtained. Normalized ground level and elevated concentrations near a highway are obtained for winds perpendicular, parallel, and at 45 degrees.…

Eulerian and Lagrangian Parameterization of the Oceanic Mixed Layer using Large Eddy Simulation and MPAS-Ocean

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

Van Roekel, Luke

We have conducted a suite of Large Eddy Simulation (LES) to form the basis of a multi-model comparison (left). The results have led to proposed model improvements. We have verified that Eulerian-Lagrangian effective diffusivity estimates of mesoscale mixing are consistent with traditional particle statistics metrics (right). LES and Lagrangian particles will be utilized to better represent the movement of water into and out of the mixed layer.

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.

On the phase lag of turbulent dissipation in rotating tidal flows

NASA Astrophysics Data System (ADS)

Zhang, Qianjiang; Wu, Jiaxue

2018-03-01

Field observations of rotating tidal flows in a shallow tidally swept sea reveal that a notable phase lag of both shear production and turbulent dissipation increases with height above the seafloor. These vertical delays of turbulent quantities are approximately equivalent in magnitude to that of squared mean shear. The shear production approximately equals turbulent dissipation over the phase-lag column, and thus a main mechanism of phase lag of dissipation is mean shear, rather than vertical diffusion of turbulent kinetic energy. By relating the phase lag of dissipation to that of the mean shear, a simple formulation with constant eddy viscosity is developed to describe the phase lag in rotating tidal flows. An analytical solution indicates that the phase lag increases linearly with height subjected to a combined effect of tidal frequency, Coriolis parameter and eddy viscosity. The vertical diffusion of momentum associated with eddy viscosity produces the phase lag of squared mean shear, and resultant delay of turbulent quantities. Its magnitude is inhibited by Earth's rotation. Furthermore, a theoretical formulation of the phase lag with a parabolic eddy viscosity profile can be constructed. A first-order approximation of this formulation is still a linear function of height, and its magnitude is approximately 0.8 times that with constant viscosity. Finally, the theoretical solutions of phase lag with realistic viscosity can be satisfactorily justified by realistic phase lags of dissipation.

Optimization of diffusion-weighted single-refocused spin-echo EPI by reducing eddy-current artifacts and shortening the echo time.

PubMed

Shrestha, Manoj; Hok, Pavel; Nöth, Ulrike; Lienerth, Bianca; Deichmann, Ralf

2018-03-30

The purpose of this work was to optimize the acquisition of diffusion-weighted (DW) single-refocused spin-echo (srSE) data without intrinsic eddy-current compensation (ECC) for an improved performance of ECC postprocessing. The rationale is that srSE sequences without ECC may yield shorter echo times (TE) and thus higher signal-to-noise ratios (SNR) than srSE or twice-refocused spin-echo (trSE) schemes with intrinsic ECC. The proposed method employs dummy scans with DW gradients to drive eddy currents into a steady state before data acquisition. Parameters of the ECC postprocessing algorithm were also optimized. Simulations were performed to obtain minimum TE values for the proposed sequence and sequences with intrinsic ECC. Experimentally, the proposed method was compared with standard DW-trSE imaging, both in vitro and in vivo. Simulations showed substantially shorter TE for the proposed method than for methods with intrinsic ECC when using shortened echo readouts. Data of the proposed method showed a marked increase in SNR. A dummy scan duration of at least 1.5 s improved performance of the ECC postprocessing algorithm. Changes proposed for the DW-srSE sequence and for the parameter setting of the postprocessing ECC algorithm considerably reduced eddy-current artifacts and provided a higher SNR.

Comparisons between thermodynamic and one-dimensional combustion models of spark-ignition engines

NASA Technical Reports Server (NTRS)

Ramos, J. I.

1986-01-01

Results from a one-dimensional combustion model employing a constant eddy diffusivity and a one-step chemical reaction are compared with those of one-zone and two-zone thermodynamic models to study the flame propagation in a spark-ignition engine. One-dimensional model predictions are found to be very sensitive to the eddy diffusivity and reaction rate data. The average mixing temperature found using the one-zone thermodynamic model is higher than those of the two-zone and one-dimensional models during the compression stroke, and that of the one-dimensional model is higher than those predicted by both thermodynamic models during the expansion stroke. The one-dimensional model is shown to predict an accelerating flame even when the front approaches the cold cylinder wall.

Dynamical influences on thermospheric composition: implications for semi-empirical models

NASA Astrophysics Data System (ADS)

Sutton, E. K.; Solomon, S. C.

2014-12-01

The TIE-GCM was recently augmented to include helium and argon, two approximately inert species that can be used as tracers of dynamics in the thermosphere. The former species is treated as a major species due to its large abundance near the upper boundary. The effects of exospheric transport are also included in order to simulate realistic seasonal and latitudinal helium distributions. The latter species is treated as a classical minor species, imparting absolutely no forces on the background atmosphere. In this study, we examine the interplay of the various dynamical terms - i.e. background circulation, molecular and Eddy diffusion - as they drive departures from the distributions that would be expected under the assumption of diffusive equilibrium. As this has implications on the formulation of all empirical thermospheric models, we use this understanding to address the following questions: (1) how do errors caused by the assumption of diffusive equilibrium manifest within empirical models of the thermosphere? and (2) where and when does an empirical model's output disagree with its underlying datasets due to the inherent limitations of said model's formulation?

Large Eddy Simulation of Entropy Generation in a Turbulent Mixing Layer

NASA Astrophysics Data System (ADS)

Sheikhi, Reza H.; Safari, Mehdi; Hadi, Fatemeh

2013-11-01

Entropy transport equation is considered in large eddy simulation (LES) of turbulent flows. The irreversible entropy generation in this equation provides a more general description of subgrid scale (SGS) dissipation due to heat conduction, mass diffusion and viscosity effects. A new methodology is developed, termed the entropy filtered density function (En-FDF), to account for all individual entropy generation effects in turbulent flows. The En-FDF represents the joint probability density function of entropy, frequency, velocity and scalar fields within the SGS. An exact transport equation is developed for the En-FDF, which is modeled by a system of stochastic differential equations, incorporating the second law of thermodynamics. The modeled En-FDF transport equation is solved by a Lagrangian Monte Carlo method. The methodology is employed to simulate a turbulent mixing layer involving transport of passive scalars and entropy. Various modes of entropy generation are obtained from the En-FDF and analyzed. Predictions are assessed against data generated by direct numerical simulation (DNS). The En-FDF predictions are in good agreements with the DNS data.

Annual Research Briefs

NASA Technical Reports Server (NTRS)

Spinks, Debra (Compiler)

1997-01-01

This report contains the 1997 annual progress reports of the research fellows and students supported by the Center for Turbulence Research (CTR). Titles include: Invariant modeling in large-eddy simulation of turbulence; Validation of large-eddy simulation in a plain asymmetric diffuser; Progress in large-eddy simulation of trailing-edge turbulence and aeronautics; Resolution requirements in large-eddy simulations of shear flows; A general theory of discrete filtering for LES in complex geometry; On the use of discrete filters for large eddy simulation; Wall models in large eddy simulation of separated flow; Perspectives for ensemble average LES; Anisotropic grid-based formulas for subgrid-scale models; Some modeling requirements for wall models in large eddy simulation; Numerical simulation of 3D turbulent boundary layers using the V2F model; Accurate modeling of impinging jet heat transfer; Application of turbulence models to high-lift airfoils; Advances in structure-based turbulence modeling; Incorporating realistic chemistry into direct numerical simulations of turbulent non-premixed combustion; Effects of small-scale structure on turbulent mixing; Turbulent premixed combustion in the laminar flamelet and the thin reaction zone regime; Large eddy simulation of combustion instabilities in turbulent premixed burners; On the generation of vorticity at a free-surface; Active control of turbulent channel flow; A generalized framework for robust control in fluid mechanics; Combined immersed-boundary/B-spline methods for simulations of flow in complex geometries; and DNS of shock boundary-layer interaction - preliminary results for compression ramp flow.

Impacts of Mesoscale Eddies on the Vertical Nitrate Flux in the Gulf Stream Region

NASA Astrophysics Data System (ADS)

Zhang, Shuwen; Curchitser, Enrique N.; Kang, Dujuan; Stock, Charles A.; Dussin, Raphael

2018-01-01

The Gulf Stream (GS) region has intense mesoscale variability that can affect the supply of nutrients to the euphotic zone (Zeu). In this study, a recently developed high-resolution coupled physical-biological model is used to conduct a 25-year simulation in the Northwest Atlantic. The Reynolds decomposition method is applied to quantify the nitrate budget and shows that the mesoscale variability is important to the vertical nitrate supply over the GS region. The decomposition, however, cannot isolate eddy effects from those arising from other mesoscale phenomena. This limitation is addressed by analyzing a large sample of eddies detected and tracked from the 25-year simulation. The eddy composite structures indicate that positive nitrate anomalies within Zeu exist in both cyclonic eddies (CEs) and anticyclonic eddies (ACEs) over the GS region, and are even more pronounced in the ACEs. Our analysis further indicates that positive nitrate anomalies mostly originate from enhanced vertical advective flux rather than vertical turbulent diffusion. The eddy-wind interaction-induced Ekman pumping is very likely the mechanism driving the enhanced vertical motions and vertical nitrate transport within ACEs. This study suggests that the ACEs in GS region may play an important role in modulating the oceanic biogeochemical properties by fueling local biomass production through the persistent supply of nitrate.

Impacts of Diffuse Radiation on Light Use Efficiency across Terrestrial Ecosystems Based on Eddy Covariance Observation in China

PubMed Central

Huang, Kun; Wang, Shaoqiang; Zhou, Lei; Wang, Huimin; Zhang, Junhui; Yan, Junhua; Zhao, Liang; Wang, Yanfen; Shi, Peili

2014-01-01

Ecosystem light use efficiency (LUE) is a key factor of production models for gross primary production (GPP) predictions. Previous studies revealed that ecosystem LUE could be significantly enhanced by an increase on diffuse radiation. Under large spatial heterogeneity and increasing annual diffuse radiation in China, eddy covariance flux data at 6 sites across different ecosystems from 2003 to 2007 were used to investigate the impacts of diffuse radiation indicated by the cloudiness index (CI) on ecosystem LUE in grassland and forest ecosystems. Our results showed that the ecosystem LUE at the six sites was significantly correlated with the cloudiness variation (0.24≤R2≤0.85), especially at the Changbaishan temperate forest ecosystem (R2 = 0.85). Meanwhile, the CI values appeared more frequently between 0.8 and 1.0 in two subtropical forest ecosystems (Qianyanzhou and Dinghushan) and were much larger than those in temperate ecosystems. Besides, cloudiness thresholds which were favorable for enhancing ecosystem carbon sequestration existed at the three forest sites, respectively. Our research confirmed that the ecosystem LUE at the six sites in China was positively responsive to the diffuse radiation, and the cloudiness index could be used as an environmental regulator for LUE modeling in regional GPP prediction. PMID:25393629

Impacts of diffuse radiation on light use efficiency across terrestrial ecosystems based on Eddy covariance observation in China.

PubMed

Huang, Kun; Wang, Shaoqiang; Zhou, Lei; Wang, Huimin; Zhang, Junhui; Yan, Junhua; Zhao, Liang; Wang, Yanfen; Shi, Peili

2014-01-01

Ecosystem light use efficiency (LUE) is a key factor of production models for gross primary production (GPP) predictions. Previous studies revealed that ecosystem LUE could be significantly enhanced by an increase on diffuse radiation. Under large spatial heterogeneity and increasing annual diffuse radiation in China, eddy covariance flux data at 6 sites across different ecosystems from 2003 to 2007 were used to investigate the impacts of diffuse radiation indicated by the cloudiness index (CI) on ecosystem LUE in grassland and forest ecosystems. Our results showed that the ecosystem LUE at the six sites was significantly correlated with the cloudiness variation (0.24 ≤ R(2) ≤ 0.85), especially at the Changbaishan temperate forest ecosystem (R(2) = 0.85). Meanwhile, the CI values appeared more frequently between 0.8 and 1.0 in two subtropical forest ecosystems (Qianyanzhou and Dinghushan) and were much larger than those in temperate ecosystems. Besides, cloudiness thresholds which were favorable for enhancing ecosystem carbon sequestration existed at the three forest sites, respectively. Our research confirmed that the ecosystem LUE at the six sites in China was positively responsive to the diffuse radiation, and the cloudiness index could be used as an environmental regulator for LUE modeling in regional GPP prediction.

Prospective and retrospective high order eddy current mitigation for diffusion weighted echo planar imaging.

PubMed

Xu, Dan; Maier, Joseph K; King, Kevin F; Collick, Bruce D; Wu, Gaohong; Peters, Robert D; Hinks, R Scott

2013-11-01

The proposed method is aimed at reducing eddy current (EC) induced distortion in diffusion weighted echo planar imaging, without the need to perform further image coregistration between diffusion weighted and T2 images. These ECs typically have significant high order spatial components that cannot be compensated by preemphasis. High order ECs are first calibrated at the system level in a protocol independent fashion. The resulting amplitudes and time constants of high order ECs can then be used to calculate imaging protocol specific corrections. A combined prospective and retrospective approach is proposed to apply correction during data acquisition and image reconstruction. Various phantom, brain, body, and whole body diffusion weighted images with and without the proposed method are acquired. Significantly reduced image distortion and misregistration are consistently seen in images with the proposed method compared with images without. The proposed method is a powerful (e.g., effective at 48 cm field of view and 30 cm slice coverage) and flexible (e.g., compatible with other image enhancements and arbitrary scan plane) technique to correct high order ECs induced distortion and misregistration for various diffusion weighted echo planar imaging applications, without the need for further image post processing, protocol dependent prescan, or sacrifice in signal-to-noise ratio. Copyright © 2013 Wiley Periodicals, Inc.

A low-dissipation monotonicity-preserving scheme for turbulent flows in hydraulic turbines

NASA Astrophysics Data System (ADS)

Yang, L.; Nadarajah, S.

2016-11-01

The objective of this work is to improve the inherent dissipation of the numerical schemes under the framework of a Reynolds-averaged Navier-Stokes (RANS) simulation. The governing equations are solved by the finite volume method with the k-ω SST turbulence model. Instead of the van Albada limiter, a novel eddy-preserving limiter is employed in the MUSCL reconstructions to minimize the dissipation of the vortex. The eddy-preserving procedure inactivates the van Albada limiter in the swirl plane and reduces the artificial dissipation to better preserve vortical flow structures. Steady and unsteady simulations of turbulent flows in a straight channel and a straight asymmetric diffuser are demonstrated. Profiles of velocity, Reynolds shear stress and turbulent kinetic energy are presented and compared against large eddy simulation (LES) and/or experimental data. Finally, comparisons are made to demonstrate the capability of the eddy-preserving limiter scheme.

A Kolmogorov-Brutsaert structure function model for evaporation into a turbulent atmosphere

NASA Astrophysics Data System (ADS)

Katul, Gabriel; Liu, Heping

2017-05-01

In 1965, Brutsaert proposed a model that predicted mean evaporation rate E¯ from rough surfaces to scale with the 3/4 power law of the friction velocity (u∗) and the square-root of molecular diffusivity (Dm) for water vapor. In arriving at these results, a number of assumptions were made regarding the surface renewal rate describing the contact durations between eddies and the evaporating surface, the diffusional mass process from the surface into eddies, and the cascade of turbulent kinetic energy sustaining the eddy renewal process itself. The working hypothesis explored here is that E¯˜Dmu∗3/4 is a direct outcome of the Kolmogorov scaling for inertial subrange eddies modified to include viscous cutoff thereby bypassing the need for a surface renewal assumption. It is demonstrated that Brutsaert's model for E¯ may be more general than its original derivation implied.

Hadley cell dynamics of a cold and virtually dry Snowball Earth atmosphere

NASA Astrophysics Data System (ADS)

Voigt, Aiko; Held, Isaac; Marotzke, Jochem

2010-05-01

We use the full-physics atmospheric general circulation model ECHAM5 to investigate a cold and virtually dry Snowball Earth atmosphere that results from specifying sea ice as the surface boundary condition everywhere, corresponding to a frozen aquaplanet, while keeping total solar irradiance at its present-day value of 1365 Wm-2. The aim of this study is the investigation of the zonal-mean circulation of a Snowball Earth atmosphere, which, due to missing moisture, might constitute an ideal though yet unexplored testbed for theories of atmospheric dynamics. To ease comparison with theories, incoming solar insolation follows permanent equinox conditions with disabled diurnal cycle. The meridional circulation consists of a thermally direct cell extending from the equator to 45 N/S with ascent in the equatorial region, and a weak thermally indirect cell with descent between 45 and 65 N/S and ascent in the polar region. The former cell corresponds to the present-day Earth's Hadley cell, while the latter can be viewed as an eddy-driven Ferrell cell; the present-day Earth's direct polar cell is missing. The Hadley cell itself is subdivided into a vigorous cell confined to the troposphere and a weak deep cell reaching well into the stratosphere. The dynamics of the vigorous Snowball Earth Hadley cell differ substantially from the dynamics of the present-day Hadley cell. The zonal momentum balance shows that in the poleward branch of the vigorous Hadley cell, mean flow meridional advection of absolute vorticity is not only balanced by eddy momentum flux convergence but also by vertical diffusion. Inside the poleward branch, eddies are more important in the upper part and vertical diffusion is more important in the lower part. Vertical diffusion also contributes to the meridional momentum balance as it decelerates the vigorous Hadley cell by downgradient momentum mixing between its poleward and equatorward branch. Zonal winds, therefore, are not in thermal wind balance in the vigorous Hadley cell. Suppressing vertical momentum diffusion above 870 hPa results in a doubling of the vigorous Hadley cell strength. Simulations where we only suppress either vertical diffusion of zonal or meridional momentum show that this doubling can be understood from the decelerating effect of vertical diffusion in the meridional momentum balance. Comparing our simulations with theories, we conclude that neither the axisymmetric Hadley cell model of Held & Hou (1980) nor the eddy-permitting model of T. Schneider et al. (2005, 2006, 2008) are applicable to a Snowball Earth atmosphere since both assume an inviscid upper Hadley cell branch.

Time evolution of the eddy viscosity in two-dimensional navier-stokes flow

PubMed

Chaves; Gama

2000-02-01

The time evolution of the eddy viscosity associated with an unforced two-dimensional incompressible Navier-Stokes flow is analyzed by direct numerical simulation. The initial condition is such that the eddy viscosity is isotropic and negative. It is shown by concrete examples that the Navier-Stokes dynamics stabilizes negative eddy viscosity effects. In other words, this dynamics moves monotonically the initial negative eddy viscosity to positive values before relaxation due to viscous term occurs.

Finite-element numerical modeling of atmospheric turbulent boundary layer

NASA Technical Reports Server (NTRS)

Lee, H. N.; Kao, S. K.

1979-01-01

A dynamic turbulent boundary-layer model in the neutral atmosphere is constructed, using a dynamic turbulent equation of the eddy viscosity coefficient for momentum derived from the relationship among the turbulent dissipation rate, the turbulent kinetic energy and the eddy viscosity coefficient, with aid of the turbulent second-order closure scheme. A finite-element technique was used for the numerical integration. In preliminary results, the behavior of the neutral planetary boundary layer agrees well with the available data and with the existing elaborate turbulent models, using a finite-difference scheme. The proposed dynamic formulation of the eddy viscosity coefficient for momentum is particularly attractive and can provide a viable alternative approach to study atmospheric turbulence, diffusion and air pollution.

Small-scale structure and turbulence observed in MAP/WINE)

NASA Technical Reports Server (NTRS)

Blix, T. A.

1989-01-01

During MAP/WINE small scale structure and turbulence in the mesosphere and lower thermosphere was studied in situ by rocket-borne instruments as well as from the ground by remote sensing techniques. The eight salvoes launched during the campaign resulted in a wealth of information on the dynamical structure of these regions. The experimental results are reviewed and their interpretation is discussed in terms of gravity waves and turbulence. It is shown that eddy diffusion coefficients and turbulent energy dissipation rates may be derived from the in situ measurements in a consistent manner. The observations are also shown to be consistent with the hypothesis that turbulence can be created by a process of gravity wave saturation.

A numerical method for the solution of internal pipe/channel flows in laminar or turbulent motion

NASA Astrophysics Data System (ADS)

Lourenco, L.; Essers, J. A.

1981-11-01

A computer program which is useful in the solution of problems of internal turbulent or laminar flow without recirculation is described. The flow is treated in terms of parabolic boundary layer differential equations. The eddy diffusivity concept is used to model turbulent stresses. Two turbulent models are available: the Prandtl mixing length model and the Nee-Kovasznay model for the effective viscosity. Fluid is considered incompressible, but little program modification is needed to treat compressible flows. Initial conditions are prescribed as well as the boundary conditions. The differencing scheme employed is fully implicit for the dependent variables. This allows the use of relatively large forward steps without stability problems.

Convection with a simple chemically reactive passive scalar

NASA Astrophysics Data System (ADS)

Herring, J. R.; Wyngaard, J. C.

Convection between horizontal stress-free perfectly conducting plates is examined in the turbulent regime for air. Results are presented for an additional scalar undergoing simple linear decay. We discuss qualitative aspects of the flow in terms of spectral and three-dimensional contour maps of the velocity and scalar fields. The horizontal mean profiles of scalar gradients and fluxes agree rather well with simple mixing-length concepts. Further, the mean profiles for a range of the destruction-rate parameter are shown to be nearly completely characterized by the boundary fluxes. Finally, we shall use the present numerical data as a basis for exploring a generalization of eddy-diffusion concepts so as to properly incorporate non-local effects.

Analysis of phase error effects in multishot diffusion-prepared turbo spin echo imaging

PubMed Central

Cervantes, Barbara; Kooijman, Hendrik; Karampinos, Dimitrios C.

2017-01-01

Background To characterize the effect of phase errors on the magnitude and the phase of the diffusion-weighted (DW) signal acquired with diffusion-prepared turbo spin echo (dprep-TSE) sequences. Methods Motion and eddy currents were identified as the main sources of phase errors. An analytical expression for the effect of phase errors on the acquired signal was derived and verified using Bloch simulations, phantom, and in vivo experiments. Results Simulations and experiments showed that phase errors during the diffusion preparation cause both magnitude and phase modulation on the acquired data. When motion-induced phase error (MiPe) is accounted for (e.g., with motion-compensated diffusion encoding), the signal magnitude modulation due to the leftover eddy-current-induced phase error cannot be eliminated by the conventional phase cycling and sum-of-squares (SOS) method. By employing magnitude stabilizers, the phase-error-induced magnitude modulation, regardless of its cause, was removed but the phase modulation remained. The in vivo comparison between pulsed gradient and flow-compensated diffusion preparations showed that MiPe needed to be addressed in multi-shot dprep-TSE acquisitions employing magnitude stabilizers. Conclusions A comprehensive analysis of phase errors in dprep-TSE sequences showed that magnitude stabilizers are mandatory in removing the phase error induced magnitude modulation. Additionally, when multi-shot dprep-TSE is employed the inconsistent signal phase modulation across shots has to be resolved before shot-combination is performed. PMID:28516049

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.

Quantitative computational infrared imaging of buoyant diffusion flames

NASA Astrophysics Data System (ADS)

Newale, Ashish S.

Studies of infrared radiation from turbulent buoyant diffusion flames impinging on structural elements have applications to the development of fire models. A numerical and experimental study of radiation from buoyant diffusion flames with and without impingement on a flat plate is reported. Quantitative images of the radiation intensity from the flames are acquired using a high speed infrared camera. Large eddy simulations are performed using fire dynamics simulator (FDS version 6). The species concentrations and temperature from the simulations are used in conjunction with a narrow-band radiation model (RADCAL) to solve the radiative transfer equation. The computed infrared radiation intensities rendered in the form of images and compared with the measurements. The measured and computed radiation intensities reveal necking and bulging with a characteristic frequency of 7.1 Hz which is in agreement with previous empirical correlations. The results demonstrate the effects of stagnation point boundary layer on the upstream buoyant shear layer. The coupling between these two shear layers presents a model problem for sub-grid scale modeling necessary for future large eddy simulations.

Spectral analysis of large-eddy advection in ET from eddy covariance towers and a large weighting lysimeter

USDA-ARS?s Scientific Manuscript database

Evapotranspiration was continuously measured by an array of eddy covariance systems and large weighting lysimeter in a cotton field in Bushland, Texas. The advective divergence from both horizontal and vertical directions were measured through profile measurements above canopy. All storage terms wer...

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.

Eddy covariance measurement of CO2 flux to the atmosphere from a area of high volcanogenic emissions, Mammoth Mountain, California

USGS Publications Warehouse

Anderson, D.E.; Farrar, C.D.

2001-01-01

Three pilot studies were performed to assess application of the eddy covariance micrometeorological method in the measurement of carbon dioxide (CO2) flux of volcanic origin. The selected study area is one of high diffuse CO2 emission on Mammoth Mountain, CA. Because terrain and source characteristics make this a complex setting for this type of measurement, added consideration was given to source area and upwind fetch. Footprint analysis suggests that the eddy covariance measurements were representative of an upwind elliptical source area (3.8 ?? 103 m2) which can vary with mean wind direction, surface roughness, and atmospheric stability. CO2 flux averaged 8-16 mg m-2 s-1 (0.7-1.4 kg m-2 day-1). Eddy covariance measurements of flux were compared with surface chamber measurements made in separate studies [Geophys. Res. Lett. 25 (1998a) 1947; EOS Trans. 79 (1998) F941.] and were found to be similar. ?? 2001 Elsevier Science B.V. All rights reserved.

Anisotropic shear dispersion parameterization for ocean eddy transport

NASA Astrophysics Data System (ADS)

Reckinger, Scott; Fox-Kemper, Baylor

2015-11-01

The effects of mesoscale eddies are universally treated isotropically in global ocean general circulation models. However, observations and simulations demonstrate that the mesoscale processes that the parameterization is intended to represent, such as shear dispersion, are typified by strong anisotropy. We extend the Gent-McWilliams/Redi mesoscale eddy parameterization to include anisotropy and test the effects of varying levels of anisotropy in 1-degree Community Earth System Model (CESM) simulations. Anisotropy has many effects on the simulated climate, including a reduction of temperature and salinity biases, a deepening of the southern ocean mixed-layer depth, impacts on the meridional overturning circulation and ocean energy and tracer uptake, and improved ventilation of biogeochemical tracers, particularly in oxygen minimum zones. A process-based parameterization to approximate the effects of unresolved shear dispersion is also used to set the strength and direction of anisotropy. The shear dispersion parameterization is similar to drifter observations in spatial distribution of diffusivity and high-resolution model diagnosis in the distribution of eddy flux orientation.

Excitation of turbulence by density waves

NASA Technical Reports Server (NTRS)

Tichen, C. M.

1985-01-01

A nonlinear system describes the microdynamical state of turbulence that is excited by density waves. It consists of an equation of propagation and a master equation. A group-scaling generates the scaled equations of many interacting groups of distribution functions. The two leading groups govern the transport processes of evolution and eddy diffusivity. The remaining sub-groups represent the relaxation for the approach of diffusivity to equilibrium. In strong turbulence, the sub-groups disperse themselves and the ensemble acts like a medium that offers an effective damping to close the hierarchy. The kinetic equation of turbulence is derived. It calculates the eddy viscosity and identifies the effective damping of the assumed medium self-consistently. It formulates the coupling mechanism for the intensification of the turbulent energy at the expense of the wave energy, and the transfer mechanism for the cascade. The spectra of velocity and density fluctuations find the power law k sup-2 and k sup-4, respectively.

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.

Gas exchange across the air-sea interface

NASA Astrophysics Data System (ADS)

Hasse, L.; Liss, P. S.

1980-10-01

The physics of gas exchange at the air-sea interface are reviewed. In order to describe the transfer of gases in the liquid near the boundary, a molecular plus eddy diffusivity concept is used, which has been found useful for smooth flow over solid surfaces. From consideration of the boundary conditions, a similar dependence of eddy diffusivity on distance from the interface can be derived for the flow beneath a gas/liquid interface, at least in the absence of waves. The influence of waves is then discussed. It is evident from scale considerations that the effect of gravity waves is small. It is known from wind tunnel work that capillary waves enhance gas transfer considerably. The existing hypotheses are apparently not sufficient to explain the observations. Examination of field data is even more frustrating since the data do not show the expected increase of gas exchange with wind speed.

An Examination of the Evolution of Radiation and Advection Fogs

DTIC Science & Technology

1993-01-01

and fog diagnostic and prediction models have developed in sophistication so that they can reproduce fairly accurate one- or two-dimensional...occurred only by molecular diffusion near the interface created between the species during the mixing process. The rate of homogenization is minimal until...of excess vapor by molecular diffusion at the interfaces of nearly saturated air mixing in eddies is faster than the relaxation time of droplet

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.

Budget of Turbulent Kinetic Energy in a Shock Wave Boundary-Layer Interaction

NASA Technical Reports Server (NTRS)

Vyas, Manan A.; Waindim, Mbu; Gaitonde, Datta V.

2016-01-01

Implicit large-eddy simulation (ILES) of a shock wave/boundary-layer interaction (SBLI) was performed. Quantities present in the exact equation of the turbulent kinetic energy transport were accumulated and used to calculate terms like production, dissipation, molecular diffusion, and turbulent transport. The present results for a turbulent boundary layer were validated by comparison with direct numerical simulation data. It was found that a longer development domain was necessary for the boundary layer to reach an equilibrium state and a finer mesh resolution would improve the predictions. In spite of these findings, trends of the present budget match closely with that of the direct numerical simulation. Budgets for the SBLI region are presented at key axial stations. These budgets showed interesting dynamics as the incoming boundary layer transforms and the terms of the turbulent kinetic energy budget change behavior within the interaction region.

Estimates of advection and diffusion in the Potomac estuary

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

Elliott, A.J.

1976-01-01

A two-layered dispersion model, suitable for application to partially-mixed estuaries, has been developed to provide hydrological interpretation of the results of biological sampling. The model includes horizontal and vertical advection plus both horizontal and vertical diffusion. A pseudo-geostrophic method, which includes a damping factor to account for internal eddy friction, is used to estimate the horizontal advective fluxes and the results are compared with field observations. A salt balance model is then used to estimate the effective diffusivities in the Potomac estuary during the Spring of 1974.

The vertical distribution of tropospheric ammonia

NASA Technical Reports Server (NTRS)

Levine, J. S.; Hoell, J. M.; Augustsson, T. R.

1980-01-01

A one-dimensional tropospheric photochemical model is used to simulate measured profiles of NH3 obtained with the Infrared Heterodyne Radiometer. The relative roles of homogeneous loss, heterogeneous loss, and vertical eddy transport are discussed in terms of selecting parameters which best fit the measurements. The best fit was obtained for a vertical eddy diffusion coefficient of 200,000/sq cm per sec or greater (corresponding to a characteristic vertical transport time in excess of about 35 days), and a characteristic heterogeneous loss time in excess of 10 days. The characteristic homogeneous chemical loss time was found to be about 40 days at the surface and decreased to about 180 days at 10 km, and not very sensitive to model chemical perturbations. Increased ground-level concentrations of NH3 to about 10 ppb, compared to background surface concentrations of about 1 ppb, were measured several weeks after application of ammonium nitrate fertilizer. This suggests that the volatilization of ammonium nitrate fertilizer is rapid, and an important source of NH3. Because of the characteristic times for the loss mechanisms, synoptic time-scale phenomena may play an important role in determining the tropospheric distribution of NH3 concentrations.

Large-eddy simulation of sand dune morphodynamics

NASA Astrophysics Data System (ADS)

Khosronejad, Ali; Sotiropoulos, Fotis; St. Anthony Falls Laboratory, University of Minnesota Team

2015-11-01

Sand dunes are natural features that form under complex interaction between turbulent flow and bed morphodynamics. We employ a fully-coupled 3D numerical model (Khosronejad and Sotiropoulos, 2014, Journal of Fluid Mechanics, 753:150-216) to perform high-resolution large-eddy simulations of turbulence and bed morphodynamics in a laboratory scale mobile-bed channel to investigate initiation, evolution and quasi-equilibrium of sand dunes (Venditti and Church, 2005, J. Geophysical Research, 110:F01009). We employ a curvilinear immersed boundary method along with convection-diffusion and bed-morphodynamics modules to simulate the suspended sediment and the bed-load transports respectively. The coupled simulation were carried out on a grid with more than 100 million grid nodes and simulated about 3 hours of physical time of dune evolution. The simulations provide the first complete description of sand dune formation and long-term evolution. The geometric characteristics of the simulated dunes are shown to be in excellent agreement with observed data obtained across a broad range of scales. This work was supported by NSF Grants EAR-0120914 (as part of the National Center for Earth-Surface Dynamics). Computational resources were provided by the University of Minnesota Supercomputing Institute.

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

Wong, May Wai San; Ovchinnikov, Mikhail; Wang, Minghuai

Potential ways of parameterizing vertical turbulent fluxes of hydrometeors are examined using a high-resolution cloud-resolving model. The cloud-resolving model uses the Morrison microphysics scheme, which contains prognostic variables for rain, graupel, ice, and snow. A benchmark simulation with a horizontal grid spacing of 250 m of a deep convection case carried out to evaluate three different ways of parameterizing the turbulent vertical fluxes of hydrometeors: an eddy-diffusion approximation, a quadrant-based decomposition, and a scaling method that accounts for within-quadrant (subplume) correlations. Results show that the down-gradient nature of the eddy-diffusion approximation tends to transport mass away from concentrated regions, whereasmore » the benchmark simulation indicates that the vertical transport tends to transport mass from below the level of maximum to aloft. Unlike the eddy-diffusion approach, the quadri-modal decomposition is able to capture the signs of the flux gradient but underestimates the magnitudes. The scaling approach is shown to perform the best by accounting for within-quadrant correlations, and improves the results for all hydrometeors except for snow. A sensitivity study is performed to examine how vertical transport may affect the microphysics of the hydrometeors. The vertical transport of each hydrometeor type is artificially suppressed in each test. Results from the sensitivity tests show that cloud-droplet-related processes are most sensitive to suppressed rain or graupel transport. In particular, suppressing rain or graupel transport has a strong impact on the production of snow and ice aloft. Lastly, a viable subgrid-scale hydrometeor transport scheme in an assumed probability density function parameterization is discussed.« less

Seismic Oceanography in the Tyrrhenian Sea: Thermohaline Staircases, Eddies, and Internal Waves

NASA Astrophysics Data System (ADS)

Buffett, G. G.; Krahmann, G.; Klaeschen, D.; Schroeder, K.; Sallarès, V.; Papenberg, C.; Ranero, C. R.; Zitellini, N.

2017-11-01

We use seismic oceanography to document and analyze oceanic thermohaline fine structure across the Tyrrhenian Sea. Multichannel seismic (MCS) reflection data were acquired during the MEDiterranean OCcidental survey in April-May 2010. We deployed along-track expendable bathythermograph probes simultaneous with MCS acquisition. At nearby locations we gathered conductivity-temperature-depth data. An autonomous glider survey added in situ measurements of oceanic properties. The seismic reflectivity clearly delineates thermohaline fine structure in the upper 2,000 m of the water column, indicating the interfaces between Atlantic Water/Winter Intermediate Water, Levantine Intermediate Water, and Tyrrhenian Deep Water. We observe the Northern Tyrrhenian Anticyclone, a near-surface mesoscale eddy, plus laterally and vertically extensive thermohaline staircases. Using MCS, we are able to fully image the anticyclone to a depth of 800 m and to confirm the horizontal continuity of the thermohaline staircases of more than 200 km. The staircases show the clearest step-like gradients in the center of the basin while they become more diffuse toward the periphery and bottom, where impedance gradients become too small to be detected by MCS. We quantify the internal wave field and find it to be weak in the region of the eddy and in the center of the staircases, while it is stronger near the coastlines. Our results indicate this is because of the influence of the boundary currents, which disrupt the formation of staircases by preventing diffusive convection. In the interior of the basin, the staircases are clearer and the internal wave field weaker, suggesting that other mixing processes such as double diffusion prevail.

Time variable eddy mixing in the global Sea Surface Salinity maxima

NASA Astrophysics Data System (ADS)

Busecke, J. J. M.; Abernathey, R.; Gordon, A. L.

2016-12-01

Lateral mixing by mesoscale eddies is widely recognized as a crucial mechanism for the global ocean circulation and the associated heat/salt/tracer transports. The Salinity in the Upper Ocean Processes Study (SPURS) confirmed the importance of eddy mixing for the surface salinity fields even in the center of the subtropical gyre of the North Atlantic. We focus on the global salinity maxima due to their role as indicators for global changes in the hydrological cycle as well as providing the source water masses for the shallow overturning circulation. We introduce a novel approach to estimate the contribution of eddy mixing to the global sea surface salinity maxima. Using a global 2D tracer experiments in a 1/10 degree MITgcm setup driven by observed surface velocities, we analyze the effect of eddy mixing using a water mass framework, thus focussing on the diffusive flux across surface isohalines. This enables us to diagnose temporal variability on seasonal to inter annual time scales, revealing regional differences in the mechanism causing temporal variability.Sensitivity experiments with various salinity backgrounds reveal robust inter annual variability caused by changes in the surface velocity fields potentially forced by large scale climate.

Modeling Postconvective Submesoscale Coherent Vortices in the Northwestern Mediterranean Sea

NASA Astrophysics Data System (ADS)

Damien, P.; Bosse, A.; Testor, P.; Marsaleix, P.; Estournel, C.

2017-12-01

For the first time, the formation of submesoscale coherent vortices (SCVs) during intermediate and deep convection events is documented in a realistic high-resolution (1 km) numerical simulation of the oceanic circulation in the northwestern Mediterranean Sea. Winter intermediate and deep convection leads to the formation of anticyclonic and cyclonic eddies with lifetimes exceeding 1 year. By focusing on three typical eddies, the main characteristics of such vortices are discussed. The anticyclonic eddies are typical of SCVs observed in deep convection areas so far. They are characterized by a small radius (˜6.5 km) and orbital peak velocities of about 7 cm/s located at great depth (˜1500 m) or intermediate depth (˜500 m). The cyclonic vortices show very similar characteristics, such as a high Rossby number (˜0.4), but with surface-intensified structures. The long lifetimes of both anticyclonic and cyclonic eddies reflect very slow diffusive processes between their core and their surroundings and a strong resistance to external perturbations. These long-lived eddies are found to participate in the spreading of a significant portion (from 15 to 35%) of the convected waters in the Gulf of Lions and contribute to the ventilation of the deep basin.

Increasing of eddy activity in the northeastern Pacific during 1993-2011

NASA Astrophysics Data System (ADS)

Ding, M.; Lin, P.; Liu, H.; Chai, F.

2017-12-01

We study the long-term behaviors of eddy activity in the northeastern Pacific (NEP) and the dynamic mechanism behind them based on the 3rd version of the mesoscale eddy trajectories dataset released by Chelton et al. (2013) combined with other observation and reanalysis datasets. Both the eddy kinetic energy (EKE) and eddy occurrence number (EON) present prominent increasing trends, with inter-annual and decadal variabilities northeast of the Hawaii-Emperor seamounts. The increasing trend of the EON is mainly due to prolongation of the eddy lifetime associated with the eddy intensification, particularly for anticyclonic eddies (AEs). Weakened surface winds tend to prolong the eddy lifetimes, as the eddy attenuation time scale is inversely proportional to the wind speed. The enhanced anticyclonic wind stress curl (WSC) anomalies inject more energy into the AE over the study region and provide a more suitable environment for AEs growth. The decadal climate modes, such as the Pacific decadal oscillation (PDO) and the North Pacific gyre oscillation (NPGO), may also modulate eddy activities in the NEP by exerting fluctuations in the surface wind system.

Vertically migrating swimmers generate aggregation-scale eddies in a stratified column.

PubMed

Houghton, Isabel A; Koseff, Jeffrey R; Monismith, Stephen G; Dabiri, John O

2018-04-01

Biologically generated turbulence has been proposed as an important contributor to nutrient transport and ocean mixing 1-3 . However, to produce non-negligible transport and mixing, such turbulence must produce eddies at scales comparable to the length scales of stratification in the ocean. It has previously been argued that biologically generated turbulence is limited to the scale of the individual animals involved 4 , which would make turbulence created by highly abundant centimetre-scale zooplankton such as krill irrelevant to ocean mixing. Their small size notwithstanding, zooplankton form dense aggregations tens of metres in vertical extent as they undergo diurnal vertical migration over hundreds of metres 3,5,6 . This behaviour potentially introduces additional length scales-such as the scale of the aggregation-that are of relevance to animal interactions with the surrounding water column. Here we show that the collective vertical migration of centimetre-scale swimmers-as represented by the brine shrimp Artemia salina-generates aggregation-scale eddies that mix a stable density stratification, resulting in an effective turbulent diffusivity up to three orders of magnitude larger than the molecular diffusivity of salt. These observed large-scale mixing eddies are the result of flow in the wakes of the individual organisms coalescing to form a large-scale downward jet during upward swimming, even in the presence of a strong density stratification relative to typical values observed in the ocean. The results illustrate the potential for marine zooplankton to considerably alter the physical and biogeochemical structure of the water column, with potentially widespread effects owing to their high abundance in climatically important regions of the ocean 7 .

Subgrid-scale models for large-eddy simulation of rotating turbulent channel flows

NASA Astrophysics Data System (ADS)

Silvis, Maurits H.; Bae, Hyunji Jane; Trias, F. Xavier; Abkar, Mahdi; Moin, Parviz; Verstappen, Roel

2017-11-01

We aim to design subgrid-scale models for large-eddy simulation of rotating turbulent flows. Rotating turbulent flows form a challenging test case for large-eddy simulation due to the presence of the Coriolis force. The Coriolis force conserves the total kinetic energy while transporting it from small to large scales of motion, leading to the formation of large-scale anisotropic flow structures. The Coriolis force may also cause partial flow laminarization and the occurrence of turbulent bursts. Many subgrid-scale models for large-eddy simulation are, however, primarily designed to parametrize the dissipative nature of turbulent flows, ignoring the specific characteristics of transport processes. We, therefore, propose a new subgrid-scale model that, in addition to the usual dissipative eddy viscosity term, contains a nondissipative nonlinear model term designed to capture transport processes, such as those due to rotation. We show that the addition of this nonlinear model term leads to improved predictions of the energy spectra of rotating homogeneous isotropic turbulence as well as of the Reynolds stress anisotropy in spanwise-rotating plane-channel flows. This work is financed by the Netherlands Organisation for Scientific Research (NWO) under Project Number 613.001.212.

Concerning the measurement of atmospheric trace gas fluxes with open- and closed-path eddy covariance systems: The density terms and spectral attenuation [Chapter 7

Treesearch

W. J. Massman

2004-01-01

Atmospheric trace gas fluxes measured with an eddy covariance sensor that detects a constituent's density fluctuations within the in situ air need to include terms resulting from concurrent heat and moisture fluxes, the so called 'density' or 'WPL corrections' (Webb et al. 1980). The theory behind these additional terms is well established. But...

Renormalization group estimates of transport coefficients in the advection of a passive scalar by incompressible turbulence

NASA Technical Reports Server (NTRS)

Zhou, YE; Vahala, George

1993-01-01

The advection of a passive scalar by incompressible turbulence is considered using recursive renormalization group procedures in the differential sub grid shell thickness limit. It is shown explicitly that the higher order nonlinearities induced by the recursive renormalization group procedure preserve Galilean invariance. Differential equations, valid for the entire resolvable wave number k range, are determined for the eddy viscosity and eddy diffusivity coefficients, and it is shown that higher order nonlinearities do not contribute as k goes to 0, but have an essential role as k goes to k(sub c) the cutoff wave number separating the resolvable scales from the sub grid scales. The recursive renormalization transport coefficients and the associated eddy Prandtl number are in good agreement with the k-dependent transport coefficients derived from closure theories and experiments.

An analysis and implications of alternative methods of deriving the density (WPL) terms for eddy covariance flux measurements

Treesearch

W. J. Massman; J. -P. Tuovinen

2006-01-01

We explore some of the underlying assumptions used to derive the density or WPL terms (Webb et al. (1980) Quart J RoyMeteorol Soc 106:85-100) required for estimating the surface exchange fluxes by eddy covariance. As part of this effort we recast the origin of the density terms as an assumption regarding the density fluctuations rather than as a (dry air) flux...

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.

Lagrangian analysis by clustering. An example in the Nordic Seas.

NASA Astrophysics Data System (ADS)

Koszalka, Inga; Lacasce, Joseph H.

2010-05-01

We propose a new method for obtaining average velocities and eddy diffusivities from Lagrangian data. Rather than grouping the drifter-derived velocities in uniform geographical bins, as is commonly done, we group a specified number of nearest-neighbor velocities. This is done via a clustering algorithm operating on the instantaneous positions of the drifters. Thus it is the data distribution itself which determines the positions of the averages and the areal extent of the clusters. A major advantage is that because the number of members is essentially the same for all clusters, the statistical accuracy is more uniform than with geographical bins. We illustrate the technique using synthetic data from a stochastic model, employing a realistic mean flow. The latter is an accurate representation of the surface currents in the Nordic Seas and is strongly inhomogeneous in space. We use the clustering algorithm to extract the mean velocities and diffusivities (both of which are known from the stochastic model). We also compare the results to those obtained with fixed geographical bins. Clustering is more successful at capturing spatial variability of the mean flow and also improves convergence in the eddy diffusivity estimates. We discuss both the future prospects and shortcomings of the new method.

Heat uptake in the Southern Ocean in a warmer, windier world: a process-based analysis using an AOGCM with an eddy-permitting ocean

NASA Astrophysics Data System (ADS)

Kuhlbrodt, T.; Gregory, J. M.

2016-02-01

About 90% of the anthropogenic increase in heat stored in the climate system is found the oceans. Therefore it is relevant to understand the details of ocean heat uptake. Here we present a detailed, process-based analysis of ocean heat uptake (OHU) processes in HiGEM1.2, an atmosphere-ocean general circulation model (AOGCM) with an eddy-permitting ocean component of 1/3° resolution. Similarly to various other models, HiGEM1.2 shows that the global heat budget is dominated by a downward advection of heat compensated by upward isopycnal diffusion. This upward isopycnal diffusion of heat is located mostly in the Southern Ocean (Fig. 1a).We compare the responses to a 4xCO2 forcing and an enhancement of the windstress forcing in the Southern Ocean. In line with the CMIP5 models, HiGEM1.2 shows a band of strong OHU in the mid-latitude Southern Ocean in the 4xCO2 run, which is mostly advective. By contrast, in the high-latitude Southern Ocean regions it is the suppression of convection that leads to OHU (Fig. 1b). In the enhanced windstress run, convection is strengthened at high Southern latitudes (Fig. 1c), leading to heat loss, while the magnitude of the OHU in the Southern mid-latitudes is very similar to the 4xCO2 results. Remarkably, there is only very small global OHU in the enhanced windstress run. The wind stress forcing just leads to a redistribution of heat. We relate the ocean changes at high southern latitudes to the effect of climate change on the Antarctic Circumpolar Current (ACC). It weakens in the 4xCO2 run and strengthens in the wind stress run. The weakening is due to a narrowing of the ACC, caused by an expansion of the Weddell Gyre, and a flattening of the isopycnals, which are explained by a combination of the wind stress forcing and increased precipitation. The presentation will also try to clarify the definitions of terms like "advective", "diffusive" and "eddy-induced" when used for observed and modelled (at various resolutions) ocean heat uptake processes. Fig. 1: Horizontally averaged temperature tendency diagnostics for the high-latitude Southern Ocean, for (a) the control run, (b) the 4xCO2 anomalies and (c) the windstress anomalies. Both axes are scaled according to a power law. "VM"- vertical mixing, which includes convection ("conv").

A LES-CMC formulation for premixed flames including differential diffusion

NASA Astrophysics Data System (ADS)

Farrace, Daniele; Chung, Kyoungseoun; Bolla, Michele; Wright, Yuri M.; Boulouchos, Konstantinos; Mastorakos, Epaminondas

2018-05-01

A finite volume large eddy simulation-conditional moment closure (LES-CMC) numerical framework for premixed combustion developed in a previous studyhas been extended to account for differential diffusion. The non-unity Lewis number CMC transport equation has an additional convective term in sample space proportional to the conditional diffusion of the progress variable, that in turn accounts for diffusion normal to the flame front and curvature-induced effects. Planar laminar simulations are first performed using a spatially homogeneous non-unity Lewis number CMC formulation and validated against physical-space fully resolved reference solutions. The same CMC formulation is subsequently used to numerically investigate the effects of curvature for laminar flames having different effective Lewis numbers: a lean methane-air flame with Leeff = 0.99 and a lean hydrogen-air flame with Leeff = 0.33. Results suggest that curvature does not affect the conditional heat release if the effective Lewis number tends to unity, so that curvature-induced transport may be neglected. Finally, the effect of turbulence on the flame structure is qualitatively analysed using LES-CMC simulations with and without differential diffusion for a turbulent premixed bluff body methane-air flame exhibiting local extinction behaviour. Overall, both the unity and the non-unity computations predict the characteristic M-shaped flame observed experimentally, although some minor differences are identified. The findings suggest that for the high Karlovitz number (from 1 to 10) flame considered, turbulent mixing within the flame weakens the differential transport contribution by reducing the conditional scalar dissipation rate and accordingly the conditional diffusion of the progress variable.

A deliberate tracer experiment in Santa Monica Basin. [for ocean density strata diffusion

NASA Technical Reports Server (NTRS)

Ledwell, J. R.; Broecker, W. S.; Watson, A. J.

1986-01-01

A tracer technique was developed for measurements of diffusion across oceanic density strata using SF6 and perfluorodecalin (PFD) tracers in the Santa Monica Basin. Fifty days after injection, the tracers were found to have mixed along the isopycnal surface to nearly every part of the basin. The diapycnal spreading of the tracer distributions yielded an apparent eddy diffusivity of 0.33 + or - 0.08 sq cm/s at the ambient density gradient of 4.0 + or - 0.5 x 10 to the -9th g/cm to the 4th.

A Fractional PDE Approach to Turbulent Mixing; Part II: Numerical Simulation

NASA Astrophysics Data System (ADS)

Samiee, Mehdi; Zayernouri, Mohsen

2016-11-01

We propose a generalizing fractional order transport model of advection-diffusion kind with fractional time- and space-derivatives, governing the evolution of passive scalar turbulence. This approach allows one to incorporate the nonlocal and memory effects in the underlying anomalous diffusion i.e., sub-to-standard diffusion to model the trapping of particles inside the eddied, and super-diffusion associated with the sudden jumps of particles from one coherent region to another. For this nonlocal model, we develop a high order numerical (spectral) method in addition to a fast solver, examined in the context of some canonical problems. PhD student, Department of Mechanical Engineering, & Department Computational Mathematics, Science, and Engineering.

Annual Research Briefs, 1998

NASA Technical Reports Server (NTRS)

Spinks, Debra (Compiler)

1998-01-01

The topics contained in this progress report are direct numerical simulation of turbulent non-premixed combustion with realistic chemistry; LES of non-premixed turbulent reacting flows with conditional source term estimation; measurements of the three-dimensional scalar dissipation rate in gas-phase planar turbulent jets; direct simulation of a jet diffusion flame; on the use of interpolating wavelets in the direct numerical simulation of combustion; on the use of a dynamically adaptive wavelet collocation algorithm in DNS (direct numerical simulation) of non-premixed turbulent combustion; 2D simulations of Hall thrusters; computation of trailing-edge noise at low mach number using LES and acoustic analogy; weakly nonlinear modeling of the early stages of bypass transition; interactions between freestream turbulence and boundary layers; interfaces at the outer boundaries of turbulent motions; largest scales of turbulent wall flows; the instability of streaks in near-wall turbulence; an implementation of the v(sup 2) - f model with application to transonic flows; heat transfer predictions in cavities; a structure-based model with stropholysis effects; modeling a confined swirling coaxial jet; subgrid-scale models based on incremental unknowns for large eddy simulations; subgrid scale modeling taking the numerical error into consideration; towards a near-wall model for LES of a separated diffuser flow; on the feasibility of merging LES with RANS (Reynolds Averaging Numerical simulation) for the near-wall region of attached turbulent flows; large-eddy simulation of a separated boundary layer; numerical study of a channel flow with variable properties; on the construction of high order finite difference schemes on non-uniform meshes with good conservation properties; development of immersed boundary methods for complex geometries; and particle methods for micro and macroscale flow simulations.

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.

Dynamic behavior of turbulent flow in a widely-spaced co-axial jet diffusion flame combustor

NASA Astrophysics Data System (ADS)

Sturgess, G. J.; Syed, S. A.

1983-01-01

Reacting flows in a bluff-body stabilized diffusion flame research combustor operated by the Wright Aeronautical Propulsion Laboratory exhibit the presence of coherent structures where, because of dynamic behavior the flame consists of large, discrete flame eddies passing down the combustion tunnel separated in time by axial regions where no flame is visible. It is proposed that the formation of these structures and their subsequent behavior are the result of vortex-shedding from the flameholder and, in the main, interaction with the organ-pipe natural frequencies of the long combustion tunnel. A simulation of the flow is made based on a finite difference solution of the time-average, steady state, elliptic form of the Reynolds equations using the two-equation turbulence model and a 'mixed is burned' combustion model for closure. The simulation of the eddies and, in conjunction with a universal Strouhal number-Reynolds number correlation, provides successful prediction of the flame frequencies.

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.