Numerical homogenization of a heat transfer problem with phase transitions in heterogeneous frozen soils

広大な凍てつく北の大地において、地面は単なる土ではありません。それは、世界最大の都市やパイプライン、道路を支える構造的な基礎なのです。何世紀もの間、この永久凍土は固まった状態を維持し、人間の建設物にとって強固な岩盤として機能してきました。しかし、気候が温暖化するにつれ、この凍った地面が解け始めています。土壌内の氷が溶けると、地面は強度を失い、不等沈下や建物のひび割れを引き起こします。この融解が正確にどの深さまで進むのかを予測することは、安全とエンジニアリングに関わる問題ですが、土壌そのものが巨大な障害となっています。実際の土壌は均一な氷の塊ではなく、層やレンズ状の構造、異なる材料のポケットが入り混じった混沌とした混合物であり、その特性はわずか数センチメロ単位で変化します。このような複雑な地面の融解を完璧にシミュレートしようとすると、コンピュータはあらゆる層を追跡する必要があり、それは現在の技術では大規模かつ長期的な予測を行うには不十分な計算能力を必要とします。

これを解決するために、研究者は「数値ホモジナイゼーション（数理的均質化）」と呼ばれる手法を用います。スポンジの中を水が流れる様子を予測することを想像してみてください。スポンジのあらゆる小さな穴や繊維をモデル化する代わりに、スポadiaンジを、水を通す「平均的な」能力を持つ単一の材料として扱います。これがホモジナイゼーションの本質です。つまり、乱雑で詳細な現実を、コンピュータで高速に実行可能な、より単純化された平均モデルに置き換えるのです。しかし、地面が溶けているとき、この平均化は罠となります。凍結した土壌と融解した土壌の境界、すなわち「フロント（前線）」は、コンピュータの計算ブロックの端で綺麗に止まることはありません。それはブロックの真ん中を真っ二つに切り裂いて進みます。ここでジレンマが生じます。もし土壌のブロックが半分は凍っており、半分は溶けている場合、コンピュータはそのブロックの特性をどのように決定すべきでしょうか？ それは氷のように振る舞うべきか、水のように振る舞うべきか、あるいはその中間であるべきでしょうか？ この選択が、モデルが「建物が沈下する」と予測するか、「建物を維持できる」と予測するかを決定します。

最近の研究において、北東連邦大学のティムール・ザハロフ氏は、これら「混合状態」のブロックを処理するための3つの異なる方法をテストすることで、この特定のジレンマに取り組みました。研究は、凍結した地面の上に置かれた熱を持つ建物というシナリオに焦点を当てました。これは、構造物からの熱によって、その下の土壌がゆっくりと溶けていく状況です。チームは、自然界の混沌とした地面を模倣するために、土壌の熱伝導率に現実的なランダムパターンを用いた、この環境のデジタルツインを作成しました。そして、各ブロックの土壌が実際の層よりもはるかに大きい「粗い格子（コース・グリッド）」上でシミュレーションを実行しました。これにより、相変化のフロントがこれらのブロックの真ん中を切り裂く状況を作り出しました。彼らは、これら切り裂かれたブロックの中で何が起きているかを決定するための、3つの異なる手法を比較しました。

第一の手法は「テーブル・アプローチ」と呼ばれ、土壌ブロックのあらゆる可能な状態に対して特性のリストを事前計算しておき、ブロックの平均温度に基づいて答えを検索するものです。第二の手法である「ノード・スキーム」は、ブロックの角における土壌の状態を決定し、それらの決定を平均化します。第三の手法は「h-スムージング」と呼ばれ、線形温度勾配に基づいて、ブロックの溶けた正確な割合を計算する、つまり、ブロックのどの部分が水に変わったかを測定する手法です。研究者たちはこれらを二次元および三次元の両方で実行し、すべての層を追跡した非常に詳細な「リファレンス（参照）」解と比較しました。

結果は、エンジニアが道具をどのように選ぶべきかに挑戦する、驚くべきトレードオフを明らかにしました。目的が土壌の温度を予測することである場合、h-スムージング法が最も正確であることが証明されました。この手法は、他の手法が苦戦した粗い格子においても、熱の流れを極めて高い精度で追跡しました。しかし、研究者が融解フロント（凍結した地面と融解した地面の実際の境界）の位置に着目すると、物語は一変しました。h-スムージング法は、温度に関しては優れていたものの、融解フロントの深さを一貫して過大評価していました。つまり、地面が実際よりも深く溶けるという予測を出していたのです。

この不一致の理由は、エネルギー保存の法則にあります。h-スムージング法は、融解ゾーンの形状を精密に捉えようとするあまり、各タイムステップにおいてシステムに微量の余分な熱を意図せず加えてしまっていました。一年間の経過の中で、この微小な余剰が蓄積し、シミュレーション上の地面は現実の地面よりも多く溶けてしまったのです。対照的に、他の2つの手法、特にテーブル・アプローチとノード・スキームは「保守的」でした。これらはエネルギーを加算も減算もしなかったため、温度の読み取りはわずかに精度が劣るものの、融解の深さについてはるかに正確に予測できました。この研究は、もし研究者が土壌の正確な温度を知りたいのであれば、h-スムージングが優れた選択肢であることを示しました。しかし、もし目的が「融解による基礎の沈下深さ」を知ることであるならば、保守的な手法の方がはるかに信頼できるということです。

また、研究では、土壌が溶ける温度が場所によって一定ではなく、層ごとに異なる場合の影響についても調査されました。このようなより複雑なシナリオでは、特性の精緻なルックアップを使用するテーブルベースの手法が、二テンソル法よりも融解フロントの位置の予測において優れた性能を発揮しました。これは、高度に変動する地面条件においては、より詳細な事前計算済みの土壌状態マップが必要であることを示唆しています。最終的に、この研究はエンジニアや科学者に明確な指針を提供しています。すなわち、「唯一の最善の手法」は存在しないということです。選択は、どのような問いに答える必要があるかによって決まります。熱的な正確さを優先するなら一つの方法が勝ち、物理的な融解の境界を優先するなら別の方法がリードします。これらの違いを明確にすることで、この研究は、北極圏の建物を守るために使用されるモデルが、単に高速であるだけでなく、信頼できるものであることを保証する助けとなります。