このコンテンツを表示するには、JoVEへの購読が必要です。 または、無料トライアルをお申し込みください。

方法論記事

水平冷面の霜のCFDシミュレーション

574 閲覧数

DOI:

10.3791/68133

2025年9月12日

この記事について

サマリー

ここでは、リー相変化アプローチを使用したオイラー多相モデルを使用して、水平冷面での霜形成をシミュレートするための数値モデルを提示します。このモデルは、霜の体積率を動的に更新して密度の変化を捉え、霜の厚さ、密度、分布に関する実験データに対して検証されます。

要約

霜の形成は、冷凍、建設、天然ガス処理など、さまざまな分野で観察される一般的な現象です。ただし、その複雑な性質により、正確で信頼性の高い数値モデルの開発は依然として大きな課題です。この問題に取り組むためのこれまでの努力にもかかわらず、現在のモデルにはまだ一定の制限があります。本稿では、霜形成の根底にある基本的なメカニズムに基づいて開発された、霜形成の修正数値モデルを紹介します。このモデルは、リー相変化モデルと組み合わせたオイラー多相流アプローチを利用しています。さらに、最大霜の体積率を決定するためのアプローチが更新され、モデルが霜取りプロセス中の密度変動を考慮できるようになりました。このモデルは、さまざまな研究からの厚さ、密度、分布に関する実験データと比較することにより、厳密に検証されています。結果は、霜の厚さの平均絶対相対偏差 (MARD) が 8.97% であるのに対し、密度の MARD は 16.06% であることを示しています。さらに、モデルによって予測された霜の形態は、参考文献で報告された実験的観察とほぼ一致します。

概要

霜の形成は、さまざまな分野で観察される一般的な現象です。熱交換器の表面に霜が蓄積すると、熱伝達効率1 が著しく損なわれ、流体の流れが妨げられ、熱交換器2 の全体的な性能が損なわれ、最終的には熱交換器の正常な動作が妨げられます3。したがって、冷凍システムにおけるこの問題に対処するには、霜形成のメカニズムと挙動を理解することが重要です4。ここ数十年で、これらのシステムにおける霜形成の原因と特徴を調査するために、かなりの研究が行われてきました。

実験的研究により、霜の形成は、気温、湿度、冷たい表面の温度などのさまざまな要因の影響を受けることが実証されています5678910。多くの実験結果は、入ってくる気温が低いほど霜層が厚くなる傾向があり6、湿度レベルが高いほど霜層が密度が高くなる傾向があることを示しています7。水平面での霜の形成を研究し、冷たい表面の周期的な温度変化が霜層界面で融解を引き起こし、霜の形成速度、霜層の厚さ、および動的霜密度に大きな影響を与える可能性があることを発見しました8。他の研究では、霜の形態と分布の両方を調べています。Jeongらは実験で、最初に入り江の近くに霜が形成され、霜の丘9として知られる現象の発生につながることを観察しました。Noorshamsらは、水平円形のチューブ表面の霜の形成を調査し、円柱の前面と背面の霜層が上面の霜層よりも厚いことを発見しました10。さらに、いくつかの研究111213141516 は、理論的および経験的アプローチの両方とともに実験的な霜形成パターンを利用して、一次元の霜層の厚さを予測するモデルを開発しました。ジョーンズとパーカーは、分子拡散理論に基づいて霜の厚さの予測モデルを開発しました11。モデルと実験データの間の不一致は、30時間にわたって3%未満のままでした。計算技術の進歩に伴い、霜の形成をシミュレートするために数値流体力学 (CFD) に注目する研究者が増えています。従来の 1 次元モデルとは異なり、CFD シミュレーションは、特に霜の厚さの分布と温度プロファイルの視覚化において大きな利点をもたらします。Cuiらは、核形成理論に基づいて霜形成のCFDシミュレーションを実行しました12。霜の厚さの予測は、Lenicらによって提供された実験データからの13%未満の偏差を示しました13。並行して、構造化チューブ内の凝縮に関するCFD研究は、ディンプル14 やらせんピッチ15 などの幾何学的特徴が局所的な熱と物質の伝達を促進することを示しています。最近、Youら16 は、霜層を成長する多孔質媒体として特徴付け、蒸気拡散を直接組み込む動的メッシュベースのCFDモデルを開発し、低い計算コストを維持しながら5%未満の相対偏差を達成しました。これらの発見は、複雑な相変化現象の解決における CFD の可能性を強調し、霜形成モデリングに貴重な洞察を提供します。

