Scribd
Upload
39 views

Transient in A Bar by Explicit Method: T 0 T 0.00048 T 0.00096 T 0.00144 T 0.00192 T 0.0024 T 0.00512

The document describes an explicit method for modeling transient one-dimensional heat conduction in a bar over time. Three simulations are run with different time step sizes (Dt) to test for stability. In the first two cases, the solution remains stable as Dt is increased. However, in the third case the solution becomes unstable when Dt is largest, indicating the importance of satisfying the stability criterion when using explicit methods.

Uploaded by

ishanv_143
Copyright
© Attribution Non-Commercial (BY-NC)
Available Formats
Download as XLS, PDF, TXT or read online on Scribd
39 views

Transient in A Bar by Explicit Method: T 0 T 0.00048 T 0.00096 T 0.00144 T 0.00192 T 0.0024 T 0.00512

The document describes an explicit method for modeling transient one-dimensional heat conduction in a bar over time. Three simulations are run with different time step sizes (Dt) to test for stability. In the first two cases, the solution remains stable as Dt is increased. However, in the third case the solution becomes unstable when Dt is largest, indicating the importance of satisfying the stability criterion when using explicit methods.

Uploaded by

ishanv_143
Copyright
© Attribution Non-Commercial (BY-NC)
Available Formats
Download as XLS, PDF, TXT or read online on Scribd
You are on page 1/ 6

Transient in a bar by explicit method

t=0

t=0.00048

t=0.00096

t=0.00144

t=0.00192

t=0.0024

t=0.00512

1
0.9
0.8
0.7
0.6

T 0.5
0.4
0.3
0.2
0.1
0
0

0.05

0.1

0.15

0.2

0.25

0.3

0.35

0.4

0.45

0.5

Explicit Formulation for One-Dimensional Transient Conduction ia a bar:


Consider the equation: qxx = qt, for 0 < x < 1 and t > 0. At t = 0 q = 1 throughout.
For t > 0, q = 0 at x = 0 and x = 1. Determine the solution by using Dx = 0.1 and
Dt = 0.0025 , 0.005 and 0.0075 for three time steps. Comment on the results.
Use an explicit method.
In the explicit method, the temperature at any interior node is given by:
q(i,j+1) = q(i,j) + (Dt/Dx2)[q(i+1,j) - 2 q(i,j) + q(i-1,j)]
Case 1: DFo = (Dt/Dx ) = 0.0025/0.1 =
t
0
0.0025
0.005
i
x
T, j = 1
T, j = 2
1
0
1
0
2
0.1
1
1
3
0.2
1
1
4
0.3
1
1
5
0.4
1
1
6
0.5
1
1
Solution is stable
Case 2: DFo = (Dt/Dx2) = 0.005/0.12 =
t
0
0.005
0.01
i
x
T, j = 1
T, j = 2
1
0
1
0
2
0.1
1
1
3
0.2
1
1
4
0.3
1
1
5
0.4
1
1
6
0.5
1
1
Solution is stable
Case 3: DFo = (Dt/Dx2) = 0.0075/0.12 =
t
0
0.0075
0.015
i
x
T, j = 1
T, j = 2
1
0
1
0
2
0.1
1
1
3
0.2
1
1
4
0.3
1
1
5
0.4
1
1
6
0.5
1
1
Solution is unstable

0.25
0.01
T, j = 3
0
0.75
1
1
1
1

0.0125
T, j = 4
0
0.625
0.938
1
1
1

0.015
T, j = 5
0
0.547
0.875
0.984
1
1

0.0175
T, j = 6
0
0.492
0.820
0.961
0.996
1

0.02
T, j = 7
0
0.451
0.773
0.935
0.988
0.998

0.5
0.015
T, j = 3
0
0.5
1
1
1
1

0.02
T, j = 4
0
0.5
0.75
1
1
1

0.025
T, j = 5
0
0.375
0.75
0.875
1
1

0.03
T, j = 6
0
0.375
0.625
0.875
0.938
1

0.035
T, j = 7
0
0.313
0.625
0.781
0.938
0.938

0.75
0.0225
T, j = 3
0
0.25
1
1
1
1

0.03
T, j = 4
0
0.625
0.438
1
1
1

