算例路径： olaFlow\tutorials\baseWaveFlume
算例描述： 一个基础的二维波浪水槽
算例快照：

文件结构：

├── 0.org
│   ├── U 
│   ├── alpha.water
│   ├── alpha.water.org
│   └── p_rgh
├── cleanCase
├── constant
│   ├── dynamicMeshDict
│   ├── g
│   ├── transportProperties
│   ├── turbulenceProperties
│   └── waveDict   --> 设置波浪要素
├── runCase
└── system├── blockMeshDict├── controlDict├── decomposeParDict├── fvSchemes├── fvSolution└── setFieldsDict  --> 设置水深

算例文件解析：

【0.org\U】

dimensions      [0 1 -1 0 0 0 0];  // 量纲 m/s
internalField   uniform (0 0 0);   // 内部速度场 均一场 0 0 0
boundaryField   //边界场
{inlet      // 造波边界{type            waveVelocity; // 波浪速度waveDictName    waveDict;     // 读取constant\waveDict中的波浪要素value           uniform (0 0 0);  //初值为 0 0 0}outlet    // 消波边界{type            waveAbsorption2DVelocity;  // 使用了二维消波理论，olaFlow采用主动消波法value           uniform (0 0 0);}bottom   // 底部边界为固壁边界，边界上速度为零{type            fixedValue;   value           uniform (0 0 0);}atmosphere // 大气边界，允许空气流出和流入{type            pressureInletOutletVelocity;value           uniform (0 0 0);}frontAndBack // 前后面，empty指示模型为二维模型{type            empty;}
}

【0.org\p_rgh】

