couette flow

16
Couette Flow Analytical & Numerical Solution RASIKH TARIQ [email protected]

Upload: rizwanherl

Post on 15-Dec-2015

36 views

Category:

Documents


1 download

DESCRIPTION

The presentation about Coutte Flow

TRANSCRIPT

Couette Flow Analytical & Numerical Solution

RASIKH [email protected]

Group Members

Number Name Registration Number

1 Khawar Shahzad ME 113 009

2 Muhammad Ali ME 113 115

3 Mohsin Raza ME 113 080

4 Zohaib Ahmad ME 113 022

5 Aoun Abbas ME 113 074

6 Omer Ayub ME 113 089

7 Muhammad Ramzan ME 113 100

8 Khuram Yousaf ME 113 107

9 Zain Talib ME 113 108

10 Hassan Aamir ME 113 053

11 Adeel Anwar ME 113 039

2

Introduction

http://upload.wikimedia.org/wikipedia/commons/thumb/9/93/Laminar_shear.svg/800px-Laminar_shear.svg.png

3

Introduction

Laminar

Viscous

Incompressible Fluid

4

Introduction

5

Analytical Solution

The governing y-direction momentum equation is:

๐›’ ๐›› ๐ฏ๐›› ๐ญ

+๐›’๐ฎ ๐››๐ฏ๐››๐ฑ

+๐›’๐ฏ ๐›› ๐ฏ๐›› ๐ฒ

+๐›’๐ฐ ๐›› ๐ฏ๐››๐ณ

=โˆ’ ๐››๐ฉ๐››๐ฒ

+๐››๐›•๐ฑ๐ฒ

๐››๐ฑ+๐››๐›•๐ฒ๐ฒ

๐›› ๐ฒ+๐››๐›•๐ณ๐ฒ

๐››๐ณ+๐›’๐Ÿ ๐ฒ

Steady State Assumption No Body Forces

๐ŸŽ=๐››๐ฉ๐››๐ฒ

6

Analytical Solution

The governing x-direction momentum equation is:

๐›’ ๐››๐ฎ๐››๐ญ

+๐›’๐ฎ ๐››๐ฎ๐››๐ฑ

+๐›’ ๐ฏ ๐››๐ฎ๐›› ๐ฒ

+๐›’๐ฐ ๐››๐ฎ๐››๐ณ

=โˆ’ ๐๐’‘๐ ๐’™

+๐๐‰๐’™๐’™

๐๐’™+๐๐‰ ๐’š๐’™

๐ ๐’š+๐๐‰๐’›๐’™

๐๐’›+๐† ๐’‡ ๐’™

Steady State Assumption Y-Direction

Velocity (v) is Zero

No beginning or end of this flow in x-direction

Final Equation

๐œ•๐œ• y ( ฮผ ๐œ•u๐œ• y )=0

๐››๐Ÿ๐’–๐››๐ฒ๐Ÿ=๐ŸŽ

7

Using boundary condition at we have, we get.

Using boundary condition at we have, we get.

๐’–๐’–๐’†

=๐’š๐‘ซ

Analytical Solution

8

Analytical Solution โ€“ Matlab Codeclcclear%Workspace and Command History Cleared

ue = 0.01 %Velocity of Upper Plate moving with 0.01m/s

D = 0.05 %Separation Distance between plate. 5cm%y is distance from lower plate to upper plate

y = linspace(0,D,500);for i=1:1:500 U(i)=(ue*y(i))/D;endplot (y,U)legend ('Couette Flow Velocity Profile','Fontsize',12')xlabel ('Width Between Plates (m)',' Fontsize ',12')ylabel ('Velocity Profile (m/s)',' Fontsize ',12')axis ([-0.01 0.06 -0.001 0.012])

-0.01 0 0.01 0.02 0.03 0.04 0.05 0.06

0

2

4

6

8

10

12x 10

-3

Width Between Plates (m)

Ve

locity P

rofile

(m

/s)

Couette Flow Velocity Profile

9

Numerical Implicit Method

Nature of Partial Differential Equation

Parabolic Nature Equation because

Conclusion: Time Marching Possible

10

The Numerical Formulation

11

The Numerical Formulation

12

The Numerical Formulation

๐ฎ ๐ฃ๐ง+๐Ÿ=๐ฎ ๐ฃ

๐ง+๐šซ๐ญ

๐Ÿ (๐šซ๐ฒ )๐Ÿ๐‘๐ž๐ƒ

๐ฎ ๐ฃ+๐Ÿ๐ง+๐Ÿ+

๐šซ๐ญ๐Ÿ (๐šซ ๐ฒ )๐Ÿ๐‘๐ž๐ƒ

๐ฎ ๐ฃ+๐Ÿ๐ง โˆ’

๐šซ๐ญ๐Ÿ (๐šซ๐ฒ )๐Ÿ๐‘๐ž๐ƒ

๐Ÿ๐ฎ ๐ฃ๐ง+๐Ÿโˆ’

๐šซ๐ญ๐Ÿ (๐šซ๐ฒ )๐Ÿ๐‘๐ž๐ƒ

๐Ÿ๐ฎ ๐ฃ๐ง+

๐šซ๐ญ๐Ÿ (๐šซ๐ฒ )๐Ÿ๐‘๐ž๐ƒ

๐ฎ ๐ฃโˆ’๐Ÿ๐ง+๐Ÿ+

๐šซ๐ญ๐Ÿ (๐šซ ๐ฒ )๐Ÿ๐‘๐ž๐ƒ

๐ฎ ๐ฃโˆ’๐Ÿ๐ง

13

The Numerical Formulation

A B ๐Š ๐ฃ

๐€๐ฎ ๐ฃโˆ’ 1๐ง+1+๐๐ฎ ๐ฃ

๐ง+1+๐€๐ฎ ๐ฃ+1๐ง+1=๐Š ๐ฃ

14

The Numerical Formulation

๐€๐ฎ ๐ฃโˆ’ 1๐ง+1+๐๐ฎ ๐ฃ

๐ง+ 1+๐€ ๐ฎ ๐ฃ+1๐ง+1=๐Š ๐ฃ

[ยฟ๐ยฟ๐€ยฟ๐ŸŽยฟ๐ŸŽยฟ๐ŸŽ

ยฟ๐€ยฟ๐ยฟ๐€ยฟ๐ŸŽยฟ๐ŸŽ

ยฟ๐ŸŽยฟ๐€ยฟ๐ยฟ๐€ยฟ๐ŸŽ

ยฟ๐ŸŽยฟ๐ŸŽยฟ๐€ยฟ๐ยฟ๐€

ยฟ๐ŸŽยฟ๐ŸŽยฟ๐ŸŽยฟ๐€ยฟ๐

] [ยฟ๐ฎ๐Ÿ

๐ง+๐Ÿ

ยฟ๐ฎ๐Ÿ‘๐ง+๐Ÿ

ยฟ๐ฎ๐Ÿ’๐ง+๐Ÿ

ยฟ๐ฎ๐Ÿ“๐ง+๐Ÿ

ยฟ๐ฎ๐Ÿ”๐ง+๐Ÿ

]=[ยฟ๐Š ๐Ÿโˆ’๐€ ๐ฎ๐Ÿ

๐ง+๐Ÿ

ยฟ๐Š๐Ÿ‘

ยฟ๐Š๐Ÿ’

ยฟ๐Š๐Ÿ“

ยฟ๐Š๐Ÿ”

ยฟ๐Š ๐Ÿ”โˆ’๐€ ๐ฎ๐Ÿ•๐ง+๐Ÿ

]15

Thank you!

16