scholarly journals Wall-bounded flow over a realistically rough superhydrophobic surface

2019 ◽  
Vol 873 ◽  
pp. 977-1019 ◽  
Author(s):  
Karim Alamé ◽  
Krishnan Mahesh

Direct numerical simulation (DNS) is performed for two wall-bounded flow configurations: laminar Couette flow at $Re=740$ and turbulent channel flow at $Re_{\unicode[STIX]{x1D70F}}=180$, where $\unicode[STIX]{x1D70F}$ is the shear stress at the wall. The top wall is smooth and the bottom wall is a realistically rough superhydrophobic surface (SHS), generated from a three-dimensional surface profile measurement. The air–water interface, which is assumed to be flat, is simulated using the volume-of-fluid (VOF) approach. The two flow cases are studied with varying interface heights $h$ to understand its effect on slip and drag reduction ($DR$). For the laminar Couette flow case, the presence of the surface roughness is felt up to $40\,\%$ of the channel height in the wall-normal direction. Nonlinear dependence of $DR$ on $h$ is observed with three distinct regions. A nonlinear curve fit is obtained for gas fraction $\unicode[STIX]{x1D719}_{g}$ as a function of $h$, where $\unicode[STIX]{x1D719}_{g}$ determines the amount of slip area exposed to the flow. A power law fit is obtained from the data for the effective slip length as a function of $\unicode[STIX]{x1D719}_{g}$ and is compared to those derived for structured geometry. For the turbulent channel flow, statistics of the flow field are compared to that of a smooth wall to understand the effects of roughness and $h$. Four cases are simulated ranging from fully wetted to fully covered and two intermediate regions in between. Scaling laws for slip length, slip velocity, roughness function and $DR$ are obtained for different penetration depths and are compared to past work for structured geometry. $DR$ is shown to depend on a competing effect between slip velocity and turbulent losses due to the Reynolds shear stress contribution. Presence of trapped air in the cavities significantly alters near-wall flow physics where we examine near-wall structures and propose a physical mechanism for their behaviour. The fully wetted roughness increases the peak value of turbulent intensities, whereas the presence of the interface suppresses them. The pressure fluctuations have competing contributions between turbulent pressure fluctuations and stagnation due to asperities, the near-wall structure is altered and breaks down with increasing slip. Overall, there exists a competing effect between the interface and the asperities, the interface suppresses turbulence whereas the asperities enhance them. The present work demonstrates DNS over a realistic multiphase SHS for the first time, to the best of our knowledge.

Optics ◽  
2020 ◽  
Vol 1 (1) ◽  
pp. 40-51
Author(s):  
Esther Mäteling ◽  
Michael Klaas ◽  
Wolfgang Schröder

An extended experimental method is presented in which the micro-pillar shear-stress sensor (MPS 3 ) and high-speed stereo particle-image velocimetry measurements are simultaneously performed in turbulent channel flow to conduct concurrent time-resolved measurements of the two-dimensional wall-shear stress (WSS) distribution and the velocity field in the outer flow. The extended experimental setup, which involves a modified MPS 3 measurement setup and data evaluation compared to the standard method, is presented and used to investigate the footprint of the outer, large-scale motions (LSM) onto the near-wall small-scale motions. The measurements were performed in a fully developed, turbulent channel flow at a friction Reynolds number R e τ = 969 . A separation between large and small scales of the velocity fluctuations and the WSS fluctuations was performed by two-dimensional empirical mode decomposition. A subsequent cross-correlation analysis between the large-scale velocity fluctuations and the large-scale WSS fluctuations shows that the streamwise inclination angle between the LSM in the outer layer and the large-scale footprint imposed onto the near-wall dynamics has a mean value of Θ ¯ x = 16.53 ∘ , which is consistent with the literature relying on direct numerical simulations and hot-wire anemometry data. When also considering the spatial shift in the spanwise direction, the mean inclination angle reduces to Θ ¯ x z = 13.92 ∘ .


