05、【算例】openFoam盖驱动空腔流动

管网:https://doc.cfd.direct/openfoam/user-guide-v9/cavity

一、算例实现

文件结构

  1. 0:存放初场
  2. constant:存放网格信息
  3. system:存放网格划分、计算等工具

800

1、画网格

blockMesh

800

2、求解

icoFoam

05、【算例】openFoam盖驱动空腔流动

3、保存文件

touch cavity.OpenFOAM

4、后处理

paraview

使用openFoamReader打开
image
显示所有边界
05、【算例】openFoam盖驱动空腔流动
apply
image

二、网格加密

1、网格划分

在icoFoam文件夹内创建cavityFine文件夹,复制cavity文件夹内的constant和system文件夹至该文件夹。
在system/blockMeshDict文件改变网格数量,由之前的(20,20,1)改为(40,40,1),之后使用blockMesh生成网格信息
image.png

2、mapFields映射结果

将粗网格中的0.5结果,映射到细网格中。
修改controlDict文件
800
映射
mapFields ../cavity/cavity -consistent

800

3、计算并输出log文件

// 计算并输出log文件
icoFoam > log
// 查看log文件
cat log

// 可以合并运行
icoFoam > log && cat log

4、calcType

可以对速度的某个量进行可视化处理并咬着计算域内的一条线来绘制一个二维的图

postProcess -func "components(U)"

此时会在时间步长文件夹下出现3个新文件
image

5、后处理

使用paraview加载数据,就会显示刚才计算出的值。
image

三、网格非均匀分布

1、修改blockMeshDict

在icoFoam文件夹内创建asMeshCavity文件夹,复制cavity文件夹内的constant和system文件夹至该文件夹。
修改constant/blockMeshDict文件,如下
image

