An Immersed Boundary Direct Numerical Simulation Study of the Wall Effect on the Dynamics of a Rigid Sphere

Author(s):  
Jason Gatewood ◽  
Zhi-Gang Feng

The presence of a wall near a rigid sphere is known to disturb the particle fore and aft flow field and thereby affect particle drag and lift. This effect has wide ranging implications in particulate flows such as the dynamics of blood cells in microvessels or the transport of particulates in channel and pipe flows. In this study, an Immersed Boundary Direct Numerical Simulation (IB-DNS) is used to predict the dynamics of a rigid spherical body in the presence of a wall at laminar flows. The wall effect is shown to be significant when the dimensionless ratio (L/D) of the particle diameter (D) to the wall distance (L) is less than 3, and when particle Reynolds number is less than 10. Based on the IB-DNS results, a correlation for the wall effect on drag coefficient is derived that can be used to predict the actual drag coefficient for rigid spheres under the influence of a wall for L/D between 0.75 and 3 and Reynolds number between 0.18 and 10. The data underlying the correlation developed herein is validated by comparison to published experimental, numerical, and analytical correlations. The application of the IB-DNS method to study the wall effect is both novel and significant. It is novel in that such an application is not yet demonstrated. It is significant in that it; (1) utilizes a uniform Cartesian fluid mesh and (2) requires no sub domains of higher grid resolution in the wall gap.

2018 ◽  
Vol 857 ◽  
pp. 878-906 ◽  
Author(s):  
T. Nagata ◽  
T. Nonomura ◽  
S. Takahashi ◽  
Y. Mizuno ◽  
K. Fukuda

In this study, direct numerical simulation of the flow around a rotating sphere at high Mach and low Reynolds numbers is conducted to investigate the effects of rotation rate and Mach number upon aerodynamic force coefficients and wake structures. The simulation is carried out by solving the three-dimensional compressible Navier–Stokes equations. A free-stream Reynolds number (based on the free-stream velocity, density and viscosity coefficient and the diameter of the sphere) is set to be between 100 and 300, the free-stream Mach number is set to be between 0.2 and 2.0, and the dimensionless rotation rate defined by the ratio of the free-stream and surface velocities above the equator is set between 0.0 and 1.0. Thus, we have clarified the following points: (1) as free-stream Mach number increased, the increment of the lift coefficient due to rotation was reduced; (2) under subsonic conditions, the drag coefficient increased with increase of the rotation rate, whereas under supersonic conditions, the increment of the drag coefficient was reduced with increasing Mach number; and (3) the mode of the wake structure becomes low-Reynolds-number-like as the Mach number is increased.


2002 ◽  
Vol 473 ◽  
pp. 379-388 ◽  
Author(s):  
P. BAGCHI ◽  
S. BALACHANDAR

The lift forces on rigid spheres entrained in a vortex and a linear shear flow are computed using a direct numerical simulation. The sphere Reynolds number is in the range 10 to 100. The lift coefficient in a vortex is shown to be nearly two orders of magnitude higher than that in a shear flow. The inviscid mechanism is shown to be inadequate to account for the enhanced lift force. The effect of free rotation of the sphere is also shown to be too small to account for the enhanced lift force. Flow structure around the sphere is studied to explain the generation of the strong lift force in a vortex.


Author(s):  
Alessandro Chiarini ◽  
Maurizio Quadrio

AbstractA direct numerical simulation (DNS) of the incompressible flow around a rectangular cylinder with chord-to-thickness ratio 5:1 (also known as the BARC benchmark) is presented. The work replicates the first DNS of this kind recently presented by Cimarelli et al. (J Wind Eng Ind Aerodyn 174:39–495, 2018), and intends to contribute to a solid numerical benchmark, albeit at a relatively low value of the Reynolds number. The study differentiates from previous work by using an in-house finite-differences solver instead of the finite-volumes toolbox OpenFOAM, and by employing finer spatial discretization and longer temporal average. The main features of the flow are described, and quantitative differences with the existing results are highlighted. The complete set of terms appearing in the budget equation for the components of the Reynolds stress tensor is provided for the first time. The different regions of the flow where production, redistribution and dissipation of each component take place are identified, and the anisotropic and inhomogeneous nature of the flow is discussed. Such information is valuable for the verification and fine-tuning of turbulence models in this complex separating and reattaching flow.


