Low-Dissipation Advection Schemes Designed for Large Eddy Simulations of Hypersonic Propulsion Systems

NASA Technical Reports Server (NTRS)

White, Jeffrey A.; Baurle, Robert A.; Fisher, Travis C.; Quinlan, Jesse R.; Black, William S.

2012-01-01

The 2nd-order upwind inviscid flux scheme implemented in the multi-block, structured grid, cell centered, finite volume, high-speed reacting flow code VULCAN has been modified to reduce numerical dissipation. This modification was motivated by the desire to improve the codes ability to perform large eddy simulations. The reduction in dissipation was accomplished through a hybridization of non-dissipative and dissipative discontinuity-capturing advection schemes that reduces numerical dissipation while maintaining the ability to capture shocks. A methodology for constructing hybrid-advection schemes that blends nondissipative fluxes consisting of linear combinations of divergence and product rule forms discretized using 4th-order symmetric operators, with dissipative, 3rd or 4th-order reconstruction based upwind flux schemes was developed and implemented. A series of benchmark problems with increasing spatial and fluid dynamical complexity were utilized to examine the ability of the candidate schemes to resolve and propagate structures typical of turbulent flow, their discontinuity capturing capability and their robustness. A realistic geometry typical of a high-speed propulsion system flowpath was computed using the most promising of the examined schemes and was compared with available experimental data to demonstrate simulation fidelity.

A Partially-Stirred Batch Reactor Model for Under-Ventilated Fire Dynamics

NASA Astrophysics Data System (ADS)

McDermott, Randall; Weinschenk, Craig

2013-11-01

A simple discrete quadrature method is developed for closure of the mean chemical source term in large-eddy simulations (LES) and implemented in the publicly available fire model, Fire Dynamics Simulator (FDS). The method is cast as a partially-stirred batch reactor model for each computational cell. The model has three distinct components: (1) a subgrid mixing environment, (2) a mixing model, and (3) a set of chemical rate laws. The subgrid probability density function (PDF) is described by a linear combination of Dirac delta functions with quadrature weights set to satisfy simple integral constraints for the computational cell. It is shown that under certain limiting assumptions, the present method reduces to the eddy dissipation concept (EDC). The model is used to predict carbon monoxide concentrations in direct numerical simulation (DNS) of a methane slot burner and in LES of an under-ventilated compartment fire.

Foucault Dissipation in a Rolling Cylinder: A Webcam Quantitative Study

ERIC Educational Resources Information Center

Bonanno, A.; Bozzo, G.; Camarca, M.; Sapia, P.

2011-01-01

In this paper we present an experimental strategy to measure the micro power dissipation due to Foucault "eddy" currents in a copper cylinder rolling on two parallel conductive rails in the presence of a magnetic field. Foucault power dissipation is obtained from kinematical measurements carried out by using a common PC webcam and video analysis…

Statistical turbulence theory and turbulence phenomenology

NASA Technical Reports Server (NTRS)

Herring, J. R.

1973-01-01

The application of deductive turbulence theory for validity determination of turbulence phenomenology at the level of second-order, single-point moments is considered. Particular emphasis is placed on the phenomenological formula relating the dissipation to the turbulence energy and the Rotta-type formula for the return to isotropy. Methods which deal directly with most or all the scales of motion explicitly are reviewed briefly. The statistical theory of turbulence is presented as an expansion about randomness. Two concepts are involved: (1) a modeling of the turbulence as nearly multipoint Gaussian, and (2) a simultaneous introduction of a generalized eddy viscosity operator.

On the Roles of Upper- versus Lower-level Thermal Forcing in Shifting the Eddy-Driven Jet

NASA Astrophysics Data System (ADS)

Zhang, Y.; Nie, Y.; Chen, G.; Yang, X. Q.

2017-12-01

One most drastic atmospheric change in the global warming scenario is the increase in temperature over tropical upper-troposphere and polar surface. The strong warming over those two area alters the spacial distributions of the baroclinicity in the upper-troposphere of subtropics and in the lower-level of subpolar region, with competing effects on the mid-latitude atmospheric circulation. The final destination of the eddy-driven jet in future climate could be "a tug of war" between the impacts of such upper- versus lower-level thermal forcing. In this study, the roles of upper- versus lower-level thermal forcing in shifting the eddy-driven jet are investigated using a nonlinear multi-level quasi-geostrophic channel model. All of our sensitivity experiments show that the latitudinal position of the eddy-driven jet is more sensitive to the upper-level thermal forcing. Such upper-level dominance over the lower-level forcing can be attributed to the different mechanisms through which eddy-driven jet responses to them. The upper-level thermal forcing induces a jet shift mainly by affecting the baroclinic generation of eddies, which supports the latitudinal shift of the eddy momentum flux convergence. The jet response to the lower-level thermal forcing, however, is strongly "eddy dissipation control". The lower-level forcing, by changing the baroclinicity in the lower troposphere, induces a direct thermal zonal wind response in the upper level thus modifies the nonlinear wave breaking and the resultant irreversible eddy mixing, which amplifies the latitudinal shift of the eddy-driven jet. Whether the eddy response is "generation control" or "dissipation control" may strongly depend on the eddy behavior in its baroclinic processes. Only the anomalous eddy generation that penetrates into the upper troposphere can have a striking impact on the eddy momentum flux, which pushes the jet shift more efficiently and dominates the eddy response.

Onshore Wind Stress and Buoyancy Flux Observed on a Dissipative Mediterranean Beach

DTIC Science & Technology

2015-12-01

a climatologically Mediterranean coastline to explore the wind stress and buoyancy flux. An eddy covariance system was deployed in the intertidal... climatologically Mediterranean coastline to explore the wind stress and buoyancy flux. An eddy covariance system was deployed in the intertidal zone

Interrelationship of Cn2 & Eddy Dissipation rate based on Scintillometer and Doppler Lidar observations in complex terrain during the Perdigao Campaign 2017

NASA Astrophysics Data System (ADS)

Creegan, E. D.; Krishnamurthy, R.; Hocut, C. M.; Pattantyus, A.; Leo, L. S.; Wang, Y.; Fernando, H. J.; Bariteau, L.

2017-12-01

The Perdigao campaign is a joint EU/US science project designed to provide information on flow field(s) over complex terrain and through wind turbines at unprecedented high spatial and temporal resolution. The goal is to improve wind energy physics and overcome the current deficiencies of wind resource models. Topographically the Perdigao location is an expansion of the "double hill in crossflow", consisting of two parallel ridges along the NW-SE direction. The site was heavily instrumented with an array of towers (with multiple transects along the valley and across two ridges) and a large suite of ground based and aerial remote sensing platforms. On the outflow side of the NW ridge a scintillometer was emplaced with the line-of-sight (LOS) running adjacent to the towers comprising the NE transect from the ridgetop down to the base. Scanning lidars were placed at both ends of this LOS. Other instruments included a tethered lifting system (TLS), sodar, microwave radiometer, an energy budget flux tower and radiosonde releases. Scintillomoter data provides a quantitative measure of the intensity of optical turbulence, through the refractive index structure parameter, Cn2, where averaged Cn2 is often determined as a function of local differences in temperature, moisture, and wind velocity at discrete points. The refractive index structure parameter is also a function of the inner (dissipation) and outer (energy producing) turbulent scales. The scintillometer directly gives path averaged Cn2 and Eddy Dissipation rate along the LOS. Coplanar scans along the same path were synchronized using two scanning coherent Doppler lidars. Algorithms have been developed to estimate both eddy dissipation rate and Cn2 from Doppler lidar data effectively creating a new lidar data product. Additionally, from TLS measurements, Cn2 and dissipation rate are calculated using the high frequency spectra of the hot-wire sensor. In this work, measurements of Cn2 and Eddy Dissipation rate between multiple Doppler lidars, scintillometer and TLS are compared and the relationship between refractive index structure parameter and turbulence is explored. The effect of optical turbulence under various atmospheric conditions in complex terrain will be investigated.

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.

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

Davidovits, Seth; Fisch, Nathaniel J.

Here we report compression of turbulent plasma can amplify the turbulent kinetic energy, if the compression is fast compared to the viscous dissipation time of the turbulent eddies. A sudden viscous dissipation mechanism is demonstrated, whereby this amplified turbulent kinetic energy is rapidly converted into thermal energy, suggesting a new paradigm for fast ignition inertial fusion.

Effect of mean velocity shear on the dissipation rate of turbulent kinetic energy

NASA Technical Reports Server (NTRS)

Yoshizawa, Akira; Liou, Meng-Sing

1992-01-01

The dissipation rate of turbulent kinetic energy in incompressible turbulence is investigated using a two-scale DIA. The dissipation rate is shown to consist of two parts; one corresponds to the dissipation rate used in the current turbulence models of eddy-viscosity type, and another comes from the viscous effect that is closely connected with mean velocity shear. This result can elucidate the physical meaning of the dissipation rate used in the current turbulence models and explain part of the discrepancy in the near-wall dissipation rates between the current turbulence models and direct numerical simulation of the Navier-Stokes equation.

Derivation of revised formulae for eddy viscous forces used in the ocean general circulation model

NASA Technical Reports Server (NTRS)

Chou, Ru Ling

1988-01-01

Presented is a re-derivation of the eddy viscous dissipation tensor commonly used in present oceanographic general circulation models. When isotropy is imposed, the currently-used form of the tensor fails to return to the laplacian operator. In this paper, the source of this error is identified in a consistent derivation of the tensor in both rectangular and earth spherical coordinates, and the correct form of the eddy viscous tensor is presented.

Eddy currents in a conducting sphere

NASA Technical Reports Server (NTRS)

Bergman, John; Hestenes, David

1986-01-01

This report analyzes the eddy current induced in a solid conducting sphere by a sinusoidal current in a circular loop. Analytical expressions for the eddy currents are derived as a power series in the vectorial displacement of the center of the sphere from the axis of the loop. These are used for first order calculations of the power dissipated in the sphere and the force and torque exerted on the sphere by the electromagnetic field of the loop.

Turbulent convection in geostrophic circulation with wind and buoyancy forcing

NASA Astrophysics Data System (ADS)

Sohail, Taimoor; Gayen, Bishakhdatta; Hogg, Andy

2017-11-01

We conduct a direct numerical simulation of geostrophic circulation forced by surface wind and buoyancy to model a circumpolar ocean. The imposed buoyancy forcing (represented by Rayleigh number) drives a zonal current and supports small-scale convection in the buoyancy destabilizing region. In addition, we observe eddy activity which transports heat southward, supporting a large amount of heat uptake. Increasing wind stress enhances the meridional buoyancy gradient, triggering more eddy activity inside the boundary layer. Therefore, heat uptake increases with higher wind stress. The majority of dissipation is confined within the surface boundary layer, while mixing is dominant inside the convective plume and the buoyancy destabilizing region of the domain. The relative strength of the mixing and dissipation in the system can be expressed by mixing efficiency. This study finds that mixing is much greater than viscous dissipation, resulting in higher values of mixing efficiency than previously used. Supported by Australian Research Council Grant DP140103706.

Implicit and explicit subgrid-scale modeling in discontinuous Galerkin methods for large-eddy simulation

NASA Astrophysics Data System (ADS)

Fernandez, Pablo; Nguyen, Ngoc-Cuong; Peraire, Jaime

2017-11-01

Over the past few years, high-order discontinuous Galerkin (DG) methods for Large-Eddy Simulation (LES) have emerged as a promising approach to solve complex turbulent flows. Despite the significant research investment, the relation between the discretization scheme, the Riemann flux, the subgrid-scale (SGS) model and the accuracy of the resulting LES solver remains unclear. In this talk, we investigate the role of the Riemann solver and the SGS model in the ability to predict a variety of flow regimes, including transition to turbulence, wall-free turbulence, wall-bounded turbulence, and turbulence decay. The Taylor-Green vortex problem and the turbulent channel flow at various Reynolds numbers are considered. Numerical results show that DG methods implicitly introduce numerical dissipation in under-resolved turbulence simulations and, even in the high Reynolds number limit, this implicit dissipation provides a more accurate representation of the actual subgrid-scale dissipation than that by explicit models.

Quasi-geostrophic free mode models of long-lived Jovian eddies: Forcing mechanisms and crucial observational tests

NASA Technical Reports Server (NTRS)

Read, P. L.

1986-01-01

Observations of Jupiter and Saturn long-lived eddies, such as Jupiter's Great Red Spot and White Ovals, are presently compared with laboratory experiments and corresponding numerical simulations for free thermal convection in a rotating fluid that is subject to horizontal differential heating and cooling. Difficulties in determining the essential processes maintaining and dissipating stable eddies, on the basis of global energy budget studies, are discussed; such difficulties do not arise in considerations of the flow's potential vorticity budget. On Jupiter, diabatically forced and transient eddy-driven flows primarily differ in the implied role of transient eddies in transporting potential vorticity across closed geostrophic streamlines in the time mean.

On the unified estimation of turbulence eddy dissipation rate using Doppler cloud radars and lidars: Radar and Lidar Turbulence Estimation

DOE PAGES

Borque, Paloma; Luke, Edward; Kollias, Pavlos

2016-05-27

Coincident profiling observations from Doppler lidars and radars are used to estimate the turbulence energy dissipation rate (ε) using three different data sources: (i) Doppler radar velocity (DRV), (ii) Doppler lidar velocity (DLV), and (iii) Doppler radar spectrum width (DRW) measurements. Likewise, the agreement between the derived ε estimates is examined at the cloud base height of stratiform warm clouds. Collocated ε estimates based on power spectra analysis of DRV and DLV measurements show good agreement (correlation coefficient of 0.86 and 0.78 for both cases analyzed here) during both drizzling and nondrizzling conditions. This suggests that unified (below and abovemore » cloud base) time-height estimates of ε in cloud-topped boundary layer conditions can be produced. This also suggests that eddy dissipation rate can be estimated throughout the cloud layer without the constraint that clouds need to be nonprecipitating. Eddy dissipation rate estimates based on DRW measurements compare well with the estimates based on Doppler velocity but their performance deteriorates as precipitation size particles are introduced in the radar volume and broaden the DRW values. And, based on this finding, a methodology to estimate the Doppler spectra broadening due to the spread of the drop size distribution is presented. Furthermore, the uncertainties in ε introduced by signal-to-noise conditions, the estimation of the horizontal wind, the selection of the averaging time window, and the presence of precipitation are discussed in detail.« less

On the unified estimation of turbulence eddy dissipation rate using Doppler cloud radars and lidars: Radar and Lidar Turbulence Estimation

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

Borque, Paloma; Luke, Edward; Kollias, Pavlos

