scholarly journals Unstructured Grid Solutions for Incompressible Laminar Flow over a Circular Cylinder Using a Particular Finite Volume-Finite Element Method

2013 ◽  
Vol 2013 ◽  
pp. 1-9 ◽  
Author(s):  
Mahdi Yousefifard ◽  
Parviz Ghadimi ◽  
Rahim Zamanian

A numerical modeling of a 2D Navier-Stokes equation by a particular vertex centered control volume framework on an unstructured grid is presented in this paper. Triangular elements are applied with an effective high performance fully coupled algorithm, to simulate incompressible laminar flow over a circular cylinder. The cell face velocities in the discretization of the continuity and momentum equations are calculated by a combined linear and momentum interpolation scheme, respectively, and their performances are compared. Flow analyses have been conducted based on various Reynolds numbers up to 200 for the steady and unsteady flows using structured and unstructured grids. The robustness and accuracy of the scheme in the unstructured mesh are proved using the benchmark problems of incompressible laminar flow over a circular cylinder at low and medium Reynolds numbers. Results have been compared with the structured grid results, both cases with equal cell numbers and same strategy for the mesh refinement. Current results display good agreement with the experimental values. Overall, it is shown that, using the suggested method for the current problem, unstructured grids are highly competitive with the structured grids.

1981 ◽  
Vol 103 (4) ◽  
pp. 785-790 ◽  
Author(s):  
J. H. Masliyah ◽  
K. Nandakumar

The Navier-Stokes equation in a rotating frame of reference is solved numerically to obtain the flow field for a steady, fully developed laminar flow of a Newtonian fluid in a twisted tube having a square cross-section. The macroscopic force and energy balance equations and the viscous dissipation term are presented in terms of variables in a rotating reference frame. The computed values of friction factor are presented for dimensionless twist ratios, (i.e., length of tube over a rotation of π radians normalized with respect to half the width of tube) of 20, 10, 5 and 2.5 and for Reynolds numbers up to 2000. The qualitative nature of the axial velocity profile was observed to be unaffected by the swirling motion. The secondary motion was found to be most important near the wall.


1957 ◽  
Vol 2 (3) ◽  
pp. 237-262 ◽  
Author(s):  
Ian Proudman ◽  
J. R. A. Pearson

This paper is concerned with the problem of obtaining higher approximations to the flow past a sphere and a circular cylinder than those represented by the well-known solutions of Stokes and Oseen. Since the perturbation theory arising from the consideration of small non-zero Reynolds numbers is a singular one, the problem is largely that of devising suitable techniques for taking this singularity into account when expanding the solution for small Reynolds numbers.The technique adopted is as follows. Separate, locally valid (in general), expansions of the stream function are developed for the regions close to, and far from, the obstacle. Reasons are presented for believing that these ‘Stokes’ and ‘Oseen’ expansions are, respectively, of the forms $\Sigma \;f_n(R) \psi_n(r, \theta)$ and $\Sigma \; F_n(R) \Psi_n(R_r, \theta)$ where (r, θ) are spherical or cylindrical polar coordinates made dimensionless with the radius of the obstacle, R is the Reynolds number, and $f_{(n+1)}|f_n$ and $F_{n+1}|F_n$ vanish with R. Substitution of these expansions in the Navier-Stokes equation then yields a set of differential equations for the coefficients ψn and Ψn, but only one set of physical boundary conditions is applicable to each expansion (the no-slip conditions for the Stokes expansion, and the uniform-stream condition for the Oseen expansion) so that unique solutions cannot be derived immediately. However, the fact that the two expansions are (in principle) both derived from the same exact solution leads to a ‘matching’ procedure which yields further boundary conditions for each expansion. It is thus possible to determine alternately successive terms in each expansion.The leading terms of the expansions are shown to be closely related to the original solutions of Stokes and Oseen, and detailed results for some further terms are obtained.


1975 ◽  
Vol 68 (3) ◽  
pp. 591-608 ◽  
Author(s):  
L. E. Seeley ◽  
R. L. Hummel ◽  
J. W. Smith

