OpenFoam 建筑物绕流计算设置

最新推荐文章于 2024-05-17 13:26:14 发布
最新推荐文章于 2024-05-17 13:26:14 发布
阅读量88 收藏 1
点赞数
文章标签： c++
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_22015481/article/details/133885075
版权

计算建筑物绕流模拟时，需要的设定文件和内容。

一，polyMesh 中 boundary文件内容。

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  5.x                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
/*   Windows 32 and 64 bit porting by blueCAPE: http://www.bluecape.com.pt   *\
|  Based on Windows porting (2.0.x v4) by Symscape: http://www.symscape.com   |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       polyBoundaryMesh;
    location    "constant/polyMesh";
    object      boundary;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

8
(
    e
    {
        type            patch;
        nFaces          808;
        startFace       3707465;
    }
    n
    {
        type            patch;
        nFaces          4756;
        startFace       3708273;
    }
    s
    {
        type            patch;
        nFaces          808;
        startFace       3713029;
    }
    w
    {
        type            patch;
        nFaces          808;
        startFace       3713837;
    }
    top
    {
        type            symmetry;
        inGroups        1(symmetry);
        nFaces          16;
        startFace       3714645;
    }
    bottom
    {
        type            wall;
        inGroups        1(wall);
        nFaces          43390;
        startFace       3714661;
    }

    inner_bottom
    {
        type            wall;
        physicalType    patch;
        inGroups        1(wall);
        nFaces          909;
        startFace       3766612;
    }

    bldg
    {
        type            wall;
        physicalType    patch;
        inGroups        1(wall);
        nFaces          8561;
        startFace       3758051;
    }
)

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

二，constant / transportProperties 和 turbulenceProperties 文件内容。

transportProperties

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  4.x                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      transportProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

transportModel  Newtonian;

nu              [0 2 -1 0 0 0 0] 1.5e-05;

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

 turbulenceProperties

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  4.x                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      RASProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

simulationType RAS;

RAS
{
    RASModel        kEpsilon;

    turbulence      on;

    printCoeffs     on;
}

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

三，初始变量 0/U，p, k, epsilon, nut 文件内容。

0/U

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  4.x                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volVectorField;
    object      U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (0 0 0);

boundaryField
{
    e
    {
        type            zeroGradient;
    }

    n
    {
	    type            atmBoundaryLayerInletVelocity;
        flowDir         (0 -1 0);
        zDir            (0 0 1);
        Uref            9.52;
        Zref            6.0;
        z0              uniform 0.01;
        zGround         uniform 0.0;
	    value         uniform (0 0 0);
	    initABL        true;
    }

    s
    {
        type            zeroGradient;
    }

    w
    {
        type            zeroGradient;
    } 
   
    top
    {
        type            symmetry;
    }

    bottom
    {
        type            noSlip;
    }

    inner_bottom
    {
        type            noSlip;
    }

    bldg
    {
        type            noSlip;
    }
}

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

 0/p

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  4.x                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 2 -2 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    e
    {
        type            zeroGradient;

    }
    n
    {
        type            zeroGradient;

    }

    s
    {
        type            fixedValue;
        value           uniform 0;
    }

    w
    {
        type            zeroGradient;

    }
       
    top
    {
        type            symmetry;

    }
    
    bottom
    {
        type            zeroGradient;

    }    
    
    inner_bottom
    {
        type            zeroGradient;
    }
    bldg
    {
        type            zeroGradient;
    }
}

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

0/k 

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  4.x                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      k;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

kInlet          1.5;   // approx k = 0.5*(I*U)^2 ; I = 0.1

dimensions      [0 2 -2 0 0 0 0];

internalField   uniform $kInlet;

boundaryField
{
    e
    {
        type            zeroGradient;
    }

    n
    {
	    type            atmBoundaryLayerInletK;
        flowDir         (0 -1 0);
        zDir            (0 0 1);
        Uref            9.52;
        Zref            6.0;
        z0              uniform 0.01;
        zGround         uniform 0.0;
	    value         uniform 0;
	    initABL        true;
    }

    s
    {
        type            zeroGradient;
    }

    w
    {
        type            zeroGradient;
    }
    
    top
    {
        type            symmetry;
    }

   bottom
    {
        type            kqRWallFunction;
        value           uniform $kInlet;
    }   

    inner_bottom
    {
        type            kqRWallFunction;
        value           uniform $kInlet;
    }

    bldg
    {
        type            kqRWallFunction;
        value           uniform $kInlet;
    }
}

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

0/epsilon 

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  4.x                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      epsilon;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

epsilonInlet  0.03; // Cmu^0.75 * k^1.5 / L ; L =10

dimensions      [0 2 -3 0 0 0 0];

internalField   uniform $epsilonInlet;

boundaryField
{
    e
    {
        type            zeroGradient;
    }

    n
    {
	    type            atmBoundaryLayerInletEpsilon;
        flowDir         (0 -1 0);
        zDir            (0 0 1);
        Uref            9.52;
        Zref            6.0;
        z0              uniform 0.01;
        zGround         uniform 0.0;
	    value         uniform 0;
	    initABL        true;
    }

    s
    {
        type            zeroGradient;
    }

    w
    {
        type            zeroGradient;
    }
    
    top
    {
        type            symmetry;
    }
    
    bottom
    {
        type            epsilonWallFunction;
        value           uniform $epsilonInlet;
	    //z0               uniform 0.01;
    }

    inner_bottom
    {
        type            epsilonWallFunction;
        value           uniform $epsilonInlet;
	    //z0               uniform 0.01;
    }

    bldg
    {
        type            epsilonWallFunction;
        value           uniform $epsilonInlet;
	    //z0               uniform 0.005;
    }

    
}

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

0/nut 

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  4.x                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      nut;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 2 -1 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    e
    {
        type            calculated;
        value           uniform 0;

    }

    n
    {
        type            calculated;
        value           uniform 0;
    }

    s
    {
        type            calculated;
        value           uniform 0;
    }

    w
    {
        type            calculated;
        value           uniform 0;
    }

    top
    {
        type            symmetry;
    }

   bottom
    {
        type            nutkWallFunction;
        value           uniform 0;
	   //z0             uniform 0.01
    }
    
    inner_bottom
    {
        type            nutkWallFunction;
        value           uniform 0;
	    //z0            uniform 0.01
    }

    bldg
    {
        type            nutkWallFunction;
        value           uniform 0;
	    //z0            uniform 0.005
    }
}

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

四，system 下文件设置

system/controlDict

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v1706+                                |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
/* Butterfly 0.0.4                https://github.com/ladybug-tools/butterfly *\
\*---------------------------------------------------------------------------*/
FoamFile
{
	version		4.0;
	format		ascii;
	class		dictionary;
	location	"system";
	object		controlDict;
}