Coincident profiling observations from Doppler lidars and radars are used to estimate the turbulence energy dissipation rate (ε) using three different data sources: (i) Doppler radar velocity (DRV), (ii) Doppler lidar velocity (DLV), and (iii) Doppler radar spectrum width (DRW) measurements. Likewise, the agreement between the derived ε estimates is examined at the cloud base height of stratiform warm clouds. Collocated ε estimates based on power spectra analysis of DRV and DLV measurements show good agreement (correlation coefficient of 0.86 and 0.78 for both cases analyzed here) during both drizzling and nondrizzling conditions. This suggests that unified (below and abovemore » cloud base) time-height estimates of ε in cloud-topped boundary layer conditions can be produced. This also suggests that eddy dissipation rate can be estimated throughout the cloud layer without the constraint that clouds need to be nonprecipitating. Eddy dissipation rate estimates based on DRW measurements compare well with the estimates based on Doppler velocity but their performance deteriorates as precipitation size particles are introduced in the radar volume and broaden the DRW values. And, based on this finding, a methodology to estimate the Doppler spectra broadening due to the spread of the drop size distribution is presented. Furthermore, the uncertainties in ε introduced by signal-to-noise conditions, the estimation of the horizontal wind, the selection of the averaging time window, and the presence of precipitation are discussed in detail.« less

Modeling and calculation of turbulent lifted diffusion flames

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

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

1994-01-01

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

Large eddy simulation of turbine wakes using higher-order methods

NASA Astrophysics Data System (ADS)

Deskos, Georgios; Laizet, Sylvain; Piggott, Matthew D.; Sherwin, Spencer

2017-11-01

Large eddy simulations (LES) of a horizontal-axis turbine wake are presented using the well-known actuator line (AL) model. The fluid flow is resolved by employing higher-order numerical schemes on a 3D Cartesian mesh combined with a 2D Domain Decomposition strategy for an efficient use of supercomputers. In order to simulate flows at relatively high Reynolds numbers for a reasonable computational cost, a novel strategy is used to introduce controlled numerical dissipation to a selected range of small scales. The idea is to mimic the contribution of the unresolved small-scales by imposing a targeted numerical dissipation at small scales when evaluating the viscous term of the Navier-Stokes equations. The numerical technique is shown to behave similarly to the traditional eddy viscosity sub-filter scale models such as the classic or the dynamic Smagorinsky models. The results from the simulations are compared to experimental data for a Reynolds number scaled by the diameter equal to ReD =1,000,000 and both the time-averaged stream wise velocity and turbulent kinetic energy (TKE) are showing a good overall agreement. At the end, suggestions for the amount of numerical dissipation required by our approach are made for the particular case of horizontal-axis turbine wakes.

Seasonality of eddy kinetic energy in an eddy permitting global climate model

NASA Astrophysics Data System (ADS)

Uchida, Takaya; Abernathey, Ryan; Smith, Shafer

2017-10-01

We examine the seasonal cycle of upper-ocean mesoscale turbulence in a high resolution CESM climate simulation. The ocean model component (POP) has 0.1° resolution, mesoscale resolving at low and middle latitudes. Seasonally and regionally resolved wavenumber power spectra are calculated for sea-surface eddy kinetic energy (EKE). Although the interpretation of the spectral slopes in terms of turbulence theory is complicated by the strong presence of dissipation and the narrow inertial range, the EKE spectra consistently show higher power at small scales during winter throughout the ocean. Potential hypotheses for this seasonality are investigated. Diagnostics of baroclinc energy conversion rates and evidence from linear quasigeostrophic stability analysis indicate that seasonally varying mixed-layer instability is responsible for the seasonality in EKE. The ability of this climate model, which is not considered submesoscale resolving, to produce mixed layer instability although damped by dissipation, demonstrates the ubiquity and robustness of this process for modulating upper ocean EKE.

A Parallel, Finite-Volume Algorithm for Large-Eddy Simulation of Turbulent Flows

NASA Technical Reports Server (NTRS)

Bui, Trong T.

1999-01-01

A parallel, finite-volume algorithm has been developed for large-eddy simulation (LES) of compressible turbulent flows. This algorithm includes piecewise linear least-square reconstruction, trilinear finite-element interpolation, Roe flux-difference splitting, and second-order MacCormack time marching. Parallel implementation is done using the message-passing programming model. In this paper, the numerical algorithm is described. To validate the numerical method for turbulence simulation, LES of fully developed turbulent flow in a square duct is performed for a Reynolds number of 320 based on the average friction velocity and the hydraulic diameter of the duct. Direct numerical simulation (DNS) results are available for this test case, and the accuracy of this algorithm for turbulence simulations can be ascertained by comparing the LES solutions with the DNS results. The effects of grid resolution, upwind numerical dissipation, and subgrid-scale dissipation on the accuracy of the LES are examined. Comparison with DNS results shows that the standard Roe flux-difference splitting dissipation adversely affects the accuracy of the turbulence simulation. For accurate turbulence simulations, only 3-5 percent of the standard Roe flux-difference splitting dissipation is needed.

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.

Large eddy simulation in a rotary blood pump: Viscous shear stress computation and comparison with unsteady Reynolds-averaged Navier-Stokes simulation.

PubMed

Torner, Benjamin; Konnigk, Lucas; Hallier, Sebastian; Kumar, Jitendra; Witte, Matthias; Wurm, Frank-Hendrik

2018-06-01

Numerical flow analysis (computational fluid dynamics) in combination with the prediction of blood damage is an important procedure to investigate the hemocompatibility of a blood pump, since blood trauma due to shear stresses remains a problem in these devices. Today, the numerical damage prediction is conducted using unsteady Reynolds-averaged Navier-Stokes simulations. Investigations with large eddy simulations are rarely being performed for blood pumps. Hence, the aim of the study is to examine the viscous shear stresses of a large eddy simulation in a blood pump and compare the results with an unsteady Reynolds-averaged Navier-Stokes simulation. The simulations were carried out at two operation points of a blood pump. The flow was simulated on a 100M element mesh for the large eddy simulation and a 20M element mesh for the unsteady Reynolds-averaged Navier-Stokes simulation. As a first step, the large eddy simulation was verified by analyzing internal dissipative losses within the pump. Then, the pump characteristics and mean and turbulent viscous shear stresses were compared between the two simulation methods. The verification showed that the large eddy simulation is able to reproduce the significant portion of dissipative losses, which is a global indication that the equivalent viscous shear stresses are adequately resolved. The comparison with the unsteady Reynolds-averaged Navier-Stokes simulation revealed that the hydraulic parameters were in agreement, but differences for the shear stresses were found. The results show the potential of the large eddy simulation as a high-quality comparative case to check the suitability of a chosen Reynolds-averaged Navier-Stokes setup and turbulence model. Furthermore, the results lead to suggest that large eddy simulations are superior to unsteady Reynolds-averaged Navier-Stokes simulations when instantaneous stresses are applied for the blood damage prediction.

Retrieving clear-air turbulence information from regular commercial aircraft using Mode-S and ADS-B broadcast

NASA Astrophysics Data System (ADS)

Kopeć, J. M.; Kwiatkowski, K.; de Haan, S.; Malinowski, S. P.

2015-11-01

Navigational information broadcast by commercial aircraft in the form of Mode-S and ADS-B messages can be considered a new and valid source of upper air turbulence measurements. A set of three processing methods is proposed and analysed using a quality record of turbulence encounters made by a research aircraft. The proposed methods are based on processing the vertical acceleration or the background wind into the eddy dissipation rate. All the necessary parameters are conveyed in the Mode-S/ADS-B messages. The comparison of the results of application of the processing against a reference eddy dissipation rate obtained using on-board accelerometer indicate a significant potential of those methods. The advantages and limitation of the presented approaches are discussed.

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

NASA Astrophysics Data System (ADS)

Silvis, Maurits; Trias, Xavier; Abkar, Mahdi; Bae, Hyunji Jane; Lozano-Duran, Adrian; Verstappen, Roel

2016-11-01

This paper discusses subgrid models for large-eddy simulation of anisotropic flows using anisotropic grids. In particular, we are looking into ways to model not only the subgrid dissipation, but also transport processes, since these are expected to play an important role in rotating turbulent flows. We therefore consider subgrid-scale models of the form τ = - 2νt S +μt (SΩ - ΩS) , where the eddy-viscosity νt is given by the minimum-dissipation model, μt represents a transport coefficient; S is the symmetric part of the velocity gradient and Ω the skew-symmetric part. To incorporate the effect of mesh anisotropy the filter length is taken in such a way that it minimizes the difference between the turbulent stress in physical and computational space, where the physical space is covered by an anisotropic mesh and the computational space is isotropic. The resulting model is successfully tested for rotating homogeneous isotropic turbulence and rotating plane-channel flows. The research was largely carried out during the CTR SP 2016. M.S, and R.V. acknowledge the financial support to attend this Summer Program.

Baroclinic Adjustment of the Eddy-Driven Jet

NASA Astrophysics Data System (ADS)

Novak, Lenka; Ambaum, Maarten H. P.; Harvey, Ben J.

2017-04-01

The prediction of poleward shift in the midlatitude eddy-driven jets due to anthropogenic climate change is now a robust feature of climate models, but the magnitude of this shift or the processes responsible for it are less certain. This uncertainty comes from the complex response in storm tracks to large-scale forcing and their nonlinear modulation of the jet. This study uses global circulation models to reveal a relationship between eddy growth rate (referred to as baroclinicity) and eddy activity, whereby baroclinicity responds most rapidly to an eddy-dissipating forcing whereas eddy activity responds most rapidly to a baroclinicity-replenishing forcing. This nonlinearity can be generally explained using a two-dimensional dynamical system essentially describing the baroclinic adjustment as a predator-prey relationship. Despite this nonlinearity, the barotropic changes in the eddy-driven jet appear to be of a comparable magnitude for the ranges of both types of forcing tested in this study. It is implied that while changes in eddy activity or baroclinicity may indicate the sign of latitudinal jet shifting, the precise magnitude of this shifting is a result of a balance between these two quantities.

Improving Parameterization of Entrainment Rate for Shallow Convection with Aircraft Measurements and Large-Eddy Simulation

DOE PAGES

Lu, Chunsong; Liu, Yangang; Zhang, Guang J.; ...

2016-02-01

This work examines the relationships of entrainment rate to vertical velocity, buoyancy, and turbulent dissipation rate by applying stepwise principal component regression to observational data from shallow cumulus clouds collected during the Routine AAF [Atmospheric Radiation Measurement (ARM) Aerial Facility] Clouds with Low Optical Water Depths (CLOWD) Optical Radiative Observations (RACORO) field campaign over the ARM Southern Great Plains (SGP) site near Lamont, Oklahoma. The cumulus clouds during the RACORO campaign simulated using a large eddy simulation (LES) model are also examined with the same approach. The analysis shows that a combination of multiple variables can better represent entrainment ratemore » in both the observations and LES than any single-variable fitting. Three commonly used parameterizations are also tested on the individual cloud scale. A new parameterization is therefore presented that relates entrainment rate to vertical velocity, buoyancy and dissipation rate; the effects of treating clouds as ensembles and humid shells surrounding cumulus clouds on the new parameterization are discussed. Physical mechanisms underlying the relationships of entrainment rate to vertical velocity, buoyancy and dissipation rate are also explored.« less

Minimum-dissipation scalar transport model for large-eddy simulation of turbulent flows

NASA Astrophysics Data System (ADS)

Abkar, Mahdi; Bae, Hyun J.; Moin, Parviz

2016-08-01

Minimum-dissipation models are a simple alternative to the Smagorinsky-type approaches to parametrize the subfilter turbulent fluxes in large-eddy simulation. A recently derived model of this type for subfilter stress tensor is the anisotropic minimum-dissipation (AMD) model [Rozema et al., Phys. Fluids 27, 085107 (2015), 10.1063/1.4928700], which has many desirable properties. It is more cost effective than the dynamic Smagorinsky model, it appropriately switches off in laminar and transitional flows, and it is consistent with the exact subfilter stress tensor on both isotropic and anisotropic grids. In this study, an extension of this approach to modeling the subfilter scalar flux is proposed. The performance of the AMD model is tested in the simulation of a high-Reynolds-number rough-wall boundary-layer flow with a constant and uniform surface scalar flux. The simulation results obtained from the AMD model show good agreement with well-established empirical correlations and theoretical predictions of the resolved flow statistics. In particular, the AMD model is capable of accurately predicting the expected surface-layer similarity profiles and power spectra for both velocity and scalar concentration.

A comparison of observed and numerically predicted eddy kinetic energy budgets for a developing extratropical cyclone

NASA Technical Reports Server (NTRS)

Dare, P. M.; Smith, P. J.

1983-01-01

The eddy kinetic energy budget is calculated for a 48-hour forecast of an intense occluding winter cyclone associated with a strong well-developed jet stream. The model output consists of the initialized (1200 GMT January 9, 1975) and the 12, 24, 36, and 48 hour forecast fields from the Drexel/NCAR Limited Area Mesoscale Prediction System (LAMPS) model. The LAMPS forecast compares well with observations for the first 24 hours, but then overdevelops the low-level cyclone while inadequately developing the upper-air wave and jet. Eddy kinetic energy was found to be concentrated in the upper-troposphere with maxima flanking the primary trough. The increases in kinetic energy were found to be due to an excess of the primary source term of kinetic energy content, which is the horizontal flux of eddy kinetic energy over the primary sinks, and the generation and dissipation of eddy kinetic energy.

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.

Dynamical analysis of a satellite-observed anticyclonic eddy in the northern Bering Sea

NASA Astrophysics Data System (ADS)

Li, Yineng; Li, Xiaofeng; Wang, Jia; Peng, Shiqiu

2016-05-01

The characteristics and evolution of a satellite-observed anticyclonic eddy in the northern Bering Sea during March and April 1999 are investigated using a three-dimensional Princeton Ocean Model (POM). The anticyclonic-like current pattern and asymmetric feature of the eddy were clearly seen in the synthetic aperture radar (SAR), sea surface temperature, and ocean color images in April 1999. The results from model simulation reveal the three-dimensional structure of the anticyclonic eddy, its movement, and dissipation. Energy analysis indicates that the barotropic instability (BTI) is the main energy source for the growth of the anticyclonic eddy. The momentum analysis further reveals that the larger magnitude of the barotropic pressure gradient in the meridional direction causes the asymmetry of the anticyclonic eddy in the zonal and meridional directions, while the different magnitudes of the meridional baroclinic pressure gradient are responsible for the different intensity of currents between the northern and southern parts of the anticyclonic eddy. This article was corrected on 23 JUL 2016. See the end of the full text for details.

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

NASA Astrophysics Data System (ADS)

Novikov, E. A.

1990-05-01

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

Linear Multivariable Regression Models for Prediction of Eddy Dissipation Rate from Available Meteorological Data

NASA Technical Reports Server (NTRS)

MCKissick, Burnell T. (Technical Monitor); Plassman, Gerald E.; Mall, Gerald H.; Quagliano, John R.

2005-01-01

Linear multivariable regression models for predicting day and night Eddy Dissipation Rate (EDR) from available meteorological data sources are defined and validated. Model definition is based on a combination of 1997-2000 Dallas/Fort Worth (DFW) data sources, EDR from Aircraft Vortex Spacing System (AVOSS) deployment data, and regression variables primarily from corresponding Automated Surface Observation System (ASOS) data. Model validation is accomplished through EDR predictions on a similar combination of 1994-1995 Memphis (MEM) AVOSS and ASOS data. Model forms include an intercept plus a single term of fixed optimal power for each of these regression variables; 30-minute forward averaged mean and variance of near-surface wind speed and temperature, variance of wind direction, and a discrete cloud cover metric. Distinct day and night models, regressing on EDR and the natural log of EDR respectively, yield best performance and avoid model discontinuity over day/night data boundaries.

