A numerical investigation of acceleration-skewed oscillatory flows

2016 ◽  
Vol 808 ◽  
pp. 576-613 ◽  
Author(s):  
Pietro Scandura ◽  
Carla Faraci ◽  
Enrico Foti
Keyword(s):  
Reynolds Number ◽  
Shear Stress ◽  
Wall Shear Stress ◽  
Relative Intensity ◽  
Half Cycle ◽  
Bed Shear Stress ◽  
Streamwise Direction ◽  
Oscillatory Flows ◽  
The Asymmetry ◽  
Wall Shear

Numerical simulations of wall-bounded acceleration-skewed oscillatory flows are here presented. The relevance of this type of boundary layer arises in connection with coastal hydrodynamics and sediment transport, as it is generated at the bottom of sea waves in shallow water. Because of the acceleration skewness, the bed shear stress during the onshore half-cycle is larger than in the offshore half-cycle. The asymmetry in the bed shear stress increases with increasing acceleration skewness, while an increase of the Reynolds number from the laminar regime causes the asymmetry first to decrease and then increase. Low- and high-speed streaks of fluid elongated in the streamwise direction emerge near the wall, shortly after the beginning of each half-cycle, at a phase that depends on the flow parameters. Such flow structures strengthen during the first part of the accelerating phase, without causing a significant deviation of the streamwise wall shear stress from the laminar values. Before the occurrence of the peak of the free stream velocity, the low-speed streaks break down into small turbulent structures causing a large increase in wall shear stress. The ratio of the root-mean-square (r.m.s.) of the fluctuations to the mean value (relative intensity) of the wall shear stress is approximately 0.4 throughout a relatively wide interval of the flow cycle that begins when breaking down of the streaks has occurred in the entire fluid domain. The acceleration skewness and the Reynolds number determine the phase at which this time interval begins. Both the skewness and the flatness coefficients of the streamwise wall shear stress are large when elongated streaks are present, while values of approximately 1.1 and 5.4 respectively occur just after breaking has occurred. The trend of both the relative intensity and the flatness of the spanwise wall shear stress are qualitatively similar to those of the wall shear in the streamwise direction. As a result of the acceleration skewness, the period-averaged Reynolds stress does not vanish. Consequently, an offshore directed steady streaming is generated which persists into the irrotational region.

Flow and Thermal Fields in a Pendant Droplet Moving on Lyophobic Surface

2010 ◽  
Author(s):  
Basant Singh Sikarwar ◽  
K. Muralidhar ◽  
Sameer Khandekar
Keyword(s):  
Heat Transfer ◽  
Contact Angle ◽  
Reynolds Number ◽  
Shear Stress ◽  
Heat Transfer Coefficient ◽  
Wall Shear Stress ◽  
Transfer Coefficient ◽  
Maximum Shear Stress ◽  
Computational Domain ◽  
Wall Shear

Clusters of liquid drops growing and moving on physically or chemically textured lyophobic surfaces are encountered in drop-wise mode of vapor condensation. As opposed to film-wise condensation, drops permit a large heat transfer coefficient and are hence attractive. However, the temporal sustainability of drop formation on a surface is a challenging task, primarily because the sliding drops eventually leach away the lyophobicity promoter layer. Assuming that there is no chemical reaction between the promoter and the condensing liquid, the wall shear stress (viscous resistance) is the prime parameter for controlling physical leaching. The dynamic shape of individual droplets, as they form and roll/slide on such surfaces, determines the effective shear interaction at the wall. Given a shear stress distribution of an individual droplet, the net effect of droplet ensemble can be determined using the time averaged population density during condensation. In this paper, we solve the Navier-Stokes and the energy equation in three-dimensions on an unstructured tetrahedral grid representing the computational domain corresponding to an isolated pendant droplet sliding on a lyophobic substrate. We correlate the droplet Reynolds number (Re = 10–500, based on droplet hydraulic diameter), contact angle and shape of droplet with wall shear stress and heat transfer coefficient. The simulations presented here are for Prandtl Number (Pr) = 5.8. We see that, both Poiseuille number (Po) and Nusselt number (Nu), increase with increasing the droplet Reynolds number. The maximum shear stress as well as heat transfer occurs at the droplet corners. For a given droplet volume, increasing contact angle decreases the transport coefficients.

A Dynamic Procedure for Wall-Modeled Large-Eddy Simulation of High Reynolds Number Flows

2011 ◽  
Author(s):  
Soshi Kawai
Keyword(s):  
Reynolds Number ◽  
Shear Stress ◽  
Wall Shear Stress ◽  
Large Eddy Simulation ◽  
High Reynolds Number ◽  
Eddy Simulation ◽  
Large Eddy ◽  
Wall Shear ◽  
High Reynolds ◽  
Near Wall

