エアリフトポンプの自由表面流解析
毛細管現象と土壌への透水現象の解析事例
水上ライドの解析事例
石油タンクのスロッシング現象
はんだの溶融解析事例
サンドセパレータの解析事例
防雪林の解析事例
噴霧エアノズルと噴霧燃焼の解析事例
噴霧・塗布ノズルと枚葉式洗浄機の解析事例
静電スプレーガンの解析事例
デフロスタの解析事例
洗濯物の浴室乾燥
天然氷の凝固・融解解析事例
冷凍ブロックの流水解凍

Analysis Case Studies

解析事例集−混相流編エアリフトポンプの自由表面流解析『STREAM®』を利用し、エアリフトポンプをMARS法で自由表面流解析
自由表面流解析

エアリフトポンプの自由表面流解析
界面捕獲法（Interface Capturing Method）

界面を表す関数を移流させることにより界面挙動を解析。
MAC（Marker and Cell）法、井戸や温泉の汲み上げ、浄化槽などの揚水で利用されるエアリフトポンプを、界面捕獲法であるMARS（Multi-interface Advection and Reconstruction Solver）
法で自由表面流解析を実行。

Level Set 法、VOF（Volume of Fluid）法、MARS 法など。

エアリフトポンプの仕組み
図1のように下部を水面下に設置した管（揚水管）に空気を流入させると、管内の水が空気と混合し、比重が小さくなり上に押し上げられる。
揚水量は空気流入量、
浸水深さ、揚水高さから経験式で求まるが、目的により曝気（エアレーション）
を促進する
場合と、促進しない場合があるため、揚水管内の気液二を把握する必要がある。
相流の流動様式（図2）
様々な制約により可視化実験ができないこともあり、流体シミュレーション解析で流動様式が把握できれば有効である。
図1 エアリフトポンプ

図2

管内の気液二相流の流動様式

解析結果

解析モデル

揚水管

送気管

揚水管

5 cm幅矩形

送気管

2 cm幅矩形

流入口

4 ヶ所、千鳥配置

空気量

流入口1 ヶ所あたり25 L/minで0.1秒毎に流入ポンプを1m浸水させて
水 面 から10 cm高 い 位
置の容器に水を汲み上
げる。

スラグ流
が発生

送気管先端の空気流入口

図3 解析したポンプ

図4 等値面図

図5 気液分布図

comment
気液界面を等値面として可視化しました（図4）。空気が混合した水が揚水され、水しぶきとなって容器に注水される様子が再現されています。図5は揚水管内の中心断面における気液分布図です。青色が水を表しています。
この図から、揚水管内がスラグ流だと推定することができます。

2 ｜ ソフトウェアクレイドル事例集【混相流編】

毛細管現象と土壌への透水現象の解析事例『STREAM®』を利用し、毛細管現象と土壌への透水現象をMARS法で解析
VOF（Volume of Fluid）法
F値の分布例

VOF法とは… ・解析領域の各要素に占める流体の

MARS法

ドナーアクセプター法

体積率をF値（0≦F≦1）
として定義
0.0

0.4

0.9

0.0

0.4

0.8

0.3

1.0

1.0

0.3

1.0

1.0

0.6

1.0

1.0

0.7

1.0

1.0

し、
この輸送方程式を解くことにより界面を求める自由表面流解析法。
・ F値が0であれば1番目の流体（例え

F 値を厳密
に保存

界面は矩形
で表現され
てしまう

ば空気）
、F値が1であれば2番目の
界面勾配を
1次関数で
近似する

流体
（例えば水）
として、それぞれの界面形状の再現性が高い

物性（密度や粘度など）を考慮する。

毛細管現象の解析事例
解析結果
接触角

毛細管現象により親水性壁面を有す

上昇る平板内では水面が上昇し、撥水性
壁面を有する平板内では水面が下降する。
このように接触角により流体の挙動は大きく変化する。

下降

図1 接触角

接触角とは… 壁面と流体の自由表面との角度（図1）。接 触 角が小さければ壁面は濡れ易い性質（親水性）、接
触角が大きければ壁は濡れ難い

