scholarly journals Numerical simulation of turbulent channel flow over a viscous hyper-elastic wall

2017 ◽  
Vol 830 ◽  
pp. 708-735 ◽  
Author(s):  
Marco E. Rosti ◽  
Luca Brandt
Keyword(s):  
Channel Flow ◽  
Solid Phase ◽  
Turbulent Channel Flow ◽  
Inertial Range ◽  
Spanwise Direction ◽  
Logarithmic Scale ◽  
Normal Velocity ◽  
Mean Velocity ◽  
Elastic Wall ◽  
Hyper Elastic

We perform numerical simulations of a turbulent channel flow over an hyper-elastic wall. In the fluid region the flow is governed by the incompressible Navier–Stokes (NS) equations, while the solid is a neo-Hookean material satisfying the incompressible Mooney–Rivlin law. The multiphase flow is solved with a one-continuum formulation, using a monolithic velocity field for both the fluid and solid phase, which allows the use of a fully Eulerian formulation. The simulations are carried out at Reynolds bulk $Re=2800$ and examine the effect of different elasticity and viscosity of the deformable wall. We show that the skin friction increases monotonically with the material elastic modulus. The turbulent flow in the channel is affected by the moving wall even at low values of elasticity since non-zero fluctuations of vertical velocity at the interface influence the flow dynamics. The near-wall streaks and the associated quasi-streamwise vortices are strongly reduced near a highly elastic wall while the flow becomes more correlated in the spanwise direction, similarly to what happens for flows over rough and porous walls. As a consequence, the mean velocity profile in wall units is shifted downwards when shown in logarithmic scale, and the slope of the inertial range increases in comparison to that for the flow over a rigid wall. We propose a correlation between the downward shift of the inertial range, its slope and the wall-normal velocity fluctuations at the wall, extending results for the flow over rough walls. We finally show that the interface deformation is determined by the fluid fluctuations when the viscosity of the elastic layer is low, while when this is high the deformation is limited by the solid properties.

scholarly journals Turbulent channel flow over an anisotropic porous wall – drag increase and reduction

2018 ◽  
Vol 842 ◽  
pp. 381-394 ◽  
Author(s):  
Marco E. Rosti ◽  
Luca Brandt ◽  
Alfredo Pinelli
Keyword(s):  
Drag Reduction ◽  
Channel Flow ◽  
High Speed ◽  
Vertical Direction ◽  
Porous Wall ◽  
Turbulent Channel Flow ◽  
Spanwise Direction ◽  
Normal Velocity ◽  
Porous Walls ◽  
The One

The effect of the variations of the permeability tensor on the close-to-the-wall behaviour of a turbulent channel flow bounded by porous walls is explored using a set of direct numerical simulations. It is found that the total drag can be either reduced or increased by more than 20 % by adjusting the permeability directional properties. Drag reduction is achieved for the case of materials with permeability in the vertical direction lower than the one in the wall-parallel planes. This configuration limits the wall-normal velocity at the interface while promoting an increase of the tangential slip velocity leading to an almost ‘one-component’ turbulence where the low- and high-speed streak coherence is strongly enhanced. On the other hand, strong drag increase is found when high wall-normal and low wall-parallel permeabilities are prescribed. In this condition, the enhancement of the wall-normal fluctuations due to the reduced wall-blocking effect triggers the onset of structures which are strongly correlated in the spanwise direction, a phenomenon observed by other authors in flows over isotropic porous layers or over ribletted walls with large protrusion heights. The use of anisotropic porous walls for drag reduction is particularly attractive since equal gains can be achieved at different Reynolds numbers by rescaling the magnitude of the permeability only.

The structure of the viscous sublayer and the adjacent wall region in a turbulent channel flow

1974 ◽  
Vol 65 (3) ◽  
pp. 439-459 ◽  
Author(s):  
Helmut Eckelmann
Keyword(s):  
Channel Flow ◽  
Turbulent Channel Flow ◽  
Streamwise Velocity ◽  
Viscous Sublayer ◽  
Wall Region ◽  
Normal Velocity ◽  
Mean Velocity ◽  
Velocity Fluctuations ◽  
Mean Values ◽  
The Mean

