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

DTIC Science & Technology

2015-08-01

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

Sediment Suspension by Straining-Induced Convection at the Head of Salinity Intrusion

NASA Astrophysics Data System (ADS)

Zhang, Qianjiang; Wu, Jiaxue

2018-01-01

The tidal straining can generate convective motions and exert a periodic modification of turbulence and sediment transport in estuarine and coastal bottom boundary layers. However, the evidence and physics of convection and sediment suspension induced by tidal straining have not been straightforward. To examine these questions, mooring and transect surveys have been conducted in September 2015 in the region of the Yangtze River plume influence. Field observations and scaling analyses indicate an occurrence of convective motions at the head of saline wedge. Theoretical analyses of stratification evolution in the saline wedge show that unstable stratification and resultant convection are induced by tidal straining. Vertical turbulent velocity and eddy viscosity at the head of saline wedge are both larger than their neutral counterparts in the main body, largely enhancing sediment suspension at the head of saline wedge. Moreover, sediment suspension in both neutral and convection-affected flows is supported by the variance of vertical turbulent velocity, rather than the shearing stress. Finally, the stability correction functions in the Monin-Obukhov similarity theory can be simply derived from the local turbulent kinetic energy balance to successfully describe the effects of tidal straining on turbulent length scale, eddy viscosity, and sediment diffusivity in the convection-affected flow. These recognitions may provide novel understanding of estuarine turbidity maxima, and the dynamical structure and processes for coastal hypoxia.

A priori testing of subgrid-scale models for the velocity-pressure and vorticity-velocity formulations

NASA Technical Reports Server (NTRS)

Winckelmans, G. S.; Lund, T. S.; Carati, D.; Wray, A. A.

1996-01-01

Subgrid-scale models for Large Eddy Simulation (LES) in both the velocity-pressure and the vorticity-velocity formulations were evaluated and compared in a priori tests using spectral Direct Numerical Simulation (DNS) databases of isotropic turbulence: 128(exp 3) DNS of forced turbulence (Re(sub(lambda))=95.8) filtered, using the sharp cutoff filter, to both 32(exp 3) and 16(exp 3) synthetic LES fields; 512(exp 3) DNS of decaying turbulence (Re(sub(Lambda))=63.5) filtered to both 64(exp 3) and 32(exp 3) LES fields. Gaussian and top-hat filters were also used with the 128(exp 3) database. Different LES models were evaluated for each formulation: eddy-viscosity models, hyper eddy-viscosity models, mixed models, and scale-similarity models. Correlations between exact versus modeled subgrid-scale quantities were measured at three levels: tensor (traceless), vector (solenoidal 'force'), and scalar (dissipation) levels, and for both cases of uniform and variable coefficient(s). Different choices for the 1/T scaling appearing in the eddy-viscosity were also evaluated. It was found that the models for the vorticity-velocity formulation produce higher correlations with the filtered DNS data than their counterpart in the velocity-pressure formulation. It was also found that the hyper eddy-viscosity model performs better than the eddy viscosity model, in both formulations.

A kinematic eddy viscosity model including the influence of density variations and preturbulence

NASA Technical Reports Server (NTRS)

Cohen, L. S.

1973-01-01

A model for the kinematic eddy viscosity was developed which accounts for the turbulence produced as a result of jet interactions between adjacent streams as well as the turbulence initially present in the streams. In order to describe the turbulence contribution from jet interaction, the eddy viscosity suggested by Prandtl was adopted, and a modification was introduced to account for the effect of density variation through the mixing layer. The form of the modification was ascertained from a study of the compressible turbulent boundary layer on a flat plate. A kinematic eddy viscosity relation which corresponds to the initial turbulence contribution was derived by employing arguments used by Prandtl in his mixing length hypothesis. The resulting expression for self-preserving flow is similar to that which describes the mixing of a submerged jet. Application of the model has led to analytical predictions which are in good agreement with available turbulent mixing experimental data.

SGS Dynamics and Modeling near a Rough Wall.

NASA Astrophysics Data System (ADS)

Juneja, Anurag; Brasseur, James G.

1998-11-01

Large-eddy simulation (LES) of the atmospheric boundary layer (ABL) using classical subgrid-scale (SGS) models is known to poorly predict mean shear at the first few grid cells near the rough surface, creating error which can propogate vertically to infect the entire ABL. Our goal was to determine the first-order errors in predicted SGS terms that arise as a consequence of necessary under-resolution of integral scales and anisotropy which exist at the first few grid levels in LES of rough wall turbulence. Analyzing the terms predicted from eddy-viscosity and similarity closures with DNS anisotropic datasets of buoyancy- and shear-driven turbulence, we uncover three important issues which should be addressed in the design of SGS closures for rough walls and we provide a priori tests for the SGS model. Firstly, we identify a strong spurious coupling between the anisotropic structure of the resolved velocity field and predicted SGS dynamics which can create a feedback loop to incorrectly enhance certain components of the predicted resolved velocity. Secondly, we find that eddy viscosity and similarity SGS models do not contain enough degrees of freedom to capture, at a sufficient level of accuracy, both RS-SGS energy flux and SGS-RS dynamics. Thirdly, to correctly capture pressure transport near a wall, closures must be made more flexible to accommodate proper partitioning between SGS stress divergence and SGS pressure gradient.

Renormalization-group theory for the eddy viscosity in subgrid modeling

NASA Technical Reports Server (NTRS)

Zhou, YE; Vahala, George; Hossain, Murshed

1988-01-01

Renormalization-group theory is applied to incompressible three-dimensional Navier-Stokes turbulence so as to eliminate unresolvable small scales. The renormalized Navier-Stokes equation now includes a triple nonlinearity with the eddy viscosity exhibiting a mild cusp behavior, in qualitative agreement with the test-field model results of Kraichnan. For the cusp behavior to arise, not only is the triple nonlinearity necessary but the effects of pressure must be incorporated in the triple term. The renormalized eddy viscosity will not exhibit a cusp behavior if it is assumed that a spectral gap exists between the large and small scales.

Solar tachocline dynamics: eddy viscosity, anti-friction, or something in between?

NASA Astrophysics Data System (ADS)

McIntyre, Michael E.

The tachocline has values of the stratification or buoyancy frequency N two or more orders of magnitude greater than the Coriolis frequency. In this and other respects it is very like the Earth's atmosphere, viewed globally, except that the Earth's solid surface is replaced by an abrupt, magnetically-constrained "tachopause" (Gough & McIntyre 1998). The tachocline is helium-poor through fast ventilation from above, down to the tachopause, on timescales of only a few million years. The corresponding sound-speed anomaly fits helioseismic data with a tachocline thickness (0.019±0.001) Rsolar, about 0.13×105km (Elliott & Gough 1999), implying large values of the gradient Richardson number such that stratification dominates vertical shear even more strongly than in the Earth's stratosphere, as earlier postulated by Spiegel & Zahn (1992). Therefore the tachocline ventilation circulation cannot be driven by vertically-transmitted frictional torques, any more than the ozone-transporting circulation and differential rotation of the Earth's stratosphere can thus be driven. Rather, the tachocline circulation must be driven mainly by the Reynolds and Maxwell stresses interior to the convection zone, through a gyroscopic pumping action and the downward-burrowing response to it. If layerwise-two-dimensional turbulence is important, then because of its potential-vorticity-transporting properties the effect will be anti-frictional rather than eddy-viscosity-like. In order to correctly predict the differential rotation of the Sun's convection zone, even qualitatively, a convection-zone model must be fully coupled to a tachocline model.

Large eddy simulation of spanwise rotating turbulent channel flow with dynamic variants of eddy viscosity model

NASA Astrophysics Data System (ADS)

Jiang, Zhou; Xia, Zhenhua; Shi, Yipeng; Chen, Shiyi

2018-04-01

A fully developed spanwise rotating turbulent channel flow has been numerically investigated utilizing large-eddy simulation. Our focus is to assess the performances of the dynamic variants of eddy viscosity models, including dynamic Vreman's model (DVM), dynamic wall adapting local eddy viscosity (DWALE) model, dynamic σ (Dσ ) model, and the dynamic volumetric strain-stretching (DVSS) model, in this canonical flow. The results with dynamic Smagorinsky model (DSM) and direct numerical simulations (DNS) are used as references. Our results show that the DVM has a wrong asymptotic behavior in the near wall region, while the other three models can correctly predict it. In the high rotation case, the DWALE can get reliable mean velocity profile, but the turbulence intensities in the wall-normal and spanwise directions show clear deviations from DNS data. DVSS exhibits poor predictions on both the mean velocity profile and turbulence intensities. In all three cases, Dσ performs the best.

Transport Coefficients in weakly compressible turbulence

NASA Technical Reports Server (NTRS)

Rubinstein, Robert; Erlebacher, Gordon

1996-01-01

A theory of transport coefficients in weakly compressible turbulence is derived by applying Yoshizawa's two-scale direct interaction approximation to the compressible equations of motion linearized about a state of incompressible turbulence. The result is a generalization of the eddy viscosity representation of incompressible turbulence. In addition to the usual incompressible eddy viscosity, the calculation generates eddy diffusivities for entropy and pressure, and an effective bulk viscosity acting on the mean flow. The compressible fluctuations also generate an effective turbulent mean pressure and corrections to the speed of sound. Finally, a prediction unique to Yoshizawa's two-scale approximation is that terms containing gradients of incompressible turbulence quantities also appear in the mean flow equations. The form these terms take is described.

On the Asymptotic Regimes and the Strongly Stratified Limit of Rotating Boussinesq Equations

NASA Technical Reports Server (NTRS)

Babin, A.; Mahalov, A.; Nicolaenko, B.; Zhou, Y.

1997-01-01

Asymptotic regimes of geophysical dynamics are described for different Burger number limits. Rotating Boussinesq equations are analyzed in the asymptotic limit, of strong stratification in the Burger number of order one situation as well as in the asymptotic regime of strong stratification and weak rotation. It is shown that in both regimes horizontally averaged buoyancy variable is an adiabatic invariant for the full Boussinesq system. Spectral phase shift corrections to the buoyancy time scale associated with vertical shearing of this invariant are deduced. Statistical dephasing effects induced by turbulent processes on inertial-gravity waves are evidenced. The 'split' of the energy transfer of the vortical and the wave components is established in the Craya-Herring cyclic basis. As the Burger number increases from zero to infinity, we demonstrate gradual unfreezing of energy cascades for ageostrophic dynamics. The energy spectrum and the anisotropic spectral eddy viscosity are deduced with an explicit dependence on the anisotropic rotation/stratification time scale which depends on the vertical aspect ratio parameter. Intermediate asymptotic regime corresponding to strong stratification and weak rotation is analyzed where the effects of weak rotation are accounted for by an asymptotic expansion with full control (saturation) of vertical shearing. The regularizing effect of weak rotation differs from regularizations based on vertical viscosity. Two scalar prognostic equations for ageostrophic components (divergent velocity potential and geostrophic departure ) are obtained.

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.

On the Subgrid-Scale Modeling of Compressible Turbulence

NASA Technical Reports Server (NTRS)

Squires, Kyle; Zeman, Otto

1990-01-01

A new sub-grid scale model is presented for the large-eddy simulation of compressible turbulence. In the proposed model, compressibility contributions have been incorporated in the sub-grid scale eddy viscosity which, in the incompressible limit, reduce to a form originally proposed by Smagorinsky (1963). The model has been tested against a simple extension of the traditional Smagorinsky eddy viscosity model using simulations of decaying, compressible homogeneous turbulence. Simulation results show that the proposed model provides greater dissipation of the compressive modes of the resolved-scale velocity field than does the Smagorinsky eddy viscosity model. For an initial r.m.s. turbulence Mach number of 1.0, simulations performed using the Smagorinsky model become physically unrealizable (i.e., negative energies) because of the inability of the model to sufficiently dissipate fluctuations due to resolved scale velocity dilations. The proposed model is able to provide the necessary dissipation of this energy and maintain the realizability of the flow. Following Zeman (1990), turbulent shocklets are considered to dissipate energy independent of the Kolmogorov energy cascade. A possible parameterization of dissipation by turbulent shocklets for Large-Eddy Simulation is also presented.

Effect of turbulent eddy viscosity on the unstable surface mode above an acoustic liner

NASA Astrophysics Data System (ADS)

Marx, David; Aurégan, Yves

2013-07-01

Lined ducts are used to reduce noise radiation from ducts in turbofan engines. In certain conditions they may sustain hydrodynamic instabilities. A local linear stability analysis of the flow in a 2D lined channel is performed using a numerical integration of the governing equations. Several model equations are used, one of them taking into account turbulent eddy viscosity, and a realistic turbulent mean flow profile is used that vanishes at the wall. The stability analysis results are compared to published experimental results. Both the model and the experiments show the existence of an unstable mode, and the importance of taking into account eddy viscosity in the model is shown. When this is done, quantities such as the growth rate and the velocity eigenfunctions are shown to agree correctly.

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.

Mesospheric heating due to intense tropospheric convection

NASA Technical Reports Server (NTRS)

Taylor, L. L.

1979-01-01

A series of rocket measurements made twice daily at Wallops Island, Va., revealed a rapid heating of the mesosphere on the order of 10 K on days when thunderstorms or squall lines were in the area. This heating is explained as the result of frictional dissipation of vertically propagating internal gravity waves generated by intense tropospheric convection. Ray-tracing theory is used to determine the spectrum of gravity wave groups that actually reach mesospheric heights. This knowledge is used in an equation describing the spectral energy density of a penetrative convective element to calculate the fraction of the total energy initially available to excite those waves that do reach the level of heating. This value, converted into a vertical velocity, is used as the lower boundary condition for a multilayer model used to determine the detailed structure of the vertically propagating waves. The amount of frictional dissipation produced by the waves is calculated from the solutions of the frictionless model by use of a vertically varying eddy viscosity coefficient. The heating produced by the dissipation is then calculated from the thermodynamic equation.

Effects of eddy viscosity and thermal conduction and Coriolis force in the dynamics of gravity wave driven fluctuations in the OH nightglow

NASA Technical Reports Server (NTRS)

Hickey, M. P.

1988-01-01

The chemical-dynamical model of Walterscheid et al. (1987), which describes wave-driven fluctuations in OH nightglow, was modified to include the effects of both eddy thermal conduction and viscosity, as well as the Coriolis force (with the shallow atmosphere approximation). Using the new model, calculations were performed for the same nominal case as used by Walterscheid et al. but with only wave periods considered. For this case, the Coriolis force was found to be unimportant at any wave period. For wave periods greater than 2 or 3 hours, the inclusion of thermal conduction alone greatly modified the results (in terms of a complex ratio 'eta' which expresses the relationship between the intensity oscillation about the time-averaged intensity and the temperature oscillation about the time-averaged temperature); this effect was reduced with the further inclusion of the eddy viscosity.

Numerical dissipation vs. subgrid-scale modelling for large eddy simulation

NASA Astrophysics Data System (ADS)

Dairay, Thibault; Lamballais, Eric; Laizet, Sylvain; Vassilicos, John Christos

2017-05-01

This study presents an alternative way to perform large eddy simulation based on a targeted numerical dissipation introduced by the discretization of the viscous term. It is shown that this regularisation technique is equivalent to the use of spectral vanishing viscosity. The flexibility of the method ensures high-order accuracy while controlling the level and spectral features of this purely numerical viscosity. A Pao-like spectral closure based on physical arguments is used to scale this numerical viscosity a priori. It is shown that this way of approaching large eddy simulation is more efficient and accurate than the use of the very popular Smagorinsky model in standard as well as in dynamic version. The main strength of being able to correctly calibrate numerical dissipation is the possibility to regularise the solution at the mesh scale. Thanks to this property, it is shown that the solution can be seen as numerically converged. Conversely, the two versions of the Smagorinsky model are found unable to ensure regularisation while showing a strong sensitivity to numerical errors. The originality of the present approach is that it can be viewed as implicit large eddy simulation, in the sense that the numerical error is the source of artificial dissipation, but also as explicit subgrid-scale modelling, because of the equivalence with spectral viscosity prescribed on a physical basis.

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.

Local dynamic subgrid-scale models in channel flow

NASA Technical Reports Server (NTRS)

Cabot, William H.

1994-01-01

The dynamic subgrid-scale (SGS) model has given good results in the large-eddy simulation (LES) of homogeneous isotropic or shear flow, and in the LES of channel flow, using averaging in two or three homogeneous directions (the DA model). In order to simulate flows in general, complex geometries (with few or no homogeneous directions), the dynamic SGS model needs to be applied at a local level in a numerically stable way. Channel flow, which is inhomogeneous and wall-bounded flow in only one direction, provides a good initial test for local SGS models. Tests of the dynamic localization model were performed previously in channel flow using a pseudospectral code and good results were obtained. Numerical instability due to persistently negative eddy viscosity was avoided by either constraining the eddy viscosity to be positive or by limiting the time that eddy viscosities could remain negative by co-evolving the SGS kinetic energy (the DLk model). The DLk model, however, was too expensive to run in the pseudospectral code due to a large near-wall term in the auxiliary SGS kinetic energy (k) equation. One objective was then to implement the DLk model in a second-order central finite difference channel code, in which the auxiliary k equation could be integrated implicitly in time at great reduction in cost, and to assess its performance in comparison with the plane-averaged dynamic model or with no model at all, and with direct numerical simulation (DNS) and/or experimental data. Other local dynamic SGS models have been proposed recently, e.g., constrained dynamic models with random backscatter, and with eddy viscosity terms that are averaged in time over material path lines rather than in space. Another objective was to incorporate and test these models in channel flow.

Prediction of Transonic Vortex Flows Using Linear and Nonlinear Turbulent Eddy Viscosity Models

NASA Technical Reports Server (NTRS)

Bartels, Robert E.; Gatski, Thomas B.

2000-01-01

Three-dimensional transonic flow over a delta wing is investigated with a focus on the effect of transition and influence of turbulence stress anisotropies. The performance of linear eddy viscosity models and an explicit algebraic stress model is assessed at the start of vortex flow, and the results compared with experimental data. To assess the effect of transition location, computations that either fix transition or are fully turbulent are performed. To assess the effect of the turbulent stress anisotropy, comparisons are made between predictions from the algebraic stress model and the linear eddy viscosity models. Both transition location and turbulent stress anisotropy significantly affect the 3D flow field. The most significant effect is found to be the modeling of transition location. At a Mach number of 0.90, the computed solution changes character from steady to unsteady depending on transition onset. Accounting for the anisotropies in the turbulent stresses also considerably impacts the flow, most notably in the outboard region of flow separation.

Convection without eddy viscosity: An attempt to model the interiors of giant planets

NASA Technical Reports Server (NTRS)

Ingersoll, A. P.

1986-01-01

In the theory of hydrostatic quasi-geostrophic flow in the Earth's atmosphere the principal results do not depend on the eddy viscosity. This contrasts with published theories of convection in deep rotating fluid spheres, where the wavelength of the fastest growing disturbance varies as E sup 1/3, where E, the Ekman number, is proportional to the eddy viscosity. A new theory of quasi-columnar motions in stably stratified fluid spheres attempts to capture the luck of the meteorologists. The theory allows one to investigate the stability of barotropic and baroclinic zonal flows that extend into the planetary interior. It is hypothesized that the internal heat Jupiter and Saturn comes out not radially but on sloping surfaces defined by the internal entropy distribution. To test the hypothesis one searches for basic states in which the wavelength of the fastest-growing disturbance remains finite as E tends to zero, and is which the heat flux vector is radially outward and poleward.

Hybrid Reynolds-Averaged/Large Eddy Simulation of a Cavity Flameholder; Assessment of Modeling Sensitivities

NASA Technical Reports Server (NTRS)

Baurle, R. A.

2015-01-01

Steady-state and scale-resolving simulations have been performed for flow in and around a model scramjet combustor flameholder. The cases simulated corresponded to those used to examine this flowfield experimentally using particle image velocimetry. A variety of turbulence models were used for the steady-state Reynolds-averaged simulations which included both linear and non-linear eddy viscosity models. The scale-resolving simulations used a hybrid Reynolds-averaged / large eddy simulation strategy that is designed to be a large eddy simulation everywhere except in the inner portion (log layer and below) of the boundary layer. Hence, this formulation can be regarded as a wall-modeled large eddy simulation. This effort was undertaken to formally assess the performance of the hybrid Reynolds-averaged / large eddy simulation modeling approach in a flowfield of interest to the scramjet research community. The numerical errors were quantified for both the steady-state and scale-resolving simulations prior to making any claims of predictive accuracy relative to the measurements. The steady-state Reynolds-averaged results showed a high degree of variability when comparing the predictions obtained from each turbulence model, with the non-linear eddy viscosity model (an explicit algebraic stress model) providing the most accurate prediction of the measured values. The hybrid Reynolds-averaged/large eddy simulation results were carefully scrutinized to ensure that even the coarsest grid had an acceptable level of resolution for large eddy simulation, and that the time-averaged statistics were acceptably accurate. The autocorrelation and its Fourier transform were the primary tools used for this assessment. The statistics extracted from the hybrid simulation strategy proved to be more accurate than the Reynolds-averaged results obtained using the linear eddy viscosity models. However, there was no predictive improvement noted over the results obtained from the explicit Reynolds stress model. Fortunately, the numerical error assessment at most of the axial stations used to compare with measurements clearly indicated that the scale-resolving simulations were improving (i.e. approaching the measured values) as the grid was refined. Hence, unlike a Reynolds-averaged simulation, the hybrid approach provides a mechanism to the end-user for reducing model-form errors.

Subregional characterization of mesoscale eddies across the Brazil-Malvinas Confluence

NASA Astrophysics Data System (ADS)

Mason, Evan; Pascual, Ananda; Gaube, Peter; Ruiz, Simón; Pelegrí, Josep L.; Delepoulle, Antoine

2017-04-01

Horizontal and vertical motions associated with coherent mesoscale structures, including eddies and meanders, are responsible for significant global transports of many properties, including heat and mass. Mesoscale vertical fluxes also influence upper ocean biological productivity by mediating the supply of nutrients into the euphotic layer, with potential impacts on the global carbon cycle. The Brazil-Malvinas Confluence (BMC) is a western boundary current region in the South Atlantic with intense mesoscale activity. This region has an active role in the genesis and transformation of water masses and thus is a critical component of the Atlantic meridional overturning circulation. The collision between the Malvinas and Brazil Currents over the Patagonian shelf/slope creates an energetic front that translates offshore to form a vigorous eddy field. Recent improvements in gridded altimetric sea level anomaly fields allow us to track BMC mesoscale eddies with high spatial and temporal resolutions using an automated eddy tracker. We characterize the eddies across fourteen 5° × 5° subregions. Eddy-centric composites of tracers and geostrophic currents diagnosed from a global reanalysis of surface and in situ data reveal substantial subregional heterogeneity. The in situ data are also used to compute the evolving quasi-geostrophic vertical velocity (QG-ω) associated with each instantaneous eddy instance. The QG-ω eddy composites have the expected dipole patterns of alternating upwelling/downwelling, however, the magnitude and sign of azimuthally averaged vertical velocity varies among subregions. Maximum eddy values are found near fronts and sharp topographic gradients. In comparison with regional eddy composites, subregional composites provide refined information about mesoscale eddy heterogeneity.

Water hammer prediction and control: the Green's function method

NASA Astrophysics Data System (ADS)

Xuan, Li-Jun; Mao, Feng; Wu, Jie-Zhi

2012-04-01

By Green's function method we show that the water hammer (WH) can be analytically predicted for both laminar and turbulent flows (for the latter, with an eddy viscosity depending solely on the space coordinates), and thus its hazardous effect can be rationally controlled and minimized. To this end, we generalize a laminar water hammer equation of Wang et al. (J. Hydrodynamics, B2, 51, 1995) to include arbitrary initial condition and variable viscosity, and obtain its solution by Green's function method. The predicted characteristic WH behaviors by the solutions are in excellent agreement with both direct numerical simulation of the original governing equations and, by adjusting the eddy viscosity coefficient, experimentally measured turbulent flow data. Optimal WH control principle is thereby constructed and demonstrated.

The formation processes of phytoplankton growth and decline in mesoscale eddies in the western North Pacific Ocean

NASA Astrophysics Data System (ADS)

Chang, Yu-Lin; Miyazawa, Yasumasa; Oey, Lie-Yauw; Kodaira, Tsubasa; Huang, Shihming

2017-05-01

In this study, we investigate the processes of phytoplankton growth and decline in mesoscale eddies in the western North Pacific Ocean based on the in situ chlorophyll data obtained from 52 cruises conducted by the Japan Meteorological Agency together with idealized numerical simulations. Both the observation and model results suggest that chlorophyll/phytoplankton concentrations are higher in cold than in warm eddies in near-surface water (z > -70 m). In the idealized simulation, the isopycnal movements associated with upwelling/downwelling transport phytoplankton and nutrients to different vertical depths during eddy formation (stage A). Phytoplankton and nutrients in cold eddies is transported toward shallower waters while those in warm eddies move toward deeper waters. In the period after the eddy has formed (stage B), sunlight and initially upwelled nutrients together promote the growth of phytoplankton in cold eddies. Phytoplankton in warm eddies decays due to insufficient sunlight in deeper waters. In stage B, upwelling and downwelling coexist in both warm and cold eddies, contributing nearly equally to vertical displacement. The upwelling/downwelling-induced nitrate flux accounts for a small percentage (˜3%) of the total nitrate flux in stage B. The vertical velocity caused by propagating eddies, therefore, is not the primary factor causing differences in phytoplankton concentrations between stage-B warm and cold eddies.

Moffatt eddies at an interface

NASA Astrophysics Data System (ADS)

Shtern, Vladimir

2014-12-01

It is shown that an infinite set of eddies can develop near the interface-wall intersection in a two-fluid flow. A striking feature is that the eddy occurrence depends on from what side of the interface the flow is driven. In air-water flows where the viscosity ratio is 0.018, the eddies develop if a driving source is located on (i) the air side for , (ii) any side for , and (iii) the water side for , where is the upper interface-wall angle.

Evaluation of the Momentum Closure Schemes in MPAS-Ocean

NASA Astrophysics Data System (ADS)

Zhao, Shimei; Liu, Yudi; Liu, Wei

2018-04-01

In order to compare and evaluate the performances of the Laplacian viscosity closure, the biharmonic viscosity closure, and the Leith closure momentum schemes in the MPAS-Ocean model, a variety of physical quantities, such as the relative reference potential energy (RPE) change, the RPE time change rate (RPETCR), the grid Reynolds number, the root mean square (RMS) of kinetic energy, and the spectra of kinetic energy and enstrophy, are calculated on the basis of results of a 3D baroclinic periodic channel. Results indicate that: 1) The RPETCR demonstrates a saturation phenomenon in baroclinic eddy tests. The critical grid Reynolds number corresponding to RPETCR saturation differs between the three closures: the largest value is in the biharmonic viscosity closure, followed by that in the Laplacian viscosity closure, and that in the Leith closure is the smallest. 2) All three closures can effectively suppress spurious dianeutral mixing by reducing the grid Reynolds number under sub-saturation conditions of the RPETCR, but they can also damage certain physical processes. Generally, the damage to the rotation process is greater than that to the advection process. 3) The dissipation in the biharmonic viscosity closure is strongly dependent on scales. Most dissipation concentrates on small scales, and the energy of small-scale eddies is often transferred to large-scale kinetic energy. The viscous dissipation in the Laplacian viscosity closure is the strongest on various scales, followed by that in the Leith closure. Note that part of the small-scale kinetic energy is also transferred to large-scale kinetic energy in the Leith closure. 4) The characteristic length scale L and the dimensionless parameter D in the Leith closure are inherently coupled. The RPETCR is inversely proportional to the product of D and L. When the product of D and L is constant, both the simulated RPETCR and the inhibition of spurious dianeutral mixing are the same in all tests using the Leith closure. The dissipative scale in the Leith closure depends on the parameter L, and the dissipative intensity depends on the parameter D. 5) Although optimal results may not be achieved by using the optimal parameters obtained from the 2D barotropic model in the 3D baroclinic simulation, the total energies are dissipative in all three closures. Dissipation is the strongest in the biharmonic viscosity closure, followed by that in the Leith closure, and that in the Laplacian viscosity closure is the weakest. Mesoscale eddies develop the fastest in the biharmonic viscosity closure after the baroclinic adjustment process finishes, and the kinetic energy reaches its maximum, which is attributed to the smallest dissipation of enstrophy in the biharmonic viscosity closure. Mesoscale eddies develop the slowest, and the kinetic energy peak value is the smallest in the Laplacian viscosity closure. Results in the Leith closure are between that in the biharmonic viscosity closure and the Laplacian viscosity closure.

Non-closure of the surface energy balance explained by phase difference between vertical velocity and scalars of large atmospheric eddies

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

Gao, Zhongming; Liu, Heping; Katul, Gabriel G.

