诈欺底下的号召创建案例文献结构
发布日期：2024-06-21 15:04 点击次数：74
Open FOAM是Open Source Field Operation and Manipulation的英文缩写，从其名字就不错看出，其是一个用来操纵和处置场数据的开源轮番库。CFD中也包含大皆的物理场，如压力场、速率场、组分场等，这些物理场有的是标量场，有的是向量场。诈欺OpenFOAM很容易收场对此类物理场的操作。
本案例演示诈欺OpenFOAM操作场数据。1 准备文献本演示案例只需要一个源文献。诈欺底下的号召创建案例文献结构。
cd $FOAM_RUNmkdir demo4 && cd demo4mkdir Maketouch demo4.C Make/files Make/optionscp -r $FOAM_TUTORIALS/incompressible/icoFoam/cavity/cavity .
文献结构如下图所示。
图片声讯系统
这里的cavity文献夹用于轮番测试。
修改files文献与options文献。
files文献的本体demo4.CEXE = demo4options文献的本体
EXE_INC = \ -I$(FOAM_SRC)/finiteVolume/lnInclude \ -I$(FOAM_SRC)/meshTools/lnInclude EXE_LIBS = \ -lfiniteVolume \ -lmeshTools2 源文献
本案例演示从文献中取得场数据，然后再对读取的数据进行操作。
轮番代码如下。
#include "fvCFD.H" // 声明一个函数，背面有该函数的具体收场经过// 输入本事t，空间坐标x，参考点x0以及缩放因子scalescalar calculatedPressure(scalar t， vector x， vector x0， scalar scale); int main(int argc， char *argv[]){#include "setRootCase.H"#include "createTime.H"#include "createMesh.H" Info << "读取transportProperties文献数据" << endl; // 从transportProperties文献中读取数据 // 先界说一个IOdictionary对象，其构造函数参数为一个IOobject对象 IOdictionary transportProperties( IOobject( "transportProperties"， // 字典文献名 runTime.constant()， // 字典文献地点旅途，这里为constant文献夹下 mesh， // 一个objectRegistry类对象，这里没什么用 IOobject::MUST_READ_IF_MODIFIED， // 若是文献被修改，则必须重新读取 IOobject::NO_WRITE //暗示文献为只读 )); // 界说一个dimensionedScalar变量nu // nu有量纲，其量纲dimViscosity等同于(0，2，-1，声讯系统0，0，0，0)，单元为m2/s dimensionedScalar nu( "nu"， //指命称呼 dimViscosity， //指定scalar的量纲 // 旧版块或org版块的写法：transportProperties.lookup("nu") transportProperties // com新版块写法，自动把柄称呼在字典中搜索 ); Info << "读取到的nu为：" << nu << endl; //-------------从p文献中读取压力场p-------- // 不需要查找关节字，因为通盘文献皆是对于压力场的数据 Info << "读取压力场数据" << endl; // 界说一个标量场p，无需指定量纲，因为其量纲一经在相应的文献中指定了 volScalarField p( IOobject( "p"， //指命称呼 runTime.timeName()， // 取恰现时本事 mesh， IOobject::MUST_READ， IOobject::AUTO_WRITE)， mesh); // -------从U文献中读取速率场--------- Info << "读取速率场数据" << endl; // 界说一个向量场U volVectorField U( IOobject( "U"， runTime.timeName()， mesh， IOobject::MUST_READ， IOobject::AUTO_WRITE)， mesh); // 界说一个场向量 const vector originVector(0.05， 0.05， 0.005); // 探求打算向量originVector与各网格中心的距离中的最大值 // 诈欺dimensionedVector将vector退换为具有长度量纲的向量 // mag函数探求打算向量的模 // value函数将数值退换为无量纲标量值 const scalar rFarCell = max( mag(dimensionedVector( "x0"， dimLength， originVector) - mesh.C())) .value(); // 在案例迭代经过中进行场变量探求打算 while (runTime.loop()) { Info << "Time = " << runTime.timeName() << nl << endl; // 在通盘的网格上轮回 for (label cellI = 0; cellI < mesh.C().size(); cellI++) { // 把柄自界说的函数探求打算压力值 p[cellI] = calculatedPressure(runTime.time().value()， mesh.C()[cellI]， originVector， rFarCell); } // 探求打算压力梯度，并退换为速率，由于量纲不一致，是以需要乘以一个本事量纲进行退换 // 预防在不成压缩求解器中，压力的量纲为[0 2 -2 0 0 0 0]，压力梯度量纲[0 1 -2 0 0 0 0] // 因此压力梯度乘上本事得到速率量纲[0 1 -1 0 0 0 0] U = fvc::grad(p) * dimensionedScalar("tmp"， dimTime， 1.0); // 将数据写入到文献中，鲜艳为AUTO_WRITE的场数据才可被写入 runTime.write(); } Info << "探求打算完成"; return 0;} scalar calculatedPressure(scalar t， vector x， vector x0， scalar scale){ scalar r(mag(x - x0) / scale); // 分母加1e-12是为了防患除以零 scalar rR(1.0 / (r + 1e-12)); scalar f(1.0); // 复返一个以x0为中心的正弦分别的压力 return Foam::sin(2.0 * Foam::constant::mathematical::pi * f * t) * rR*100;} 在demo4旅途诈欺wmake编译轮番。
图片
3 测试轮番轮番需要调解OpenFOAM案例初始。本案例一经事前准备了cavity案例，诈欺底下的号召初始轮番：
cd $FOAM_RUN/demo4/cavityblockMesh../demo4
初始效用如下图所示。
图片
不错参预ParaView稽查数据。
图片
本站仅提供存储职业，通盘本体均由用户发布，如发现存害或侵权本体，请点击举报。