Normal and tangential velocities in the boundary layer and out into the free stream have been obtained using a non-disturbing flow visualization technique for uniform laminar flow around a sphere. The non-similar data are available in tables at 2.5° intervals from 20° from the front to about 15° past the separation point a t Reynolds numbers of 290, 750, 1300 and 3000. Stream functions calculated by LeClair using a numerical solution of the Navier-Stokes equation at Re 21 300 are not in good agreement with measured values from 30° to 60°, but are in much better agreement around the separation point. Too few grid points near the sphere where the tangential velocities rise to a maximum above free-stream values may account for the difference.


Author(s):  
Lunji Song

To simulate incompressible Navier–Stokes equation, a temporal splitting scheme in time and high-order symmetric interior penalty Galerkin (SIPG) method in space discretization are employed, while the local Lax-Friedrichs flux is applied in the discretization of the nonlinear term. Under a constraint of the Courant–Friedrichs–Lewy (CFL) condition, two benchmark problems in 2D are simulated by the fully discrete SIPG method. One is a lid-driven cavity flow and the other is a circular cylinder flow. For the former, we compute velocity field, pressure contour and vorticity contour. In the latter, while the von Kármán vortex street appears with Reynolds number 50≤Re≤400, we simulate different dynamical behavior of circular cylinder flows, and numerically estimate the Strouhal numbers comparable to the existing experimental results. The calculations on vortex dominated flows are carried out to investigate the potential application of the SIPG method.


Aerospace ◽  
2021 ◽  
Vol 8 (8) ◽  
pp. 216
Author(s):  
Emanuel A. R. Camacho ◽  
Fernando M. S. P. Neves ◽  
André R. R. Silva ◽  
Jorge M. M. Barata

Natural flight has consistently been the wellspring of many creative minds, yet recreating the propulsive systems of natural flyers is quite hard and challenging. Regarding propulsive systems design, biomimetics offers a wide variety of solutions that can be applied at low Reynolds numbers, achieving high performance and maneuverability systems. The main goal of the current work is to computationally investigate the thrust-power intricacies while operating at different Reynolds numbers, reduced frequencies, nondimensional amplitudes, and mean angles of attack of the oscillatory motion of a NACA0012 airfoil. Simulations are performed utilizing a RANS (Reynolds Averaged Navier-Stokes) approach for a Reynolds number between 8.5×103 and 3.4×104, reduced frequencies within 1 and 5, and Strouhal numbers from 0.1 to 0.4. The influence of the mean angle-of-attack is also studied in the range of 0∘ to 10∘. The outcomes show ideal operational conditions for the diverse Reynolds numbers, and results regarding thrust-power correlations and the influence of the mean angle-of-attack on the aerodynamic coefficients and the propulsive efficiency are widely explored.


2007 ◽  
Vol 573 ◽  
pp. 171-190 ◽  
Author(s):  
A. DIPANKAR ◽  
T. K. SENGUPTA ◽  
S. B. TALLA

Vortex shedding behind a cylinder can be controlled by placing another small cylinder behind it, at low Reynolds numbers. This has been demonstrated experimentally by Strykowski & Sreenivasan (J. Fluid Mech. vol. 218, 1990, p. 74). These authors also provided preliminary numerical results, modelling the control cylinder by the innovative application of boundary conditions on some selective nodes. There are no other computational and theoretical studies that have explored the physical mechanism. In the present work, using an over-set grid method, we report and verify numerically the experimental results for flow past a pair of cylinders. Apart from providing an accurate solution of the Navier–Stokes equation, we also employ an energy-based receptivity analysis method to discuss some aspects of the physical mechanism behind vortex shedding and its control. These results are compared with the flow picture developed using a dynamical system approach based on the proper orthogonal decomposition (POD) technique.


1992 ◽  
Vol 114 (4) ◽  
pp. 936-943 ◽  
Author(s):  
Z. F. Dong ◽  
M. A. Ebadian

