[論文レビュー] Refined radial basis function-generated finite difference analysis of non-Newtonian natural convection

Miha Rot, Gregor Kosec|arXiv (Cornell University)|Feb 16, 2022

Numerical methods in engineering被引用数 7

ひとこと要約

本稿では、2次元および3次元の複雑な幾何形状における非ニュートン対流のシミュレーションを目的とした、洗練された径数基底関数生成有限差分（RBF-FD）メッシュレス法を提示する。散乱ノードを用い、局所的勾配に基づく区分線形ノード密度関数による適応的グリッド細分化を採用することで、均一な離散化と比較して同等の精度を達成するためのノード数を最大5倍まで削減した。また、人工拡散を低減するための数値安定化手法を用いない。

ABSTRACT

Code and data for https://doi.org/10.1063/5.0257896 published in Physics of Fluids 37, 033130 (2025) Refined radial basis function-generated finite difference analysis of non-Newtonian natural convection In this paper, we present a refined radial basis function-generated finite difference solution for a non-Newtonian fluid in a closed differentially heated cavity. The non-Newtonian behavior is modeled with the Ostwald–de Waele power law and the buoyancy with the Boussinesq approximation. The problem domain is discretized with scattered nodes without any requirement for a topological relation between them. This allows a trivial generalization of the solution procedure to complex irregular three dimensional (3D) domains, which is also demonstrated by solving the problem in a two dimensional (2D) and 3D geometry mimicking a porous filter. The results in 2D are compared with two reference solutions that use the finite volume method in a conjunction with two different stabilization techniques, where we achieved good agreement with the reference data. The refinement is implemented on top of a dedicated meshless node positioning algorithm using piecewise linear node density function that ensures sufficient node density in the center of the domain while maximizing the node density in a boundary layer where the most intense dynamic is expected. The results show that with a refined approach, more than five times fewer nodes are required to obtain the results with the same accuracy compared to the regular discretization. The paper also discusses the convergence with refined discretization for different scenarios for up to 2e5 nodes, the impact of method parameters, the behavior of the flow in the boundary layer, the behavior of the viscosity, and the geometric flexibility of the proposed solution procedure.

研究の動機と目的

• 複雑な幾何形状における非ニュートン対流を、次元に依存しないメッシュレス数値的手法として開発すること。
• 計算領域におけるメッシュ生成およびトポロジカル接続の必要性を排除すること。
• 上流化やQUICKなどの安定化手法を避けることで、数値拡散を低減すること。
• 局所的なスティンプおよび基底関数制御を可能にすることで、効率的なhp適応性を実現すること。
• ベンチマーク用2次元および複雑な3次元多孔質フィルタ類似領域において、本手法の精度と幾何的柔軟性を実証すること。

提案手法

• 本手法は、トポロジカル制約なしに、散乱された非構造的ノードを用いた洗練されたRBF-FD離散化を採用する。
• 専用のメッシュレスノード配置アルゴリズムが、区分線形ノード密度関数を用い、境界層および勾配が高い領域にノードを集中させる。
• ナビエ–ストークス方程式にオストワルド＝ド・ワールの指数則モデルとボウシネスク近似を適用し、圧力・速度結合にはチョリンの投影法を用いる。
• 明示的時間刻み法を用い、物理的精度を保つために人工安定化（例：上流化やQUICK）を一切適用しない。
• ノードごとに変動するスティンプサイズ、多項式補完次数、RBFタイプを許容することで、hp適応性を実現する。
• 定式化は次元に依存せず、2次元から3次元へのスムーズな拡張が可能であり、コード変更を最小限に抑える。

実験結果

リサーチクエスチョン

• RQ1 洗練されたRBF-FDメッシュレス法は、数値安定化を用いずに非ニュートン対流において高い精度を達成できるか？
• RQ2 局所的勾配に基づく適応的ノード細分化は、均一なノード配置と比較して、どのように効率性と精度を向上させるか？
• RQ3 多孔質フィルタのような複雑な3次元幾何形状を、メッシュ生成なしにどの程度処理できるか？
• RQ4 せん断薄体挙動に起因する粘度変化が、閉じた領域内の流れ構造および速度プロファイルにどのように影響を与えるか？
• RQ5 スケーラブルで保守性の高いコードを実現するため、次元に依存しないオブジェクト指向的な実装が可能か？

主な発見

• 洗練されたRBF-FD法は、2次元ベンチマークケースにおいて、均一な離散化と同等の精度を維持しながら、最大5倍のノード数削減を達成した。
• 2次元ケースの結果は、上流化やQUICK安定化を用いた参考のFVM解と非常に良好に一致し、本手法の精度を裏付けた。
• 数値安定化を用いないことで、数値拡散が低減され、流れおよび温度勾配がより鋭い形状を保った。
• 3次元多孔質フィルターケースでは、高流量チャネルや壁および細い通路付近における強いせん断薄体効果といった複雑な流れ構造を的確に捉えた。
• せん断薄体に起因する粘度変化が明確に可視化され、最も低い粘度が壁付近および細い通路内などの高せん断領域に現れた。
• 本手法は強力な幾何的柔軟性を示し、最小限のコード変更で複雑な3次元領域の高精度なシミュレーションが可能となった。