Hot-film anemometer measurements have been carried out in a fully developed turbulent channel flow. An oil channel with a thick viscous sublayer was used, which permitted measurements very close to the wall. In the viscous sublayer between y+ ≃ 0·1 and y+ = 5, the streamwise velocity fluctuations decreased at a higher rate than the mean velocity; in the region y+ [lsim ] 0·1, these fluctuations vanished at the same rate as the mean velocity.The streamwise velocity fluctuations u observed in the viscous sublayer and the fluctuations (∂u/∂y)0 of the gradient at the wall were almost identical in form, but the fluctuations of the gradient at the wall were found to lag behind the velocity fluctuations with a lag time proportional to the distance from the wall. Probability density distributions of the streamwise velocity fluctuations were measured. Furthermore, measurements of the skewness and flatness factors made by Kreplin (1973) in the same flow channel are discussed. Measurements of the normal velocity fluctuations v at the wall and of the instantaneous Reynolds stress −ρuv were also made. Periods of quiescence in the − ρuv signal were observed in the viscous sublayer as well as very active periods where ratios of peak to mean values as high as 30:1 occurred.

Fluid particle dynamics and the non-local origin of the Reynolds shear stress

2018 ◽  
Vol 847 ◽  
pp. 520-551 ◽  
Author(s):  
Peter S. Bernard ◽  
Martin A. Erinin
Keyword(s):  
Shear Stress ◽  
Channel Flow ◽  
Reynolds Shear Stress ◽  
Mixing Time ◽  
Turbulent Channel Flow ◽  
Structural Features ◽  
Fluid Particle ◽  
Normal Velocity ◽  
Mean Velocity ◽  
Non Local

The causative factors leading to the Reynolds shear stress distribution in turbulent channel flow are analysed via a backward particle path analysis. It is found that the classical displacement transport mechanism, by which changes in the mean velocity field over a mixing time correlate with the wall-normal velocity, is the dominant source of Reynolds shear stress. Approximately 20 % of channel flow at any given time contains fluid motions that contribute to displacement transport. Much rarer events provide a small but non-negligible contribution to the Reynolds shear stress due to fluid particle accelerations and long-lived correlations deriving from structural features of the near-wall flow. The Reynolds shear stress in channel flow is shown to be a non-local phenomenon that is not conducive to description via a local model and particularly one depending directly on the local mean velocity gradient.

Direct numerical simulation of turbulent channel flow up to

2015 ◽  
Vol 774 ◽  
pp. 395-415 ◽  
Author(s):  
Myoungkyu Lee ◽  
Robert D. Moser
Keyword(s):  
Numerical Simulation ◽  
Reynolds Number ◽  
Direct Numerical Simulation ◽  
Turbulent Flows ◽  
Channel Flow ◽  
Turbulent Channel Flow ◽  
Streamwise Velocity ◽  
Mean Velocity ◽  
Layer Structures ◽  
High Reynolds

A direct numerical simulation of incompressible channel flow at a friction Reynolds number ($\mathit{Re}_{{\it\tau}}$) of 5186 has been performed, and the flow exhibits a number of the characteristics of high-Reynolds-number wall-bounded turbulent flows. For example, a region where the mean velocity has a logarithmic variation is observed, with von Kármán constant ${\it\kappa}=0.384\pm 0.004$. There is also a logarithmic dependence of the variance of the spanwise velocity component, though not the streamwise component. A distinct separation of scales exists between the large outer-layer structures and small inner-layer structures. At intermediate distances from the wall, the one-dimensional spectrum of the streamwise velocity fluctuation in both the streamwise and spanwise directions exhibits $k^{-1}$ dependence over a short range in wavenumber $(k)$. Further, consistent with previous experimental observations, when these spectra are multiplied by $k$ (premultiplied spectra), they have a bimodal structure with local peaks located at wavenumbers on either side of the $k^{-1}$ range.

scholarly journals Coherent structures in the linearized impulse response of turbulent channel flow

2019 ◽  
Vol 863 ◽  
pp. 1190-1203 ◽  
Author(s):  
Sabarish B. Vadarevu ◽  
Sean Symon ◽  
Simon J. Illingworth ◽  
Ivan Marusic
Keyword(s):  
Channel Flow ◽  
Coherent Structures ◽  
Large Scale ◽  
Body Force ◽  
Stokes Equations ◽  
Turbulent Channel Flow ◽  
Kinetic Energy Density ◽  
Reynolds Stresses ◽  
Mean Velocity ◽  
Velocity Fluctuations

