論 文

随伴変数法および有限要素法に基づく熱対流場における熱源境界形状の最適化
青木 崇1・倉橋 貴彦1・片峯 英次2・河田 剛毅3
1 長岡技術科学大学大学院 機械創造工学専攻（Department of Mechanical Engineering, Graduate School of Nagaoka University of Technology）
2 岐阜工業高等専門学校 機械工学科（Department of Mechanical Engineering, National Institute of Technology, Gifu College）
3 機械工学科（Department of Mechanical Engineering, National Institute of Technology, Nagaoka College）

Optimization of Heat Source Boundary Shape in Thermal Convection Field Based on
the Adjoint Variable and the Finite Element Methods
Takashi AOKI 1, Takahiko KURAHASHI 1, Eiji KATAMINE 2
and Yoshitaka KAWADA 3
AbstractIn this paper, we present a shape identification analysis about the cooling efficiency of a heater based on
the adjoint variable and the finite element methods. The purpose of this study is to obtain an identified
shape so as to approach the target radiation amount of the arbitrary boundary in the target shape. Also, theshape gradient for this identification problem is derived by the first variation of the Lagrange function.
Reshaping is accomplished using a traction method that is proposed as a solution to the domain optimization
problems. The traction method is an algorithm for non-parametric shape optimization problems, and thesmooth shape can be obtained by using this method. Based on these considerations, we have constructed aprogram code for shape identification analysis in thermal convection field. FreeFem++ was used tocalculate the identified shape.
Key Words : shape identification, finite element method, adjoint variable method, thermal convection field,
natural convection field

１．はじめに昇が問題となっており，そのために冷却効率を高くする，つまり放熱量を大きくすることが求められる．

電子機器や熱交換器などにおいて，性能改善を目

一般的に，電子機器においては温度上昇を抑えるた的とした形状設計は非常に重要である．例えば，電めに撹拌機などの冷却機構が設置されている．しか子機器においては発熱による電子機器本体の温度上し，機器のメンテナンス性やコンパクト性を考慮す
17

ると，そのような冷却機構の設置はあまり望ましく

𝑞𝑞(𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡) = ℎ��𝜃𝜃(𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡) − 𝜃𝜃�𝑓𝑓 �ないと言える．そこで，上記の性能を向上させるような例として，熱対流の中でも特に攪拌する機構などを必要としない，自然対流による冷却を利用したである．ここで，𝑄𝑄𝑖𝑖 は重み定数，𝑞𝑞 は各計算ステッ
プにおける放熱量，𝑞𝑞(𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡) は目標とする放熱量，ℎ� は熱伝達係数，𝜃𝜃 は温度，𝜃𝜃(𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡) は目標とする温機器の設計について取り上げる．
自然対流を利用した電子機器の例の 1 つとして，「クライストロン電源」と呼ばれるものがある．こ度，𝜃𝜃�𝑓𝑓 は外気温度を表している．𝑡𝑡0 , 𝑡𝑡𝑓𝑓 , 𝛤𝛤 はそれぞれについて，近藤らの 1)研究においては発熱する電れ数値解析における初期時刻，最終時刻，数値計算源をヒーターと見做し，ヒーターの位置や形状につ境界を表している．重み定数 𝑄𝑄𝑖𝑖 は，放熱量を規定いて実験的に検討を行っている．近年では，数値計する境界では 1 とし，他の境界では 0 として与える．

算によるシミュレーションが発展してきており，前述した検討は数理設計の立場から位置および形状を２．２決定するような形状決定問題として捉えることがで

支配方程式

非定常自然対流場における支配方程式として，きる．本研究では，ヒーターの形状設計に応用する

Boussinesq 近似された Navier-Stokes 方程式，連続のかたちで，自然対流場の形状設計問題について取り式，熱伝達方程式を考える．Navier-Stokes 方程式，上げる．これまでに，自然対流場の形状設計問題に連続の式は非圧縮性とする．これらの支配方程式は関して，トポロジー最適化における流速最大化問題
2)(3)