図2 毛細管 平板の接触角60°の親水性壁面（左）平板の接触角120°の撥水性壁面（右）性質（撥水性）になる。

土壌への透水現象
解析モデル（土壌断面）
初期水位0.5 mm
凹み

砂地
30 x 70 mm

深さ

30 mm

幅

70 mm

中央部の凹み深さ1 mm

空隙率0.15（砂地、樹脂）

接触角°90（樹脂）

ダルシー係数

5 × 10-10 m2（空隙率を0とした場合の浸透率に相当）

初期条件

一様な高さ0.5 mmの水位（中央部の高さは1.5 mm）

図3 解析対象

comment

解析結果

図4は2秒後の解析結果です。見かけの透水速度は中央部で約10 ㎜/sになります。図5は砂地の保水性を高めるためにダルシー係数を2桁小さくした樹脂を深さ5 ㎜に一部設置した解析結果です。樹脂により透水が妨げられている様子が再図4 砂地への透水

図5 砂地への透水（樹脂設置）

現されています。

ソフトウェアクレイドル事例集【混相流編】 ｜3

水上ライドの解析事例『SCRYU/Tetra®』を利用し、水上ライドの移動をVOF法と重合格子で解析
VOF法と重合格子
・ 最新の流体シミュレーション解析ソフトウェアでは、VOF法と他の解析機能との組み合せが可能となり、物体の移動を考慮した自由表面流などが解析できるようになった。
・ 重合格子（overset grid）
は、移動領域と静止領域の各要素を重ね合せる方法。要素を再生成する必要がな
いため、
プログラムが複雑にならないことや計算安定性が維持できる。

図1 重合格子

解析内容
解析モデル
乗車人数

大人4名

イス

4脚（乗客の体重を考慮し1脚あたり86 kg）

直径

3m

密度
移動

360 kg/m3（下部円柱とフェンス）

6 DOF（6 Degrees of freedom）

水流からライドが受ける力を考慮し、並進および回転させる。
図2 曲がり流路

図3 ライド解析結果
速度の違いによって高低差が大きくなる

高さ位置[m]
移動速度

6 km/h
移動速度

7 km/h移動距離[m]
点線は流路の曲がり位置

図4 ライドの水上移動3秒後（速度6 ㎞/h）

図5 ライドの水上移動3秒後（速度7 ㎞/h）

図6 移動距離に伴うライドの高さ位置

comment
図4は水上移動するライドの速度が約6 ㎞/hの場合の解析結果です。
ライドは水上を揺れながら移動します。図5は水上移動するライドの速度が約7 ㎞/hの場合の解析結果です。
ライドは転覆こそしませんが、傾きが大きく乗客の安全が保障できません。図6は移動距離に伴うライドの高さ位置を表しています。
ライドの移動速度が6 ㎞/hの場合では10 mで0.54 m、7 ㎞/hの場合では10 mで1.2 mの高低差となります。解析結果から、ライドの移動速度が6 ㎞/hの場合でも、かなり大きなアップダウンのあるウォーターアトラクションであることが分かります。

4 ｜ ソフトウェアクレイドル事例集【混相流編】

石油タンクのスロッシング現象『STREAM®』を利用し、石油タンクのスロッシング現象を解析
石油タンクの解析事例
解析モデル
26,000 kL

石油貯蔵量（密度740 kg/m3、粘度0.0026 Pa·s、
表面張力係数0.020 N/m）
横揺れ・100 gal（震度5相当、加速度1 m/s2）
・10秒間継続・周期は長周期地震動として5秒

要素数

1,306,800（直交格子）

・MARS法（Multi-interface Advection and Reconstruction Solver法）で実時間60秒を解析。
・時間刻みはクーラン数0.9を超えないように自動設定。
図1 30,000 kl級石油タンク

計算時間の短縮を検討

A1(1) 対称面設定

30%
up

A2

(2) 並列計算

図のように対称面を設定
し要素数を653,400に軽減。
これにより計算時間は52 ％短縮。

(3) 部分的な要素粗化

2並列計算を行った場
合、計算時間はさらに
47 ％短縮。

