Prediction of Turbulent Rough-Wall Skin Friction Using a Discrete Element Approach

1985 ◽  
Vol 107 (2) ◽  
pp. 251-257 ◽  
Author(s):  
R. P. Taylor ◽  
H. W. Coleman ◽  
B. K. Hodge

A discrete element model for turbulent flow over rough surfaces has been derived from basic principles. This formulation includes surface roughness form drag and blockage effects as a constituent part of the partial differential equations and does not rely on a single-length-scale concept such as equivalent sandgrain roughness. The roughness model includes the necessary empirical information on the interaction between three-dimensional roughness elements and the flow in a general way which does not require experimental data on each specific surface. This empirical input was determined using data from well-accepted experiments. Predictions using the model are compared with additional data for fully-developed and boundary layer flows. The predictions are shown to compare equally well with both transitionally rough and fully rough turbulent flows without modification of the roughness model.

2015 ◽  
Vol 9 (2) ◽  
pp. 1425-1460 ◽  
Author(s):  
P. Hagenmuller ◽  
G. Chambon ◽  
M. Naaim

Abstract. Rapid and large deformations of snow are mainly controlled by grain rearrangements, which occur through the failure of cohesive bonds and the creation of new contacts. We exploit a granular description of snow to develop a discrete element model based on the full three-dimensional microstructure captured by microtomography. The model assumes that snow is composed of rigid grains interacting through localized contacts accounting for cohesion and friction. The geometry of the grains and of the intergranular bonding system are explicitly defined from microtomographic data using geometrical criteria based on curvature and contiguity. Single grains are represented as rigid clumps of spheres. The model is applied to different snow samples subjected to confined compression tests. A detailed sensitivity analysis shows that artifacts introduced by the modeling approach and the influence of numerical parameters are limited compared to variations due to the geometry of the microstructure. The model shows that the compression behavior of snow is mainly controlled by the density of the samples, but that deviations from a pure density parameterization are not insignificant during the first phase of deformation. In particular, the model correctly predicts that, for a given density, faceted crystals are less resistant to compression than rounded grains or decomposed snow. For larger compression strains, no clear differences between snow types are observed.


2009 ◽  
Vol 131 (2) ◽  
Author(s):  
Stephen T. McClain ◽  
Jason M. Brown

The discrete-element model for flows over rough surfaces was recently modified to predict drag and heat transfer for flow over randomly rough surfaces. However, the current form of the discrete-element model requires a blockage fraction and a roughness-element diameter distribution as a function of height to predict the drag and heat transfer of flow over a randomly rough surface. The requirement for a roughness-element diameter distribution at each height from the reference elevation has hindered the usefulness of the discrete-element model and inhibited its incorporation into a computational fluid dynamics (CFD) solver. To incorporate the discrete-element model into a CFD solver and to enable the discrete-element model to become a more useful engineering tool, the randomly rough surface characterization must be simplified. Methods for determining characteristic diameters for drag and heat transfer using complete three-dimensional surface measurements are presented. Drag and heat transfer predictions made using the model simplifications are compared to predictions made using the complete surface characterization and to experimental measurements for two randomly rough surfaces. Methods to use statistical surface information, as opposed to the complete three-dimensional surface measurements, to evaluate the characteristic dimensions of the roughness are also explored.


1991 ◽  
Vol 113 (1) ◽  
pp. 34-41 ◽  
Author(s):  
G. J. Yoo ◽  
R. M. C. So ◽  
B. C. Hwang

Internal rotating boundary-layer flows are strongly influenced by large circumferential strain and the turbulence field is anisotropic. This is especially true in the entry region of a rotating pipe where the flow is three dimensional, the centrifugal force due to fluid rotation is less important, and the circumferential strain created by surface rotation has a significant effect on the turbulence field near the wall. Consequently, viscous effects cannot be neglected in the near-wall region. Several low-Reynolds-number turbulence closures are proposed for the calculation of developing rotating pipe flows. Some are two-equation closures with and without algebraic stress correction, while others are full Reynolds-stress closures. It is found that two-equation closures with and without algebraic stress correction are totally inadequate for this three-dimensional flow, while Reynolds-stress closures give results that are in good agreement with measurements over a wide range of rotation numbers.


Author(s):  
Stephen T. McClain ◽  
Jason M. Brown