式(4)～(6)のように表される．

や温度分布を規定する問題 3)が行われている．しかし，これらの問題は定常問題，あるいは断熱境界1𝑢𝑢̇𝚤𝚤 + 𝑢𝑢𝑗𝑗 𝑢𝑢𝑖𝑖,𝑗𝑗 + 𝑝𝑝,𝑖𝑖 − 𝜈𝜈𝑢𝑢𝑖𝑖,𝑗𝑗𝑗𝑗 + 𝛽𝛽𝛽𝛽(𝜃𝜃 − 𝜃𝜃0 )𝑒𝑒𝑖𝑖
𝜌𝜌= 0 in 𝑡𝑡 ∈ �𝑡𝑡0 , 𝑡𝑡𝑓𝑓 � on 𝛺𝛺

を設計境界として取り扱う問題であった．
本論文では，ある任意の形状における境界の放熱量を目標値として設定し，熱源（ヒーター）の境界形状を変化させることで，境界の放熱量との 2 乗誤𝑢𝑢𝑖𝑖,𝑖𝑖 = 0 in 𝑡𝑡 ∈ �𝑡𝑡0 , 𝑡𝑡𝑓𝑓 � on 𝛺𝛺

差を最小化することを目的とした放熱量規定問題について検討する．まず，非定常自然対流場の支配方
程式を示し，随伴変数法および有限要素法を用いて，形状更新に必要な形状勾配関数を導出する．次に，導出した形状勾配関数に基づいて力法 4)を適用し，ヒーターの形状を決定する．

𝜃𝜃̇ + 𝑢𝑢𝑖𝑖 𝜃𝜃,𝑖𝑖 − 𝛼𝛼𝜃𝜃,𝑖𝑖𝑖𝑖 = 0 in 𝑡𝑡 ∈ �𝑡𝑡0 , 𝑡𝑡𝑓𝑓 � on 𝛺𝛺

(4)(5)(6)

ここで，𝑢𝑢𝑖𝑖 は流速，𝑝𝑝 は圧力である．𝜌𝜌, 𝜈𝜈, 𝛽𝛽, 𝑔𝑔, 𝜃𝜃0 , 𝛼𝛼はそれぞれ密度，動粘性係数，体積膨張率，重力加速度，基準温度，温度拡散率を表している．また，

２．放熱量規定問題の定式化
２．１

𝑒𝑒1 = 𝑒𝑒3 = 0, 𝑒𝑒2 = −1 である．

目的汎関数

２．３

放熱量規定問題の目的汎関数 𝐽𝐽 を，式(1)のように境界条件および初期条件

非定常自然対流場における境界条件，初期条件を定義する．この問題では，ある任意の境界における式(7)～(14)のように定義する．

放熱量を目標の放熱量に近づけることを考える．

𝐽𝐽 =1 𝑡𝑡𝑓𝑓2� � 𝑄𝑄 �𝑞𝑞 − 𝑞𝑞(𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡) � 𝑑𝑑𝛤𝛤𝑑𝑑𝑑𝑑
2 𝑡𝑡0 𝛤𝛤 𝑖𝑖(1)ただし，𝑞𝑞 = ℎ��𝜃𝜃 − 𝜃𝜃�𝑓𝑓 �𝑢𝑢𝑖𝑖 = 𝑢𝑢�𝑖𝑖 at 𝑡𝑡 = 𝑡𝑡0 in 𝛺𝛺(7)𝑢𝑢𝑖𝑖 = 𝑢𝑢�𝑖𝑖 in 𝑡𝑡 ∈ �𝑡𝑡0 , 𝑡𝑡𝑓𝑓 � on 𝛤𝛤1𝑢𝑢 and 𝛤𝛤𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑(8)

