Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Three Days Workshop On Computational Fluid Dynamics 7 To9 February 2014

Download as doc, pdf, or txt
Download as doc, pdf, or txt
You are on page 1of 4

Page 1 of 4

Three Days Workshop on Computational Fluid Dynamics 7th to 9th February 2014

Course Coordinator:

Prof. Vikas Lakhera Department of Mech. Engg. NIRMA University

Course Instructor:

Prof Atul Sharma Department of Mech. Engg. IIT Bombay

Lab # 1: Computational Heat Conduction Date: 08/02/20 ! Time: 2 P"# $%$0 P"

Note: &oth pro'lem an( ans)er sheet are pro*i(e( 'elo). To sa*e figure% +o to ,+raphi- )in(o) num'er./ -li-k on ,0ile./ then -li-k on ,12port to./ sele-t ,Win(o)s &"P image. in the ,0iles of type.. "ake sure to sa*e the file in the same lo-ation )here you ha*e this file. "ore (etails are gi*en in the ne2t page3 after the pro'lems.

. 1-D Computationa !eat Conduction "C!C# %


4onsi(er a D 4artesian -omputational (omain of si5e L=1 m/ for 464 transient simulation in a thin 7in x(ire-tion8 long steel plate 7(ensity% 99:0 kg/m$/ spe-ifi-#heat% :00 ;/<g </ thermal#-on(u-ti*ity% =.2 W/m#<. The plate is taken from a furna-e at a temperature of :0 04 7initial#-on(ition8 an( is su'>e-te( to *arious thermal 'oun(ary -on(itions an( *olumetri- heat generation3 sho)n in a ta'le 'elo). ?sing the program 7A @ D-on(u-tion.s-i8 for the a'o*e pro'lem/ run the -o(e for four (ifferent 'oun(ary#-on(itions 74ases8/ gi*en as% 4ases &oun(ary 4on(itions Volumetri- 6eat Left Aight +eneration 7W/m$8 A 00C 1000C 0 0 & 100 C 0 h=100W/m2.K, T=300C 0 0 4 0C 100 C 30,000 2 0 D 1000C 30,000 h=100W/m .K, T=30 C Take the ma2imum num'er of gri( points in 2#an( y#(ire-tion as max=12 an( -on*ergen-e -riteria as 0.00000 . Aeport the results as a8 Plot the stea(y state temperature profiles for the (ifferent -ases 7! figures8. '8 Dis-uss the effe-t of *olumetri- heat generation on the results for 'oth the types of &4s.

2. $# D Computationa !eat Conduction "C!C# %


4onsi(er a 2D 4artesian -omputational x-! (omain of si5e L=1 m an( "=1 m/ for 464 transient simulation in a long 7in #-(ire-tion8 steel plate 7(ensity% 99:0 kg/m$/ spe-ifi-#heat% :00 ;/<g </ thermal#-on(u-ti*ity% =.2 W/m#<. The plate is taken from a furna-e at a temperature of :0 04 7initial#-on(ition8 an( is su'>e-te( to *arious thermal 'oun(ary -on(itions an( *olumetri- heat generation3 sho)n in a ta'le 'elo). ?sing the program 7A @2D-on(u-tion.s-i8 for the a'o*e pro'lem/ run the -o(e for four (ifferent 'oun(ary#-on(itions 74ases8/ gi*en as% 4ases &oun(ary 4on(itions Volumetri- 6eat Left &ottom Aight Top +eneration 7W/m$8 A 1000C 2000C 3000C $000C 0 0 2 2 0 & 100 C Bnsulate( %W=10,000 W/m 0 h=100W/m .K, T=30 C 0 0 0 0 4 100 C 200 C 300 C $00 C &0,000 D 1000C Bnsulate( %W=10,000 W/m2 10,000 h=100W/m2.K, T=300C

Page $ of 4
Take the ma2imum num'er of gri( points in 2#an( y#(ire-tion as max='max=12 an( -on*ergen-e -riteria as 0.00000 . Aeport the results as a8 Plot the stea(y state temperature -ontours for the (ifferent -ases 7! figures8. '8 Dis-uss the effe-t of *olumetri- heat generation on the results for 'oth the types of &4s.

CDT1% 4DD1 1E14?TBDC BC S4BLA& To open a Scilab console% Linu2% Appli-ationsFProgrammingFS-ila' or Appli-ationsFS-ien-eFS-ila' Win(o)s% StartFAll programsFS-ila'FS-ila'.e2e To oad a source code 7or a (c note file8% Bn the S-ila' -onsole/ go to top menu 'ar 0ileFDpen a fileGF7&ro)se for the H.s-i file path8 The sour-e -o(e opens in a ne) (c note )in(o). N%T&: Dn-e this )in(o) is open/ su'seIuent loa(ing of a ne) H.s-i file )ill open a ne) ta' in the same (c note )in(o).

MU*TI,*& -CIN%T& TAB-

To e'ecute a pro(ram )ritten in Scinote% Bn the S-inote )in(o) 7)ith the (esire( program ta' open8/ go to top menu 'ar 12e-uteFGfile )ith no e-ho 7please (o not sele-t Gfile )ith e-ho8 4lear 4onsole Sele-t ,JesK The e2e-ution 'egins in the (c )a* con+o)e )in(o).

N%T&: %n y one Scilab code can be e'ecuted at a time. Dn-e an e2e-ute -omman( is gi*en an( the -o(e is to 'e stoppe( at an interme(iate stage/ use the metho( gi*en 'elo). Bf another e2e-ute -omman( is gi*en )ithout -ompleting/a'orting the pre*ious run/ erroneous results may 'e pro(u-e(. To abort a runnin( pro(ram% +o to the (c )a* con+o)e )in(o) an( press CTR*+C. This interrupts the -o(e e2e-ution. A prompt appears asking for user input. 1nter ,a'ort. here to stop the -o(e e2e-ution. <eep Playing )ith the DP1C#SD?A41 4o(es 7,e-e)ope, *! .rof. /t0) (harma, 11T 2om*a!8 in 0uture also.

Page . of 4

Answer Sheet
Pro'lem L % #D Computationa !eat Conduction "C!C# %
a8 Plot the stea(y#state temperature profile 7 T3x48 for the (ifferent -ases 7! figures8.

7a8

7'8

7-8

7(8

0ig. . % Stea(y state temperature profile T3x4 in the D (omain for 7a8 4ase A/ 7'8 4ase &/ 7-8 4ase 4 an( 7(8 4ase D '8 Dis-uss the effe-t of *olumetri- heat generation on the results for 'oth the types of &4s.

Write your ans)er/ limite( insi(e this te2t 'o2 only

Pro'lem L 2% 2#D Computationa !eat Conduction "C!C# %

Page 4 of 4
a8 Plot the stea(y state temperature -ontours for the (ifferent -ases 7! figures8.

7a8

7'8

7-8

7(8

0ig. .2% Stea(y state temperature -ontours in the 2D (omain for 7a8 4ase A/ 7'8 4ase &/ 7-8 4ase 4 an( 7(8 4ase D

'8 Dis-uss the effe-t of *olumetri- heat generation on the results for 'oth the types of &4s. Write your ans)er/ limite( insi(e this te2t 'o2 only

You might also like