結論として、かなりの数の研究5678910111213 が冷たい表面での霜の形成を調査しており、さまざまなパラメーターの下での霜の形成パターンの理解の進化に貢献しています。さまざまな次元やメカニズムを組み込んだ複数の数値モデルが開発されていますが、包括的な検証が欠けていることがよくあります。ほとんどの研究678910111213 は主に霜の厚さを使用してモデルを検証しており、より広範な適用性が制限されています16。これらの限界を克服するために、本稿では、霜形成の根底にあるコアメカニズムに重点を置き、リー相変化モデルとオイラー多相流モデルを統合した数値モデルを紹介します。さらに、霜密度の時間依存の変化を考慮して霜体積率の上限を計算する新しい方法が導入され、以前のモデルの欠点に対処します。提案モデルの精度と信頼性は、さまざまな実験条件にわたる霜の厚さ、密度変動、霜の形成パターンなど、さまざまな角度から評価されます。この広範な検証により、実際のアプリケーションにおける霜の挙動をより正確に予測するための堅牢な理論的枠組みが提供されます。

アクセスが制限されています。このコンテンツを表示するにはログインするか、トライアルを開始してください。

プロトコル

1. 物理モデルとメッシュ

  1. SpaceClaimを開き、[ スケッチ] タブを選択し、[作成]機能の[ 長方形 ]オプションを選択します。
  2. XOY 平面上に、X 軸に沿って長さが 500 mm、y 軸に沿って幅が 15 mm の 2D 幾何学モデルを作成します。
  3. ICEM を開き、[ファイル] タブに移動し、[ ジオメトリ ] タブを選択し、[ ジオメトリを開く ] を選択して 2D ジオメトリ モデルをインポートします。
  4. 2D モデルの境界を選択し、[部品] タブの [部品を作成] 機能を開き、さまざまな境界に名前を割り当てます。
  5. [ ブロック] タブを選択し、[ ブロックの作成] をクリックします。[ブロックを作成]ウィンドウで、[ パーツ名を継承 ]オプションをオンにします。次に、[ブロックの作成] タブに移動し、[ ブロックの初期化] を選択します。
  6. [ブロックの初期化]タブで、タイプとして [2Dサーフェスブロッキング ]を選択し、[ 設定で初期化]をオンにします。
  7. [選択]タブで、[方法]オプションで [フォームサーフェス ]を選択します。次に、「サーフェス」(Surfaces) で「 サーフェスを選択」(Select Surface(s) ) アイコンをクリックし、グラフィックスウィンドウで 2D 平面 を選択します。
  8. [サーフェス ブロッキング]タブで、[方法]オプションに[ ほとんどマッピング] を選択し、[フリー フェース メッシュ タイプ]に[ すべてのクワッド ]を選択し、[フリー フェース メッシュ方法]に[ ICEM CFD クワッド ]を選択します。
  9. [カーブをまたいでブロックをマージ]タブで、[方法]で [すべて ]を選択し、[サイズを無視]を 0.0 に設定します。
  10. [Pre-Mesh Params] タブの [Meshing Parameters] で、[ Edge Params] をクリックします。[Edge] で [ Select Edge(s)] アイコンをクリックし、グラフィカル インターフェイスで X 方向の Edge を選択します。
  11. [メッシュ パラメータ] で、[ノード] を 500、[間隔 1 から 1e+10]、[比率 1 から 2]、[間隔 2 から 1e+10]、[比率 2 から 2]、[最大スペース] を 1e+10 に設定します。次に、[ パラメータのコピー] オプションをオンにします。
  12. メッシュパラメータオプションで、メッシュ法則をバイジオメトリックに設定し、 適用をクリックします。[Pre-Mesh Params] タブの [Meshing Parameters] オプションで、[ Edge Params] をクリックします。
  13. メッシュパラメータ(Meshing Parameters)のエッジ(Edge)オプションで、 エッジの選択(Select Edge(s)) アイコンをクリックし、グラフィカルインターフェイスでY方向のエッジを選択します。
  14. [メッシュ パラメータ] で、[ノード] を 150、間隔 1 から 1e+10、比率 1 から 2、間隔 2 から 1e+10、比率 2 から 2、比率 2 から 2、最大スペースを 1e+10 に設定します。次に、[ パラメータのコピー] オプションをオンにします。
  15. [メッシュ パラメータ]オプションで、[メッシュ ロー]を[BiGeometric]に設定し、[適用]をクリックします。[ブロック]オプションで、[プレメッシュ]機能をアクティブにし、プロンプトが表示されたら[ はい ]をクリックします。
  16. Pre-Mesh を右クリックし、コンテキストメニューから Convert to Unstructured Mesh を選択します。
  17. [出力メッシュ]タブに移動し、[ ソルバー設定]をクリックします。ソルバー設定で、ソルバーとして ANSYS Fluent を選択し、 適用をクリックします。
  18. [出力] タブで、[ 入力の書き込み] をクリックします。新しいウィンドウ 1 で、[ 保存] をクリックします。新しいウィンドウ 2 で、[ はい]、[ 保存] の順にクリックします。新しいウィンドウ 3 で、[ 開く] をクリックします。をクリックし、新しいウィンドウ 4 で、[ 完了] をクリックします。