//#include		"probes";

application		simpleFoam;

startFrom		latestTime;

startTime		0;

stopAt		endTime;

endTime		3000;

deltaT		1;

writeControl		timeStep;

writeInterval		200;

purgeWrite		0;

writeFormat		ascii;

writePrecision		8;

writeCompression		off;

timeFormat		general;

timePrecision		6;

runTimeModifiable		true;

functions{

 #includeFunc  probes
 #includeFunc  probes_lines_0.1m
 #includeFunc  probes_lines_0.01m
 #includeFunc  probes_lines_0m
 #includeFunc  residuals

}

system/fvSchemes

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v1706+                                |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
/* Butterfly 0.0.4                https://github.com/ladybug-tools/butterfly *\
\*---------------------------------------------------------------------------*/
FoamFile
{
	version		4.0;
	format		ascii;
	class		dictionary;
	location	"system";
	object		fvSchemes;
}

ddtSchemes
{

    default		steadyState;

}

gradSchemes
{

    default		cellLimited leastSquares 1;

}

divSchemes
{

    default		none;

    div(phi,epsilon)		bounded Gauss linearUpwind grad(epsilon);

    div(phi,U)		bounded Gauss linearUpwindV grad(U);

    div((nuEff*dev2(T(grad(U)))))		Gauss linear;

    div(phi,k)		bounded Gauss linearUpwind grad(k);

}

laplacianSchemes
{

    default		Gauss linear limited corrected 0.333;

}

interpolationSchemes
{

    default		linear;

}

snGradSchemes
{

    default		limited corrected 0.333;

}

fluxRequired
{

    default		no;

}

system/fvSolution

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v1706+                                |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
/* Butterfly 0.0.4                https://github.com/ladybug-tools/butterfly *\
\*---------------------------------------------------------------------------*/
FoamFile
{
	version		4.0;
	format		ascii;
	class		dictionary;
	location	"system";
	object		fvSolution;
}

