【OpenFOAM】-olaFlow-算例5- oppositeSolitariesFlume

news2025/1/11 23:48:14

算例路径： olaFlow\tutorials\oppositeSolitariesFlume
算例描述： 两列反向的孤立波相互作用
学习目标： 熟练掌握olaFlow的造波设置，波浪方向与消波方向设置
算例快照：
在这里插入图片描述

图1 两列反向孤立波相互作用

文件结构：

.
├── 0.org
│   ├── U
│   ├── alpha.water
│   ├── alpha.water.org
│   └── p_rgh
├── cleanCase
├── constant
│   ├── dynamicMeshDict
│   ├── g
│   ├── transportProperties
│   ├── turbulenceProperties
│   ├── waveDict
│   └── waveDictOpposite
├── runCase
└── system
    ├── blockMeshDict
    ├── controlDict
    ├── decomposeParDict
    ├── fvSchemes
    ├── fvSolution
    └── setFieldsDict

算例文件解析：

【0.org\U】

dimensions      [0 1 -1 0 0 0 0];
internalField   uniform (0 0 0);
boundaryField
{
    inlet // 造波边界 1
    {
        type            waveVelocity;
        waveDictName    waveDict;    // 读取constant/waveDict文件
        value           uniform (0 0 0);
    }
    outlet // 造波边界 2
    {
        type            waveVelocity;
        waveDictName    waveDictOpposite;    // 读取constant/waveDictOpposite文件
        value           uniform (0 0 0);
    }
    bottom
   ...
    atmosphere
   ...
    frontAndBack
   ...
}

【0.org\alpha.water】

dimensions      [0 0 0 0 0 0 0];
internalField   uniform 0;
boundaryField
{
    inlet
    {
        type            waveAlpha;
        waveDictName    waveDict;    // 读取constant/waveDict文件
        value           uniform 0;
    }
    outlet
    {
        type            waveAlpha;
        waveDictName    waveDictOpposite;    // 读取constant/waveDictOpposite文件
        value           uniform 0;
    }
    bottom
   ...
    atmosphere
   ...
    frontAndBack
   ...
}

【0.org\p_rgh】
参考【OpenFOAM】-olaFlow-算例1- baseWaveFlume

【constant\dynamicMeshDict】
【constant\g】
【constant\transportProperties】
【constant\turbulenceProperties】
以上均与一般水槽设置相同，参考【OpenFOAM】-olaFlow-算例1- baseWaveFlume

【constant\waveDict】

waveType        solitary; // 孤立波
waveTheory      Grimshaw; // 孤立波理论
genAbs          1;        // 消波
absDir          0.0;      // 消波方向。 修正速度设置方向(角度为消波边界指向域内的法向量与x轴正方向的夹角)。与waveDir相同，单位为度。方向须指向域内，否则不起作用或适得其反。若大于360°(默认值)，则边界条件选择与造波板垂直的方向。
nPaddles        1;        // 
waveHeight      0.10;
waveDir         0.0;      // 波浪方向。X轴正方向为零度，角度按逆时针(转向Y轴)计。
// Change both entries to true to re-read this dictionary upon restart.
rereadAlpha     false;
rereadU         false;

【constant\waveDictOpposite】

waveType        solitary;
waveTheory      Grimshaw;
genAbs          1;
absDir          180.0;     // 消波方向
nPaddles        1;
waveHeight      0.10;
waveDir         180.0;     // 波浪方向
// Change both entries to true to re-read this dictionary upon restart.
rereadAlpha     false;
rereadU         false;

【system\blockMeshDict】
【system\controlDict】
【system\decomposeParDict】
【system\fvSchemes】
【system\fvSolution】
【system\setFieldsDict】
参考【OpenFOAM】-olaFlow-算例1- baseWaveFlume