i need some help in fluid mechanics lab report with CFD

profilemoadjalsaljhc8
Untitled2.pdf

K B

50 05

: F lu

id S

ys te

m s

C om

pu ta

tio na

l F lu

id D

yn am

ic s

M oh

am m

ad .R

ah m

at i@

no rt

hu m

br ia

.a c.

uk R

oo m

:W J3

02 1

W ha

t y ou

le ar

n in

th is

le ct

ur e

• T

o gi

ve y

ou a

n in

tr od

uc tio

n on

C om

pu ta

tio na

l F lu

id

D yn

am ic

s (C

F D

)

• T

o pr

ov id

e an

a pp

re ci

at io

n fo

r th

e co

m pl

ex ity

o f C

F D

m

et ho

ds

• T

o gi

ve a

n id

ea a

bo ut

th e

po te

nt ia

l o f C

F D

in fl

ui d

m ec

ha ni

cs a

nd a

pp lic

at io

ns

2 M

T R

ah m

at i,

Fl ui

d S

ys te

m 2

01 8

M od

el s

in e

ng in

ee rin

g

• E

ng in

ee rs

cr ea

te m

od el

s of

sy st

em s

be ca

us e

th is

pr

oc es

s sa

ve s

m on

ey a

nd r

es ul

ts in

b et

te r

de si

gn s.

• M

od el

s ar

e to

ol s

w hi

ch r

ep re

se nt

s im

pl ifi

ed v

er si

on o

f re

al ity

• M

od el

s ar

e va

lu ab

le w

he n

th ey

h el

p en

gi ne

er s

an d

ot he

r pr

of es

si on

al s

re ac

h go

al s

in an

e co

no m

ic al

w ay

.

3 M

T R

ah m

at i,

Fl ui

d S

ys te

m 2

01 8

A na

ly tic

al F

lu id

D yn

am ic

s (A

F D

)

E xp

er im

en ta

l F

lu id

D yn

am ic

s (E

F D

)

C om

pu ta

tio na

l F

lu id

D yn

am ic

s (C

F D

)

W ha

t i s

C om

pu ta

tio na

l F lu

id D

yn am

ic s

(C F

D )

?

A na

ly tic

al F

lu id

D yn

am ic

s (A

F D

)

W ha

t i s

C om

pu ta

tio na

l F lu

id D

yn am

ic s

(C F

D )

?

Li m

ite d

kn ow

le dg

e

O nl

y so

lv ab

le fo

r si

m pl

e ca

se s

E xa

m pl

es th

at y

ou d

id : H

el e-

S ha

w fl

ow , f

lo w

in p

ip es

, C

ou et

te fl

ow e

tc .

y

x

) .

( 21

2 y

d y

xp u

GG P

H el

e- S

ha w

fl ow

E xp

er im

en ta

l F

lu id

D yn

am ic

s (E

F D

)

W ha

t i s

C om

pu ta

tio na

l F lu

id D

yn am

ic s

(C F

D )

?

E xp

en si

ve Ti

m e

co ns

um in

g D

iff ic

ul t t

o is

ol at

e su

b pr

ob le

m s

P ro

hi bi

tiv e

to p

er fo

rm

fu ll-

sc al

e te

st s

Im po

ss ib

le to

m ea

su re

in

ag gr

es si

ve e

nv ir

on m

en ts

E xa

m pl

es th

at y

ou h

av e

do ne

o r

yo u

ar e

do in

g th

e la

b: fl

ow in

a p

ip e

an d

flo w

o ve

r an

a irf

oi l

F lo

w o

ve r

an a

irf oi

l

C om

pu ta

tio na

l F

lu id

D yn

am ic

s (C

F D

)

W ha

t i s

C om

pu ta

tio na

l F lu

id D

yn am

ic s

(C F

D )

?

R ed

uc ed

ti m

e- to

-m ar

ke t

Fu ll-

sc al

e te

st in

g is

p os

si bl

e A

gr es

si ve

e nv

ir on

m en

ts C

om pl

ex c

om pl

ex p

ro bl

em s

th at

c an

no t e

ff ec

tiv el

y be

m od

el ed

w ith

an al

yt ic

al or

ex pe

ri m

en ta

l f lu

id m

ec ha

ni cs

. C

FD is

w id

el y

us ed

in in

du st

ry .

E xa

m pl

es th

at h

op ef

ul ly

y ou

w ill

d o:

F lo

w o

ve r

an a

irf oi

l a nd

fl ow

in a

p ip

e

R ol

e of

C F

D in

d es

ig n

of

B oe

in g

77 7

C F

D E

vo lu

tio n

at A

irb us

o ve

r 40

y ea

rs

Ir ro

ta tio

na l

In vi

sc id

R ot

at io

na l

In vi

sc id

R ot

at io

na l

V is

co us

C F

D P

ro ce

ss

1 2

3 4

5

Th e

S tr

ea m

F un

ct io

n an

d ve

lo ci

ty

dy dx

v

u

V )