Computational fluid dynamic on the temperature simulation of air preheat effect combustion in propane turbulent flame

NASA Astrophysics Data System (ADS)

Elwina; Yunardi; Bindar, Yazid

2018-04-01

this paper presents results obtained from the application of a computational fluid dynamics (CFD) code Fluent 6.3 to modelling of temperature in propane flames with and without air preheat. The study focuses to investigate the effect of air preheat temperature on the temperature of the flame. A standard k-ε model and Eddy Dissipation model are utilized to represent the flow field and combustion of the flame being investigated, respectively. The results of calculations are compared with experimental data of propane flame taken from literature. The results of the study show that a combination of the standard k-ε turbulence model and eddy dissipation model is capable of producing reasonable predictions of temperature, particularly in axial profile of all three flames. Both experimental works and numerical simulation showed that increasing the temperature of the combustion air significantly increases the flame temperature.

Large Eddy simulation of compressible flows with a low-numerical dissipation patch-based adaptive mesh refinement method

NASA Astrophysics Data System (ADS)

Pantano, Carlos

2005-11-01

We describe a hybrid finite difference method for large-eddy simulation (LES) of compressible flows with a low-numerical dissipation scheme and structured adaptive mesh refinement (SAMR). Numerical experiments and validation calculations are presented including a turbulent jet and the strongly shock-driven mixing of a Richtmyer-Meshkov instability. The approach is a conservative flux-based SAMR formulation and as such, it utilizes refinement to computational advantage. The numerical method for the resolved scale terms encompasses the cases of scheme alternation and internal mesh interfaces resulting from SAMR. An explicit centered scheme that is consistent with a skew-symmetric finite difference formulation is used in turbulent flow regions while a weighted essentially non-oscillatory (WENO) scheme is employed to capture shocks. The subgrid stresses and transports are calculated by means of the streched-vortex model, Misra & Pullin (1997)

Time tracking and interaction of energy-eddies at different scales

NASA Astrophysics Data System (ADS)

Cardesa, Jose I.; Vela-Martin, Alberto; Jimenez, Javier

2016-11-01

We study the energy cascade through coherent structures obtained in time-resolved simulations of incompressible, statistically steady isotropic turbulence. The structures are defined as geometrically connected regions of the flow with high kinetic energy. We compute the latter by band-pass filtering the velocity field around a scale r. We analyse the dynamics of structures extracted with different r, which are a proxy for eddies containing energy at those r. We find that the size of these "energy-eddies" scales with r, while their lifetime scales with the local eddy-turnover r 2 / 3ɛ - 1 / 3 , where ɛ is the energy dissipation averaged over all space and time. Furthermore, a statistical analysis over the lives of the eddies shows a slight predominance of the splitting over the merging process. When we isolate the eddies which do not interact with other eddies of the same scale, we observe a parent-child dependence by which, on average, structures are born at scale r during the decaying part of the life of a structure at scale r' > r . The energy-eddy at r' lives in the same region of space as that at r. Finally, we investigate how interactions between eddies at the same scale are echoed across other scales. Funded by the ERC project Coturb.

Comparison of Forecast and Observed Energetics

NASA Technical Reports Server (NTRS)

Baker, W. E.; Brin, Y.

1985-01-01

An energetics analysis scheme was developed to compare the observed kinetic energy balance over North America with that derived from forecast cyclone case. It is found that: (1) the observed and predicted kinetic energy and eddy conversion are in good qualitative agreement, although the model eddy conversion tends to be 2 to 3 times stronger than the observed values. The eddy conversion which is stronger in the 12 h forecast than in observations and may be due to several factors is studied; (2) vertical profiles of kinetic energy generation and dissipation exhibit lower and upper tropospheric maxima in both the forecast and observations; and (3) a lag in the observational analysis with the maximum in the observed kinetic energy occurring at 0000 GMT 14 January over the same region as the maximum Eddy conversion 12 h earlier is noted.

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

A New Eddy Dissipation Rate Formulation for the Terminal Area PBL Prediction System(TAPPS)

NASA Technical Reports Server (NTRS)

Charney, Joseph J.; Kaplan, Michael L.; Lin, Yuh-Lang; Pfeiffer, Karl D.

2000-01-01

The TAPPS employs the MASS model to produce mesoscale atmospheric simulations in support of the Wake Vortex project at Dallas Fort-Worth International Airport (DFW). A post-processing scheme uses the simulated three-dimensional atmospheric characteristics in the planetary boundary layer (PBL) to calculate the turbulence quantities most important to the dissipation of vortices: turbulent kinetic energy and eddy dissipation rate. TAPPS will ultimately be employed to enhance terminal area productivity by providing weather forecasts for the Aircraft Vortex Spacing System (AVOSS). The post-processing scheme utilizes experimental data and similarity theory to determine the turbulence quantities from the simulated horizontal wind field and stability characteristics of the atmosphere. Characteristic PBL quantities important to these calculations are determined based on formulations from the Blackadar PBL parameterization, which is regularly employed in the MASS model to account for PBL processes in mesoscale simulations. The TAPPS forecasts are verified against high-resolution observations of the horizontal winds at DFW. Statistical assessments of the error in the wind forecasts suggest that TAPPS captures the essential features of the horizontal winds with considerable skill. Additionally, the turbulence quantities produced by the post-processor are shown to compare favorably with corresponding tower observations.

Assessment of dynamic closure for premixed combustion large eddy simulation

NASA Astrophysics Data System (ADS)

Langella, Ivan; Swaminathan, Nedunchezhian; Gao, Yuan; Chakraborty, Nilanjan

2015-09-01

Turbulent piloted Bunsen flames of stoichiometric methane-air mixtures are computed using the large eddy simulation (LES) paradigm involving an algebraic closure for the filtered reaction rate. This closure involves the filtered scalar dissipation rate of a reaction progress variable. The model for this dissipation rate involves a parameter βc representing the flame front curvature effects induced by turbulence, chemical reactions, molecular dissipation, and their interactions at the sub-grid level, suggesting that this parameter may vary with filter width or be a scale-dependent. Thus, it would be ideal to evaluate this parameter dynamically by LES. A procedure for this evaluation is discussed and assessed using direct numerical simulation (DNS) data and LES calculations. The probability density functions of βc obtained from the DNS and LES calculations are very similar when the turbulent Reynolds number is sufficiently large and when the filter width normalised by the laminar flame thermal thickness is larger than unity. Results obtained using a constant (static) value for this parameter are also used for comparative evaluation. Detailed discussion presented in this paper suggests that the dynamic procedure works well and physical insights and reasonings are provided to explain the observed behaviour.

Computational Analysis of Static and Dynamic Behaviour of Magnetic Suspensions and Magnetic Bearings

NASA Technical Reports Server (NTRS)

Britcher, Colin P. (Editor); Groom, Nelson J.

1996-01-01

Static modelling of magnetic bearings is often carried out using magnetic circuit theory. This theory cannot easily include nonlinear effects such as magnetic saturation or the fringing of flux in air-gaps. Modern computational tools are able to accurately model complex magnetic bearing geometries, provided some care is exercised. In magnetic suspension applications, the magnetic fields are highly three-dimensional and require computational tools for the solution of most problems of interest. The dynamics of a magnetic bearing or magnetic suspension system can be strongly affected by eddy currents. Eddy currents are present whenever a time-varying magnetic flux penetrates a conducting medium. The direction of flow of the eddy current is such as to reduce the rate-of-change of flux. Analytic solutions for eddy currents are available for some simplified geometries, but complex geometries must be solved by computation. It is only in recent years that such computations have been considered truly practical. At NASA Langley Research Center, state-of-the-art finite-element computer codes, 'OPERA', 'TOSCA' and 'ELEKTRA' have recently been installed and applied to the magnetostatic and eddy current problems. This paper reviews results of theoretical analyses which suggest general forms of mathematical models for eddy currents, together with computational results. A simplified circuit-based eddy current model proposed appears to predict the observed trends in the case of large eddy current circuits in conducting non-magnetic material. A much more difficult case is seen to be that of eddy currents in magnetic material, or in non-magnetic material at higher frequencies, due to the lower skin depths. Even here, the dissipative behavior has been shown to yield at least somewhat to linear modelling. Magnetostatic and eddy current computations have been carried out relating to the Annular Suspension and Pointing System, a prototype for a space payload pointing and vibration isolation system, where the magnetic actuator geometry resembles a conventional magnetic bearing. Magnetostatic computations provide estimates of flux density within airgaps and the iron core material, fringing at the pole faces and the net force generated. Eddy current computations provide coil inductance, power dissipation and the phase lag in the magnetic field, all as functions of excitation frequency. Here, the dynamics of the magnetic bearings, notably the rise time of forces with changing currents, are found to be very strongly affected by eddy currents, even at quite low frequencies. Results are also compared to experimental measurements of the performance of a large-gap magnetic suspension system, the Large Angle Magnetic Suspension Test Fixture (LAMSTF). Eddy current effects are again shown to significantly affect the dynamics of the system. Some consideration is given to the ease and accuracy of computation, specifically relating to OPERA/TOSCA/ELEKTRA.

Copepod behavior response to Burgers' vortex treatments mimicking turbulent eddies

NASA Astrophysics Data System (ADS)

Elmi, D.; Webster, D. R.; Fields, D. M.

2017-11-01

Copepods detect hydrodynamic cues in the water by their mechanosensory setae. We expect that copepods sense the flow structure of turbulent eddies in order to evoke behavioral responses that lead to population-scale distribution patterns. In this study, the copepods' response to the Burgers' vortex is examined. The Burgers' vortex is a steady-state solution of three-dimensional Navier-Stokes equations that allows us to mimic turbulent vortices at the appropriate scale and eliminate the stochastic nature of turbulence. We generate vortices in the laboratory oriented in the horizontal and vertical directions each with four intensity levels. The objective of including vortex orientation as a parameter in the study is to quantify directional responses that lead to vertical population distribution patterns. The four intensity levels correspond to target vortex characteristics of eddies corresponding to the typical dissipative vortices in isotropic turbulence with mean turbulent dissipation rates in the range of 0.002 to 0.25 cm2/s3. These vortices mimic the characteristics of eddies that copepods most likely encounter in coastal zones. We hypothesize that the response of copepods to hydrodynamic features depends on their sensory architecture and relative orientation with respect to gravity. Tomo-PIV is used to quantify the vortex circulation and axial strain rate for each vortex treatment. Three-dimensional trajectories of the copepod species Calanus finmarchicus are analyzed to examine their swimming kinematics in and around the vortex to quantify the hydrodynamic cues that trigger their behavior.

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.

A Comparison of Martian Transient Wave Energetics in High and Low Optical Depth Environments

NASA Astrophysics Data System (ADS)

Battalio, J. M.; Szunyogh, I.; Lemmon, M. T.

2016-12-01

The local energetics of individual transient eddies from the Mars Analysis Correction Data Assimilation (MACDA) is compared between a year with a global-scale dust storm (MY 25) and two years of relatively low optical depth conditions. Eddies in each year are considered from a period of strong wave activity in the northern hemisphere before the winter solstice (Ls=170-240°). The local growth of eddies is typically triggered by geopotential flux convergence. While all waves exhibit some baroclinic growth, baroclinic energy conversion is weaker in the waves that occur during the global-scale dust storm. The weaker baroclinic energy conversion in these waves, however, is compensated by a more intense barotropic transfer of the kinetic energy from the mean flow to the waves: the contribution from barotropic energy conversion allows eddies during the global-scale dust storm to attain roughly the same maximum eddy kinetic energy as eddies during the low optical depth years. Individual eddies in the waves decay through a combination of barotropic conversion of the kinetic energy from the waves to the mean flow, geopotential flux divergence, and dissipation in both the high- and the low-optical-depth years.

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.

Parametric study of power absorption from electromagnetic waves by small ferrite spheres

NASA Technical Reports Server (NTRS)

Englert, Gerald W.

1989-01-01

Algebraic expressions in terms of elementary mathematical functions are derived for power absorption and dissipation by eddy currents and magnetic hysteresis in ferrite spheres. Skin depth is determined by using a variable inner radius in descriptive integral equations. Numerical results are presented for sphere diameters less than one wavelength. A generalized power absorption parameter for both eddy currents and hysteresis is expressed in terms of the independent parameters involving wave frequency, sphere radius, resistivity, and complex permeability. In general, the hysteresis phenomenon has a greater sensitivity to these independent parameters than do eddy currents over the ranges of independent parameters studied herein. Working curves are presented for obtaining power losses from input to the independent parameters.

A numerical comparison of cold flow and combustion characteristics for GCH4/GO2 splash platelet injector

NASA Astrophysics Data System (ADS)

Yin, Liang; Liu, Weiqiang

2018-04-01

The differences between cold flow and combustion under the same condition were investigated by the numerical simulations, an eddy dissipation concept (EDC) with 16 species and 41 reactions is considered for the CH4/O2 combustion. Three configurations of the splash platelet injector were selected for these simulations. Results show that cold flow and combustion have evident differences. Compared with cold flow, CH4 mole fraction was more evenly distributed in the combustion chamber head, and the mixing of propellants was lagged by the combustion of multi-elements. However, this conclusion is contrary for the single element. The recirculation zones were observable near the injector faceplate at the combustion condition. Moreover, the cold flow simulation cannot reflect the actual combustion but can provide a reference value for experimental research.

Air motions accompanying the development of a planetary wave critical layer

NASA Technical Reports Server (NTRS)

Salby, Murry L.; O'Sullivan, Donal; Callaghan, Patrick; Garcia, Rolando R.

1990-01-01

The horizontal air motions accompanying the development of a planetary wave critical layer are presently investigated on the sphere, in terms of wave amplitude, the characteristics of the zonal flow, and dissipation. While attention is given to adiabatic motions, which should furnish an upper bound on the redistribution of conserved quantities by eddy stirring, nonconservative processes may be important in determining how large a role eddy stirring actually plays in the redistribution of atmospheric constituents. Nonconservative processes may also influence tracer distributions by directly affecting dynamics.

Large-eddy simulation of flow past a circular cylinder

NASA Technical Reports Server (NTRS)

Mittal, R.

1995-01-01

Some of the most challenging applications of large-eddy simulation are those in complex geometries where spectral methods are of limited use. For such applications more conventional methods such as finite difference or finite element have to be used. However, it has become clear in recent years that dissipative numerical schemes which are routinely used in viscous flow simulations are not good candidates for use in LES of turbulent flows. Except in cases where the flow is extremely well resolved, it has been found that upwind schemes tend to damp out a significant portion of the small scales that can be resolved on the grid. Furthermore, it has been found that even specially designed higher-order upwind schemes that have been used successfully in the direct numerical simulation of turbulent flows produce too much dissipation when used in conjunction with large-eddy simulation. The objective of the current study is to perform a LES of incompressible flow past a circular cylinder at a Reynolds number of 3900 using a solver which employs an energy-conservative second-order central difference scheme for spatial discretization and compare the results obtained with those of Beaudan & Moin (1994) and with the experiments in order to assess the performance of the central scheme for this relatively complex geometry.