2. 霜形成シミュレーションソフトウェアの操作

  1. Ansys Fluentを開きます。[ファイル] タブに移動し、[読み取り] の下の [メッシュ] オプションを選択します。スケールメッシュに移動し、メッシュが作成された場所でmmを設定します。ここで紹介する手順は、表1のケース1で定義された実験条件に基づいています。
  2. ソルバー設定で、 圧力ベースの タイプ、 絶対速度 定式化、および 過渡時間を選択します。
  3. Y方向の重力加速度(Gravitational Acceleration in the Y-direction)を-9.81に設定します。[モデル]で、[エネルギー]をクリックし、[ エネルギー方程式]を有効にします。
  4. [モデル] で [Viscous ] を有効にし、[ k-epsilon (2 eqn)] を選択します。[k-epsilon モデル] セクションで、[ 標準] を選択します。
  5. [壁面近くの処理]で、[ 標準壁機能]を選択します。[乱流多相モデル] で [混合] を選択します。
  6. モデル定数オプションで、Cmu、C1-イプシロン、C2-イプシロン、TKEプラントル数、TDRプラントル数、分散プラントル数、エネルギープラントル数、壁プラントル数、乱流シュミット数をそれぞれ0.09、1.44、1.92、1、1.3、0.75、0.85、0.85、0.75、0.7に設定します。
  7. ユーザー定義関数で、乱流粘度混合物、フェーズ 1、フェーズ 2 をなしに設定します。[モデル] で [種 ] をクリックし、[ 種の輸送] を有効にします。
  8. [種輸送モデル] ウィンドウで、[ 拡散エネルギー源 ] をオンにして [ フェーズ 1] を選択します。
  9. [モデル] で [ マルチフェーズ] を有効にし、 オイラー モデルを選択します。「多相モデル」ウィンドウで、「 体積分率パラメーター」「定式化」「暗黙的」を選択し、「 オイラー相の数 」を 「2」に設定します。
  10. [多相モデル]ウィンドウで、[ フェーズ] タブを選択し、[ フェーズ-1 - プライマリフェーズ]を選択し、[ 名前 ]を [フェーズ-1]に設定し、[ フェーズ材料] [混合テンプレート]に設定します。
  11. [多相モデル]ウィンドウで、[ フェーズ] タブを選択し、[ フェーズ-2 - 二次フェーズ]を選択し、[ 名前 ]を フェーズ-2に設定し、[ フェーズ材料][氷]に設定して、 粒度を有効にします。
  12. [フェーズ]タブの[フェーズ設定]ウィンドウで、[ フェーズ プロパティ][粒度温度モデル]として選択し、[ 直径] 0.0001 に設定します。
  13. [相] タブの [粒状プロパティ] ウィンドウで、[粒状粘度] を 1e-05 に設定し、[粒状体積粘度] を 30 に設定し、[固体圧力] を lun-et-al に設定し、[粒状温度] を代数的に設定し、[摩擦粘度] を [なし] に設定し、[保圧限界] を [ユーザー定義] に設定し、[半径分布] を lun-et-al に設定し、[弾性係数] を [導出] に設定します。
  14. [多相モデル] ウィンドウの [フェーズ相互作用] タブで、[ フォース] タグを選択します。次に 、フェーズ 1、フェーズ 2 を選択し、係数を wen-yu に設定します。
  15. [多相モデル] ウィンドウの [フェーズ相互作用] タブで、[ フォース] タグを選択します。次に 、フェーズ 2、フェーズ 2 を選択し、返還係数を 0.9 に設定します。
  16. [多相モデル] ウィンドウの [相相互作用] タブで、[ 界面面積 ] タブを選択します。次に、 ia-symmetric オプションを選択します。
  17. [境界条件]で、[ 入口 ]をクリックし、[ フェーズ 1]タブを選択し、[フェーズ]オプションを[フェーズ 1]に設定します。
  18. フェーズ 1 タブの [速度入口] ウィンドウで、[ 運動量] タブを選択し、[速度指定方法] を [大きさ] に設定し、[垂直] を [境界] に設定し、[参照系] を [絶対] に設定し、[速度の大きさ] を 0.6 に設定します。
  19. フェーズ 1 の Velocity Inlet ウィンドウで、 Thermal タグを選択し、Temperature を 292.8 に設定します。phase-1タグのVelocity Inletウィンドウで、 Species タグを選択し、h2oを0.008202に設定します。
  20. [境界条件]で、[ 入口 ]をクリックし、 phase-2 タグを選択し、[フェーズ]オプションをphase-2に設定します。
  21. phase-2タグのVelocity Inletウィンドウで、 Momentum タグを選択し、Velocity Specification Methodを'Magnitude, Normal to Boundary'に設定し、Reference FrameをAbsoluteに設定し、Velocity Magnitudeを0に設定し、Granular Temperatureを0.0001に設定します。
  22. phase-2タグのVelocity Inletウィンドウで、 Thermal タブを選択し、Temperatureを273に設定します。phase-2タグのVelocity Inletウィンドウで、 Multiphase タブを選択し、Volume Fractionを0に設定します。
  23. [境界条件]で、[ アウトレット ]をクリックし、[ phase-1 ]タグを選択し、[フェーズオプション]を[phase-1]に設定します。
  24. phase-1タブの圧力出口ウィンドウで、 Thermal タブを選択し、Backflow Total Temperatureを300に設定します。
  25. phase-1タブの圧力出口ウィンドウで、 タブを選択し、h2oを0に設定します。[境界条件] で [ アウトレット ] をクリックし、[phase-2] タブを選択し、[Phase] オプションを phase-2 に設定します。
  26. phase-2 タブの Pressure Outlet ウィンドウで、 Thermal タブを選択し、Backflow Total Temperature を 300 に設定します。
  27. phase-2タブの圧力出口ウィンドウで、 Multiphase タブを選択し、Backflow Granular Temperatureを0.0001に設定し、Volume Fraction Specification MethodをBackflow Volume Fractionに設定し、Backflow Volume Fractionを0に設定します。
  28. 2.28 境界条件で、 をクリックし、 コールドウォール タブを選択し、フェーズオプションを混合に設定します。
  29. コールドウォールタブの壁ウィンドウで、運動 タブを選択し、壁のモーションオプションを固定壁に設定し、壁粗さモデルオプションを標準に設定し、粗さの高さを0に設定し、粗さ定数を0.5に設定します。
  30. コールドウォールタブの壁ウィンドウで、[ 熱] タブを選択し、[ 熱条件]として[温度]を選択し、[温度]を 252.65 に設定し、[材料]を[鋼]に設定します。
  31. [ソリューション] で [メソッド] をクリックし、[ソリューション メソッド] ウィンドウを開きます。[解法] ウィンドウで、[ 圧力-速度結合スキーム] として [位相結合 SIMPLE] を選択し、[ 勾配空間離散化] として [最小二乗セル ベース] を選択し、[ 圧力空間離散化] として [2 次] を選択し、[密度空間離 散化] として [1 次上風] を選択し、[ 運動量空間離散化] として [1 次上風] を選択し、[体 積分率空間離散化] として [1 次上風] を選択し、 として [1 次上風] を選択します乱 流運動エネルギー空間離散化、乱 流散逸率空間離散化として 1 次風上を選択し、 エネルギー空間離散化として風上第 1 次を選択し、 フェーズ 1 h2o 空間離散化として風上第 1 次上風を選択し、 過渡定式化として第 1 次陰的を選択します。
  32. [ソリューション] で [コントロール] をクリックし、[ソリューション コントロール] ウィンドウを開きます。[圧力緩和係数] を 0.4 に設定し、[密度緩和係数] を 1 に設定し、[身体力緩和係数過少] を 1 に設定し、[運動量緩和係数] を 0.4 に設定し、[体積分率緩和係数] を 0.4 に設定し、[粒状温度緩和係数] を 0.3 に設定し、[乱流運動エネルギー緩和係数] を 0.3 に設定し、[乱流散逸率緩和係数] を 0.3 に設定し、[乱流粘度緩和係数] を 0.3 に設定します。 Energy Under-Relaxation Factorを0.4に設定し、phase-1 h2o Under-Relaxation Factorを0.4に設定します。
  33. [ソリューション] で [初期化] をクリックし、[ソリューションの初期化] ウィンドウを開きます。初期化方法として 標準初期化 を選択し、参照フレームとしてセル ゾーンに相対的 を選択します。
  34. [溶液の初期化]ウィンドウで、[ ゲージ圧 ]を0に設定し、[乱流運動エネルギー]を0.00135に設定し、[乱流散逸率]を0.001143987に設定し、[phase-1 X Velocity]を0に設定し、phase-1 y速度を0に設定し、phase-1 h2oを0.008202に設定し、phase-2 X速度を0に設定し、phase-2 Y速度を0に設定し、phase-2 Volume Fractionを0に設定します。 phase-2 の粒度温度を 0.0001 に設定し、phase-2 の温度を 273 に設定します。次に、[ 初期化] をクリックします。