Recent data obtained in an unsteady turbulent channel flow is reviewed. Results concerning the modulation characteristics of the Reynolds shear stresses, of the structural parameters and of the length scales inferred from unsteady spatial correlations are discussed. The close examination of both the amplitude and the phase shifts of the Reynolds shear stresses confirms the existence of three distinct inposed frequency regimes, namely the quasi-steady regime, the relaxation regime, in which the amplitudes decrease and which is accompanied by large time lags, and a subsequent third regime wherein the modulation characteristics change considerably. The fine structure of the near-wall turbulence response through quadrant analysis reveals large cyclic variations of the contributions of ejections and sweeps to the Reynolds shear stress. The reaction of the spanwise extent of the near-wall structures is investigated through the spanwise correlation coefficient between the wall shear stress and the streamwise velocity, and the resulting length scales. A temporal filtering of both signals shows that the inactive motions respond uniformly in the whole imposed frequency regime. A strong correlation is found between the modulation characteristics of the streak spacing and the ejection frequency.


2014 ◽  
Vol 756 ◽  
pp. 30-53 ◽  
Author(s):  
Euiyoung Kim ◽  
Haecheon Choi

AbstractThe space–time characteristics of a compliant wall in a turbulent channel flow are investigated using direct numerical simulation (DNS). The compliant wall is modelled as a homogeneous plane supported by spring-and-damper arrays and is passively driven by wall-pressure fluctuations. The frequency/wavenumber spectra and convection velocities of the wall-pressure fluctuations, wall displacement and wall velocity are obtained from the present simulation. As the spring, damping, and tension coefficients decrease, the wall becomes softer and the wall displacement and velocity fluctuations increase. For a relatively stiff compliant wall (i.e. large spring, damping and streamwise tension coefficients), there are few changes in the skin-friction drag and near-wall turbulence structures. However, when a compliant wall is soft (i.e. small spring, damping and streamwise tension coefficients), the wall moves in the form of a large-amplitude quasi-two-dimensional wave travelling in the downstream direction. This wave is generated by the resonance of the wall property and the near-wall flow is significantly activated by this wall motion. The power spectra of wall variables show distinct peaks near the resonance frequencies. The convection velocities of the wall motion and wall-pressure fluctuations become smaller with a softer wall.


2009 ◽  
Vol 631 ◽  
pp. 281-309 ◽  
Author(s):  
D. CHUNG ◽  
D. I. PULLIN

We report large-eddy simulation (LES) of turbulent channel flow. This LES neither resolves nor partially resolves the near-wall region. Instead, we develop a special near-wall subgrid-scale (SGS) model based on wall-parallel filtering and wall-normal averaging of the streamwise momentum equation, with an assumption of local inner scaling used to reduce the unsteady term. This gives an ordinary differential equation (ODE) for the wall shear stress at every wall location that is coupled with the LES. An extended form of the stretched-vortex SGS model, which incorporates the production of near-wall Reynolds shear stress due to the winding of streamwise momentum by near-wall attached SGS vortices, then provides a log relation for the streamwise velocity at the top boundary of the near-wall averaged domain. This allows calculation of an instantaneous slip velocity that is then used as a ‘virtual-wall’ boundary condition for the LES. A Kármán-like constant is calculated dynamically as part of the LES. With this closure we perform LES of turbulent channel flow for Reynolds numbers Reτ based on the friction velocity uτ and the channel half-width δ in the range 2 × 103 to 2 × 107. Results, including SGS-extended longitudinal spectra, compare favourably with the direct numerical simulation (DNS) data of Hoyas & Jiménez (2006) at Reτ = 2003 and maintain an O(1) grid dependence on Reτ.


2016 ◽  
Vol 790 ◽  
Author(s):  
Taeyong Jung ◽  
Haecheon Choi ◽  
John Kim