Malvinas Current variations: dissipation of mesoscale activity over the Malvinas Plateau and recurrent blocking events in the Argentine Basin.

NASA Astrophysics Data System (ADS)

Provost, C.; Artana, C.; Ferrari, R.; Koenig, Z.; Saraceno, M.; Piola, A. R.

2016-12-01

The Malvinas Current (MC) is an offshoot of the Antarctic Circumpolar Current (ACC). Downstream of Drake Passage, the northern fronts of the ACC veer northward, cross over the North Scotia Ridge (NSR) and the Malvinas Plateau and enter the Argentine Basin. We investigate the variations of the MC circulation between the NSR and 41°S and their possible relations with the ACC circulation using data from Argo floats and satellite altimetry. The data depict meandering and eddy-shedding of the northern ACC jets as they cross the NSR. The satellite fields (altimetry and high resolution sea surface temperature images) show that these eddies are trapped, break down and dissipate over the Malvinas Plateau, suggesting that this region is a hot spot for dissipation of mesoscale variability. Variations of sea level anomalies (SLA) across the NSR do not impact the MC further north, except for intra-seasonal variability associated with coastal trapped waves. Altimetry and float trajectories show events during which a large fraction of the MC is cut off from the ACC. During these blocking events, the MC does not collapse as a robust cyclonic cell is established to the north of the cut-off. The MC becomes the western boundary current of the cell and small cyclonic eddies locally reinforce the circulation. Blocking events at around 48.5°S are a recurrent feature of the MC circulation. Over the 23 year altimetry record, we detected 26 events during which the MC surface transport at 48.5°S was reduced to less than half its long term mean. Blocking events last from 10 to 35 days and do not present any significant trend. These events were tracked back to positive SLA that built up over the Argentine Abyssal Plain.

Numerical modeling of separated flows at moderate Reynolds numbers appropriate for turbine blades and unmanned aero vehicles

NASA Astrophysics Data System (ADS)

Castiglioni, Giacomo

Flows over airfoils and blades in rotating machinery, for unmanned and micro-aerial vehicles, wind turbines, and propellers consist of a laminar boundary layer near the leading edge that is often followed by a laminar separation bubble and transition to turbulence further downstream. Typical Reynolds averaged Navier-Stokes turbulence models are inadequate for such flows. Direct numerical simulation is the most reliable, but is also the most computationally expensive alternative. This work assesses the capability of immersed boundary methods and large eddy simulations to reduce the computational requirements for such flows and still provide high quality results. Two-dimensional and three-dimensional simulations of a laminar separation bubble on a NACA-0012 airfoil at Rec = 5x104 and at 5° of incidence have been performed with an immersed boundary code and a commercial code using body fitted grids. Several sub-grid scale models have been implemented in both codes and their performance evaluated. For the two-dimensional simulations with the immersed boundary method the results show good agreement with the direct numerical simulation benchmark data for the pressure coefficient Cp and the friction coefficient Cf, but only when using dissipative numerical schemes. There is evidence that this behavior can be attributed to the ability of dissipative schemes to damp numerical noise coming from the immersed boundary. For the three-dimensional simulations the results show a good prediction of the separation point, but an inaccurate prediction of the reattachment point unless full direct numerical simulation resolution is used. The commercial code shows good agreement with the direct numerical simulation benchmark data in both two and three-dimensional simulations, but the presence of significant, unquantified numerical dissipation prevents a conclusive assessment of the actual prediction capabilities of very coarse large eddy simulations with low order schemes in general cases. Additionally, a two-dimensional sweep of angles of attack from 0° to 5° is performed showing a qualitative prediction of the jump in lift and drag coefficients due to the appearance of the laminar separation bubble. The numerical dissipation inhibits the predictive capabilities of large eddy simulations whenever it is of the same order of magnitude or larger than the sub-grid scale dissipation. The need to estimate the numerical dissipation is most pressing for low-order methods employed by commercial computational fluid dynamics codes. Following the recent work of Schranner et al., the equations and procedure for estimating the numerical dissipation rate and the numerical viscosity in a commercial code are presented. The method allows for the computation of the numerical dissipation rate and numerical viscosity in the physical space for arbitrary sub-domains in a self-consistent way, using only information provided by the code in question. The method is first tested for a three-dimensional Taylor-Green vortex flow in a simple cubic domain and compared with benchmark results obtained using an accurate, incompressible spectral solver. Afterwards the same procedure is applied for the first time to a realistic flow configuration, specifically to the above discussed laminar separation bubble flow over a NACA 0012 airfoil. The method appears to be quite robust and its application reveals that for the code and the flow in question the numerical dissipation can be significantly larger than the viscous dissipation or the dissipation of the classical Smagorinsky sub-grid scale model, confirming the previously qualitative finding.

A mixed multiscale model better accounting for the cross term of the subgrid-scale stress and for backscatter

NASA Astrophysics Data System (ADS)

Thiry, Olivier; Winckelmans, Grégoire

2016-02-01

In the large-eddy simulation (LES) of turbulent flows, models are used to account for the subgrid-scale (SGS) stress. We here consider LES with "truncation filtering only" (i.e., that due to the LES grid), thus without regular explicit filtering added. The SGS stress tensor is then composed of two terms: the cross term that accounts for interactions between resolved scales and unresolved scales, and the Reynolds term that accounts for interactions between unresolved scales. Both terms provide forward- (dissipation) and backward (production, also called backscatter) energy transfer. Purely dissipative, eddy-viscosity type, SGS models are widely used: Smagorinsky-type models, or more advanced multiscale-type models. Dynamic versions have also been developed, where the model coefficient is determined using a dynamic procedure. Being dissipative by nature, those models do not provide backscatter. Even when using the dynamic version with local averaging, one typically uses clipping to forbid negative values of the model coefficient and hence ensure the stability of the simulation; hence removing the backscatter produced by the dynamic procedure. More advanced SGS model are thus desirable, and that better conform to the physics of the true SGS stress, while remaining stable. We here investigate, in decaying homogeneous isotropic turbulence, and using a de-aliased pseudo-spectral method, the behavior of the cross term and of the Reynolds term: in terms of dissipation spectra, and in terms of probability density function (pdf) of dissipation in physical space: positive and negative (backscatter). We then develop a new mixed model that better accounts for the physics of the SGS stress and for the backscatter. It has a cross term part which is built using a scale-similarity argument, further combined with a correction for Galilean invariance using a pseudo-Leonard term: this is the term that also does backscatter. It also has an eddy-viscosity multiscale model part that accounts for all the remaining phenomena (also for the incompleteness of the cross term model), that is dynamic and that adjusts the overall dissipation. The model is tested, both a priori and a posteriori, and is compared to the direct numerical simulation and to the exact SGS terms, also in time. The model is seen to provide accurate energy spectra, also in comparison to the dynamic Smagorinsky model. It also provides significant backscatter (although four times less than the real SGS stress), while remaining stable.

A novel eddy current damper: theory and experiment

NASA Astrophysics Data System (ADS)

Ebrahimi, Babak; Khamesee, Mir Behrad; Golnaraghi, Farid

2009-04-01

A novel eddy current damper is developed and its damping characteristics are studied analytically and experimentally. The proposed eddy current damper consists of a conductor as an outer tube, and an array of axially magnetized ring-shaped permanent magnets separated by iron pole pieces as a mover. The relative movement of the magnets and the conductor causes the conductor to undergo motional eddy currents. Since the eddy currents produce a repulsive force that is proportional to the velocity of the conductor, the moving magnet and the conductor behave as a viscous damper. The eddy current generation causes the vibration to dissipate through the Joule heating generated in the conductor part. An accurate, analytical model of the system is obtained by applying electromagnetic theory to estimate the damping properties of the proposed eddy current damper. A prototype eddy current damper is fabricated, and experiments are carried out to verify the accuracy of the theoretical model. The experimental test bed consists of a one-degree-of-freedom vibration isolation system and is used for the frequency and transient time response analysis of the system. The eddy current damper model has a 0.1 m s-2 (4.8%) RMS error in the estimation of the mass acceleration. A damping coefficient as high as 53 Ns m-1 is achievable with the fabricated prototype. This novel eddy current damper is an oil-free, inexpensive damper that is applicable in various vibration isolation systems such as precision machinery, micro-mechanical suspension systems and structure vibration isolation.

On the properties of energy stable flux reconstruction schemes for implicit large eddy simulation

NASA Astrophysics Data System (ADS)

Vermeire, B. C.; Vincent, P. E.

2016-12-01

We begin by investigating the stability, order of accuracy, and dispersion and dissipation characteristics of the extended range of energy stable flux reconstruction (E-ESFR) schemes in the context of implicit large eddy simulation (ILES). We proceed to demonstrate that subsets of the E-ESFR schemes are more stable than collocation nodal discontinuous Galerkin methods recovered with the flux reconstruction approach (FRDG) for marginally-resolved ILES simulations of the Taylor-Green vortex. These schemes are shown to have reduced dissipation and dispersion errors relative to FRDG schemes of the same polynomial degree and, simultaneously, have increased Courant-Friedrichs-Lewy (CFL) limits. Finally, we simulate turbulent flow over an SD7003 aerofoil using two of the most stable E-ESFR schemes identified by the aforementioned Taylor-Green vortex experiments. Results demonstrate that subsets of E-ESFR schemes appear more stable than the commonly used FRDG method, have increased CFL limits, and are suitable for ILES of complex turbulent flows on unstructured grids.

Extension of a Kolmogorov Atmospheric Turbulence Model for Time-Based Simulation Implementation

NASA Technical Reports Server (NTRS)

McMinn, John D.

1997-01-01

The development of any super/hypersonic aircraft requires the interaction of a wide variety of technical disciplines to maximize vehicle performance. For flight and engine control system design and development on this class of vehicle, realistic mathematical simulation models of atmospheric turbulence, including winds and the varying thermodynamic properties of the atmosphere, are needed. A model which has been tentatively selected by a government/industry group of flight and engine/inlet controls representatives working on the High Speed Civil Transport is one based on the Kolmogorov spectrum function. This report compares the Dryden and Kolmogorov turbulence forms, and describes enhancements that add functionality to the selected Kolmogorov model. These added features are: an altitude variation of the eddy dissipation rate based on Dryden data, the mapping of the eddy dissipation rate database onto a regular latitude and longitude grid, a method to account for flight at large vehicle attitude angles, and a procedure for transitioning smoothly across turbulence segments.

DNS/LES Simulations of Separated Flows at High Reynolds Numbers

NASA Technical Reports Server (NTRS)

Balakumar, P.

2015-01-01

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

A Computational Fluid Dynamic Model for a Novel Flash Ironmaking Process

NASA Astrophysics Data System (ADS)

Perez-Fontes, Silvia E.; Sohn, Hong Yong; Olivas-Martinez, Miguel

A computational fluid dynamic model for a novel flash ironmaking process based on the direct gaseous reduction of iron oxide concentrates is presented. The model solves the three-dimensional governing equations including both gas-phase and gas-solid reaction kinetics. The turbulence-chemistry interaction in the gas-phase is modeled by the eddy dissipation concept incorporating chemical kinetics. The particle cloud model is used to track the particle phase in a Lagrangian framework. A nucleation and growth kinetics rate expression is adopted to calculate the reduction rate of magnetite concentrate particles. Benchmark experiments reported in the literature for a nonreacting swirling gas jet and a nonpremixed hydrogen jet flame were simulated for validation. The model predictions showed good agreement with measurements in terms of gas velocity, gas temperature and species concentrations. The relevance of the computational model for the analysis of a bench reactor operation and the design of an industrial-pilot plant is discussed.

Implicit LES using adaptive filtering

NASA Astrophysics Data System (ADS)

Sun, Guangrui; Domaradzki, Julian A.

2018-04-01

In implicit large eddy simulations (ILES) numerical dissipation prevents buildup of small scale energy in a manner similar to the explicit subgrid scale (SGS) models. If spectral methods are used the numerical dissipation is negligible but it can be introduced by applying a low-pass filter in the physical space, resulting in an effective ILES. In the present work we provide a comprehensive analysis of the numerical dissipation produced by different filtering operations in a turbulent channel flow simulated using a non-dissipative, pseudo-spectral Navier-Stokes solver. The amount of numerical dissipation imparted by filtering can be easily adjusted by changing how often a filter is applied. We show that when the additional numerical dissipation is close to the subgrid-scale (SGS) dissipation of an explicit LES the overall accuracy of ILES is also comparable, indicating that periodic filtering can replace explicit SGS models. A new method is proposed, which does not require any prior knowledge of a flow, to determine the filtering period adaptively. Once an optimal filtering period is found, the accuracy of ILES is significantly improved at low implementation complexity and computational cost. The method is general, performing well for different Reynolds numbers, grid resolutions, and filter shapes.

Nondestructive Testing Eddy Current Basic Principles RQA/M1-5330.12 (V-I).

ERIC Educational Resources Information Center

National Aeronautics and Space Administration, Huntsville, AL. George C. Marshall Space Flight Center.

As one in the series of programmed instruction handbooks, prepared by the U.S. space program, home study material is presented in this volume concerning familiarization and orientation on basic eddy current principles. The subject is presented under the following headings: Basic Eddy Current Concepts, Eddy Current Generation and Distribution,…

78 FR 69928 - 35th Meeting: RTCA Special Committee 206, Aeronautical Information and Meteorological Data Link...

Federal Register 2010, 2011, 2012, 2013, 2014

2013-11-21

...-Groups meetings SG4: SE2020 Eddy Dissipation Rate (EDR) Turbulence Project Plenary--SG3 Architecture Document FRAC Resolution 11 December Plenary--SG3 Architecture Document FRAC Resolution Sub-Group Meetings 12 December Plenary--SG3 Architecture Document FRAC Resolution Sub-Group Meetings 13 December Closing...

Effect of gravity waves on the North Atlantic circulation

NASA Astrophysics Data System (ADS)

Eden, Carsten

2017-04-01

The recently proposed IDEMIX (Internal wave Dissipation, Energy and MIXing) parameterisation for the effect of gravity waves offers the possibility to construct consistent ocean models with a closed energy cycle. This means that the energy available for interior mixing in the ocean is only controlled by external energy input from the atmosphere and the tidal system and by internal exchanges. A central difficulty is the unknown fate of meso-scale eddy energy. In different scenarios for that eddy dissipation, the parameterized internal wave field provides between 2 and 3 TW for interior mixing from the total external energy input of about 4 TW, such that a transfer between 0.3 and 0.4 TW into mean potential energy contributes to drive the large-scale circulation in the model. The impact of the different mixing on the meridional overturning in the North Atlantic is discussed and compared to hydrographic observations. Furthermore, the direct energy exchange of the wave field with the geostrophic flow is parameterized in extended IDEMIX versions and the sensitivity of the North Atlantic circulation by this gravity wave drag is discussed.

Large-eddy simulation of a spatially-evolving turbulent mixing layer

NASA Astrophysics Data System (ADS)

Capuano, Francesco; Catalano, Pietro; Mastellone, Andrea

2015-11-01

