# Chad Baker # Created on 2011 Feb 10 # Distribution Modules import matplotlib.pyplot as plt import os, sys # User Defined Modules cmd_folder = os.path.dirname(os.path.abspath('../Modules/hx.py')) if cmd_folder not in sys.path: sys.path.insert(0, cmd_folder) import real_hx reload(real_hx) hx_val = real_hx.get_hx() k_gypsum = 0.17e-3 # thermal conductivity (kW / (m * K)) of gypsum board from Incropera # and DeWitt Intro. to Heat Transfer 5th ed., Table A.3 thickness_gypsum = 0.25 * 2.54e-2 # thickness (m) of gypsum board hx_val.R_extra = thickness_gypsum / k_gypsum hx_val.exh.T_inlet = 700. hx_val.cool.T_outlet = 300. hx_val.exh.mdot = 0.221 hx_val.cool.mdot = 0.5 hx_val.solve_hx() print "\nPlotting..."
# hx_exp.file = 'copper' fit_status = "_fit" # fit_status = '_exp' hx_exp.import_data() hx_exp.get_Qdot_fit() SIZE = 10 mdot_surf = np.linspace(hx_exp.exh.mdot.min(), hx_exp.exh.mdot.max(), SIZE) T_in_surf = np.linspace(hx_exp.exh.T_in.min(), hx_exp.exh.T_in.max(), SIZE) mdot2d, T_in2d = np.meshgrid(mdot_surf, T_in_surf) hx_exp.rep_Qdot_surf(mdot_surf, T_in_surf) hx_mod = real_hx.get_hx() if hx_exp.file == "gypsum": k_gypsum = 0.17e-3 # # thermal conductivity (kW / (m * K)) of gypsum board from Incropera # # and DeWitt Intro. to Heat Transfer 5th ed., Table A.3 thickness_gypsum = 0.25 * 2.54e-2 # thickness (m) of gypsum board hx_mod.R_extra = thickness_gypsum / k_gypsum elif hx_exp.file == "copper": # the following expressions for k come from Li and Peterson, 2006. The # effective thermal conductivity of wire screen k_copper = 400.0e-3 # thermal conductivity of copper (kW / (m * K)) # from Wolfram Alpha