The SGS kinetic energy and the viscous dissipation equations as closure relations in LES

Author(s):  
F. Gallerano ◽  
L. Melilla ◽  
E. Pasero
2021 ◽  
Vol 7 (1) ◽  
Author(s):  
Antoine Riaud ◽  
Cui Wang ◽  
Jia Zhou ◽  
Wanghuai Xu ◽  
Zuankai Wang

AbstractElectric energy generation from falling droplets has seen a hundred-fold rise in efficiency over the past few years. However, even these newest devices can only extract a small portion of the droplet energy. In this paper, we theoretically investigate the contributions of hydrodynamic and electric losses in limiting the efficiency of droplet electricity generators (DEG). We restrict our analysis to cases where the droplet contacts the electrode at maximum spread, which was observed to maximize the DEG efficiency. Herein, the electro-mechanical energy conversion occurs during the recoil that immediately follows droplet impact. We then identify three limits on existing droplet electric generators: (i) the impingement velocity is limited in order to maintain the droplet integrity; (ii) much of droplet mechanical energy is squandered in overcoming viscous shear force with the substrate; (iii) insufficient electrical charge of the substrate. Of all these effects, we found that up to 83% of the total energy available was lost by viscous dissipation during spreading. Minimizing this loss by using cascaded DEG devices to reduce the droplet kinetic energy may increase future devices efficiency beyond 10%.


In this problem a mean turbulent shear layer originally exists, homogeneous in the streamwise direction, formed perhaps by previous instabilities, but in equilibrium with the fine-grained turbulence. At a given time, a large eddy of a fixed horizontal wavenumber is initiated. We study the subsequent time development of the non-equilibrium interactions between the three components of flow as they adjust towards ultimate simultaneous equilibrium, using the integrated energy-balance conservation equations to derive the amplitude equations. This necessarily involves the usual averaging procedure and a conditional or phase-averaging procedure by which the large structure motion is educed from the total fluctuations. In general, the mean flow growth is due to the energy transfer to both fluctuating components, the large eddy gains energy from the mean motion and exchanges energy with the fine-grained turbulence, while the fine-grained turbulence gains energy from the mean flow and exchanges with the large eddy and converts its energy to heat through viscous dissipation of the smallest scales. The closure problem is obtained via the shape assumptions which enter into the interaction integrals. The situation in which the fine-grained turbulent kinetic energy production and viscous dissipation are in local balance is considered, the displacement from equilibrium being due only to the energy transfer from the large eddy. The large eddy shape is taken to be two-dimensional, instability-wavelike, with its vorticity axis perpendicular to the direction of the mean outer stream. Prior to averaging, detailed but approximate calculations of the wave-induced turbulent Reynolds stresses are obtained; the product of these stresses with the appropriate large-eddy rates of strain give the energy transfer mechanism between the two disparate scales of fluctuations. Coupled, nonlinear amplitude or energy density equations for the three components of motion are obtained, the coefficients of which are the interaction integrals guided by the shape assumptions. It is found that for the special case of parallel flow, the energy of the large eddy first undergoes a hydrodynamic-instability type of amplification but eventually decays due to the energy transfer to the fine-grained turbulence, while the turbulent kinetic energy is displaced from an original level of equilibrium to a new one because of the ability of the large eddy to negotiate an indirect energy transfer from the mean flow. For the growing shear layer, approximate considerations show that if the mechanism of energy transfer from the large to the small scale is eventually weakened by the shear layer growth compared to the large-eddy production mechanism so that the amplification and decay process repeats, ‘bursts’ of the remnant of the same large eddy will occur repeatedly until an ultimate equilibrium is reached among the three interacting components of motion. However, for the large eddy whose wavenumber corresponds to that of the initially most amplified case, the ‘bursting’ phenomenon is much less pronounced and equilibrium is very nearly reached at the end of the very first ‘burst’.


2015 ◽  
Vol 773 ◽  
pp. 395-417 ◽  
Author(s):  
K. Petschel ◽  
S. Stellmach ◽  
M. Wilczek ◽  
J. Lülff ◽  
U. Hansen

