Tutorial 4

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
"""
This example shows how to simulate heat diffusion
in a hollow sphere made out of ATS58 (PCM material).
The sphere has an inner radius of 5 mm and an outer
radius of 25 mm. The initial temperature is 40℃.
The sphere is exposed to a convection heat transfer
with a heat transfer coefficient of 50 W/m²K. The
simulation time is 100 minutes, and the output is
saved every 10 minutes.
"""

import openterrace
import matplotlib.pyplot as plt
import numpy as np

def main():
    Ri = 0.005
    Ro = 0.025
    T_init = 40+273.15
    T_inf = 80+273.15
    h = 50

    ot = openterrace.Simulate(t_end=6000, dt=0.05)

    bed = ot.create_phase(n=30, type='bed')
    bed.select_substance(substance='ATS58')
    bed.select_domain_shape(domain='hollow_sphere_1d', Rinner=Ri, Router=Ro)
    bed.select_schemes(diff='central_difference_1d')
    bed.select_initial_conditions(T=T_init)
    bed.select_bc(bc_type='zero_gradient', 
                parameter='T',
                position=(slice(None, None, None), 0))
    bed.select_bc(bc_type='zero_gradient', 
                parameter='T',
                position=(slice(None, None, None), -1))

    # Initialise array of thermal resistances
    R = np.inf*np.ones_like(bed.T)

    # Set thermal resistance for the surface
    R[0][-1] = 1/(h*4*np.pi*Ro**2)

    bed.add_sourceterm_thermal_resistance(R=R, T_inf=T_inf)
    bed.select_output(times=range(0, 6000+600, 600))

    ot.run_simulation()

    plt.plot(bed.node_pos,bed.data.T[:,0,:].T-273.15, label=bed.data.time)
    plt.legend(title='Simulation time (s)')
    plt.show()
    plt.xlabel(u'Sphere radial position (m)')
    plt.ylabel(u'Temperature (℃)')
    plt.grid()
    plt.grid(which='major', color='#DDDDDD', linewidth=1)
    plt.grid(which='minor', color='#EEEEEE', linestyle=':', linewidth=0.8)
    plt.minorticks_on()
    plt.savefig('ot_plot_tutorial4.svg')

if __name__ == "__main__":
    main()