, ,

( z

y x

\

B y

de fin

iti on

th e

ve lo

ci ty

o f a

fl ui

d pa

rt ic

le is

a lw

ay s

ta ng

en t t

o th

e st

re am

lin e

al on

g w

hi ch

it tr

av el

s. In

o th

er w

or ds

, a s

tr ea

m li

ne is

a

cu rv

e th

at is

e ve

ry w

he re

ta ng

en t t

o th

e lo

ca l v

el oc

ity v

ec to

r.

y

x st

re am

lin e

w k

vj ui

V �

U

M T

R ah

m at

i, Fl

ui d

S ys

te m

2 01

8

3D N

av ie

r- st

ok es

e qu

at io

ns fo

r In

co m

pr es

si bl

e flo

w

12 M

T R

ah m

at i,

Fl ui

d S

ys te

m 2

01 8

¸̧ ¹· ¨̈ ©§

ww �

ww �

ww �

ww �

ww

� ww

� ww

� ww

¸̧ ¹· ¨̈ ©§

ww �

ww �

ww �

ww �

ww

� ww

� ww

� ww

¸̧ ¹· ¨̈ ©§

ww �

ww �

ww �

ww �

ww

� ww

� ww

� ww

ww

� ww

� ww

2

2

2

2

2

2

2

2

2

2

2

2

2

2

2

2

2

2

111

0

zw yw

xw zp

zw w

yw v

xw u

tw

zv yv

xv yp

zv w

yv v

xv u

tv

zu yu

xu xp

zu w

yu v

xu u

tu

zw yv

xu

UP U

UP U

U U

UP U

M om

en tu

m e

qu at

io n

al on

g z

M om

en tu

m e

qu at

io n

al on

g y

M om

en tu

m e

qu at

io n

al on

g x

C on

tin ui

ty e

qu at

io n

A cc

el er

at io

n; a

Fo rc

e du

e to

p re

ss ur

e Fo

rc e

du e

to v

is co

si ty

3D N

av ie

r- st

ok es

e qu

at io

ns fo

r In

co m

pr es

si bl

e flo

w

13 M

T R

ah m

at i,

Fl ui

d S

ys te

m 2

01 8

T he

se P

D E

s re

pr es

en t t

he sc

ie nt

ifi c

la w

s th

at go

ve rn

flo w

in g

flu id

s. S

ol vi

ng th

es e

eq ua

tio ns

gi ve

s nu

m er

ic al

va lu

es fo

r th

e pr

es su

re fi

el d

an d

th e

ve lo

ci ty

T hu

s, so

lv in

g N

av ie

r- st

ok es

e qu

at io

n is

th e

ul tim

at e

so lu

tio n

te ch

ni qu

e.

A n

ex ac

t s ol

ut io

n in

vo lv

es a

p hy

si ca

ls itu

at io

n in

w hi

ch th

e eq

ua tio

ns of

m ot

io n

re du

ce

to e

qu at

io ns

th at

c an

b e

so lv

ed . T

he re

a re

a bo

ut fe

w s

uc h

so lu

tio ns

in

ex is

te nc

e. E

xa m

pl es

in cl

ud e

C ou

et te

flo w

a nd

H el

e- S

ha w

fl ow

:

La m

in ar

F lo

w B

et w

ee n

S ta

tio na

ry P

ar al

le l P

la te

s (H

el e-

S ha

w F

lo w

)

14 M

T R

ah m

at i,

Fl ui

d S

ys te

m 2

01 8

y

x

z

yG xGzG

3D v

ie w

yGyG

3D v

ie w ¸̧ ¹·

¨̈ ©§ ww

� ww

� ww

� ww

ww �

ww �

ww �

ww

¸̧ ¹· ¨̈ ©§

ww �

ww �

ww �

ww �

ww

� ww

� ww

� ww

¸̧ ¹· ¨̈ ©§

ww �

ww �

ww �

ww �

ww

� ww

� ww

� ww

2

2

2

2

2

2

2

2

2

2

2

2

2

2

2

2

2

2

111

zw yw

xw zp

zw w

yw v

xw u

tw

zv yv

xv yp

zv w

yv v

xv u

tv

zu yu

xu xp

zu w

yu v

xu u

tu

UP U

UP U

U U

UP U

B ec

au se

w =0

B ec

au se

v =0

B ec

au se

v =0

B ec

au se

w =0

B ec

au se

fl ow

is s

te ad

y B

ec au

se V

el oc

ity c

ha ng

es o

nl y

in o

ne d

ire ct

io n;

y

B ec

au se

V el

oc ity

c ha

ng es

o nl

y in

o ne

d ire

ct io

n; y

Fl ow

is :

-S te

ad y

- O

ne d

im en

si on

al

u

¸̧ ¹· ¨̈ ©§ ww

� ww

2

2

0 yu

xp P

=

0

0

0 0

0 0

0

0

3D N

av ie

r- st

ok es

e qu

at io

n fo

r In

co m

pr es

si bl

e flo

w

15 M

T R

ah m

at i,

Fl ui

d S

ys te

m 2

01 8

¸̧ ¹· ¨̈ ©§

ww �

ww �

ww �

ww �

ww

� ww

� ww

� ww

¸̧ ¹· ¨̈ ©§

ww �

ww �

ww �

ww �

ww

� ww

� ww

� ww

¸̧ ¹· ¨̈ ©§

ww �

ww �

ww �

ww �

ww

� ww

� ww

� ww

ww

� ww

� ww

2

2

2

2

2

2

2

2

2

2

2

2

2

2

2

2

2

2

111

0

zw yw

xw zp

zw w

yw v

xw u

tw

zv yv

xv yp

zv w

yv v

xv u

tv

zu yu

xu xp

zu w

yu v

xu u

tu

zw yv

xu

UP U

UP U

U U

UP U

M om

en tu

m e

qu at

io n

al on

g z

M om

en tu

m e

qu at

io n

al on

g y

M om

en tu

m e

qu at

io n

al on

g x

C on

tin ui

ty e

qu at

io n

A cc

el er

at io

n; a

Fo rc

e du

e to

p re

ss ur

e Fo

rc e

du e

to v

is co

si ty

N

ot e:

if y

ou n

eg le

ct th

e vi

sc ou

s te

rm s

th e

eq ua

tio ns

ar

e ca

lle d

E ul

er e

qu at

io ns

.

R A

N S

e qu

at io

ns

T ur

bu le

nt fl

ow

• flu

id m

ov es

in v

er y

irr eg

ul ar

p at

hs

• ve

lo ci

ty a

t a p

oi nt

fl uc

tu at

es

T he

m et

ho d

us ed

is th

e m

os ta

cc ur

at e

ap pr

oa ch

fo r

m od

el lin

g tu

rb ul

en t f

lo w

. I t i

s ca

lle d

ca lle

d di

re ct

n um

er ic

al s

im ul

at io

n (D

N S

}, in

vo lv

es se

tti ng

th e

gr id

a nd

tim

e st

ep s

fin e

en ou

gh to

re so

lv e

th e

fe at

ur es

o f t

he tu

rb ul

en tf

lo w

. T he

D

N S

so lu

tio n

is ve

ry a

cc ur

at e

bu tf

or f

or 9

9. 9%

of C

F D

pr ob

le m

s is

c om

pu ta

tio na

lly

un re

al is

tic .

T hi

s is

be ca

us e

th e

re qu

ire d

co m

pu te

r tim

e is

to o

la rg

e fo

r to

da y'

s co

m pu

te rs

.

Tu rb

ul en

ce M

od el

in g.

Tu rb

ul en

ce m

od el

in g

in vo

lv es

th e

pr ed

ic tio

n of

e ffe

ct s

of

tu rb

ul en

ce by

ap pl

yi ng

si m

pl ifi

ed eq

ua tio

ns . T

he se

eq ua

tio ns

ar e

si m

pl er

th at

th

e fu

ll tim

e- de

pe nd

en t N

av ie

r- S

to ke

s eq

ua tio

ns .O

ne o

f t he

m os

t w id

el y

us ed

tu

rb ul

en ce

m od

el s

is ca

lle d

th e

k- ep

si lo

n m

od el

o r

k- ep

sl io

n. m

od el

.T he

se eq

ua tio

ns ar

e us

ed to

ge th

er w

ith th

e R

ey no

ld s-

av er

ag ed

N

av ie

r- S

to ke

s (R

A N

S )

eq ua

tio ns

.

D is

cr et

iz at

io n:

To s

ol ve

N S

E , w

e m

us t c

on ve

rt g

ov er

ni ng

P D

E ’s

to a

lg eb

ra ic

e qu

at io

ns .

T he

p ro

ce ss

o fd

is cr

et iz

at io

n in

vo lv

es d

ev el

op in

g a

se t o

f a lg

eb ra

ic

eq ua

tio ns

w hi

ch is

u se

d in

p la

ce o

f t he

p ar

tia l d

iff er

en tia

l e qu

at io

ns . O

f th

e va

rio us

d is

cr et

iz at

io n

te ch

ni qu

es a

va ila

bl e

fo r

th e

nu m

er ic

al s

ol ut

io n

of th

e go

ve rn

in g

di ffe

re nt

ia l e

qu at

io ns

, t he

fo llo

w in

g th

re e

F in

ite d

iff er

en ce

m et

ho ds

( F

D M

) E

ac h

te rm

in N

S E

a pp

ro xi

m at

ed u

si ng

T ay

lo r

se rie

s, e

.g .,

F in

ite v

ol um

e m

et ho

ds (

F V

M )

U se

C on

tr ol

V ol

um e

fo rm

o f N

S E

e qu

at io

ns o

n ea

ch g

rid c

el l

M os

t p op

ul ar

a pp

ro ac

h, e

sp ec

ia lly

fo r

co m

m er

ci al

c od

es

F in

ite e

le m

en t m

et ho

ds (

F E

M )

S ol

ve P

D E

’s b

y re

pl ac

in g

co nt

in uo

us fu

nc tio

ns b

y pi

ec ew

is e

ap pr

ox im

at io

ns d

ef in

ed o

n po

ly go

ns , w

hi ch

a re

r ef

er re

d to

a s

el em

en ts

.

� �

� �

� �

1

2 2

1 1

2 2

2

i i

i i

i

U U

U O

x x

x U

U U

U O

x x

x

� � �

w �

' w

' w

� �

' w

'

0

� ¸ ¹·

¨ ©§ S

dxdT k

dxd

C on

si de

r 1 -d

im en

si on

al c

on du

ct io

n of

h ea

t:

w he

re :

k =

th er

m al

c on

du ct

iv ity

T =

te m

pe ra

tu re

S =

r at

e of

h ea

t g en

er at

io n

pe r

un it

vo lu

m e

D is

cr et

iz at

io n:

a n

ex am

pl e

C on

si de

r t he

1 -d

im en

si on

al g

rid s

ys te

m ( '

y an

d '

z =

1) :

W P

E

w e

w x)