1
𝑡𝑡𝑖𝑖 = 𝑡𝑡𝑖𝑖̂ = �− 𝑝𝑝𝛿𝛿𝑖𝑖𝑖𝑖 + 𝑢𝑢𝑖𝑖,𝑗𝑗 � 𝑛𝑛𝑗𝑗
𝜌𝜌
in 𝑡𝑡 ∈ �𝑡𝑡0 , 𝑡𝑡𝑓𝑓 � on 𝛤𝛤2𝑢𝑢

(2)

18(9)

𝜃𝜃 = 𝜃𝜃� at 𝑡𝑡 = 𝑡𝑡0 in 𝛺𝛺𝜃𝜃 = 𝜃𝜃� in 𝑡𝑡 ∈ �𝑡𝑡0 , 𝑡𝑡𝑓𝑓 � on 𝛤𝛤1𝜃𝜃

𝑏𝑏 = 𝑏𝑏� = −𝛼𝛼𝜃𝜃,𝑖𝑖 𝑛𝑛𝑖𝑖 in 𝑡𝑡 ∈ �𝑡𝑡0 , 𝑡𝑡𝑓𝑓 � on 𝛤𝛤2𝜃𝜃

𝑞𝑞 = ℎ��𝜃𝜃 − 𝜃𝜃�𝑓𝑓 � = −𝛼𝛼𝜃𝜃,𝑖𝑖 𝑛𝑛𝑖𝑖
in 𝑡𝑡 ∈ �𝑡𝑡0 , 𝑡𝑡𝑓𝑓 � on 𝛤𝛤3𝜃𝜃
𝑝𝑝 = 𝑝𝑝̂ at 𝑡𝑡 = 𝑡𝑡0 in 𝛺𝛺

1(10)∗− 𝑢𝑢𝑖𝑖,𝑖𝑖= 0 in 𝑡𝑡 ∈ �𝑡𝑡0 , 𝑡𝑡𝑓𝑓 � on 𝛺𝛺(11)∗−𝜃𝜃 ∗̇ − (𝜃𝜃 ∗ 𝑢𝑢𝑖𝑖 ),𝑖𝑖 − 𝛼𝛼𝜃𝜃,𝑖𝑖𝑖𝑖+ 𝛽𝛽𝛽𝛽𝑢𝑢𝑖𝑖 ∗ 𝑒𝑒𝑖𝑖
𝑡𝑡+ ∫ 𝑓𝑓 ∫ 𝑄𝑄ℎ�2 �𝜃𝜃 − 𝜃𝜃(𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡𝑡) � 𝑑𝑑𝛤𝛤𝑑𝑑𝑑𝑑 = 0𝜌𝜌

𝑡𝑡0(12)

(18)

(13)

(19)

(14)

𝜃𝜃 ∗ = 0 at 𝑡𝑡 = 𝑡𝑡𝑓𝑓 in 𝛺𝛺

(20)𝑢𝑢𝑖𝑖 ∗ = 0 in 𝑡𝑡 ∈ �𝑡𝑡0 , 𝑡𝑡𝑓𝑓 � on 𝛤𝛤1𝑢𝑢 and 𝛤𝛤𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑(21)𝜃𝜃 ∗ = 0 in 𝑡𝑡 ∈ �𝑡𝑡0 , 𝑡𝑡𝑓𝑓 � on 𝛤𝛤1𝜃𝜃(22)いる．また，𝛤𝛤1𝑢𝑢 , 𝛤𝛤2𝑢𝑢 , 𝛤𝛤1𝜃𝜃 , 𝛤𝛤2𝜃𝜃 , 𝛤𝛤3𝜃𝜃 はそれぞれ流れ場