The anisotropy of the slip length and its effect on the skin-friction drag are numerically investigated for a turbulent channel flow with an idealized superhydrophobic surface having an air layer, where the idealized air–water interface is flat and does not contain the surface-tension effect. Inside the air layer, both the shear-driven flow and recirculating flow with zero net mass flow rate are considered. With increasing air-layer thickness, the slip length, slip velocity and percentage of drag reduction increase. It is shown that the slip length is independent of the water flow and depends only on the air-layer geometry. The amount of drag reduction obtained is in between those by the empirical formulae from the streamwise slip only and isotropic slip, indicating that the present air–water interface generates an anisotropic slip, and the streamwise slip length ($b_{x}$) is larger than the spanwise one ($b_{z}$). From the joint probability density function of the slip velocities and velocity gradients at the interface, we confirm the anisotropy of the slip lengths and obtain their relative magnitude ($b_{x}/b_{z}=4$) for the present idealized superhydrophobic surface. It is also shown that the Navier slip model is valid only in the mean sense, and it is generally not applicable to fluctuating quantities.


Author(s):  
Derek C. Tretheway ◽  
Luoding Zhu ◽  
Linda Petzold ◽  
Carl D. Meinhart

This work examines the slip boundary condition by Lattice Boltzmann simulations, addresses the validity of the Navier’s hypothesis that the slip velocity is proportional to the shear rate and compares the Lattice Boltzmann simulations to the experimental results of Tretheway and Meinhart (Phys. of Fluids, 14, L9–L12). The numerical simulation models the boundary condition as the probability, P, of a particle to bounce-back relative to the probability of specular reflection, 1−P. For channel flow, the numerically calculated velocity profiles are consistent with the experimental profiles for both the no-slip and slip cases. No-slip is obtained for a probability of 100% bounce-back, while a probability of 0.03 is required to generate a slip length and slip velocity consistent with the experimental results of Tretheway and Meinhart for a hydrophobic surface. The simulations indicate that for microchannel flow the slip length is nearly constant along the channel walls, while the slip velocity varies with wall position as a results of variations in shear rate. Thus, the resulting velocity profile in a channel flow is more complex than a simple combination of the no-slip solution and slip velocity as is the case for flow between two infinite parallel plates.


2018 ◽  
Vol 15 (2) ◽  
pp. 75-89
Author(s):  
Muhammad Saiful Islam Mallik ◽  
Md. Ashraf Uddin

A large eddy simulation (LES) of a plane turbulent channel flow is performed at a Reynolds number Re? = 590 based on the channel half width, ? and wall shear velocity, u? by approximating the near wall region using differential equation wall model (DEWM). The simulation is performed in a computational domain of 2?? x 2? x ??. The computational domain is discretized by staggered grid system with 32 x 30 x 32 grid points. In this domain the governing equations of LES are discretized spatially by second order finite difference formulation, and for temporal discretization the third order low-storage Runge-Kutta method is used. Essential turbulence statistics of the computed flow field based on this LES approach are calculated and compared with the available Direct Numerical Simulation (DNS) and LES data where no wall model was used. Comparing the results throughout the calculation domain we have found that the LES results based on DEWM show closer agreement with the DNS data, especially at the near wall region. That is, the LES approach based on DEWM can capture the effects of near wall structures more accurately. Flow structures in the computed flow field in the 3D turbulent channel have also been discussed and compared with LES data using no wall model.


2019 ◽  
Vol 4 (7) ◽  
Author(s):  
Pierre-Alain Gubian ◽  
Jordan Stoker ◽  
James Medvescek ◽  
Laurent Mydlarski ◽  
B. Rabi Baliga

Author(s):  
Boris Arcen ◽  
Anne Tanie`re ◽  
Benoiˆt Oesterle´

The importance of using the lift force and wall-corrections of the drag coefficient for modeling the motion of solid particles in a fully-developed channel flow is investigated by means of direct numerical simulation (DNS). The turbulent channel flow is computed at a Reynolds number based on the wall-shear velocity and channel half-width of 185. Contrary to most of the numerical simulations, we consider in the present study a lift force formulation that accounts for the weak and strong shear as well as for the wall effects (hereinafter referred to as optimum lift force), and the wall-corrections of the drag force. The DNS results show that the optimum lift force and the wall-corrections of the drag together have little influence on most of the statistics (particle concentration, mean velocities, and mean relative and drift velocities), even in the near wall region.


Sign in / Sign up

Export Citation Format

Share Document