Direct simulation of initial value problems for the motion of solid bodies in a Newtonian fluid Part 1. Sedimentation

1994 ◽  
Vol 261 ◽  
pp. 95-134 ◽  
Author(s):  
J. Feng ◽  
H. H. Hu ◽  
D. D. Joseph
Keyword(s):  
Reynolds Number ◽  
Stokes Equations ◽  
Steady Motion ◽  
Vertical Channel ◽  
Periodic Oscillation ◽  
Reynolds Numbers ◽  
Two Dimensions ◽  
Fluid Particle ◽  
Initial Value ◽  
Circular Particle

This paper reports the result of direct simulations of fluid–particle motions in two dimensions. We solve the initial value problem for the sedimentation of circular and elliptical particles in a vertical channel. The fluid motion is computed from the Navier–Stokes equations for moderate Reynolds numbers in the hundreds. The particles are moved according to the equations of motion of a rigid body under the action of gravity and hydrodynamic forces arising from the motion of the fluid. The solutions are as exact as our finite-element calculations will allow. As the Reynolds number is increased to 600, a circular particle can be said to experience five different regimes of motion: steady motion with and without overshoot and weak, strong and irregular oscillations. An elliptic particle always turn its long axis perpendicular to the fall, and drifts to the centreline of the channel during sedimentation. Steady drift, damped oscillation and periodic oscillation of the particle are observed for different ranges of the Reynolds number. For two particles which interact while settling, a steady staggered structure, a periodic wake-action regime and an active drafting–kissing–tumbling scenario are realized at increasing Reynolds numbers. The non-linear effects of particle–fluid, particle–wall and interparticle interactions are analysed, and the mechanisms controlling the simulated flows are shown to be lubrication, turning couples on long bodies, steady and unsteady wakes and wake interactions. The results are compared to experimental and theoretical results previously published.

Capsular flow in pipelines

1972 ◽  
Vol 56 (1) ◽  
pp. 49-59 ◽  
Author(s):  
A. E. Vardy ◽  
M. I. G. Bloor ◽  
J. A. Fox
Keyword(s):  
Reynolds Number ◽  
Pressure Gradient ◽  
Numerical Solutions ◽  
Stokes Equations ◽  
Steady Motion ◽  
Reynolds Numbers ◽  
Navier Stokes ◽  
Hydraulic Pressure ◽  
Central Difference ◽  
Navier Stokes Equations

The problem considered is that of the steady motion of a series of neutrally buoyant, flat-faced, rigid, cylindrical capsules along the axis of a pipeline under the influence of a hydraulic pressure gradient. The Navier-Stokes equations are non-dimensionalized and expressed in central-difference form. Numerical solutions are found by the method of relaxation for Reynolds numbers up to 20 000 and a close agreement is obtained with readings from a laboratory apparatus for Reynolds numbers up to 2200.The flow is examined in detail and the existence of toroidal vortices between successive capsules is demonstrated. Their shape is shown to be increasingly influenced by inertial forces as the Reynolds number increases, but the overall pressure gradient is not greatly dependent on the Reynolds number.

Entrainment by a Refrigerated Air Curtain Down a Wall

2004 ◽  
Vol 126 (5) ◽  
pp. 871-879 ◽  
Author(s):  
Pratik Bhattacharjee ◽  
Eric Loth
Keyword(s):  
Reynolds Number ◽  
Stokes Equations ◽  
Reynolds Numbers ◽  
Two Dimensions ◽  
Transitional Flow ◽  
Navier Stokes ◽  
Wall Jets ◽  
Adiabatic Wall ◽  
Air Curtain ◽  
Thermal Entrainment

Thermal entrainment is important as it adversely affects energy consumption and evaporator humidity levels of refrigerated air curtain display cases, often at transitional Reynolds numbers. In order to get a more fundamental understanding of the mean and unsteady thermal entrainment processes, the shelf structure of a display case has been idealized to that of a plane, adiabatic wall subjected to refrigerated wall jets at laminar and transitional flow conditions. The wall jets are studied at different inflow profiles, Reynolds numbers, and Richardson numbers to investigate the effect on thermal entrainment rates. The primary simulation technique was direct numerical simulation of the Navier–Stokes equations in two dimensions for the low and moderate Reynolds numbers (though three-dimensional simulations were also conducted). At higher Reynolds numbers, a conventional Reynolds averaged Navier–Stokes approach was employed, which was found to give reasonable agreement with the above approach at a wall jet (early-transitional) Reynolds number of 2000. In general, the results yielded a significant variation in entrainment as a function of Reynolds number, with a minimum occurring at flow speeds immediately prior to transition. The entrainment rates were also sensitive to the initial velocity distribution, whereby a constant gradient profile (where any local velocity-gradient peaks were minimized) provided the least entrainment. Entrainment was also found to decrease with increasing Richardson number.