(4) 1流体設定

液面の高低差が最大になる時の界面付近より下部の要素を粗く設定
すると、
要素数は
505,296となり、計算時間を23 ％短縮。

プログラムコードによっては1流体設定が可能。
スロッシング解析で1流体設定すると計算時間を30 ％短縮できる。

計算時間が4時間から約30分に86％ の短縮

初期設定から

解析結果
赤 線は初 期 の 液面高さから3.5 m
高い位置を表示

図2 最大液面高さ 周期5秒・13秒後（左）、周期4秒・15秒後（中央）、周期2秒・17秒後（右）comment
図2は100 galの横揺れ周期を変更して計算した結果です。左は周期5秒、中央は周期4秒、右は周期2秒の最大高さの図です。
周期5秒では、スロッシングの固有周期と近いために共振現象が起きていることがわかります。

ソフトウェアクレイドル事例集【混相流編】 ｜5

はんだの溶融解析事例『SCRYU/Tetra®』を利用し、リフローはんだ付けの作業プロセスを解析
はんだの溶融解析事例
解析モデル
基板上の2 ヵ所のランド部にクリームはんだを塗布し、VOF法を使って角形チップ抵抗器を配置。

はんだ

はんだ

基板

図1

ランド

チップ抵抗器長さ 0.4 mm
幅0.2 mm

溶接はんだの粘性係数

100 〜 0.020 〜 100 Pa·sに変化させることで溶融を表現

溶接はんだの密度

8000 kg/m3

表面張力0.40 N/m

ランドやチップ
抵抗器との接触角

30°0402サイズ チップ抵抗器• チップ抵抗器を囲む移動領域と、静止領域の各要素と重ね合せる重合格子（overset grid）
を設定。
• チップ抵抗器の移動は6DOF（6 Degrees of freedom）。運動方程式を解析する
ことにより、溶融はんだから受ける力を考慮して並進および回転。
図2 重合格子

解析結果
フィレット

溶接はんだ表面

チップ抵抗器の両端に溶融はんだ表面のフィレットが形成されていく様子を再現できた

ズレなし（ズレ幅小）
横ズレ条件1

comment
図3の横ズレ条件1の結果から、多少の横ズレならフィレット形式に問題ないことがわかります。横ズレ
条件2の結果のように、電子部品が立ってしまう現象は、米国ニューヨーク市のマンハッタン区の高層ビルに例えてマンハッタン現象と呼ばれます。

横ズレ条件2（ズレ幅大）
図3 溶接はんだ解析結果図 （ズレなし[上]、横ズレ条件1[中]、横ズレ条件2[下]）6 ｜ ソフトウェアクレイドル事例集【混相流編】

サンドセパレータの解析事例『SCRYU/Tetra®』を利用し、粒子追跡法で水と砂の分離現象を解析
粒子追跡法（Particle Tracking Method）

図1 Two way coupling

粒子追跡法は、図1のように格子で区切られた要素を用いて流体を連続相として解析（Euler法）
し、流体中の粒子を分散 相として個 別に追 跡（Lagrange法）
し、流体と粒子を連成解析する方法。

図2 粒子追跡法の適用

図3 流体中の粒子にかかる力

粒子追跡法では粒子を質点として考慮する。そのため、図2のように粒子の接触や衝突が支配的な対象には適用できない。

粒子密度が流体密度よりも大きい場合は抗力と重力が主な力になる。
ただし、
粒子の幾
何学形状を球として抗力を算出するため、注意が必要。
例えば、
必ずしも球形ではない
砂や埃、液滴などは形状が変化する。

サンドセパレータの解析事例
解析モデル

解析結果

図7 粒径による除去率の変化

図4 サンドセパレータ

サンドセパレータ
遠心力を利用して水と砂を分離する装 置。水と砂 の 密 度 差は約2.5で、空気と砂に比べると密度差が小さく分離し難い。

図5 粒子挙動

図6

流線