(G e

x) (G

x '

W e

in te

gr at

e th

e he

at e

qu at

io n

ov er

th is

v ol

um e:

0

� ¸ ¹·

¨ ©§ �

¸ ¹· ¨ ©§

³ dx

S dxdT

k dxdT

k e w

w e

D is

cr et

iz at

io n:

0

� ¸ ¹·

¨ ©§ �

¸ ¹· ¨ ©§

³ dx

S dxdT

k dxdT

k e w

w e

W e

ne ed

to a

ss um

e a

pr of

ile to

e va

lu at

e th

e gr

ad ie

nt te

rm s:

W w

P e

E W

w P

e E

dT /d

x N

O T

de fin

ed fo

r th

is p

ro fil

e!

D is

cr et

iz at

io n:

If w

e us

e th

e lin

ea r p

ro fil

e as

su m

pt io

n:

W w

P e

E

W P

E

w e

w x)

(G e

x) (G

0 )

( )

( )

( )

(

' �

� �

� V

S x T

T k

x T

T k

w

W P

w

e

P E

e

G G

W e

ar riv

e at

th e

fo llo

w in

g di

sc re

tiz ed

e qu

at io

n (w

ith u

ni fo

rm c

el ls

):

0

� ¸ ¹·

¨ ©§ �

¸ ¹· ¨ ©§

³ dx

S dxdT

k dxdT

k e w

w e

D is

cr et

iz at

io n:

T he

d is

cr et

iz ed

e qu

at io

n is

w rit

te n

as :

¦ �

b T

a b

T a

T a

T a

nb nb

W W

E E

P P

w he

re :

( )

( )

e w

E W

P E

W e

w

k k

a a

a a

a b

S x

x x

d d

= =

= +

= D

In 2

d im

en si

on s,

th is

is o

fte n

ill us

tr at

ed a

s an

a m

oe ba

:

D is

cr et

iz at

io n:

S tr

uc tu

re d

G rid

s

xT he

gr id

is la

id ou

ti n

a re

gu la

r re

pe at

in g

pa tte

rn c

al le

d a

bl oc

k x

D et

ai ls

ar e

sp ec

ifi ed

by th

e us

er .

xT he

ad va

nt ag

e: th

e us

er ca

n se

tu p

th e

gr id

to m

ax im

iz e

ac cu

ra cy

w hi

le

ac hi

ev in

g ac

ce pt

ab le

r un

ti m

e.

xD ra

w ba

ck :i

tc an

ta ke

si gn

ifi ca

nt tim

e fo

r th

e us

er to

in pu

t t he

pa ra

m et

er s

ne ed

ed

to cr

ea te

th e

gr id

.

U ns

tr uc

tu re

d G

rid s

U ns

tr uc

tu re

d gr

id a

ro un

d an

a irf

oi l,

R ah

m at

i, 20

06

U ns

tr uc

tu re

d gr

id a

ro un

d a

bl ad

e,

R ah

m at

i 2 00

6

x It

is ba

se d

on a

c om

pu te

r al

go rit

hm th

at se

le ct

s an

a rb

itr ar

y co

lle ct

io n

of

el em

en ts

to fil

lt he

so lu

tio n

do m

ai n.

B e

ca us

e th

e el

em en

ts la

ck a

pa tte

rn ,

th e

gr id

is ca

lle d

un st

ru ct

ur ed

. x

It is

w el

ls ui

te d

fo r

no vi

ce s

be ca

us e

th e

gr id

ca n

be se

tu p

ea si

ly

an d

qu ic

kl y

an d

do es

n ot

r eq

ui re

m uc

h us

er ex

pe rt

is e.

x T

he d

ra w

ba ck

s ar

e th

at th

e gr

id m

ay

no t b

e go

od as

a st

ru ct

ur ed

g rid

in

te rm

s of

a cc

ur ac

y an

d so

lu tio

n tim

e.

H yb

rid G

rid s

A n

ax ia

l t ur

bo m

ac hi

ne ry

b la

de h

yb rid

m es

h, R

ah m

at i 2

00 8

x H

yb rid

gr id

m et

ho ds

ar e

de si

gn ed

to ta

ke ad

va nt

ag e

of th

e po

si tiv

e as

pe ct

s of

bo th

st ru

ct ur

ed a

nd un

st ru

ct ur

ed gr

id s.

H yb

rid gr

id s

us e

a st

ru ct

ur ed

g rid

in lo

ca l r

eg io

ns w

hi le

us in

g an

un st

ru ct

ur ed

gr id

in th

e bu

lk

of th

e do

m ai

n.

B ou

nd ar

y co

nd iti

on

P er

io di

c P er

io di

c

W al

l N

o- sl

ip

V el

oc ity

in le

t

P re

ss ur

e ou

tle t

Ty pi

ca l c

on di

tio ns

W al

l N

o- sl

ip (

u =

v =

w =

0 )

S lip

( ta

ng en

tia l s

tr es

s =

0, n

or m

al v

el oc

ity

= 0)

V el

oc ity

in le

t P

re ss

ur e

ou tle

t S

ym m

et ry

a nd

P er

io di

ci ty

R ah

m at

i 2 01

3

-A s

ol ve

r is

th e

co m

pu te

r al

go rit

hm th

at s

ol ve

s th

e al

ge br

ai c

eq ua

tio ns

u se

d by

th e

C F

D c

od e.

-T he

o ut

pu ts

o f t

he s

ol ve

r ar

e th

e va

lu es

o f t

he

ve lo

ci ty

, p re

ss ur

e, a

nd o

th er

r el

ev an

t f ie

ld s.

-A n

in iti

al v

el oc

ity a

nd p

re ss

ur e

fie ld

a re

as

su m

ed a

nd th

ey a

re c

or re

ct ed

u nt

il th

e go

ve rn

in g

eq ua

tio ns

a re

s at

is fie

d fo

r al

l c el

ls .

-T he

d et

ai ls

b ey

on d

th e

sc op

e of

th is

m od

ul e.

S ol

ve r

V er

ifi ca

tio n

ex am

in es

th e

de gr

ee to

w hi

ch th

e nu

m er

ic al

m et

ho ds

u se

d by

th e

co de

re

su lt

in ac

cu ra

te an

sw er

s. V

er ifi

ca tio

n ca

n in

vo lv

e va

ry in

g th

e sp

ac in

g in

th e

gr id

an d

en su

rin g

th at

th e

pr ed

ic te

d re

su lts

a re

no td

ep en

de nt

on th

e gr

id s

pa ci

ng .

V er

ifi ca

tio n:

S ol

vi ng

th e

eq ua

tio ns

r ig

ht

Va lid

at io

n ex

am in

es th

e de

gr ee

to w

hi ch

C F

D p

re di

ct io

ns a

gr ee

w ith

re al

-w or

ld o

bs er

va tio

ns .

A c

om m

on va

lid at

io n

st ra

te gy

is to

sy st

em at

ic al

ly co

m pa

re C

F D

pr ed

ic tio

ns to

e xp

er im

en ta

l da

ta or

to so

lu tio

ns to

w el

l-k no

w n

pr ob

le m

s, c

al le

d be

nc hm

ar k

so lu

tio ns

.

Va lid

at io

n: S

ol vi

ng th

e rig

ht e

qu at

io ns

C F

D P

ro ce

ss

1 2

3 4

5

P hy

si cs

: F lo

w o

ve r

a C

yl in

de r

S te

p 1:

S pe

ci fy

th e

G ov

er ni

ng e

qu at

io ns

- C

on se

rv at

io n

of m

as s,

m om

en tu

m a

nd e

ne rg

y

S te

p 2:

S pe

ci fy

th e

do m

ai n

- F

lo w

a ro

un d

th e

cy lin

de r

S te

p 3:

A ss

um pt

io n

an d

si m

pl ifi

ca tio

ns -

Ir ro

ta tio

na l

- In

vi sc

id

P ro

bl em

: M od

el th

e flo

w a

ro un

d a

cy lin

de r

w ith

d ia

m et

er o

f D .

(R e:

is k

no w

n)

0 2

2

2

2

ww

� ww

y x

M M

La pl

ac e

E qu

at io

n

co nt

in ui

ty P

ot en

tia l f

lo w

0 2

2

2

2

ww

� ww

y x

\ \

Ir ro

ta tio

na l

1

G rid

2 3

D is

cr et

is e

G rid

2 3

D is

cr et

is e

Ta yl

or s

er ie

s

... )

