解説

重み付き残差式の右辺分布熱流束項は、 （式 1.17 参照）

$$\int_{S^{dist}} T^* {}^{t} q^{dist} \mathrm{d} S$$ (5.33)

要素 $Ie$ 境界表面 $Iface$ の要素分布熱流束項ベクトルは、 以下のように求まる。 なお、以下では、表面 $Iface = 0$ の場合を考える。

重み関数 $T^*$ の式 5.10 を分布熱流束項に代入すると、

$$\int_{S^{dist}} T^* {}^{t} q^{dist} \mathrm{d} S$$

$$= \sum_{Ie} \sum_{Iface} \int_{S_{(Ie Iface)}} { ( \langle \mathbf{ N } \rangle \{ \mathbf{ T^* } \} ^{(Ie)} ) } ^ { T } {}^{t} q^{dist} \mathrm{d} S$$

$$= \sum_{Ie} \sum_{Iface} \langle \mathbf{ T^* } \rangle ^{(Ie)} {}^{t} q^{dist} \int_{S_{(Ie Iface)}} \{ \mathbf{ N } \} \mathrm{d} S$$

$$= \langle \mathbf{ T^* } \rangle ( \sum_{Ie} \sum_{Iface} {}^{t} q^... ...t\{ \begin{array}{c} 0 \\ L1 \\ L2 \\ L3 \end{array} \right\} \mathrm{d} S )$$

$$= \langle \mathbf{ T^* } \rangle ( \sum_{Ie} \sum_{Iface} 1 / 3 S {}^{t} q^{dist} \left\{ \begin{array}{c} 0 \\ 1 \\ 1 \\ 1 \end{array} \right\} )$$ (5.34)