The discrete-element model for flows over rough surfaces was recently modified to predict drag and heat transfer for flow over randomly-rough surfaces. However, the current form of the discrete-element model requires a blockage fraction and a roughness-element diameter distribution as a function of height to predict the drag and heat transfer of flow over a randomly-rough surface. The requirement for a roughness element-diameter distribution at each height from the reference elevation has hindered the usefulness of the discrete-element model and inhibited its incorporation into a computational fluid dynamics (CFD) solver. To incorporate the discrete-element model into a CFD solver and to enable the discrete-element model to become a more useful engineering tool, the randomly-rough surface characterization must be simplified. Methods for determining characteristic diameters for drag and heat transfer using complete three-dimensional surface measurements are presented. Drag and heat transfer predictions made using the model simplifications are compared to predictions made using the complete surface characterization and to experimental measurements for two randomly-rough surfaces. Methods to use statistical surface information, as opposed to the complete three-dimensional surface measurements, to evaluate the characteristic dimensions of the roughness are also explored.


2013 ◽  
Vol 135 (12) ◽  
Author(s):  
S. Cherubini ◽  
M. D. de Tullio ◽  
P. De Palma ◽  
G. Pascazio

This work provides a three-dimensional energy optimization analysis, looking for perturbations inducing the largest energy growth at a finite time in a boundary-layer flow in the presence of roughness elements. The immersed boundary technique has been coupled with a Lagrangian optimization in a three-dimensional framework. Four roughness elements with different heights have been studied, inducing amplification mechanisms that bypass the asymptotical growth of Tollmien–Schlichting waves. The results show that even very small roughness elements, inducing only a weak deformation of the base flow, can strongly localize the optimal disturbance. Moreover, the highest value of the energy gain is obtained for a varicose perturbation. This result demonstrates the relevance of varicose instabilities for such a flow and shows a different behavior with respect to the secondary instability theory of boundary layer streaks.


2021 ◽  
Vol 2127 (1) ◽  
pp. 012006
Author(s):  
M V Philippov ◽  
A V Zolotukhin ◽  
I A Chokhar ◽  
V V Terekhov ◽  
V I Terekhov ◽  
...  

Abstract The work is aimed at creating a technique for reconstructing the three-dimensional structure of turbulent flows using data of a two-component laser-Doppler anemometer at the example of a flow in flat coplanar channels. Data of paired measurements by a two-component LDA at two different angles relative to the measurement point are used to obtain information about the three-dimensional structure of the flow. Further, all three components of the velocity vector were calculated using transformations. The developed method is used to measure the velocity field in a coplanar channel.


Author(s):  
S. Cherubini ◽  
M. D. de Tullio ◽  
P. De Palma ◽  
G. Pascazio

This work provides a three-dimensional energy optimization analysis, looking for perturbations inducing the largest energy growth at a finite time in a boundary-layer flow in the presence of roughness elements. Amplification mechanisms are described which by-pass the asymptotical growth of Tollmien–Schlichting waves. The immersed boundary technique has been coupled with a Lagrangian optimization in a three-dimensional framework. Two types of roughness elements have been studied, characterized by a different height. The results show that even very small roughness elements, inducing only a weak deformation of the base flow, can strongly localize the optimal disturbance. Moreover, the highest value of the energy gain is obtained for a varicose perturbation, pointing out the importance of varicose instabilities for such a flow and a different behavior with respect to the secondary instability theory of boundary layer streaks.


2022 ◽  
Vol 9 ◽  
Author(s):  
Kai Wu ◽  
Zan Li ◽  
Zhibin Liu ◽  
Songyu Liu

This work provides a three-dimensional discrete element simulation (DEM) model to study the air sparging technology. The simulations have taken into account the multi-phases of bubble (gas) - fluid (water) - soil (solid) particles. Bubbles are treated as discrete individual particles, with buoyancy and drag forces applied to bubbles and soil particles. The trajectory of each discrete bubble particle can be tracked using the discrete element model. It is found that the diffusion of the whole bubble is inverted conical though the motion behavior of a single bubble particle is random. Furthermore, the distribution of the radius of influence (ROI) is not uniform. The bubbles become more concentrated as in the center of the inverted cone. The number of bubbles dissipated from the water surface is normally distributed. The DEM simulation is a novel approach to studying air sparging technology that can provide us a deeper insight into bubble migration at the microscopic level.


Sign in / Sign up

Export Citation Format

Share Document