3. 後処理とデータエクスポートの構成

  1. 「結果」(Results) で「 等高線」(Contours ) をクリックして「等高線」(Contours) ウィンドウを開きます。Contoursウィンドウで、 Filled、Node Values、Boundary Values、Global Range、 および Auto Range オプションを有効にします。[タイプ]の[輪郭]として [フェーズ] を選択し、[ 体積率]を選択します。次に、フェーズとしてフェーズ 2 を選択し、[ 保存/表示] をクリックします。
  2. 「計算アクティビティ」で、「 ソリューション・アニメーション」 をクリックして、「アニメーション定義」ウィンドウを開きます。
  3. [アニメーション定義]ウィンドウで、[間隔の後に記録]を 1 に設定し、[ タイム ステップ]を選択します。[ストレージタイプ]として[ HSFファイル ]を選択します。アニメーション オブジェクト(Animation Object)オプションで contour-1 を選択し、[ OK]をクリックします。
  4. 結果(Results)でサー フェス(Surfaces)オプションを選択し、 新しいライン/レーキサーフェス(New Line/Rake Surface )オプションをクリックして、ライン/レーキサーフェス(Line/Rake Surface)ウィンドウを開きます。
  5. [ライン/レーキ サーフェス]ウィンドウで、[ ライン]を有効にし、[x0 [m]を0.21に設定し、[x1 [m]を0.21に設定し、[y0 [m]を0に設定し、y1 [m]を0.015に設定して、[ 作成]をクリックします。
  6. [ファイル] タブで、[エクスポート] の下の [ 計算中 ] オプションを選択します。次に、「 ソリューション・データ」 をクリックして「自動エクスポート」ウィンドウを開きます。
  7. [自動エクスポート] ウィンドウで、[ファイルの種類] として [ASCII] を選択します。[場所] として [セルの中心] を選択し、区切り文字として [スペース ] を選択します。 [データのエクスポート間隔] 1 に設定し、[ タイム ステップ] を選択します。「サーフェス」(Surfaces) オプションで、「 ライン 1」(line-1) を選択します。[数量]オプションで、[ 密度(フェーズ-2) ]と [体積分率(フェーズ-2)]を選択します。次に、[ 参照] をクリックし、[ファイルの選択] ウィンドウを開きます。[ファイルの選択] ウィンドウで、[ OK] をクリックします。[自動エクスポート] ウィンドウで、[ OK] をクリックします。
  8. 「計算アクティビティ」で、「 自動保存(フロー時間ごと) 」をクリックして、「自動保存」ウィンドウを開きます。次に、「データ・ファイルの保存間隔」を「100」に設定し、「 フロー時間」を選択し、「関連ケース・ファイルの保存タイプ」として 「変更された場合にのみ」 を選択して、「 OK」をクリックします。
  9. 「ソリューション」で、「 計算の実行 」をクリックして、「計算の実行」ウィンドウを開きます。次に、[時間の進行タイプ]として[固定]を選択し、[時間の進行方法]として [ユーザー指定 ]を選択します。
  10. [計算の実行] ウィンドウで、[時間ステップ数] を 7200 に設定し、[時間ステップ サイズ] を 1 に設定し、[最大反復/時間ステップ] を 20 に設定し、[レポート間隔] を 1 に設定し、[プロファイル更新間隔] を 2 に設定します。次に、[ 計算] をクリックします。

アクセスが制限されています。このコンテンツを表示するにはログインするか、トライアルを開始してください。

結果

提案された修正数値モデルは、霜形成の主要な特徴を効果的に捉えています。これは、霜の成長の根底にある基本的なメカニズムに基づいており、オイラー多相流アプローチとリー相変化モデルを組み合わせています。このアプローチにより、モデルは最大霜量率を更新できるため、霜取りプロセス全体の密度変動を考慮できます。シミュレーションは、霜の厚さ、密度、分布の予測をさまざまな研究の実験データと比較することにより、モデルが厳密に検証されていることを示しています。さらに、モデルによって予測された霜の形態は、文献で報告されている実験的観察と密接に一致しています。

図1は、シミュレートされた物理モデルの概略図を示しています。計算領域には、前面に完全に発達した入口領域が含まれ、上壁は断熱に設定されています。速度入口は左側の境界に定義され、右側の境界は圧力出口として機能します。この設定により、対流下での霜形成の現実的な境界条件が確保されます。

図2

アクセスが制限されています。このコンテンツを表示するにはログインするか、トライアルを開始してください。

ディスカッション

本研究では、低温水平冷地面において、時間と運転条件の両面に応じて体積分率の上限を動的に調整することで、霜の発生をシミュレーションし、霜密度の変動を再現できる数値モデルを開発します。ここで提示する検証は従来の冷地温度に限定されていますが、このモデルは、シミュレーションと実験の霜の厚さの間の最大偏差が-20%を超えないこと(MARD = 8.97%)、および霜密度の最大偏差が29.72%以内であることを保証します(MARD = 16.06%)。さらに、シミュレートされた霜の分布は、実験的観測値と密接に一致しています。これらの結果は、高精度の数値霜モデルの構築にとって非常に重要です。

文献11121314

アクセスが制限されています。このコンテンツを表示するにはログインするか、トライアルを開始してください。

開示事項

著者らは、この論文で報告された研究に影響を与えると思われる競合する経済的利益や個人的な関係が知られていないことを宣言します。

謝辞

この研究は、(XLYC2203184)、(U23A20657)、および(LJ222410153082)によってサポートされています。

アクセスが制限されています。このコンテンツを表示するにはログインするか、トライアルを開始してください。

材料

この記事で使用された材料の一覧
名前会社カタログ番号コメント
ぺらぺらアンシス
ICEMのダッソーシステム
スペースクレームアンシス

参考文献

  1. Saygin, A., Basol, A. M., Arik, M. An experimental study on the frost formation over a flat plate: Effect of frosting on heat transfer. Exp Therm Fluid Sci. 144, 110862(2023).
  2. Fang, X., et al. A new frictional pressure drop correlation based on flow patterns for hydrocarbon refrigerants condensation flow. Int J Refrig. 170, 214-223 (2025).
  3. Rong, X., et al. Experimental study on a multi-evaporator mutual defrosting system for air source heat pumps. Appl Energy. 332, 120528(2023).
  4. Jia, Y., Xu, X., Li, Y., Liang, X., Yao, M. Experimental studies on frost and defrost of fine tube bundles under coolant temperature between −20 and −5 °C. Int J Heat Mass Transf. 116, 617-620 (2018).
  5. Song, M., Dang, C. Review on the measurement and calculation of frost characteristics. Int. J. Heat Mass Transf. 124, 586-614 (2018).
  6. Lee, J., Lee, K. -S. The behavior of frost layer growth under conditions favorable for desublimation. Int J Heat Mass Transf. 120, 259-266 (2018).
  7. Lee, Y. B., Ro, S. T. Frost formation on a vertical plate in simultaneously developing flow. Exp Therm Fluid Sci. 26 (8), 939-945 (2002).
  8. Mengjie, S., Shangwen, L., Hosseini, S. H., Xiaoyan, L., Zhihua, W. An experimental study on the effect of horizontal cold plate surface temperature on frosting characteristics under natural convection. Appl Therm Eng. 211, 118416(2022).
  9. Jeong, H., Byun, S., Kim, D. R., Lee, K. S. Frost growth mechanism and its behavior under ultra-low temperature conditions. Int J Heat Mass Transf. 169, 120941(2021).
  10. Barzanoni, Y., Noorshams, O., Basirat Tabrizi, H., Damangir, E. Experimental investigation of frost formation on a horizontal cold cylinder under cross flow. Int J Refrig. 34 (4), 1174-1180 (2011).
  11. Jones, B. W., Parker, J. D. Frost formation with varying environmental parameters. J Heat Transf. 97 (2), 255-259 (1975).
  12. Cui, J., Li, W. Z., Liu, Y., Jiang, Z. Y. A new time- and space-dependent model for predicting frost formation. Appl Therm Eng. 31 (4), 447-457 (2011).
  13. Lenic, K., Trp, A., Frankovic, B. Transient two-dimensional model of frost formation on a fin-and-tube heat exchanger. Int J Heat Mass Transf. 52 (1-2), 22-32 (2009).
  14. Yu, J., Huo, R., Shen, H., Li, X., Zhu, Z. A simulation study on the condensation flow and thermal control characteristics of mixed refrigerant in a dimpled tube. Appl Therm Eng. 231, 120889(2023).
  15. Yu, J., Jiang, Y., Cai, W., Li, X., Zhu, Z. Condensation flow patterns and heat transfer correction for zeotropic hydrocarbon mixtures in a helically coiled tube. Int J Heat Mass Transf. 143, 118500(2019).
  16. You, Y., Wang, S., Lv, W., Chen, Y., Gross, U. A CFD model of frost formation based on dynamic meshes technique via secondary development of ANSYS Fluent. Int J Heat Fluid Flow. 89, 108807(2021).
  17. Cai, W., Fang, X., Li, S., Qiu, G. A modified CFD model for frosting on a horizontal plate. Int J Heat Mass Transf. 229, 125726(2024).
  18. Boyina, K. S., et al. Condensation frosting on meter-scale superhydrophobic and superhydrophilic heat exchangers. Int J Heat Mass Transf. 145, 118694(2019).

アクセスが制限されています。このコンテンツを表示するにはログインするか、トライアルを開始してください。

再版と許可

タグ

Lee