伝熱学会北陸信越支部春季セミナーで研究成果を発表しました．

2022年5月14日（土）にオンラインで開催された日本伝熱学会北陸信越支部春季セミナーで研究成果を発表しました．

演題：熱・化学的非平衡流れのOpenFOAMシミュレーション
著者：坂村芳孝（登壇者）・中山勝之・大嶋元啓

Teamsでのオンライン会議でしたが，iPadのTeamsアプリの使用して画面をブロードキャストしながら，PCを使って発表しました．iPadを使うとレーザーポインター風のカーソルが使えて便利です．

謝辞：本研究で使用したソルバーhy2Foamの開発者・Vincent Casseau博士に深く感謝します．なお，本研究はJSPS科学研究費JP21K03880の助成を受けたものです．

日本機械学会北陸信越支部2022年合同講演会で本研究室の研究成果を発表しました．

2022年3月4〜5日にオンラインで開催された標記講演会におきまして，研究成果を発表しました．

演題：スペクトルフィッティング法による大気圧誘電体バリア放電プラズマ中の窒素分子の回転および振動温度の推算
著者：坂村芳孝（登壇者），大嶋元啓，Joshua Lim，小林龍司
発表日：2022年3月5日

オープンCAEシンポジウム2021で研究成果を発表しました．

オンラインで開催されたオープンCAEシンポジウム2021（2021年12月2日〜4日）において研究成果を発表しました．

講演題目：OpenFOAMを用いた熱化学的非平衡流れの数値シミュレーション
著者：坂村芳孝（登壇者）・中山勝之・大嶋元啓
発表日：2021年12月4日

本研究で使用したソルバーhy2Foamの開発者であり，GitHubリポジトリ上で情報を公開しているVincent Casseau博士に深く感謝します．なお，本研究はJSPS科学研究費JP21K03880の助成を受けたものです．

JSME第99期流体工学部門講演会で本研究室の研究成果を発表しました．

2021年11月9日（火）にオンラインで開催されている日本機械学会(JSME)第99期流体工学部門講演会で本研究室の研究成果を発表しました．

題目：OpenFOAMを用いた極超音速飛行体周りの熱化学的非平衡流の数値解析（Parkの2温度モデルに基づく数値解析コードhy2Foamの検証）
著者：坂村芳孝（発表者）・中山勝之・大嶋元啓

JSME北陸信越支部第58期総会・講演会で研究成果を発表しました．

2021年3月6日（土）にオンラインで開催された標記講演会で，われわれの研究室の成果を発表しました．

演題：OpenFOAMに実装された重合格子法による造波抵抗の予測精度
著者：後藤啓太（登壇者），中山勝之，大嶋元啓，坂村芳孝

12月に発表した内容をさらに発展させて，運動する物体にはたらく造波抗力の予測精度を検証しました．

なぜ気体の粘性係数は密度に依存しないのか？（気体分子運動論による説明）

はじめに

先日開催された「プレゼンテーション演習」の合同発表会で，本研究室のプレ配属生が次のような質問を受けていました．

火星の大気密度が小さいのであれば，密度に応じて粘性係数も変化するのでは？

私自身，温度が決まれば気体の粘性係数が決まることは，知識として知っていましたが，密度に対する依存性については意識の外にありました（ウッカリしていました）．

理科年表⁽¹⁾を紐解いても

気体の粘度（粘性係数　著者注）は数十Paより数気圧に至る広い範囲において圧力にはほとんど無関係である．

とは書かれていますが，その理由までは示されていません．

そこで，慌ててVincenti & Krugerのテキスト“Introduction to Physical Gas Dynamics”⁽²⁾を確認して整理してみましたので，備忘録として掲載しておきます．

輸送現象

気体の巨視的な特性量（速度，温度，濃度）が空間的に一様でないとき，一様な状態に戻す現象（粘性，熱伝導，拡散）が生じることは経験的によく知られているところです．長時間経つと，最終的に一様な状態に達し，これらの現象は消滅します．

これらの物理現象は，気体を構成する分子のランダムな運動（熱運動）によって運動量やエネルギー，質量が輸送されることで生じるものであり，輸送現象と呼ばれています．巨視的な特性量の空間的な偏りが，分子運動によって輸送される微視的な物理量の流れ（流束）にアンバランスをもたらし，巨視的な特性量の輸送を生み出すわけです．

本稿では気体の粘性についてのみ考えていきますが，この場合考察の対象となるのは気体の速度の非一様性です．例えば，速度の $x$ 方向成分 $u$ が $y$ 方向にだけ非一様に分布している場合， $y=$ 一定となる平面上には速度成分 $u$ を一様にする方向に単位面積あたり $\tau$ の大きさをもつせん断力（せん断応力）が働きます．これはニュートンの粘性法則として知られていて，次式のように表すことができます．

$\displaystyle \tau = \mu \frac{d u}{dy} \hfill (1)$

ここで現れた $\mu$ が粘性係数（粘度）と呼ばれる物性値です．次節では，気体分子運動論に基づき，粘性係数が密度や圧力に依らず温度だけで決まることを説明したいと思います．

気体分子運動論による輸送現象の説明

気体分子がもつ物理量（運動量，エネルギー，質量）の巨視的平均量 $\tilde{a}$ を考えましょう．ここでは $\tilde{a}$ は $y$ 方向にだけ変化するものとします．つまり

$\tilde{a} = \tilde{a}(y)$

