<jupyter_start><jupyter_text>放射熱伝達率のテスト 無限平行面と2次元空間を比較する<jupyter_code># 表面温度の設定
theta_1 = 20.0
theta_2 = 15.0
# 面1の放射率
eps1 = 0.9

# 面2の放射率
eps2 = np.arange(1, 11) / 10.0
print(eps2)
# 無限平行平板の場合の有効放射率
infinite_emissivity = htc.effective_emissivity_parallel(eps1, eps2)
print(infinite_emissivity)
# 有限平行平板の場合の有効放射率
# 通気層厚さをパラメータに設定する
l_d = np.arange(15, 50) / 1000.0
print(l_d)

# 有効放射率を計算
finete_emissivity = [[htc.effective_emissivity_two_dimension(
    emissivity_1=eps1,
    emissivity_2=eps,
    l_d=l,
    l_s=0.455) for eps in eps2.tolist()] for l in l_d.tolist()]
# 放射熱伝達率の計算
infinite_rad_heat_transfer = [htc.radiative_heat_transfer_coefficient_detailed(
    theta_1=theta_1,
    theta_2=theta_2,
    effective_emissivity=emissivity) for emissivity in infinite_emissivity]
print(infinite_rad_heat_transfer)

finite_rad_heat_transfer = (htc.radiative_heat_transfer_coefficient_detailed(
    theta_1=theta_1,
    theta_2=theta_2,
    effective_emissivity=np.array(finete_emissivity)))
# グラフの描画（放射率別）
fig1 = plt.figure()
ax1 = fig1.subplots()

ax1.set_xlabel('Emissibity_2[-]')
ax1.set_ylabel('radiative heat transfer[W/(m2 K)]')

plt.grid(which='major', color='gray', linestyle='-')
plt.grid(which='minor', color='gray', linestyle='-')

plt.plot(eps2, infinite_rad_heat_transfer,linestyle='solid',
         color='r', label='infinite')
plt.plot(eps2, finite_rad_heat_transfer[0,:],linestyle='solid',
         color='g', label='finite thickness:0.015m')
plt.plot(eps2, finite_rad_heat_transfer[15,:],linestyle='solid',
         color='b', label='finite thickness:0.03m')
plt.plot(eps2, finite_rad_heat_transfer[34,:],linestyle='solid',
         color='y', label='finite thickness:0.049m')

plt.legend()
plt.show()
# グラフの描画（中空層厚さ別）
fig1 = plt.figure()
ax1 = fig1.subplots()

ax1.set_xlabel('thickness[m]')
ax1.set_ylabel('radiative heat transfer[W/(m2 K)]')

plt.grid(which='major', color='gray', linestyle='-')
plt.grid(which='minor', color='gray', linestyle='-')

plt.plot(l_d, np.full(35,infinite_rad_heat_transfer[0]),linestyle='solid',
         color='r', label='infinite emissivity=0.1')
plt.plot(l_d, np.full(35,infinite_rad_heat_transfer[4]),linestyle='solid',
         color='g', label='infinite emissivity=0.5')
plt.plot(l_d, np.full(35,infinite_rad_heat_transfer[9]),linestyle='solid',
         color='b', label='infinite emissivity=0.9')
plt.plot(l_d, finite_rad_heat_transfer[:,0],linestyle='dashed',
         color='r', label='finite emissivity=0.1')
plt.plot(l_d, finite_rad_heat_transfer[:,4],linestyle='dashed',
         color='g', label='finite emissivity=0.5')
plt.plot(l_d, finite_rad_heat_transfer[:,9],linestyle='dashed',
         color='b', label='finite emissivity=0.9')

plt.legend(bbox_to_anchor=(1.05, 1), loc='upper left')
plt.show()
# グラフの描画（中空層厚さ別）

#fig1 = plt.figure()
#ax1 = fig1.subplots()
# グラフサイズを指定する　デフォルト：figsize=(6.4, 4.8)、小サイズ：figsize=(4.0, 3.5)
fig, ax1 = plt.subplots(figsize=(4.5, 3.5))
    
# フォントサイズを設定
plt.rcParams["font.size"] = 14

ax1.set_xlabel('Thickness of vented air-layer[m]')
ax1.set_ylabel('Radiative heat transfer\n[W/(m2 K)]')

plt.grid(which='major', color='gray', linestyle='-')
plt.grid(which='minor', color='gray', linestyle='-')

plt.plot(l_d, np.full(35,infinite_rad_heat_transfer[0]),linestyle='solid',
         color='r', label='infinite emissivity=0.1')
plt.plot(l_d, np.full(35,infinite_rad_heat_transfer[4]),linestyle='solid',
         color='g', label='infinite emissivity=0.5')
plt.plot(l_d, np.full(35,infinite_rad_heat_transfer[9]),linestyle='solid',
         color='b', label='infinite emissivity=0.9')
plt.plot(l_d, finite_rad_heat_transfer[:,0],linestyle='dashed',
         color='r', label='finite emissivity=0.1')
plt.plot(l_d, finite_rad_heat_transfer[:,4],linestyle='dashed',
         color='g', label='finite emissivity=0.5')
plt.plot(l_d, finite_rad_heat_transfer[:,9],linestyle='dashed',
         color='b', label='finite emissivity=0.9')

# 軸の描画範囲設定
ax1.set_xlim(0.012, 0.052)
ax1.set_ylim(0, 6)
    
plt.legend(bbox_to_anchor=(1.05, 1), loc='upper left')
plt.show()<jupyter_output><empty_output>
{'fraction_non_alphanumeric': 0.17552019583843328, 'fraction_numerical': 0.04186046511627907, 'mean_word_length': 6.24468085106383, 'pattern_counts': {'":': 0, '<': 5, '<?xml version=': 0, '>': 5, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
{'chain_length': '1', 'dir': 'jupyter-structured-clean-dedup', 'id': '124113', 'license': 'MIT', 'n_tokens_mistral': 1895, 'n_tokens_neox': 1817, 'n_words': 217, 'path': 'valid_radiative_heat_transfer.ipynb', 'repo_name': 'BRI-EES-House/03_ventilation_layer', 'size': '4085'}