This paper addresses the error in large-eddy simulation with wall-modeling (i.e., when the wall shear stress is modeled and the viscous near-wall layer is not resolved): the error in estimating the wall shear stress from a given outer-layer velocity field using auxiliary near-wall RANS equations where convection is not neglected. By considering the behavior of turbulence length scales near a wall, the cause of the errors is diagnosed and solutions that remove the errors are proposed based solidly on physical reasoning. The resulting method is shown to accurately predict equilibrium boundary layers at very high Reynolds number, with both realistic instantaneous fields (without overly elongated unphysical near-wall structures) and accurate statistics (both skin friction and turbulence quantities).

Evolution of wall shear stress with Reynolds number in fully developed turbulent channel flow experiments

2019 ◽  
Vol 4 (7) ◽  
Author(s):  
Pierre-Alain Gubian ◽  
Jordan Stoker ◽  
James Medvescek ◽  
Laurent Mydlarski ◽  
B. Rabi Baliga
Keyword(s):  
Reynolds Number ◽  
Shear Stress ◽  
Wall Shear Stress ◽  
Channel Flow ◽  
Turbulent Channel Flow ◽  
Wall Shear ◽  
Flow Experiments

scholarly journals Physiological non-Newtonian blood flow through single stenosed artery

2016 ◽  
Vol 43 (1) ◽  
pp. 99-115 ◽  
Author(s):  
Khairuzzaman Mamun ◽  
Most. Akhter ◽  
Mohammad Ali
Keyword(s):  
Reynolds Number ◽  
Shear Stress ◽  
Wall Shear Stress ◽  
Three Dimensional ◽  
Cross Sectional ◽  
Number Range ◽  
Newtonian Model ◽  
Early Systole ◽  
Wall Shear ◽  
Downstream Flow

A numerical simulation to investigate the Non-Newtonian modeling effects on physiological flows in a three dimensional idealized artery with a single stenosis of 85% severity is given. The wall vessel is considered to be rigid. Oscillatory physiological and parabolic velocity profile has been imposed for inlet boundary condition. Determination of the physiological waveform is performed using a Fourier series with sixteen harmonics. The investigation has a Reynolds number range of 96 to 800. Low Reynolds number k ? w model is used as governing equation. The investigation has been carried out to characterize two Non-Newtonian constitutive equations of blood, namely, (i) Carreau and (ii) Cross models. The Newtonian model has also been investigated to study the physics of fluid. The results of Newtonian model are compared with the Non-Newtonian models. The numerical results are presented in terms of velocity, pressure, wall shear stress distributions and cross sectional velocities as well as the streamlines contour. At early systole pressure differences between Newtonian and Non-Newtonian models are observed at pre-stenotic, throat and immediately after throat regions. In the case of wall shear stress, some differences between Newtonian and Non-Newtonian models are observed when the flows are minimum such as at early systole or diastole. In general, the velocities at throat regions are highest at all-time phase. However, at pick systole higher velocities are observed at post-stenotic region. Downstream flow of all models creates some recirculation regions at diastole.

Very Large-Scale Structures and Their Effects on the Wall Shear-Stress Fluctuations in a Turbulent Channel Flow up to Reτ=640

2004 ◽  
Vol 126 (5) ◽  
pp. 835-843 ◽  
Author(s):  
Hiroyuki Abe ◽  
Hiroshi Kawamura ◽  
Haecheon Choi
Keyword(s):  
Reynolds Number ◽  
Shear Stress ◽  
Wall Shear Stress ◽  
High Reynolds Number ◽  
Large Scale ◽  
Turbulent Channel Flow ◽  
Large Scale Structures ◽  
Scale Structures ◽  
Wall Shear ◽  
High Reynolds

Direct numerical simulation of a fully developed turbulent channel flow has been carried out at three Reynolds numbers, 180, 395, and 640, based on the friction velocity and the channel half width, in order to investigate very large-scale structures and their effects on the wall shear-stress fluctuations. It is shown that very large-scale structures exist in the outer layer and that they certainly contribute to inner layer structures at high Reynolds number. Moreover, it is revealed that very large-scale structures exist even in the wall shear-stress fluctuations at high Reynolds number, which are essentially associated with the very large-scale structures in the outer layer.

Characterization of a Sudden Expansion Flow Chamber to Study the Response of Endothelium to Flow Recirculation