It is now accepted that large-scale turbulent eddies impact the widely reported non-closure of the surface energy balance when latent and sensible heat fluxes are measured using the eddy covariance method in the atmospheric surface layer (ASL). However, a mechanistic link between large eddies and non-closure of the surface energy balance remains a subject of inquiry. Here, measured 10 Hz time series of vertical velocity, air temperature, and water vapor density collected in the ASL are analyzed for conditions where entrainment and/or horizontal advection separately predominate. The series are decomposed into small- and large- eddies based on a frequency cutoffmore » and their contributions to turbulent fluxes are analyzed. Phase difference between vertical velocity and water vapor density associated with large eddies reduces latent heat fluxes, especially in conditions where advection prevails. Furthermore, enlarged phase difference of large eddies linked to entrainment or advection occurrence leads to increased residuals of the surface energy balance.« less

Non-closure of the surface energy balance explained by phase difference between vertical velocity and scalars of large atmospheric eddies

DOE PAGES

Gao, Zhongming; Liu, Heping; Katul, Gabriel G.; ...

2017-03-16

It is now accepted that large-scale turbulent eddies impact the widely reported non-closure of the surface energy balance when latent and sensible heat fluxes are measured using the eddy covariance method in the atmospheric surface layer (ASL). However, a mechanistic link between large eddies and non-closure of the surface energy balance remains a subject of inquiry. Here, measured 10 Hz time series of vertical velocity, air temperature, and water vapor density collected in the ASL are analyzed for conditions where entrainment and/or horizontal advection separately predominate. The series are decomposed into small- and large- eddies based on a frequency cutoffmore » and their contributions to turbulent fluxes are analyzed. Phase difference between vertical velocity and water vapor density associated with large eddies reduces latent heat fluxes, especially in conditions where advection prevails. Furthermore, enlarged phase difference of large eddies linked to entrainment or advection occurrence leads to increased residuals of the surface energy balance.« less

Development of a two-dimensional zonally averaged statistical-dynamical model. III - The parameterization of the eddy fluxes of heat and moisture

NASA Technical Reports Server (NTRS)

Stone, Peter H.; Yao, Mao-Sung

1990-01-01

A number of perpetual January simulations are carried out with a two-dimensional zonally averaged model employing various parameterizations of the eddy fluxes of heat (potential temperature) and moisture. The parameterizations are evaluated by comparing these results with the eddy fluxes calculated in a parallel simulation using a three-dimensional general circulation model with zonally symmetric forcing. The three-dimensional model's performance in turn is evaluated by comparing its results using realistic (nonsymmetric) boundary conditions with observations. Branscome's parameterization of the meridional eddy flux of heat and Leovy's parameterization of the meridional eddy flux of moisture simulate the seasonal and latitudinal variations of these fluxes reasonably well, while somewhat underestimating their magnitudes. New parameterizations of the vertical eddy fluxes are developed that take into account the enhancement of the eddy mixing slope in a growing baroclinic wave due to condensation, and also the effect of eddy fluctuations in relative humidity. The new parameterizations, when tested in the two-dimensional model, simulate the seasonal, latitudinal, and vertical variations of the vertical eddy fluxes quite well, when compared with the three-dimensional model, and only underestimate the magnitude of the fluxes by 10 to 20 percent.

Existence of k⁻¹ power-law scaling in the equilibrium regions of wall-bounded turbulence explained by Heisenberg's eddy viscosity.

PubMed

Katul, Gabriel G; Porporato, Amilcare; Nikora, Vladimir

2012-12-01

The existence of a "-1" power-law scaling at low wavenumbers in the longitudinal velocity spectrum of wall-bounded turbulence was explained by multiple mechanisms; however, experimental support has not been uniform across laboratory studies. This letter shows that Heisenberg's eddy viscosity approach can provide a theoretical framework that bridges these multiple mechanisms and explains the elusiveness of the "-1" power law in some experiments. Novel theoretical outcomes are conjectured about the role of intermittency and very-large scale motions in modifying the k⁻¹ scaling.

A Nonlinear Interactions Approximation Model for Large-Eddy Simulation

NASA Astrophysics Data System (ADS)

Haliloglu, Mehmet U.; Akhavan, Rayhaneh

2003-11-01

A new approach to LES modelling is proposed based on direct approximation of the nonlinear terms \\overlineu_iuj in the filtered Navier-Stokes equations, instead of the subgrid-scale stress, τ_ij. The proposed model, which we call the Nonlinear Interactions Approximation (NIA) model, uses graded filters and deconvolution to parameterize the local interactions across the LES cutoff, and a Smagorinsky eddy viscosity term to parameterize the distant interactions. A dynamic procedure is used to determine the unknown eddy viscosity coefficient, rendering the model free of adjustable parameters. The proposed NIA model has been applied to LES of turbulent channel flows at Re_τ ≈ 210 and Re_τ ≈ 570. The results show good agreement with DNS not only for the mean and resolved second-order turbulence statistics but also for the full (resolved plus subgrid) Reynolds stress and turbulence intensities.

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.

Eddy damping, backscatter, and subgrid stresses in subgrid modeling of turbulence

NASA Technical Reports Server (NTRS)

Zhou, YE

1991-01-01

The Brief Report demonstrates the relationship of eddy-viscosity models to subgrid stresses. A formula that determines the relative importance of the cross stress and the Reynolds stress for the net eddy-damping and backscatter contributions is derived. The cross-stress term with sharp-cut filtering is identified as an important quantity to model. These concepts could prove useful as a basis for constructing specific models for the Reynolds and cross stresses.

Enhanced Nitrogen Loss by Eddy-Induced Vertical Transport in the Offshore Peruvian Oxygen Minimum Zone

PubMed Central

Callbeck, Cameron M.; Lavik, Gaute; Stramma, Lothar; Kuypers, Marcel M. M.; Bristow, Laura A.

2017-01-01

The eastern tropical South Pacific (ETSP) upwelling region is one of the ocean’s largest sinks of fixed nitrogen, which is lost as N2 via the anaerobic processes of anammox and denitrification. One-third of nitrogen loss occurs in productive shelf waters stimulated by organic matter export as a result of eastern boundary upwelling. Offshore, nitrogen loss rates are lower, but due to its sheer size this area accounts for ~70% of ETSP nitrogen loss. How nitrogen loss and primary production are regulated in the offshore ETSP region where coastal upwelling is less influential remains unclear. Mesoscale eddies, ubiquitous in the ETSP region, have been suggested to enhance vertical nutrient transport and thereby regulate primary productivity and hence organic matter export. Here, we investigated the impact of mesoscale eddies on anammox and denitrification activity using 15N-labelled in situ incubation experiments. Anammox was shown to be the dominant nitrogen loss process, but varied across the eddy, whereas denitrification was below detection at all stations. Anammox rates at the eddy periphery were greater than at the center. Similarly, depth-integrated chlorophyll paralleled anammox activity, increasing at the periphery relative to the eddy center; suggestive of enhanced organic matter export along the periphery supporting nitrogen loss. This can be attributed to enhanced vertical nutrient transport caused by an eddy-driven submesoscale mechanism operating at the eddy periphery. In the ETSP region, the widespread distribution of eddies and the large heterogeneity observed in anammox rates from a compilation of stations suggests that eddy-driven vertical nutrient transport may regulate offshore primary production and thereby nitrogen loss. PMID:28122044

Some effects of swirl on turbulent mixing and combustion

NASA Technical Reports Server (NTRS)

Rubel, A.

1972-01-01

A general formulation of some effects of swirl on turbulent mixing is given. The basis for the analysis is that momentum transport is enhanced by turbulence resulting from rotational instability of the fluid field. An appropriate form for the turbulent eddy viscosity is obtained by mixing length type arguments. The result takes the form of a corrective factor that is a function of the swirl and acts to increase the eddy viscosity. The factor is based upon the initial mixing conditions implying that the rotational turbulence decays in a manner similar to that of free shear turbulence. Existing experimental data for free jet combustion are adequately matched by using the modifying factor to relate the effects of swirl on eddy viscosity. The model is extended and applied to the supersonic combustion of a ring jet of hydrogen injected into a constant area annular air stream. The computations demonstrate that swirling the flow could: (1) reduce the burning length by one half, (2) result in more uniform burning across the annulus width, and (3) open the possibility of optimization of the combustion characteristics by locating the fuel jet between the inner wall and center of the annulus width.

Induction logging device

DOEpatents

Koelle, A.R.; Landt, J.A.

An instrument is disclosed for mapping vertical conductive fractures in a resistive bedrock, magnetically inducing eddy currents by a pair of vertically oriented, mutually perpendicular, coplanar coils. The eddy currents drive magnetic fields which are picked up by a second, similar pair of coils.

A Nonlinear Dynamic Subscale Model for Partially Resolved Numerical Simulation (PRNS)/Very Large Eddy Simulation (VLES) of Internal Non-Reacting Flows

NASA Technical Reports Server (NTRS)

Shih, Tsan-Hsing; Liu, nan-Suey

2010-01-01

A brief introduction of the temporal filter based partially resolved numerical simulation/very large eddy simulation approach (PRNS/VLES) and its distinct features are presented. A nonlinear dynamic subscale model and its advantages over the linear subscale eddy viscosity model are described. In addition, a guideline for conducting a PRNS/VLES simulation is provided. Results are presented for three turbulent internal flows. The first one is the turbulent pipe flow at low and high Reynolds numbers to illustrate the basic features of PRNS/VLES; the second one is the swirling turbulent flow in a LM6000 single injector to further demonstrate the differences in the calculated flow fields resulting from the nonlinear model versus the pure eddy viscosity model; the third one is a more complex turbulent flow generated in a single-element lean direct injection (LDI) combustor, the calculated result has demonstrated that the current PRNS/VLES approach is capable of capturing the dynamically important, unsteady turbulent structures while using a relatively coarse grid.

Application of renormalization group theory to the large-eddy simulation of transitional boundary layers

NASA Technical Reports Server (NTRS)

Piomelli, Ugo; Zang, Thomas A.; Speziale, Charles G.; Lund, Thomas S.

1990-01-01

An eddy viscosity model based on the renormalization group theory of Yakhot and Orszag (1986) is applied to the large-eddy simulation of transition in a flat-plate boundary layer. The simulation predicts with satisfactory accuracy the mean velocity and Reynolds stress profiles, as well as the development of the important scales of motion. The evolution of the structures characteristic of the nonlinear stages of transition is also predicted reasonably well.

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.

Improving a two-equation eddy-viscosity turbulence model to predict the aerodynamic performance of thick wind turbine airfoils

NASA Astrophysics Data System (ADS)

Bangga, Galih; Kusumadewi, Tri; Hutomo, Go; Sabila, Ahmad; Syawitri, Taurista; Setiadi, Herlambang; Faisal, Muhamad; Wiranegara, Raditya; Hendranata, Yongki; Lastomo, Dwi; Putra, Louis; Kristiadi, Stefanus

2018-03-01

Numerical simulations for relatively thick airfoils are carried out in the present studies. An attempt to improve the accuracy of the numerical predictions is done by adjusting the turbulent viscosity of the eddy-viscosity Menter Shear-Stress-Transport (SST) model. The modification involves the addition of a damping factor on the wall-bounded flows incorporating the ratio of the turbulent kinetic energy to its specific dissipation rate for separation detection. The results are compared with available experimental data and CFD simulations using the original Menter SST model. The present model improves the lift polar prediction even though the stall angle is still overestimated. The improvement is caused by the better prediction of separated flow under a strong adverse pressure gradient. The results show that the Reynolds stresses are damped near the wall causing variation of the logarithmic velocity profiles.

Evaluation of machine learning algorithms for prediction of regions of high Reynolds averaged Navier Stokes uncertainty

DOE PAGES

Ling, Julia; Templeton, Jeremy Alan

2015-08-04

Reynolds Averaged Navier Stokes (RANS) models are widely used in industry to predict fluid flows, despite their acknowledged deficiencies. Not only do RANS models often produce inaccurate flow predictions, but there are very limited diagnostics available to assess RANS accuracy for a given flow configuration. If experimental or higher fidelity simulation results are not available for RANS validation, there is no reliable method to evaluate RANS accuracy. This paper explores the potential of utilizing machine learning algorithms to identify regions of high RANS uncertainty. Three different machine learning algorithms were evaluated: support vector machines, Adaboost decision trees, and random forests.more » The algorithms were trained on a database of canonical flow configurations for which validated direct numerical simulation or large eddy simulation results were available, and were used to classify RANS results on a point-by-point basis as having either high or low uncertainty, based on the breakdown of specific RANS modeling assumptions. Classifiers were developed for three different basic RANS eddy viscosity model assumptions: the isotropy of the eddy viscosity, the linearity of the Boussinesq hypothesis, and the non-negativity of the eddy viscosity. It is shown that these classifiers are able to generalize to flows substantially different from those on which they were trained. As a result, feature selection techniques, model evaluation, and extrapolation detection are discussed in the context of turbulence modeling applications.« less

LARGE EDDY SIMULATION OF TURBULENT FLOW OVER MARGINALLY RESOLVED THREE BLUFF BODIES USING AN IMMERSED BOUNDARY METHOD AND LAGRANGIAN DYNAMIC EDDY-VISCOSITY MODELS. (R828771C004)

EPA Science Inventory

The perspectives, information and conclusions conveyed in research project abstracts, progress reports, final reports, journal abstracts and journal publications convey the viewpoints of the principal investigator and may not represent the views and policies of ORD and EPA. Concl...

Eddy-induced Sea Surface Salinity changes in the tropical Pacific

NASA Astrophysics Data System (ADS)

Delcroix, T. C.; Chaigneau, A.; Soviadan, D.; Boutin, J.

2017-12-01

We analyse the Sea Surface Salinity (SSS) signature of westward propagating mesoscale eddies in the tropical Pacific by collocating 5 years (2010-2015) of SMOS (Soil Moisture and Ocean Salinity) SSS and altimetry-derived sea level anomalies. The main characteristics of mesoscale eddies are first identified in SLA maps. Composite analyses in the Central and Eastern ITCZ regions then reveal regionally dependent impacts with opposite SSS anomalies for the cyclonic and anticyclonic eddies. In the Central region (where we have the largest meridional SSS gradient), we found dipole-like SSS changes with maximum anomalies on the leading edge of the eddy. In the Eastern region (where we have the largest near-surface vertical salinity gradient) we found monopole-like SSS changes with maximum anomalies in the eddy centre. These dipole/monopole patterns and the rotational sense of eddies suggest the dominant role of horizontal and vertical advection in the Central and Eastern ITCZ regions, respectively.

A study of overflow simulations using MPAS-Ocean: Vertical grids, resolution, and viscosity

NASA Astrophysics Data System (ADS)

Reckinger, Shanon M.; Petersen, Mark R.; Reckinger, Scott J.

2015-12-01

MPAS-Ocean is used to simulate an idealized, density-driven overflow using the dynamics of overflow mixing and entrainment (DOME) setup. Numerical simulations are carried out using three of the vertical coordinate types available in MPAS-Ocean, including z-star with partial bottom cells, z-star with full cells, and sigma coordinates. The results are first benchmarked against other models, including the MITgcm's z-coordinate model and HIM's isopycnal coordinate model, which are used to set the base case used for this work. A full parameter study is presented that looks at how sensitive overflow simulations are to vertical grid type, resolution, and viscosity. Horizontal resolutions with 50 km grid cells are under-resolved and produce poor results, regardless of other parameter settings. Vertical grids ranging in thickness from 15 m to 120 m were tested. A horizontal resolution of 10 km and a vertical resolution of 60 m are sufficient to resolve the mesoscale dynamics of the DOME configuration, which mimics real-world overflow parameters. Mixing and final buoyancy are least sensitive to horizontal viscosity, but strongly sensitive to vertical viscosity. This suggests that vertical viscosity could be adjusted in overflow water formation regions to influence mixing and product water characteristics. Lastly, the study shows that sigma coordinates produce much less mixing than z-type coordinates, resulting in heavier plumes that go further down slope. Sigma coordinates are less sensitive to changes in resolution but as sensitive to vertical viscosity compared to z-coordinates.

Electromagnetic Modeling of the Passive Stabilization Loop at EAST

NASA Astrophysics Data System (ADS)

Ji, Xiang; Song, Yuntao; Wu, Songtao; Wang, Zhibin; Shen, Guang; Liu, Xufeng; Cao, Lei; Zhou, Zibo; Peng, Xuebing; Wang, Chenghao

2012-09-01

A passive stabilization loop (PSL) has been designed and manufactured in order to enhance the control of vertical instability and accommodate the new stage for high-performance plasma at EAST. Eddy currents are induced by vertical displacement events (VDEs) and disruption, which can produce a magnetic field to control the vertical instability of the plasma in a short timescale. A finite element model is created and meshed using ANSYS software. Based on the simulation of plasma VDEs and disruption, the distribution and decay curve of the eddy currents on the PSL are obtained. The largest eddy current is 200 kA and the stress is 68 MPa at the outer current bridge, which is the weakest point of the PSL because of the eddy currents and the magnetic fields. The analysis results provide the supporting data for the structural design.

Induction logging device with a pair of mutually perpendicular bucking coils

DOEpatents

Koelle, Alfred R.; Landt, Jeremy A.

1981-01-01

An instrument is disclosed for mapping vertical conductive fractures in a resistive bedrock, magnetically inducing eddy currents by a pair of vertically oriented, mutually perpendicular, coplanar coils. The eddy currents drive magnetic fields which are picked up by a second, similar pair of coils.

Atmospheric flow over two-dimensional bluff surface obstructions

NASA Technical Reports Server (NTRS)

Bitte, J.; Frost, W.

1976-01-01

The phenomenon of atmospheric flow over a two-dimensional surface obstruction, such as a building (modeled as a rectangular block, a fence or a forward-facing step), is analyzed by three methods: (1) an inviscid free streamline approach, (2) a turbulent boundary layer approach using an eddy viscosity turbulence model and a horizontal pressure gradient determined by the inviscid model, and (3) an approach using the full Navier-Stokes equations with three turbulence models; i.e., an eddy viscosity model, a turbulence kinetic-energy model and a two-equation model with an additional transport equation for the turbulence length scale. A comparison of the performance of the different turbulence models is given, indicating that only the two-equation model adequately accounts for the convective character of turbulence. Turbulence flow property predictions obtained from the turbulence kinetic-energy model with prescribed length scale are only insignificantly better than those obtained from the eddy viscosity model. A parametric study includes the effects of the variation of the characteristics parameters of the assumed logarithmic approach velocity profile. For the case of the forward-facing step, it is shown that in the downstream flow region an increase of the surface roughness gives rise to higher turbulence levels in the shear layer originating from the step corner.

Modifications of the law of the wall and algebraic turbulence modelling for separated boundary layers

NASA Technical Reports Server (NTRS)

Baldwin, B. S.; Maccormack, R. W.

1976-01-01

Various modifications of the conventional algebraic eddy viscosity turbulence model are investigated for application to separated flows. Friction velocity is defined in a way that avoids singular behavior at separation and reattachment but reverts to the conventional definition for flows with small pressure gradients. This leads to a modified law of the wall for separated flows. The effect on the calculated flow field of changes in the model that affect the eddy viscosity at various distances from the wall are determined by (1) switching from Prandtl's form to an inner layer formula due to Clauser at various distances from the wall, (2) varying the constant in the Van Driest damping factor, (3) using Clauser's inner layer formula all the way to the wall, and (4) applying a relaxation procedure in the evaluation of the constant in Clauser's inner layer formula. Numerical solutions of the compressible Navier-Stokes equations are used to determine the effects of the modifications. Experimental results from shock-induced separated flows at Mach numbers 2.93 and 8.45 are used for comparison. For these cases improved predictions of wall pressure distribution and positions of separation and reattachment are obtained from the relaxation version of the Clauser inner layer eddy viscosity formula.

Crossover between two- and three-dimensional turbulence in spatial mixing layers

NASA Astrophysics Data System (ADS)

Biancofiore, Luca

2016-11-01

We investigate how the domain depth affects the turbulent behaviour in spatially developing mixing layers by means of large-eddy simulations (LES) based on a spectral vanishing viscosity technique. Analyses of spectra of the vertical velocity, of Lumley's diagrams, of the turbulent kinetic energy and of the vortex stretching show that a two-dimensional behaviour of the turbulence is promoted in spatial mixing layers by constricting the fluid motion in one direction. This finding is in agreement with previous works on turbulent systems constrained by a geometric anisotropy, pioneered by Smith, Chasnov & Waleffe. We observe that the growth of the momentum thickness along the streamwise direction is damped in a confined domain. A full two-dimensional turbulent behaviour is observed when the momentum thickness is of the same order of magnitude as the confining scale.

Learning an Eddy Viscosity Model Using Shrinkage and Bayesian Calibration: A Jet-in-Crossflow Case Study

DOE PAGES

Ray, Jaideep; Lefantzi, Sophia; Arunajatesan, Srinivasan; ...

2017-09-07

In this paper, we demonstrate a statistical procedure for learning a high-order eddy viscosity model (EVM) from experimental data and using it to improve the predictive skill of a Reynolds-averaged Navier–Stokes (RANS) simulator. The method is tested in a three-dimensional (3D), transonic jet-in-crossflow (JIC) configuration. The process starts with a cubic eddy viscosity model (CEVM) developed for incompressible flows. It is fitted to limited experimental JIC data using shrinkage regression. The shrinkage process removes all the terms from the model, except an intercept, a linear term, and a quadratic one involving the square of the vorticity. The shrunk eddy viscositymore » model is implemented in an RANS simulator and calibrated, using vorticity measurements, to infer three parameters. The calibration is Bayesian and is solved using a Markov chain Monte Carlo (MCMC) method. A 3D probability density distribution for the inferred parameters is constructed, thus quantifying the uncertainty in the estimate. The phenomenal cost of using a 3D flow simulator inside an MCMC loop is mitigated by using surrogate models (“curve-fits”). A support vector machine classifier (SVMC) is used to impose our prior belief regarding parameter values, specifically to exclude nonphysical parameter combinations. The calibrated model is compared, in terms of its predictive skill, to simulations using uncalibrated linear and CEVMs. Finally, we find that the calibrated model, with one quadratic term, is more accurate than the uncalibrated simulator. The model is also checked at a flow condition at which the model was not calibrated.« less

Learning an Eddy Viscosity Model Using Shrinkage and Bayesian Calibration: A Jet-in-Crossflow Case Study

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

Ray, Jaideep; Lefantzi, Sophia; Arunajatesan, Srinivasan

In this paper, we demonstrate a statistical procedure for learning a high-order eddy viscosity model (EVM) from experimental data and using it to improve the predictive skill of a Reynolds-averaged Navier–Stokes (RANS) simulator. The method is tested in a three-dimensional (3D), transonic jet-in-crossflow (JIC) configuration. The process starts with a cubic eddy viscosity model (CEVM) developed for incompressible flows. It is fitted to limited experimental JIC data using shrinkage regression. The shrinkage process removes all the terms from the model, except an intercept, a linear term, and a quadratic one involving the square of the vorticity. The shrunk eddy viscositymore » model is implemented in an RANS simulator and calibrated, using vorticity measurements, to infer three parameters. The calibration is Bayesian and is solved using a Markov chain Monte Carlo (MCMC) method. A 3D probability density distribution for the inferred parameters is constructed, thus quantifying the uncertainty in the estimate. The phenomenal cost of using a 3D flow simulator inside an MCMC loop is mitigated by using surrogate models (“curve-fits”). A support vector machine classifier (SVMC) is used to impose our prior belief regarding parameter values, specifically to exclude nonphysical parameter combinations. The calibrated model is compared, in terms of its predictive skill, to simulations using uncalibrated linear and CEVMs. Finally, we find that the calibrated model, with one quadratic term, is more accurate than the uncalibrated simulator. The model is also checked at a flow condition at which the model was not calibrated.« 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.

Numerical modeling of turbulent swirling flow in a multi-inlet vortex nanoprecipitation reactor using dynamic DDES

NASA Astrophysics Data System (ADS)

Hill, James C.; Liu, Zhenping; Fox, Rodney O.; Passalacqua, Alberto; Olsen, Michael G.

2015-11-01

The multi-inlet vortex reactor (MIVR) has been developed to provide a platform for rapid mixing in the application of flash nanoprecipitation (FNP) for manufacturing functional nanoparticles. Unfortunately, commonly used RANS methods are unable to accurately model this complex swirling flow. Large eddy simulations have also been problematic, as expensive fine grids to accurately model the flow are required. These dilemmas led to the strategy of applying a Delayed Detached Eddy Simulation (DDES) method to the vortex reactor. In the current work, the turbulent swirling flow inside a scaled-up MIVR has been investigated by using a dynamic DDES model. In the DDES model, the eddy viscosity has a form similar to the Smagorinsky sub-grid viscosity in LES and allows the implementation of a dynamic procedure to determine its coefficient. The complex recirculating back flow near the reactor center has been successfully captured by using this dynamic DDES model. Moreover, the simulation results are found to agree with experimental data for mean velocity and Reynolds stresses.

Distribution of the near-inertial kinetic energy inside mesoscale eddies: Observations in the Gulf of Mexico

NASA Astrophysics Data System (ADS)

Ixetl Garcia Gomez, Beatriz; Pallas Sanz, Enric; Candela Perez, Julio

2017-04-01

The near-inertial oscillations (NIOs), generated by the wind stress on the surface mixed layer, are the inertia gravity waves with the lowest frequency and the highest kinetic energy. NIOs are important because they drive vertical mixing in the interior ocean during wave breaking events. Although the interaction between NIOs and mesoscale eddies has been reported by several authors, these studies are mostly analytical and numerical, and only few observational studies have attempted to show the differences in near-inertial kinetic energy (KEi) between anticyclonic and cyclonic eddies. In this work the spatial structure of the KEi inside the mesoscale eddies is computed using daily satellite altimetry and observations of horizontal velocity from 23 moorings equipped with acoustic Doppler current profilers in the western Gulf of Mexico. Consistent to theory, the obtained four-year KEi-composites show two times more KEi inside the anticyclonic eddies than inside the cyclonic ones. The vertical and horizontal cross-sections of the KEi-composites show that the KEi is mainly located near to the surface of the cyclonic eddies (positive vorticity), whereas the KEi in anticyclonic eddies (negative vorticity) is maximum in the eddy's center near to the base of the eddy where the NIOs become more inertial, are trapped, and amplified. The mean vertical profiles show that the cyclonic eddies present a maximum of KEi near to the surface at 50, while the maximum of KEi in the anticyclonic eddies occurs between 900 and 1100 m. Inside anticyclonic eddies another two relative maximums are observed, one in the mixed layer and the second at 300 m. In contrast, the mean profile of KEi outside the mesoscale eddies has the maximum value at the surface ( 50 m), with high values of KEi in the first 200 m and negligible energy beneath that depth. A different mean distribution of the KEi is observed depending on the type of wind generator: tropical storms or unidirectional wind.