We study the evolution of velocity fluctuations due to an isolated spatio-temporal impulse using the linearized Navier–Stokes equations. The impulse is introduced as an external body force in incompressible channel flow at $Re_{\unicode[STIX]{x1D70F}}=10\,000$. Velocity fluctuations are defined about the turbulent mean velocity profile. A turbulent eddy viscosity is added to the equations to fix the mean velocity as an exact solution, which also serves to model the dissipative effects of the background turbulence on large-scale fluctuations. An impulsive body force produces flow fields that evolve into coherent structures containing long streamwise velocity streaks that are flanked by quasi-streamwise vortices; some of these impulses produce hairpin vortices. As these vortex–streak structures evolve, they grow in size to be nominally self-similar geometrically with an aspect ratio (streamwise to wall-normal) of approximately 10, while their kinetic energy density decays monotonically. The topology of the vortex–streak structures is not sensitive to the location of the impulse, but is dependent on the direction of the impulsive body force. All of these vortex–streak structures are attached to the wall, and their Reynolds stresses collapse when scaled by distance from the wall, consistent with Townsend’s attached-eddy hypothesis.

Relationship between the heat transfer law and the scalar dissipation function in a turbulent channel flow

2017 ◽  
Vol 830 ◽  
pp. 300-325 ◽  
Author(s):  
Hiroyuki Abe ◽  
Robert Anthony Antonia
Keyword(s):  
Heat Transfer ◽  
Prandtl Number ◽  
Channel Flow ◽  
Dissipation Rate ◽  
Turbulent Channel Flow ◽  
Scalar Dissipation ◽  
Scalar Dissipation Rate ◽  
Mean Velocity ◽  
The Mean ◽  
Logarithmic Dependence

Integration across a fully developed turbulent channel flow of the transport equations for the mean and turbulent parts of the scalar dissipation rate yields relatively simple relations for the bulk mean scalar and wall heat transfer coefficient. These relations are tested using direct numerical simulation datasets obtained with two isothermal boundary conditions (constant heat flux and constant heating source) and a molecular Prandtl number Pr of 0.71. A logarithmic dependence on the Kármán number $h^{+}$ is established for the integrated mean scalar in the range $h^{+}\geqslant 400$ where the mean part of the total scalar dissipation exhibits near constancy, whilst the integral of the turbulent scalar dissipation rate $\overline{\unicode[STIX]{x1D700}_{\unicode[STIX]{x1D703}}}$ increases logarithmically with $h^{+}$. This logarithmic dependence is similar to that established in a previous paper (Abe & Antonia, J. Fluid Mech., vol. 798, 2016, pp. 140–164) for the bulk mean velocity. However, the slope (2.18) for the integrated mean scalar is smaller than that (2.54) for the bulk mean velocity. The ratio of these two slopes is 0.85, which can be identified with the value of the turbulent Prandtl number in the overlap region. It is shown that the logarithmic $h^{+}$ increase of the integrated mean scalar is intrinsically associated with the overlap region of $\overline{\unicode[STIX]{x1D700}_{\unicode[STIX]{x1D703}}}$, established for $h^{+}$ (${\geqslant}400$). The resulting heat transfer law also holds at a smaller $h^{+}$ (${\geqslant}200$) than that derived by assuming a log law for the mean temperature.

NUMERICAL RESEARCH ON THE EFFECT OF FIBERS ON THE PROPERTY OF TURBULENT FIBER SUSPENSION IN A CHANNEL

2009 ◽  
Vol 23 (03) ◽  
pp. 509-512 ◽  
Author(s):  
SUHUA SHEN ◽  
JIANZHONG LIN
Keyword(s):  
Reynolds Number ◽  
Kinetic Energy ◽  
Turbulent Kinetic Energy ◽  
Channel Flow ◽  
Turbulent Channel Flow ◽  
Fiber Suspension ◽  
Additional Stress ◽  
Mean Velocity ◽  
Navier Stokes Equation ◽  
Fiber Concentration