Large-eddy simulations of a spatially-evolving turbulent mixing layer have been performed. The flow conditions correspond to those of a documented experimental campaign (Delville, Appl. Sci. Res. 1994). The flow evolves downstream of a splitter plate separating two fully turbulent boundary layers, with Reθ = 2900 on the high-speed side and Reθ = 1200 on the low-speed side. The computational domain starts at the trailing edge of the splitter plate, where experimental mean velocity profiles are prescribed; white-noise perturbations are superimposed to mimic turbulent fluctuations. The fully compressible Navier-Stokes equations are solved by means of a finite-volume method implemented into the in-house code SPARK-LES. The results are mainly checked in terms of the streamwise evolution of the vorticity thickness and averaged velocity profiles. The combined effects of inflow perturbations, numerical accuracy and subgrid-scale model are discussed. It is found that excessive levels of dissipation may damp inlet fluctuations and delay the virtual origin of the turbulent mixing layer. On the other hand, non-dissipative, high-resolution computations provide results that are in much better agreement with experimental data.

High-resolution multi-code implementation of unsteady Navier-Stokes flow solver based on paralleled overset adaptive mesh refinement and high-order low-dissipation hybrid schemes

NASA Astrophysics Data System (ADS)

Li, Gaohua; Fu, Xiang; Wang, Fuxin

2017-10-01

The low-dissipation high-order accurate hybrid up-winding/central scheme based on fifth-order weighted essentially non-oscillatory (WENO) and sixth-order central schemes, along with the Spalart-Allmaras (SA)-based delayed detached eddy simulation (DDES) turbulence model, and the flow feature-based adaptive mesh refinement (AMR), are implemented into a dual-mesh overset grid infrastructure with parallel computing capabilities, for the purpose of simulating vortex-dominated unsteady detached wake flows with high spatial resolutions. The overset grid assembly (OGA) process based on collection detection theory and implicit hole-cutting algorithm achieves an automatic coupling for the near-body and off-body solvers, and the error-and-try method is used for obtaining a globally balanced load distribution among the composed multiple codes. The results of flows over high Reynolds cylinder and two-bladed helicopter rotor show that the combination of high-order hybrid scheme, advanced turbulence model, and overset adaptive mesh refinement can effectively enhance the spatial resolution for the simulation of turbulent wake eddies.

Large Eddy Simulation study of the development of finite-channel lock-release currents at high Grashof numbers

NASA Astrophysics Data System (ADS)

Ooi, Seng-Keat

2005-11-01

Lock-exchange gravity current flows produced by the instantaneous release of a heavy fluid are investigated using 3-D well resolved Large Eddy Simulation simulations at Grashof numbers up to 8*10^9. It is found the 3-D simulations correctly predict a constant front velocity over the initial slumping phase and a front speed decrease proportional to t-1/3 (the time t is measured from the release) over the inviscid phase, in agreement with theory. The evolution of the current in the simulations is found to be similar to that observed experimentally by Hacker et al. (1996). The effect of the dynamic LES model on the solutions is discussed. The energy budget of the current is discussed and the contribution of the turbulent dissipation to the total dissipation is analyzed. The limitations of less expensive 2D simulations are discussed; in particular their failure to correctly predict the spatio-temporal distributions of the bed shear stresses which is important in determining the amount of sediment the gravity current can entrain in the case in advances of a loose bed.

Turbulence: does vorticity affect the structure and shape of body and fin propulsors?

PubMed

Webb, P W; Cotel, A J

2010-12-01

Over the past century, many ideas have been developed on the relationships between water flow and the structure and shape of the body and fins of fishes, largely during swimming in relatively steady flows. However, both swimming by fishes and the habitats they occupy are associated with vorticity, typically concentrated as eddies characteristic of turbulent flow. Deployment of methods to examine flow in detail suggests that vorticity impacts the lives of fishes. First, vorticity near the body and fins can increase thrust and smooth variations in thrust that are a consequence of using oscillating and undulating propulsors to swim. Second, substantial mechanical energy is dissipated in eddies in the wake and adaptations that minimize these losses would be anticipated. We suggest that such mechanisms may be found in varying the length of the propulsive wave, stiffening propulsive surfaces, and shifting to using median and paired fins when swimming at low speeds. Eddies in the flow encountered by fishes may be beneficial, but when eddy radii are of the order of 0.25 of the fish's total length, negative impacts occur due to greater difficulties in controlling stability. The archetypal streamlined "fish" shape reduces destabilizing forces for fishes swimming into eddies.

Integrated and spectral energy flows of the GLAS GCM

NASA Technical Reports Server (NTRS)

Tennebaum, J.

1981-01-01

Methods to analyze the generation, transport, and dissipation of energy to study geophysical fluid flows are discussed. Energetics analyses are pursued in several directions: (1) the longitudinal and time dependence on the energy flow to the stratosphere was examined as a function of geographical sector; (2) strong and weak energy flows were correlated by medium range forecasts; (3) the one dimensional spectral results (Fourier services around latitude circles) were extended to spherical harmonics over a global domain; (4) the validity of vertical velocities derived from mass convergence was examined for their effect on the conversion of eddy available potential energy to eddy kinetic energy.

Dissipation in the Baltic proper during winter stratification

NASA Astrophysics Data System (ADS)

Lass, Hans Ulrich; Prandke, Hartmut; Liljebladh, Bengt

2003-06-01

Profiles of dissipation rates and stratification between 10 and 120 m depth were measured with a loosely tethered profiler over a 9-day winter period in the Gotland Basin of the Baltic Sea. Supplementary measurements of current profiles were made with moored ADCPs. Temporal and spatial patterns of the stratification were observed by means of towed CTD. Shallow freshwater lenses in the surface mixed layer, mesoscale eddies, inertial oscillations, and inertial waves as part of the internal wave spectrum provided the marine physical environment for the small-scale turbulence. Two well-separated turbulence regimes were detected. The turbulence in the surface mixed layer was well correlated with the wind. The majority of the energy flux from the wind to the turbulent kinetic energy was dissipated within the surface mixed layer. A minor part of this flux was consumed by changes of the potential energy of the fresh water lenses. The penetration depth Hpen of the wind-driven turbulence into the weakly stratified surface mixed layer depended on the local wind speed (W10) as Hpen = cW103/2 Active erosion of the Baltic halocline by wind-driven turbulence is expected for wind speeds greater than 14 m/s. The turbulence in the strongly stratified interior of the water column was quite independent of the meteorological forcing at the sea surface. The integrated production of turbulent kinetic energy exceeded the energy loss of inertial oscillations in the surface layer suggesting additional energy sources which might have been provided by inertial wave radiation during geostrophic adjustment of coastal jets and mesoscale eddies. The averaged dissipation rate profile in the stratified part of the water column, best fitted by ɛ ∝ EN, was different from the scaling of the dissipation in the thermocline of the ocean [, 1986]. The diapycnical mixing coefficient (Kv) was best fit by Kv = a0/N according to [1987] with a0 ≈ 0.87 × 10-7 m2/s2. The diapycnal diffusivity estimated from the dissipation rate was lower than those estimated by the bulk method.

Thermal Analysis of the PediaFlow pediatric ventricular assist device.

PubMed

Gardiner, Jeffrey M; Wu, Jingchun; Noh, Myounggyu D; Antaki, James F; Snyder, Trevor A; Paden, David B; Paden, Brad E

2007-01-01

Accurate modeling of heat dissipation in pediatric intracorporeal devices is crucial in avoiding tissue and blood thermotrauma. Thermal models of new Maglev ventricular assist device (VAD) concepts for the PediaFlow VAD are developed by incorporating empirical heat transfer equations with thermal finite element analysis (FEA). The models assume three main sources of waste heat generation: copper motor windings, active magnetic thrust bearing windings, and eddy currents generated within the titanium housing due to the two-pole motor. Waste heat leaves the pump by convection into blood passing through the pump and conduction through surrounding tissue. Coefficients of convection are calculated and assigned locally along fluid path surfaces of the three-dimensional pump housing model. FEA thermal analysis yields a three-dimensional temperature distribution for each of the three candidate pump models. Thermal impedances from the motor and thrust bearing windings to tissue and blood contacting surfaces are estimated based on maximum temperature rise at respective surfaces. A new updated model for the chosen pump topology is created incorporating computational fluid dynamics with empirical fluid and heat transfer equations. This model represents the final geometry of the first generation prototype, incorporates eddy current heating, and has 60 discrete convection regions. Thermal analysis is performed at nominal and maximum flow rates, and temperature distributions are plotted. Results suggest that the pump will not exceed a temperature rise of 2 degrees C during normal operation.

Improved turbulence models based on large eddy simulation of homogeneous, incompressible turbulent flows

NASA Technical Reports Server (NTRS)

Bardino, J.; Ferziger, J. H.; Reynolds, W. C.

1983-01-01

The physical bases of large eddy simulation and subgrid modeling are studied. A subgrid scale similarity model is developed that can account for system rotation. Large eddy simulations of homogeneous shear flows with system rotation were carried out. Apparently contradictory experimental results were explained. The main effect of rotation is to increase the transverse length scales in the rotation direction, and thereby decrease the rates of dissipation. Experimental results are shown to be affected by conditions at the turbulence producing grid, which make the initial states a function of the rotation rate. A two equation model is proposed that accounts for effects of rotation and shows good agreement with experimental results. In addition, a Reynolds stress model is developed that represents the turbulence structure of homogeneous shear flows very well and can account also for the effects of system rotation.

Reduction of magneto rheological dampers stiffness by incorporating of an eddy current damper

NASA Astrophysics Data System (ADS)

Asghar Maddah, Ali; Hojjat, Yousef; Reza Karafi, Mohammad; Reza Ashory, Mohammad

2017-05-01

In this paper, a hybrid damper is developed to achieve lower stiffness compared to magneto rheological dampers. The hybrid damper consists of an eddy current damper (ECD) and a Magneto Rheological Damper (MRD). The aim of this research is to reduce the stiffness of MRDs with equal damping forces. This work is done by adding an eddy current passive damper to a semi-active MRD. The ECDs are contactless dampers which show an almost viscous damping behavior without increasing the stiffness of a system. However, MRDs increase damping and stiffness of a system simultaneously, when a magnetic field is applied. Damping of each part is studied theoretically and experimentally. A semi-empirical model is developed to explain the viscoelastic behavior of the damper. The experimental results showed that the hybrid damper is able to dissipate energy as much as those of MRDs while its stiffness is 12% lower at a zero excitation current.

Comparison of Forecast and Observed Energetics

NASA Technical Reports Server (NTRS)

Baker, W. E.; Brin, Y.

1984-01-01

An energetics analysis scheme was developed to compare the observed kinetic energy balance over North America with that derived from forecast fields of the GLAS fourth order model for the 13 to 15 January 1979 cyclone case. It is found that: (1) the observed and predicted kinetic energy and eddy conversion are in good qualitative agreement, although the model eddy conversion tends to be 2 to 3 times stronger than the observed values. The eddy conversion which is stronger in the 12 h forecast than in observations and may be due to several factors is studied; (2) vertical profiles of kinetic energy generation and dissipation exhibit lower and upper tropospheric maxima in both the forecast and observations; (3) a lag in the observational analysis with the maximum in the observed kinetic energy occurring at 0000 GMT 14 January over the same region as the maximum ddy conversion 12 h earlier is noted.

A comparison of observed and forecast energetics over North America

NASA Technical Reports Server (NTRS)

Baker, W. E.; Brin, Y.

1985-01-01

The observed kinetic energy balance is calculated over North America and compared with that computed from forecast fields for the 13-15 January 1979 cyclone. The FGGE upper-air rawinsonde network serves as the observational database while the forecast energetics are derived from a numerical integration with the GLAS fourth-order general circulation model initialized at 00 GMT 13 January. Maps of the observed and predicted kinetic energy and eddy conversion are in good qualitative agreement, although the model eddy conversion tends to be 2 to 3 times stronger than the observed values. Both the forecast and observations exhibit the lower and upper tropospheric maxima in vertical profiles of kinetic energy generation and dissipation typically found in cyclonic disturbances. An interesting time lag is noted in the observational analysis with the maximum observed kinetic energy occurring 12 h later than the maximum eddy conversion over the same region.

The effect of the polymer relaxation time on the nonlinear energy cas- cade and dissipation of statistically steady and decaying homogeneous isotropic turbulence

NASA Astrophysics Data System (ADS)

Valente, Pedro C.; da Silva, Carlos B.; Pinho, Fernando T.

2013-11-01

We report a numerical study of statistically steady and decaying turbulence of FENE-P fluids for varying polymer relaxation times ranging from the Kolmogorov dissipation time-scale to the eddy turnover time. The total turbulent kinetic energy dissipation is shown to increase with the polymer relaxation time in both steady and decaying turbulence, implying a ``drag increase.'' If the total power input in the statistically steady case is kept equal in the Newtonian and the viscoelastic simulations the increase in the turbulence-polymer energy transfer naturally lead to the previously reported depletion of the Newtonian, but not the overall, kinetic energy dissipation. The modifications to the nonlinear energy cascade with varying Deborah/Weissenberg numbers are quantified and their origins investigated. The authors acknowledge the financial support from Fundação para a Ciência e a Tecnologia under grant PTDC/EME-MFE/113589/2009.

Global transport calculations with an equivalent barotropic system

NASA Technical Reports Server (NTRS)

Salby, Murry L.; O'Sullivan, Donal; Garcia, Rolando R.; Tribbia, Joseph

1990-01-01

Transport properties of the two-dimensional equations governing equivalent barotropic motion are investigated on the sphere. This system has ingredients such as forcing, equivalent depth, and thermal dissipation explicitly represented, and takes into account compression effects associated with vertical motion along isentropic surfaces. Horizontal transport properties of this system are investigated under adiabatic and diabatic conditions for different forms of dissipation, and over a range of resolutions. It is shown that forcing represetative of time-mean and amplified conditions at 10 mb leads to the behavior typical of observations at this level. The displacement of the polar night vortex and its distortion into a comma shape are evident, as is irreversible mixing under sufficiently strong forcing amplitude. It is shown that thermal dissipation influences the behavior significantly by inhibiting the amplification of unstable eddies and thereby the horizontal stirring of air.

From mesoscale eddies to small-scale turbulence in the Antarctic Circumpolar Current

NASA Astrophysics Data System (ADS)

Naveira Garabato, A.; Brearley, J. A.; Sheen, K. L.; Waterman, S. N.

2012-12-01

A foremost question in physical oceanography is that of how the oceanic mesoscale dissipates. The Antarctic Circumpolar Current (ACC), in the Southern Ocean, is forced strongly by the wind and hosts a vigorous mesoscale eddy field. It has been recently suggested that substantial dampening of mesoscale flows in the region may occur through interactions with topography, on the basis of a number of indirect approaches. Here, we present the first direct evidence of a transfer of energy between mesoscale eddies and small-scale turbulence in the ACC, via the radiation, instability and breaking of internal waves generated as mesoscale flows impinge on rough topography. The evidence is provided by analysis of two data sets gathered by the DIMES (Diapycnal and Isopycnal Experiment in the Southern Ocean) experiment: (1) the observations of a mooring cluster, specifically designed to measure dynamical exchanges between the mesoscale eddy and internal wave fields in Drake Passage over a 2-year deployment; and (2) an extensive fine- and microstructure survey of the region. The physical mechanisms implicated in the cascade of energy across scales will be discussed.