Calculation of eddy viscosity in a compressible turbulent boundary layer with mass injection and chemical reaction, volume 2. [computer programs

NASA Technical Reports Server (NTRS)

Omori, S.

1973-01-01

As described in Vol. 1, the eddy viscosity is calculated through the turbulent kinetic energy, in order to include the history of the flow and the effect of chemical reaction on boundary layer characteristics. Calculations can be performed for two different cooling concepts; that is, transpiration and regeneratively cooled wall cases. For the regenerative cooling option, coolant and gas side wall temperature and coolant bulk temperature in a rocket engine can be computed along the nozzle axis. Thus, this computer program is useful in designing coolant flow rate and cooling tube geometry, including the tube wall thickness as well as in predicting the effects of boundary layers along the gas side wall on thrust performances.

Characteristics of subgrid-resolved-scale dynamics in anisotropic turbulence, with application to rough-wall boundary layers

NASA Astrophysics Data System (ADS)

Juneja, Anurag; Brasseur, James G.

1999-10-01

Large-eddy simulation (LES) of the atmospheric boundary layer (ABL) using eddy viscosity subgrid-scale (SGS) models is known to poorly predict mean shear at the first few grid cells near the ground, a rough surface with no viscous sublayer. It has recently been shown that convective motions carry this localized error vertically to infect the entire ABL, and that the error is more a consequence of the SGS model than grid resolution in the near-surface inertial layer. Our goal was to determine what first-order errors in the predicted SGS terms lead to spurious expectation values, and what basic dynamics in the filtered equation for resolved scale (RS) velocity must be captured by SGS models to correct the deficiencies. Our analysis is of general relevance to LES of rough-wall high Reynolds number boundary layers, where the essential difficulty in the closure is the importance of the SGS acceleration terms, a consequence of necessary under-resolution of relevant energy-containing motions at the first few grid levels, leading to potentially strong couplings between the anisotropies in resolved velocity and predicted SGS dynamics. We analyze these two issues (under-resolution and anisotropy) in the absence of a wall using two direct numerical simulation datasets of homogeneous turbulence with very different anisotropic structure characteristic of the near-surface ABL: shear- and buoyancy-generated turbulence. We uncover three important issues which should be addressed in the design of SGS closures near rough walls and we provide a priori tests for the SGS model. First, we identify a strong spurious coupling between the anisotropic structure of the resolved velocity field and predicted SGS dynamics which can create a feedback loop to incorrectly enhance certain components of the predicted velocity field. Second, we find that eddy viscosity and "similarity" SGS models do not contain enough degrees of freedom to capture, at a sufficient level of accuracy, both RS-SGS energy flux and SGS-RS dynamics. Third, to correctly capture pressure transport near a wall, closures must be made more flexible to accommodate proper partitioning between SGS stress divergence and SGS pressure gradient.

Calculation of eddy viscosity in a compressible turbulent boundary layer with mass injection and chemical reaction, volume 1. [theoretical analysis

NASA Technical Reports Server (NTRS)

Omori, S.

1973-01-01

The turbulent kinetic energy equation is coupled with boundary layer equations to solve the characteristics of compressible turbulent boundary layers with mass injection and combustion. The Reynolds stress is related to the turbulent kinetic energy using the Prandtl-Wieghardt formulation. When a lean mixture of hydrogen and nitrogen is injected through a porous plate into the subsonic turbulent boundary layer of air flow and ignited by external means, the turbulent kinetic energy increases twice as much as that of noncombusting flow with the same mass injection rate of nitrogen. The magnitudes of eddy viscosity between combusting and noncombusting flows with injection, however, are almost the same due to temperature effects, while the distributions are different. The velocity profiles are significantly affected by combustion; that is, combustion alters the velocity profile as if the mass injection rate is increased, reducing the skin-friction as a result of a smaller velocity gradient at the wall. If pure hydrogen as a transpiration coolant is injected into a rocket nozzle boundary layer flow of combustion products, the temperature drops significantly across the boundary layer due to the high heat capacity of hydrogen. At a certain distance from the wall, hydrogen reacts with the combustion products, liberating an extensive amount of heat. The resulting large increase in temperature reduces the eddy viscosity in this region.

Overflow Simulations using MPAS-Ocean in Idealized and Realistic Domains

NASA Astrophysics Data System (ADS)

Reckinger, S.; Petersen, M. R.; Reckinger, S. J.

2016-02-01

MPAS-Ocean is used to simulate an idealized, density-driven overflow using the dynamics of overflow mixing and entrainment (DOME) setup. Numerical simulations are benchmarked against other models, including the MITgcm's z-coordinate model and HIM's isopycnal coordinate model. A full parameter study is presented that looks at how sensitive overflow simulations are to vertical grid type, resolution, and viscosity. Horizontal resolutions with 50 km grid cells are under-resolved and produce poor results, regardless of other parameter settings. Vertical grids ranging in thickness from 15 m to 120 m were tested. A horizontal resolution of 10 km and a vertical resolution of 60 m are sufficient to resolve the mesoscale dynamics of the DOME configuration, which mimics real-world overflow parameters. Mixing and final buoyancy are least sensitive to horizontal viscosity, but strongly sensitive to vertical viscosity. This suggests that vertical viscosity could be adjusted in overflow water formation regions to influence mixing and product water characteristics. Also, the study shows that sigma coordinates produce much less mixing than z-type coordinates, resulting in heavier plumes that go further down slope. Sigma coordinates are less sensitive to changes in resolution but as sensitive to vertical viscosity compared to z-coordinates. Additionally, preliminary measurements of overflow diagnostics on global simulations using a realistic oceanic domain are presented.

Prediction of the rate of the rise of an air bubble in nanofluids in a vertical tube.

PubMed

Cho, Heon Ki; Nikolov, Alex D; Wasan, Darsh T

2018-04-19

Our recent experiments have demonstrated that when a bubble rises through a nanofluid (a liquid containing dispersed nanoparticles) in a vertical tube, a nanofluidic film with several particle layers is formed between the gas bubble and the glass tube wall, which significantly changes the bubble velocity due to the nanoparticle layering phenomenon in the film. We calculated the structural nanofilm viscosity as a function of the number of particle layers confined in it and found that the film viscosity increases rather steeply when the film contains only one or two particle layers. The nanofilm viscosity was found to be several times higher than the bulk viscosity of the fluid. Consequently, the Bretherton equation cannot accurately predict the rate of the rise of a slow-moving long bubble in a vertical tube in a nanofluid because it is valid only for very thick films and uses the bulk viscosity of the fluid. However, in this brief note, we demonstrate that the Bretherton equation can indeed be used for predicting the rate of the rise of a long single bubble through a vertical tube filled with a nanofluid by simply replacing the bulk viscosity with the proper structural nanofilm viscosity of the fluid. Copyright © 2018. Published by Elsevier Inc.

Inference and Biogeochemical Response of Vertical Velocities inside a Mode Water Eddy

NASA Astrophysics Data System (ADS)

Barceló-Llull, B.; Pallas Sanz, E.; Sangrà, P.

2016-02-01

With the aim to study the modulation of the biogeochemical fluxes by the ageostrophic secondary circulation in anticyclonic mesoscale eddies, a typical eddy of the Canary Eddy Corridor was interdisciplinary surveyed on September 2014 in the framework of the PUMP project. The eddy was elliptical shaped, 4 month old, 110 km diameter and 400 m depth. It was an intrathermocline type often also referred as mode water eddy type. We inferred the mesoscale vertical velocity field resolving a generalized omega equation from the 3D density and ADCP velocity fields of a five-day sampled CTD-SeaSoar regular grid centred on the eddy. The grid transects where 10 nautical miles apart. Although complex, in average, the inferred omega velocity field (hereafter w) shows a dipolar structure with downwelling velocities upstream of the propagation path (west) and upwelling velocities downstream. The w at the eddy center was zero and maximum values were located at the periphery attaining ca. 6 m day-1. Coinciding with the occurrence of the vertical velocities cells a noticeable enhancement of phytoplankton biomass was observed at the eddy periphery respect to the far field. A corresponding upward diapycnal flux of nutrients was also observed at the periphery. As minimum velocities where reached at the eddy center, lineal Ekman pumping mechanism was discarded. Minimum values of phytoplankton biomass where also observed at the eddy center. The possible mechanisms for such dipolar w cell are still being investigated, but an analysis of the generalized omega equation forcing terms suggest that it may be a combination of horizontal deformation and advection of vorticity by the ageostrophic current (related to nonlinear Ekman pumping). As expected for Trades, the wind was rather constant and uniform with a speed of ca. 5 m s-1. Diagnosed nonlinear Ekman pumping leaded also to a dipolar cell that mirrors the omega w dipolar cell.

Study of the far wake vortex field generated by a rectangular airfoil in a water tank

NASA Technical Reports Server (NTRS)

Lezius, D. K.

1973-01-01

Underwater towing experiments were carried out with a rectangular airfoil of aspect ratio 5.3 at 4 and 8 deg angles of attack and at chord-based Reynolds numbers between 2 x 100,000 and 7.5 x 100,000. Quantitative measurements by means of the hydrogen bubble technique indicated lower peak swirl velocities in the range of 100 to 1000 lenghts downstream than have been measured in wind tunnel of flight tests. The maximum circumferential velocity decayed whereas the turbulent eddy viscosity increased. This behavior and other known rates of vortex decay are explained in terms of an analytical solution for the vortex problem with time varying eddy viscosity. It is shown that this case corresponds to nonequilibrium turbulent vortex flow.

An Improved K-Epsilon Model for Near-Wall Turbulence and Comparison with Direct Numerical Simulation

NASA Technical Reports Server (NTRS)

Shih, T. H.

1990-01-01

An improved k-epsilon model for low Reynolds number turbulence near a wall is presented. The near-wall asymptotic behavior of the eddy viscosity and the pressure transport term in the turbulent kinetic energy equation is analyzed. Based on this analysis, a modified eddy viscosity model, having correct near-wall behavior, is suggested, and a model for the pressure transport term in the k-equation is proposed. In addition, a modeled dissipation rate equation is reformulated. Fully developed channel flows were used for model testing. The calculations using various k-epsilon models are compared with direct numerical simulations. The results show that the present k-epsilon model performs well in predicting the behavior of near-wall turbulence. Significant improvement over previous k-epsilon models is obtained.

PRATHAM: Parallel Thermal Hydraulics Simulations using Advanced Mesoscopic Methods

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

Joshi, Abhijit S; Jain, Prashant K; Mudrich, Jaime A

2012-01-01

At the Oak Ridge National Laboratory, efforts are under way to develop a 3D, parallel LBM code called PRATHAM (PaRAllel Thermal Hydraulic simulations using Advanced Mesoscopic Methods) to demonstrate the accuracy and scalability of LBM for turbulent flow simulations in nuclear applications. The code has been developed using FORTRAN-90, and parallelized using the message passing interface MPI library. Silo library is used to compact and write the data files, and VisIt visualization software is used to post-process the simulation data in parallel. Both the single relaxation time (SRT) and multi relaxation time (MRT) LBM schemes have been implemented in PRATHAM.more » To capture turbulence without prohibitively increasing the grid resolution requirements, an LES approach [5] is adopted allowing large scale eddies to be numerically resolved while modeling the smaller (subgrid) eddies. In this work, a Smagorinsky model has been used, which modifies the fluid viscosity by an additional eddy viscosity depending on the magnitude of the rate-of-strain tensor. In LBM, this is achieved by locally varying the relaxation time of the fluid.« less

[CO2 turbulent exchange in a broadleaved Korean pine forest in Changbai Mountains].

PubMed

Wu, Jia-bing; Guan, De-xin; Sun, Xiao-min; Shi, Ting-ting; Han, Shi-jie; Jin, Chang-jie

2007-05-01

The measurement of CO2 turbulent exchange in a broadleaved Korean pine forest in Changbai Mountains by an open-path eddy covariance system showed that with near neutral atmospheric stratification, the CO2 and vertical wind components over canopy in inertial subrange followed the expected -2/3 power law, and the dominant vertical eddy scale was about 40 m. The frequency ranges of eddy contributions to CO2 fluxes were mostly within 0.01-2.0 Hz, and the eddy translated by low frequency over canopy contributed more of CO2 fluxes. The open-path eddy covariance system could satisfy the estimation of turbulent fluxes over canopy, but the CO2 fluxes between forest and atmosphere were generally underestimated at night because the increment of non turbulent processes, suggesting that the CO2 fluxes estimated under weak turbulence needed to revise correspondingly.

Coherent mesoscale eddies in the North Atlantic subtropical gyre: 3-D structure and transport with application to the salinity maximum

NASA Astrophysics Data System (ADS)

Amores, Angel; Melnichenko, Oleg; Maximenko, Nikolai

2017-01-01

The mean vertical structure and transport properties of mesoscale eddies are investigated in the North Atlantic subtropical gyre by combining historical records of Argo temperature/salinity profiles and satellite sea level anomaly data in the framework of the eddy tracking technique. The study area is characterized by a low eddy kinetic energy and sea surface salinity maximum. Although eddies have a relatively weak signal at surface (amplitudes around 3-7 cm), the eddy composites reveal a clear deep signal that penetrates down to at least 1200 m depth. The analysis also reveals that the vertical structure of the eddy composites is strongly affected by the background stratification. The horizontal patterns of temperature/salinity anomalies can be reconstructed by a linear combination of a monopole, related to the elevation/depression of the isopycnals in the eddy core, and a dipole, associated with the horizontal advection of the background gradient by the eddy rotation. A common feature of all the eddy composites reconstructed is the phase coherence between the eddy temperature/salinity and velocity anomalies in the upper ˜300 m layer, resulting in the transient eddy transports of heat and salt. As an application, a box model of the near-surface layer is used to estimate the role of mesoscale eddies in maintaining a quasi-steady state distribution of salinity in the North Atlantic subtropical salinity maximum. The results show that mesoscale eddies are able to provide between 4 and 21% of the salt flux out of the area required to compensate for the local excess of evaporation over precipitation.

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

The role of bulk viscosity on the decay of compressible, homogeneous, isotropic turbulence

NASA Astrophysics Data System (ADS)

Johnsen, Eric; Pan, Shaowu

2016-11-01

The practice of neglecting bulk viscosity in studies of compressible turbulence is widespread. While exact for monatomic gases and unlikely to strongly affect the dynamics of fluids whose bulk-to-shear viscosity ratio is small and/or of weakly compressible turbulence, this assumption is not justifiable for compressible, turbulent flows of gases whose bulk viscosity is orders of magnitude larger than their shear viscosities (e.g., CO2). To understand the mechanisms by which bulk viscosity and the associated phenomena affect compressible turbulence, we conduct DNS of freely decaying compressible, homogeneous, isotropic turbulence for ratios of bulk-to-shear viscosity ranging from 0-1000. Our simulations demonstrate that bulk viscosity increases the decay rate of turbulent kinetic energy; while enstrophy exhibits little sensitivity to bulk viscosity, dilatation is reduced by an order of magnitude within the two eddy turnover time. Via a Helmholtz decomposition of the flow, we determined that bulk viscosity damps the dilatational velocity and reduces dilatational-solenoidal exchanges, as well as pressure-dilatation coupling. In short, bulk viscosity renders compressible turbulence incompressible by reducing energy transfer between translational and internal modes.

RAS one-equation turbulence model with non-singular eddy-viscosity coefficient

NASA Astrophysics Data System (ADS)

Rahman, M. M.; Agarwal, R. K.; Siikonen, T.

2016-02-01

A simplified consistency formulation for Pk/ε (production to dissipation ratio) is devised to obtain a non-singular Cμ (coefficient of eddy-viscosity) in the explicit algebraic Reynolds stress model of Gatski and Speziale. The coefficient Cμ depends non-linearly on both rotational/irrotational strains and is used in the framework of an improved RAS (Rahman-Agarwal-Siikonen) one-equation turbulence model to calculate a few well-documented turbulent flows, yielding predictions in good agreement with the direct numerical simulation and experimental data. The strain-dependent Cμ assists the RAS model in constructing the coefficients and functions such as to benefit complex flows with non-equilibrium turbulence. Comparisons with the Spalart-Allmaras one-equation model and the shear stress transport k-ω model demonstrate that Cμ improves the response of RAS model to non-equilibrium effects.

Subduction at upper ocean fronts by baroclinic instability

NASA Astrophysics Data System (ADS)

Verma, Vicky; Pham, Hieu T.; Radhakrishnan, Anand; Sarkar, Sutanu

2017-11-01

Large eddy simulations of upper ocean fronts that are initially in geostrophic balance show that the linear and subsequent nonlinear evolution of baroclinic intability are effective in restratifying the front. During the growth of baroclinic instability, the front develops thin regions with enhanced vertical vorticity, i.e., vorticity filaments. Moreover, the vorticity filaments organize into submesoscale eddies. The subsequent frontal dynamics is dominated by the vorticity filaments and the submesoscale eddies. Diagnosis of the horizontal force balance reveals that the regions occupied by these coherent structures have significantly large imbalance, and are characterized by large vertical velocity. High density fluid from the heavier side of the front is subducted by the vertical velocity to the bottom of the mixed layer. The process of subduction is illustrated by Lagrangian tracking of fluid particles released at a fixed depth.

Observations and analysis of a stratification-destratification event in a tropical estuary

NASA Astrophysics Data System (ADS)

Uncles, R. J.; Ong, J. E.; Gong, W. K.

1990-11-01

A data set comprising 31 continuous tidal cycles was collected in the Sungai Merbok Estuary, Malaysia, in June 1987 as part of an ecological study of nutrient fluxes from a tropical mangrove estuary. Currents, salinity and salinity stratification at a deep-channel (15 m) station near the mouth of the Merbok Estuary showed a pronounced spring-neap variability. The slow currents and weak vertical mixing at neap tides favoured the formation of a stratified water column and generated a neap-spring cycle of water column stabilization and destabilization. A strong stratification event occurred during the period of observations. This was partly driven by a modest freshwater spate which coincided with neap tides. An eddy viscosity-diffusivity model of the stratification, which assumed a constant, longitudinal salinity gradient, demonstrated a pronounced stratification-destratification cycle due to neap-spring variations in vertical mixing. Larger and more realistic stratification was modelled when the estimated, time-varying longitudinal salinity gradient was incorporated. This gradient maximized in response to the peak in freshwater runoff. The measured and modelled density-driven circulations showed qualitative similarities and were of the order of 10 cm s -1 at neap tides. The circulation was weaker during spring tides. The tidally-filtered salt transport due to vertical shear was directed up-estuary and was an order of magnitude smaller during spring tides. The results are discussed in terms of their relevance to mangrove system oceanography.

right-sized dimple evaluator

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

Rodriguez, Sal

2017-08-24

The code (aka computer program written as a Matlab script) uses a unique set of n independent equations to solve for n turbulence variables. The code requires the input of a characteristic dimension, a characteristic fluid velocity, the fluid dynamic viscosity, and the fluid density. Most importantly, the code estimates the size of three key turbulent eddies: Kolmogorov, Taylor, and integral. Based on the eddy sizes, dimples dimensions are prescribed such that the key eddies (principally Taylor, and sometimes Kolmogorov), can be generated by the dimple rim and flow unimpeded through the dimple’s concave cavity. It is hypothesized that turbulentmore » eddies are generated by the dimple rim at the dimple-surface interface. The newly-generated eddies in turn entrain the movement of surrounding regions of fluid, creating more mixing. The eddies also generate lift near the wall surrounding the dimple, as they accelerate and reduce pressure in the regions near and at the dimple cavity, thereby minimizing the fluid drag.« less

Impacts of the Mesoscale Ocean-Atmosphere Coupling on the Peru-Chile Ocean Dynamics: The Current-Induced Wind Stress Modulation

NASA Astrophysics Data System (ADS)

Oerder, V.; Colas, F.; Echevin, V.; Masson, S.; Lemarié, F.

2018-02-01

The ocean dynamical responses to the surface current-wind stress interaction at the oceanic mesoscale are investigated in the South-East Pacific using a high-resolution regional ocean-atmosphere coupled model. Two simulations are compared: one includes the surface current in the wind stress computation while the other does not. In the coastal region, absolute wind velocities are different between the two simulations but the wind stress remains very similar. As a consequence, the mean regional oceanic circulation is almost unchanged. On the contrary, the mesoscale activity is strongly reduced when taking into account the effect of the surface current on the wind stress. This is caused by a weakening of the eddy kinetic energy generation near the coast by the wind work and to intensified offshore eddy damping. We show that, above coherent eddies, the current-stress interaction generates eddy damping through Ekman pumping and eddy kinetic energy dissipation through wind work. This alters significantly the coherent eddy vertical structures compared with the control simulation, weakening the temperature and vorticity anomalies and increasing strongly the vertical velocity anomalies associated to eddies.

Very-high-Reynolds-number vortex dynamics via Coherent-vorticity-Preserving (CvP) Large-eddy simulations

NASA Astrophysics Data System (ADS)

Chapelier, Jean-Baptiste; Wasistho, Bono; Scalo, Carlo

2017-11-01

A new approach to Large-Eddy Simulation (LES) is introduced, where subgrid-scale (SGS) dissipation is applied proportionally to the degree of local spectral broadening, hence mitigated in regions dominated by large-scale vortical motion. The proposed CvP-LES methodology is based on the evaluation of the ratio of the test-filtered to resolved (or grid-filtered) enstrophy: σ = ξ ∧ / ξ . Values of σ = 1 indicate low sub-test-filter turbulent activity, justifying local deactivation of any subgrid-scale model. Values of σ < 1 span conditions ranging from incipient spectral broadening σ <= 1 , to equilibrium turbulence σ =σeq < 1 , where σeq is solely as a function of the test-to-grid filter-width ratio Δ ∧ / Δ , derived assuming a Kolmogorov's spectrum. Eddy viscosity is fully restored for σ <=σeq . The proposed approach removes unnecessary SGS dissipation, can be applied to any eddy-viscosity model, is algorithmically simple and computationally inexpensive. A CvP-LES of a pair of unstable helical vortices, representative of rotor-blade wake dynamics, show the ability of the method to sort the coherent motion from the small-scale dynamics. This work is funded by subcontract KSC-17-001 between Purdue University and Kord Technologies, Inc (Huntsville), under the US Navy Contract N68335-17-C-0159 STTR-Phase II, Purdue Proposal No. 00065007, Topic N15A-T002.

Interactions between a fractal tree-like object and hydrodynamic turbulence: flow structure and characteristic mixing length

NASA Astrophysics Data System (ADS)

Meneveau, C. V.; Bai, K.; Katz, J.

2011-12-01

The vegetation canopy has a significant impact on various physical and biological processes such as forest microclimate, rainfall evaporation distribution and climate change. Most scaled laboratory experimental studies have used canopy element models that consist of rigid vertical strips or cylindrical rods that can be typically represented through only one or a few characteristic length scales, for example the diameter and height for cylindrical rods. However, most natural canopies and vegetation are highly multi-scale with branches and sub-branches, covering a wide range of length scales. Fractals provide a convenient idealization of multi-scale objects, since their multi-scale properties can be described in simple ways (Mandelbrot 1982). While fractal aspects of turbulence have been studied in several works in the past decades, research on turbulence generated by fractal objects started more recently. We present an experimental study of boundary layer flow over fractal tree-like objects. Detailed Particle-Image-Velocimetry (PIV) measurements are carried out in the near-wake of a fractal-like tree. The tree is a pre-fractal with five generations, with three branches and a scale reduction factor 1/2 at each generation. Its similarity fractal dimension (Mandelbrot 1982) is D ~ 1.58. Detailed mean velocity and turbulence stress profiles are documented, as well as their downstream development. We then turn attention to the turbulence mixing properties of the flow, specifically to the question whether a mixing length-scale can be identified in this flow, and if so, how it relates to the geometric length-scales in the pre-fractal object. Scatter plots of mean velocity gradient (shear) and Reynolds shear stress exhibit good linear relation at all locations in the flow. Therefore, in the transverse direction of the wake evolution, the Boussinesq eddy viscosity concept is appropriate to describe the mixing. We find that the measured mixing length increases with increasing streamwise locations. Conversely, the measured eddy viscosity and mixing length decrease with increasing elevation, which differs from eddy viscosity and mixing length behaviors of traditional boundary layers or canopies studied before. In order to find an appropriate length for the flow, several models based on the notion of superposition of scales are proposed and examined. One approach is based on spectral distributions. Another more practical approach is based on length-scale distributions evaluated using fractal geometry tools. These proposed models agree well with the measured mixing length. The results indicate that information about multi-scale clustering of branches as it occurs in fractals has to be incorporated into models of the mixing length for flows through canopies with multiple scales. The research is supported by National Science Foundation grant ATM-0621396 and AGS-1047550.

Foam-machining tool with eddy-current transducer

NASA Technical Reports Server (NTRS)

Copper, W. P.

1975-01-01

Three-cutter machining system for foam-covered tanks incorporates eddy-current sensor. Sensor feeds signal to numerical controller which programs rotational and vertical axes of sensor travel, enabling cutterhead to profile around tank protrusions.

Eddy-driven nutrient transport and associated upper-ocean primary production along the Kuroshio

NASA Astrophysics Data System (ADS)

Uchiyama, Yusuke; Suzue, Yota; Yamazaki, Hidekatsu

2017-06-01

The Kuroshio is one of the most energetic western boundary currents accompanied by vigorous eddy activity both on mesoscale and submesoscale, which affects biogeochemical processes in the upper ocean. We examine the primary production around the Kuroshio off Japan using a climatological ocean modeling based on the Regional Oceanic Modeling System (ROMS) coupled with a nitrogen-based nutrient, phytoplankton and zooplankton, and detritus (NPZD) biogeochemical model in a submesoscale eddy-permitting configuration. The model indicates significant differences of the biogeochemical responses to eddy activities in the Kuroshio Region (KR) and Kuroshio Extension Region (KE). In the KR, persisting cyclonic eddies developed between the Kuroshio and coastline are responsible for upwelling-induced eutrophication. However, the eddy-induced vertical nutrient flux counteracts and promotes pronounced southward and downward diapycnal nutrient transport from the mixed-layer down beneath the main body of the Kuroshio, which suppresses the near-surface productivity. In contrast, the KE has a 23.5% higher productivity than the KR, even at comparable eddy intensity. Upward nutrient transport prevails near the surface due to predominant cyclonic eddies, particularly to the north of the KE, where the downward transport barely occurs, except at depths deeper than 400 m and to a much smaller degree than in the KR. The eddy energy conversion analysis reveals that the combination of shear instability around the mainstream of the Kuroshio with prominent baroclinic instability near the Kuroshio front is essential for the generation of eddies in the KR, leading to the increase of the eddy-induced vertical nitrate transport around the Kuroshio.

Tropical waves and the quasi-biennial oscillation in the lower stratosphere

NASA Technical Reports Server (NTRS)

Miller, A. J.; Angell, J. K.; Korshover, J.

1976-01-01

By means of spectrum analysis of 11 years of lower stratospheric daily winds and temperatures at Balboa, Ascension and Canton-Singapore, evidence is presented supporting the existence of two principal wave modes with periods of about 11-17 days (Kelvin waves) and about 4-5 days (mixed Rossby-gravity waves). The structure of the two wave modes, as well as the vertical eddy momentum flux by the waves, is shown to be related to the quasi-biennial cycle, although for the mixed Rossby-gravity waves this is obvious only at Ascension. In addition, the Coriolis term, suggested as a source of vertical easterly momentum flux for the mixed Rossby-gravity waves, is investigated and found to be of the same magnitude as the vertical eddy flux term. Finally, we have examined the mean meridional motion and the meridional eddy momentum flux for its possible association with the quasi- biennial variation.

Polar symmetric flow of a viscous compressible atmosphere; an application to Mars

NASA Technical Reports Server (NTRS)

Pirraglia, J. A.

1974-01-01

The atmosphere is assumed to be driven by a polar symmetric temperature field and the equations of motion in pressure ratio coordinates are linearized by considering the zero order in terms of a thermal Rossby number R delta I/(2a omega) sq where delta T is a measure of the latitudinal temperature gradient. When the eddy viscosity is greater than 1 million sq cm/sec, the boundary layer extends far up into the atmosphere, making the geostrophic approximation invalid for the bulk of the atmosphere. A temperature model for Mars was used which was based on Mariner 9 infrared spectral data with a 30% increase in the depth averaged temperature from the winter pole to the subsolar point. The results obtained for the increase in surface pressure from the subsolar point to the winter pole, as a function of eddy viscosity and with no-slip conditions imposed at the surface, are given.

Reference Solutions for Benchmark Turbulent Flows in Three Dimensions

NASA Technical Reports Server (NTRS)

Diskin, Boris; Thomas, James L.; Pandya, Mohagna J.; Rumsey, Christopher L.

2016-01-01

A grid convergence study is performed to establish benchmark solutions for turbulent flows in three dimensions (3D) in support of turbulence-model verification campaign at the Turbulence Modeling Resource (TMR) website. The three benchmark cases are subsonic flows around a 3D bump and a hemisphere-cylinder configuration and a supersonic internal flow through a square duct. Reference solutions are computed for Reynolds Averaged Navier Stokes equations with the Spalart-Allmaras turbulence model using a linear eddy-viscosity model for the external flows and a nonlinear eddy-viscosity model based on a quadratic constitutive relation for the internal flow. The study involves three widely-used practical computational fluid dynamics codes developed and supported at NASA Langley Research Center: FUN3D, USM3D, and CFL3D. Reference steady-state solutions computed with these three codes on families of consistently refined grids are presented. Grid-to-grid and code-to-code variations are described in detail.

Reynolds averaged turbulence modelling using deep neural networks with embedded invariance

DOE PAGES

Ling, Julia; Kurzawski, Andrew; Templeton, Jeremy

2016-10-18

There exists significant demand for improved Reynolds-averaged Navier–Stokes (RANS) turbulence models that are informed by and can represent a richer set of turbulence physics. This paper presents a method of using deep neural networks to learn a model for the Reynolds stress anisotropy tensor from high-fidelity simulation data. A novel neural network architecture is proposed which uses a multiplicative layer with an invariant tensor basis to embed Galilean invariance into the predicted anisotropy tensor. It is demonstrated that this neural network architecture provides improved prediction accuracy compared with a generic neural network architecture that does not embed this invariance property.more » Furthermore, the Reynolds stress anisotropy predictions of this invariant neural network are propagated through to the velocity field for two test cases. For both test cases, significant improvement versus baseline RANS linear eddy viscosity and nonlinear eddy viscosity models is demonstrated.« less

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

NASA Astrophysics Data System (ADS)

Cook, Andrew W.

2007-05-01

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

Wave attenuation in the marginal ice zone during LIMEX

NASA Technical Reports Server (NTRS)

Liu, Antony K.; Vachon, Paris W.; Peng, Chih Y.; Bhogal, A. S.

1992-01-01

The effect of ice cover on ocean-wave attenuation is investigated for waves under flexure in the marginal ice zone (MIZ) with SAR image spectra and the results of models. Directional wavenumber spectra are taken from the SAR image data, and the wave-attenuation rate is evaluated with SAR image spectra and by means of the model by Liu and Mollo-Christensen (1988). Eddy viscosity is described by means of dimensional analysis as a function of ice roughness and wave-induced velocity, and comparisons are made with the remotely sensed data. The model corrects the open-water model by introducing the effects of a continuous ice sheet, and turbulent eddy viscosity is shown to depend on ice thickness, floe sizes, significant wave height, and wave period. SAR and wave-buoy data support the trends described in the model results, and a characteristic rollover is noted in the model and experimental wave-attenuation rates at high wavenumbers.

Reynolds averaged turbulence modelling using deep neural networks with embedded invariance

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

Ling, Julia; Kurzawski, Andrew; Templeton, Jeremy

There exists significant demand for improved Reynolds-averaged Navier–Stokes (RANS) turbulence models that are informed by and can represent a richer set of turbulence physics. This paper presents a method of using deep neural networks to learn a model for the Reynolds stress anisotropy tensor from high-fidelity simulation data. A novel neural network architecture is proposed which uses a multiplicative layer with an invariant tensor basis to embed Galilean invariance into the predicted anisotropy tensor. It is demonstrated that this neural network architecture provides improved prediction accuracy compared with a generic neural network architecture that does not embed this invariance property.more » Furthermore, the Reynolds stress anisotropy predictions of this invariant neural network are propagated through to the velocity field for two test cases. For both test cases, significant improvement versus baseline RANS linear eddy viscosity and nonlinear eddy viscosity models is demonstrated.« less

Computation of turbulent rotating channel flow with an algebraic Reynolds stress model

NASA Technical Reports Server (NTRS)

Warfield, M. J.; Lakshminarayana, B.

1986-01-01

An Algebraic Reynolds Stress Model has been implemented to modify the Kolmogorov-Prandtl eddy viscosity relation to produce an anisotropic turbulence model. The eddy viscosity relation becomes a function of the local turbulent production to dissipation ratio and local turbulence/rotation parameters. The model is used to predict fully-developed rotating channel flow over a diverse range of rotation numbers. In addition, predictions are obtained for a developing channel flow with high rotation. The predictions are compared with the experimental data available. Good predictions are achieved for mean velocity and wall shear stress over most of the rotation speeds tested. There is some prediction breakdown at high rotation (rotation number greater than .10) where the effects of the rotation on turbulence become quite complex. At high rotation and low Reynolds number, the laminarization on the trailing side represents a complex effect of rotation which is difficult to predict with the described models.

Prediction of High-Lift Flows using Turbulent Closure Models

NASA Technical Reports Server (NTRS)

Rumsey, Christopher L.; Gatski, Thomas B.; Ying, Susan X.; Bertelrud, Arild

1997-01-01

The flow over two different multi-element airfoil configurations is computed using linear eddy viscosity turbulence models and a nonlinear explicit algebraic stress model. A subset of recently-measured transition locations using hot film on a McDonnell Douglas configuration is presented, and the effect of transition location on the computed solutions is explored. Deficiencies in wake profile computations are found to be attributable in large part to poor boundary layer prediction on the generating element, and not necessarily inadequate turbulence modeling in the wake. Using measured transition locations for the main element improves the prediction of its boundary layer thickness, skin friction, and wake profile shape. However, using measured transition locations on the slat still yields poor slat wake predictions. The computation of the slat flow field represents a key roadblock to successful predictions of multi-element flows. In general, the nonlinear explicit algebraic stress turbulence model gives very similar results to the linear eddy viscosity models.

A simple method for simulating wind profiles in the boundary layer of tropical cyclones

DOE PAGES

Bryan, George H.; Worsnop, Rochelle P.; Lundquist, Julie K.; ...

2016-11-01

A method to simulate characteristics of wind speed in the boundary layer of tropical cyclones in an idealized manner is developed and evaluated. The method can be used in a single-column modelling set-up with a planetary boundary-layer parametrization, or within large-eddy simulations (LES). The key step is to include terms in the horizontal velocity equations representing advection and centrifugal acceleration in tropical cyclones that occurs on scales larger than the domain size. Compared to other recently developed methods, which require two input parameters (a reference wind speed, and radius from the centre of a tropical cyclone) this new method alsomore » requires a third input parameter: the radial gradient of reference wind speed. With the new method, simulated wind profiles are similar to composite profiles from dropsonde observations; in contrast, a classic Ekman-type method tends to overpredict inflow-layer depth and magnitude, and two recently developed methods for tropical cyclone environments tend to overpredict near-surface wind speed. When used in LES, the new technique produces vertical profiles of total turbulent stress and estimated eddy viscosity that are similar to values determined from low-level aircraft flights in tropical cyclones. Lastly, temporal spectra from LES produce an inertial subrange for frequencies ≳0.1 Hz, but only when the horizontal grid spacing ≲20 m.« less

A Simple Method for Simulating Wind Profiles in the Boundary Layer of Tropical Cyclones

NASA Astrophysics Data System (ADS)

Bryan, George H.; Worsnop, Rochelle P.; Lundquist, Julie K.; Zhang, Jun A.

2017-03-01

A method to simulate characteristics of wind speed in the boundary layer of tropical cyclones in an idealized manner is developed and evaluated. The method can be used in a single-column modelling set-up with a planetary boundary-layer parametrization, or within large-eddy simulations (LES). The key step is to include terms in the horizontal velocity equations representing advection and centrifugal acceleration in tropical cyclones that occurs on scales larger than the domain size. Compared to other recently developed methods, which require two input parameters (a reference wind speed, and radius from the centre of a tropical cyclone) this new method also requires a third input parameter: the radial gradient of reference wind speed. With the new method, simulated wind profiles are similar to composite profiles from dropsonde observations; in contrast, a classic Ekman-type method tends to overpredict inflow-layer depth and magnitude, and two recently developed methods for tropical cyclone environments tend to overpredict near-surface wind speed. When used in LES, the new technique produces vertical profiles of total turbulent stress and estimated eddy viscosity that are similar to values determined from low-level aircraft flights in tropical cyclones. Temporal spectra from LES produce an inertial subrange for frequencies ≳ 0.1 Hz, but only when the horizontal grid spacing ≲ 20 m.

Wave energetics of the southern hemisphere of Mars

NASA Astrophysics Data System (ADS)

Battalio, Michael; Szunyogh, Istvan; Lemmon, Mark

2018-07-01

An assessment of the energetics of transient waves in the southern hemisphere of Mars is presented using the Mars Analysis Correction Data Assimilation (MACDA) dataset (v1.0) and the eddy kinetic energy equation. The dataset is divided into four representative periods covering the summer and winter solstices, a late fall period, and an early spring period for three Mars years. Spring eddies are the most intense, with eddies during the fall being less intense due to a marginally more stable mean-temperature profile and reduced recirculation of ageostrophic geopotential fluxes compared to the spring. Eddy kinetic energy during winter is reduced in intensity as a result of the winter solstitial pause in wave activity, and eddy kinetic energy during the summer is limited. Baroclinic energy conversion acts as a source in fall and spring but disappears during the winter as a result of a stabilized vertical temperature profile. Barotropic energy conversion acts as both a source and a sink of eddy kinetic energy, being most positive during the solstitial pause. Eddies take a northwest to southeast track across the southern highlands in the fall but have a more zonal track in the spring due to stronger eddy kinetic energy advection. Wave energetics is less intense in the southern compared to the northern hemisphere as a result of a shallower baroclinically unstable vertical profile.

Complementary Use of Glider Data, Altimetry, and Model for Exploring Mesoscale Eddies in the Tropical Pacific Solomon Sea

NASA Astrophysics Data System (ADS)

Gourdeau, L.; Verron, J.; Chaigneau, A.; Cravatte, S.; Kessler, W.

2017-11-01

Mesoscale activity is an important component of the Solomon Sea circulation that interacts with the energetic low-latitude western boundary currents of the South Tropical Pacific Ocean carrying waters of subtropical origin before joining the equatorial Pacific. Mixing associated with mesoscale activity could explain water mass transformation observed in the Solomon Sea that likely impacts El Niño Southern Oscillation dynamics. This study makes synergetic use of glider data, altimetry, and high-resolution model for exploring mesoscale eddies, especially their vertical structures, and their role on the Solomon Sea circulation. The description of individual eddies observed by altimetry and gliders provides the first elements to characterize the 3-D structure of these tropical eddies, and confirms the usefulness of the model to access a more universal view of such eddies. Mesoscale eddies appear to have a vertical extension limited to the Surface Waters (SW) and the Upper Thermocline Water (UTW), i.e., the first 140-150 m depth. Most of the eddies are nonlinear, meaning that eddies can trap and transport water properties. But they weakly interact with the deep New Guinea Coastal Undercurrent that is a key piece of the equatorial circulation. Anticyclonic eddies are particularly efficient to advect salty and warm SW coming from the intrusion of equatorial Pacific waters at Solomon Strait, and to impact the characteristics of the New Guinea Coastal Current. Cyclonic eddies are particularly efficient to transport South Pacific Tropical Water (SPTW) anomalies from the North Vanuatu Jet and to erode by diapycnal mixing the high SPTW salinity.

Contribution of mesoscale eddies to Black Sea ventilation

NASA Astrophysics Data System (ADS)

Capet, Arthur; Mason, Evan; Pascual, Ananda; Grégoire, Marilaure

2017-04-01

The shoaling of the Black Sea oxycline is one of the most urgent environmental issues in the Black Sea. The permanent oxycline derives directly from the Black Sea permanent stratification and has shoaled alarmingly in the last decades, due to a shifting balance between oxygen consumption and ventilation processes (Capet et al. 2016). The understanding of this balance is thus of the utmost importance and requires to quantify 1) the export of nutrients and organic materials from the shelf regions to the open sea and 2) the ventilation processes. These two processes being influenced by mesoscale features, it is critical to understand the role of the semi-permanent mesoscale structures in horizontal (center/periphery) and vertical (diapycnal and isopycnal) exchanges. A useful insight can be obtained by merging observations from satellite altimeter and in situ profilers (ARGO). In such composite analyses, eddies are first automatically identified and tracked from altimeter data (Mason et al. 2014, py-eddy-tracker). Vertical ARGO profiles are then expressed in terms of their position relative to eddy centers and radii. Derived statistics indicate how consistently mesoscale eddies alter the vertical structure, and provide a deeper understanding of the associated horizontal and vertical fluxes. However, this data-based approach is limited in the Black Sea due to the lower quality of gridded altimetric products in the vicinity of the coast, where semi-permanent mesoscale structures prevail. To complement the difficult analysis of this sparse dataset, a compositing methodology. is also applied to model outputs from the 5km GHER-BHAMBI Black Sea implementation (CMEMS BS-MFC). Characteristic biogeochemical anomalies associated with eddies in the model are analyzed per se, and compared to the observation-based analysis. Capet, A., Stanev, E. V., Beckers, J.-M., Murray, J. W., and Grégoire, M.: Decline of the Black Sea oxygen inventory, Biogeosciences, 13, 1287-1297, doi:10.5194/bg-13-1287-2016, 2016. Mason, Evan, Ananda Pascual, and James C. McWilliams. "A new sea surface height-based code for oceanic mesoscale eddy tracking." Journal of Atmospheric and Oceanic Technology 31.5 (2014): 1181-1188.

