OpenFoam--UserGuide--记录(二)

2 案例网格细化

2.1 拷贝信息

 

cd $FOAM_RUN  //  run
mkdir cavityFine
cp -r cavity/constant cavityFine
cp -r cavity/system cavityFine

 

2.2 映射算例结果

mapFields 程序把一个给定几何的场信息，映射到另一个几何的对应场。

在例子中，场被认为是相同的，因为被映射场和原始场的几何，边界类型和条件都是一样的。因此在执行mapFields 程序时，可以使用-consistent 附加命令选项。

mapFields 首先读取某个时间步的场信息，这通过被映射场算例的controlDict 中的startFrom/startTime（例如需要被映射的时间步）来指定，在这个算例中，我们打算把粗网格的最终结果映射到cavityFine 的细网格中，由于cavity 算例的最终结果存储在0.5 文件夹中， 因此controlDict 文件中的startTime 应该设定为0.5，startFrom 设定为startTime。现在这个算例就可以运行mapFields 了，键入mapFields -help 可以显示mapFields 的命令参数。可以看出，提供源算例文件的地址是必须的，另外我们还需要使用-consistent 选项， 因此在cavityFine 终端里输入：

mapFields ../cavity -consistent

得到以下信息：

 

 

2.3 控制参数

在算例求解时，最好保持Co<1 。 网格大小减半后时间步长必须减半, 故 controlDict 字典中deltaT需要设定为0.0025。目前场信息设置为每多少个固定时间步写一次。

固定的时间间隔来写结果文件

在controlDict 中的writeControl 中，之前设置为每个固定timeStep 输出，也可以指定为runTime 来指定每多少个物理时间间隔来写结果文件。在这个算例中我们打算每0.1 秒写一次结果。因此我们设定writeControl 为runTime，writeInterval 为0.1。最后，由于算例从一个粗网格的结果来进行计算，因此我们只需要运行短暂的时间就可以达到稳态。我们将endTime 设置为0.7 秒，确认所有的设置正确后我们保存文件。

2.4 绘制数据图

postProcess 中包含的函数功能可以通过键入：

postProcess -list

components 以及mag 函数可用来对矢量进行操作来提取矢量分量以及模的大小

在案例目录下执行命令：

postProcess -func “components(U)”

在Filter 选项中选择Data Analysis 下的Plot Over Line。