Tutorial 4

"""
This example shows how to simulate heat diffusion
in a hollow sphere made out of ATS58 (PCM material).
"""

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()