における Dirichlet 境界，流れ場における Neumann 境∗∗𝑡𝑡𝑖𝑖 ∗ = 𝑡𝑡𝑖𝑖̂ = �𝑢𝑢𝑗𝑗 𝑢𝑢𝑖𝑖 ∗ + 𝑝𝑝 ∗ 𝛿𝛿𝑖𝑖𝑖𝑖 + 𝜈𝜈𝑢𝑢𝑖𝑖,𝑗𝑗�𝑛𝑛𝑗𝑗
in 𝑡𝑡 ∈ �𝑡𝑡0 , 𝑡𝑡𝑓𝑓 � on 𝛤𝛤2𝑢𝑢

界，温度場における Dirichlet 境界，温度場におけるNeumann 境界，Robin 境界を表している．

�∗ = �𝜃𝜃 ∗ 𝑢𝑢𝑖𝑖 + 𝛼𝛼𝜃𝜃,𝑖𝑖∗ �𝑛𝑛𝑖𝑖
𝑏𝑏 ∗ = 𝑏𝑏
in 𝑡𝑡 ∈ �𝑡𝑡0 , 𝑡𝑡𝑓𝑓 � on 𝛤𝛤2𝜃𝜃

随伴方程式および形状勾配関数の導出

任意の随伴変数 𝑢𝑢𝑖𝑖 ∗ , 𝑝𝑝 ∗ , 𝜃𝜃 ∗ を用いることで，目的汎関数 𝐽𝐽 はLagrange関数 𝐽𝐽∗ に拡張され，式(15)で表される．

∗in 𝑡𝑡 ∈ �𝑡𝑡0 , 𝑡𝑡𝑓𝑓 � on 𝛺𝛺𝑢𝑢𝑖𝑖 ∗ = 0 at 𝑡𝑡 = 𝑡𝑡𝑓𝑓 in 𝛺𝛺� はそれぞれクロネッカーのデルタ，
ここで，𝛿𝛿𝑖𝑖𝑖𝑖 , 𝑛𝑛𝑗𝑗 , (∙)
境界における単位法線ベクトル，既知関数を表して２．４

𝛤𝛤(17)(23)(24)

さらに，設計境界の座標に対するLagrange関数の勾配 𝐺𝐺𝑖𝑖 は，式(25)から導出される．

𝐽𝐽 = 𝐽𝐽

𝑡𝑡𝑓𝑓

𝑡𝑡𝑓𝑓
1+ � � 𝑢𝑢𝑖𝑖 ∗ �𝑢𝑢̇𝚤𝚤 + 𝑢𝑢𝑗𝑗 𝑢𝑢𝑖𝑖,𝑗𝑗 + 𝑝𝑝,𝑖𝑖 − 𝜈𝜈𝑢𝑢𝑖𝑖,𝑗𝑗𝑗𝑗
𝜌𝜌𝑡𝑡0 𝛺𝛺𝑡𝑡𝑓𝑓

𝑡𝑡0+ 𝛽𝛽𝛽𝛽(𝜃𝜃 − 𝜃𝜃0 )𝑒𝑒𝑖𝑖 � 𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑

∗

𝛺𝛺𝑡𝑡0

𝛺𝛺

𝑡𝑡𝑓𝑓

𝑡𝑡0

𝛤𝛤𝑡𝑡𝑓𝑓𝛤𝛤= � (𝐹𝐹 𝑏𝑏 + 𝐹𝐹 𝑑𝑑 + 𝐻𝐻𝑏𝑏 + 𝐻𝐻𝑑𝑑 ) 𝑑𝑑𝑑𝑑
𝑡𝑡0+ � � 𝑝𝑝 𝑢𝑢𝑖𝑖,𝑖𝑖 𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑
𝑡𝑡0𝑡𝑡𝑓𝑓

� � 𝑡𝑡𝑖𝑖 ∗ 𝛿𝛿𝑢𝑢𝑖𝑖 𝑑𝑑𝛤𝛤𝑑𝑑𝑑𝑑 + � � 𝑏𝑏 ∗ 𝛿𝛿𝛿𝛿 𝑑𝑑𝛤𝛤𝑑𝑑𝑑𝑑