The kinetic energy balance in Rayleigh–Bénard convection is investigated by means of direct numerical simulations for the Prandtl number range $0.01\leqslant \mathit{Pr}\leqslant 150$ and for fixed Rayleigh number $\mathit{Ra}=5\times 10^{6}$. The kinetic energy balance is divided into a dissipation, a production and a flux term. We discuss the profiles of all the terms and find that the different contributions to the energy balance can be spatially separated into regions where kinetic energy is produced and where kinetic energy is dissipated. By analysing the Prandtl number dependence of the kinetic energy balance, we show that the height dependence of the mean viscous dissipation is closely related to the flux of kinetic energy. We show that the flux of kinetic energy can be divided into four additive contributions, each representing a different elementary physical process (advection, buoyancy, normal viscous stresses and viscous shear stresses). The behaviour of these individual flux contributions is found to be surprisingly rich and exhibits a pronounced Prandtl number dependence. Different flux contributions dominate the kinetic energy transport at different depths, such that a comprehensive discussion requires a decomposition of the domain into a considerable number of sublayers. On a less detailed level, our results reveal that advective kinetic energy fluxes play a key role in balancing the near-wall dissipation at low Prandtl number, whereas normal viscous stresses are particularly important at high Prandtl number. Finally, our work reveals that classical velocity boundary layers are deeply connected to the kinetic energy transport, but fail to correctly represent regions of enhanced viscous dissipation.


2018 ◽  
Vol 23 (3) ◽  
pp. 787-801
Author(s):  
B. Zigta

Abstract This study examines the effect of thermal radiation, chemical reaction and viscous dissipation on a magnetohydro- dynamic flow in between a pair of infinite vertical Couette channel walls. The momentum equation accounts the effects of both the thermal and the concentration buoyancy forces of the flow. The energy equation addresses the effects of the thermal radiation and viscous dissipation of the flow. Also, the concentration equation includes the effects of molecular diffusivity and chemical reaction parameters. The gray colored fluid considered in this study is a non-scattering medium and has the property of absorbing and emitting radiation. The Roseland approximation is used to describe the radiative heat flux in the energy equation. The velocity of flow transforms kinetic energy into heat energy. The increment of the velocity due to internal energy results in heating up of the fluid and consequently it causes increment of the thermal buoyancy force. The Eckert number being the ratio of the kinetic energy of the flow to the temperature difference of the channel walls is directly proportional to the thermal energy dissipation. It can be observed that increasing the Eckert number results in increasing velocity. A uniform magnetic field is applied perpendicular to the channel walls. The temperature of the moving wall is high enough due to the presence of thermal radiation. The solution of the governing equations is obtained using regular perturbation techniques. These techniques help to convert partial differential equations to a set of ordinary differential equations in dimensionless form and thus they are solved analytically. The following results are obtained: from the simulation study it is observed that the flow pattern of the fluid is affected due to the influence of the thermal radiation, the chemical reaction and viscous dissipation. The increment in the Hartmann number results in the increment of the Lorentz force but a decrement in velocity of the flow. An increment in the radiative parameter results in a decrement in temperature. An increment in the Prandtl number results in a decrement in thermal diffusivity. An increment in both the chemical reaction parameter and molecular diffusivity results in a decrement in concentration.


2015 ◽  
Vol 767 ◽  
pp. 342-363 ◽  
Author(s):  
J. Bertram

AbstractFollowing Malkus’s (J. Fluid Mech., vol. 1, 1956, pp. 521–539) proposal that turbulent Poiseuille channel flow maximises total viscous dissipation $D$, a variety of variational procedures have been explored involving the maximisation of different flow quantities under different constraints. However, the physical justification for these variational procedures has remained unclear. Here we address more recent claims that mean flow viscous dissipation $D_{m}$ should be maximised on the basis of a statistical stability argument, and that maximising $D_{m}$ yields realistic mean velocity profiles (Malkus, J. Fluid Mech., vol. 489, 2003, pp. 185–198). We clarify the connection between maximising $D_{m}$ and other flow quantities, verify Malkus & Smith’s, (J. Fluid Mech., vol. 208, 1989, pp. 479–507) claim that maximising the ‘efficiency’ yields realistic profiles and show that, in contrast, maximising $D_{m}$ does not yield realistic mean velocity profiles as recently claimed. This leads us to revisit Malkus’s statistical stability argument for maximising $D_{m}$ and to address some of its limitations. We propose an alternative statistical stability argument leading to a principle of minimum kinetic energy for fixed pressure gradient, which suggests a principle of maximum $D$ for fixed Reynolds number under certain conditions. We discuss possible ways to reconcile these conflicting results, focusing on the choice of constraints.