() (

!31 )

() (

!21 )

() (

!11 )

( )

( 3

0 3

3 2

0 2

2

0 0

� �

w w

� �

ww �

� w

w �

x

x x x

f x

x x x

f x

x xx f

x f

x f

) (

0x x

x �

'

... )

() (

!31 )

() (

!21 )

() (

!11 )

( )

( 3

3

3 2

2

2

0 �

' w

w �

' ww

� '

w w

x x

x x

f x

x x

f x

xx f

x f

x f

xx f

x x

f x

f w

w

'� )

( )

( )

( )

( 0

G rid

2 3

D is

cr et

is e

x x

j i

j i

'�

ww �

, ,1

\ \

\

y y

j i

j i

'�

ww �

, 1

, \

\ \

2 ,1

, ,1

2

2 2 x

x j

i j

i j

i

' �

ww �

� \

\ \

\

2 1

, ,

1 ,

2

2 2 y

y j

i j

i j

i

' �

ww �

� \

\ \

\

0 2

2

2

2

ww

� ww

y x

\ \

0 2

2 2

1 ,

, 1

, 2

,1 ,

,1

' �

� �

' �

� �

� �

y x

j i

j i

j i

j i

j i

j i

\ \

\ \

\ \

G rid

2 3

D is

cr et

is e

S ol

ve S

pe ci

fy th

e bo

un da

ry c

on di

tio ns

S

ol ve

it u

si ng

a c

om pu

te r

4 5 A

na ly

se :

V er

ify th

e co

de ,

Va lid

at e

th e

re su

lts u

si ng

th e

nu m

er ic

al o

r an

al yt

ic al

r es

ul ts

V is

ua liz

e th

e re

su lts

S ol

ve S

pe ci

fy th

e bo

un da

ry c

on di

tio ns

S

ol ve

it u

si ng

a c

om pu

te r

4 5 A

na ly

se :

V er

ify th

e co

de ,

Va lid

at e

th e

re su

lts u

si ng

th e

nu m

er ic

al o

r an

al yt

ic al

r es

ul ts

V is

ua liz

e th

e re

su lts

C F

D P

ro ce

ss

1 2

3 4

5

P hy

si cs

: F lo

w o

ve r

a C

yl in

de r

1

P ro

bl em

: M od

el th

e flo

w a

ro un

d a

cy lin

de r

w ith

d ia

m et

er o

f D .

(R e:

is k

no w

n)

¸̧ ¹· ¨̈ ©§

ww �

ww �

ww �

ww �

ww

� ww

� ww

� ww

¸̧ ¹· ¨̈ ©§

ww �

ww �

ww �

ww �

ww

� ww

� ww

� ww

¸̧ ¹· ¨̈ ©§

ww �

ww �

ww �

ww �

ww

� ww

� ww

� ww

ww

� ww

� ww

2

2

2

2

2

2

2

2

2

2

2

2

2

2

2

2

2

2

111

0

zw yw

xw zp

zw w

yw v

xw u

tw

zv yv

xv yp

zv w

yv v

xv u

tv

zu yu

xu xp

zu w

yu v

xu u

tu

zw yv

xu

UP U

UP U

U U

UP U

S te

p 1:

S pe

ci fy

th e

G ov

er ni

ng e

qu at

io ns

- C

on se

rv at

io n

of m

as s,

m om

en tu

m a

nd e

ne rg

y

3D N

av ie

r- st

ok es

eq

ua tio

ns fo

r In

co m

pr es

si bl

e flo

w

S te

p 2:

S pe

ci fy

th e

do m

ai n

- F

lo w

a ro

un d

th e

cy lin

de r

S te

p 3:

A ss

um pt

io n

an d

si m

pl ifi

ca tio

ns -

U si

ng a

tu rb

ul en

ce m

od el

to ta

ke in

to

ac co

un t t

he e

ffe ct

o f v

is co

si ty

A pp

lic at

io ns

o f C

F D

C om

pu ta

tio na

l d om

ai n

& B

ou nd

ar y

co nd

iti on

s

W al

l, N

o- sl

ip (

u =

v =

w =

0 )

S lip

( ta

ng en

tia l s

tr es

s =

0, n

or m

al v

el oc

ity =

0 )

V el

oc ity

in le

t P

re ss

ur e

ou tle

t S

ym m

et ry

a nd

P er

io di

ci ty

N on

-s lip

c on

di tio

n

V is

co si

ty (μ

) c an

b e

co ns

id er

ed a

s th

e st

ic ki

ne ss

of

fl ui

ds .

B ec

au se

o f v

is co

si ty

, a fl

ui d

in d

ir ec

t c on

ta ct

w ith

a s

ol id

“s

tic ks

” to

th e

su rf

ac e

an d

th er

e is

n o

sl ip

. T hi

s is

kn

ow n

as th

e no

-s lip

co

nd iti

on w

hi ch

m ea

ns th

e lo

ca l f

lo w

v el

oc ity

a t s

ol id

w

al ls

is z

er o.

In vi

sc id

fl ow

V is

co us

fl ow

37 M

T R

ah m

at i,

Fl ui

d S

ys te

m 2

01 8

P er

io di

c B

ou nd

ar ie

s

P er

io di

c

P er

io di

c

B y

us in

g pe

rio di

c bo

un da

ry c

on di

tio n

on ly

a s

in gl

e m

od el

is

m od

el le

d us

in g

C F

D

S ym

m et

ry B

ou nd

ar y

co nd

iti on

T he

s ym

m et

ry b

ou nd

ar y

co nd

iti on

fo rc

es fl

ow fi

el d

va ria

bl es

to b

e m

irr or

-im ag

ed

ac ro

ss a

s ym

m et

ry p

la ne

. T o

im po

se th

at b

ou nd

ar y

co nd

iti on

in a

c od

e, th

e gr

ad ie

nt s

of fl

ow fi

el d

va ria

bl es

in th

e di

re ct

io n

no rm

al to

th e

sy m

m et

ry p

la ne

a re

s et

to z

er o

ac ro

ss th

e pl

an e

of s

ym m

et ry

. B y

us in

g th

is b

ou nd

ar y

co nd

iti on

it is

p os

si bl

e to

m od

el a

po

rt io

n of

th e

ph ys

ic al

fl ow

d om

ai n

in o

rd er

to s

av e

co m

pu te

r re

so ur

ce s.

G rid

2 3

D is

cr et

is e

T he

N -S

e qu

at io

ns a

re a

re c

on ve

rt ed

to to

a lg

eb ra

ic e

qu at

io ns

. T

he p

ro ce

ss o

fd is

cr et

iz at

io n

in vo

lv es

d ev

el op

in g

a se

t o f a

lg eb

ra ic

e qu

at io

ns

w hi

ch is

u se

d in

p la

ce o

f t he

p ar

tia l d

iff er

en tia

l e qu

at io

ns .

S ol

ve S

pe ci

fy th

e bo

un da

ry

co nd

iti on

s S

ol ve

it u

si ng

a

co m

pu te

r

4

5 A

na ly

se :

V er

ify th

e co

de ,

Va lid

at e

th e

re su

lts u

si ng

th e

nu m

er ic

al o

r an

al yt

ic al

r es

ul ts

V is

ua liz

e th

e re

su lts

S ea

n D

al to

n, s

tu de

nt p

ro je

ct , 2

01 7

A pp

lic at

io ns

T

he se

w ak

es c

an le

d to

V or

te x

in du

ce d

vi br

at io

n (V

IV )

w hi

ch is

a v

er y

im po

rt an

t p he

no m

en on

to c

on si

de r

in m

an y

fie ld

s of

e ng

in ee

rin g

e. g.

c iv

il an

d m

ar in

e en

gi ne

er in

g. T

he y

ar e

kn ow

n to

h av

e ca

us ed

c om

po ne

nt

fa ilu

re in

a r

an ge

o f i

nd us

tr ia

l c om

po ne

nt s,

s uc

h as

b rid

ge s

up po

rt

st ru

ct ur

es , m

ar in

e ca

bl es

. R

is er

p ip

es u

se d

to e

xt ra

ct o

il fo

rm d

ee p

oc ea

ns a

nd p

ile s

ar e

an ot

he r

ex am

pl e

of s

tr uc

tu re

s ex

po se

d to

ti da

l w

av es

o r

riv er

fl ow

in w

hi ch

V IV

n ee

ds to

b e

co ns

id er

ed .

F ix

ed a

nd fl

ex ib

le r

is er

p ip

es u

se d

to e

xt ra

ct o

il fo

rm th

e de

ep s

ea s.

V or

te x

ca n

le d

to

vi br

at io

n of

p ip

es

an d

st ru

ct ur

al

fa ilu

re

La bs

e xp

er im

en t l

ea rn

in g

ou tc

om es

:

1. To

c om

m un

ic at

e te

ch ni

ca l i

nf or

m at

io n

us in

g w

rit te

n an

d gr

ap hi

ca l t

ec hn

iq ue

s. 2.

To a

na ly

ze th

e re

su lts

a nd

d is

cu ss

a pp

ro pr

ia te

ly

3. T

o be

a bl

e to

c rit

ic al

ly a

ss es

s th

e ex

pe rim

en t.

4. T

o co

m pa

re th

e ex

pe rim

en ta

l r es

ul ts

w ith

n um

er ic

al r

es ul

ts

T hi

s w

rit te

n re

po rt

is w

or th

3 0%

o f t

he m

od ul

e (o

r 10

0% o

f t he

la b

co m

po ne

nt ).

Yo u

w ill

p er

fo rm

b ot

h la

bs b

ut w

ill s

ub m

it a

w rit

te n

re po

rt o

n on

ly o

ne o

f t he

tw

o la

b ex

pe rim

en ts

id en

tif ie

d by

L ab

1- La

b2 :

La b1

: F lo

w in

si de

a p

ip e

La b2

: F lo

w o

ve r

an A

er of

oi l

Yo u

m ay

s ub

m it

a re

po rt

o n

La b1

o r

La b2

d ep

en di

ng o

n yo

ur p

re fe

re nc

es . Y

ou

sh ou

ld u

se y

ou r

ow n

da ta

g en

er at

ed d

ur in

g th

e la

b se

ss io

n.

In te

rm s

of m

ar ki

ng th

e la

b re

po rt

s, th

er e

ar e

no p

re fe

re nc

es b

et w

ee n

La b1

an

d La

b2 ; b

ot h

la bs

h av

e eq

ua l m

ar ks

a nd

im po

rt an

ce .

Ite m

M

ar k

1 S

ho rt

in tr

od uc

tio n

an d

ai m

s an

d ov

er al

l p re

se nt

at io

n

10

2 E

xp er

im en

ta l R

es ul

ts : T

ab le

s an

d G

ra ph

s 20

3 N

um er

ic al

m et

ho d

(C F

D )

an d

re su

lts

50

4 C

om pa

ris on

s an

d D

is cu

ss io

ns

20

To ta

l 1 00

( 3

0% o

f t he

to ta

l m od

ul e)

) 1(

4 4

22 2

2 2

Rr dxdp

R dxdp

R r

u �

P P

m ax

u av g

u

La b1

: F lo

w in

si de

a p

ip e

(N ot

th e

sa m

e co

nd iti

on a

s th

e la

b ex

pe rim

en t)

F lo

w in

si de

a p

ip e

La b2

: F lo

w o

ve r

an A

er of

oi l(

N ot

th e

sa m

e co

nd iti

on a

s th

e la

b ex

pe rim

en t)