原水の流入量を2.8 ㎥ /hで一定と設定すると、約3 %がドレン水として流出する。
この原水と併せて粒径一定で密度2500 ㎏/㎥ の粒子100個を0.1秒毎に流入させ、30秒にわたり粒子追跡解析を行う。
ドレン水と共に流出する粒子数を、
ドレン水および処理水と共に流出する粒子数で除算して除去率とする。

comment
・図5は粒径100 ㎛の粒子挙動の解析結果です。
速度ベクトルを付与した粒子の挙動を表示しています。
原水と共に流入する粒子は壁面に沿いながら下方に移動し、ドレン水と共に流出します。
・図6は水の流線を矢印で表現しています。
サンドセパレータ内に生じる旋回流の様子がよくわかります。
・図7は粒径5、25、50、75、100 ㎛と粒径を変えた場合の除去率の変化です。粒径が小さくなるにつれて分離できなくなることがわかります。

ソフトウェアクレイドル事例集【混相流編】 ｜7

防雪林の解析事例『STREAM®』を使い、防雪林を粒子追跡法で解析
解析内容
道路域

防雪域

吹雪域

流入空気

流速8 m/s（基準高さ10 m）で氷点下20 ℃の空気がべき乗則（平坦地）で流入

雪粒子・密度200 ㎏/㎥ 、粒径100㎛
・地上から高さ3 mまで均一に流入・有効個数を1,000 としてパーセル近似

反発係数0.5（地面や防風柵など境界面に適用）

計算時間

8秒間

秒速8 m/s
防雪柵（開口率30%）
二車線道路

盛土

防風柵
防雪林

図1 道路防雪林

針葉樹の解析モデル
葉面積密度8 ㎡ /㎥

抵抗係数0.7 の圧損領域

雪粒子の反発係数0.1

針葉樹の配置

4 m間隔の2列の千鳥配置で狭林帯とする

針葉樹の設定・保育期（植樹から15年程度の期間）
・育成期（保育期からさらに15年程度の期間）
・維持期（植樹から30年以上）

図3 保育期（樹高4.5 〜 6.0 m）
防風柵2 m 、防雪柵2 m

図2 針葉樹（断面図）

図4

育成期（樹高6.5 〜 8.0 m）
防風柵2 m

図5 維持期（樹高9.0 〜 10 m）
防風柵2 m

解析結果

保育期

育成期

維持期

図6 雪粒子挙動

comment
防雪林が成長するに従い、防雪域での雪粒子数が増え、道路域での雪粒子数が減り、防雪効果が高くなっていることがわか
ります。
防災林なし

防災林あり
（維持期）
図7 道路視程

8 ｜ ソフトウェアクレイドル事例集【混相流編】

噴霧エアノズルと噴霧燃焼の解析事例『STREAM®』を利用し、粒子追跡法で蒸発や揮発を考慮した現象を解析
噴霧エアノズルの解析事例
粒子追跡法を用い、高温ガスを冷却する水滴の噴霧エアノズルを粒子追跡法で解析。
2パターンの噴霧条件を比較する。

水滴の蒸発モデル

解析結果
蒸発成分が生成される
ため、化学種を変数として取り扱う解析設定
流入する高温ガス

500 ℃

噴霧エア

27 ℃

形状

軸対称、軸対称2次元（冷却容器・ノズル）

高温ガス流速25 m/s

噴霧エア流速40 m/s壁面

静止壁、断熱条件

図1 噴霧エアノズル

図2 温度分布（7秒後）
フルコーン
（左）、
ホローコーン（右）燃料滴の噴霧燃焼の解析事例
燃料滴の噴霧燃焼を粒子追跡法で解析。

燃料滴の噴霧燃焼モデル

解析結果

燃料滴の噴霧燃焼モデル

揮発により流体中に生成する有機化合物は燃焼し、二酸化炭素と水蒸気が生成される

解析設定
燃料滴
噴霧流量
噴霧速度
噴霧パターン

ザウター平均粒径で20㎛（粒径分布に抜山・棚沢分布を適用）

0.005 kg/s
（パーセル近似10,000 s-1）

図3 燃焼容器

15 m/s
ホローコーン
広がり角110 〜 120°・1秒間は空気と噴霧エアだけを流入し、1秒後か
ら燃料滴を噴霧・燃料滴の潜熱は200 kJ/kg、Antoine式の定数はn‐デカンを用いる

