def T_gamma(T, C, shape, xx, yy, zz, t):
temp = np.ones(shape)
# Symmetric B.C.
temp[-1, :, :] = 0
temp[:, -1, :] = 0
T_shape = T.reshape(shape)
# Air facing edges compare to room temperature
temp[:, :, 1:] *= (1 - f_func(T_shape[:, :, 1:])) * \
c.h * (c.T_room - T_shape[:, :, 1:])
# Bottom facing edge compare to pan temperature
temp[:, :, 0] = (1-f_func(T_shape[:, :, 0])) * \
c.h * (c.T_pan - T_shape[:, :, 0])
# TODO: Check if h should be different for the pan (or sous vide)
return temp.ravel()
def T_gamma(T, C, shape, xx, yy, zz, t):
temp = np.ones(shape)
# Bottom has a different heat transfer than the rest
# temp[:,:,0] *= c.h_plate / c.h_air ## No; assume heat transfer is material independent
temp[-1, :, :] = 0
temp[:, -1, :] = 0
return temp.ravel() * (1 - f_func(T)) * c.h * (c.T_oven - T)
def C_gamma(T, C, shape, xx, yy, zz, t):
temp = np.ones(shape)
# No flux through bottom
temp[:, :, 0] = 0
# Symmetric B.C.
temp[-1, :, :] = 0
temp[:, -1, :] = 0
return temp.ravel() * f_func(T) * c.h * (c.T_oven - T)/(c.H_evap * c.rho_w) * (C - C_eq(T))