𝑡𝑡𝑓𝑓+ � � 𝜃𝜃 ∗ �𝜃𝜃̇ + 𝑢𝑢𝑖𝑖 𝜃𝜃,𝑖𝑖 − 𝛼𝛼𝜃𝜃,𝑖𝑖𝑖𝑖 � 𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑= � (𝐹𝐹 𝑏𝑏 + 𝐻𝐻𝑏𝑏 ) 𝑑𝑑𝑑𝑑 + 𝐺𝐺
𝑡𝑡0(15)(25)

ここで，𝐹𝐹 𝑏𝑏 , 𝐹𝐹 𝑑𝑑 , 𝐻𝐻𝑏𝑏 , 𝐻𝐻𝑑𝑑 , 𝐺𝐺 はそれぞれ式(26)～(30)で表される．

ここで，𝑢𝑢𝑖𝑖 ∗ は随伴流速，𝑝𝑝∗ は随伴圧力，𝜃𝜃 ∗ は随伴温度である．紙面の都合上，詳細は省略するが，Lagrange 関数の停留状態（𝛿𝛿𝐽𝐽∗ = 0）を考慮すると，随伴方程式および随伴条件は次のように決定される．

∗−𝑢𝑢𝚤𝚤̇ ∗ + 𝑢𝑢𝑗𝑗,𝑖𝑖 𝑢𝑢𝑗𝑗 − �𝑢𝑢𝑗𝑗 𝑢𝑢𝑖𝑖 �,𝑗𝑗 − 𝑝𝑝,𝑖𝑖∗ − 𝜈𝜈𝑢𝑢𝑖𝑖,𝑗𝑗𝑗𝑗∗∗+𝜃𝜃 ∗ 𝜃𝜃,𝑖𝑖 = 0 in 𝑡𝑡 ∈ �𝑡𝑡0 , 𝑡𝑡𝑓𝑓 � on 𝛺𝛺𝐹𝐹 𝑏𝑏 = �𝑡𝑡𝑖𝑖 ∗ 𝛿𝛿𝑢𝑢𝑖𝑖 𝑑𝑑𝛤𝛤(26)

𝐹𝐹 𝑑𝑑 = �𝑡𝑡𝑖𝑖 ∗ 𝑢𝑢𝑖𝑖,𝑗𝑗 𝛿𝛿𝑥𝑥𝑗𝑗 𝑑𝑑𝛤𝛤(27)

𝛤𝛤1𝑢𝑢 +𝛤𝛤2𝑢𝑢

𝛤𝛤𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑

(16)

19

Step8. 力法に基づいて形状更新を行い，Step 4に戻𝑏𝑏∗

𝐻𝐻 = �𝑏𝑏 𝛿𝛿𝛿𝛿 𝑑𝑑𝛤𝛤

(28)

る．

𝐻𝐻𝑑𝑑 = �𝑏𝑏 ∗ 𝜃𝜃,𝑖𝑖 𝛿𝛿𝑥𝑥𝑖𝑖 𝑑𝑑𝛤𝛤

(29)

３．数値解析例

�𝐺𝐺𝑗𝑗 𝛿𝛿𝑥𝑥𝑗𝑗 + 𝐺𝐺𝑖𝑖 𝛿𝛿𝑥𝑥𝑖𝑖 � 𝑑𝑑𝛤𝛤(30)

𝛤𝛤1𝜃𝜃 +𝛤𝛤2𝜃𝜃𝛤𝛤𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑

𝐺𝐺 = �

𝛤𝛤𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑

３．１

放熱量規定問題の解析モデルとして，図-1 のようなヒーターを設定した．本論文のはじめに述べた
ように，冷却効率の向上を目的とするヒーターの形

勾配 𝐺𝐺𝑖𝑖 は，線形弾性体の外力として考え，変位値状の一例として，赤線に示す長方形を初期形状とし，を修正勾配 𝐺𝐺𝑖𝑖 ∗ として用いる手法，力法を適用する