図4

燃料滴挙動と温度分布（8秒後）

図5 燃焼ガスの質量分率分布（8秒後）有機化合物（左）、二酸化炭素（右）comment
上の噴霧エアノズルの解析結果から、
ホローコーンの方が比較的早く冷却されることがわかります。下の噴霧燃焼の解析結果は、燃料滴の噴霧燃焼により容器内の燃焼ガス温度が高くなることを示しています。

ソフトウェアクレイドル事例集【混相流編】 ｜9

噴霧・塗布ノズルと枚葉式洗浄機の解析事例『STREAM®』と『SCRYU/Tetra®』を活用し、粒子の液化挙動を解析
噴霧・塗布ノズルの解析（SCRYU/Tetra®）
解析結果

液膜モデル

板

50㎜ × 40㎜、
斜めに設置

噴霧粒子

粒径100㎛、ノズルで噴霧し板上で液化させ塗布している

塗布液膜の厚さ

1㎜オーダー、
粒子は長軸30°、短のフラット形
軸10°で噴霧している

粒子が壁面で液化して生じる液膜の物性（密度と粘度）と液膜の高さを考慮し、壁面上を移動する液膜を解析する方法。表
面張力や接触角を考慮しないため、計算負荷が小さい。

図1 解析結果（左：ノズル横からの視点

右：ノズル上部からの視点）

液体をほぼ均一に塗布できていることが確認できる半導体ウェハの枚葉式洗浄機の解析（STREAM®）
自由表面流解析と組み合わせる方法
粒子が壁面や液面に付着した場合に、粒子から液に変換する。液に変換した粒子は消滅させ、それ以上は追跡しない。表面張力や接触角まで考慮するため、液の分裂や結合を解析できるが、液膜モデルと比較して計算負荷が大きくなる。

解析結果
噴霧ノズル

粒径50 ㎛（1流体ノズル）

噴霧流量0.50 ㎏/s（30 L/ ）噴霧速度

15 m/s

物性

純水（洗浄水）接触角50°

図2 枚葉式洗浄機

図3

噴霧の様子

円錐噴霧広がり角30（左）° 、5（右）°楕円錐噴霧の場合は、洗浄効率を高めるためにウェハを高速回転すると、洗浄水が周方向にも広がることが確認できる
図4 洗浄水の様子（100秒後） 円錐噴霧（左）、楕円錐噴霧（中）、楕円錐噴霧2500 rpm（右）comment
STREAM®とSCRYU/Tetra®を用い、粒子追跡法で粒子の液化を再現することができました。

10 ｜ ソフトウェアクレイドル事例集【混相流編】

静電スプレーガンの解析事例『STREAM®』を利用し、静電スプレーガンを粒子追跡法で解析
静電スプレーガンの解析事例
解析モデル
塗装ブース

500 ㎜ × 300 ㎜ × 250 ㎜
・幅20 ㎜ × 高さ200 ㎜ の20面体柱

20面体柱・静電スプレーガン先端との距離が100 ㎜となる位置に置いて塗装（接地電位）
・電位は0 V

空気速度・ブース天井面から一様の速度0.3 m/sで流入塗着されない塗料粒子を排気（グレーチングのような床面）
・ブース床面は排気面

電束密度
空気の比誘電率

0 C/㎡（上下面を含む塗装ブースの全ての壁面）
空気の比誘電率 1.000586

図1 塗装ブース

解析結果
電荷なし

図2 静電スプレーガン

図3 中央断面での電位分布

空気流量40 L/min（ノズル径10 ㎜）

塗料吐出量

100 g/min（塗料密度は1000 ㎏/㎥）

塗料粒子の粒径50 ㎛

ノズル先端電位-50,000 V

電荷あり

20面体柱に付着した塗料粒子は消滅させ、それ以上追跡しない
図4 解析結果（0.05秒後） 正面の塗膜分布（上）、背面の塗膜分布（下）堆積させ塗膜の厚さに換算する