Turbulence Characteristics of Swirling Flowfields. Ph.D. Thesis

NASA Technical Reports Server (NTRS)

Jackson, T. W.

1983-01-01

Combustor design phenomena; recirculating flows research; single-wire, six-orientation, eddy dissipation rate, and turbulence modeling measurement; directional sensitivity (DS); calibration equipment, confined jet facility, and hot-wire instrumentation; effects of swirl, strong contraction nozzle, and expansion ratio; and turbulence parameters; uncertain; and DS in laminar jets; turbulent nonswirling jets, and turbulent swirling jets are discussed.

Cascades and dissipation ratio in rotating magnetohydrodynamic turbulence at low magnetic Prandtl number.

PubMed

Plunian, Franck; Stepanov, Rodion

2010-10-01

A phenomenology of isotropic magnetohydrodynamic (MHD) turbulence subject to both rotation and applied magnetic field is presented. It is assumed that the triple correlation decay time is the shortest between the eddy turn-over time and the ones associated to the rotating frequency and the Alfvén wave period. For Pm=1 it leads to four kinds of piecewise spectra, depending on four parameters: injection rate of energy, magnetic diffusivity, rotation rate, and applied field. With a shell model of MHD turbulence (including rotation and applied magnetic field), spectra for Pm ≤ 1 are presented, together with the ratio between magnetic and viscous dissipations.

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.

The Death of Two Eddies, Against the Shelf

NASA Astrophysics Data System (ADS)

Zavala-Trujillo, B.; Badan, A.; Rivas, D.; Ochoa, J.; Sheinbaum, J.; Candela, J.

2007-05-01

A set of five moorings deployed in front of the coast of Tamaulipas, western Gulf of Mexico, provided fourteen months (from August 2004 to November 2005) of surface to bottom observations of currents and temperature that document the processes associated with the collision and dissipation of two warm mesoscale eddies with the continental slope. Two Loop Current eddies (Titanic and Ulysses) were identified reaching the study area during the observation period. On September 2004, the two southernmost 2000-m moorings show that temperature and salinity increases throughout the entire water column, related to eddy Titanic; similarily; on April 2005, eddy Ulysses caused a strong increase of temperature in the 3500-m mooring. The velocity field suggests three different régimes: a coastal region, the continental slope currents, and the abyssal circulation. Over the slope, three different layers can be identified: a surface layer (above 500 m depth), influenced by eddies and transients, a deep layer (under de 1900 m) with a persistent southerly current and a transition layer (from 500 to 1900 m) that separates them. The variance ellipses at ~ 700 m at the 3500-m mooring have no a predominant orientation of the mayor axis. At the northernmost 2000-m mooring, the axis of maximum variation is oriented with the bathymetry, but at the southernmost 2000-m mooring it is perpendicular to the coast. The spectral characteristics of the measurements are also discussed.

Eddy-Current Reference Standard

NASA Technical Reports Server (NTRS)

Ambrose, H. H., Jr.

1985-01-01

Magnetic properties of metallic reference standards duplicated and stabilized for eddy-current coil measurements over long times. Concept uses precisely machined notched samples of known annealed materials as reference standards.

On the spatial distribution of the transpiration and soil moisture of a Mediterranean heterogeneous ecosystem in water-limited conditions.

NASA Astrophysics Data System (ADS)

Curreli, Matteo; Corona, Roberto; Montaldo, Nicola; Albertson, John D.; Oren, Ram

2014-05-01

Mediterranean ecosystems are characterized by a strong heterogeneity, and often by water-limited conditions. In these conditions contrasting plant functional types (PFT, e.g. grass and woody vegetation) compete for the water use. Both the vegetation cover spatial distribution and the soil properties impact the soil moisture (SM) spatial distribution. Indeed, vegetation cover density and type affects evapotranspiration (ET), which is the main lack of the soil water balance in these ecosystems. With the objective to carefully estimate SM and ET spatial distribution in a Mediterranean water-limited ecosystem and understanding SM and ET relationships, an extended field campaign is carried out. The study was performed in a heterogeneous ecosystem in Orroli, Sardinia (Italy). The experimental site is a typical Mediterranean ecosystem where the vegetation is distributed in patches of woody vegetation (wild olives mainly) and grass. Soil depth is low and spatially varies between 10 cm and 40 cm, without any correlation with the vegetation spatial distribution. ET, land-surface fluxes and CO2 fluxes are estimated by an eddy covariance technique based micrometeorological tower. But in heterogeneous ecosystems a key assumption of the eddy covariance theory, the homogeneity of the surface, is not preserved and the ET estimate may be not correct. Hence, we estimate ET of the woody vegetation using the thermal dissipation method (i.e. sap flow technique) for comparing the two methodologies. Due the high heterogeneity of the vegetation and soil properties of the field a total of 54 sap flux sensors were installed. 14 clumps of wild olives within the eddy covariance footprint were identified as the most representative source of flux and they were instrumented with the thermal dissipation probes. Measurements of diameter at the height of sensor installation (height of 0.4 m above ground) were recorded in all the clumps. Bark thickness and sapwood depth were measured on several trees to obtain a generalized estimates of sapwood depth. The known of allometric relationships between sapwood area, diameter and canopy cover area within the eddy covariance footprint helped for the application of a reliable scaling procedure of the local sap flow estimates which are in a good agreement with the estimates of ET eddy covariance based. Soil moisture were also extensively monitored through 25 probes installed in the eddy covariance footprint. Results show that comparing eddy covariance and sap flow ET estimates eddy covariance technique is still accurate in this heterogeneous field, whereas the key assumption, surface homogeneity, is not preserved. Furthermore, interestingly wild olives still transpire at higher rates for the driest soil moisture conditions, confirming the hydraulic redistribution from soil below the roots, and from roots penetrating deep cracks in the underlying basalt parent rock.

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

Retrieving atmospheric turbulence information from regular commercial aircraft using Mode-S and ADS-B

NASA Astrophysics Data System (ADS)

Kopeć, Jacek M.; Kwiatkowski, Kamil; de Haan, Siebren; Malinowski, Szymon P.

2016-05-01

Navigational information broadcast by commercial aircraft in the form of Mode-S EHS (Mode-S Enhanced Surveillance) and ADS-B (Automatic Dependent Surveillance-Broadcast) messages can be considered a new source of upper tropospheric and lower stratospheric turbulence estimates. A set of three processing methods is proposed and analysed using a quality record of turbulence encounters made by a research aircraft.The proposed methods are based on processing the vertical acceleration or the background wind into the eddy dissipation rate. Turbulence intensity can be estimated using the standard content of the Mode-S EHS/ADS-B.The results are based on a Mode-S EHS/ADS-B data set generated synthetically based on the transmissions from the research aircraft. This data set was validated using the overlapping record of the Mode-S EHS/ADS-B received from the same research aircraft. The turbulence intensity, meaning the eddy dissipation rate, obtained from the proposed methods based on the Mode-S EHS/ADS-B is compared with the value obtained using on-board accelerometer. The results of the comparison indicate the potential of the methods. The advantages and limitation of the presented approaches are discussed.

Capturing remote mixing due to internal tides using multi-scale modeling tool: SOMAR-LES

NASA Astrophysics Data System (ADS)

Santilli, Edward; Chalamalla, Vamsi; Scotti, Alberto; Sarkar, Sutanu

2016-11-01

Internal tides that are generated during the interaction of an oscillating barotropic tide with the bottom bathymetry dissipate only a fraction of their energy near the generation region. The rest is radiated away in the form of low- high-mode internal tides. These internal tides dissipate energy at remote locations when they interact with the upper ocean pycnocline, continental slope, and large scale eddies. Capturing the wide range of length and time scales involved during the life-cycle of internal tides is computationally very expensive. A recently developed multi-scale modeling tool called SOMAR-LES combines the adaptive grid refinement features of SOMAR with the turbulence modeling features of a Large Eddy Simulation (LES) to capture multi-scale processes at a reduced computational cost. Numerical simulations of internal tide generation at idealized bottom bathymetries are performed to demonstrate this multi-scale modeling technique. Although each of the remote mixing phenomena have been considered independently in previous studies, this work aims to capture remote mixing processes during the life cycle of an internal tide in more realistic settings, by allowing multi-level (coarse and fine) grids to co-exist and exchange information during the time stepping process.

Magnetic Reconnection in MHD and Kinetic Turbulence

NASA Astrophysics Data System (ADS)

Loureiro, Nuno; Boldyrev, Stanislav

2017-10-01

Recent works have revisited the current understanding of Alfvénic turbulence to account for the role of magnetic reconnection. Theoretical arguments suggest that reconnection inevitably becomes important in the inertial range, at the scale where it becomes faster than the eddy turnover time. This leads to a transition to a new sub-inertial interval, suggesting a route to energy dissipation that is fundamentally different from that envisioned in the usual Kolmogorov-like phenomenology. These concepts can be extended to collisionless plasmas, where reconnection is enabled by electron inertia rather than resistivity. Although several different cases must then be considered, a common result is that the energy spectrum exhibits a scaling with the perpendicular wave number that scales between k⊥- 8 / 3 and k⊥- 3 , in favourable agreement with many numerical results and observations. Work supported by NSF-DOE Partnership in Basic Plasma Science and Engineering, Award No. DE-SC0016215, and by NSF CAREER Award No. 1654168 (NFL); and by NSF Grant NSF AGS- 1261659 and by the Vilas Associates Award of UWM (SB).

Progress in turbulence modeling for complex flow fields including effects of compressibility

NASA Technical Reports Server (NTRS)

Wilcox, D. C.; Rubesin, M. W.

1980-01-01

Two second-order-closure turbulence models were devised that are suitable for predicting properties of complex turbulent flow fields in both incompressible and compressible fluids. One model is of the "two-equation" variety in which closure is accomplished by introducing an eddy viscosity which depends on both a turbulent mixing energy and a dissipation rate per unit energy, that is, a specific dissipation rate. The other model is a "Reynolds stress equation" (RSE) formulation in which all components of the Reynolds stress tensor and turbulent heat-flux vector are computed directly and are scaled by the specific dissipation rate. Computations based on these models are compared with measurements for the following flow fields: (a) low speed, high Reynolds number channel flows with plane strain or uniform shear; (b) equilibrium turbulent boundary layers with and without pressure gradients or effects of compressibility; and (c) flow over a convex surface with and without a pressure gradient.

ATMOSPHERIC CIRCULATION OF BROWN DWARFS: JETS, VORTICES, AND TIME VARIABILITY

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

Zhang, Xi; Showman, Adam P., E-mail: xiz@lpl.arizona.edu

2014-06-10

A variety of observational evidence demonstrates that brown dwarfs exhibit active atmospheric circulations. In this study we use a shallow-water model to investigate the global atmospheric dynamics in the stratified layer overlying the convective zone on these rapidly rotating objects. We show that the existence and properties of the atmospheric circulation crucially depend on key parameters including the energy injection rate and radiative timescale. Under conditions of strong internal heat flux and weak radiative dissipation, a banded flow pattern comprised of east-west jet streams spontaneously emerges from the interaction of atmospheric turbulence with the planetary rotation. In contrast, when themore » internal heat flux is weak and/or radiative dissipation is strong, turbulence injected into the atmosphere damps before it can self-organize into jets, leading to a flow dominated by transient eddies and isotropic turbulence instead. The simulation results are not very sensitive to the form of the forcing. Based on the location of the transition between jet-dominated and eddy-dominated regimes, we suggest that many brown dwarfs may exhibit atmospheric circulations dominated by eddies and turbulence (rather than jets) due to the strong radiative damping on these worlds, but a jet structure is also possible under some realistic conditions. Our simulated light curves capture important features from observed infrared light curves of brown dwarfs, including amplitude variations of a few percent and shapes that fluctuate between single-peak and multi-peak structures. More broadly, our work shows that the shallow-water system provides a useful tool to illuminate fundamental aspects of the dynamics on these worlds.« less

The equilibrium tide in stars and giant planets. I. The coplanar case

NASA Astrophysics Data System (ADS)

Remus, F.; Mathis, S.; Zahn, J.-P.

2012-08-01

Context. Since 1995, more than 500 extrasolar planets have been discovered orbiting very close to their parent star, where they experience strong tidal interactions. Their orbital evolution depends on the physical mechanisms that cause tidal dissipation, which remain poorly understood. Aims: We refine the theory of the equilibrium tide in fluid bodies that are partly or entirely convective, to predict the dynamical evolution of the systems. In particular, we examine the validity of modeling the tidal dissipation using the quality factor Q, which is commonly done. We consider here the simplest case where the considered star or planet rotates uniformly, all spins are aligned, and the companion is reduced to a point mass. Methods: We expand the tidal potential as a Fourier series, and express the hydrodynamical equations in the reference frame, which rotates with the corresponding Fourier component. The results are cast in the form of a complex disturbing function, which may be implemented directly in the equations governing the dynamical evolution of the system. Results: The first manifestation of the tide is to distort the shape of the star or planet adiabatically along the line of centers. This generates the divergence-free velocity field of the adiabatic equilibrium tide, which is stationary in the frame rotating with the considered Fourier component of the tidal potential; this large-scale velocity field is decoupled from the dynamical tide. The tidal kinetic energy is dissipated into heat by means of turbulent friction, which is modeled here as an eddy-viscosity acting on the adiabatic tidal flow. This dissipation induces a second velocity field, the dissipative equilibrium tide, which is in quadrature with the exciting potential; this field is responsible for the imaginary part of the disturbing function, which is implemented in the dynamical evolution equations, from which one derives the characteristic evolutionary times. Conclusions: The rate at which the system evolves depends on the physical properties of the tidal dissipation, and specifically on both how the eddy viscosity varies with tidal frequency and the thickness of the convective envelope for the fluid equilibrium tide. At low frequency, this tide is retarded by a constant time delay, whereas it lags behind by a constant angle when the tidal frequency exceeds the convective turnover rate.

Wind turbine wakes in forest and neutral plane wall boundary layer large-eddy simulations

NASA Astrophysics Data System (ADS)

Schröttle, Josef; Piotrowski, Zbigniew; Gerz, Thomas; Englberger, Antonia; Dörnbrack, Andreas

2016-09-01

Wind turbine wake flow characteristics are studied in a strongly sheared and turbulent forest boundary layer and a neutral plane wall boundary layer flow. The reference simulations without wind turbine yield similar results as earlier large-eddy simulations by Shaw and Schumann (1992) and Porte-Agel et al. (2000). To use the fields from the homogeneous turbulent boundary layers on the fly as inflow fields for the wind turbine wake simulations, a new and efficient methodology was developed for the multiscale geophysical flow solver EULAG. With this method fully developed turbulent flow fields can be achieved upstream of the wind turbine which are independent of the wake flow. The large-eddy simulations reproduce known boundary-layer statistics as mean wind profile, momentum flux profile, and eddy dissipation rate of the plane wall and the forest boundary layer. The wake velocity deficit is more asymmetric above the forest and recovers faster downstream compared to the velocity deficit in the plane wall boundary layer. This is due to the inflection point in the mean streamwise velocity profile with corresponding turbulent coherent structures of high turbulence intensity in the strong shear flow above the forest.

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.

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