赤色破線に示す三角形を目標形状として設定した．

ことで変更される．式(31)に示すように，修正され初期形状および目標形状において，Robin 境界 𝛤𝛤3𝜃𝜃

た勾配 𝐺𝐺𝑖𝑖 ∗ を用いて，目標境界上の座標は更新されにおける目的汎関数を計算する．また，温度既知境界 𝛤𝛤1𝜃𝜃 を設計境界 𝛤𝛤𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑𝑑 として 𝜃𝜃 = 303.15 [K]を与る．ここで，𝑙𝑙 は形状更新ステップ回数，𝜂𝜂はステップ幅である．

𝑥𝑥𝑖𝑖 (𝑙𝑙+1) = 𝑥𝑥𝑖𝑖 (𝑙𝑙) − 𝜂𝜂𝐺𝐺𝑖𝑖
２．５

解析条件

え，すべての境界に対してノンスリップ境界条件∗ (𝑙𝑙)

𝑢𝑢𝑖𝑖 = 0 [m/s]を与えた．また，Neumann 境界 𝛤𝛤2𝜃𝜃 においては断熱境界として 𝑞𝑞 = 0 [W/m2]を与えた．

(31)

図-1 中の領域の寸法は，𝐿𝐿 = 1 [m], 𝑙𝑙 = 0.05 [m]である．表-1 に，図-1 の解析モデルにおける計算条計算の流れ件を示す．

本研究における形状最適化は，以下のステップを繰り返すことによって行われる．
Step1. 目標形状を与える．
Step2. 状態方程式を用いて，時間 𝑡𝑡 = 𝑡𝑡0 から 𝑡𝑡 = 𝑡𝑡𝑓𝑓の方向へ，流速分布 𝑢𝑢𝑖𝑖 ，圧力分布 𝑝𝑝 ，温度分布 𝜃𝜃 を解析する．このとき，目的汎関数の計算のために目標形状における温度分布 𝜃𝜃 を記録する．
Step3. 初期形状を与える．

Step4. 状態方程式を用いて，時間 𝑡𝑡 = 𝑡𝑡0 から 𝑡𝑡 = 𝑡𝑡𝑓𝑓の方向へ，流速分布 𝑢𝑢𝑖𝑖 ，圧力分布 𝑝𝑝 ，温度

図-1 数値解析モデル

収束判定条件式を満たした場合，解析を終了

表-1 数値解析条件

分布 𝜃𝜃 を解析する．
Step5. 目的汎関数の収束判定を行う．式(32)に示す
する．
�𝐽𝐽(𝑘𝑘) − 𝐽𝐽(𝑘𝑘−1) �⁄𝐽𝐽(0) < 𝜀𝜀(32)Step6. 上記のStep 4で得られた流速分布 𝑢𝑢𝑖𝑖 ，圧力分

布 𝑝𝑝 ，温度分布 𝜃𝜃 および随伴方程式を用いて，
時間 𝑡𝑡 = 𝑡𝑡𝑓𝑓 から 𝑡𝑡 = 𝑡𝑡0 の方向へ，随伴流速分

布 𝑢𝑢𝑖𝑖 ∗ ，随伴圧力分布 𝑝𝑝∗ ，随伴温度分布 𝜃𝜃 ∗を解析する．

Step7. これらの結果を用いて，形状勾配関数 𝐺𝐺𝑖𝑖 を導出する．

20

節点数
要素数
時間ステップ 𝑡𝑡𝑓𝑓 [s]
時間刻み幅 ∆𝑡𝑡 [s]収束条件値 𝜀𝜀 [-]動粘度 𝜈𝜈 [m2/s]
密度 𝜌𝜌 [kg/m3]
体積膨張率 𝛽𝛽 [1/K]
重力加速度 𝑔𝑔 [m/s2]
基準温度 𝜃𝜃0 [K]
温度拡散率 𝛼𝛼 [m2/s]

初期形状

目標形状

2769
5284