comment
20面体柱に付着した塗料粒子数と噴霧した塗料粒子数から塗着効率を算出します。塗料粒子の電荷がない場合は59.6 %でしたが、塗料粒子の電荷がある場合は84.5 %となり、静電塗装の効果が再現されました。

ソフトウェアクレイドル事例集【混相流編】 ｜11

デフロスタの解析事例『STREAM®』
で室内の窓や壁で生じる結露・蒸発を計算
電熱式デフロスタの解析事例
デフロスタとは…

車両の窓ガラスの結露による曇りを除去する機能

結露（蒸発）解析とは… ・気相側の流体解析で水蒸気を変数として取り扱う
・図1のように、空気中の絶対湿度が壁面温度の飽和絶
対湿度より大きい場合には結露速度[ kg/(m2 s) ]が生じ、
潜熱分2500 kJ/kgを発熱する
・空気中の絶対湿度が壁面温度の飽和絶対湿度より小さい場合は蒸発速度（負の結露速度）
が生じ、
潜熱分を
吸熱する

図1 結露速度解析モデル

運転席ガラス枠

運転席ガラス
車輛外側-5 ℃

運転手36 ℃
水分発生量

電熱線

100 W/本・二重ガラス構造と仮定。

運転席扉
20 ℃

電熱式デフロスタ

50 mL/h

枠ガラス

運転席
初期温度 20 ℃
初期温度 55 %（相対湿度）

運転車両外観

電熱式デフロスタ・電熱線を4本として開始時から90秒後に発熱し、10秒間で１本あたり100 Wの発熱 量に到 達、そ の 後はその発熱量を維持するもの
とする。

運転席

図2 鉄道車両

解析結果
計算開始10秒後

90秒後

300秒後

900秒後

解析対象

運転手席

表面温度

36 ℃

水分発生量

9.0 × 10-6 kg/(m2 s)（50 mL/h）

表面温度

・−5 ℃
（車両外側運転席ガラス）（運転席扉）
・20 ℃
・その他は断熱条件と想定

初期温度

20 ℃（運転席）

初期湿度

55 ％（相対湿度）
温度差による浮力を考慮

図3 運転手ガラス面の温度分布（上）と結露量分布（下）comment
開始時から90秒後には運転席ガラス面に結露が生じますが、その後は電熱線の発熱に伴い、結露が除去されていく様子が確認できました。

12 ｜ ソフトウェアクレイドル事例集【混相流編】

洗濯物の浴室乾燥『STREAM®』
で固体の吸放湿特性を考慮した結露・蒸発解析を実行
吸放湿特性の設定
吸放湿特性とは…
固体の湿度（絶対湿度）を変化させるのに必要な水分量と、固体の温度を変化させるのに必
要な水分量。

図1 固体中の湿気

図2

固体中の吸放湿図1のように乾燥時に空隙のある固体（多孔質体）を考える。固体中の湿気は水分の浸透や水蒸気の拡散で移動する。
同時に、図2のように固体中では水蒸気が凝縮して水分になり、水分が蒸発して水蒸気になる相変化が起きる。相変化に伴う潜熱分の発熱と吸熱が影響するので、固体中の湿気と熱の移動は連成して解析する必要がある。
この連成解析を実施する場合に吸放湿特性を設定する。

洗濯物の浴室乾燥の解析事例
解析結果

解析モデル

洗濯物表面と中心断面での温度分布

図3 洗濯物の浴室乾燥
洗濯物

長さ126 cm×幅60 cm×厚さ1cm
天井近傍の横棒にかけてあると仮定

空隙率0.4（乾燥時）密度1.30 kg/m3

比熱1.15 kJ/(kg·K)

熱伝導度0.086 W/(m·K)

湿気伝導度

5.0×10-6 kg/(m·s (kg/kgDA))*1

絶対湿度を変化させるのに必要な水分量

5.0×103 kg/(m3 (kg/kgDA))*1

温度を変化させるのに必要な水分量2.5 kg/(m3 K)*1

DA…Dry Airの略(kg/kgDA)…絶対湿度の単位

洗濯物表面での相対湿度分布

吹出口から吹き出す空気流量50m3/h

