輻射を考慮した直方体領域内の自然対流（fvDOM モデル）

OpenFOAM 4.x

ケース

$FOAM_TUTORIALS/heatTransfer/buoyantSimpleFoam/hotRadiationRoomFvDOM

概要

直方体領域に置かれた温度 500 K の発熱体（領域 box）の周りの定常流れを計算します。下面（領域 floor）、上面（領域 ceilling）は温度 300 K、側面（領域 fixedWalls）は熱的に解放され熱が逃げられるものとします。ただし解析領域への流体の流入出は無いものとします。

モデル形状 モデル形状

解析では輻射（放射）が考慮されます。輻射の設定はファイル constant/radiationProperties で以下の様にして指定します。輻射のモデルとして finite volume Discrete Ordinary Method（fvDOM）を使用します。

radiation on;
radiationModel  fvDOM;

fvDOMCoeffs
{
    nPhi        3;
    nTheta      5;
    convergence 1e-3;
    maxIter     10;
    cacheDiv    false;
}

solverFreq 10;

absorptionEmissionModel constantAbsorptionEmission;
constantAbsorptionEmissionCoeffs
{
   absorptivity    absorptivity    [ m^-1 ]         0.5;
   emissivity      emissivity      [ m^-1 ]         0.5;
   E               E               [ kg m^-1 s^-3 ] 0;
}

scatterModel    none;
sootModel       none;

また重力加速度の大きさ、方向はファイル constant/g で以下の様にして Z 軸負方向を指定します。

dimensions      [0 1 -2 0 0 0 0];
value           (0 0 -9.81);

メッシュは以下の通りで、メッシュ数は104000です。

メッシュ メッシュ
メッシュ（発熱体付近） メッシュ（発熱体付近）

計算結果は以下の通りです。

温度（T） 温度（T）
XZ 面（Y=1）での流速（U） XZ 面（Y=1）での流速（U）
XZ 面（Y=1）での温度（T） XZ 面（Y=1）での温度（T）
XY 面（Z=0.25）での流速（U） XY 面（Z=0.25）での流速（U）
XY 面（Z=0.25）での温度（T） XY 面（Z=0.25）での温度（T）

計算結果から発熱体部分で上昇気流が起き、自然対流が発生していることがわかります。

実行コマンド

cp -r $FOAM_TUTORIALS/heatTransfer/buoyantSimpleFoam/hotRadiationRoomFvDOM hotRadiationRoomFvDOM
cd hotRadiationRoomFvDOM

blockMesh
buoyantSimpleFoam

paraFoam

計算時間

16分12.51秒　※シングル、Inter(R) Core(TM) i7-2600 CPU @ 3.40GHz 3.40GHz

参照