A Numerical Study of Unsteady Laminar Forced Convection From a Circular Cylinder

1976 ◽  
Vol 98 (2) ◽  
pp. 303-307 ◽  
Author(s):  
P. C. Jain ◽  
B. S. Goel

A numerical investigation of an unsteady laminar forced convection from a circular cylinder is presented. The Navier-Stokes equations and the energy equation for an unsteady incompressible fluid flow are solved by the finite difference method. The results are obtained at Reynolds numbers 100 and 200. The temperature field around the cylinder is obtained throughout the region of computation and is shown by isotherms at different times. The variations of the local Nusselt number around the cylinder at different times are computed and shown by graphs. The mean Nusselt number and the Strouhal number are also calculated. The computed results are compared with the other available experimental and theoretical results and are found to be in good agreement with them.

1991 ◽  
Vol 225 ◽  
pp. 557-574 ◽  
Author(s):  
Saul S. Abarbanel ◽  
Wai Sun Don ◽  
David Gottlieb ◽  
David H. Rudy ◽  
James C. Townsend

A detailed numerical study of two-dimensional flow past a circular cylinder at moderately low Reynolds numbers has been conducted using three different numerical algorithms for solving the time-dependent compressible Navier–Stokes equations. It was found that if the algorithm and associated boundary conditions were consistent and stable, then the major features of the unsteady wake were well predicted. However, it was also found that even stable and consistent boundary conditions could introduce additional periodic phenomena reminiscent of the type seen in previous wind-tunnel experiments. However, these additional frequencies were eliminated by formulating the boundary conditions in terms of the characteristic variables. An analysis based on a simplified model provides an explanation for this behaviour.


1974 ◽  
Vol 65 (3) ◽  
pp. 461-480 ◽  
Author(s):  
W. M. Collins ◽  
S. C. R. Dennis

The flow normal to an infinite circular cylinder which is uniformly accelerated from rest in a viscous fluid is considered. The flow is assumed to remain symmetrical about the direction of motion of the cylinder. Two types of solution are presented. In the first an expansion in powers of the time from the start of the motion is given which extends the results of boundary-layer theory by taking into account corrections for finite Reynolds numbers. Physical properties of the flow for small times and finite but large Reynolds numbers are calculated from this expansion. In the second method of solution the Navier-Stokes equations are integrated by an accurate procedure which is a logical extension of the solution in powers of the time. Results are obtained forR2= 97·5, 5850, 122 × 103and ∞, whereRis the Reynolds number. This is defined asR= 2a(ab)½/v, whereais the radius of the cylinder,bthe uniform acceleration andvthe kinematic viscosity of the fluid. The methods are in good agreement for small times.The numerical method of integration has been carried to moderate times and various flow properties have been calculated. The growth of the length of the separated wake behind the cylinder forR2= 97·5, 5850 and 122 × 103is compared with the results of recent experimental measurements. The agreement is only moderate forR2= 97·5 but it improves greatly asRincreases. The numerical integrations were continued in each case until the implicit method of integration failed to converge, which terminated the procedure. A secondary vortex appeared on the surface of the cylinder for the caseR2= 122 × 103.


2015 ◽  
Vol 3 (2) ◽  
pp. 28-49
Author(s):  
Ridha Alwan Ahmed

       In this paper, the phenomena of vortex shedding from the circular cylinder surface has been studied at several Reynolds Numbers (40≤Re≤ 300).The 2D, unsteady, incompressible, Laminar flow, continuity and Navier Stokes equations have been solved numerically by using CFD Package FLUENT. In this package PISO algorithm is used in the pressure-velocity coupling.        The numerical grid is generated by using Gambit program. The velocity and pressure fields are obtained upstream and downstream of the cylinder at each time and it is also calculated the mean value of drag coefficient and value of lift coefficient .The results showed that the flow is strongly unsteady and unsymmetrical at Re>60. The results have been compared with the available experiments and a good agreement has been found between them


1998 ◽  
Vol 14 (1) ◽  
pp. 23-29
Author(s):  
Robert R. Hwang ◽  
Sheng-Yuh Jaw

ABSTRACTThis paper presents a numerical study on turbulent vortex shedding flows past a square cylinder. The 2D unsteady periodic shedding motion was resolved in the calculation and the superimposed turbulent fluctuations were simulated with a second-order Reynolds-stress closure model. The calculations were carried out by solving numerically the fully elliptic ensemble-averaged Navier-Stokes equations coupled with the turbulence model equations together with the two-layer approach in the treatment of the near-wall region. The performance of the computations was evaluated by comparing the numerical results with data from available experiments. Results indicate that the present study gives good agreement in the shedding frequency and mean drag as well as in some phase profiles of the mean velocity.