scholarly journals Finite-amplitude steady waves in plane viscous shear flows

1985 ◽  
Vol 160 ◽  
pp. 281-295 ◽  
Author(s):  
F. A. Milinazzo ◽  
P. G. Saffman
Keyword(s):  
Boundary Layer ◽  
Reynolds Number ◽  
Stokes Equations ◽  
Reynolds Numbers ◽  
Finite Amplitude ◽  
Linear Instability ◽  
Navier Stokes ◽  
Unidirectional Flow ◽  
Viscous Shear ◽  
Finite Amplitude Waves

Computations of two-dimensional solutions of the Navier–Stokes equations are carried out for finite-amplitude waves on steady unidirectional flow. Several cases are considered. The numerical method employs pseudospectral techniques in the streamwise direction and finite differences on a stretched grid in the transverse direction, with matching to asymptotic solutions when unbounded. Earlier results for Poiseuille flow in a channel are re-obtained, except that attention is drawn to the dependence of the minimum Reynolds number on the physical constraint of constant flux or constant pressure gradient. Attempts to calculate waves in Couette flow by continuation in the velocity of a channel wall fail. The asymptotic suction boundary layer is shown to possess finite-amplitude waves at Reynolds numbers orders of magnitude less than the critical Reynolds number for linear instability. Waves in the Blasius boundary layer and unsteady Rayleigh profile are calculated by employing the artifice of adding a body force to cancel the spatial or temporal growth. The results are verified by comparison with perturbation analysis in the vicinity of the linear-instability critical Reynolds numbers.

Lid-Driven Square Cavity Flow: A Benchmark Solution with an 8192x8192 Grid

2021 ◽  
Author(s):  
Carlos Marchi ◽  
Cosmo D. Santiago ◽  
Carlos Alberto Rezende de Carvalho Junior
Keyword(s):  
Reynolds Number ◽  
Stokes Equations ◽  
Accurate Determination ◽  
Stable Solution ◽  
Reynolds Numbers ◽  
Accurate Solution ◽  
Discretization Error ◽  
Square Cavity ◽  
Solution Vector ◽  
Wide Range

Abstract The incompressible steady-state fluid flow inside a lid-driven square cavity was simulated using the mass conservation and Navier-Stokes equations. This system of equations is solved for Reynolds numbers of up to 10,000 to the accuracy of the computational machine round-off error. The computational model used was the second-order accurate finite volume method. A stable solution is obtained using the iterative multigrid methodology with 8192 × 8192 volumes, while degree-10 interpolation and Richardson extrapolation were used to reduce the discretization error. The solution vector comprised five entries of velocities, pressure, and location. For comparison purposes, 65 different variables of interest were chosen, such as velocity profile, its extremum values and location, extremum values and location of the stream function. The discretization error for each variable of interest was estimated using two types of estimators and their apparent order of accuracy. The variations of the 11 selected variables are shown across 38 Reynolds number values between 0.0001 and 10,000. In this study, we provide a more accurate determination of the Reynolds number value at which the upper secondary vortex appears. The results of this study were compared with those of several other studies in the literature. The current solution methodology was observed to produce the most accurate solution till date for a wide range of Reynolds numbers.

Effect of wall-mounted V-baffle position in a turbulent flow through a channel

2019 ◽  
Vol 29 (10) ◽  
pp. 3908-3937 ◽  
Author(s):  
Younes Menni ◽  
Ahmed Azzi ◽  
Ali J. Chamkha ◽  
Souad Harmand
Keyword(s):  
Heat Transfer ◽  
Reynolds Number ◽  
Stokes Equations ◽  
Numerical Study ◽  
Rectangular Channel ◽  
Simple Algorithm ◽  
Two Dimensions ◽  
Large Reynolds Number ◽  
Constant Property ◽  
Content Type