Benchmarking the mesoscale variability in global ocean eddy-permitting numerical systems

NASA Astrophysics Data System (ADS)

Cipollone, Andrea; Masina, Simona; Storto, Andrea; Iovino, Doroteaciro

2017-10-01

The role of data assimilation procedures on representing ocean mesoscale variability is assessed by applying eddy statistics to a state-of-the-art global ocean reanalysis (C-GLORS), a free global ocean simulation (performed with the NEMO system) and an observation-based dataset (ARMOR3D) used as an independent benchmark. Numerical results are computed on a 1/4 ∘ horizontal grid (ORCA025) and share the same resolution with ARMOR3D dataset. This "eddy-permitting" resolution is sufficient to allow ocean eddies to form. Further to assessing the eddy statistics from three different datasets, a global three-dimensional eddy detection system is implemented in order to bypass the need of regional-dependent definition of thresholds, typical of commonly adopted eddy detection algorithms. It thus provides full three-dimensional eddy statistics segmenting vertical profiles from local rotational velocities. This criterion is crucial for discerning real eddies from transient surface noise that inevitably affects any two-dimensional algorithm. Data assimilation enhances and corrects mesoscale variability on a wide range of features that cannot be well reproduced otherwise. The free simulation fairly reproduces eddies emerging from western boundary currents and deep baroclinic instabilities, while underestimates shallower vortexes that populate the full basin. The ocean reanalysis recovers most of the missing turbulence, shown by satellite products , that is not generated by the model itself and consistently projects surface variability deep into the water column. The comparison with the statistically reconstructed vertical profiles from ARMOR3D show that ocean data assimilation is able to embed variability into the model dynamics, constraining eddies with in situ and altimetry observation and generating them consistently with local environment.

The Influence of Slope Breaks on Lava Flow Surface Disruption

NASA Technical Reports Server (NTRS)

Glaze, Lori S.; Baloga, Stephen M.; Fagents, Sarah A.; Wright, Robert

2014-01-01

Changes in the underlying slope of a lava flow impart a significant fraction of rotational energy beyond the slope break. The eddies, circulation and vortices caused by this rotational energy can disrupt the flow surface, having a significant impact on heat loss and thus the distance the flow can travel. A basic mechanics model is used to compute the rotational energy caused by a slope change. The gain in rotational energy is deposited into an eddy of radius R whose energy is dissipated as it travels downstream. A model of eddy friction with the ambient lava is used to compute the time-rate of energy dissipation. The key parameter of the dissipation rate is shown to be rho R(sup 2/)mu, where ? is the lava density and mu is the viscosity, which can vary by orders of magnitude for different flows. The potential spatial disruption of the lava flow surface is investigated by introducing steady-state models for the main flow beyond the steepening slope break. One model applies to slow-moving flows with both gravity and pressure as the driving forces. The other model applies to fast-moving, low-viscosity, turbulent flows. These models provide the flow velocity that establishes the downstream transport distance of disrupting eddies before they dissipate. The potential influence of slope breaks is discussed in connection with field studies of lava flows from the 1801 Hualalai and 1823 Keaiwa Kilauea, Hawaii, and 2004 Etna eruptions.

Comparative Study Of Four Models Of Turbulence

NASA Technical Reports Server (NTRS)

Menter, Florian R.

1996-01-01

Report presents comparative study of four popular eddy-viscosity models of turbulence. Computations reported for three different adverse pressure-gradient flowfields. Detailed comparison of numerical results and experimental data given. Following models tested: Baldwin-Lomax, Johnson-King, Baldwin-Barth, and Wilcox.

Enhanced vertical mixing within mesoscale eddies due to high frequency winds in the South China Sea

NASA Astrophysics Data System (ADS)

Cardona, Yuley; Bracco, Annalisa

The South China Sea is a marginal basin with a complex circulation influenced by the East Asian Monsoon, river discharge and intricate bathymetry. As a result, both the mesoscale eddy field and the near-inertial energy distribution display large spatial variability and they strongly influence the oceanic transport and mixing. With an ensemble of numerical integrations using a regional ocean model, this work investigates how the temporal resolution of the atmospheric forcing fields modifies the horizontal and vertical velocity patterns and impacts the transport properties in the basin. The response of the mesoscale circulation in the South China Sea is investigated under three different forcing conditions: monthly, daily and 6-hourly momentum and heat fluxes. While the horizontal circulation does not display significant differences, the representation of the vertical velocity field displays high sensitivity to the frequency of the wind forcing. If the wind field contains energy at the inertial frequency or higher (daily and 6-hourly cases), then submesoscale fronts, vortex Rossby waves and near inertial waves are excited as ageostrophic expression of the vigorous eddy field. Those quasi- and near-inertial waves dominate the vertical velocity field in the mixed layer (vortex Rossby waves) and below the first hundred meters (near inertial waves) and they are responsible for the differences in the vertical transport properties under the various forcing fields as quantified by frequency spectra, vertical velocity profiles and vertical dispersion of Lagrangian tracers.

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.

Acoustic measurements of the spatial and temporal structure of the near-bottom boundary layer in the 1990-1991 STRESS experiment

NASA Astrophysics Data System (ADS)

Lynch, James F.; Irish, James D.; Gross, Thomas F.; Wiberg, Patricia L.; Newhall, Arthur E.; Traykovski, Peter A.; Warren, Joseph D.

1997-08-01

As part of the 1990-1991 Sediment TRansport Events on Shelves and Slopes (STRESS) experiment, a 5 MHz Acoustic BackScatter System (ABSS) was deployed in 90 m of water to measure vertical profiles of near-bottom suspended sediment concentration. By looking at the vertical profile of concentration from 0 to 50 cm above bottom (cmab) with 1 cm vertical resolution, the ABSS was able to examine the detailed structure of the bottom boundary layer created by combined wave and current stresses. The acoustic profiles clearly showed the wave-current boundary layer, which extends to (order) 10 cmab. The profiles also showed evidence of an "intermediate" boundary layer, also influenced by combined wave and current stresses, just above the wave-current boundary layer. This paper examines the boundary-layer structure by comparing acoustic data obtained by the authors to a 1-D eddy viscosity model formulation. Specifically, these data are compared to a simple extension of the Grant-Glenn-Madsen model formulation. Also of interest is the appearance of apparently 3-D "advective plume" structures in these data. This is an interesting feature in a site which was initially chosen to be a good example of (temporally averaged) 1-D bottom boundary-layer dynamics. Computer modeling and sector-scanning sonar images are presented to justify the plausibility of observing 3-D structure at the STRESS site. 1997 Elsevier Science Ltd

Statistical characterization of high-to-medium frequency mesoscale gravity waves by lidar-measured vertical winds and temperatures in the MLT

NASA Astrophysics Data System (ADS)

Lu, Xian; Chu, Xinzhao; Li, Haoyu; Chen, Cao; Smith, John A.; Vadas, Sharon L.

2017-09-01

We present the first statistical study of gravity waves with periods of 0.3-2.5 h that are persistent and dominant in the vertical winds measured with the University of Colorado STAR Na Doppler lidar in Boulder, CO (40.1°N, 105.2°W). The probability density functions of the wave amplitudes in temperature and vertical wind, ratios of these two amplitudes, phase differences between them, and vertical wavelengths are derived directly from the observations. The intrinsic period and horizontal wavelength of each wave are inferred from its vertical wavelength, amplitude ratio, and a designated eddy viscosity by applying the gravity wave polarization and dispersion relations. The amplitude ratios are positively correlated with the ground-based periods with a coefficient of 0.76. The phase differences between the vertical winds and temperatures (φW -φT) follow a Gaussian distribution with 84.2±26.7°, which has a much larger standard deviation than that predicted for non-dissipative waves ( 3.3°). The deviations of the observed phase differences from their predicted values for non-dissipative waves may indicate wave dissipation. The shorter-vertical-wavelength waves tend to have larger phase difference deviations, implying that the dissipative effects are more significant for shorter waves. The majority of these waves have the vertical wavelengths ranging from 5 to 40 km with a mean and standard deviation of 18.6 and 7.2 km, respectively. For waves with similar periods, multiple peaks in the vertical wavelengths are identified frequently and the ones peaking in the vertical wind are statistically longer than those peaking in the temperature. The horizontal wavelengths range mostly from 50 to 500 km with a mean and median of 180 and 125 km, respectively. Therefore, these waves are mesoscale waves with high-to-medium frequencies. Since they have recently become resolvable in high-resolution general circulation models (GCMs), this statistical study provides an important and timely reference for them.

Study of Near-Surface Models in Large-Eddy Simulations of a Neutrally Stratified Atmospheric Boundary Layer

NASA Technical Reports Server (NTRS)

Senocak, I.; Ackerman, A. S.; Kirkpatrick, M. P.; Stevens, D. E.; Mansour, N. N.

2004-01-01

Large-eddy simulation (LES) is a widely used technique in armospheric modeling research. In LES, large, unsteady, three dimensional structures are resolved and small structures that are not resolved on the computational grid are modeled. A filtering operation is applied to distinguish between resolved and unresolved scales. We present two near-surface models that have found use in atmospheric modeling. We also suggest a simpler eddy viscosity model that adopts Prandtl's mixing length model (Prandtl 1925) in the vicinity of the surface and blends with the dynamic Smagotinsky model (Germano et al, 1991) away from the surface. We evaluate the performance of these surface models by simulating a neutraly stratified atmospheric boundary layer.

Francis-99 turbine numerical flow simulation of steady state operation using RANS and RANS/LES turbulence model

NASA Astrophysics Data System (ADS)

Minakov, A.; Platonov, D.; Sentyabov, A.; Gavrilov, A.

2017-01-01

We performed numerical simulation of flow in a laboratory model of a Francis hydroturbine at three regimes, using two eddy-viscosity- (EVM) and a Reynolds stress (RSM) RANS models (realizable k-ɛ, k-ω SST, LRR) and detached-eddy-simulations (DES), as well as large-eddy simulations (LES). Comparison of calculation results with the experimental data was carried out. Unlike the linear EVMs, the RSM, DES, and LES reproduced well the mean velocity components, and pressure pulsations in the diffusor draft tube. Despite relatively coarse meshes and insufficient resolution of the near-wall region, LES, DES also reproduced well the intrinsic flow unsteadiness and the dominant flow structures and the associated pressure pulsations in the draft tube.

Nonlinear closures for scale separation in supersonic magnetohydrodynamic turbulence

NASA Astrophysics Data System (ADS)

Grete, Philipp; Vlaykov, Dimitar G.; Schmidt, Wolfram; Schleicher, Dominik R. G.; Federrath, Christoph

2015-02-01

Turbulence in compressible plasma plays a key role in many areas of astrophysics and engineering. The extreme plasma parameters in these environments, e.g. high Reynolds numbers, supersonic and super-Alfvenic flows, however, make direct numerical simulations computationally intractable even for the simplest treatment—magnetohydrodynamics (MHD). To overcome this problem one can use subgrid-scale (SGS) closures—models for the influence of unresolved, subgrid-scales on the resolved ones. In this work we propose and validate a set of constant coefficient closures for the resolved, compressible, ideal MHD equations. The SGS energies are modeled by Smagorinsky-like equilibrium closures. The turbulent stresses and the electromotive force (EMF) are described by expressions that are nonlinear in terms of large scale velocity and magnetic field gradients. To verify the closures we conduct a priori tests over 137 simulation snapshots from two different codes with varying ratios of thermal to magnetic pressure ({{β }p}=0.25,1,2.5,5,25) and sonic Mach numbers ({{M}s}=2,2.5,4). Furthermore, we make a comparison to traditional, phenomenological eddy-viscosity and α -β -γ closures. We find only mediocre performance of the kinetic eddy-viscosity and α -β -γ closures, and that the magnetic eddy-viscosity closure is poorly correlated with the simulation data. Moreover, three of five coefficients of the traditional closures exhibit a significant spread in values. In contrast, our new closures demonstrate consistently high correlations and constant coefficient values over time and over the wide range of parameters tested. Important aspects in compressible MHD turbulence such as the bi-directional energy cascade, turbulent magnetic pressure and proper alignment of the EMF are well described by our new closures.

VALIDITY OF A TWO-DIMENSIONAL MODEL FOR VARIABLE-DENSITY HYDRODYNAMIC CIRCULATION

EPA Science Inventory