1990 ◽  
Vol 220 ◽  
pp. 459-484 ◽  
Author(s):  
H. M. Badr ◽  
M. Coutanceau ◽  
S. C. R. Dennis ◽  
C. Ménard

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.


Author(s):  
Wei Ning ◽  
Li He

A numerical study has been carried out to investigate modelling issues on trailing edge vortex shedding. The vortex shedding from a circular cylinder and a VKI turbine blade is calculated using a 2-D unsteady multi-block Navier-Stokes solver. The unsteady stresses are calculated from the unsteady solutions. The distributions of the unsteady stresses are analysed and compared for the cylinder case and the cascade case, respectively. The time-averaged equations are then solved and the effectiveness of the “unsteady stresses” in suppressing trailing edge vortex shedding is checked. Finally, the time-independent solution produced by solving the time-averaged equations is compared with the time-averaged solution obtained by integrating the unsteady solutions. The numerical results have demonstrated that a time-independent vortex shedding solution can be achieved by solving the Navier-Stokes equations with the unsteady stresses and the time-averaged effects of the vortex shedding can be included.


1997 ◽  
Vol 119 (3) ◽  
pp. 556-567 ◽  
Author(s):  
Mihai Arghir ◽  
Jean Freˆne

The paper presents a method to calculate the rotordynamic coefficients of circumferentially-grooved liquid seals operating in centered position and turbulent flow regimes. The method is based on the integration of the averaged Navier-Stokes equations and uses a coordinate transformation proposed by Dietzen and Nordmann (1987). The effect of the coordinate transformation on the components of the stress tensor is included in the first order transport equations. To ensure grid independent solutions, numerical boundary conditions for the first-order velocities were formulated using the logarithmic law. The perturbation of the turbulent viscosity was also considered. A pressure recovery effect at the exit section was included in the first order mathematical model. The method is validated by calculations for straight and circumferentially-grooved seals. Comparisons with experimental and theoretical results show a good agreement for straight seals and for seals with few grooves, and a reasonable agreement for severe industrial cases (high Reynolds numbers and large number of grooves).


1997 ◽  
Vol 57 (1) ◽  
pp. 195-201 ◽  
Author(s):  
D. O. MARTÍNEZ ◽  
S. CHEN ◽  
G. D. DOOLEN ◽  
R. H. KRAICHNAN ◽  
L.-P. WANG ◽  
...  

High-resolution, direct numerical simulations of three-dimensional incompressible Navier–Stokes equations are carried out to study the energy spectrum in the dissipation range. An energy spectrum of the form A(k/kd)α exp[−βk/kd] is confirmed. The possible values of the parameters α and β, as well as their dependence on Reynolds numbers and length scales, are investigated, showing good agreement with recent theoretical predictions. A ‘bottleneck’-type effect is reported at k/kd≈4, exhibiting a possible transition from near-dissipation to far-dissipation.


1990 ◽  
Vol 210 ◽  
pp. 281-302 ◽  
Author(s):  
D. B. Ingham ◽  
T. Tang ◽  
B. R. Morton

A numerical and experimental study is described for the two-dimensional steady flow through a uniform cascade of normal flat plates. The Navier–Stokes equations are written in terms of the stream function and vorticity and are solved using a second-order-accurate finite-difference scheme which is based on a modified procedure to preserve accuracy and iterative convergence at higher Reynolds numbers. The upstream and downstream boundary conditions are discussed and an asymptotic solution is employed both upstream and downstream. A frequently used method for dealing with corner singularities is shown to be inaccurate and a method for overcoming this problem is described. Numerical solutions have been obtained for blockage ratio of 50 % and Reynolds numbers in the range 0 [les ]R[les ] 500 and results for both the lengths of attached eddies and the drag coefficients are presented. The calculations indicate that the eddy length increases linearly withR, at least up toR= 500, and that the multiplicative constant is in very good agreement with the theoretical prediction of Smith (1985a), who considered a related problem. In the case ofR= 0 the Navier–Stokes equations are solved using the finite-difference scheme and a modification of the boundary-element method which treats the corner singularities. The solutions obtained by the two methods are compared and the results are shown to be in good agreement. An experimental investigation has been performed at small and moderate values of the Reynolds numbers and there is excellent agreement with the numerical results both for flow streamlines and eddy lengths.


Sign in / Sign up

Export Citation Format

Share Document