2599
4923

3000
3
10−7
1.16 × 10−2860
7.5 × 10−49.8
293.158.48 × 10−5

3.2

解析結果

析については今後の課題とする．

図-2，図-3，図-4 にそれぞれ初期形状における有限要素メッシュ，温度分布，流速・圧力分布，図-5，図-6，図-7 にそれぞれ目標形状における有限謝辞：本論文を執筆するにあたり，理化学研究所放要素メッシュ，温度分布，流速・圧力分布，図-8，射光科学研究センターの近藤力様の助言を受けた．

図-9，図-10 にそれぞれ同定形状における有限要素記して深く謝意を表する．

メッシュ，温度分布，流速・圧力分布を示す．図11 に目的汎関数の収束履歴を示す．形状更新の際には，FreeFem++上の関数「adaptmesh」を用いてリメッシュを行った．
図-5 と図-8 を比較すると，得られた同定形状と目標形状はおおよそ一致していることが確認できる．
図-6，図-7 と図-9，図-10 を比較すると，対流の発生によって得られた温度分布も一致するような結果が得られていることが確認できる．同定形状において，三角形の下部が膨らむような形状になっているが，これは領域内でスムーズに対流が起こりやすくなるように変化したからであると考えられる．また，
図-11 より，初期値を 100 [%]とした目的汎関数が十図-2 初期形状における有限要素メッシュ分に収束していることが確認できる．形状更新回数
100 回目付近で目的汎関数は既に収束しているように見えるが，図-12 において形状更新回数 100 回目と 150 回目の形状を比較すると，形状が変化していることが確認できる．そのため，目的汎関数の収束に伴う形状の変化も考慮して，適切な収束条件値 𝜀𝜀を決定する必要があると言える．なお，本解析では，
形状更新回数 232 回目において随伴解析が発散してしまったため，形状更新回数 231 回目で最適化解析を終了した．しかし，このまま形状更新を続けていくことができれば，より目標形状に近い形状が得られると考えられる．

図-3 初期形状における温度分布４．結論
本論文では，ある目標形状において放熱量を規定し，形状を同定する問題について取り組んだ．熱源を設計境界とした最適化問題として，予め目標とする熱源形状を設定し，最終的に得られた熱源形状との比較を行った．熱源を設計境界とした場合，場の状態が常に変化し続けることで目標とする放熱量も変化し続けてしまい，解の収束が安定しないことが考えられる．そのため，目標値を設定し，規定問題とすることによって，良好な結果が得られていることが確認できた．本検討では，Navier-Stokes 方程式に Boussinesq 近似を用いており，温度変化が小さい問題を対象としている．温度変化が大きい問題の解図-4 初期形状における圧力・流速分布21図-5 目標形状における有限要素メッシュ

図-8 同定形状における有限要素メッシュ

図-6 目標形状における温度分布

図-9 同定形状における温度分布

図-7 目標形状における圧力・流速分布

図-10 同定形状における圧力・流速分布22図-11 目的汎関数の収束履歴

図-12 形状更新 100 回目と 150 回目における形状の比較

参考文献
1) 近藤力 他，“クライストロン電源における絶縁油の冷却効率の向上”，第 4 回加速器学会，TP42，2003．
2) Joe Alexandersen, Niels Aage, Casper Schousboe Andreasen,Ole Sigmund, “Topology optimisation for natural convectionproblems”, International Journal for Numerical Methods inFluids, Vol.76, No.10, pp.699-721, 2014．
3) 今井伸哉，片峯英次，“温度分布を規定する非定常自然対流場の形状同定問題の解法”，日本機械学会論文
集，Vol.82，No.833，pp.15-00578，2016．
https://doi.org/10.1299/transjsme.15-00578.
4) 畔上秀幸，“領域最適化問題の一解法”，日本機械学会論文集 A 編，Vol.60，No.574，pp.1479-1486，1994．
(2020 8. 1 受付)

23