A three-dimensional model of temperatures and currents has been formulated to assist in the analysis and interpretation of the dynamics of stratified lakes. In this model, nonlinear eddy coefficients for viscosity and conductivities are included. A two-dimensional model (one vert...

A novel VLES model accounting for near-wall turbulence: physical rationale and applications

NASA Astrophysics Data System (ADS)

Jakirlic, Suad; Chang, Chi-Yao; Kutej, Lukas; Tropea, Cameron

2014-11-01

A novel VLES (Very Large Eddy Simulation) model whose non-resolved residual turbulence is modelled by using an advanced near-wall eddy-viscosity model accounting for the near-wall Reynolds stress anisotropy influence on the turbulence viscosity by modelling appropriately the velocity scale in the relevant formulation (Hanjalic et al., 2004) is proposed. It represents a variable resolution Hybrid LES/RANS (Reynolds-Averaged Navier-Stokes) computational scheme enabling a seamless transition from RANS to LES depending on the ratio of the turbulent viscosities associated with the unresolved scales corresponding to the LES cut-off and the `unsteady' scales pertinent to the turbulent properties of the VLES residual motion, which varies within the flow domain. The VLES method is validated interactively in the process of the model derivation by computing fully-developed flow in a plane channel (important representative of wall-bounded flows, underlying the log-law for the velocity field, for studying near-wall Reynolds stress anisotropy) and a separating flow over a periodic arrangement of smoothly-contoured 2-D hills. The model performances are also assessed in capturing the natural decay of the homogeneous isotropic turbulence. The model is finally applied to swirling flow in a vortex tube, flow in an IC-engine configuration and flow past a realistic car model.

Simulations of Turbulent Flow Over Complex Terrain Using an Immersed-Boundary Method

NASA Astrophysics Data System (ADS)

DeLeon, Rey; Sandusky, Micah; Senocak, Inanc

2018-02-01

We present an immersed-boundary method to simulate high-Reynolds-number turbulent flow over the complex terrain of Askervein and Bolund Hills under neutrally-stratified conditions. We reconstruct both the velocity and the eddy-viscosity fields in the terrain-normal direction to produce turbulent stresses as would be expected from the application of a surface-parametrization scheme based on Monin-Obukhov similarity theory. We find that it is essential to be consistent in the underlying assumptions for the velocity reconstruction and the eddy-viscosity relation to produce good results. To this end, we reconstruct the tangential component of the velocity field using a logarithmic velocity profile and adopt the mixing-length model in the near-surface turbulence model. We use a linear interpolation to reconstruct the normal component of the velocity to enforce the impermeability condition. Our approach works well for both the Askervein and Bolund Hills when the flow is attached to the surface, but shows slight disagreement in regions of flow recirculation, despite capturing the flow reversal.

Simulations of Turbulent Flow Over Complex Terrain Using an Immersed-Boundary Method

NASA Astrophysics Data System (ADS)

DeLeon, Rey; Sandusky, Micah; Senocak, Inanc

2018-06-01

We present an immersed-boundary method to simulate high-Reynolds-number turbulent flow over the complex terrain of Askervein and Bolund Hills under neutrally-stratified conditions. We reconstruct both the velocity and the eddy-viscosity fields in the terrain-normal direction to produce turbulent stresses as would be expected from the application of a surface-parametrization scheme based on Monin-Obukhov similarity theory. We find that it is essential to be consistent in the underlying assumptions for the velocity reconstruction and the eddy-viscosity relation to produce good results. To this end, we reconstruct the tangential component of the velocity field using a logarithmic velocity profile and adopt the mixing-length model in the near-surface turbulence model. We use a linear interpolation to reconstruct the normal component of the velocity to enforce the impermeability condition. Our approach works well for both the Askervein and Bolund Hills when the flow is attached to the surface, but shows slight disagreement in regions of flow recirculation, despite capturing the flow reversal.

A new time scale based k-epsilon model for near wall turbulence

NASA Technical Reports Server (NTRS)

Yang, Z.; Shih, T. H.

1992-01-01

A k-epsilon model is proposed for wall bonded turbulent flows. In this model, the eddy viscosity is characterized by a turbulent velocity scale and a turbulent time scale. The time scale is bounded from below by the Kolmogorov time scale. The dissipation equation is reformulated using this time scale and no singularity exists at the wall. The damping function used in the eddy viscosity is chosen to be a function of R(sub y) = (k(sup 1/2)y)/v instead of y(+). Hence, the model could be used for flows with separation. The model constants used are the same as in the high Reynolds number standard k-epsilon model. Thus, the proposed model will be also suitable for flows far from the wall. Turbulent channel flows at different Reynolds numbers and turbulent boundary layer flows with and without pressure gradient are calculated. Results show that the model predictions are in good agreement with direct numerical simulation and experimental data.

Near-wall k-epsilon turbulence modeling

NASA Technical Reports Server (NTRS)

Mansour, N. N.; Kim, J.; Moin, P.

1987-01-01

The flow fields from a turbulent channel simulation are used to compute the budgets for the turbulent kinetic energy (k) and its dissipation rate (epsilon). Data from boundary layer simulations are used to analyze the dependence of the eddy-viscosity damping-function on the Reynolds number and the distance from the wall. The computed budgets are used to test existing near-wall turbulence models of the k-epsilon type. It was found that the turbulent transport models should be modified in the vicinity of the wall. It was also found that existing models for the different terms in the epsilon-budget are adequate in the region from the wall, but need modification near the wall. The channel flow is computed using a k-epsilon model with an eddy-viscosity damping function from the data and no damping functions in the epsilon-equation. These computations show that the k-profile can be adequately predicted, but to correctly predict the epsilon-profile, damping functions in the epsilon-equation are needed.

New time scale based k-epsilon model for near-wall turbulence

NASA Technical Reports Server (NTRS)

Yang, Z.; Shih, T. H.

1993-01-01

A k-epsilon model is proposed for wall bonded turbulent flows. In this model, the eddy viscosity is characterized by a turbulent velocity scale and a turbulent time scale. The time scale is bounded from below by the Kolmogorov time scale. The dissipation equation is reformulated using this time scale and no singularity exists at the wall. The damping function used in the eddy viscosity is chosen to be a function of R(sub y) = (k(sup 1/2)y)/v instead of y(+). Hence, the model could be used for flows with separation. The model constants used are the same as in the high Reynolds number standard k-epsilon model. Thus, the proposed model will be also suitable for flows far from the wall. Turbulent channel flows at different Reynolds numbers and turbulent boundary layer flows with and without pressure gradient are calculated. Results show that the model predictions are in good agreement with direct numerical simulation and experimental data.

A defect stream function, law of the wall/wake method for compressible turbulent boundary layers

NASA Technical Reports Server (NTRS)

Barnwell, Richard W.; Dejarnette, Fred R.; Wahls, Richard A.

1989-01-01

The application of the defect stream function to the solution of the two-dimensional, compressible boundary layer is examined. A law of the wall/law of the wake formulation for the inner part of the boundary layer is presented which greatly simplifies the computational task near the wall and eliminates the need for an eddy viscosity model in this region. The eddy viscosity model in the outer region is arbitrary. The modified Crocco temperature-velocity relationship is used as a simplification of the differential energy equation. Formulations for both equilibrium and nonequilibrium boundary layers are presented including a constrained zero-order form which significantly reduces the computational workload while retaining the significant physics of the flow. A formulation for primitive variables is also presented. Results are given for the constrained zero-order and second-order equilibrium formulations and are compared with experimental data. A compressible wake function valid near the wall has been developed from the present results.

A normal stress subgrid-scale eddy viscosity model in large eddy simulation

NASA Technical Reports Server (NTRS)

Horiuti, K.; Mansour, N. N.; Kim, John J.

1993-01-01

The Smagorinsky subgrid-scale eddy viscosity model (SGS-EVM) is commonly used in large eddy simulations (LES) to represent the effects of the unresolved scales on the resolved scales. This model is known to be limited because its constant must be optimized in different flows, and it must be modified with a damping function to account for near-wall effects. The recent dynamic model is designed to overcome these limitations but is compositionally intensive as compared to the traditional SGS-EVM. In a recent study using direct numerical simulation data, Horiuti has shown that these drawbacks are due mainly to the use of an improper velocity scale in the SGS-EVM. He also proposed the use of the subgrid-scale normal stress as a new velocity scale that was inspired by a high-order anisotropic representation model. The testing of Horiuti, however, was conducted using DNS data from a low Reynolds number channel flow simulation. It was felt that further testing at higher Reynolds numbers and also using different flows (other than wall-bounded shear flows) were necessary steps needed to establish the validity of the new model. This is the primary motivation of the present study. The objective is to test the new model using DNS databases of high Reynolds number channel and fully developed turbulent mixing layer flows. The use of both channel (wall-bounded) and mixing layer flows is important for the development of accurate LES models because these two flows encompass many characteristic features of complex turbulent flows.

Analysis of turbulent heat transfer, mass transfer, and friction in smooth tubes at high Prandtl and Schmidt numbers

NASA Technical Reports Server (NTRS)

Deissler, Robert G

1955-01-01

The expression for eddy diffusivity from a previous analysis was modified in order to account for the effect of kinematic viscosity on the turbulence in the region close to a wall. By using the modified expression, good agreement was obtained between predicted and experimental results for heat and mass transfer at Prandtl and Schmidt numbers between 0.5 and 3000. The effects of length-to-diameter ratio and of variable viscosity were also investigated for a wide range of Prandtl numbers.

Nonlinear Eddy-Eddy Interactions in Dry Atmospheres Macroturbulence

NASA Astrophysics Data System (ADS)

Ait Chaalal, F.; Schneider, T.

2012-12-01

The statistical moment equations derived from the atmospheric equation of motions are not closed. However neglecting the large-scale eddy-eddy nonlinear interactions in an idealized dry general circulation model (GCM), which is equivalent to truncating the moment equations at the second order, can reproduce some of the features of the general circulation ([1]), highlighting the significance of eddy-mean flow interactions and the weakness of eddy-eddy interactions in atmospheric macroturbulence ([2]). The goal of the present study is to provide new insight into the rôle of these eddy-eddy interactions and discuss the relevance of a simple stochastic parametrization to represent them. We investigate in detail the general circulation in an idealized dry GCM, comparing full simulations with simulations where the eddy-eddy interactions are removed. The radiative processes are parametrized through Newtonian relaxation toward a radiative-equilibrium state with a prescribed equator to pole temperature contrast. A convection scheme relaxing toward a prescribed convective vertical lapse rate mimics some aspects of moist convection. The study is performed over a wide range of parameters covering the planetary rotation rate, the equator to pole temperature contrast and the vertical lapse rate. Particular attention is given to the wave-mean flow interactions and to the spectral budget. It is found that the no eddy-eddy simulations perform well when the baroclinic activity is weaker, for example for lower equator to pole temperature contrasts or higher rotation rates: the mean meridional circulation is well reproduced, with realistic eddy-driven jets and energy-containing eddy length scales of the order of the Rossby deformation radius. For a stronger baroclinic activity the no eddy-eddy model does not achieve a realistic isotropization of the eddies, the meridional circulation is compressed in the meridional direction and secondary eddy-driven jets emerge. In addition, the baroclinic wave activity does not reach the upper troposphere in association with a very weak or absent Rossby wave absorption in the upper subtropical troposphere. Understanding these deficiencies and the rôle of the eddy-eddy nonlinear interactions in determining the mean meridional circulation paves the way to the development of stochastic third order moments parametrizations, to eventually build GCMs that directly solve for the flow statistics and that could provide a deeper understanding of anthropogenic and natural climate changes. [1] O'Gorman, P. A., & Schneider, T. 2007, Geophysical Research Letters, 34, 22801 [2] Schneider, T., and C. C. Walker, 2006, Journal of the Atmospheric Sciences, 63, 1569-1586.

Characterizing frontal eddies along the East Australian Current from HF radar observations

NASA Astrophysics Data System (ADS)

Schaeffer, Amandine; Gramoulle, A.; Roughan, M.; Mantovanelli, A.

2017-05-01

The East Australian Current (EAC) dominates the ocean circulation along south-eastern Australia, however, little is known about the submesoscale frontal instabilities associated with this western boundary current. One year of surface current measurements from HF radars, in conjunction with mooring and satellite observations, highlight the occurrence and propagation of meanders and frontal eddies along the inshore edge of the EAC. Eddies were systematically identified using the geometry of the high spatial resolution (˜1.5 km) surface currents, and tracked every hour. Cyclonic eddies were observed irregularly, on average every 7 days, with inshore radius ˜10 km. Among various forms of structures, frontal eddies associated with EAC meanders were characterized by poleward advection speeds of ˜0.3-0.4 m/s, migrating as far as 500 km south, based on satellite imagery. Flow field kinematics show that cyclonic eddies have high Rossby numbers (0.6-1.9) and enhance particle dispersion. Patches of intensified surface divergence at the leading edge of the structures are expected to generate vertical uplift. This is confirmed by subsurface measurements showing temperature uplift of up to 55 m over 24 h and rough estimates of vertical velocities of 10s of meters per day. While frontal eddies propagate through the radar domain independently of local wind stress, upfront wind can influence their stalling and growth, and can also generate large cold core eddies through intense shear. Such coherent structures are a major mechanism for the transport and entrainment of nutrient rich coastal or deep waters, influencing physical and biological dynamics, and connectivity over large distances.

A singularity free approach to post glacial rebound calculations

NASA Technical Reports Server (NTRS)

Fang, Ming; Hager, Bradford H.

1994-01-01

Calculating the post glacial response of a viscoelastic Earth model using the exponential decay normal mode technique leads to intrinsic singularities if viscosity varies continuously as a function of radius. We develop a numerical implementation of the Complex Real Fourier transform (CRFT) method as an accurate and stable procedure to avoid these singularities. Using CRFT, we investigate the response of a set of Maxwell Earth models to surface loading. We find that the effect of expanding a layered viscosity structure into a continuously varying structure is to destroy the modes associated with the boundary between layers. Horizontal motion is more sensitive than vertical motion to the viscosity structure just below the lithosphere. Horizontal motion is less sensitive to the viscosity of the lower mantle than the vertical motion is. When the viscosity increases at 670 km depth by a factor of about 60, the response of the lower mantle is close to its elastic limit. Any further increase of the viscosity contrast at 670 km depth or further increase of viscosity as a continuous function of depth starting from 670 km depth is unlikely to be resolved.

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 .

Inexpensive Eddy-Current Standard

NASA Technical Reports Server (NTRS)

Berry, Robert F., Jr.

1985-01-01

Radial crack replicas serve as evaluation standards. Technique entails intimately joining two pieces of appropriate aluminum alloy stock and centering drilled hole through and along interface. Bore surface of hole presents two vertical stock interface lines 180 degrees apart. These lines serve as radial crack defect replicas during eddy-current technique setup and verification.

Assessment of benthic flux of dissolved organic carbon in wetland and estuarine sediments using the eddy-correlation technique

NASA Astrophysics Data System (ADS)

Swett, M. P.; Amirbahman, A.; Boss, E.

2009-12-01

Wetland and estuarine sediments release significant amounts of dissolved organic carbon (DOC) due to high levels of microbial activity, particularly sulfate reduction. Changes in climate and hydrologic conditions have a potential to alter DOC release from these systems as well. This is a concern, as high levels of DOC can lead to mobilization of toxic metals and organics in natural waters. In addition, source waters high in DOC produce undesirable disinfection byproducts in water treatment. Various in situ methods, such as peepers and sediment core centrifugation, exist to quantify vertical benthic fluxes of DOC and other dissolved species from the sediment-water interface (SWI). These techniques, however, are intrusive and involve disturbance of the sediment environment. Eddy-correlation allows for real-time, non-intrusive, in situ flux measurement of important analytes, such as O2 and DOC. An Acoustic Doppler Velocimeter (ADV) is used to obtain three-dimensional fluid velocity measurements. The eddy-correlation technique employs the mathematical separation of fluid velocity into mean velocity and fluctuating velocity components, with the latter representing turbulent eddy velocity. DOC concentrations are measured using a colored dissolved organic matter (CDOM) fluorometer, and instantaneous vertical flux is determined from the correlated data. This study assesses DOC flux at three project sites: a beaver pond in the Lower Penobscot Watershed, Maine; a mudflat in Penobscot River, Maine; and a mudflat in Great Bay, New Hampshire. Eddy flux values are compared with results obtained using peepers and centrifugation, as well as vertical profiling.

Energetics of eddy-mean flow interactions in the Brazil current between 20°S and 36°S

NASA Astrophysics Data System (ADS)

Magalhães, F. C.; Azevedo, J. L. L.; Oliveira, L. R.

2017-08-01

The energetics of eddy-mean flow interactions in the Brazil Current (BC) between 20°S and 36°S are investigated in 19 transects perpendicular to the 200 m isobath. Ten years (2000-2009) of output data from the Hybrid Coordinate Ocean Model (HYCOM) NCODA reanalysis, with a spatial resolution of 1/12.5° and 5 day averages, are used. The mean kinetic energy (MKE) and eddy kinetic energy (EKE) fields presented the same subsurface spatial pattern but with reduced values. The EKE increases southward, with high values along the BC path and the offshore portion of the jet. The values of the barotropic conversion term (BTC) are highest in the surface layers and decreased with depth, whereas the values of the baroclinic conversion term (BCC) and the vertical eddy heat flux (VEHF) are highest in the subsurface. Despite the vertical thickening of the BC, the highest energy conversion rates are confined to the upper 700 m of the water column. The energetic analysis showed that the current features mixed instability processes. The vertical weighted mean of the BTC and BCC presented an oscillatory pattern related to the bathymetry. The eddy field accelerates the time-mean flow upstream and downstream of bathymetric features and drains energy from the time-mean flow over the features. The BC is baroclinically unstable south of 28°S, and the highest energy conversion rates occur in Cabo de São Tomé, Cabo Frio, and the Cone do Rio Grande.

Simulation of Atmospheric-Entry Capsules in the Subsonic Regime

NASA Technical Reports Server (NTRS)

Murman, Scott M.; Childs, Robert E.; Garcia, Joseph A.

2015-01-01

The accuracy of Computational Fluid Dynamics predictions of subsonic capsule aerodynamics is examined by comparison against recent NASA wind-tunnel data at high-Reynolds-number flight conditions. Several aspects of numerical and physical modeling are considered, including inviscid numerical scheme, mesh adaptation, rough-wall modeling, rotation and curvature corrections for eddy-viscosity models, and Detached-Eddy Simulations of the unsteady wake. All of these are considered in isolation against relevant data where possible. The results indicate that an improved predictive capability is developed by considering physics-based approaches and validating the results against flight-relevant experimental data.

Subgrid or Reynolds stress-modeling for three-dimensional turbulence computations

NASA Technical Reports Server (NTRS)

Rubesin, M. W.

1975-01-01

A review is given of recent advances in two distinct computational methods for evaluating turbulence fields, namely, statistical Reynolds stress modeling and turbulence simulation, where large eddies are followed in time. It is shown that evaluation of the mean Reynolds stresses, rather than use of a scalar eddy viscosity, permits an explanation of streamline curvature effects found in several experiments. Turbulence simulation, with a new volume averaging technique and third-order accurate finite-difference computing is shown to predict the decay of isotropic turbulence in incompressible flow with rather modest computer storage requirements, even at Reynolds numbers of aerodynamic interest.

Effect of mesoscale oceanic eddies on mid-latitude storm-tracks.

NASA Astrophysics Data System (ADS)

Foussard, Alexis; Lapeyre, Guillaume; Plougonven, Riwal

2017-04-01

Sharp sea surface temperature (SST) gradients associated with oceanic western boundary currents (WBC) exert an influence on the position and intensity of mid-latitude storm-tracks. This occurs through strong surface baroclinicity maintained by cross frontal SST gradient and deep vertical atmospheric motion due to convection on the warm flank of the WBC. However the additional role of mesoscale oceanic structures (30-300km) has not yet been explored although they have a non-negligible influence on surface heat fluxes. Using the Weather Research and Forecasting model, we investigate the potential role of these oceanic eddies in the case of an idealized atmospheric mid-latitude storm track forced by a mesoscale oceanic eddy field superposed with a large-scale SST gradient. Surface latent and sensible fluxes are shown to react with a non-linear response to the SST variations, providing additional heat and moisture supply at large scales. The atmospheric response is not restricted to the boundary layer but reaches the free troposphere, especially through increased water vapor vertical transport and latent heat release. This additional heating in presence of eddies is balanced by a shift of the storm-track and its poleward heat flux toward high latitudes, with amplitude depending on atmospheric configuration and eddies amplitude. We also explore how this displacement of perturbations changes the position and structure of the mid-latitude jet through eddy momentum fluxes.

On Eddy Viscosity, Energy Cascades, and the Horizontal Resolution of Gridded Satellite Altimeter Products

DTIC Science & Technology

2013-02-01

Manuscript received 21 December 2011, in final form 23 August 2012) ABSTRACT Motivated by the recent interest in ocean energetics, the widespread use...Inhomogeneous two-dimensional turbu- lence in the atmosphere. Advances in Turbulence, G. Comte - Bellot and J. Mathieu, Eds., Springer-Verlag, 269-278

Numerical simulation of turbulence and sediment transport of medium sand

NASA Astrophysics Data System (ADS)

Schmeeckle, M. W.

2012-12-01

Eleven numerical simulations, ranging from no transport to bedload to vigorous suspension transport, are presented of a combined large eddy simulation (LES) and distinct element model (DEM) of an initially flat bed of medium sand. The fluid and particles are fully coupled in momentum. The friction coefficient, defined here as the squared ratio of the friction velocity to the depth-averaged velocity, is in good agreement with well-known rough bed relations at no transport and increases with the intensity of bedload transport. The friction coefficient nearly doubles in value at the onset of sediment suspension owing to a rapid increase of the depth over which particles and fluid exchange momentum. The friction coefficient decreases with increasing suspension intensity because of increasingly stable stratification. Fluid Reynolds stress and time-averaged velocity profiles in the bedload regime agree well with previous experiments and simulations. Also consistent with previous studies of suspended sediment, there is an increase in slope of the lower portion of the velocity profile that has been modeled in the past using stably stratified eddy viscosity closures or an adjusted von Karman constant. Stokes numbers in the simulations, using an estimated lagrangian integral time scale, are less than unity. As such, particles faithfully follow the fluid, except for particle settling and grain-grain interactions near the bed. Fluid-particle velocity correlation coefficients approach one in portions of the flow where volumetric sediment concentrations are below about ten percent. Bedload entrainment is critically connected to vertical velocity fluctuations. When a fluid packet approaches the bed from the interior of the flow (i.e. a sweep), fluid is forced into the bed, and at the edges of the sweep, fluid is forced out of the bed. Much of the particle entrainment occurs at these sweep edges. Fluid velocity statistics following the particles reveal that moving bedload particles are preferentially concentrated in zones of upward fluid velocity. This may explain previous observations noting a rapid vertical rise at the beginning of saltation trajectories. The simulations described here have no lift forces. Because of the short particle time scales relative to that of the turbulent structures, high transport stage bedload entrainment zones involve mutual interaction between turbulence structures and bed deformation. These deformation structures appear as depressed areas of the bed at the center of the sweep and raised areas of entraining particles at the edges of the sweep penetration. Suspended sediment entrainment structures are similar to these bedload entrainment structures but have much larger scales. Preferential concentration of suspended grains in zones of upward moving fluid dampens turbulence intensities and momentum transport. Much of the suspended transport takes place within this highly concentrated near-bed zone of damped turbulence. Particle-fluid correlation coefficients are relatively low in the lower portion of this highly concentrated suspended sediment zone, owing to particle-particle interactions. As such, Rouse-like profiles utilizing eddy viscosity closures, adjusted according to flux Richardson numbers, do not adequately describe the physics of this zone.

The effect of artificial bulk viscosity in simulations of forced compressible turbulence

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

Campos, A.; Morgan, B.

The use of an artificial bulk viscosity for shock stabilization is a common approach employed in turbulence simulations with high-order numerics. The effect of the artificial bulk viscosity is analyzed in the context of large eddy simulations by using as a test case simulations of linearly-forced compressible homogeneous turbulence (Petersen and Livescu, 2010 [12]). This case is unique in that it allows for the specification of a priori target values for total dissipation and ratio of solenoidal to dilatational dissipation. A comparison between these target values and the true predicted levels of dissipation is thus used to investigate the performancemore » of the artificial bulk viscosity. Results show that the artificial bulk viscosity is effective at achieving stable solutions, but also leads to large values of artificial dissipation that outweigh the physical dissipation caused by fluid viscosity. An alternate approach, which employs the artificial thermal conductivity only, shows that the dissipation of dilatational modes is entirely due to the fluid viscosity. However, this method leads to unwanted Gibbs oscillations around the shocklets. The use of shock sensors that further localize the artificial bulk viscosity did not reduce the amount of artificial dissipation introduced by the artificial bulk viscosity. Finally, an improved forcing function that explicitly accounts for the role of the artificial bulk viscosity in the budget of turbulent kinetic energy was explored.« less

The effect of artificial bulk viscosity in simulations of forced compressible turbulence

DOE PAGES

Campos, A.; Morgan, B.

2018-05-17

The use of an artificial bulk viscosity for shock stabilization is a common approach employed in turbulence simulations with high-order numerics. The effect of the artificial bulk viscosity is analyzed in the context of large eddy simulations by using as a test case simulations of linearly-forced compressible homogeneous turbulence (Petersen and Livescu, 2010 [12]). This case is unique in that it allows for the specification of a priori target values for total dissipation and ratio of solenoidal to dilatational dissipation. A comparison between these target values and the true predicted levels of dissipation is thus used to investigate the performancemore » of the artificial bulk viscosity. Results show that the artificial bulk viscosity is effective at achieving stable solutions, but also leads to large values of artificial dissipation that outweigh the physical dissipation caused by fluid viscosity. An alternate approach, which employs the artificial thermal conductivity only, shows that the dissipation of dilatational modes is entirely due to the fluid viscosity. However, this method leads to unwanted Gibbs oscillations around the shocklets. The use of shock sensors that further localize the artificial bulk viscosity did not reduce the amount of artificial dissipation introduced by the artificial bulk viscosity. Finally, an improved forcing function that explicitly accounts for the role of the artificial bulk viscosity in the budget of turbulent kinetic energy was explored.« less

Scale-Similar Models for Large-Eddy Simulations

NASA Technical Reports Server (NTRS)

Sarghini, F.

1999-01-01

Scale-similar models employ multiple filtering operations to identify the smallest resolved scales, which have been shown to be the most active in the interaction with the unresolved subgrid scales. They do not assume that the principal axes of the strain-rate tensor are aligned with those of the subgrid-scale stress (SGS) tensor, and allow the explicit calculation of the SGS energy. They can provide backscatter in a numerically stable and physically realistic manner, and predict SGS stresses in regions that are well correlated with the locations where large Reynolds stress occurs. In this paper, eddy viscosity and mixed models, which include an eddy-viscosity part as well as a scale-similar contribution, are applied to the simulation of two flows, a high Reynolds number plane channel flow, and a three-dimensional, nonequilibrium flow. The results show that simulations without models or with the Smagorinsky model are unable to predict nonequilibrium effects. Dynamic models provide an improvement of the results: the adjustment of the coefficient results in more accurate prediction of the perturbation from equilibrium. The Lagrangian-ensemble approach [Meneveau et al., J. Fluid Mech. 319, 353 (1996)] is found to be very beneficial. Models that included a scale-similar term and a dissipative one, as well as the Lagrangian ensemble averaging, gave results in the best agreement with the direct simulation and experimental data.

Numerical study of large-eddy breakup and its effect on the drag characteristics of boundary layers

NASA Technical Reports Server (NTRS)

Kinney, R. B.; Taslim, M. E.; Hung, S. C.

1985-01-01

The break-up of a field of eddies by a flat-plate obstacle embedded in a boundary layer is studied using numerical solutions to the two-dimensional Navier-Stokes equations. The flow is taken to be incompressible and unsteady. The flow field is initiated from rest. A train of eddies of predetermined size and strength are swept into the computational domain upstream of the plate. The undisturbed velocity profile is given by the Blasius solution. The disturbance vorticity generated at the plate and wall, plus that introduced with the eddies, mix with the background vorticity and is transported throughout the entire flow. All quantities are scaled by the plate length, the unidsturbed free-stream velocity, and the fluid kinematic viscosity. The Reynolds number is 1000, the Blasius boundary layer thickness is 2.0, and the plate is positioned a distance of 1.0 above the wall. The computational domain is four units high and sixteen units long.

Noncontact Measurement Of Sizes And Eccentricities Of Holes

NASA Technical Reports Server (NTRS)

Chern, Engmin J.

1993-01-01

Semiautomatic eddy-current-probe apparatus makes noncontact measurements of nominally round holes in electrically conductive specimens and processes measurement data into diameters and eccentricities of holes. Includes x-y translation platform, which holds specimen and moves it horizontally. Probe mounted on probe scanner, positioning probe along vertical (z) direction and rotates probe about vertical axis at preset low speed. Eddy-current sensing coil mounted in side of probe near tip. As probe rotates, impedance analyzer measures electrical impedance (Z) of coil as function of instantaneous rotation angle. Translation and rotation mechanisms and impedance analyzer controlled by computer, which also processes impedance-measurement data.

Wind Effects on Flow Patterns and Net Fluxes in Density-Driven High-Latitude Channel Flow

NASA Astrophysics Data System (ADS)

Huntley, Helga S.; Ryan, Patricia

2018-01-01

A semianalytic two-dimensional model is used to analyze the interplay between the different forces acting on density-driven flow in high-latitude channels. In particular, the balance between wind stress, viscous forces, baroclinicity, and sea surface slope adjustments under specified flux conditions is examined. Weak winds are found not to change flow patterns appreciably, with minimal (<7%) adjustments to horizontal velocity maxima. In low-viscosity regimes, strong winds change the flow significantly, especially at the surface, by either strengthening the dual-jet pattern, established without wind, by a factor of 2-3 or initiating return flow at the surface. A nonzero flux does not result in the addition of a uniform velocity throughout the channel cross section, but modifies both along-channel and cross-channel velocities to become more symmetric, dominated by a down-channel jet centered in the domain and counter-clockwise lateral flow. We also consider formulations of the model that allow adjustments of the net flux in response to the wind. Flow patterns change, beyond uniform intensification or weakening, only for strong winds and high Ekman number. Comparisons of the model results to observational data collected in Nares Strait in the Canadian Archipelago in the summer of 2007 show rough agreement, but the model misses the upstream surface jet on the east side of the strait and propagates bathymetric effects too strongly in the vertical for this moderately high eddy viscosity. Nonetheless, the broad strokes of the observed high-latitude flow are reproduced.

Fronts and eddies: Engines for biogeochemical variability of the Central Red Sea during winter-spring periods

NASA Astrophysics Data System (ADS)

Zarokanellos, Nikolaos; Jones, Burton

2017-04-01

The central Red Sea (CRS) has been shown to be characterized by significant eddy activity throughout the year. In winter, weakened stratification may lead to enhanced vertical exchange contributing to physical and biogeochemical processes. In winter 2014-2015 we began an extended glider time series to monitor a region in the northern CRS where eddy activity is significant. Remote sensing and glider observations that include CTD, oxygen, CDOM and chlorophyll fluorescence, and multi-wavelength optical backscatter, have been used to characterize the effects of winter mixing and eddy activity in this region. During winter, deep mixing driven by surface cooling and strong winds combined with eddy features, can supply nutrients into the upper layer dramatically modifies the environment from its typically stratified conditions. These mixing events disperse the phytoplankton from the deep chlorophyll maximum throughout the upper mixed layer, and increase the chlorophyll signature detected by ocean color imagery. In addition to the mixing, cyclonic eddies in the region can enhance the vertical displacement of deeper, nutrient containing water toward the euphotic zone contributing to increased chlorophyll concentration and biological productivity. Remote sensing analyses indicate that these eddies also contribute to significant horizontal dispersion including the exchange between the open sea and coastal coral reef ecosystems. During the winter mixing periods, diel fluctuations in phytoplankton biomass have been observed indicative of solar driven plankton dynamics. The biogeochemical response to the subsurface physical processes provides a sensitive indicator to the processes that result from the mixing and eddy dynamics - processes that are not necessarily detectable via remote sensing. In order to understand the seasonal responses, but also the interannual influences on these processes, sustained in situ autonomous platform measurements are essential.

Impacts of mesoscale eddies in the South China Sea on biogeochemical cycles

NASA Astrophysics Data System (ADS)

Guo, Mingxian; Chai, Fei; Xiu, Peng; Li, Shiyu; Rao, Shivanesh

2015-09-01

Biogeochemical cycles associated with mesoscale eddies in the South China Sea (SCS) were investigated. The study was based on a coupled physical-biogeochemical Pacific Ocean model (Regional Ocean Model System-Carbon, Silicate, and Nitrogen Ecosystem, ROMS-CoSiNE) simulation for the period from 1991 to 2008. A total of 568 mesoscale eddies with lifetime longer than 30 days were used in the analysis. Composite analysis revealed that the cyclonic eddies were associated with abundance of nutrients, phytoplankton, and zooplankton while the anticyclonic eddies depressed biogeochemical cycles, which are generally controlled by the eddy pumping mechanism. In addition, diatoms were dominant in phytoplankton species due to the abundance of silicate. Dipole structures of vertical fluxes with net upward motion in cyclonic eddies and net downward motion in anticyclonic eddies were revealed. During the lifetime of an eddy, the evolutions of physical, biological, and chemical structures were not linearly coupled at the eddy core where plankton grew, and composition of the community depended not only on the physical and chemical processes but also on the adjustments by the predator-prey relationship.

Seasonal overturning circulation in the Red Sea: 2. Winter circulation

NASA Astrophysics Data System (ADS)

Yao, Fengchao; Hoteit, Ibrahim; Pratt, Larry J.; Bower, Amy S.; Köhl, Armin; Gopalakrishnan, Ganesh; Rivas, David

2014-04-01

The shallow winter overturning circulation in the Red Sea is studied using a 50 year high-resolution MITgcm (MIT general circulation model) simulation with realistic atmospheric forcing. The overturning circulation for a typical year, represented by 1980, and the climatological mean are analyzed using model output to delineate the three-dimensional structure and to investigate the underlying dynamical mechanisms. The horizontal model circulation in the winter of 1980 is dominated by energetic eddies. The climatological model mean results suggest that the surface inflow intensifies in a western boundary current in the southern Red Sea that switches to an eastern boundary current north of 24°N. The overturning is accomplished through a cyclonic recirculation and a cross-basin overturning circulation in the northern Red Sea, with major sinking occurring along a narrow band of width about 20 km along the eastern boundary and weaker upwelling along the western boundary. The northward pressure gradient force, strong vertical mixing, and horizontal mixing near the boundary are the essential dynamical components in the model's winter overturning circulation. The simulated water exchange is not hydraulically controlled in the Strait of Bab el Mandeb; instead, the exchange is limited by bottom and lateral boundary friction and, to a lesser extent, by interfacial friction due to the vertical viscosity at the interface between the inflow and the outflow.

Calculations of High-Temperature Jet Flow Using Hybrid Reynolds-Average Navier-Stokes Formulations

NASA Technical Reports Server (NTRS)

Abdol-Hamid, Khaled S.; Elmiligui, Alaa; Giriamaji, Sharath S.

2008-01-01

Two multiscale-type turbulence models are implemented in the PAB3D solver. The models are based on modifying the Reynolds-averaged Navier Stokes equations. The first scheme is a hybrid Reynolds-averaged- Navier Stokes/large-eddy-simulation model using the two-equation k(epsilon) model with a Reynolds-averaged-Navier Stokes/large-eddy-simulation transition function dependent on grid spacing and the computed turbulence length scale. The second scheme is a modified version of the partially averaged Navier Stokes model in which the unresolved kinetic energy parameter f(sub k) is allowed to vary as a function of grid spacing and the turbulence length scale. This parameter is estimated based on a novel two-stage procedure to efficiently estimate the level of scale resolution possible for a given flow on a given grid for partially averaged Navier Stokes. It has been found that the prescribed scale resolution can play a major role in obtaining accurate flow solutions. The parameter f(sub k) varies between zero and one and is equal to one in the viscous sublayer and when the Reynolds-averaged Navier Stokes turbulent viscosity becomes smaller than the large-eddy-simulation viscosity. The formulation, usage methodology, and validation examples are presented to demonstrate the enhancement of PAB3D's time-accurate turbulence modeling capabilities. The accurate simulations of flow and turbulent quantities will provide a valuable tool for accurate jet noise predictions. Solutions from these models are compared with Reynolds-averaged Navier Stokes results and experimental data for high-temperature jet flows. The current results show promise for the capability of hybrid Reynolds-averaged Navier Stokes and large eddy simulation and partially averaged Navier Stokes in simulating such flow phenomena.

A Coherent vorticity preserving eddy-viscosity correction for Large-Eddy Simulation

NASA Astrophysics Data System (ADS)

Chapelier, J.-B.; Wasistho, B.; Scalo, C.

2018-04-01

This paper introduces a new approach to Large-Eddy Simulation (LES) where subgrid-scale (SGS) dissipation is applied proportionally to the degree of local spectral broadening, hence mitigated or deactivated in regions dominated by large-scale and/or laminar vortical motion. The proposed coherent-vorticity preserving (CvP) LES methodology is based on the evaluation of the ratio of the test-filtered to resolved (or grid-filtered) enstrophy, σ. Values of σ close to 1 indicate low sub-test-filter turbulent activity, justifying local deactivation of the SGS dissipation. The intensity of the SGS dissipation is progressively increased for σ < 1 which corresponds to a small-scale spectral broadening. The SGS dissipation is then fully activated in developed turbulence characterized by σ ≤σeq, where the value σeq is derived assuming a Kolmogorov spectrum. The proposed approach can be applied to any eddy-viscosity model, is algorithmically simple and computationally inexpensive. LES of Taylor-Green vortex breakdown demonstrates that the CvP methodology improves the performance of traditional, non-dynamic dissipative SGS models, capturing the peak of total turbulent kinetic energy dissipation during transition. Similar accuracy is obtained by adopting Germano's dynamic procedure albeit at more than twice the computational overhead. A CvP-LES of a pair of unstable periodic helical vortices is shown to predict accurately the experimentally observed growth rate using coarse resolutions. The ability of the CvP methodology to dynamically sort the coherent, large-scale motion from the smaller, broadband scales during transition is demonstrated via flow visualizations. LES of compressible channel are carried out and show a good match with a reference DNS.

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.

Eddy Correlation Flux Measurement System (ECOR) Handbook

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

Cook, DR

2011-01-31

The eddy correlation (ECOR) flux measurement system provides in situ, half-hour measurements of the surface turbulent fluxes of momentum, sensible heat, latent heat, and carbon dioxide (CO2) (and methane at one Southern Great Plains extended facility (SGP EF) and the North Slope of Alaska Central Facility (NSA CF). The fluxes are obtained with the eddy covariance technique, which involves correlation of the vertical wind component with the horizontal wind component, the air temperature, the water vapor density, and the CO2 concentration.

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.

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.

Wind-Wave Effects on Vertical Mixing in Chesapeake Bay, USA: comparing observations to second-moment closure predictions.

NASA Astrophysics Data System (ADS)

Fisher, A. W.; Sanford, L. P.; Scully, M. E.

2016-12-01

Coherent wave-driven turbulence generated through wave breaking or nonlinear wave-current interactions, e.g. Langmuir turbulence (LT), can significantly enhance the downward transfer of momentum, kinetic energy, and dissolved gases in the oceanic surface layer. There are few observations of these processes in the estuarine or coastal environments, where wind-driven mixing may co-occur with energetic tidal mixing and strong density stratification. This presents a major challenge for evaluating vertical mixing parameterizations used in modeling estuarine and coastal dynamics. We carried out a large, multi-investigator study of wind-driven estuarine dynamics in the middle reaches of Chesapeake Bay, USA, during 2012-2013. The center of the observational array was an instrumented turbulence tower with both atmospheric and marine turbulence sensors as well as rapidly sampled temperature and conductivity sensors. For this paper, we examined the impacts of surface gravity waves on vertical profiles of turbulent mixing and compared our results to second-moment turbulence closure predictions. Wave and turbulence measurements collected from the vertical array of Acoustic Doppler Velocimeters (ADVs) provided direct estimates of the dominant terms in the TKE budget and the surface wave field. Observed dissipation rates, TKE levels, and turbulent length scales are compared to published scaling relations and used in the calculation of second-moment nonequilibrium stability functions. Results indicate that in the surface layer of the estuary, where elevated dissipation is balanced by vertical divergence in TKE flux, existing nonequilibrium stability functions underpredict observed eddy viscosities. The influences of wave breaking and coherent wave-driven turbulence on modeled and observed stability functions will be discussed further in the context of turbulent length scales, TKE and dissipation profiles, and the depth at which the wave-dominated turbulent transport layer transitions to a turbulent log layer. The influences of fetch-limited wind waves, density stratification, and surface buoyancy fluxes will also be discussed.

A comparison of the structure, properties, and water mass composition of quasi-isotropic eddies in western boundary currents in an eddy-resolving ocean model

NASA Astrophysics Data System (ADS)

Rykova, Tatiana; Oke, Peter R.; Griffin, David A.

2017-06-01

Using output from a near-global eddy-resolving ocean model, we analyse the properties and characteristics of quasi-isotropic eddies in five Western Boundary Current (WBC) regions, including the extensions of the Agulhas, East Australian Current (EAC), Brazil-Malvinas Confluence (BMC), Kuroshio and Gulf Stream regions. We assess the model eddies by comparing to satellite and in situ observations, and show that most aspects of the model's representation of eddies are realistic. We find that the mean eddies differ dramatically between these WBC regions - all with some unique and noteworthy characteristics. We find that the vertical displacement of isopycnals of Agulhas eddies is the greatest, averaging 350-450 m at depths of over 800-900 m. EAC (BMC) eddies are the least (most) barotropic, with only 50% (85-90%) of the velocity associated with the barotropic mode. Kuroshio eddies are the most stratified, resulting in small isopycnal displacement, even for strong eddies; and Gulf Stream eddies carry the most heat. Despite their differences, we explicitly show that the source waters for anticyclonic eddies are a mix of the WBC water (from the boundary current itself) and water that originates equatorward of the WBC eddy-field; and cyclonic eddies are a mix of WBC water and water that originates poleward of the WBC eddy-field.

Observed and modeled mesoscale variability near the Gulf Stream and Kuroshio Extension

NASA Astrophysics Data System (ADS)

Schmitz, William J.; Holland, William R.

1986-08-01

Our earliest intercomparisons between western North Atlantic data and eddy-resolving two-layer quasi-geostrophic symmetric-double-gyre steady wind-forced numerical model results focused on the amplitudes and largest horizontal scales in patterns of eddy kinetic energy, primarily abyssal. Here, intercomparisons are extended to recent eight-layer model runs and new data which allow expansion of the investigation to the Kuroshio Extension and throughout much of the water column. Two numerical experiments are shown to have realistic zonal, vertical, and temporal eddy scales in the vicinity of the Kuroshio Extension in one case and the Gulf Stream in the other. Model zonal mean speeds are larger than observed, but vertical shears are in general agreement with the data. A longitudinal displacement between the maximum intensity in surface and abyssal eddy fields as observed for the North Atlantic is not found in the model results. The numerical simulations examined are highly idealized, notably with respect to basin shape, topography, wind-forcing, and of course dissipation. Therefore the zero-order agreement between modeled and observed basic characteristics of mid-latitude jets and their associated eddy fields suggests that such properties are predominantly determined by the physical mechanisms which dominate the models, where the fluctuations are the result of instability processes. The comparatively high vertical resolution of the model is needed to compare with new higher-resolution data as well as for dynamical reasons, although the precise number of layers required either kinematically or dynamically (or numerically) has not been determined; we estimate four to six when no attempt is made to account for bottom- or near-surface-intensified phenomena.

Distribution of living radiolarians and its response on the environments in spring from the section South China Sea

NASA Astrophysics Data System (ADS)

Zhang, L.; Hu, W.; Chen, M.; Zeng, L.; Xiang, R.; Zhou, W.

2013-12-01

The composition and spatial (horizontal and vertical) distribution of living radiolarians in spring was firstly studied in the section (18°N and 113°E) South China Sea. Vertical plankton tows were collected at depth-intervals from 0 to 300 m in spring using a closing-type net with 62 um mesh size. And we distinguished the living specimens by staining with Rose Bengal. It dominated by tropical-subtropical warm species in spring from the studied areas. The abundance of nassellarians was the almost same as that of spumellarians in the upper-surface waters (0-25m). In the below-surface waters (25-50m), nassellarian abundance was the almost twice that of spumellarians. And the abundances generally decreased with depth (more than 50m), but nasselarian abundance reduced more quickly. The results showed that the horizontal and vertical distribution patterns of living radiolarians were closely related to the mesoscale eddies. The horizontal distributions of radiolarian abundance were uneven and pachy, which may be related to the complicated mecoscale eddies during the sampling period. That is, there were comparatively high abundances in the upper-surface waters where had the cold eddies development. But in the cold eddies of Meigong River mouth, radiolarian abundance was low due to the large input of fresh water, suggesting that low salinity had more important influence than the nutrient on the radiolarian development and reproduction. Vertically, the highest abundances occurred at the mixed layer in the cold eddies, and gradually decreased with depth. However, in the warm eddies, the maximum abundances were in the thermocline layers, where had an abundant supply of nutrients for radiolarians. This study showed that Didymocyrtis tetrathalamus tetrathalamus mostly occurred at the mixed layer, which should be closely related to the cold eddies and rich nutrition and be limited by the fresh water. Based on the distribution of Didymocyrtis tetrathalamus tetrathalamus, we concluded that the influence of west Pacific waters was obviously weak on the northwestern Luzon Island during the sampling period. As a tropical surface warm species, Tetrapyle octacantha was also found to be indicator of tropical upwelling eutrophication water. Acanthodesmia vinculata was mainly living in the mixed layer, and had a good response to the cold eddies far away the continental shelf. Besides, we also concluded that Siphonosphaera polysiphonia should be tropical surface warm species, having a gregarious life, which had a closely related to the warm eddies. Interestingly, the typical deep-dwellers (Cornutella profunda and Cyrtopera laguncula) occurred in the different depth intervals, even in the upper-surface waters, which suggested that the temperature might not be the mostly one of factors to control their living-depth. This study was funded by the following research programs: the National Natural Science Foundation of China (Nos. 41276051, 91228207, 40906030).

Glider and satellite high resolution monitoring of a mesoscale eddy in the Algerian basin: effects on the mixed layer depth and biochemistry

NASA Astrophysics Data System (ADS)

Cotroneo, Yuri; Aulicino, Giuseppe; Ruiz, Simón; Pascual, Ananda; Budillon, Giorgio; Fusco, Giannetta; Tintoré, Joaquin

2016-04-01

Despite of the extensive bibliography about the circulation of the Mediterranean Sea and its sub-basins, the debate on mesoscale dynamics and its impacts on biochemical processes is still open because of their intrinsic time scales and of the difficulties in sampling. In order to clarify some of these processes, the "Algerian BAsin Circulation Unmanned Survey - ABACUS" project was proposed and realized through access to JERICO Trans National Access (TNA) infrastructures between September and December 2014. In this framework, a deep glider cruise was carried out in the area between Balearic Islands and Algerian coasts to establish an endurance line for monitoring the basin circulation. During the mission, a mesoscale eddy, identified on satellite altimetry maps, was sampled at high-spatial horizontal resolution (4 km) along its main axes and from surface to 1000 m depth. Data were collected by a Slocum glider equipped with a pumped CTD and biochemical sensors that collected about 100 complete casts inside the eddy. In order to describe the structure of the eddy, in situ data were merged with new generation remotely sensed data as daily synoptic sea surface temperature (SST) and chlorophyll concentration (Chl-a) images from MODIS satellites as well as sea surface height and geostrophic velocities from AVISO. From its origin along the Algerian coast in the eastern part of the basin, the eddy propagated to north-west at a mean speed of about 4 km/day with a mean diameter of 112/130 km, a mean elevation of 15.7 cm and clearly distinguished by the surrounding waters thanks to its higher SST and Chl-a values. Temperature and salinity values along the water column confirm the origin of the eddy from the AC showing the presence of recent Atlantic water in the surface layer and Levantine Intermediate Water (LIW) in the deeper layer. Eddy footprint is clearly evident in the multiparametric vertical sections conducted along its main axes. Deepening of temperature, salinity and density isolines at the center of the eddy is associated with variations in the Chl-a, oxygen concentration and turbidity pattern. In particular at 50 m depth, along the eddy borders, Chl-a values are higher (1.1-5.2 μg/l) than in correspondence of the eddy center (0.5-0.7 μg/l) with maxima values registered in the southeastern sector of the eddy. Calculation of geostrophic velocities along transects and vertical quasi geostrophic velocities (QG-w) over a regular 5 km grid from glider data, helped in describing the mechanism and functioning of the eddy. QG-w presents an asymmetric pattern, with associated relatively strong downwelling in the western part of the eddy and upwelling in the southeastern part of it. This asymmetry in the vertical velocity pattern, bringing LIW in the euphotic layer, as well as eventual advection from the northeastern sector of the eddy may justify the observed increase in Chl-a values.

A Particle Representation Model for the Deformation of Homogeneous Turbulence

NASA Technical Reports Server (NTRS)

Kassinos, S. C.; Reynolds, W. C.

1996-01-01

In simple flows, where the mean deformation rates are mild and the turbulence has time to come to equilibrium with the mean flow, the Reynolds stresses are determined by the applied strain rate. Hence in these flows, it is often adequate to use an eddy-viscosity representation. The modern family of kappa-epsilon models has been very useful in predicting near equilibrium turbulent flows, where the rms deformation rate S is small compared to the reciprocal time scale of the turbulence (epsilon/kappa). In modern engineering applications, turbulence models are quite often required to predict flows with very rapid deformations (large S kappa/epsilon). In these flows, the structure takes some time to respond and eddy viscosity models are inadequate. The response of turbulence to rapid deformations is given by rapid distortion theory (RDT). Under RDT the nonlinear effects due to turbulence-turbulence interactions are neglected in the governing equations, but even when linearized in this fashion, the governing equations are unclosed at the one-point level due to the non-locality of the pressure fluctuations.

Relaxation of an unsteady turbulent boundary layer on a flat plate in an expansion tube

NASA Technical Reports Server (NTRS)

Gurta, R. N.; Trimpi, R. L.

1974-01-01

An analysis is presented for the relaxation of a turbulent boundary layer on a semi-infinite flat plate after passage of a shock wave and a trailing driver gas-driven gas interface. The problem has special application to expansion-tube flows. The flow-governing equations have been transformed into the Crocco variables, and a time-similar solution is presented in terms of the dimensionless distance-time variable alpha and the dimensionless velocity variable beta. An eddy-viscosity model, similar to that of time-steady boundary layers, is applied to the inner and outer regions of the boundary layer. A turbulent Prandtl number equal to the molecular Prandtl number is used to relate the turbulent heat flux to the eddy viscosity. The numerical results, obtained by using the Gauss-Seidel line-relaxation method, indicate that a fully turbulent boundary layer relaxes faster to the final steady-state values of heat transfer and skin friction than a laminar boundary layer. The results also give a fairly good estimate of the local skin friction and heat transfer for near steady-flow conditions.

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

NASA Technical Reports Server (NTRS)

Kim, Sang-Wook; Chen, Yen-Sen

1986-01-01

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

An experimental investigation of turbulent boundary layers along curved surfaces

NASA Technical Reports Server (NTRS)

So, R. M. C.; Mellor, G. L.

1972-01-01

A curved wall tunnel was designed, and an equilibrium turbulent boundary layer was set up on the straight section preceding the curved test section. Turbulent boundary layer flows with uniform and adverse pressure distributions along convex and concave walls were investigated. Hot-wire measurements along the convex surface indicated that turbulent mixing between fluid layers was very much reduced. However, the law of the wall held and the skin friction, thus determined, correlated well with other measurements. Hot-wire measurements along the concave test wall revealed a system of longitudinal vortices inside the boundary layer and confirmed that concave curvature enhances mixing. A self-consistent set of turbulent boundary layer equations for flows along curved surfaces was derived together with a modified eddy viscosity. Solution of these equations together with the modified eddy viscosity gave results that correlated well with the present data on flows along the convex surface with arbitrary pressure distribution. However, it could only be used to predict the mean characteristics of the flow along concave walls because of the existence of the system of longitudinal vortices inside the boundary layer.

A low dimensional dynamical system for the wall layer

NASA Technical Reports Server (NTRS)

Aubry, N.; Keefe, L. R.

1987-01-01

Low dimensional dynamical systems which model a fully developed turbulent wall layer were derived.The model is based on the optimally fast convergent proper orthogonal decomposition, or Karhunen-Loeve expansion. This decomposition provides a set of eigenfunctions which are derived from the autocorrelation tensor at zero time lag. Via Galerkin projection, low dimensional sets of ordinary differential equations in time, for the coefficients of the expansion, were derived from the Navier-Stokes equations. The energy loss to the unresolved modes was modeled by an eddy viscosity representation, analogous to Heisenberg's spectral model. A set of eigenfunctions and eigenvalues were obtained from direct numerical simulation of a plane channel at a Reynolds number of 6600, based on the mean centerline velocity and the channel width flow and compared with previous work done by Herzog. Using the new eigenvalues and eigenfunctions, a new ten dimensional set of ordinary differential equations were derived using five non-zero cross-stream Fourier modes with a periodic length of 377 wall units. The dynamical system was integrated for a range of the eddy viscosity prameter alpha. This work is encouraging.

Mesoscale Eddies Are Oases for Higher Trophic Marine Life

PubMed Central

Godø, Olav R.; Samuelsen, Annette; Macaulay, Gavin J.; Patel, Ruben; Hjøllo, Solfrid Sætre; Horne, John; Kaartvedt, Stein; Johannessen, Johnny A.

2012-01-01

Mesoscale eddies stimulate biological production in the ocean, but knowledge of energy transfers to higher trophic levels within eddies remains fragmented and not quantified. Increasing the knowledge base is constrained by the inability of traditional sampling methods to adequately sample biological processes at the spatio-temporal scales at which they occur. By combining satellite and acoustic observations over spatial scales of 10 s of km horizontally and 100 s of m vertically, supported by hydrographical and biological sampling we show that anticyclonic eddies shape distribution and density of marine life from the surface to bathyal depths. Fish feed along density structures of eddies, demonstrating that eddies catalyze energy transfer across trophic levels. Eddies create attractive pelagic habitats, analogous to oases in the desert, for higher trophic level aquatic organisms through enhanced 3-D motion that accumulates and redistributes biomass, contributing to overall bioproduction in the ocean. Integrating multidisciplinary observation methodologies promoted a new understanding of biophysical interaction in mesoscale eddies. Our findings emphasize the impact of eddies on the patchiness of biomass in the sea and demonstrate that they provide rich feeding habitat for higher trophic marine life. PMID:22272294

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.

The structure of turbulent flow around vertical plates containing holes and attached to a channel bed

NASA Astrophysics Data System (ADS)

Basnet, K.; Constantinescu, G.

2017-11-01

High-resolution, 3-D large eddy simulations are conducted to study the physics of flow past 2-D solid and porous vertical plates of height H mounted on a horizontal surface (no bottom gap) with a fully developed, turbulent incoming flow. The porous plate consists of an array of spanwise-oriented, identical solid cylinders of rectangular cross section. The height of the solid cylinders and the spacing between the solid cylinders, corresponding to the plate's "holes," are kept constant for any given configuration, as the present study considers only plates of uniform porosity. The paper discusses how the mean flow and turbulence structure around the vertical plate, the unsteady forces acting on the plate, the dynamics of the large-scale turbulent eddies, the spectral content of the wake, and the distribution of the bed friction velocity on the horizontal channel bed vary as a function of the plate porosity (0% < P < 36%), the relative spacing between the solid elements of the porous plate (d/H), and the roughness of the channel bed surface. Simulation results are used to explain how the bleeding flow affects the dynamics on the larger billow eddies advected in the separated shear layer (SSL) forming at the top of the plate and the wake structure. It is found that the main recirculation eddy in the wake remains attached to the plate for P < 30%. For larger porosities, the main recirculation eddy forms away from the porous plate. The energy of the billows advected in the SSL decays monotonically with increasing plate porosity. For cases when the recirculation eddy remains attached to the plate, the larger billows advected in the downstream part of the SSL are partially reinjected inside the main recirculation eddy as a result of their interaction with the channel bed. This creates a feedback mechanism that induces large-scale disturbances of the spanwise-oriented vortex tubes advected inside the upstream part of the SSL. Results also show that the mean drag coefficient and the root-mean-square of the drag coefficient fluctuations increase mildly with increasing d/H. Meanwhile, varying d/H has a negligible effect on the position and size of the main recirculation eddy. The presence of large-scale roughness elements (2-D ribs) at the bed results in the decrease of the mean drag coefficient of the plate and, in the case of a solid plate, in a large decrease of the frequency of the large-scale eddies advected in the SSL.

Navier-Stokes Computations With One-Equation Turbulence Model for Flows Along Concave Wall Surfaces

NASA Technical Reports Server (NTRS)

Wang, Chi R.

2005-01-01

This report presents the use of a time-marching three-dimensional compressible Navier-Stokes equation numerical solver with a one-equation turbulence model to simulate the flow fields developed along concave wall surfaces without and with a downstream extension flat wall surface. The 3-D Navier- Stokes numerical solver came from the NASA Glenn-HT code. The one-equation turbulence model was derived from the Spalart and Allmaras model. The computational approach was first calibrated with the computations of the velocity and Reynolds shear stress profiles of a steady flat plate boundary layer flow. The computational approach was then used to simulate developing boundary layer flows along concave wall surfaces without and with a downstream extension wall. The author investigated the computational results of surface friction factors, near surface velocity components, near wall temperatures, and a turbulent shear stress component in terms of turbulence modeling, computational mesh configurations, inlet turbulence level, and time iteration step. The computational results were compared with existing measurements of skin friction factors, velocity components, and shear stresses of the developing boundary layer flows. With a fine computational mesh and a one-equation model, the computational approach could predict accurately the skin friction factors, near surface velocity and temperature, and shear stress within the flows. The computed velocity components and shear stresses also showed the vortices effect on the velocity variations over a concave wall. The computed eddy viscosities at the near wall locations were also compared with the results from a two equation turbulence modeling technique. The inlet turbulence length scale was found to have little effect on the eddy viscosities at locations near the concave wall surface. The eddy viscosities, from the one-equation and two-equation modeling, were comparable at most stream-wise stations. The present one-equation turbulence model is an effective approach for turbulence modeling in the near solid wall surface region of flow over a concave wall.

A frequency domain linearized Navier-Stokes method including acoustic damping by eddy viscosity using RANS

NASA Astrophysics Data System (ADS)

Holmberg, Andreas; Kierkegaard, Axel; Weng, Chenyang

2015-06-01

In this paper, a method for including damping of acoustic energy in regions of strong turbulence is derived for a linearized Navier-Stokes method in the frequency domain. The proposed method is validated and analyzed in 2D only, although the formulation is fully presented in 3D. The result is applied in a study of the linear interaction between the acoustic and the hydrodynamic field in a 2D T-junction, subject to grazing flow at Mach 0.1. Part of the acoustic energy at the upstream edge of the junction is shed as harmonically oscillating disturbances, which are conveyed across the shear layer over the junction, where they interact with the acoustic field. As the acoustic waves travel in regions of strong shear, there is a need to include the interaction between the background turbulence and the acoustic field. For this purpose, the oscillation of the background turbulence Reynold's stress, due to the acoustic field, is modeled using an eddy Newtonian model assumption. The time averaged flow is first solved for using RANS along with a k-ε turbulence model. The spatially varying turbulent eddy viscosity is then added to the spatially invariant kinematic viscosity in the acoustic set of equations. The response of the 2D T-junction to an incident acoustic field is analyzed via a plane wave scattering matrix model, and the result is compared to experimental data for a T-junction of rectangular ducts. A strong improvement in the agreement between calculation and experimental data is found when the modification proposed in this paper is implemented. Discrepancies remaining are likely due to inaccuracies in the selected turbulence model, which is known to produce large errors e.g. for flows with significant rotation, which the grazing flow across the T-junction certainly is. A natural next step is therefore to test the proposed methodology together with more sophisticated turbulence models.

A Large-eddy Simulation Study of Vertical Axis Wind Turbine Wakes in the Atmospheric Boundary Layer

NASA Astrophysics Data System (ADS)

Shamsoddin, Sina; Porté-Agel, Fernando

2016-04-01

Vertical axis wind turbines (VAWTs) offer some advantages over their horizontal axis counterparts, and are being considered as a viable alternative to conventional horizontal axis wind turbines (HAWTs). Nevertheless, a relative shortage of scientific, academic and technical investigations of VAWTs is observed in the wind energy community with respect to HAWTs. Having this in mind, in this work, we aim to study the wake of a single VAWT, placed in the atmospheric boundary layer, using large-eddy simulation (LES) coupled with actuator line model (ALM). It is noteworthy that this is the first time that such a study is being performed. To do this, for a typical 1 MW VAWT design, first, the variation of power coefficient with both the chord length of the blades and the tip-speed ratio is analyzed using LES-ALM, and an optimum combination of chord length and tip-speed ratio is obtained. Subsequently, the wake of a VAWT with these optimum specifications is thoroughly examined by showing different relevant mean and turbulent wake flow statistics. Keywords: vertical axis wind turbine (VAWT); VAWT wake; Atmospheric Boundary Layer (ABL); large eddy simulation (LES); actuator line model (ALM); turbulence.

The Energetics of Transient Eddies in the Martian Northern Hemisphere

NASA Astrophysics Data System (ADS)

Battalio, Joseph Michael; Szunyogh, Istvan; Lemmon, Mark T.

2016-10-01

The energetics of northern hemisphere transient waves in the Mars Analysis Correction Data Assimilation is analyzed. Three periods between the fall and spring equinoxes (Ls=200°-230°, 255°-285°, and 330°-360°) during three Mars Years are selected to exemplify the fall, winter, and spring wave activity. Fall and spring eddy energetics is similar with some inter-annual and inter-seasonal variability, but winter eddy kinetic energy and its transport are strongly reduced in intensity as a result of the solsticial pause in eddy activity. Barotropic energy conversion acts as a sink of eddy kinetic energy throughout the northern hemisphere eddy period with little reduction in amplitude during the solsticial pause. Baroclinic energy conversion acts as a source in fall and spring but disappears during the winter period as a result of the stabilized vertical shear profile of the westerly jet around winter solstice.

Role of eddy pumping in enhancing primary production in the ocean

NASA Technical Reports Server (NTRS)

Falkowski, Paul G.; Kolber, Zbigniew; Ziemann, David; Bienfang, Paul K.

1991-01-01

Eddy pumping is considered to explain the disparity between geochemical estimates and biological measurements of exported production. Episodic nutrient injections from the ocean into the photic zone can be generated by eddy pumping, which biological measurements cannot sample accurately. The enhancement of production is studied with respect to a cyclonic eddy in the subtropical Pacific. A pump-and-probe fluorimeter generates continuous vertical profiles of primary productivity from which the contributions of photochemical and nonphotochemical processes to fluorescence are derived. A significant correlation is observed between the fluorescence measurements and radiocarbon measurements. The results indicate that eddy pumping has an important effect on phytoplankton production and that this production is near the maximum relative specific growth rates. Based on the production enhancement observed in this case, eddy pumping increases total primary production by only 20 percent and does not account for all enhancement.

Heating, moisture, and water budgets of tropical and midlatitude squall lines - Comparisons and sensitivity to longwave radiation

NASA Technical Reports Server (NTRS)

Tao, W.-K.; Simpson, J.; Sui, C.-H.; Ferrier, B.; Lang, S.; Scala, J.; Chou, M.-D.; Pickering, K.

1993-01-01

A 2D time-dependent and nonhydrostatic numerical cloud model is presently used to estimate the heating, moisture, and water budgets in the convective and stratiform regions for both a tropical and a midlatitude squall line. The model encompasses a parameterized, three-class ice phase microphysical scheme and longwave radiative transfer process. It is noted that the convective region plays an important role in the generation of stratiform rainfall for both cases. While a midlevel minimum in the moisture profile for the tropical case is due to vertical eddy transport in the convective region, the contribution to the heating budget by the cloud-scale fluxes is minor; by contrast, the vertical eddy heat-flux is relatively important for the midlatitude case due to the stronger vertical velocities present in the convective cells.

Counter-rotating accretion discs

NASA Astrophysics Data System (ADS)

Dyda, S.; Lovelace, R. V. E.; Ustyugova, G. V.; Romanova, M. M.; Koldoba, A. V.

2015-01-01

Counter-rotating discs can arise from the accretion of a counter-rotating gas cloud on to the surface of an existing corotating disc or from the counter-rotating gas moving radially inwards to the outer edge of an existing disc. At the interface, the two components mix to produce gas or plasma with zero net angular momentum which tends to free-fall towards the disc centre. We discuss high-resolution axisymmetric hydrodynamic simulations of viscous counter-rotating discs for the cases where the two components are vertically separated and radially separated. The viscosity is described by an isotropic α-viscosity including all terms in the viscous stress tensor. For the vertically separated components, a shear layer forms between them and the middle part of this layer free-falls to the disc centre. The accretion rates are increased by factors of ˜102-104 over that for a conventional disc rotating in one direction with the same viscosity. The vertical width of the shear layer and the accretion rate are strongly dependent on the viscosity and the mass fraction of the counter-rotating gas. In the case of radially separated components where the inner disc corotates and the outer disc rotates in the opposite direction, a gap between the two components opens and closes quasi-periodically. The accretion rates are ≳25 times larger than those for a disc rotating in one direction with the same viscosity.

Turbulent fluxes by "Conditional Eddy Sampling"

NASA Astrophysics Data System (ADS)

Siebicke, Lukas

2015-04-01

Turbulent flux measurements are key to understanding ecosystem scale energy and matter exchange, including atmospheric trace gases. While the eddy covariance approach has evolved as an invaluable tool to quantify fluxes of e.g. CO2 and H2O continuously, it is limited to very few atmospheric constituents for which sufficiently fast analyzers exist. High instrument cost, lack of field-readiness or high power consumption (e.g. many recent laser-based systems requiring strong vacuum) further impair application to other tracers. Alternative micrometeorological approaches such as conditional sampling might overcome major limitations. Although the idea of eddy accumulation has already been proposed by Desjardin in 1972 (Desjardin, 1977), at the time it could not be realized for trace gases. Major simplifications by Businger and Oncley (1990) lead to it's widespread application as 'Relaxed Eddy Accumulation' (REA). However, those simplifications (flux gradient similarity with constant flow rate sampling irrespective of vertical wind velocity and introduction of a deadband around zero vertical wind velocity) have degraded eddy accumulation to an indirect method, introducing issues of scalar similarity and often lack of suitable scalar flux proxies. Here we present a real implementation of a true eddy accumulation system according to the original concept. Key to our approach, which we call 'Conditional Eddy Sampling' (CES), is the mathematical formulation of conditional sampling in it's true form of a direct eddy flux measurement paired with a performant real implementation. Dedicated hardware controlled by near-real-time software allows full signal recovery at 10 or 20 Hz, very fast valve switching, instant vertical wind velocity proportional flow rate control, virtually no deadband and adaptive power management. Demonstrated system performance often exceeds requirements for flux measurements by orders of magnitude. The system's exceptionally low power consumption is ideal for the field (one to two orders of magnitude lower compared to current closed-path laser based eddy covariance systems). Potential applications include fluxes of CO2, CH4, N2O, VOCs and other tracers. Finally we assess the flux accuracy of the Conditional Eddy Sampling (CES) approach as in our real implementation relative to alternative techniques including eddy covariance (EC) and relaxed eddy accumulation (REA). We further quantify various sources of instrument and method specific measurement errors. This comparison uses real measurements of 20 Hz turbulent time series of 3D wind velocity, sonic temperature and CO2 mixing ratio over a mixed decidious forest at the 'ICOS' flux tower site 'Hainich', Germany. Results from a simulation using real wind and CO2 timeseries from the Hainich site from 30 April to 3 November 2014 and real instrument performance suggest that the maximum flux estimates error (50% and 75% error quantiles) from Conditional Eddy Sampling (CES) relative to the true flux is 1.3% and 10%, respectively for monthly net fluxes, 1.6% and 7%, respectively for daily net fluxes and 8% and 35%, respectively for 30-minute CO2 flux estimates. Those results from CES are promising and outperform our REA estimates by about a factor of 50 assuming REA with constant b value. Results include flux time series from the EC, CES and REA approaches from 30-min to annual resolution.

Large eddy simulation of incompressible turbulent channel flow

NASA Technical Reports Server (NTRS)

Moin, P.; Reynolds, W. C.; Ferziger, J. H.

1978-01-01

The three-dimensional, time-dependent primitive equations of motion were numerically integrated for the case of turbulent channel flow. A partially implicit numerical method was developed. An important feature of this scheme is that the equation of continuity is solved directly. The residual field motions were simulated through an eddy viscosity model, while the large-scale field was obtained directly from the solution of the governing equations. An important portion of the initial velocity field was obtained from the solution of the linearized Navier-Stokes equations. The pseudospectral method was used for numerical differentiation in the horizontal directions, and second-order finite-difference schemes were used in the direction normal to the walls. The large eddy simulation technique is capable of reproducing some of the important features of wall-bounded turbulent flows. The resolvable portions of the root-mean square wall pressure fluctuations, pressure velocity-gradient correlations, and velocity pressure-gradient correlations are documented.

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.

LES FOR SIMULATING THE GAS EXCHANGE PROCESS IN A SPARK IGNITION ENGINE

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

Ameen, Muhsin M; yang, xiaofeng; kuo, tang-wei

2015-01-01

The gas exchange process is known to be a significant source of cyclic variability in Internal Combustion Engines (ICE). Traditionally, Large Eddy Simulations (LES) are expected to capture these cycle-to-cycle variations. This paper reports a numerical effort to establish best practices for capturing cyclic variability with LES tools in a Transparent Combustion Chamber (TCC) spark ignition engine. The main intention is to examine the sensitivity of cycle averaged mean and Root Mean Square (RMS) flow fields and Proper Orthogonal Decomposition (POD) modes to different computational hardware, adaptive mesh refinement (AMR) and LES sub-grid scale (SGS) models, since these aspects havemore » received little attention in the past couple of decades. This study also examines the effect of near-wall resolution on the predicted wall shear stresses. LES is pursued with commercially available CONVERGE code. Two different SGS models are tested, a one-equation eddy viscosity model and dynamic structure model. The results seem to indicate that both mean and RMS fields without any SGS model are not much different than those with LES models, either one-equation eddy viscosity or dynamic structure model. Computational hardware results in subtle quantitative differences, especially in RMS distributions. The influence of AMR on both mean and RMS fields is negligible. The predicted shear stresses near the liner walls is also found to be relatively insensitive to near-wall resolution except in the valve curtain region.« less

Quantifying mesoscale eddies in the Lofoten Basin

NASA Astrophysics Data System (ADS)

Raj, R. P.; Johannessen, J. A.; Eldevik, T.; Nilsen, J. E. Ø.; Halo, I.

2016-07-01

The Lofoten Basin is the most eddy rich region in the Norwegian Sea. In this paper, the characteristics of these eddies are investigated from a comprehensive database of nearly two decades of satellite altimeter data (1995-2013) together with Argo profiling floats and surface drifter data. An automated method identified 1695/1666 individual anticyclonic/cyclonic eddies in the Lofoten Basin from more than 10,000 altimeter-based eddy observations. The eddies are found to be predominantly generated and residing locally. The spatial distributions of lifetime, occurrence, generation sites, size, intensity, and drift of the eddies are studied in detail. The anticyclonic eddies in the Lofoten Basin are the most long-lived eddies (>60 days), especially in the western part of the basin. We reveal two hotspots of eddy occurrence on either side of the Lofoten Basin. Furthermore, we infer a cyclonic drift of eddies in the western Lofoten Basin. Barotropic energy conversion rates reveals energy transfer from the slope current to the eddies during winter. An automated colocation of surface drifters trapped inside the altimeter-based eddies are used to corroborate the orbital speed of the anticyclonic and cyclonic eddies. Moreover, the vertical structure of the altimeter-based eddies is examined using colocated Argo profiling float profiles. Combination of altimetry, Argo floats, and surface drifter data is therefore considered to be a promising observation-based approach for further studies of the role of eddies in transport of heat and biomass from the slope current to the Lofoten Basin.

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

A High Frequency Response Relaxed Eddy Accumulation Flux Measurement System for Sampling Short-Lived Biogenic Volatile Organic Compounds

EPA Science Inventory

A second-generation relaxed eddy accumulation system was built and tested with the capability to measure vertical biogenic volatile organic compound (VOC) fluxes at levels as low as 10 µg C m⁻² hr⁻¹. The system features a continuous, integrated gas-phase ozo...

Southern Ocean eddy compensation in a forced eddy-resolving GCM

NASA Astrophysics Data System (ADS)

Bruun Poulsen, Mads; Jochum, Markus; Eden, Carsten; Nuterman, Roman

2017-04-01

Contemporary eddy-resolving model studies have demonstrated that the common parameterisation of isopycnal mixing in the ocean is subject to limitations in the Southern Ocean where the mesoscale eddies are of leading order importance to the dynamics. We here present forced simulations from the Community Earth System Model on a global {1/10}° and 1° horizontal grid, the latter employing an eddy parameterisation, where the strength of the zonal wind stress south of 25°S has been varied. With a 50% zonally symmetric increase of the wind stress, we show that the two models arrive at two radically different solutions in terms of the large-scale circulation, with an increase of the deep inflow of water to the Southern Ocean at 40°S by 50% in the high resolution model against 20% at coarse resolution. Together with a weaker vertical displacement of the pycnocline in the 1° model, these results suggest that the parameterised eddies have an overly strong compensating effect on the water mass transformation compared to the explicit eddies. Implications for eddy mixing parameterisations will be discussed.

Numerical solutions of the complete Navier-Strokes equations. no. 27

NASA Technical Reports Server (NTRS)

Hassan, H. A.

1996-01-01

This report describes the development of an enstrophy model capable of predicting turbulence separation and its application to two airfoils at various angles of attack and Mach numbers. In addition, a two equation kappa-xi model with a tensor eddy viscosity was developed. Plans call for this model to be used in calculating three dimensional turbulent flows.

Modeling the dispersal of Levantine Intermediate Water and its role in Mediterranean deep water formation

NASA Astrophysics Data System (ADS)

Wu, Peili; Haines, Keith

1996-03-01

This paper demonstrates the importance of Levantine Intermediate Water (LIW) in the deep water formation process in the Mediterranean using the modular ocean general circulation model at 0.25° resolution, 19 vertical levels, over the entire Mediterranean with an open Gibraltar strait. LIW formation is strongly prescribed in the Rhodes Gyre region by Haney [1971] relaxation, while in other regions, surface salinity relaxation is much reduced by applying the `mixed' thermohaline surface boundary conditions. Isopycnal diagnostics are used to trace water mass movements, and volume fluxes are monitored at straits. Low viscosity and diffusion are used to permit baroclinic eddies to play a role in water mass dispersal. The overall water budget is measured by an average flux at Gibraltar of 0.8 Sv, of which 0.7 Sv is exchanged with the eastern basin at Sicily. LIW (density around 28.95) spreads rapidly after formation throughout the entire Levantine due to baroclinic eddies. Toward the west, LIW accumulates in the northern and central Ionian, with some entering the Adriatic through Otranto and some mixing southward in eddies and exiting to the western Mediterranean through Sicily. LIW is converted to deep water in the south Adriatic at an average rate of 0.4 Sv. Water exchange through the Otranto strait appears to be buoyancy driven, with a strong bias to the end of winter (March-April), while at Sicily the exchange has a strong symmetric seasonal cycle, with maximum transport of 1.1 Sv in December indicating the effects of wind driving. LIW pathways in the west are complex and variable. In the Tyrrhenian, intermediate water becomes uniform on isopycnal surfaces due to eddy stirring. West of Sardinia, two LIW boundary currents are formed in the Balearic basin; one flows northward up the west coast of Sardinia and Corsica, and one westward along the northern African coast. The northward current is consistent with observations, while the westward current is intermittent for the first 10 years, often breaking up into eddies which enter the basin interior. Some observations of high-salinity waters near the African coast may support this interpretation. LIW retains a subsurface salinity maximum of 38.4-38.5 practical salinity units (psu) when reaching the northwestern Mediterranean, contrasting with surface waters fresher than 38.0 psu. West Mediterranean deep water is formed below 1500 m depth with climatological characteristics, when it is mixed and cooled during winter convection in Lions Gyre.

The Drainage of Thin, Vertical, Model Polyurethane Liquid Films

NASA Astrophysics Data System (ADS)

Snow, Steven; Pernisz, Udo; Braun, Richard; Naire, Shailesh

1999-11-01

We have successfully measured the drainage rate of thin, vertically-aligned, liquid films prepared from model polyurethane foam formulations. The pattern of interference fringes in these films was consistent with a wedge-shaped film profile. The time evolution of this wedge shape (the ``collapsing wedge") obeyed a power law relationship between fringe density s and time t of s = k t^m. Experimentally, m ranged from -0.47 to -0.92. The lower bound for m represented a case where the surface viscosity of the film was very high (a ``rigid" surface). Theoretical modeling of this case yielded m = -0.5, in excellent agreement with experiment. Instantaneous film drainage rate (dV/dt) could be extracted from the ``Collapsing Wedge" model. As expected, dV/dt scaled inversely with bulk viscosity. As surfactant concentration was varied at constant bulk viscosity, dV/dt passed through a maximum value, consistent with a model where the rigidity of the surface was a function of both the intensity of surface tension gradients and the surface viscosity of the film. The influence of surface viscosity on dV/dt was also modeled theoretically.

An interior penalty stabilised incompressible discontinuous Galerkin-Fourier solver for implicit large eddy simulations

NASA Astrophysics Data System (ADS)

Ferrer, Esteban

2017-11-01

We present an implicit Large Eddy Simulation (iLES) h / p high order (≥2) unstructured Discontinuous Galerkin-Fourier solver with sliding meshes. The solver extends the laminar version of Ferrer and Willden, 2012 [34], to enable the simulation of turbulent flows at moderately high Reynolds numbers in the incompressible regime. This solver allows accurate flow solutions of the laminar and turbulent 3D incompressible Navier-Stokes equations on moving and static regions coupled through a high order sliding interface. The spatial discretisation is provided by the Symmetric Interior Penalty Discontinuous Galerkin (IP-DG) method in the x-y plane coupled with a purely spectral method that uses Fourier series and allows efficient computation of spanwise periodic three-dimensional flows. Since high order methods (e.g. discontinuous Galerkin and Fourier) are unable to provide enough numerical dissipation to enable under-resolved high Reynolds computations (i.e. as necessary in the iLES approach), we adapt the laminar version of the solver to increase (controllably) the dissipation and enhance the stability in under-resolved simulations. The novel stabilisation relies on increasing the penalty parameter included in the DG interior penalty (IP) formulation. The latter penalty term is included when discretising the linear viscous terms in the incompressible Navier-Stokes equations. These viscous penalty fluxes substitute the stabilising effect of non-linear fluxes, which has been the main trend in implicit LES discontinuous Galerkin approaches. The IP-DG penalty term provides energy dissipation, which is controlled by the numerical jumps at element interfaces (e.g. large in under-resolved regions) such as to stabilise under-resolved high Reynolds number flows. This dissipative term has minimal impact in well resolved regions and its implicit treatment does not restrict the use of large time steps, thus providing an efficient stabilization mechanism for iLES. The IP-DG stabilisation is complemented with a Spectral Vanishing Viscosity (SVV) method, in the z-direction, to enhance stability in the continuous Fourier space. The coupling between the numerical viscosity in the DG plane and the SVV damping, provides an efficient approach to stabilise high order methods at moderately high Reynolds numbers. We validate the formulation for three turbulent flow cases: a circular cylinder at Re = 3900, a static and pitch oscillating NACA 0012 airfoil at Re = 10000 and finally a rotating vertical-axis turbine at Re = 40000, with Reynolds based on the circular diameter, airfoil chord and turbine diameter, respectively. All our results compare favourably with published direct numerical simulations, large eddy simulations or experimental data. We conclude that the DG-Fourier high order solver, with IP-SVV stabilisation, proves to be a valuable tool to predict turbulent flows and associated statistics for both static and rotating machinery.

On the role of the transient eddies in maintaining the seasonal mean circulation

NASA Technical Reports Server (NTRS)

White, G. H.; Hoskins, B. J.

1984-01-01

The role of transient eddies in maintaining the observed local seasonal mean atmospheric circulation was investigated by examining the time-averaged momentum balances and omega equation, using seasonal statistics calculated from daily operational analyses by the European Centre for Medium Range Weather Forecasts. While both the Northern and Southern Hemispheres and several seasons were studied, emphasis was placed upon the Northern Hemisphere during December 1981-February 1982. The results showed that transient eddies played a secondary role in the seasonal mean zonal momentum budget and in the forcing of seasonal mean vertical and a geostrophic motion.

Effect of Turbulence Modeling on Hovering Rotor Flows

NASA Technical Reports Server (NTRS)

Yoon, Seokkwan; Chaderjian, Neal M.; Pulliam, Thomas H.; Holst, Terry L.

2015-01-01

The effect of turbulence models in the off-body grids on the accuracy of solutions for rotor flows in hover has been investigated. Results from the Reynolds-Averaged Navier-Stokes and Laminar Off-Body models are compared. Advection of turbulent eddy viscosity has been studied to find the mechanism leading to inaccurate solutions. A coaxial rotor result is also included.

A computer program for the calculation of laminar and turbulent boundary layer flows

NASA Technical Reports Server (NTRS)

Dwyer, H. A.; Doss, E. D.; Goldman, A. L.

1972-01-01

The results are presented of a study to produce a computer program to calculate laminar and turbulent boundary layer flows. The program is capable of calculating the following types of flow: (1) incompressible or compressible, (2) two dimensional or axisymmetric, and (3) flows with significant transverse curvature. Also, the program can handle a large variety of boundary conditions, such as blowing or suction, arbitrary temperature distributions and arbitrary wall heat fluxes. The program has been specialized to the calculation of equilibrium air flows and all of the thermodynamic and transport properties used are for air. For the turbulent transport properties, the eddy viscosity approach has been used. Although the eddy viscosity models are semi-empirical, the model employed in the program has corrections for pressure gradients, suction and blowing and compressibility. The basic method of approach is to put the equations of motion into a finite difference form and then solve them by use of a digital computer. The program is written in FORTRAN 4 and requires small amounts of computer time on most scientific machines. For example, most laminar flows can be calculated in less than one minute of machine time, while turbulent flows usually require three or four minutes.

A Galilean and tensorial invariant k-epsilon model for near wall turbulence

NASA Technical Reports Server (NTRS)

Yang, Z.; Shih, T. H.

1993-01-01

A k-epsilon model is proposed for wall bounded turbulent flows. In this model, the eddy viscosity is characterized by a turbulent velocity scale and a turbulent time scale. The time scale is bounded from below by the Kolmogorov time scale. The dissipation rate equation is reformulated using this time scale and no singularity exists at the wall. A new parameter R = k/S(nu) is introduced to characterize the damping function in the eddy viscosity. This parameter is determined by local properties of both the mean and the turbulent flow fields and is free from any geometry parameter. The proposed model is then Galilean and tensorial invariant. The model constants used are the same as in the high Reynolds number Standard k-epsilon Model. Thus, the proposed model will also be suitable for flows far from the wall. Turbulent channel flows and turbulent boundary layer flows with and without pressure gradients are calculated. Comparisons with the data from direct numerical simulations and experiments show that the model predictions are excellent for turbulent channel flows and turbulent boundary layers with favorable pressure gradients, good for turbulent boundary layers with zero pressure gradients, and fair for turbulent boundary layer with adverse pressure gradients.

Thermo-fluid-dynamics of turbulent boundary layer over a moving continuous flat sheet in a parallel free stream

NASA Astrophysics Data System (ADS)

Afzal, Bushra; Noor Afzal Team; Bushra Afzal Team

2014-11-01

The momentum and thermal turbulent boundary layers over a continuous moving sheet subjected to a free stream have been analyzed in two layers (inner wall and outer wake) theory at large Reynolds number. The present work is based on open Reynolds equations of momentum and heat transfer without any closure model say, like eddy viscosity or mixing length etc. The matching of inner and outer layers has been carried out by Izakson-Millikan-Kolmogorov hypothesis. The matching for velocity and temperature profiles yields the logarithmic laws and power laws in overlap region of inner and outer layers, along with friction factor and heat transfer laws. The uniformly valid solution for velocity, Reynolds shear stress, temperature and thermal Reynolds heat flux have been proposed by introducing the outer wake functions due to momentum and thermal boundary layers. The comparison with experimental data for velocity profile, temperature profile, skin friction and heat transfer are presented. In outer non-linear layers, the lowest order momentum and thermal boundary layer equations have also been analyses by using eddy viscosity closure model, and results are compared with experimental data. Retired Professor, Embassy Hotel, Rasal Ganj, Aligarh 202001 India.

Impacts of mesoscale eddies on biogeochemical cycles in the South China Sea

NASA Astrophysics Data System (ADS)

Xiu, P.; Chai, F.; Guo, M.

2016-02-01

Biogeochemical cycles associated with mesoscale eddies in the South China Sea (SCS) are investigated by using satellite surface chlorophyll concentration, altimeter data, satellite sea surface temperature, and a coupled physical-biogeochemical Pacific Ocean model (ROMS-CoSiNE) simulation for the period from 1991 to 2007. Considering the annual mean, composite analysis reveals that cyclonic eddies are associated with higher concentrations of nutrients, phytoplankton and zooplankton while the anticyclonic eddies are with lower concentrations compared with surrounding waters, which is generally controlled by the eddy pumping mechanism. Dipole structures of vertical fluxes with net upward motion in cyclonic eddies and net downward motion in anticyclonic eddies are also revealed. During the lifetime of an eddy, the evolutions of physical, biological, and chemical structures are not linearly coupled at the eddy core where plankton grow and composition of the community depend not only on the physical and chemical processes but also on the adjustments by the predator-prey relationship. Considering the seasonal variability, we find eddy pumping mechanisms are generally dominant in winter and eddy advection effects are dominant in summer. Over the space, variability of chlorophyll to the west of Luzon Strait and off northwest of Luzon Island are mainly controlled by eddy pumping mechanism. In regions off the Vietnam coast, chlorophyll distributions are generally associated with horizontal eddy advection. This research highlights different mesoscale mechanisms affecting biological structures that can potentially disturb ocean biogeochemical cycling processes in the South China Sea.

South-Eastern Bay of Biscay eddy-induced anomalies and their effect on chlorophyll distribution

NASA Astrophysics Data System (ADS)

Caballero, Ainhoa; Rubio, Anna; Ruiz, Simón; Le Cann, Bernard; Testor, Pierre; Mader, Julien; Hernández, Carlos

2016-10-01

The analysis of deep-water glider hydrographic and fluorescence data, together with satellite measurements provides a new insight into eddy-induced anomalies within the South-Eastern Bay of Biscay, during summer. Two cyclonic eddies and a SWODDY have been observed in different glider transects and by means of different sources of satellite data. Vertical profiles reveal complex structures (characteristic of the second baroclinic mode): upward/downward displacement of the seasonal/permanent thermocline in the case of X13 and the opposite thermocline displacements in the case of the cyclones. This is a typical behaviour of mode-water and "cyclonic thinny" eddies. A qualitative analysis of the vertical velocities in the anticyclone indicates that though geostrophy dominates the main water column, depressing the isopycnals, near the sea-surface the eddy-wind interaction affects the vertical currents, favouring Ekman pumping and upwelling. The analysis of the Θ-S properties corroborates that inside cyclones and between the 26 and 27 isopycnals, net downwelling occurs. These two types of intra-thermocline lenses appear to deeply impact the Chl-a fluorescence profiles, since the maximum Chl-a fluorescence is located just below the seasonal thermocline. The mean Chl-a fluorescence was higher in the anticyclone than within the cyclones and the mean for the entire study period; the highest values were observed in the centre of the anticyclone. These results are in agreement with previous findings concerning the SWODDY F90 and surrounding cyclones, located in the South-Western Bay of Biscay. Significant differences in the Θ-S properties of the two cyclonic mesoscale structures have been observed: higher temperatures and lower salinity in the easternmost cyclone. Finally, time variation of the salinity content of the shallowest water masses of the anticyclone (salinity decreasing over time), probably indicates advective mixing processes occurred during the mission.

Measurement of eddy-current distribution in the vacuum vessel of the Sino-UNIted Spherical Tokamak.

PubMed

Li, G; Tan, Y; Liu, Y Q

2015-08-01

Eddy currents have an important effect on tokamak plasma equilibrium and control of magneto hydrodynamic activity. The vacuum vessel of the Sino-UNIted Spherical Tokamak is separated into two hemispherical sections by a toroidal insulating barrier. Consequently, the characteristics of eddy currents are more complex than those found in a standard tokamak. Thus, it is necessary to measure and analyze the eddy-current distribution. In this study, we propose an experimental method for measuring the eddy-current distribution in a vacuum vessel. By placing a flexible printed circuit board with magnetic probes onto the external surface of the vacuum vessel to measure the magnetic field parallel to the surface and then subtracting the magnetic field generated by the vertical-field coils, the magnetic field due to the eddy current can be obtained, and its distribution can be determined. We successfully applied this method to the Sino-UNIted Spherical Tokamak, and thus, we obtained the eddy-current distribution despite the presence of the magnetic field generated by the external coils.

Global Ocean Vertical Velocity From a Dynamically Consistent Ocean State Estimate

NASA Astrophysics Data System (ADS)

Liang, Xinfeng; Spall, Michael; Wunsch, Carl

2017-10-01

Estimates of the global ocean vertical velocities (Eulerian, eddy-induced, and residual) from a dynamically consistent and data-constrained ocean state estimate are presented and analyzed. Conventional patterns of vertical velocity, Ekman pumping, appear in the upper ocean, with topographic dominance at depth. Intense and vertically coherent upwelling and downwelling occur in the Southern Ocean, which are likely due to the interaction of the Antarctic Circumpolar Current and large-scale topographic features and are generally canceled out in the conventional zonally averaged results. These "elevators" at high latitudes connect the upper to the deep and abyssal oceans and working together with isopycnal mixing are likely a mechanism, in addition to the formation of deep and abyssal waters, for fast responses of the deep and abyssal oceans to the changing climate. Also, Eulerian and parameterized eddy-induced components are of opposite signs in numerous regions around the global ocean, particularly in the ocean interior away from surface and bottom. Nevertheless, residual vertical velocity is primarily determined by the Eulerian component, and related to winds and large-scale topographic features. The current estimates of vertical velocities can serve as a useful reference for investigating the vertical exchange of ocean properties and tracers, and its complex spatial structure ultimately permits regional tests of basic oceanographic concepts such as Sverdrup balance and coastal upwelling/downwelling.

Computer program for solving laminar, transitional, or turbulent compressible boundary-layer equations for two-dimensional and axisymmetric flow

NASA Technical Reports Server (NTRS)

Harris, J. E.; Blanchard, D. K.

1982-01-01

A numerical algorithm and computer program are presented for solving the laminar, transitional, or turbulent two dimensional or axisymmetric compressible boundary-layer equations for perfect-gas flows. The governing equations are solved by an iterative three-point implicit finite-difference procedure. The software, program VGBLP, is a modification of the approach presented in NASA TR R-368 and NASA TM X-2458, respectively. The major modifications are: (1) replacement of the fourth-order Runge-Kutta integration technique with a finite-difference procedure for numerically solving the equations required to initiate the parabolic marching procedure; (2) introduction of the Blottner variable-grid scheme; (3) implementation of an iteration scheme allowing the coupled system of equations to be converged to a specified accuracy level; and (4) inclusion of an iteration scheme for variable-entropy calculations. These modifications to the approach presented in NASA TR R-368 and NASA TM X-2458 yield a software package with high computational efficiency and flexibility. Turbulence-closure options include either two-layer eddy-viscosity or mixing-length models. Eddy conductivity is modeled as a function of eddy viscosity through a static turbulent Prandtl number formulation. Several options are provided for specifying the static turbulent Prandtl number. The transitional boundary layer is treated through a streamwise intermittency function which modifies the turbulence-closure model. This model is based on the probability distribution of turbulent spots and ranges from zero to unity for laminar and turbulent flow, respectively. Several test cases are presented as guides for potential users of the software.

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.

Eddy-driven stratification initiates North Atlantic spring phytoplankton blooms.

PubMed

Mahadevan, Amala; D'Asaro, Eric; Lee, Craig; Perry, Mary Jane

2012-07-06

Springtime phytoplankton blooms photosynthetically fix carbon and export it from the surface ocean at globally important rates. These blooms are triggered by increased light exposure of the phytoplankton due to both seasonal light increase and the development of a near-surface vertical density gradient (stratification) that inhibits vertical mixing of the phytoplankton. Classically and in current climate models, that stratification is ascribed to a springtime warming of the sea surface. Here, using observations from the subpolar North Atlantic and a three-dimensional biophysical model, we show that the initial stratification and resulting bloom are instead caused by eddy-driven slumping of the basin-scale north-south density gradient, resulting in a patchy bloom beginning 20 to 30 days earlier than would occur by warming.

A Dynamic Eddy Viscosity Model for the Shallow Water Equations Solved by Spectral Element and Discontinuous Galerkin Methods

NASA Astrophysics Data System (ADS)

Marras, Simone; Suckale, Jenny; Giraldo, Francis X.; Constantinescu, Emil

2016-04-01

We present the solution of the viscous shallow water equations where viscosity is built as a residual-based subgrid scale model originally designed for large eddy simulation of compressible [1] and stratified flows [2]. The necessity of viscosity for a shallow water model not only finds motivation from mathematical analysis [3], but is supported by physical reasoning as can be seen by an analysis of the energetics of the solution. We simulated the flow of an idealized wave as it hits a set of obstacles. The kinetic energy spectrum of this flow shows that, although the inviscid Galerkin solutions -by spectral elements and discontinuous Galerkin [4]- preserve numerical stability in spite of the spurious oscillations in the proximity of the wave fronts, the slope of the energy cascade deviates from the theoretically expected values. We show that only a sufficiently small amount of dynamically adaptive viscosity removes the unwanted high-frequency modes while preserving the overall sharpness of the solution. In addition, it yields a physically plausible energy decay. This work is motivated by a larger interest in the application of a shallow water model to the solution of tsunami triggered coastal flows. In particular, coastal flows in regions around the world where coastal parks made of mitigation hills of different sizes and configurations are considered as a means to deviate the power of the incoming wave. References [1] M. Nazarov and J. Hoffman (2013) "Residual-based artificial viscosity for simulation of turbulent compressible flow using adaptive finite element methods" Int. J. Numer. Methods Fluids, 71:339-357 [2] S. Marras, M. Nazarov, F. X. Giraldo (2015) "Stabilized high-order Galerkin methods based on a parameter-free dynamic SGS model for LES" J. Comput. Phys. 301:77-101 [3] J. F. Gerbeau and B. Perthame (2001) "Derivation of the viscous Saint-Venant system for laminar shallow water; numerical validation" Discrete Contin. Dyn. Syst. Ser. B, 1:89?102 [4] F. X. Giraldo and M. Restelli (2010) "High-order semi-implicit time-integrators for a triangular discontinuous Galerkin oceanic shallow water model. Int. J. Numer. Methods Fluids, 63:1077-1102

A new class of finite element variational multiscale turbulence models for incompressible magnetohydrodynamics

DOE PAGES

Sondak, D.; Shadid, J. N.; Oberai, A. A.; ...

2015-04-29

New large eddy simulation (LES) turbulence models for incompressible magnetohydrodynamics (MHD) derived from the variational multiscale (VMS) formulation for finite element simulations are introduced. The new models include the variational multiscale formulation, a residual-based eddy viscosity model, and a mixed model that combines both of these component models. Each model contains terms that are proportional to the residual of the incompressible MHD equations and is therefore numerically consistent. Moreover, each model is also dynamic, in that its effect vanishes when this residual is small. The new models are tested on the decaying MHD Taylor Green vortex at low and highmore » Reynolds numbers. The evaluation of the models is based on comparisons with available data from direct numerical simulations (DNS) of the time evolution of energies as well as energy spectra at various discrete times. Thus a numerical study, on a sequence of meshes, is presented that demonstrates that the large eddy simulation approaches the DNS solution for these quantities with spatial mesh refinement.« less

Large-Eddy Simulation of Internal Flow through Human Vocal Folds

NASA Astrophysics Data System (ADS)

Lasota, Martin; Šidlof, Petr

2018-06-01

The phonatory process occurs when air is expelled from the lungs through the glottis and the pressure drop causes flow-induced oscillations of the vocal folds. The flow fields created in phonation are highly unsteady and the coherent vortex structures are also generated. For accuracy it is essential to compute on humanlike computational domain and appropriate mathematical model. The work deals with numerical simulation of air flow within the space between plicae vocales and plicae vestibulares. In addition to the dynamic width of the rima glottidis, where the sound is generated, there are lateral ventriculus laryngis and sacculus laryngis included in the computational domain as well. The paper presents the results from OpenFOAM which are obtained with a large-eddy simulation using second-order finite volume discretization of incompressible Navier-Stokes equations. Large-eddy simulations with different subgrid scale models are executed on structured mesh. In these cases are used only the subgrid scale models which model turbulence via turbulent viscosity and Boussinesq approximation in subglottal and supraglottal area in larynx.

A family of dynamic models for large-eddy simulation

NASA Technical Reports Server (NTRS)

Carati, D.; Jansen, K.; Lund, T.

1995-01-01

Since its first application, the dynamic procedure has been recognized as an effective means to compute rather than prescribe the unknown coefficients that appear in a subgrid-scale model for Large-Eddy Simulation (LES). The dynamic procedure is usually used to determine the nondimensional coefficient in the Smagorinsky (1963) model. In reality the procedure is quite general and it is not limited to the Smagorinsky model by any theoretical or practical constraints. The purpose of this note is to consider a generalized family of dynamic eddy viscosity models that do not necessarily rely on the local equilibrium assumption built into the Smagorinsky model. By invoking an inertial range assumption, it will be shown that the coefficients in the new models need not be nondimensional. This additional degree of freedom allows the use of models that are scaled on traditionally unknown quantities such as the dissipation rate. In certain cases, the dynamic models with dimensional coefficients are simpler to implement, and allow for a 30% reduction in the number of required filtering operations.

Estimation of strait transport in the East China Sea

NASA Astrophysics Data System (ADS)

Moon, J.; Hirose, N.; Usui, N.; Tsujino, H.

2010-12-01

Volume transport through the major channels is still diverse in realistic eddy-resolving models. For instance, time-mean transport through the Tokara Strait is predicted as 16.9Sv by Maltrud and McClean (2005) or 36-72Sv by Hurlburt et al. (1996). However the difference may be decreased by constraining measurement data, i.e., data assimilation. The assimilated estimates from two different systems of Meteorological Research Institute and Kyushu University (MOVE-WNP and DREAMS_B) show realistic averages of 22-23Sv through the Tokara Strait. Inverse estimation of adjustable parameters implies that reduction of wind stress and strong vertical viscosity are crucial to prevent excessive transport and associated instabilities in a forward model. It is noted that both of the assimilated results show a deep northward flow of 3-4Sv through the Kerama Gap in Ryukyu Islands. The core depth (~500m) of this subsurface current is similar to Ryukyu Current. Further analysis shows coherent changes of Soya and Tsushima Warm Currents, which is consistent to the Okhotsk wind theory of Tsujino et al. (2008). On the other hand, the changes of Tsushima Strait transport are nearly independent from the Kuroshio or the Taiwan Warm Current.

Turbulent Mixing in Gravity Currents with Transverse Shear

NASA Astrophysics Data System (ADS)

White, Brian; Helfrich, Karl; Scotti, Alberto

2010-11-01

A parallel flow with horizontal shear and horizontal density gradient undergoes an intensification of the shear by gravitational tilting and stretching, rapidly breaking down into turbulence. Such flows have the potential for substantial mixing in estuaries and the coastal ocean. We present high-resolution numerical results for the mixing efficiency of these flows, which can be viewed as gravity currents with transverse shear, and contrast them with the well-studied case of stably stratified, homogeneous turbulence (uniform vertical density and velocity gradients). For a sheared gravity current, the buoyancy flux, turbulent Reynolds stress, and dissipation are well out of equilibrium. The total kinetic energy first increases as potential energy is transferred to the gravity current, but rapidly decays once turbulence sets in. Despite the non-equilibrium character, mixing efficiencies are slightly higher but qualitatively similar to homogeneous stratified turbulence. Efficiency decreases in the highly energetic regime where the dissipation rate is large compared with viscosity and stratification, ɛ/(νN^2)>100, further declining as turbulence decays and kinetic energy dissipation dominates the buoyancy flux. In general, the mixing rate, parameterized by a turbulent eddy diffusivity, increases with the strength of the transverse shear.

The Vertical Transport in the Ocean: a Pump Driven by Meso and Submesoscale Structures

NASA Astrophysics Data System (ADS)

Rosso, I.; Hogg, A. M.; Strutton, P. G.; Kiss, A. E.

2012-04-01

The upper ocean can be considered as a vehicle for the exchange of gases between air and the deep ocean. Furthermore, the transport of nutrients through the mixed layer occurs via a combination of biogeochemical and physical pumps; both of these mechanisms play a fundamental role in the carbon cycle. In the surface layer phytoplankton convert carbon dioxide into organic compounds using nutrients and light. Nutrients, which are depleted at the surface, can be transported into the mixed layer by vertical motion; recently, it has been discovered that this vertical transport is more often associated with submesoscale fronts of O(10) km (rather than inside mesoscale structures, of O(100) km, like eddies). At the submesoscale fronts, rates of O(100) m day-1 can emerge, particularly high compared to values of 10 m/day found at the mesoscales [M. Lèvy, et al. J. Mar. Res., 2001]. At this fine scale, the vertical transport of nutrients is highly effective, upwelling waters from the depth rich of nutrients and downwelling depleted waters from the surface. The fine-scale vertical transport mechanism has recently become of great interest, though is not completely understood. We investigate the dynamics and the transport of tracers at the meso and sub-mesoscales by running numerical simulations with a domain of 1024 km x 512 km x 1600 m, at 3 different resolutions: 8 km, 4km and 1km. We use the MIT general circulation model with free surface, linear bottom drag and free slip condition at the north and south walls. Non-linear 3rd order advection scheme and biharmonic viscosity are applied. Furthermore, the fluid is forced by a constant zonal wind stress. The flow is zonally periodic and presents an idealized topography. We started from an initial vertical stratification and run the model to reach an equilibrium flow state. A passive tracer is released after the equilibrium is reached. We investigate how the fine scales are affecting the rate of vertical transport and the distribution of the tracer. We show that also the presence of the topography has an impact in driving this transport. Furthermore, the model can run in a non-hydrostatic configuration, allowing us to investigate the effect of this parameterization on the transport.

Finite-mode spectral model of homogeneous and isotropic Navier-stokes turbulence: a rapidly depleted energy cascade.

PubMed

Lévêque, E; Koudella, C R

2001-04-30

An eddy-viscous term is added to Navier-Stokes dynamics at wave numbers k greater than the inflection point kc of the energy flux F(log(k)). The eddy viscosity is fixed so that the energy spectrum satisfies E(k) = E(kc) (k/kc)(-3) for k>kc. This resulting forcing induces a rapid depletion of the energy cascade at k>kc. It is observed numerically that the model reproduces turbulence energetics at k< or =kc and statistics of two-point velocity correlations at scales r>lambda (Taylor microscale). Compared to a direct numerical simulation of R(lambda) = 130 an equivalent run with the present model results in a gain of a factor 20 in CPU time.

A comparison of three approaches to compute the effective Reynolds number of the implicit large-eddy simulations

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

Zhou, Ye; Thornber, Ben

2016-04-12

Here, the implicit large-eddy simulation (ILES) has been utilized as an effective approach for calculating many complex flows at high Reynolds number flows. Richtmyer–Meshkov instability (RMI) induced flow can be viewed as a homogeneous decaying turbulence (HDT) after the passage of the shock. In this article, a critical evaluation of three methods for estimating the effective Reynolds number and the effective kinematic viscosity is undertaken utilizing high-resolution ILES data. Effective Reynolds numbers based on the vorticity and dissipation rate, or the integral and inner-viscous length scales, are found to be the most self-consistent when compared to the expected phenomenology andmore » wind tunnel experiments.« less

Image-based modelling of lateral magma flow: the Basement Sill, Antarctica

PubMed Central

Mirhadizadeh, Seyed

2017-01-01

The McMurdo Dry Valleys magmatic system, Antarctica, provides a world-class example of pervasive lateral magma flow on a continental scale. The lowermost intrusion (Basement Sill) offers detailed sections through the now frozen particle microstructure of a congested magma slurry. We simulated the flow regime in two and three dimensions using numerical models built on a finite-element mesh derived from field data. The model captures the flow behaviour of the Basement Sill magma over a viscosity range of 1–104 Pa s where the higher end (greater than or equal to 102 Pa s) corresponds to a magmatic slurry with crystal fractions varying between 30 and 70%. A novel feature of the model is the discovery of transient, low viscosity (less than or equal to 50 Pa s) high Reynolds number eddies formed along undulating contacts at the floor and roof of the intrusion. Numerical tracing of particle orbits implies crystals trapped in eddies segregate according to their mass density. Recovered shear strain rates (10−3–10−5 s−1) at viscosities equating to high particle concentrations (around more than 40%) in the Sill interior point to shear-thinning as an explanation for some types of magmatic layering there. Model transport rates for the Sill magmas imply a maximum emplacement time of ca 105 years, consistent with geochemical evidence for long-range lateral flow. It is a theoretically possibility that fast-flowing magma on a continental scale will be susceptible to planetary-scale rotational forces. PMID:28573002

Image-based modelling of lateral magma flow: the Basement Sill, Antarctica.

PubMed

Petford, Nick; Mirhadizadeh, Seyed

2017-05-01

The McMurdo Dry Valleys magmatic system, Antarctica, provides a world-class example of pervasive lateral magma flow on a continental scale. The lowermost intrusion (Basement Sill) offers detailed sections through the now frozen particle microstructure of a congested magma slurry. We simulated the flow regime in two and three dimensions using numerical models built on a finite-element mesh derived from field data. The model captures the flow behaviour of the Basement Sill magma over a viscosity range of 1-10 4  Pa s where the higher end (greater than or equal to 10 2  Pa s) corresponds to a magmatic slurry with crystal fractions varying between 30 and 70%. A novel feature of the model is the discovery of transient, low viscosity (less than or equal to 50 Pa s) high Reynolds number eddies formed along undulating contacts at the floor and roof of the intrusion. Numerical tracing of particle orbits implies crystals trapped in eddies segregate according to their mass density. Recovered shear strain rates (10 -3 -10 -5  s -1 ) at viscosities equating to high particle concentrations (around more than 40%) in the Sill interior point to shear-thinning as an explanation for some types of magmatic layering there. Model transport rates for the Sill magmas imply a maximum emplacement time of ca 10 5 years, consistent with geochemical evidence for long-range lateral flow. It is a theoretically possibility that fast-flowing magma on a continental scale will be susceptible to planetary-scale rotational forces.

Open ocean dead zones in the tropical North Atlantic Ocean

NASA Astrophysics Data System (ADS)

Karstensen, J.; Fiedler, B.; Schütte, F.; Brandt, P.; Körtzinger, A.; Fischer, G.; Zantopp, R.; Hahn, J.; Visbeck, M.; Wallace, D.

2015-04-01

Here we present first observations, from instrumentation installed on moorings and a float, of unexpectedly low (<2 μmol kg-1) oxygen environments in the open waters of the tropical North Atlantic, a region where oxygen concentration does normally not fall much below 40 μmol kg-1. The low-oxygen zones are created at shallow depth, just below the mixed layer, in the euphotic zone of cyclonic eddies and anticyclonic-modewater eddies. Both types of eddies are prone to high surface productivity. Net respiration rates for the eddies are found to be 3 to 5 times higher when compared with surrounding waters. Oxygen is lowest in the centre of the eddies, in a depth range where the swirl velocity, defining the transition between eddy and surroundings, has its maximum. It is assumed that the strong velocity at the outer rim of the eddies hampers the transport of properties across the eddies boundary and as such isolates their cores. This is supported by a remarkably stable hydrographic structure of the eddies core over periods of several months. The eddies propagate westward, at about 4 to 5 km day-1, from their generation region off the West African coast into the open ocean. High productivity and accompanying respiration, paired with sluggish exchange across the eddy boundary, create the "dead zone" inside the eddies, so far only reported for coastal areas or lakes. We observe a direct impact of the open ocean dead zones on the marine ecosystem as such that the diurnal vertical migration of zooplankton is suppressed inside the eddies.

Classical closure theory and Lam's interpretation of epsilon-RNG

NASA Technical Reports Server (NTRS)

Zhou, YE

1995-01-01

Lam's phenomenological epsilon-renormalization group (RNG) model is quite different from the other members of that group. It does not make use of the correspondence principle and the epsilon-expansion procedure. We demonstrate that Lam's epsilon-RNG model is essentially the physical space version of the classical closure theory in spectral space and consider the corresponding treatment of the eddy viscosity and energy backscatter.

An ocean large-eddy simulation of Langmuir circulations and convection in the surface mixed layer

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

Skyllingstad, E.D.; Denbo, D.W.

Numerical experiments were performed using a three-dimensional large-eddy simulation model of the ocean surface mixed layer that includes the Craik-Leibovich vortex force to parameterize the interaction of surface waves with mean currents. Results from the experiments show that the vortex force generates Langmuir circulations that can dominate vertical mixing. The simulated vertical velocity fields show linear, small-scale, coherent structures near the surface that extend downwind across the model domain. In the interior of the mixed layer, scales of motion increase to eddy sizes that are roughly equivalent to the mixed-layer depth. Cases with the vortex force have stronger circulations nearmore » the surface in contrast to cases with only heat flux and wind stress, particularly when the heat flux is positive. Calculations of the velocity variance and turbulence dissipation rates for cases with and without the vortex force, surface cooling, and wind stress indicate that wave-current interactions are a dominant mixing process in the upper mixed layer. Heat flux calculations show that the entrainment rate at the mixed-layer base can be up to two times greater when the vortex force is included. In a case with reduced wind stress, turbulence dissipation rates remained high near the surface because of the vortex force interaction with preexisting inertial currents. In deep mixed layers ({approximately}250 m) the simulations show that Langmuir circulations can vertically transport water 145 m during conditions of surface heating. Observations of turbulence dissipation rates and the vertical temperature structure support the model results. 42 refs., 20 figs., 21 tabs.« less

Quasi-laminar stability and sensitivity analyses for turbulent flows: Prediction of low-frequency unsteadiness and passive control

NASA Astrophysics Data System (ADS)

Mettot, Clément; Sipp, Denis; Bézard, Hervé

2014-04-01

This article presents a quasi-laminar stability approach to identify in high-Reynolds number flows the dominant low-frequencies and to design passive control means to shift these frequencies. The approach is based on a global linear stability analysis of mean-flows, which correspond to the time-average of the unsteady flows. Contrary to the previous work by Meliga et al. ["Sensitivity of 2-D turbulent flow past a D-shaped cylinder using global stability," Phys. Fluids 24, 061701 (2012)], we use the linearized Navier-Stokes equations based solely on the molecular viscosity (leaving aside any turbulence model and any eddy viscosity) to extract the least stable direct and adjoint global modes of the flow. Then, we compute the frequency sensitivity maps of these modes, so as to predict before hand where a small control cylinder optimally shifts the frequency of the flow. In the case of the D-shaped cylinder studied by Parezanović and Cadot [J. Fluid Mech. 693, 115 (2012)], we show that the present approach well captures the frequency of the flow and recovers accurately the frequency control maps obtained experimentally. The results are close to those already obtained by Meliga et al., who used a more complex approach in which turbulence models played a central role. The present approach is simpler and may be applied to a broader range of flows since it is tractable as soon as mean-flows — which can be obtained either numerically from simulations (Direct Numerical Simulation (DNS), Large Eddy Simulation (LES), unsteady Reynolds-Averaged-Navier-Stokes (RANS), steady RANS) or from experimental measurements (Particle Image Velocimetry - PIV) — are available. We also discuss how the influence of the control cylinder on the mean-flow may be more accurately predicted by determining an eddy-viscosity from numerical simulations or experimental measurements. From a technical point of view, we finally show how an existing compressible numerical simulation code may be used in a black-box manner to extract the global modes and sensitivity maps.

Water Vapor Feedbacks to Climate Change

NASA Technical Reports Server (NTRS)

Rind, David

1999-01-01

The response of water vapor to climate change is investigated through a series of model studies with varying latitudinal temperature gradients, mean temperatures, and ultimately, actual climate change configurations. Questions to be addressed include: what role does varying convection have in water vapor feedback; do Hadley Circulation differences result in differences in water vapor in the upper troposphere; and, does increased eddy energy result in greater eddy vertical transport of water vapor in varying climate regimes?

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

Atmospheric Moisture Budget and Spatial Resolution Dependence of Precipitation Extremes in Aquaplanet Simulations

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

Yang, Qing; Leung, Lai-Yung R.; Rauscher, Sara

This study investigates the resolution dependency of precipitation extremes in an aqua-planet framework. Strong resolution dependency of precipitation extremes is seen over both tropics and extra-tropics, and the magnitude of this dependency also varies with dynamical cores. Moisture budget analyses based on aqua-planet simulations with the Community Atmosphere Model (CAM) using the Model for Prediction Across Scales (MPAS) and High Order Method Modeling Environment (HOMME) dynamical cores but the same physics parameterizations suggest that during precipitation extremes moisture supply for surface precipitation is mainly derived from advective moisture convergence. The resolution dependency of precipitation extremes mainly originates from advective moisturemore » transport in the vertical direction. At most vertical levels over the tropics and in the lower atmosphere over the subtropics, the vertical eddy transport of mean moisture field dominates the contribution to precipitation extremes and its resolution dependency. Over the subtropics, the source of moisture, its associated energy, and the resolution dependency during extremes are dominated by eddy transport of eddies moisture at the mid- and upper-troposphere. With both MPAS and HOMME dynamical cores, the resolution dependency of the vertical advective moisture convergence is mainly explained by dynamical changes (related to vertical velocity or omega), although the vertical gradients of moisture act like averaging kernels to determine the sensitivity of the overall resolution dependency to the changes in omega at different vertical levels. The natural reduction of variability with coarser resolution, represented by areal data averaging (aggregation) effect, largely explains the resolution dependency in omega. The thermodynamic changes, which likely result from non-linear feedback in response to the large dynamical changes, are small compared to the overall changes in dynamics (omega). However, after excluding the data aggregation effect in omega, thermodynamic changes become relatively significant in offsetting the effect of dynamics leading to reduce differences between the simulated and aggregated results. Compared to MPAS, the simulated stronger vertical motion with HOMME also results in larger resolution dependency. Compared to the simulation at fine resolution, the vertical motion during extremes is insufficiently resolved/parameterized at the coarser resolution even after accounting for the natural reduction in variability with coarser resolution, and this is more distinct in the simulation with HOMME. To reduce uncertainties in simulated precipitation extremes, future development in cloud parameterizations must address their sensitivity to spatial resolution as well as dynamical cores.« less

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.

Interface instabilities during displacements of two miscible fluids in a vertical pipe

NASA Astrophysics Data System (ADS)

Scoffoni, J.; Lajeunesse, E.; Homsy, G. M.

2001-03-01

We study experimentally the downward vertical displacement of one miscible fluid by another in a vertical pipe at sufficiently high velocities for diffusive effects to be negligible. For certain viscosity ratios and flow rates, the interface between the two fluids can destabilize. We determine the dimensionless flow rate Uc above which the instability is triggered and its dependence on the viscous ratio M, resulting in a stability map Uc=Uc(M). Two different instability modes have been observed: an asymmetric "corkscrew" mode and an axisymmetric one. We remark that the latter is always eventually disturbed by "corkscrew" type instabilities. We speculate that these instabilities are driven by the viscosity stratification and are analogous to those already observed in core annular flows of immiscible fluids.

Local atmospheric response to warm mesoscale ocean eddies in the Kuroshio-Oyashio Confluence region.

PubMed

Sugimoto, Shusaku; Aono, Kenji; Fukui, Shin

2017-09-19

In the extratropical regions, surface winds enhance upward heat release from the ocean to atmosphere, resulting in cold surface ocean: surface ocean temperature is negatively correlated with upward heat flux. However, in the western boundary currents and eddy-rich regions, the warmer surface waters compared to surrounding waters enhance upward heat release-a positive correlation between upward heat release and surface ocean temperature, implying that the ocean drives the atmosphere. The atmospheric response to warm mesoscale ocean eddies with a horizontal extent of a few hundred kilometers remains unclear because of a lack of observations. By conducting regional atmospheric model experiments, we show that, in the Kuroshio-Oyashio Confluence region, wintertime warm eddies heat the marine atmospheric boundary layer (MABL), and accelerate westerly winds in the near-surface atmosphere via the vertical mixing effect, leading to wind convergence around the eastern edge of eddies. The warm-eddy-induced convergence forms local ascending motion where convective precipitation is enhanced, providing diabatic heating to the atmosphere above MABL. Our results indicate that warm eddies affect not only near-surface atmosphere but also free atmosphere, and possibly synoptic atmospheric variability. A detailed understanding of warm eddy-atmosphere interaction is necessary to improve in weather and climate projections.

Three-dimensional structure of a swarm of the salp Thalia democratica within a cold-core eddy off southeast Australia

NASA Astrophysics Data System (ADS)

Everett, J. D.; Baird, M. E.; Suthers, I. M.

2011-12-01

Swarms of the salp Thalia democratica periodically occur off southeast Australia following the austral spring bloom of phytoplankton. In October 2008 a filament of upwelled water was advected south by the adjacent East Australian Current and formed a 30 km diameter cold-core eddy (CCE). The three-dimensional structure of a subsurface swarm of T. democratica within the eddy was examined using both oblique and vertical hauls and an optical plankton counter (OPC) deployed on a towed body. The CCE displayed distinct uplift of the nutricline and elevated fluorescence. Net samples show the zooplankton community was dominated by T. democratica, comprising 73%-88% of zooplankton abundance. The size distribution of T. democratica measured from net samples was 0.5-5 mm and was used to interpret the OPC transects, which showed the swarm formed a 15 km diameter disc located 20-40 m deep in the center of the eddy. The maximum salp abundance was in the pycnocline and coincided with the subsurface fluorescence maximum. The mean abundance of T. democratica size particles within the disc was 5003 individuals m-3 (ind. m-3), contrasted with only 604 ind. m-3 at the outer edge of the eddy. The vertically concentrated and horizontally constrained disc-shaped salp swarm occurred at the interface of salp-bearing inner shelf water and nutrient-rich upwelled water in a CCE. The physical processes that formed the CCE on the inshore edge of the western boundary current led to the largest density of salps recorded.

Simulation study of disruption characteristics in KSTAR

NASA Astrophysics Data System (ADS)

Lee, Jongkyu; Kim, J. Y.; Kessel, C. E.; Poli, F.

2012-10-01

A detailed simulation study of disruption in KSTAR had been performed using the Tokamak Simulation Code(TSC) [1] during the initial design phase of KSTAR [2]. Recently, however, a partial modification in the structure of passive plate was made in relation to reduce eddy current and increase the efficiency of control of vertical position. A substantial change can then occur in disruption characteristics and plasma behavior during disruption due to changes in passive plate structure. Because of this, growth rate of vertical instability is expected to be increased and eddy current and its associated electomagnetic force are expected to be reduced. To check this in more detail, a new simulation study is here given with modified passive plate structure of KSTAR. In particular, modeling of vertical disruption that is vertical displacement event (VDE) was carried out. We calculated vertical growth rate for a drift phase of plasma and electromagnetic force acting on PFC structures and compared the results between in a new model and an old model. [4pt] [1] S.C. Jardin, N. Pomphrey and J. Delucia, J. Comp. Phys. 66, 481 (1986).[0pt] [2] J.Y. Kim, S.Y. Cho and KSTAR Team, Disruption load analysis on KSTAR PFC structures, J. Accel. Plasma Res. 5, 149 (2000).

A numerical modeling study of the East Australian Current encircling and overwashing a warm-core eddy

NASA Astrophysics Data System (ADS)

MacDonald, H. S.; Roughan, M.; Baird, M. E.; Wilkin, J.

2013-01-01