0.0375
T, j = 5
0
0.016
1.000
0.578
1
1

0.045
T, j = 6
0
0.742
-0.055
1.211
0.684
1

0.0525
T, j = 7
0
-0.412
1.492
-0.134
1.316
0.525

Case 2: DFo = (Dt/Dx ) = 0.005/0.1 =


t
0
50
100
i
x
T, j = 1
T, j = 2
1
0
100
100.0
2
0.1
30
65.0
3
0.2
30
30.0
4
0.3
30
30.0
5
0.4
30
30.0
6
0.5
30
30.0
Solution is stable

0.5
150
T, j = 3
100.0
65.0
47.5
30.0
30.0
30.0

200
T, j = 4
100.0
73.8
47.5
38.8
30.0
30.0

250
T, j = 5
100.0
73.8
56.3
38.8
34.4
30.0

300
T, j = 6
100.0
78.1
56.3
45.3
34.4
34.4

350
T, j = 7
100.0
78.1
61.7
45.3
39.8
34.4

DFo=0.25

A better solution with a larger number of nodes:


h=
0.025 Dt=
0.000156
x
0
0.025
0.05
0.075
0.1
0.125
0.15
0.175
0.2
0.225
0.25
0.275
0.3
0.325
0.35
0.375
0.4
0.425
0.45
0.475
0.5

400
T, j =8
100.0
80.9
61.7
50.8
39.8
39.8

0
0
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1

$j\downarrow
1
2
3
4
5
6
7
8
9

0.00016
0
0.75
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1

0.00032 0.00048 0.00064


0.0008 0.00096
0
0
0
0
0
0.625 0.546875 0.492188 0.451172 0.418945
0.9375
0.875 0.820313 0.773438 0.733154
1 0.984375 0.960938 0.93457 0.907715
1
1 0.996094 0.988281 0.977539
1
1
1 0.999023 0.996582
1
1
1
1 0.999756
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1

i\rightarrow$
1
2
65
30
100
47.5
100
65.0
100
69.4
100
73.8
100
75.9
100
78.1
100
79.5
100
80.9

3
30
30.0
38.8
47.5
51.9
56.3
59.0
61.7
63.6

4
30
30.0
30.0
34.4
38.8
42.0
45.3
47.8
50.2

5
30
30.0
30.0
30.0
32.2
34.4
36.6
38.8
40.5

6
30
30.0
30.0
30.0
30.0
31.1
32.2
33.3
34.4

0.00112
0
0.392761
0.698242
0.881531
0.964844
0.992615
0.999023
0.999939
1
1
1
1
1
1
1
1
1
1
1
1
1

30
30
30
30
30
30
30
30
30

0.00128
0
0.370941
0.667694
0.856537
0.950958
0.987274
0.99765
0.999725
0.999985
1
1
1
1
1
1
1
1
1
1
1
1

0.00144
0
0.352394
0.640717
0.832932
0.936432
0.980789
0.995575
0.999271
0.999924
0.999996
1
1
1
1
1
1
1
1
1
1
1

0.0016
0
0.336376
0.61669
0.810753
0.921646
0.973396
0.992803
0.99851
0.999779
0.999979
0.999999
1
1
1
1
1
1
1
1
1
1

0.00176
0
0.322361
0.595127
0.78996
0.90686
0.96531
0.989378
0.997401
0.999512
0.999934
0.999994
1
1
1
1
1
1
1
1
1
1

0.00192
0
0.309962
0.575644
0.770477
0.892248
0.956715
0.985367
0.995923
0.999089
0.999843
0.999981
0.999998
1
1
1
1
1
1
1
1
1

0.00208
0
0.298892
0.557932
0.752211
0.877922
0.947761
0.980843
0.994075
0.998486
0.999689
0.999951
0.999994
1
1
1
1
1
1
1
1
1

0.00224
0
0.288929
0.541742
0.735069
0.863954
0.938572
0.97588
0.99187
0.997684
0.999454
0.999896
0.999985
0.999998
1
1
1
1
1
1
1
1

0.0024
0
0.2799
0.52687
0.718958
0.850387
0.929244
0.970551
0.989326
0.996673
0.999122
0.999808
0.999966
0.999995
1
1
1
1
1
1
1
1