// p_rgh = p - rgh，实际压力减去静水压力
dimensions      [1 -1 -2 0 0 0 0];  // M(1) L(-1) T(-2)
internalField   uniform 0;
boundaryField
{frontAndBack{type            empty;}outlet{type            fixedFluxPressure; //将压力梯度设置为0，边界上的通量由速度边界条件指定value           uniform 0;}inlet {type            fixedFluxPressure;value           uniform 0;}bottom{type            fixedFluxPressure;value           uniform 0;}atmosphere {type            totalPressure; //总压条件:流出 p = p0; 流入 p = p0 - 0.5|U|^2U               U;phi             phi;rho             rho;psi             none;gamma           1;p0              uniform 0;value           uniform 0;}
}

【0.org\alpha.water.org】

// 设置流体体积分数
dimensions      [0 0 0 0 0 0 0]; 
internalField   uniform 0;       
boundaryField
{inlet{type            waveAlpha;  // 根据波浪条件设置waveDictName    waveDict;value           uniform 0;}frontAndBack{type            empty;}outlet{type            zeroGradient; // 零梯度}bottom{type            zeroGradient;}atmosphere{type            inletOutlet; // 当流体流出时，α的梯度为0；如果流入，α为0，100%的空气流入inletValue      uniform 0;value           uniform 0;}
}

【constant\dynamicMeshDict】

dynamicFvMesh        staticFvMesh;

【constant\g】

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

【constant\transportProperties】

phases (water air);water
{transportModel  Newtonian;nu              [0 2 -1 0 0 0 0] 1e-06;  // 流体运动粘度rho             [1 -3 0 0 0 0 0] 1000;   // 流体密度
}air
{transportModel  Newtonian;nu              [0 2 -1 0 0 0 0] 1.48e-05;rho             [1 -3 0 0 0 0 0] 1;
}sigma           [1 0 -2 0 0 0 0] 0.07;     // 水和空气之间的表面张力参数

【constant\turbulenceProperties】

simulationType  laminar;     // 设置为层流模型

【constant\waveDict】

waveType        regular;    // 规则波
waveTheory      cnoidal;    // 椭圆余弦波
genAbs          1;          // 考虑造波边界的消波性能  1/0
absDir          0.0;        // 造波边界的消波方向
nPaddles        1;          // 主动消波的Paddles数量设置
waveHeight      0.10;       // 波高
wavePeriod      3;          // 波周期
waveDir         0.0;        // 波向
wavePhase       1.57079633; // 初始相位
// Change both entries to true to re-read this dictionary upon restart.
rereadAlpha     false;
rereadU         false;

【system\blockMeshDict】

scale 1;
vertices        
((0.0     -0.02   0.0)(10.0    -0.02   0.0)(10.0    -0.02   0.7)(0.0     -0.02   0.7)(0.0     0.0     0.0)(10.0    0.0     0.0)(10.0    0.0     0.7)(0.0     0.0     0.7)
);
blocks          
(hex (0 1 5 4 3 2 6 7)   (500 1 70) simpleGrading (1 1 1)
);
edges           
(
);
patches         
(patch inlet  // 造波边界((0 4 7 3))patch outlet  // 消波边界 ((1 5 6 2))wall bottom   // 水槽底部边界((0 1 5 4))patch atmosphere  // 大气边界((3 2 6 7))empty frontAndBack // 水槽侧面边界((0 1 2 3)(4 5 6 7))
);
mergePatchPairs
(
);

【system\controlDict】

application     olaFlow;      // olaFlow求解器
startFrom       latestTime;
startTime       0;
stopAt          endTime;
endTime         60;
deltaT          0.001;       // 计算时间步
writeControl    adjustableRunTime;
writeInterval   0.05;        // 写出时间步
purgeWrite      0;
writeFormat     ascii;
writePrecision  6;
compression     off;         // 是否压缩格式写出，可节约硬盘空间, on/off
timeFormat      general;
timePrecision   6;
runTimeModifiable yes;
adjustTimeStep  yes;         // 采用自适应时间步，可能会加速计算，也可能造成时间步极小
maxCo           0.5;         // CFL条件的Courant数, 一般<1, 设置一个小值会使计算结果更精确，但也减小了时间步长，增加了计算成本
maxAlphaCo	    0.5;         // 两相交界面上的最大Courant数
maxDeltaT       0.025;

【system\decomposeParDict】

numberOfSubdomains 2;        // 并行区域分解数目
method          scotch;      // 区域分解方法
...

【system\fvSchemes】

// 指定控制方程中各项的有限体积法的离散格式
ddtSchemes  // 指定时间离散格式
{default         Euler;  // Euler法，一阶精度，条件稳定
} gradSchemes  // 梯度项离散格式
{default         Gauss linear; // 高斯定理，将网格中心的量插值到网格面上
}// olaFlow 的算例中给出了几乎所有可能出现的项，具体算例可能不会包含全部项
divSchemes  // 对流项与散度项的离散格式， 将网格中心的量插值到网格面上，因此实际上选用的是interpolationSchemes
{div(rhoPhi,U)  Gauss limitedLinearV 1;  // Guass limitedLinear(一种TVD格式，使同时满足精度和有界) V类(采用限制器时考虑了流动方向)div(U)  Gauss linear; // 二阶精度，无界div((rhoPhi|interpolate(porosity)),U)  Gauss limitedLinearV 1;div(rhoPhiPor,UPor)  Gauss limitedLinearV 1;div(rhoPhi,UPor)  Gauss limitedLinearV 1;div(rhoPhiPor,U)   Gauss limitedLinearV 1; div(phi,alpha)  Gauss vanLeer; // Gauss vanLeer(一种TVD格式，使同时满足精度和有界)div(phirb,alpha) Gauss interfaceCompression; // 界面压缩格式，基于一般限制格式div((muEff*dev(T(grad(U))))) Gauss linear;div(phi,k)      Gauss upwind; // 一阶迎风格式，有界div(phi,epsilon) Gauss upwind;div((phi|interpolate(porosity)),k)      Gauss upwind;div((phi|interpolate(porosity)),epsilon) Gauss upwind;div(phi,omega) Gauss upwind;div((phi|interpolate(porosity)),omega) Gauss upwind;
}laplacianSchemes  // 拉普拉斯项离散格式
{default         Gauss linear corrected; // Guass线性插值，corrected(显式的非正交网格修正)
}interpolationSchemes
{default         linear;  // 线性插值格式
}snGradSchemes  // 面法向梯度格式
{default         corrected;
}fluxRequired  // 
{default         no;p_rgh;pcorr;alpha.water;
}

【system\fvSolution】

// 指定方程组矩阵求解器、残差以及其他算法控制
solvers
{"alpha.water.*"{nAlphaCorr      1;  //nAlphaSubCycles 2;alphaOuterCorrectors yes;cAlpha          1;MULESCorr       no;nLimiterIter    3;solver          smoothSolver;    // 求解器：光顺求解器。对称和非对称矩阵均适用smoother        symGaussSeidel;  // 光顺器：对称Gauss-Seidel方法tolerance       1e-8;            // 标准化残差relTol          0;               // 相对残差，表征残差的递减量}"pcorr.*"  // 压强校正量{solver          PCG;            // 求解器：预处理共轭梯度法。仅适用于对称矩阵    preconditioner  DIC;            // 预处理器：基于对角的不完全Cholesky预处理器tolerance       1e-5;relTol          0;}p_rgh{solver          PCG;preconditioner  DIC;tolerance       1e-07;relTol          0.05;}p_rghFinal  // 在求解p_rgh时可设置较大的残差，在最后一步设置严格的误差{$p_rgh;relTol          0;}U{solver          smoothSolver; smoother        symGaussSeidel;tolerance       1e-06;relTol          0;}
}// 离散方程组的分离式解法
PIMPLE   // SIMPLE 与 PISO 算法的耦合
{momentumPredictor   no;nOuterCorrectors    1;  // 值不大于1时，算法变为纯PISO算法；若大于1，为非定常SIMPLE算法nCorrectors         3;nNonOrthogonalCorrectors 0;
}relaxationFactors // 松弛方法，为降低数值波动
{fields        // 表示需要使用松弛方法(显式)的变量{}equations     // 表示需要使用松弛方法(隐式)的方程{".*" 1;}
}

【system\setFieldsDict】

defaultFieldValues
(volScalarFieldValue alpha.water 0   // 初始化流体体积分数
);
regions
(boxToCell{box (-10 -1 -1) (30 1 0.4);fieldValues(volScalarFieldValue alpha.water 1  // 将box范围内的流体设置为a = 1的相，即水深);}
);