図1に示されるように，平面 $y = y_0$ の下側に存在する分子がこの平面を横切るとき，この分子は，最後に他の分子と衝突した場所 $y = y_0 - \delta y$ でもつ平均量 $\tilde{a}(y_0 - \delta y)$ をもって移動すると考えます．ここで $\delta y$ は平均自由行程 $\lambda$ と同程度の大きさをもつものと仮定し， $\delta y = \alpha \lambda$ と表します．ここで， $\alpha$ は1程度の大きさをもつ比例定数です．平面 $y = y_0$ の上側に存在する分子は，上と同様に考えて，平均量 $\tilde{a}(y_0 + \delta y)$ をもって移動することになります．

また，平面を単位時間・単位面積あたりに横切る平均的な分子数は， $n \overline{C}$ に比例するものと考えられます．ここで， $n$ は分子数密度， $\overline{C}$ はランダムな分子運動の平均速さであり， $y = y_0$ で評価されたものです．

以上のことから， $y$ 軸の正の方向に向かって輸送される，単位時間・単位面積あたりの平均量の流束 $\Lambda_{\tilde{a}}$ は

$\Lambda_{\tilde{a}} = \eta n \overline{C} \left[ \tilde{a}(y_0-\alpha \lambda) - \tilde{a}(y_0+\alpha \lambda) \right]$

として計算することができます．ここで $\eta$ は比例定数です．

さらに， $\tilde{a}$ を $y = y_0$ のまわりでTayler展開して，1次項までを残すと次式が得られます．

$\displaystyle \Lambda_{\tilde{a}} = - \beta n \overline{C} \lambda \frac{d \tilde{a} }{dy} \hfill (2)$

ただし， $\beta = 2 \alpha \eta$ としました．

運動量の輸送と粘性

分子の質量を $m$ とすると，式(2)において $\tilde{a}=mu$ とおけば気体分子がもつ $x$ 方向の運動量流束が求まります．

$\displaystyle \Lambda_{mu} = - \beta m n \overline{C} \lambda \frac{d u }{dy} = - \beta \rho \overline{C} \lambda \frac{d u }{dy} \hfill (3)$

ここで $\rho = m n$ は気体の密度です．

次に，Newtonの運動法則

ある時間内に生じた物体の運動量の変化はその間に物体に作用した力積に等しい．

を適用すれば，図2に示すような平面に作用するせん断応力 $\tau$ は，この面を単位時間に通過する $x$ 方向の運動量流束と等しくなることがわかります．ただし，せん断応力の向きについては，図2に示した平面の場合， $x$ 軸の方向を正と定義していますので，この面を通過して流入する正味の運動量が正となる場合に対応しています．一方，前節で定義した運動量流束は $y$ 軸の正の向きが正となるように，つまり流出する方向が正となるように定義されていますので，次式のように符号を変更する必要があります．

$\tau = - \Lambda_{mu} \hfill (4)$

以上で求めてきた式(1), (3), (4)を比較することで，粘性係数が次のように与えられます．

$\mu = \beta \rho \overline{C} \lambda$

ここで，平均自由行程は密度に反比例⁽³⁾し，熱運動の平均速さが絶対温度の平方根に比例する⁽⁴⁾ことを思い出せば，粘性係数は密度には依存せず，温度のみの関数となることが理解できます．

おわりに

3年生の授業（「プレゼンテーション演習」合同発表会）での質問をきっかけにして，気体の粘性係数が密度に依存しないことを気体分子運動論的に改めて考えてみました．以上の内容が記載されていたVincentiとKrugerのテキストは研究室の輪読で何度か読んでいたのですが，最近取り上げることはなかったので，完全に記憶から抹消されていました．やはり定期的に読み返す必要がありますね．

なお，同書によれば，Maxwellが1860年にこの関係を発見したとき，あまりに驚いて，彼自身が実験で確認するまでの間，眠れなかったとのことです．

参考文献

国立天文台編，理科年表2019，丸善出版（2019）, p. 395.
W. G. Vincenti and C. H. Kruger, Jr., Introduction to Physical Gas Dynamics, Krieger Pub. Co. (1986), pp. 15-18.
ibid., p. 14.
ibid., p. 9.

流体工学シンポジウム2020で研究成果を発表しました．

2020年12月12日（土）にオンラインで開催された標記シンポジウム（北陸流体工学研究会主催）で本研究室の研究成果を発表しました．

演題：OpenFOAMの重合格子機能を用いて評価した造波抗力の精度について
著者：後藤啓太（登壇者）・中山勝之・大嶋元啓・坂村芳孝

小杉高校科学探求力アップ講座2020が終了しました！

2020年12月10日（木）(9〜12時）に標記講座が本学（L205会議室）で開催されました．

本講座は，高大連携教育の一環として毎年開催しているもので，今年で8回目となります．受講生は小杉高校の2年生（19名）で，11月16日（月）に小杉高校で行った授業と今回本学で実施した授業との2コマで構成される講座です．

本講座では，スターリングエンジンの動作原理を学んでもらった後に，エンジン模型の組み立てと動作実験を行いました．例年そうなのですが，組み立ててからすぐに動くことはほとんどなくて，皆さんとても苦労していました．それでも最後まで粘り強く取り組んでいる姿勢に，いつも感動させてもらっています（ありがとうございます！）．

最後に実施したコンテストの優勝タイムは，以下の通りです．優勝されたチームの皆さん，おめでとうございます．

優勝チーム：チーム上田（タイム：1分24秒）（参考　予選タイム9分30秒）

共有:

いいね:

共有:

いいね:

共有:

いいね:

共有:

いいね:

共有:

いいね:

はじめに

輸送現象

気体分子運動論による輸送現象の説明

運動量の輸送と粘性

おわりに

参考文献

共有:

いいね:

共有:

いいね:

共有:

いいね: