dimViscosity

灯具配件

本案例演示应用OpenFOAM创建一个标量输运方程求解器。

本案例要求解的方程为：

莫得瞬态项和源项。

1 文献框架

应用foamNewApp快速创建文献框架。

runfoamNewApp demo10 && cd demo10

扩充终了后的文献框架如下图所示。

图片

本案例相比浅易，只需要一个源文献demo10.C即可。

2 源代码

源文献demo10.C的实质如下。

#include "fvCFD.H" int main(int argc， char *argv[]){#include "setRootCase.H"#include "createTime.H"#include "createMesh.H"     Info << "读取物理场beta" << endl;    volScalarField beta    (        IOobject        (            "beta"，            runTime.timeName()，            mesh，            IOobject::MUST_READ，            IOobject::AUTO_WRITE        )，        mesh    );     Info << "读取物理场U" << endl;    volVectorField U    (        IOobject        (            "U"，            runTime.timeName()，            mesh，            IOobject::MUST_READ，            IOobject::AUTO_WRITE        )，        mesh    );     // 读取字典文献    Info << "读取字典文献transportProperties" << endl;     IOdictionary transportProperties    (        IOobject        (            "transportProperties"，            runTime.constant()， //- 此文献在constant文献夹底下            mesh，            IOobject::MUST_READ_IF_MODIFIED，            IOobject::NO_WRITE        )    );     Info << "从字典文献中读取扩散整个gamma" << endl;    dimensionedScalar gamma("gamma"，
首页-湖名宝香精有限公司 dimViscosity，
首页-和盛佳染料有限公司 transportProperties);     Info << "读取/筹划面通量场phi" << endl;    surfaceScalarField phi    (        IOobject        (            "phi"，固原市尊洲蚕丝有限公司            runTime.timeName()，            mesh，            IOobject::READ_IF_PRESENT， //若phi文献存在，则读取，不然弃取速率与面积相乘筹划取得            IOobject::AUTO_WRITE       // 保存该变量        )，        fvc::interpolate(U) & mesh.Sf() //速率向量与面积向量点积    );     //求解轨则方程    solve(fvm::div(phi， beta) - fvm::laplacian(gamma， beta));     // 创建一个标量场用于存储效果，也不错凯旋用beta    volScalarField result    (        IOobject        (            "result"，灯具配件            runTime.timeName()，            mesh，            IOobject::NO_READ，            IOobject::AUTO_WRITE        )，        beta     //将beta的值拷贝给result    );     result.write();     Info << nl << "筹划终了" << endl;    runTime.printExecutionTime(Info);    return 0;}

编译并测试运行，如下图所示。

图片

3 测试筹划

案例文献需要阐述求解器的需求来准备。上头准备的求解器中需要求解物理量beta。

constant/transportProperties文献中需要指定扩散整个gamma

FoamFile{    version     2.0;    format      ascii;    class       dictionary;    location    "constant";    object      transportProperties;}// * * * * * * * * * * * * * * * * * //gamma            0.01;

由于在读取字典文献时指定了量纲，因此在字典文献中无需叠加指定量纲。

system/fvSolution文献中需要指定beta的求解形状

FoamFile{    version     2.0;    format      ascii;    class       dictionary;    location    "system";    object      fvSolution;}// * * * * * * * * * * * * * * * * * //solvers{    beta    {        solver        GAMG;        smoother    DILUGaussSeidel;        tolerance    1e-6;        relTol        0;    }}

system/fvSchemes文献需要指定轨则方程中对流项、梯度项以及laplacian项的闹翻形状

FoamFile{    version     2.0;    format      ascii;    class       dictionary;    location    "system";    object      fvSchemes;}// * * * * * * * * * * * * * * // ddtSchemes{    default none;}// 梯度项gradSchemes{    default none;    grad(beta) Gauss linear;}//散度项divSchemes{    default none;    div(phi，beta) bounded Gauss upwind;}// laplacian项laplacianSchemes{    default none;    laplacian(gamma，beta) Gauss linear corrected;}// 通量插值interpolationSchemes{    default none;    interpolate(U) linear;} snGradSchemes{    default none;} 

system/controlDict按惯例成就

FoamFile{    version     2.0;    format      ascii;    class       dictionary;    location    "system";    object      controlDict;}// * * * * * * * * * * * * * * * * //startFrom       startTime;startTime       0;stopAt          endTime;endTime         1;deltaT          1;writeControl    timeStep;writeInterval   1;purgeWrite      0;writeFormat     ascii;writePrecision  6;writeCompression off;timeFormat      general;timePrecision   6;runTimeModifiable true;

0/beta文献指定物理场beta的限度条款与运转条款

FoamFile{    version     2.0;    format      ascii;    class       volScalarField;    object      beta;}// 选藏beta是标量dimensions      [0 0 0 0 0 0 0];internalField   uniform 0;boundaryField{    left    {        type            fixedValue;        value           uniform 1;    }     lower    {        type            fixedValue;        value           uniform 0;    }     "(upper|right)"    {        type zeroGradient;    }     frontAndBack    {        type            empty;    }}

0/U文献

FoamFile{    version     2.0;    format      ascii;    class       volVectorField;    object      U;}// * * * * * * * * * * * * * // dimensions      [0 1 -1 0 0 0 0]; internalField   uniform (1 1 0); boundaryField{    "(lower|upper|left|right)"    {        type            fixedValue;        value           $internalField;    }     frontAndBack    {        type            empty;    }}

在案例根旅途下扩充底下的呐喊进行筹划。

blockMeshdemo10

很快就筹划完结。

图片

筹划效果如下图所示。

图片

首页-盛 吉玉棉类有限公司

（终了）

本站仅提供存储就业，通盘实质均由用户发布，如发现存害或侵权实质，请点击举报。

• 灯具配件
• dimViscosity