This paper numerically investigates the effects of buoyancy on fully developed laminar flow in a curved duct with an elliptic cross section. The flow of Newtonian fluids is assumed steady in terms of Boussinesq approximation. The curved elliptic duct is subjected to thermal boundary conditions of axially uniform heat flux and peripherally uniform wall temperature. The numerically generated boundary-fitted coordinate system is applied to discretize the solution domain of the elliptic duct, and the Navier-Stokes equations and the energy equation, including the curvature ratio, are solved by use of the control volume-based finite difference method. The solution covers a wide range of curvature ratios, and Dean and Grashof numbers. The results presented are displayed graphically and in tabular form to illustrate the buoyancy effect. It is further shown that buoyancy acts to increase both the Nusselt number and the friction factor and changes the distribution of the velocity and the temperature. The results for the curved circular duct with and without buoyancy are compared with the data available in the open literature for all cases. Also compared with the published data are the results of laminar flow in a curved elliptic duct, and very good agreement is obtained.


2013 ◽  
Vol 730 ◽  
pp. 145-161 ◽  
Author(s):  
Qiqi Wang ◽  
Jun-Hui Gao

AbstractThis paper analyses the adjoint solution of the Navier–Stokes equation. We focus on flow across a circular cylinder at three Reynolds numbers, ${\mathit{Re}}_{D} = 20, 100$ and $500$. The quantity of interest in the adjoint formulation is the drag on the cylinder. We use classical fluid mechanics approaches to analyse the adjoint solution, which is a vector field similar to a flow field. Production and dissipation of kinetic energy of the adjoint field is discussed. We also derive the evolution of circulation of the adjoint field along a closed material contour. These analytical results are used to explain three numerical solutions of the adjoint equations presented in this paper. The adjoint solution at ${\mathit{Re}}_{D} = 20$, a viscous steady state flow, exhibits a downstream suction and an upstream jet, the opposite of the expected behaviour of a flow field. The adjoint solution at ${\mathit{Re}}_{D} = 100$, a periodic two-dimensional unsteady flow, exhibits periodic, bean-shaped circulation in the near-wake region. The adjoint solution at ${\mathit{Re}}_{D} = 500$, a turbulent three-dimensional unsteady flow, has complex dynamics created by the shear layer in the near wake. The magnitude of the adjoint solution increases exponentially at the rate of the first Lyapunov exponent. These numerical results correlate well with the theoretical analysis presented in this paper.


2012 ◽  
Vol 2012 ◽  
pp. 1-10
Author(s):  
Hyun Jun Jeong ◽  
Wook Ryol Hwang ◽  
Chongyoup Kim

We present two-dimensional numerical simulations of the impact and spreading of a droplet containing a number of small particles on a flat solid surface, just after hitting the solid surface, to understand particle effects on spreading dynamics of a particle-laden droplet for the application to the industrial inkjet printing process. The Navier-Stokes equation is solved by a finite-element-based computational scheme that employs the level-set method for the accurate interface description between the drop fluid and air and a fictitious domain method for suspended particles to account for full hydrodynamic interaction. Focusing on the particle effect on droplet spreading and recoil behaviors, we report that suspended particles suppress the droplet oscillation and deformation, by investigating the drop deformations for various Reynolds numbers. This suppressed oscillatory behavior of the particulate droplet has been interpreted with the enhanced energy dissipation due to the presence of particles.


Author(s):  
Carl E. Rathmann

For well over 150 years now, theoreticians and practitioners have been developing and teaching students easily visualized models of fluid behavior that distinguish between the laminar and turbulent fluid regimes. Because of an emphasis on applications, perhaps insufficient attention has been paid to actually understanding the mechanisms by which fluids transition between these regimes. Summarized in this paper is the product of four decades of research into the sources of these mechanisms, at least one of which is a direct consequence of the non-linear terms of the Navier-Stokes equation. A scheme utilizing chaotic dynamic effects that become dominant only for sufficiently high Reynolds numbers is explored. This paper is designed to be of interest to faculty in the engineering, chemistry, physics, biology and mathematics disciplines as well as to practitioners in these and related applications.


Sign in / Sign up

Export Citation Format

Share Document