To explore the rheological property in turbulent channel flow of fiber suspensions, the equation of probability distribution function for mean fiber orientation and the Reynolds averaged Navier-Stokes equation with the term of additional stress resulted from fibers were solved with numerical methods to get the distributions of the mean velocity and turbulent kinetic energy. The simulation results show that the effect of fibers on turbulent channel flow is equivalent to an additional viscosity. The turbulent velocity profiles of fiber suspension become gradually sharper by increasing the fiber concentration and/or decreasing the Reynolds number. The turbulent kinetic energy will increase with increasing Reynolds number and fiber concentration.

scholarly journals Resolvent analysis of turbulent channel flow with manipulated mean velocity profile

2020 ◽  
Vol 15 (3) ◽  
pp. JFST0014-JFST0014
Author(s):  
Riko UEKUSA ◽  
Aika KAWAGOE ◽  
Yusuke NABAE ◽  
Koji FUKAGATA
Keyword(s):  
Velocity Profile ◽  
Channel Flow ◽  
Turbulent Channel Flow ◽  
Mean Velocity

Roll-cell instabilities in rotating laminar and trubulent channel flows

1976 ◽  
Vol 77 (1) ◽  
pp. 153-174 ◽  
Author(s):  
Dietrich K. Lezius ◽  
James P. Johnston
Keyword(s):  
Channel Flow ◽  
Poiseuille Flow ◽  
Equations Of Motion ◽  
Turbulent Channel Flow ◽  
Reynolds Numbers ◽  
Marginal Stability ◽  
Mean Flow ◽  
Unstable Flow ◽  
Mean Velocity ◽  
Longitudinal Instabilities

The stability of laminar and turbulent channel flow is examined for cases where Coriolis forces are introduced by steady rotation about an axis perpendicular to the plane of mean flow. Linearized equations of motion are derived for small disturbances of the Taylor type. Conditions for marginal stability in laminar Couette and Poiseuille flow correspond, in part, to the analogous solutions of buoyancy-driven convection instabilities in heated fluid layers, and to those of Taylor instabilities in the flow between rotating cylinders. In plane Poiseuille flow with rotation, the critical disturbance mode occurs at a Reynolds number of Rec= 88.53 and rotation number Ro= 0.5. At higher Reynolds numbers, unstable conditions canexist over the range of rotation numbers given by 0 < Ro< 3, provided the undisturbed flow remains laminar. A two-layer model is devised to investigate the onset of longitudinal instabilities in turbulent flow. The linear disturbance equations are solved essentially in their laminar form, whereby the velocity gradient of laminar flow is replaced by a numerically computed profile for the gradient of the turbulent mean velocity. The turbulent stress levels in the stable and unstable flow regions are represented by integrated averages of the eddy viscosity. Onset of instability for Reynolds numbers between 6000 and 35 000 is predicted to occur at Ro= 0.022, a value in remarkable agreement with the experimentally observed appearance of roll instabilities in rotating turbulent channel flow.

Shear-Improved Smagorinsky Modeling of Turbulent Channel Flow Using Generalized Lattice Boltzmann Equation

2010 ◽  
Author(s):  
Saeed Jafari ◽  
Mohammad Rahnama
Keyword(s):  
Boltzmann Equation ◽  
Channel Flow ◽  
Lattice Boltzmann ◽  
Computational Cost ◽  
Turbulent Channel Flow ◽  
Lattice Boltzmann Equation ◽  
Wall Region ◽  
Smagorinsky Model ◽  
Mean Velocity ◽  
Turbulent Statistics

Generalized Lattice Boltzmann Equation (GLBE) was used for computation of turbulent channel flow for which Large Eddy Simulation (LES) was employed as a turbulence model. The subgrid-Scale turbulence effects were simulated through a Shear-Improved Smagorinsky Model (SISM) which is capable of predicting turbulent near wall region accurately without any wall function. Computations were done for a relatively coarse grid with shear Reynolds number of 180 in a parallelized code. Good numerical stability was observed for this computational framework. Results of mean velocity distribution across the channel showed good correspondence with Direct Numerical Simulation (DNS) data. Negligible discrepancies were observed for computed turbulent statistics between present computations and those reported from DNS. Three-dimensional instantaneous vorticity contours showed complex vortical structures appeared in such flow geometries. It is concluded that such framework is capable of predicting accurate results for turbulent channel flow without adding significant complication and computational cost to the standard Smagorinsky model.

Sign in / Sign up

Export Citation Format

Share Document