Parametric Study of Decay of Homogeneous Isotropic Turbulence Using Large Eddy Simulation

NASA Technical Reports Server (NTRS)

Swanson, R. C.; Rumsey, Christopher L.; Rubinstein, Robert; Balakumar, Ponnampalam; Zang, Thomas A.

2012-01-01

Numerical simulations of decaying homogeneous isotropic turbulence are performed with both low-order and high-order spatial discretization schemes. The turbulent Mach and Reynolds numbers for the simulations are 0.2 and 250, respectively. For the low-order schemes we use either second-order central or third-order upwind biased differencing. For higher order approximations we apply weighted essentially non-oscillatory (WENO) schemes, both with linear and nonlinear weights. There are two objectives in this preliminary effort to investigate possible schemes for large eddy simulation (LES). One is to explore the capability of a widely used low-order computational fluid dynamics (CFD) code to perform LES computations. The other is to determine the effect of higher order accuracy (fifth, seventh, and ninth order) achieved with high-order upwind biased WENO-based schemes. Turbulence statistics, such as kinetic energy, dissipation, and skewness, along with the energy spectra from simulations of the decaying turbulence problem are used to assess and compare the various numerical schemes. In addition, results from the best performing schemes are compared with those from a spectral scheme. The effects of grid density, ranging from 32 cubed to 192 cubed, on the computations are also examined. The fifth-order WENO-based scheme is found to be too dissipative, especially on the coarser grids. However, with the seventh-order and ninth-order WENO-based schemes we observe a significant improvement in accuracy relative to the lower order LES schemes, as revealed by the computed peak in the energy dissipation and by the energy spectrum.

Dynamics, Stability, and Evolutionary Patterns of Mesoscale Intrathermocline Vortices

DTIC Science & Technology

2016-12-01

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

On the use of kinetic energy preserving DG-schemes for large eddy simulation

NASA Astrophysics Data System (ADS)

Flad, David; Gassner, Gregor

2017-12-01

Recently, element based high order methods such as Discontinuous Galerkin (DG) methods and the closely related flux reconstruction (FR) schemes have become popular for compressible large eddy simulation (LES). Element based high order methods with Riemann solver based interface numerical flux functions offer an interesting dispersion dissipation behavior for multi-scale problems: dispersion errors are very low for a broad range of scales, while dissipation errors are very low for well resolved scales and are very high for scales close to the Nyquist cutoff. In some sense, the inherent numerical dissipation caused by the interface Riemann solver acts as a filter of high frequency solution components. This observation motivates the trend that element based high order methods with Riemann solvers are used without an explicit LES model added. Only the high frequency type inherent dissipation caused by the Riemann solver at the element interfaces is used to account for the missing sub-grid scale dissipation. Due to under-resolution of vortical dominated structures typical for LES type setups, element based high order methods suffer from stability issues caused by aliasing errors of the non-linear flux terms. A very common strategy to fight these aliasing issues (and instabilities) is so-called polynomial de-aliasing, where interpolation is exchanged with projection based on an increased number of quadrature points. In this paper, we start with this common no-model or implicit LES (iLES) DG approach with polynomial de-aliasing and Riemann solver dissipation and review its capabilities and limitations. We find that the strategy gives excellent results, but only when the resolution is such, that about 40% of the dissipation is resolved. For more realistic, coarser resolutions used in classical LES e.g. of industrial applications, the iLES DG strategy becomes quite inaccurate. We show that there is no obvious fix to this strategy, as adding for instance a sub-grid-scale models on top doesn't change much or in worst case decreases the fidelity even more. Finally, the core of this work is a novel LES strategy based on split form DG methods that are kinetic energy preserving. The scheme offers excellent stability with full control over the amount and shape of the added artificial dissipation. This premise is the main idea of the work and we will assess the LES capabilities of the novel split form DG approach when applied to shock-free, moderate Mach number turbulence. We will demonstrate that the novel DG LES strategy offers similar accuracy as the iLES methodology for well resolved cases, but strongly increases fidelity in case of more realistic coarse resolutions.

Biological consequences of a recurrent eddy off Point Conception, California

NASA Technical Reports Server (NTRS)

Haury, Loren R.; Simpson, James J.; Pelaez, Jose; Wisenhahn, David; Koblinsky, Chester J.

1986-01-01

The biological effects on three different time scales (100-day mesoscale, annual, and several-year) of a mesoscale anticyclonic eddy consistently found in shipboard surveys and satellite-sensed data several hundred kilometers southwest of Point Conception, CA, are described. A detailed shipboard study of the eddy in January 1981 found a complex system of fronts in surface chlorophyll at the northern edge of the eddy; microplankton and zooplankton distributions were strongly affected by entrainment processes at the surface and, apparently, at depth. Concurrent satellite coastal zone color scanner ocean color images show agreement with the general surface characteristics of the eddy chlorophyll field but do not reflect features deeper than about 25 m, including the contribution of the deep chlorophyll maximum to the integrated chlorophyll values. Satellite data for the period October 1980 through October 1981 and shipboard data from California Cooperative Oceanic Fisheries Investigations (CalCOFI) for December 1980 to July 1981 show the continued presence of the eddy in the sea surface temperature and color field and in the distributions of surface chlorophyll and zooplankton displacement volume. A review of the CalCOFI survey results from 1949 to the present time demonstrates the recurrent nature of the eddy system on a year-to-year basis. The eddy system appears to have a significant effect on the distribution of both oceanic and nearshore organisms. Offshore transport of coastal species occurs in the form of large entrained plumes or filaments.

Numerical study on the combustion characteristics of a fuel-centered pintle injector for methane rocket engines

NASA Astrophysics Data System (ADS)

Son, Min; Radhakrishnan, Kanmaniraja; Yoon, Youngbin; Koo, Jaye

2017-06-01

A pintle injector is a movable injector capable of controlling injection area and velocities. Although pintle injectors are not a new concept, they have become more notable due to new applications such as planet landers and low-cost engines. However, there has been little consistent research on pintle injectors because they have many design variations and mechanisms. In particular, simulation studies are required for bipropellant applications. In this study, combustion simulation was conducted using methane and oxygen to determine the effects of injection condition and geometries upon combustion characteristics. Steady and two-dimensional axisymmetric conditions were assumed and a 6-step Jones-Lindstedt mechanism with an eddy-dissipation concept model was used for turbulent kinetic reaction. As a result, the results with wide flame angles showed good combustion performances with a large recirculation under the pintle tip. Under lower mass flow-rate conditions, the combustion performance got worse with lower flame angles. To solve this problem, decreasing the pintle opening distance was very effective and the flame angle recovered. In addition, a specific recirculation zone was observed near the post, suggesting that proper design of the post could increase the combustion performance, while the geometry without a recirculation zone had the poor performance.

Physical and biological response of mesoscale eddies to wintertime forcing in the north central Red Sea (22˚N-25.5˚N)

NASA Astrophysics Data System (ADS)

Zarokanellos, N.; Jones, B. H.

2016-02-01

Red Sea is one of the saltiest and warmer seas in the world and acts as inverted estuary. Until recently, the Red Sea has been relatively underexplored. The limited observations that exist and results from various modeling exercises for the Red Sea have indicated that the sea has a complex mesoscale circulation often dominated by eddies. These mesoscale eddies are often visible in satellite imagery of sea surface height, temperature or chlorophyll, but only the surface expression of them. Because of previously limited in situ observations, the processes that drive the physical dynamics and the coupled biological responses have been poorly understood. To resolve and understand the role of these eddies in the dynamics of the north-central Red Sea during the wintertime, we used a combination of approaches that include remote sensing and autonomous underwater gliders equipped with physical, chemical, and bio-optical sensors. Remote sensing analyses of these eddies has shown that these eddies not only affect the physical circulation, but modify and disperse the phytoplankton populations and enhance exchange between the open sea and coastal coral reef ecosystems. During winter 2015, we observed deeper mixing driven by surface cooling and strong winds. As of result of the deeper mixing, phytoplankton populations became well mixed such that the ocean color imagery now reflected the integrated vertical processes. Localized diel fluctuations in phytoplankton are clearly evident during these well mixed periods. The mixing likely contributes to enhanced nutrient fluxes as well. Through sustained AUV observations, we have better understand the development, evolution, and dissipation of eddies. We also now have a better understanding of the mixing of source water from both the northern and southern Red Sea in this region of the north central Red Sea.

The solsticial pause on Mars: 2 modelling and investigation of causes

NASA Astrophysics Data System (ADS)

Mulholland, David P.; Lewis, Stephen R.; Read, Peter L.; Madeleine, Jean-Baptiste; Forget, Francois

2016-01-01

The martian solsticial pause, presented in a companion paper (Lewis et al., 2016), was investigated further through a series of model runs using the UK version of the LMD/UK Mars Global Climate Model. It was found that the pause could not be adequately reproduced if radiatively active water ice clouds were omitted from the model. When clouds were used, along with a realistic time-dependent dust opacity distribution, a substantial minimum in near-surface transient eddy activity formed around solstice in both hemispheres. The net effect of the clouds in the model is, by altering the thermal structure of the atmosphere, to decrease the vertical shear of the westerly jet near the surface around solstice, and thus reduce baroclinic growth rates. A similar effect was seen under conditions of large dust loading, implying that northern midlatitude eddy activity will tend to become suppressed after a period of intense flushing storm formation around the northern cap edge. Suppression of baroclinic eddy generation by the barotropic component of the flow and via diabatic eddy dissipation were also investigated as possible mechanisms leading to the formation of the solsticial pause but were found not to make major contributions. Zonal variations in topography were found to be important, as their presence results in weakened transient eddies around winter solstice in both hemispheres, through modification of the near-surface flow. The zonal topographic asymmetry appears to be the primary reason for the weakness of eddy activity in the southern hemisphere relative to the northern hemisphere, and the ultimate cause of the solsticial pause in both hemispheres. The meridional topographic gradient was found to exert a much weaker influence on near-surface transient eddies.

Wave Energetics of the Atmosphere of Mars

NASA Astrophysics Data System (ADS)

Battalio, Joseph Michael

A comprehensive assessment of the energetics of transient waves is presented for the atmosphere of Mars using the Mars Analysis Correction Data Assimilation (MACDA) dataset (v1.0) and the eddy kinetic energy equation. Each hemisphere is divided into four representative periods covering the summer and winter solstices, a late fall period, and an early spring period for each of the three Mars years available. Northern hemisphere 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 winter solstitial pause in wave activity. Barotropic energy conversion acts as a sink of eddy kinetic energy throughout each year with little reduction in amplitude during the solstitial 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 temperature profile around winter solstice. Traveling waves are typically triggered by geopotential flux convergence. Individual waves decay through a combination of barotropic conversion of the kinetic energy from the waves to the mean flow, geopotential flux divergence, and dissipation. The southern hemisphere energetics is similar to the northern hemisphere in timing, but wave energetics is much weaker as a result of the high and zonally asymmetric topography. The effect of dust on baroclinic instability is examined by comparing a year with a global-scale dust storm (GDS) to two years without a GDS. In the GDS year, waves develop a mixed baroclinic/barotropic growth phase before decaying barotropically. Though the total amount of eddy kinetic energy generated by baroclinic energy conversion is lower during the GDS year, the maximum eddy intensity is not diminished. Instead, the number of intense eddies is reduced by about 50%.

A Modulated-Gradient Parametrization for the Large-Eddy Simulation of the Atmospheric Boundary Layer Using the Weather Research and Forecasting Model

NASA Astrophysics Data System (ADS)

Khani, Sina; Porté-Agel, Fernando

2017-12-01

The performance of the modulated-gradient subgrid-scale (SGS) model is investigated using large-eddy simulation (LES) of the neutral atmospheric boundary layer within the weather research and forecasting model. Since the model includes a finite-difference scheme for spatial derivatives, the discretization errors may affect the simulation results. We focus here on understanding the effects of finite-difference schemes on the momentum balance and the mean velocity distribution, and the requirement (or not) of the ad hoc canopy model. We find that, unlike the Smagorinsky and turbulent kinetic energy (TKE) models, the calculated mean velocity and vertical shear using the modulated-gradient model, are in good agreement with Monin-Obukhov similarity theory, without the need for an extra near-wall canopy model. The structure of the near-wall turbulent eddies is better resolved using the modulated-gradient model in comparison with the classical Smagorinsky and TKE models, which are too dissipative and yield unrealistic smoothing of the smallest resolved scales. Moreover, the SGS fluxes obtained from the modulated-gradient model are much smaller near the wall in comparison with those obtained from the regular Smagorinsky and TKE models. The apparent inability of the LES model in reproducing the mean streamwise component of the momentum balance using the total (resolved plus SGS) stress near the surface is probably due to the effect of the discretization errors, which can be calculated a posteriori using the Taylor-series expansion of the resolved velocity field. Overall, we demonstrate that the modulated-gradient model is less dissipative and yields more accurate results in comparison with the classical Smagorinsky model, with similar computational costs.

Numerical Performance Prediction of a Miniature Ramjet at Mach 4

DTIC Science & Technology

2012-09-01

with the computational fluids dynamic (CFD) code from ANSYS - CFX . The nozzle-throat area was varied to increase the backpressure and this pushed the...normal shock that was sitting within the inlet, out to the lip of the inlet cowl. Using the eddy dissipation combustion model in ANSYS - CFX , a...improved accuracy in turbulence modeling. 14. SUBJECT TERMS Mach 4, Ramjet, Drag, Turbulence Modeling, Simulation, ANSYS CFX 15. NUMBER

On the Use of Linearized Euler Equations in the Prediction of Jet Noise

NASA Technical Reports Server (NTRS)

Mankbadi, Reda R.; Hixon, R.; Shih, S.-H.; Povinelli, L. A.

1995-01-01

Linearized Euler equations are used to simulate supersonic jet noise generation and propagation. Special attention is given to boundary treatment. The resulting solution is stable and nearly free from boundary reflections without the need for artificial dissipation, filtering, or a sponge layer. The computed solution is in good agreement with theory and observation and is much less CPU-intensive as compared to large-eddy simulations.

Damping and scattering of electromagnetic waves by small ferrite spheres suspended in an insulator

NASA Technical Reports Server (NTRS)

Englert, Gerald W.

1992-01-01

The intentional degradation of electromagnetic waves by their penetration into a media comprised of somewhat sparsely distributed energy absorbing ferrite spheres suspended in an electrical insulator is investigated. Results are presented in terms of generalized parameters involving wave length and sphere size, sphere resistivity, permeability, and spacing; their influence on dissipation of wave power by eddy currents, magnetic hysteresis, and scattering is shown.

An improved k-epsilon model for near wall turbulence

NASA Technical Reports Server (NTRS)

Shih, T. H.; Hsu, Andrew T.

1991-01-01