1995 ◽  
Vol 117 (2) ◽  
pp. 203-210 ◽  
Author(s):  
George A. Truskey ◽  
Kevin M. Barber ◽  
Thomas C. Robey ◽  
Lauri A. Olivier ◽  
Marty P. Combs
Keyword(s):  
Reynolds Number ◽  
Shear Stress ◽  
Wall Shear Stress ◽  
Recirculation Zone ◽  
Numerical Solutions ◽  
Flow Chamber ◽  
Sudden Expansion ◽  
Shear Stresses ◽  
Stress Gradients ◽  
Wall Shear

In order to simulate regions of flow separation observed in vivo, a conventional parallel plate flow chamber was modified to produce an asymmetric sudden expansion. The flow field was visualized using light reflecting particles and the size of the recirculation zone was measured by image analysis of the particles. Finite element numerical solutions of the two and three-dimensional forms of the Navier-Stokes equation were used to determine the wall shear stress distribution and predict the location of reattachment. For two different size expansions, numerical estimates of the reattachment point along the centerline of the flow chamber agreed well with experimental values for Reynolds numbers below 473. Even at a Reynolds number of 473, the flow could be approximated as two-dimensional for 80 percent of the chamber width. Peak shear stresses in the recirculation zone as high as 80 dyne/cm2 and shear stress gradients of 2500 (dyne/cm2)/cm were produced. As an application of this flow chamber, subconfluent bovine aortic endothelial cell shape and orientation were examined in the zone of recirculation during a 24 h exposure to flow at a Reynolds number of 267. After 24 h, gradients in cell orientation and shape were observed within the recirculation zone. At the location of reattachment, where the wall shear stress was zero but the shear stress gradients were large, cells plated at low density were still aligned with the direction of flow. No preferred orientation was observed at the gasket edge where the wall shear stress and shear stress gradients were zero. At higher cell densities, no alignment was observed at the separation point. The results suggest that endothelial cells can respond to spatial gradients of wall shear stress.

Turbulent dynamics of sinusoidal oscillatory flow over a wavy bottom

2018 ◽  
Vol 858 ◽  
pp. 264-314 ◽  
Author(s):  
Asim Önder ◽  
Jing Yuan
Keyword(s):  
Shear Stress ◽  
Wall Shear Stress ◽  
Oscillatory Flow ◽  
Free Stream ◽  
Stream Velocity ◽  
Half Cycle ◽  
Primary Vortex ◽  
Vortex Dipole ◽  
Coherent Vortices ◽  
Wall Shear

A direct numerical simulation study is conducted to investigate sinusoidal oscillatory flow over a two-dimensional wavy wall. The height and wavelength of the bottom profile, and the period and amplitude of the free-stream oscillation, are selected to mimic a wave-driven boundary layer over vortex ripples on a sandy seabed. Two cases with different Reynolds numbers$(Re)$are considered, and the higher-$Re$case achieves a fully developed turbulent state with a wide separation between the energy-containing and dissipative scales. The oscillatory flow is characterized by coherent columnar vortices, which are the main transport agents of turbulent kinetic energy and enstrophy. Two classes of coherent vortices are observed: (i) a primary vortex formed at the lee side of the ripple by flow separation at the crest; (ii) a secondary vortex formed beneath the primary vortex by vortex-induced separation. When the free-stream velocity weakens, these vortices form a counter-rotating vortex dipole and eject themselves over the crest with their mutual induction. Turbulence production peaks twice in a half-cycle; during the formation of the primary vortex and during the ejection of the vortex dipole. The intensity of the former peak remains low in the lower-$Re$case, as the vortex dipole follows a higher altitude trajectory limiting its interactions with the bottom, and leaving minimal residual turbulence around the ripples for the subsequent half-cycle. Flow snapshots and spectral analysis reveal two dominant three-dimensional features: (i) an energetic vortex mode with a preferred spanwise wavelength close to the ripple wavelength; (ii) streamwise vortical structures in near-wall regions with a relatively shorter spanwise spacing influenced by viscous effects. The vortex mode becomes strong when the cores of the vortices are strained to an elliptical form while moving towards the crest. Following the detachment of the vortices from the ripple, the vortex mode in the higher-$Re$case breaks down the spanwise coherence of the columnar vortices and decomposes them into intermittent patches of turbulent vortex clusters. The distribution of wall shear stress over the ripple is also analysed in detail. The peak values are observed near the ripple crest around the ejection of the vortex dipole and the maximum free-stream velocity. In the former, both the vortex mode and streamwise vortices have strong footprints on the wall, yielding a bimodal wall-shear-stress spectrum with two distinctive peaks. In the second high-stress regime, decaying coherent vortices impose strong inhomogeneity on the wall shear stress as their wall-attached parts sweep the ripples. These spanwise variations in the wall shear provide insights into the instability of two-dimensional sand ripples.

Sign in / Sign up

Export Citation Format

Share Document