/*--------------------------------*- C++ -*----------------------------------*\
 =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  9
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    format      ascii;
    class       dictionary;
    object      blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

convertToMeters 0.1;

vertices
(
    (0 0 0)
    (0.5 0 0)
    (1 0 0)
    (0 0.5 0)
    (0.5 0.5 0)
    (1 0.5 0)
    (0 1 0)
    (0.5 1 0)
    (1 1 0)
    (0 0 0.1)
    (0.5 0 0.1)
    (1 0 0.1)
    (0 0.5 0.1)
    (0.5 0.5 0.1)
    (1 0.5 0.1)
    (0 1 0.1)
    (0.5 1 0.1)
    (1 1 0.1)
);

blocks
(
    hex (0 1 4 3 9 10 13 12) (10 10 1) simpleGrading (2 2 1)
    hex (1 2 5 4 10 11 14 13) (10 10 1) simpleGrading (0.5 2 1)
    hex (3 4 7 6 12 13 16 15) (10 10 1) simpleGrading (2 0.5 1)
    hex (4 5 8 7 13 14 17 16) (10 10 1) simpleGrading (0.5 0.5 1)
);

edges
(
);

boundary
(
    movingWall
    {
        type wall;
        faces
        (
            (6 15 16 7)
            (7 16 17 8)
        );
    }
    fixedWalls
    {
        type wall;
        faces
        (
            (3 12 15 6)
            (0 9 12 3)
            (0 1 10 9)
            (1 2 11 10)
            (2 5 14 11)
            (5 8 17 14)
        );
    }
    frontAndBack
    {
        type empty;
        faces
        (
            (0 3 4 1)
            (1 4 5 2)
            (3 6 7 4)
            (4 7 8 5)
            (9 10 13 12)
            (10 11 14 13)
            (12 13 16 15)
            (13 14 17 16)
        );
    }
);

mergePatchPairs
(
);

// ************************************************************************* //

2、画网格

blockMesh

3、生成asyMeshCavity.OpenFOAM文件

touch asyMeshCavity.OpenFOAM

4、后处理

paraview

image.png

5、修改时间步长

image
image

6、映射场

将cavityFine中的0.7文件映射过来

mapFields -consistent ../cavityFine

7、计算输出

icoFoam > log

image

四、大雷诺数

1、创建文件夹并拷贝

创建文件夹

mkdir cavityHighRe

拷贝cavity文件夹至该文件夹

cp -r ../cavity/cavity/ .

2、修改雷诺数

修改constant/transportProperties文件
image
即可提高雷诺数10倍

3、修改控制文件

从0.5开始计算到2s
image

4、计算

nohup nice -19 icoFoam > log && cat log
  1. nohup:党用户退出登录时,程序依然执行;
  2. nice:调整程序优先级,-20对应优先级最高,19对应优先级最低进程。

可以看出1.4s时Ux结束迭代,No Iterations 0表示速度求解停止。
image

五、RAS模型

进入pisoFoam/RAS/cavity
image

可以看到在0文件夹下有许多有关湍流模型的参数
image
在momentumTransport文件里定义模型

/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  9
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    format      ascii;
    class       dictionary;
    location    "constant";
    object      momentumTransport;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

simulationType  RAS;

RAS
{
    model           kEpsilon;

    turbulence      on;

    printCoeffs     on;
}

// ************************************************************************* //

simulationType:laminar, RAS and LES
printCoeffs:为on时,这些参数会被输出到终端

六、CavityClipped

1、划分网格

blockMesh

2、复制0文件夹内容至0.5

cp -r 0 0.5

3、设置

controlDict文件开始时间设置为0.5,结束时间设置为0.6
system/mapFieldsDict文件
有两个参数

patchMap        (lid movingWall);

cuttingPatches  ();
  1. 当用户向把原始场的patches条件投影到背投影场的patches的时候,我们使用patchMap。
  2. 当用户打算把原始场内的数据投影到被投影场的边界时,我们使用cuttingPatches列表

这里lib边界条件根cavity算例的movingWall边界条件是一致的。
image

4、计算

icoFoam > log

5、后处理

image
修改样式
image
image
修改Legend
image
绘制矢量图
image
image

原文链接:https://www.cnblogs.com/dbai/p/17267165.html

本站文章如无特殊说明,均为本站原创,如若转载,请注明出处:05、【算例】openFoam盖驱动空腔流动 - Python技术站

(0)
上一篇 2023年4月18日
下一篇 2023年4月18日

相关文章

  • Qt源码阅读(四) 事件循环

    事件系统 文章为本人理解,如有理解不到位之处,烦请各位指正。 @ 目录 事件系统 什么是事件循环? 事件是如何产生的? sendEvent postEvent 事件是如何处理的? 事件循环是怎么遍历的? 事件过滤器 event 夹带私货时间 Qt的事件循环,应该是所有Qter都避不开的一个点,所以,这篇博客,咱们来了解源码中一些关于Qt中事件循环的部分。先抛…

    C++ 2023年4月18日
    00
  • C++ 入门

    001 c++ 如何工作 任何以 # 开头的语句,都是预处理语句,所谓的预处理语句,在编译之前,就已经被处理了 关键字 include:找到 <> 文件(通常称为“头文件”),然后将 <> 中的所有内容拷贝到现在的文件里 main()比较特殊,虽然它的返回值类型是 int,但它不一定需要返回值,如果不设置返回值,默认返回 0 <…

    C++ 2023年5月11日
    00
  • C++重载的奥义之运算符重载

    0、引言         重载,顾名思义从字面上理解就是重复装载,打一个不恰当的比方,你可以用一个篮子装蔬菜,也可以装水果或者其它,使用的是同一个篮子,但是可以用篮子重复装载的东西不一样。         正如在之前的文章《重载的奥义之函数重载》中介绍的类似,函数的重载是指利用相同的函数名设计一系列功能相近,但是功能细节不一样的函数接口;因此运算符重载也是指…

    C++ 2023年4月18日
    00
  • 【Visual Leak Detector】源码下载

    说明 使用 VLD 内存泄漏检测工具辅助开发时整理的学习笔记。本篇介绍 VLD 源码的下载。同系列文章目录可见 《内存泄漏检测工具》目录 目录 说明 1. 下载途径 2. 不同下载途径的源文件差异 1. 下载途径 以 v2.5.1 版本为例,可以到 Github-KindDragon-vld 页面下载 master 的 zip 源码包,如下所示: 也可以到 …

    C++ 2023年4月22日
    00
  • 【Visual Leak Detector】源码文件概览

    说明 使用 VLD 内存泄漏检测工具辅助开发时整理的学习笔记。本篇对 VLD 源码包中的各文件用途做个概述。同系列文章目录可见 《内存泄漏检测工具》目录 目录 说明 1. 整体概览 2. 文件夹 .teamcity 3 文件夹 lib 3.1 文件夹 cppformat(生成 libformat) 3.2 文件夹 dbghelp 3.3 文件夹 gtest(…

    C++ 2023年4月24日
    00
  • STL容器之queue

    是什么 循环队列, FIFO先进先出 怎么用 初始化 //C11 deque<int> deq{1,2,3,4,5}; //拷贝构造,可以拷贝deque queue<int> que(deq); //100个5 queue<int> que2(100,5); //运算符重载 que2 = que; 操作 //队尾添加元素 …

    C++ 2023年4月17日
    00
  • 驱动开发:通过MDL映射实现多次通信

    在前几篇文章中LyShark通过多种方式实现了驱动程序与应用层之间的通信,这其中就包括了通过运用SystemBuf缓冲区通信,运用ReadFile读写通信,运用PIPE管道通信,以及运用ASYNC反向通信,这些通信方式在应对一收一发模式的时候效率极高,但往往我们需要实现一次性吐出多种数据,例如ARK工具中当我们枚举内核模块时,往往应用层例程中可以返回几条甚至…

    C++ 2023年4月30日
    00
  • 【Visual Leak Detector】VS 中 VLD 输出解析

    说明 使用 VLD 内存泄漏检测工具辅助开发时整理的学习笔记。同系列文章目录可见 《内存泄漏检测工具》目录 目录 说明 1. 使用方式 2. 输出报告 1. 使用方式 在 VS 中使用 VLD 的方法可以查看另外一篇博客:在 VS 2015 中使用 VLD。 2. 输出报告 在 VS 中使用 VLD 时的输出报告,与在 QT 中使用时是一致的,输出内容的解析…

    C++ 2023年4月17日
    00
合作推广
合作推广
分享本页
返回顶部