An improved k-epsilon model for low Reynolds number turbulence near a wall is presented. In the first part of this work, the near-wall asymptotic behavior of the eddy viscosity and the pressure transport term in the turbulent kinetic energy equation are analyzed. Based on these analyses, a modified eddy viscosity model with the 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, and a boundary condition for the dissipation rate is suggested. In the second part of the work, one of the deficiencies of the existing k-epsilon models, namely, the wall distance dependency of the equations and the damping functions, is examined. An improved model that does not depend on any wall distance is introduced. Fully developed turbulent channel flows and turbulent boundary layers over a flat plate are studied as validations for the proposed new models. Numerical results obtained from the present and other previous k-epsilon models are compared with data from direct numerical simulation. The results show that the present k-epsilon model, with added robustness, performs as well as or better than other existing models in predicting the behavior of near-wall turbulence.

Physical and biological characteristics of the winter-summer transition in the Central Red Sea

NASA Astrophysics Data System (ADS)

Zarokanellos, Nikolaos D.; Papadopoulos, Vassilis P.; Sofianos, Sarantis. S.; Jones, Burton H.

2017-08-01

The Central Red Sea (CRS) lies between two distinct hydrographic and atmospheric regimes. In the southern Red Sea, seasonal monsoon reversal regulates the exchange of water between the Red Sea and the Indian Ocean. In the northern Red Sea, intermediate and occasionally deep water are formed during winter to sustain the basin's overturning circulation. Highly variable mesoscale eddies and the northward flowing eastern boundary current (EBC) determine the physical and biogeochemical characteristics of the CRS. Ship-based and glider observations in the CRS between March and June 2013 capture key features of the transition from winter to summer and depict the impact of the eddy activity on the EBC flow. Less saline and relatively warmer water of Indian Ocean origin reaches the CRS via the EBC. Initially, an anticyclonic eddy with diameter of 140 km penetrating to 150m depth with maximum velocities up to 30-35 cm s-1 prevails in the CRS. This anticyclonic eddy appears to block or at least redirect the northward flow of the EBC. Dissipation of the eddy permits the near-coastal, northward flow of the EBC and gives place to a smaller cyclonic eddy with a diameter of about 50 km penetrating to 200 m depth. By the end of May, as the northerly winds become stronger and persistent throughout the basin, characteristic of the summer southwest monsoon wind regime, the EBC, and its associated lower salinity water became less evident, replaced by the saltier surface water that characterizes the onset of the summer stratification in the CRS.

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.

Advanced k-epsilon modeling of heat transfer

NASA Technical Reports Server (NTRS)

Kwon, Okey; Ames, Forrest E.

1995-01-01

This report describes two approaches to low Reynolds-number k-epsilon turbulence modeling which formulate the eddy viscosity on the wall-normal component of turbulence and a length scale. The wall-normal component of turbulence is computed via integration of the energy spectrum based on the local dissipation rate and is bounded by the isotropic condition. The models account for the anisotropy of the dissipation and the reduced mixing length due to the high strain rates present in the near-wall region. The turbulent kinetic energy and its dissipation rate were computed from the k and epsilon transport equations of Durbin. The models were tested for a wide range of turbulent flows and proved to be superior to other k-epsilon models, especially for nonequilibrium anisotropic flows. For the prediction of airfoil heat transfer, the models included a set of empirical correlations for predicting laminar-turbulent transition and laminar heat transfer augmentation due to the presence of freestream turbulence. The predictions of surface heat transfer were generally satisfactory.

LLWBCS changes through surface mesoscale activity and baroclinic tides in the Solomon Sea

NASA Astrophysics Data System (ADS)

Gourdeau, L.; Djath, B.; Ganachaud, A. S.; Tchilibou, M. L.; Verron, J. A.; Jouanno, J.

2016-02-01

In the south west Pacific, the Solomon Sea is on the pathway of the Low Latitudes Western Boundary Currents that connect the subtropics to the equator. Changes in their strengths, or in their water mass properties may have implication for ENSO and its low frequency modulation. During their transit in the Solomon Sea, the salinity maximum at thermocline level, characteristic of the South Pacific Tropical Waters (SPTW), is largely eroded. Different mechanisms could explain such salt erosion whose current/bathymetry interaction, internal tides, eddy activity. The Solomon Sea is an area of high level of eddy kinetic energy (EKE), especially in the surface layers, and its complex bathymetry is favourable for generation and dissipation of internal tides. Based on high resolution modelling, glider, and altimetric data mesoscale eddies observed at the surface are analysed in their 4D aspects. Their role on water mass transformation is explored. These eddies may affect the surface layers (σ<23.3) and the upper thermocline waters (23.3< σ <24.3), but they cannot explained the erosion of the salinity maximum below. Simulations with and without explicit tides provide a description of baroclinic tides in the Solomon Sea. Their role on water mixing is evaluated, especially for the SPTW.

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.

Non-Markovian closure models for large eddy simulations using the Mori-Zwanzig formalism

NASA Astrophysics Data System (ADS)

Parish, Eric J.; Duraisamy, Karthik

2017-01-01

This work uses the Mori-Zwanzig (M-Z) formalism, a concept originating from nonequilibrium statistical mechanics, as a basis for the development of coarse-grained models of turbulence. The mechanics of the generalized Langevin equation (GLE) are considered, and insight gained from the orthogonal dynamics equation is used as a starting point for model development. A class of subgrid models is considered which represent nonlocal behavior via a finite memory approximation [Stinis, arXiv:1211.4285 (2012)], the length of which is determined using a heuristic that is related to the spectral radius of the Jacobian of the resolved variables. The resulting models are intimately tied to the underlying numerical resolution and are capable of approximating non-Markovian effects. Numerical experiments on the Burgers equation demonstrate that the M-Z-based models can accurately predict the temporal evolution of the total kinetic energy and the total dissipation rate at varying mesh resolutions. The trajectory of each resolved mode in phase space is accurately predicted for cases where the coarse graining is moderate. Large eddy simulations (LESs) of homogeneous isotropic turbulence and the Taylor-Green Vortex show that the M-Z-based models are able to provide excellent predictions, accurately capturing the subgrid contribution to energy transfer. Last, LESs of fully developed channel flow demonstrate the applicability of M-Z-based models to nondecaying problems. It is notable that the form of the closure is not imposed by the modeler, but is rather derived from the mathematics of the coarse graining, highlighting the potential of M-Z-based techniques to define LES closures.

Fluxes of carbon, water and energy over Brazilian cerrado: an analysis using eddy covariance and stable isotopes

Treesearch

A. C. Miranda; H. S. Miranda; J. Lloyd; J. Grace; R. J. Francey; J. A. Mcintyre; P. Meir; P. Riggan; R. Lockwood; J. Brass

1997-01-01

We present the energy and mass balance of cerrado sensu stricto (a Brazilian form of savanna), in which a mixture of shrubs, trees and grasses forms a vegetation with a leaf area index of 1·0 in the wet season and 0·4 in the dry season. In the wet season the available energy was equally dissipated between sensible heat and...

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.

The relationship between a deformation-based eddy parameterization and the LANS-α turbulence model

NASA Astrophysics Data System (ADS)

Bachman, Scott D.; Anstey, James A.; Zanna, Laure

2018-06-01

A recent class of ocean eddy parameterizations proposed by Porta Mana and Zanna (2014) and Anstey and Zanna (2017) modeled the large-scale flow as a non-Newtonian fluid whose subgridscale eddy stress is a nonlinear function of the deformation. This idea, while largely new to ocean modeling, has a history in turbulence modeling dating at least back to Rivlin (1957). The new class of parameterizations results in equations that resemble the Lagrangian-averaged Navier-Stokes-α model (LANS-α, e.g., Holm et al., 1998a). In this note we employ basic tensor mathematics to highlight the similarities between these turbulence models using component-free notation. We extend the Anstey and Zanna (2017) parameterization, which was originally presented in 2D, to 3D, and derive variants of this closure that arise when the full non-Newtonian stress tensor is used. Despite the mathematical similarities between the non-Newtonian and LANS-α models which might provide insight into numerical implementation, the input and dissipation of kinetic energy between these two turbulent models differ.

Results from flamelet and non-flamelet models for supersonic combustion

NASA Astrophysics Data System (ADS)

Ladeinde, Foluso; Li, Wenhai

2017-11-01

Air-breathing propulsion systems (scramjets) have been identified as a viable alternative to rocket engines for improved efficiency. A scramjet engine, which operates at flight Mach numbers around 7 or above, is characterized by the existence of supersonic flow conditions in the combustor. In a dual-mode scramjet, this phenomenon is possible because of the relatively low value of the equivalence ratio and high stagnation temperature, which, together, inhibits thermal choking downstream of transverse injectors. The flamelet method has been our choice for turbulence-combustion interaction modeling and we have extended the basic approach in several dimensions, with a focus on the way the pressure and progress variable are modeled. Improved results have been obtained. We have also examined non-flamelet models, including laminar chemistry (QL), eddy dissipation concept (EDC), and partially-stirred reactor (PaSR). The pressure/progress variable-corrected simulations give better results compared with the original model, with reaction rates that are lower than those from EDC and PaSR. In general, QL tends to over-predict the reaction rate for the supersonic combustion problems investigated in our work.

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.

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.

Surrogates for numerical simulations; optimization of eddy-promoter heat exchangers

NASA Technical Reports Server (NTRS)

Patera, Anthony T.; Patera, Anthony

1993-01-01

Although the advent of fast and inexpensive parallel computers has rendered numerous previously intractable calculations feasible, many numerical simulations remain too resource-intensive to be directly inserted in engineering optimization efforts. An attractive alternative to direct insertion considers models for computational systems: the expensive simulation is evoked only to construct and validate a simplified, input-output model; this simplified input-output model then serves as a simulation surrogate in subsequent engineering optimization studies. A simple 'Bayesian-validated' statistical framework for the construction, validation, and purposive application of static computer simulation surrogates is presented. As an example, dissipation-transport optimization of laminar-flow eddy-promoter heat exchangers are considered: parallel spectral element Navier-Stokes calculations serve to construct and validate surrogates for the flowrate and Nusselt number; these surrogates then represent the originating Navier-Stokes equations in the ensuing design process.

Large eddy simulation of cavitating flows

NASA Astrophysics Data System (ADS)

Gnanaskandan, Aswin; Mahesh, Krishnan

2014-11-01

Large eddy simulation on unstructured grids is used to study hydrodynamic cavitation. The multiphase medium is represented using a homogeneous equilibrium model that assumes thermal equilibrium between the liquid and the vapor phase. Surface tension effects are ignored and the governing equations are the compressible Navier Stokes equations for the liquid/vapor mixture along with a transport equation for the vapor mass fraction. A characteristic-based filtering scheme is developed to handle shocks and material discontinuities in non-ideal gases and mixtures. A TVD filter is applied as a corrector step in a predictor-corrector approach with the predictor scheme being non-dissipative and symmetric. The method is validated for canonical one dimensional flows and leading edge cavitation over a hydrofoil, and applied to study sheet to cloud cavitation over a wedge. This work is supported by the Office of Naval Research.

Large-eddy simulations of compressible convection on massively parallel computers. [stellar physics

NASA Technical Reports Server (NTRS)

Xie, Xin; Toomre, Juri

1993-01-01

We report preliminary implementation of the large-eddy simulation (LES) technique in 2D simulations of compressible convection carried out on the CM-2 massively parallel computer. The convective flow fields in our simulations possess structures similar to those found in a number of direct simulations, with roll-like flows coherent across the entire depth of the layer that spans several density scale heights. Our detailed assessment of the effects of various subgrid scale (SGS) terms reveals that they may affect the gross character of convection. Yet, somewhat surprisingly, we find that our LES solutions, and another in which the SGS terms are turned off, only show modest differences. The resulting 2D flows realized here are rather laminar in character, and achieving substantial turbulence may require stronger forcing and less dissipation.

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.

Large Eddy Simulations of Supercritical Mixing Layers for Air Force Applications

DTIC Science & Technology

2010-05-01

Taskinoglu and J. Bellan Species m( gmo \\-]) 7). (K) P,(MPa) N2 28.013 126.3 3.399 C7H16 100.205...equation is that its source term, called the irreversible en- tropy production, which is by definition is the dissipation [3|, con - tains the full extent...with sim- ilar or even larger gradient magnitudes under fully turbulent con - ditions (the experimental data was for Re = O(104) - O(105)). Thus

An innovative experimental sequence on electromagnetic induction and eddy currents based on video analysis and cheap data acquisition

NASA Astrophysics Data System (ADS)

Bonanno, A.; Bozzo, G.; Sapia, P.

2017-11-01

In this work, we present a coherent sequence of experiments on electromagnetic (EM) induction and eddy currents, appropriate for university undergraduate students, based on a magnet falling through a drilled aluminum disk. The sequence, leveraging on the didactical interplay between the EM and mechanical aspects of the experiments, allows us to exploit the students’ awareness of mechanics to elicit their comprehension of EM phenomena. The proposed experiments feature two kinds of measurements: (i) kinematic measurements (performed by means of high-speed video analysis) give information on the system’s kinematics and, via appropriate numerical data processing, allow us to get dynamic information, in particular on energy dissipation; (ii) induced electromagnetic field (EMF) measurements (by using a homemade multi-coil sensor connected to a cheap data acquisition system) allow us to quantitatively determine the inductive effects of the moving magnet on its neighborhood. The comparison between experimental results and the predictions from an appropriate theoretical model (of the dissipative coupling between the moving magnet and the conducting disk) offers many educational hints on relevant topics related to EM induction, such as Maxwell’s displacement current, magnetic field flux variation, and the conceptual link between induced EMF and induced currents. Moreover, the didactical activity gives students the opportunity to be trained in video analysis, data acquisition and numerical data processing.

Eddy-Kuroshio Interactions: Local and Remote Effects

NASA Astrophysics Data System (ADS)

Jan, Sen; Mensah, Vigan; Andres, Magdalena; Chang, Ming-Huei; Yang, Yiing Jang

2017-12-01

Quasi-geostrophic mesoscale eddies regularly impinge on the Kuroshio in the western North Pacific, but the processes underlying the evolution of these eddy-Kuroshio interactions have not yet been thoroughly investigated in the literature. Here this interaction is examined with results from a semi-idealized three-dimensional numerical model and observations from four pressure-sensor equipped inverted echo sounders (PIESs) in a zonal section east of Taiwan and satellite altimeters. Both the observations and numerical simulations suggest that, during the interaction of a cyclonic eddy with the Kuroshio, the circular eddy is deformed into an elliptic shape with the major axis in the northwest-southeast direction, before being dissipated; the poleward velocity and associated Kuroshio transport decrease and the sea level and pycnocline slopes across the Kuroshio weaken. In contrast, for an anticyclonic eddy during the eddy-Kuroshio interaction, variations in the velocity, sea level, and isopycnal depth are reversed; the circular eddy is also deformed to an ellipse but with the major axis parallel to the Kuroshio. The model results also demonstrate that the velocity field is modified first and consequently the SSH and isopycnal depth evolve during the interaction. Furthermore, due to the combined effect of impingement latitude and realistic topography, some eddy-Kuroshio interactions east of Taiwan are found to have remote effects, both in the Luzon Strait and on the East China Sea shelf northeast of Taiwan.