2012 ◽  
Vol 2012 ◽  
pp. 1-13 ◽  
Author(s):  
Mohit Katragadda ◽  
Nilanjan Chakraborty ◽  
R. S. Cant

A direct numerical simulation (DNS) database of freely propagating statistically planar turbulent premixed flames with a range of different turbulent Reynolds numbers has been used to assess the performance of algebraic flame surface density (FSD) models based on a fractal representation of the flame wrinkling factor. The turbulent Reynolds number Rethas been varied by modifying the Karlovitz number Ka and the Damköhler number Da independently of each other in such a way that the flames remain within the thin reaction zones regime. It has been found that the turbulent Reynolds number and the Karlovitz number both have a significant influence on the fractal dimension, which is found to increase with increasing Retand Ka before reaching an asymptotic value for large values of Retand Ka. A parameterisation of the fractal dimension is presented in which the effects of the Reynolds and the Karlovitz numbers are explicitly taken into account. By contrast, the inner cut-off scale normalised by the Zel’dovich flame thicknessηi/δzdoes not exhibit any significant dependence on Retfor the cases considered here. The performance of several algebraic FSD models has been assessed based on various criteria. Most of the algebraic models show a deterioration in performance with increasing the LES filter width.


2016 ◽  
Vol 30 (8) ◽  
pp. 6727-6737 ◽  
Author(s):  
Cong Xu ◽  
Zhihua Wang ◽  
Wubin Weng ◽  
Kaidi Wan ◽  
Ronald Whiddon ◽  
...  

Author(s):  
Zhi-gang Feng ◽  
Jason Gatewood ◽  
E.E. Michaelides

Abstract The presence of a wall near a rigid sphere in motion is known to disturb the particle fore and aft flow field symmetry and to affect the hydrodynamic force. An Immersed Boundary Direct Numerical Simulation (IB-DNS) is used in this study to determine the wall effects on the dynamics of a free-falling sphere and the drag of a sphere moving at a constant velocity. The numerical results are validated by comparison to the published experimental, numerical, and analytical data. The pressure and velocity fields are numerically computed when the particle is in the vicinity of the wall; the transverse (lift) and longitudinal (drag) parts of the hydrodynamic force are calculated; its rotational velocity is also investigated in the case of a free-falling sphere. The flow asymmetry also causes the particle to rotate. The wall effect is shown to be significant when the dimensionless ratio of the wall distance to the particle diameter, L/D, is less than 3. The wall effects are more pronounced and when the particle Reynolds number, Re, is less than 10. Based on the computational results, a useful correlation for the wall effects on the drag coefficients spheres is derived in the range 0.75 < L/D < 3 and 0.18 < Re < 10.


2021 ◽  
Author(s):  
Silin Jing ◽  
Xianzhi Song ◽  
Zhaopeng Zhu ◽  
Buwen Yu ◽  
Shiming Duan

Abstract Accurate description of cuttings slippage in the gas-liquid phase is of great significance for wellbore cleaning and the control accuracy of bottom hole pressure during MPD. In this study, the wellbore bubble flow environment was simulated by a constant pressure air pump and the transparent wellbore, and the settling characteristics of spherical particles under different gas volume concentrations were recorded and analyzed by highspeed photography. A total of 225 tests were conducted to analyze the influence of particle diameter (1–12mm), particle density (2700–7860kg/m^3), liquid viscosity and bubble volume concentration on particle settling velocity. Gas drag force is defined to quantitatively evaluate the bubble’s resistance to particle slippage. The relationship between bubble drag coefficient and particle Reynolds number is obtained by fitting the experimental results. An explicit settling velocity equation is established by introducing Archimedes number. This explicit equation with an average relative error of only 8.09% can directly predict the terminal settling velocity of the sphere in bubble containing Newtonian fluids. The models for predicting bubble drag coefficient and the terminal settling velocity are valid with particle Reynolds number ranging from 0.05 to 167 and bubble volume concentration ranging from 3.0% to 20.0%. Besides, a trial-and-error procedure and an illustrative example are presented to show how to calculate bubble drag coefficient and settling velocity in bubble containing fluids. The results of this study will provide the theoretical basis for wellbore cleaning and accurate downhole pressure to further improve the performance of MPD in treating gas influx.


Sign in / Sign up

Export Citation Format

Share Document