Purpose The purpose of this paper is to carry out a numerical study on the dynamic and thermal behavior of a fluid with a constant property and flowing turbulently through a two-dimensional horizontal rectangular channel. The upper surface was put in a constant temperature condition, while the lower one was thermally insulated. Two transverse, solid-type obstacles, having different shapes, i.e. flat rectangular and V-shaped, were inserted into the channel and fixed to the top and bottom walls of the channel, in a periodically staggered manner to force vortices to improve the mixing, and consequently the heat transfer. The flat rectangular obstacle was put in the first position and was placed on the hot top wall of the channel. However, the second V-shaped obstacle was placed on the insulated bottom wall, at an attack angle of 45°; its position was varied to find the optimum configuration for optimal heat transfer. Design/methodology/approach The fluid is considered Newtonian, incompressible with constant properties. The Reynolds averaged Navier–Stokes equations, along with the standard k-epsilon turbulence model and the energy equation, are used to control the channel flow model. The finite volume method is used to integrate all the equations in two-dimensions; the commercial CFD software FLUENT along with the SIMPLE-algorithm is used for pressure-velocity coupling. Various values of the Reynolds number and obstacle spacing were selected to perform the numerical runs, using air as the working medium. Findings The channel containing the flat fin and the 45° V-shaped baffle with a large Reynolds number gave higher heat transfer and friction loss than the one with a smaller Reynolds number. Also, short separation distances between obstacles provided higher values of the ratios Nu/Nu0 and f/f0 and a larger thermal enhancement factor (TEF) than do larger distances. Originality/value This is an original work, as it uses a novel method for the improvement of heat transfer in completely new flow geometry.

scholarly journals Nonlinear amplification in hydrodynamic turbulence

2021 ◽  
Vol 930 ◽  
Author(s):  
Kartik P. Iyer ◽  
Katepalli R. Sreenivasan ◽  
P.K. Yeung
Keyword(s):  
Reynolds Number ◽  
Isotropic Turbulence ◽  
Stokes Equations ◽  
Three Dimensional ◽  
Reynolds Numbers ◽  
Navier Stokes ◽  
Navier Stokes Equations ◽  
Developed Turbulence ◽  
Advection Term ◽  
Nonlinear Amplification

Using direct numerical simulations performed on periodic cubes of various sizes, the largest being $8192^3$ , we examine the nonlinear advection term in the Navier–Stokes equations generating fully developed turbulence. We find significant dissipation even in flow regions where nonlinearity is locally absent. With increasing Reynolds number, the Navier–Stokes dynamics amplifies the nonlinearity in a global sense. This nonlinear amplification with increasing Reynolds number renders the vortex stretching mechanism more intermittent, with the global suppression of nonlinearity, reported previously, restricted to low Reynolds numbers. In regions where vortex stretching is absent, the angle and the ratio between the convective vorticity and solenoidal advection in three-dimensional isotropic turbulence are statistically similar to those in the two-dimensional case, despite the fundamental differences between them.

Jet-to-Plate Distance Effect on Heat Transfer From a Flat Plate to an Impinging Jet at Various Reynolds Numbers

2012 ◽  
Author(s):  
Patricia Streufert ◽  
Terry X. Yan ◽  
Mahdi G. Baygloo
Keyword(s):  
Heat Transfer ◽  
Experimental Data ◽  
Reynolds Number ◽  
Flat Plate ◽  
Numerical Solutions ◽  
Stokes Equations ◽  
Local Heat Transfer ◽  
Reynolds Numbers ◽  
Air Jet ◽  
Plate Distance

Local turbulent convective heat transfer from a flat plate to a circular impinging air jet is numerically investigated. The jet-to-plate distance (L/D) effect on local heat transfer is the main focus of this study. The eddy viscosity V2F turbulence model is used with a nonuniform structured mesh. Reynolds-Averaged Navier-Stokes equations (RANS) and the energy equation are solved for axisymmetric, three-dimensional flow. The numerical solutions obtained are compared with published experimental data. Four jet-to-plate distances, (L/D = 2, 4, 6 and 10) and seven Reynolds numbers (Re = 7,000, 15,000, 23,000, 50,000, 70,000, 100,000 and 120,000) were parametrically studied. Local and average heat transfer results are analyzed and correlated with Reynolds number and the jet-to-plate distance. Results show that the numerical solutions matched experimental data best at low jet-to-plate distances and lower Reynolds numbers, decreasing in ability to accurately predict the heat transfer as jet-to-plate distance and Reynolds number was increased.

Prediction of Reynolds Number Effects on Low-Pressure Turbines Using a High-Order ILES Method

2019 ◽  
Author(s):  
Marc Bolinches-Gisbert ◽  
David Cadrecha Robles ◽  
Roque Corral ◽  
Fernando Gisbert
Keyword(s):  
Boundary Layer ◽  
Reynolds Number ◽  
Fourth Order ◽  
Stokes Equations ◽  
Reynolds Numbers ◽  
Low Pressure ◽  
High Order ◽  
Experimental Results ◽  
Numerical Code ◽  
Design System

