非保存形式保存形式山口大学大学院坪乡浩一＜特性曲线法.ppt

「移流輸送の高次精度かつ高解像度スキームの開発」 山口大学大学院　理工学研究科 システム工学専攻　社会システム工学大講座 坪郷　浩一 山口大学大学院　　 坪郷　浩一　 ＜移流項計算スキームの検証＞ 水 密度：1000kg/? 粘性係数：1.00×10-3Pa?ｓ 表面張力なし 空気 密度：1.25kg/? 粘性係数：1.50×10-5Pa?ｓ Ｌ＝0.15ｍ，0.6ｍ×0.6m 　計算領域40×40 計算格子間隔Δx=Δy = 0.015m 時間刻みは　Δt = 0.0001sec 山口大学大学院　　 坪郷　浩一　 ＜移流項計算スキームの検証＞ 体積補正法なしで最も数値拡散が少ない組み合わせの検討 山口大学大学院　　 坪郷　浩一　 ＜全質量の相対誤差＞ 山口大学大学院　　 坪郷　浩一　 ＜保存形式高次精度6-point schemeへの拡張 ＞ 保存形式の移流方程式を有限体積法で離散化 離散化 整理すれば はコントロールボリュームの界面濃度だから5点の濃度で以下のように近似できる． まとめると a(=uDt/Dx)はクーラン数である． u　：流速一定 山口大学大学院　　 坪郷　浩一　 a(=uDt/Dx)はクーラン数である． ＜本論文で提案する保存形式高次精度6-point scheme＞ 山口大学大学院　　 坪郷　浩一　 計算格子幅 Δx=200m 計算時間間隔 Δt=100sec 流速 0.5m/sec クーラン数　0.25 中心位置1400m，標準偏差264m，ピーク値10 中心位置2400m，標準偏差264m，ピーク値6.5 の二種類のガウス分布の重ね合わせ 幅2000mの矩形分布 中心位置5800mで最大値10の半楕円分布（空間方向の半径2000m） ＜モデル計算＞ 初期分布 山口大学大学院　　 坪郷　浩一　 ＜universal limiterの導入＞ Δx=200m ， Δｔ= 100sec 特性曲線形式1次精度6-point scheme 保存形式1次精度6-point scheme universal limiter適用 山口大学大学院　　 坪郷　浩一　 ＜universal limiterの導入＞ Δx=200m ， Δｔ= 100sec 特性曲線形式1次精度6-point scheme 保存形式1次精度6-point scheme universal limiter適用 数値振動だけを除去でき，物理的に意味のある極値近傍ではlimiterを機能させない朝位ら(１９９９）のdiscriminatorを導入する． 数値振動は除去されているが，ガウスの極値の再現性の低下が問題 なお，極値の再現性の低下をクリッピングエラーと呼ぶ． 山口大学大学院　　 坪郷　浩一　 ＜朝位らのdiscriminatorの導入＞ 朝位らのdiscriminatorのアルゴリズム 朝位らのdiscriminatorの概念図 山口大学大学院　　 坪郷　浩一　 ＜朝位らのdiscriminatorの導入＞ Δx=200m ， Δｔ= 100sec 保存形式1次精度6-point scheme universal limiter適用 保存形式1次精度6-point scheme universal limiterとdiscriminatorの適用 山口大学大学院　　 坪郷　浩一　 ＜朝位らのdiscriminatorの導入＞ Δx=200m ， Δｔ= 100sec 保存形式1次精度6-point scheme universal limiter適用 保存形式1次精度6-point scheme universal limiterとdiscriminatorの適用 ガウスの極値の再現性は向上したが矩形分布の裾の近傍で数値振動が発生している． 山口大学大学院　　 坪郷　浩一　 ＜朝位らのdiscriminatorの誤作動＞ 朝位らのdiscriminatorが「極値がある」と誤作動したために数値振動が発生した． 改善 「勾配D1とD3が大きく異なる場合」においてlimiterを作用させない． 山口大学大学院　　 坪郷　浩一　 ＜修正された朝位らのdiscriminator＞ 修正された朝位らの discriminatorのアルゴリズム Maxは2個の引数の中で大きい方を選択する演算子である．dは分母が0にならないようにするために導入される非常に小さな値である．eは任意の閾値である． 山口大学大学院　　 坪郷　浩一　 ＜修正された朝位らのdiscriminator＞ 閾値とガウス分布のピーク値の関係 閾値と最小値の関係 数値実験によりeは，3≦e ≦10が推奨される．本論文では， e=3を用いる． 移流問題を