[deepconv/arare/sample]

温位の移流計算

概要

水平 2 次元な面内において温位の移流を計算する. 具体的には, 水平風速 u を与え, さらに温位偏差の初期分布として ガウシアン分布を与え, その時間変化を計算する. ここでチェックすることは,

である.

計算結果

2 次精度中心差分演算ルーチンを用いた結果

温位の移流(1) x 方向の温位の移流. 長い時間ステップと短いタイムステップの クーラン数を同じにする 下流に中心差分を用いたことによるノイズが現われる. 空間分解能を上げるとノイズは小さくなる. 空間分解能を固定しクーラン数を大きくするとノイズは小さくなる.
温位の移流(2a) x 方向の温位の移流. 長い時間ステップと短いタイムステップの クーラン数を変える. タイムフィルターなし. 下流に中心差分を用いたことによるノイズが現われる. 空間分解能を固定し長い時間ステップに対するクーラン数を大きくする とノイズは小さくなる.
温位の移流(2b) x 方向の温位の移流. 長い時間ステップと短いタイムステップの クーラン数を変える. タイムフィルターあり.  
温位の移流(3) 定常な剛体回転流を与えた場合の温位の移流
温位の移流(3β) 定常な剛体回転流を与えた場合の温位の移流 下流に中心差分を用いたことによるノイズが現われる. 初期に与えたトレーサの 振幅は次第に減衰する.
温位の移流(3γ) 定常な剛体回転流を与えた場合の温位の移流, 温位の移流(3β)と同じ設定で数値粘性を入れる プリュームの上昇(4) を安定に計算するのに必要な数値粘性係数よりも小さな値を用いても, 振幅はかなり減少してしまう. 中心差分を用いたことにより下流に生じるノイズを完全に消すことはできない.
温位の移流(4a) xz 方向の温位の移流. 格子点数を固定して, 長い時間ステップの値を変えて計算.
温位の移流(4b) xz 方向の温位の移流. クーラン数を固定して, 格子点数を変えて計算.
温位の移流(5) 定常な剛体回転流を与えた場合の温位の移流 (MPDATA) MPDATA スキームを使用. 下流に差分にともなうノイズは生じないが, 空間分解能が低いと振幅は急激に減衰する.

4 次精度中心差分演算ルーチンを用いた結果

温位の移流(1a) x 方向の温位の移流. 空間解像度を固定して時間ステップを変えて計算する. タイムフィルタなし.  
温位の移流(1b) x 方向の温位の移流. 空間解像度を固定して時間ステップを変えて計算する. タイムフィルタあり.  
温位の移流(1c) x 方向の温位の移流. クーラン数を固定して時間ステップ/空間解像度を変えて計算する. タイムフィルタなし.  
温位の移流(1d) x 方向の温位の移流. クーラン数を固定して時間ステップ/空間解像度を変えて計算する. タイムフィルタあり.  
温位の移流(2a) z 方向の温位の移流. 空間解像度を固定して時間ステップを変えて計算する. タイムフィルタなし.  
温位の移流(2b) z 方向の温位の移流. 空間解像度を固定して時間ステップを変えて計算する. タイムフィルタあり.  
温位の移流(2c) z 方向の温位の移流. クーラン数を固定して時間ステップ/空間解像度を変えて計算する. タイムフィルタなし.  
温位の移流(2d) z 方向の温位の移流. クーラン数を固定して時間ステップ/空間解像度を変えて計算する. タイムフィルタあり.  

失敗例

温位の移流(a) z 方向の温位の移流. 境界条件の影響で境界付近から波が出る
温位の移流(b) x 方向の温位の移流. 短いタイムステップはクーラン数を満たすが, 長い時間ステップはクーラン数を満たさない場合.



Last Update: 2005/04/13 (北守太一)