solvers
{

    p
    {

        agglomerator		faceAreaPair;

        relTol		0.1;

        tolerance		1e-7;

        nCellsInCoarsestLevel		10;

        smoother		GaussSeidel;

        solver		GAMG;

        cacheAgglomeration		on;

        nPostSweeps		2;

        nPreSweepsre		0;

        mergeLevels		1;

    }

    U
    {

        relTol		0.1;

        tolerance		1e-8;

        nSweeps		1;

        smoother		GaussSeidel;

        solver		smoothSolver;

    }

    k
    {

        relTol		0.1;

        tolerance		1e-8;

        nSweeps		1;

        smoother		GaussSeidel;

        solver		smoothSolver;

    }

    epsilon
    {

        relTol		0.1;

        tolerance		1e-8;

        nSweeps		1;

        smoother		GaussSeidel;

        solver		smoothSolver;

    }

}

SIMPLE
{

    nNonOrthogonalCorrectors		2;

    residualControl
    {

        nut		1e-8;

        k		1e-8;

        U		1e-8;

        p		1e-8;

        epsilon	1e-8;

    }

}

relaxationFactors
{

    k		0.7;

    U		0.7;

    epsilon		0.7;

    p		0.3;

}

不知魏晋CFD
关注
  • 0
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
OpenFOAM-圆柱绕流
asd165654的博客
08-12 1998
原版视频下载地址：https://yunpan.cn/c64yrdt9J5LmQ 访问密码 0128 首先进行建模操作，任何建模软件均可，本教程采用ICEM直接建模，模型尺寸如下： 建成的模型如下： 对建好的模型进行网格划分，划分完成的网格如下： 对圆柱的近壁面进行了加密处理，将划分的网格导出为ASCII的.msh格式（注：二进制的.msh格式OpenFO...
python做圆柱绕流_OpenFOAM-圆柱绕流
weixin_39721953的博客
11-30 1205
首先进行建模操作，任何建模软件均可，本教程采用ICEM直接建模，模型尺寸如下：建成的模型如下：对建好的模型进行网格划分，划分完成的网格如下：对圆柱的近壁面进行了加密处理，将划分的网格导出为ASCII的.msh格式（注：二进制的.msh格式OpenFOAM是不支持网格转换的）接下来转入OpenFOAM的操作：首先新建一个文件夹，名字任取，本教程中我将该文件夹命名为：train然后进入OpenFOAM...
OpenFOAM的圆柱绕流算例记录
热门推荐
陈十七的博客
09-08 1万+
圆柱绕流是流体力学里最经典的算例之一，受到众多学者的“喜爱”。初学CFD的人除了空腔、管道流之外，最先接触的莫过于圆柱绕流了。圆柱绕流能反映流体流过障碍物时的各种情况：层流向湍流的转捩、流体的剪切与分离、边界层、卡门涡街现象等等，是深入理解流动机理的重要算例。因为经常要算圆柱绕流，在这里做一些记录，记录一些细节，以防以后忘掉。这次做Re=200时的圆柱绕流算例，画网格用ICEM，CFD软件用开源软件OpenFOAM，后处理用tecplot。 1. 前处理：画网格 用ICEM画网格。圆柱直径D=1m，计算
OpenFoam 并行计算设置
qq_22015481的博客
10-17 349
初始值设定时，如果使用了入口大气边界条件，应该设定 value 值。一， 准备 decomposeParDict 文件。OpenFoam 并行计算设置需要的步骤。
openfoam大规模数据计算效率提升的分析
lijianhua1205的专栏,微信号 lijianhua1205
11-24 530
openfoam大规模数据计算效率提升的分析 openfoam大规模数据计算的效率提升，主要是两个方面 1 大规模数据文件的读取效率 2 大规模线性方程组计算的效率 提升效率主要有并行化和分布式两个手法，线性方程组并行化计算的开源数学库有一些，但是和openfoam并不是简单的集成就一定能提高效率 大规模线性方程组分布式计算的研究要少些，但是基于云服务器的分布式计算平台是非常廉价的（和高性能计算服务器或超算比较），也是一个很有价值的技术方向 ...
OpenFOAM——圆柱绕流对流换热
asd165654的博客
01-15 2130
本算例来自《ANSYS FLUENT技术基础与工程应用：流动传热与环境污染控制领域》 TOP和DOWN为对称边界（symmetry），入口速度为0.01m/s，入口温度为300K，圆柱温度为350K 流体的物性参数： 密度：1.225kg/m3 导热系数：0.0242W/(m·K) 动力黏度：1.7894×10-5kg/(m·s) 运动黏度：1.4607...
OpenFoam 计算过程中数据监测设置
qq_22015481的博客
10-17 754
三，计算开始后，postProcessing 文件夹 出现 residuals 和 probes 文件夹记录监测数据。使用gnu plot 或者 python 处理数据，进行可视化。二，在system/controlDict 中添加功能项，链接 residuals 和 probes。使用OpenFoam计算时，过程中的残差，自定义点位置变量值的监测，需要进行的设置。一，在system下创建文件。
OpenFOAM提取等值面并计算面积
李文韬的博客
08-02 934
OpenFOAM的isoSurface类简介
OpenFOAM分段边界条件的设置
CloudBird07的博客
10-23 3869
问题描述 在实际使用OpenFOAM中，我们经常会遇到根据平面上点的坐标给出边界条件的情况，例如一个三维方腔，在距离墙壁底端3m的位置上有一个直径0.5m的管道入口，流体从其中以恒定的流量流入计算域。在这样的背景下，我们就需要对这个壁面给定分段的边界条件。在本文中我们就具体讨论一下如何设置这样的边界条件。 具体步骤 我们用到的两个工具是creatPatch和topoSet。 creatPatch需要在算例的system文件下调用creatPatchDict文件，其原本用来调整周期性边界条件的位置，在这里我
OpenFOAM低雷诺数圆柱绕流算例
09-08
用OpenFOAM进行低雷诺数圆柱绕流计算，里面是一整个算例文件
cylinder_openfoam_cylinder_圆柱绕流涡量_transportcylinder_
10-01
openfoam计算二维圆柱绕流的代码
实践作业之OpenFOAM实现cavity算例温度场计算.rar
06-28
文件实现在cavity算例中增加温度场，在模型左右边界分别施加固定温度，求解模型温度场分布
openfoam-v2012
12-25
openfoam最新版，供流体力学科研人员下载使用
openFoam 3D方柱绕流仿真说明文档+源代码
09-26
openFoam是一款开源CFD平台，本资源适合初学OF的CFD研究人员
初识C语言——第二十一天
硬件菜鸟的成长之路
05-15 620
1.rand()函数是只能生成伪随机数，并且无返回值，使用前需调用srand()伪随机数生成器作为起点2.time（）函数是结构体函数，返回值是64为整数类型，可以通过time函数生成时间戳时间戳-函数调用的时间点和计算机的初始时间点之间的时间间隔3.void（）函数-无返回值类型函数4.while(1)-一直循环，直到break跳出。
C语言.动态内存管理
Joseit的博客
05-16 590
栈区（stack）：在执行函数时，函数内局部变量的存储单元都可以在栈上创建，函数执行结束时这些存储单元自动被释放。栈内存分配运算内置于处理器的指令集中，效率很高，但是分配的内存容量有限。栈区主要存放运行函数而分配的局部变量、函数参数、返回数据、返回地址等。堆区（heap）：一般由程序员分配释放， 若程序员不释放，程序结束时可能由OS(操作系统)回收。分配方式类似于链表。数据段（静态区）：（static）存放全局变量、静态数据。程序结束后由系统释放。
【CCF-CSP】202312-4 宝藏
whaoe_m的博客
05-14 408
对于所有类型为 2的事件，输出一行四个非负整数 C1,1​,C1,2​,C2,1​,C2,2​，表示该时刻的密码 C。西西艾弗岛上埋藏着一份宝藏，小 C 根据藏宝图找到了宝藏的位置。由于小 C 并不会这个问题，他向你发起了求助。你需要帮助小 C 求出所有类型为 22 的事件所对应的密码。小 C 将所有的时刻发生的事件均记录了下来。输入的第一行包含两个正整数n,m。从标准输入读入数据。
C++ 11
最新发布
qq_65112813的博客
05-17 658
对于类类型，如果定义了接受。
openfoam圆柱绕流
05-13
OpenFOAM是一款自由、开源的计算流体力学（CFD）软件，可以实现包括圆柱绕流在内的各种数值模拟。 圆柱绕流是CFD中的一个经典问题，也称为柱面绕流。该问题指的是，将一个圆柱体置于流体中，当流体流过圆柱时，会形成一定的涡流和压力分布，经常用来做为CFD算例测试。 对于圆柱绕流，在OpenFOAM中可以使用不同的求解器和网格进行计算。其中，最常用的求解器是ICOFoam，它是OpenFOAM自带的求解器之一，可以求解不可压缩流体的Navier-Stokes方程。 在网格选择方面，可以使用结构网格或非结构网格。 关于模拟结果，可以显示出在圆柱表面靠近尾部会出现小型的环流，越往靠近半径中心处环流越大，而在流体的底部和顶部也会出现涡流。除此之外，还会显示出压力分布的规律，随着圆柱与流体的相对速度的变化，压力分布也有所不同。 总之，OpenFOAM可以通过不同的求解器和网格来模拟圆柱绕流现象，并给出与实际相符的模拟结果。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值