16
0
0.271668
0.51315
0.703794
0.837244
0.919857
0.964918
0.986469
0.995449
0.998681
0.999676
0.999934
0.999989
0.999999
1
1
1
1
1
1
1

17
0
0.264121
0.50044
0.689495
0.824535
0.910469
0.95904
0.983326
0.994012
0.998122
0.999492
0.999883
0.999978
0.999997
1
1
1
1
1
1
1

18
0
0.257171
0.488624
0.675991
0.812258
0.901128
0.952969
0.979926
0.992368
0.997437
0.999247
0.999809
0.999959
0.999993
0.999999
1
1
1
1
1
1

19
0
0.250741
0.477603
0.663216
0.800409
0.891871
0.946748
0.976297
0.990525
0.996622
0.998935
0.999706
0.99993
0.999986
0.999998
1
1
1
1
1
1

20
0
0.244771
0.467291
0.651111
0.788976
0.882725
0.940416
0.972467
0.988492
0.995676
0.99855
0.999569
0.999888
0.999975
0.999995
0.999999
1
1
1
1
1

21
0
0.239208
0.457616
0.639622
0.777947
0.873711
0.934006
0.968461
0.986282
0.994598
0.998086
0.999394
0.99983
0.999958
0.999991
0.999998
1
1
1
1
1

22
0
0.234008
0.448516
0.628702
0.767307
0.864844
0.927546
0.964302
0.983906
0.993391
0.997541
0.999176
0.999753
0.999934
0.999985
0.999997
0.999999
1
1
1
1

23
0
0.229133
0.439935
0.618307
0.75704
0.856135
0.921059
0.960014
0.981376
0.992057
0.996912
0.998911
0.999654
0.999902
0.999975
0.999994
0.999999
1
1
1
1

24
0
0.22455
0.431828
0.608397
0.74713
0.847592
0.914567
0.955616
0.978706
0.990601
0.996198
0.998597
0.99953
0.999858
0.999962
0.999991
0.999998
1
1
1
1

25
0
0.220232
0.424151
0.598938
0.737562
0.83922
0.908085
0.951126
0.975907
0.989026
0.995399
0.998231
0.999379
0.999802
0.999943
0.999985
0.999997
0.999999
1
1
1

26
0
0.216154
0.416868
0.589897
0.728321
0.831022
0.901629
0.946561
0.972992
0.98734
0.994514
0.99781
0.999198
0.999731
0.999918
0.999978
0.999994
0.999999
1
1
1

27
0
0.212294
0.409947
0.581246
0.71939
0.822999
0.895211
0.941936
0.969971
0.985546
0.993544
0.997333
0.998984
0.999645
0.999886
0.999967
0.999991
0.999998
1
1
1

28
0
0.208634
0.403358
0.572957
0.710756
0.81515
0.888839
0.937263
0.966856
0.983652
0.992492
0.996798
0.998736
0.99954
0.999846
0.999953
0.999987
0.999997
0.999999
1
1

29
0
0.205156
0.397077
0.565007
0.702405
0.807474
0.882523
0.932555
0.963657
0.981663
0.991359
0.996206
0.998453
0.999416
0.999796
0.999935
0.999981
0.999995
0.999999
1
1

30
0
0.201847
0.391079
0.557374
0.694323
0.799969
0.876269
0.927823
0.960383
0.979585
0.990147
0.995556
0.998132
0.99927
0.999736
0.999912
0.999973
0.999992
0.999998
1
1

31
0
0.198694
0.385345
0.550037
0.686497
0.792632
0.870082
0.923074
0.957043
0.977425
0.988859
0.994848
0.997772
0.999102
0.999663
0.999883
0.999962
0.999989
0.999997
0.999999
1

0.00512
0
0.195683
0.379855
0.542979
0.678916
0.785461
0.863968
0.918318
0.953647
0.975188
0.987497
0.994082
0.997374
0.99891
0.999578
0.999848
0.999949
0.999984
0.999996
0.999999
0.999999

0.00528
0
0.192805
0.374593
0.536182
0.671568
0.778451
0.857929
0.913563
0.9502
0.97288
0.986066
0.993259
0.996935
0.998693
0.999478
0.999806
0.999933
0.999978
0.999994
0.999998
0.999999

You might also like