吹出口から吹き出す空気流量60m3/h
図4

解析結果（1時間30分後） 吹出口からの空気流量50m3/h（上）、60m3/h（下）

comment
浴室乾燥により洗濯物が乾いていく様子が再現されています。空気流量を大きくすることにより洗濯物の乾燥が早まることがわかります。

ソフトウェアクレイドル事例集【混相流編】 ｜13

天然氷の凝固・融解解析事例『STREAM®』を使い、凝固・融解解析をマクロにアプローチ
温度回復法
・流体中における固相の体積率を固相率と定義する。
・固液界面では平衡であると仮定し、温度回復法により固相率の変化を求める。
・図1のように流体要素の温度を解き、温度が液相線温度（水などの純物質は平衡であれば固相線温度と一致する）
よりも低い場合は潜熱と比熱から固相率を算出する。
図1 温度回復法

・次に、固相率に相当する潜熱を放出（発熱）
して再び流体要素の温度を解くこと
により、流体要素の温度が回復する。
・これを反復すると流体要素の固相率と温度が求まる。

天然氷の凝固・融解解析
解析結果

解析モデル

水面温度分布（2 時間後）

固相率分布（20 日後）

・淡水： 固相率分布から、水面から氷が成長していく様子がわかる。

図2 水池

水池

長さ 24 m×幅 14 m×深さ 0.5 m
縁石で囲まれ、その周囲は有機質土

池の上面

外気 −8 ℃
熱伝達係数 10 W/(m2K)

池の下面

外気 4 ℃( 熱伝導条件 )

解析手法

非定常解析

備考

水の流れは解かない（質量分率）の混合液で −1.8 ℃で水分が凍り始・海水（塩水）： 海水は塩分濃度が3.4 %める（液相線温度）。氷の成長はなくなり、塩分による凝固点降下が再現されている。
※赤線は水面から深さ15 cmの位置
図3

解析結果

淡水（上）、海水（下）comment
図3は2時間後の水面温度分布と20日後の固相率分布です。淡水の結果から、水面はほぼ均一に温度が下がり、水面から氷が成長していく様子がわかります。水池を海水のような塩水に変えると、氷の成長はなくなり、塩分による凝固点降下が再現されています。

14 ｜ ソフトウェアクレイドル事例集【混相流編】

冷凍ブロックの流水解凍『STREAM®』を使い、流水の解凍現象を凝固・融解解析で検証
流動限界固相率によるモデル化

図1 流動限界固相率によるモデル化

固液共存状態における液相の流れは固相の体積率（固相率）により影響を受ける。氷を水で溶かす時のように固液界面が平滑な場合には、流体要素の固相率が流動限界固相率以上で圧力と速度を0にして、流体要素は温度と固相率のみ解くことでモデル化することが可能。

図2 解析事例の速度ベクトル

図2は固体の周囲に付着している氷を流水で溶かす解析事例の速度ベクトル。流体要素がわかるように格子を表示し、一定長の速度ベクトルを要素ごとに表示している。

冷凍ブロックの流水解凍解析事例
解析モデル

解析結果
水産物の断面温度

図3 冷凍ブロック

水産物

水産物18尾（各段に9尾）
長さ266 ㎜、幅170 ㎜、高さ50 ㎜

物性値

密度900 kg/㎥、比熱2000 J/(kgK)、
熱伝導率1.40W/(m·K)

温度−20 ℃（冷凍ブロックの初期温度）
4 ℃（流水する水）

流水時間

解析

流動限界固相率

固相率の等値面

秒速50 ㎜で流水した場合

20 分間
非定常解析

0.9

秒速25 ㎜で流水した場合
図4 解析結果 断面温度と固相率の等値面（20分後）
流水秒速50 ㎜（上）、秒速25 ㎜（下）

comment
断面温度を見ますと、手前の1尾以外は0℃以下になっていますので、大部分は半解凍状態と考えられます。
この時点で水産物を流水から取り出し、自然解凍すれば水分のドリップもなく短時間で解凍できることになります。

ソフトウェアクレイドル事例集【混相流編】 ｜15