2017 ◽  
Vol 140 (1) ◽  
Author(s):  
Sanghun Choi ◽  
Jiwoong Choi ◽  
Ching-Long Lin

The aim of this study was to investigate and quantify contributions of kinetic energy and viscous dissipation to airway resistance during inspiration and expiration at various flow rates in airway models of different bifurcation angles. We employed symmetric airway models up to the 20th generation with the following five different bifurcation angles at a tracheal flow rate of 20 L/min: 15 deg, 25 deg, 35 deg, 45 deg, and 55 deg. Thus, a total of ten computational fluid dynamics (CFD) simulations for both inspiration and expiration were conducted. Furthermore, we performed additional four simulations with tracheal flow rate values of 10 and 40 L/min for a bifurcation angle of 35 deg to study the effect of flow rate on inspiration and expiration. Using an energy balance equation, we quantified contributions of the pressure drop associated with kinetic energy and viscous dissipation. Kinetic energy was found to be a key variable that explained the differences in airway resistance on inspiration and expiration. The total pressure drop and airway resistance were larger during expiration than inspiration, whereas wall shear stress and viscous dissipation were larger during inspiration than expiration. The dimensional analysis demonstrated that the coefficients of kinetic energy and viscous dissipation were strongly correlated with generation number. In addition, the viscous dissipation coefficient was significantly correlated with bifurcation angle and tracheal flow rate. We performed multiple linear regressions to determine the coefficients of kinetic energy and viscous dissipation, which could be utilized to better estimate the pressure drop in broader ranges of successive bifurcation structures.


2010 ◽  
Vol 132 (6) ◽  
Author(s):  
Ri Li ◽  
Nasser Ashgriz ◽  
Sanjeev Chandra

This theoretical study proposes an analytical model to predict the maximum spread of single droplets on solid surfaces with zero or low Weber and Reynolds numbers. The spreading droplet is assumed as a spherical cap considering low impact velocities. Three spreading states are considered, which include equilibrium spread, maximum spontaneous spread, and maximum spread. Energy conservation is applied to the droplet as a control volume. The model equation contains two viscous dissipation terms, each of which has a defined coefficient. One term is for viscous dissipation in spontaneous spreading and the other one is for viscous dissipation of the initial kinetic energy of the droplet. The new model satisfies the fundamental physics of drop-surface interaction and can be used for droplets impacting on solid surfaces with or without initial kinetic energy.


2019 ◽  
Vol 141 (12) ◽  
Author(s):  
T.-W. Lee ◽  
J. E. Park

We have used the integral form of the conservation equations, to find a cubic formula for the drop size during in liquid sprays in coflow of air (air-blast atomization). Similar to our previous work, the energy balance dictates that the initial kinetic energy of the gas and injected liquid will be distributed into the final surface tension energy, kinetic energy of the gas and droplets, and viscous dissipation. Using this approach, the drop size can be determined based on the basic injection and fluid parameters for “air-blast” atomization, where the injected liquid is atomized by high-speed coflow of air. The viscous dissipation term is estimated using appropriate velocity and length scales of liquid–air coflow breakup. The mass and energy balances for the spray flows render to an expression that relates the drop size to all of the relevant parameters, including the gas- and liquid-phase velocities and fluid properties. The results agree well with experimental data and correlations for the drop size. The solution also provides for drop size–velocity cross-correlation, leading to computed drop size distributions based on the gas-phase velocity distribution. This approach can be used in the estimation of the drop size for practical sprays and also as a primary atomization module in computational simulations of air-blast atomization over a wide range of injection and fluid conditions, the only caveat being that a parameter to account for the viscous dissipation needs to be calibrated with a minimal set of observational data.


Sign in / Sign up

Export Citation Format

Share Document