Abstract A comprehensive comparison between Implicit Large Eddy Simulations (ILES) and experimental results of a modern highlift low-pressure turbine airfoil has been carried out for an array of Reynolds numbers (Re). Experimental data were obtained in a low-speed linear cascade at the Polithecnic University of Madrid using hot-wire anemometry and LDV. The numerical code is fourth order accurate, both in time and space. The spatial discretization of the compressible Navier-Stokes equations is based on a high-order Flux Reconstruction approach while a fourth order Runge-Kutta method is used to march in time the simulations. The losses, pressure coefficient distributions, and boundary layer and wake velocity profiles have been compared for an array of realistic Reynolds numbers. Moreover, boundary layer and wake velocity fluctuations are compared for the first time with experimental results. It is concluded that the accuracy of the numerical results is comparable to that of the experiments, especially for integral quantities such as the losses or exit angle. Turbulent fluctuations in the suction side boundary layer and the wakes are well predicted also. The elapsed time of the is about 140 hours on 40 Graphics Processor Units. The numerical tool is integrated within an industrial design system and reuses pre- and post-processing tools previously developed for another kind of applications. The trend of the losses with the Reynolds number has a sub-critical regime, where the losses scale with Re−1, and a supercrital regime, where the losses scale with Re−1/2. This trend can be seen both, in the simulations and the experiments.

Unsteady flow past a rotating circular cylinder at Reynolds numbers 103 and 104

1990 ◽  
Vol 220 ◽  
pp. 459-484 ◽  
Author(s):  
H. M. Badr ◽  
M. Coutanceau ◽  
S. C. R. Dennis ◽  
C. Ménard
Keyword(s):  
Reynolds Number ◽  
Circular Cylinder ◽  
Unsteady Flow ◽  
Rotation Rate ◽  
Stokes Equations ◽  
Reynolds Numbers ◽  
Navier Stokes ◽  
Navier Stokes Equations ◽  
Number Range ◽  
Flow Past

The unsteady flow past a circular cylinder which starts translating and rotating impulsively from rest in a viscous fluid is investigated both theoretically and experimentally in the Reynolds number range 103 [les ] R [les ] 104 and for rotational to translational surface speed ratios between 0.5 and 3. The theoretical study is based on numerical solutions of the two-dimensional unsteady Navier–Stokes equations while the experimental investigation is based on visualization of the flow using very fine suspended particles. The object of the study is to examine the effect of increase of rotation on the flow structure. There is excellent agreement between the numerical and experimental results for all speed ratios considered, except in the case of the highest rotation rate. Here three-dimensional effects become more pronounced in the experiments and the laminar flow breaks down, while the calculated flow starts to approach a steady state. For lower rotation rates a periodic structure of vortex evolution and shedding develops in the calculations which is repeated exactly as time advances. Another feature of the calculations is the discrepancy in the lift and drag forces at high Reynolds numbers resulting from solving the boundary-layer limit of the equations of motion rather than the full Navier–Stokes equations. Typical results are given for selected values of the Reynolds number and rotation rate.

Analytical and numerical studies of the structure of steady separated flows

1966 ◽  
Vol 24 (1) ◽  
pp. 113-151 ◽  
Author(s):  
Odus R. Burggraf
Keyword(s):  
Reynolds Number ◽  
Numerical Solutions ◽  
Stokes Equations ◽  
Low Reynolds Number ◽  
Reynolds Numbers ◽  
Total Power ◽  
Large Reynolds Number ◽  
Boundary Velocity ◽  
Displacement Thickness ◽  
Low Reynolds

The viscous structure of a separated eddy is investigated for two cases of simplified geometry. In § 1, an analytical solution, based on a linearized model, is obtained for an eddy bounded by a circular streamline. This solution reveals the flow development from a completely viscous eddy at low Reynolds number to an inviscid rotational core at high Reynolds number, in the manner envisaged by Batchelor. Quantitatively, the solution shows that a significant inviscid core exists for a Reynolds number greater than 100. At low Reynolds number the vortex centre shifts in the direction of the boundary velocity until the inviscid core develops; at large Reynolds number, the inviscid vortex core is symmetric about the centre of the circle, except for the effect of the boundary-layer displacement-thickness. Special results are obtained for velocity profiles, skin-friction distribution, and total power dissipation in the eddy. In addition, results of the method of inner and outer expansions are compared with the complete solution, indicating that expansions of this type give valid results for separated eddies at Reynolds numbers greater than about 25 to 50. The validity of the linear analysis as a description of separated eddies is confirmed to a surprising degree by numerical solutions of the full Navier–Stokes equations for an eddy in a square cavity driven by a moving boundary at the top. These solutions were carried out by a relaxation procedure on a high-speed digital computer, and are described in § 2. Results are presented for Reynolds numbers from 0 to 400 in the form of contour plots of stream function, vorticity, and total pressure. At the higher values of Reynolds number, an inviscid core develops, but secondary eddies are present in the bottom corners of the square at all Reynolds numbers. Solutions of the energy equation were obtained also, and isotherms and wall heat-flux distributions are presented graphically.

Sign in / Sign up

Export Citation Format

Share Document