Scribd
Upload
93 views10 pages

ME608 - 2D SIMPLE Solver For Laminar Flow Over A Square Cylinder - Spring 2010

This document summarizes a numerical simulation of laminar 2D flow over a square cylinder using a SIMPLE solver. It details the geometry, mesh, boundary conditions, and discrete equations solved. Results are presented for velocity, pressure, vorticity, and streamlines at various Reynolds numbers. Validation is shown through good agreement between computed drag and recirculation length compared to published results. The solver uses a staggered mesh with Gauss-Seidel and LBL-TDMA methods to solve the discrete momentum and pressure correction equations respectively.

Uploaded by

saravana pandian
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
93 views10 pages

ME608 - 2D SIMPLE Solver For Laminar Flow Over A Square Cylinder - Spring 2010

This document summarizes a numerical simulation of laminar 2D flow over a square cylinder using a SIMPLE solver. It details the geometry, mesh, boundary conditions, and discrete equations solved. Results are presented for velocity, pressure, vorticity, and streamlines at various Reynolds numbers. Validation is shown through good agreement between computed drag and recirculation length compared to published results. The solver uses a staggered mesh with Gauss-Seidel and LBL-TDMA methods to solve the discrete momentum and pressure correction equations respectively.

Uploaded by

saravana pandian
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 10

ME608 - 2D SIMPLE Solver for Laminar Flow Over a

Square Cylinder - Spring 2010

Kameswararao Anupindi and Jeffrey R. Kennington

School of Mechanical Engineering


Purdue University
West Lafayette, IN 47907

April 22, 2010


Outline

Geometry, Mesh, and BCs


Solver details
Results
Validation
Summary
Details of Flow Domain

Figure: Geometry of the flow over 2D square cylinder

Side of the square D = 1 is the length scale.


H = 8D; L = 50D; Lin = 12.5D;
Channel is meshed with 500 x 80 control volumes with 10 x 10
across the square cylinder.
Boundary Conditions:
1 Inlet: Parabolic velocity inlet
2 Outlet: Upwinded mass flow outlet
3 Walls: No slip wall
Details of Solver

SIMPLE solver on staggered mesh


Gauss-Seidel for momentum equations
LBL-TDMA for pressure correction equation
Discrete u and v momentum equations
X
aeu ue = u
anb unb + beu + ∆y (PP − PE ) (1)
nb

X
anv vn = v
anb vnb + bnv + ∆x(PP − PN ) (2)
nb

Pressure correction equation


0 0
X
aP pP = anb pnb + (Fw∗ − Fe∗ + Fs∗ − Fn∗ ) (3)
nb
Results - Velocity

(a) Re 1 (b) Re 10

(c) Re 20 (d) Re 30

(e) Re 40 (f) Re 50

Figure: Velocity Profiles


Results - Pressure

(a) Re 1 (b) Re 10

(c) Re 20 (d) Re 30

(e) Re 40 (f) Re 50

Figure: Pressure Profiles


Results - Vorticity

Figure: Vorticity Profiles


Results - Streamlines

(a) Re 1 (b) Re 10 (c) Re 20

(d) Re 30 (e) Re 40 (f) Re 50

Figure: Streamlines
Validation

(a) Coefficient of Drag (b) Recirculation Length

Figure: Reynolds Number vs Coefficient of Drag and Recirculation Length


Summary

Closure
Laminar 2D flow over a square cylinder was simulated using SIMPLE
solver
Validated against published results with good agreement
References
ME 608 Class Notes
Accurate computations of the laminar flow past a square cylinder
based on two different methods: lattice-Boltzmann and finite-volume
by: M. Breuer, J. Bernsdorf, T. Zeiser, F. Durst
International Journal of Heat and Fluid Flow

You might also like

pFad - Phonifier reborn

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.


Alternative Proxies:

Alternative Proxy

pFad Proxy

pFad v3 Proxy

pFad v4 Proxy