Skip to content

OpenFOAM

OpenFOAM交流区

5.2k 主题 31.0k 帖子
  • Openfoam中求偏导数问题

    2
    2 帖子
    1k 浏览
    李东岳
    fvc::grad(T).component(2);
  • 3 帖子
    2k 浏览
    T

    @李东岳 多谢李老师!

  • 标量输运方程添加位置

    3
    3 帖子
    2k 浏览
    田畔的风

    一般放外面就行,如果要放pimple.loop()里的话,标量方程的时间步长需要除以nOuterCorrector。

  • cent升级g++9

    1
    1 帖子
    894 浏览
    李东岳
    yum install centos-release-scl -y yum install devtoolset-9 -y source /opt/rh/devtoolset-9/enable

    OpenMPI aclocal

    cd /usr/local wget https://ftp.gnu.org/gnu/automake/automake-1.15.tar.gz --no-check-certificate tar -xzf automake-1.15.tar.gz cd aclocal-1.15 ./configure --prefix=/usr/local/automake make make install export PATH=/usr/local/automake/bin:$PATH cd openmpi-4.1.2 autogen.pl ./makeOPENMPI yum install blas-devel lapack-devel wget https://ftp.mcs.anl.gov/pub/petsc/release-snapshots/petsc-lite-3.17.2.tar.gz --no-check-certificate ./makePETSC -hypre=http://dyfluid.com/hyper-2.24.0.tar.gz
  • 关于icoFOAM的qDNS案例的问题求助

    3
    3 帖子
    2k 浏览
    U

    @李东岳 谢谢李老师,问题已经解决了,对于DNS的槽道流需要用统计数据,流场计算时间不够就很难对得上,还需要等到湍流完全发展后再开始计算平均场。

  • 如何提取时间变量实现特定时间点文件输出?

    4
    4 帖子
    2k 浏览
    chengan.wangC

    @李东岳 我找到如何输出endTime方法了,用this->owner().db().time().endTime().value()就可以。但是放在判断语句中出现如下错误:

    note: ‘Foam::word’ is not derived from ‘const std::fpos<_StateT>’ if ((this->owner().db().time().timeName()) == (this->owner().db().time().endTime().value())) /usr/include/c++/5/bits/stl_iterator.h:349:5: note: candidate: template<class _IteratorL, class _IteratorR> bool std::operator==(const std::reverse_iterator<_Iterator>&, const std::reverse_iterator<_IteratorR>&) operator==(const reverse_iterator<_IteratorL>& __x, ^ /usr/include/c++/5/bits/stl_iterator.h:349:5: note: template argument deduction/substitution failed:

    李老师能指点一下么?

  • openFoam如何直接在整个计算域内注满粒子

    3
    3 帖子
    2k 浏览
    李东岳

    感谢分享!!

  • 看个乐:renumberMesh后的魔性云图

    13
    13 帖子
    8k 浏览
    Y

    @winsway_zero 嗯嗯,我懂了,谢谢,renumberMesh是对网格的操作,setFields也要在renumberMesh之后操作,不然也会乱掉:146:

  • interFoam标量输运

    10
    10 帖子
    5k 浏览
    李东岳

    你这个哪个是进出口,重力是向下的么,

    alpha场也发一下

  • 对边界场进行赋值

    7
    7 帖子
    4k 浏览
    ZY-GONGZ

    @李东岳 :xiexie: 感谢李老师

  • 双流体方法颗粒堆积模拟

    7
    7 帖子
    4k 浏览
    Y

    @xpqiu 好的,谢谢老师,我学习一下

  • 请教圆柱绕流一维能谱制作方法

    1
    1 帖子
    933 浏览
    yfleungY

    各位带佬大家好,小弟最近在做一个雷诺数为3900的圆柱绕流。看到文献中有取圆柱后面一条线上的点做能谱分析的。

    (如取下图X/D=5;或X/D=10处的一条直线上的点的速度取平均值)

    5c64b8f6-a5e1-4ad9-a193-62c0e672ec64-image.png

    最后得出了这样子的一维能谱图

    393a09bf-c9a0-4b64-a2e5-055baca6e431-image.png

    小弟对该图片的横坐标和纵坐标表达的含义都不是很清楚。只知道在红色部分区域好像是符合湍流的-5/3幂律关系。所以想向各位带佬请教该图片横、纵坐标代表的含义,以及制作的方法。:chouchou: :chigua:

  • 论坛无法搜索

    6
    6 帖子
    3k 浏览
    L

    @李东岳 好的,感谢东岳老师:papa:

  • coded边界条件中变量的输出

    3
    3 帖子
    2k 浏览
    M

    @李东岳 感谢李老师回复。是一个volScalarField,其实就是grad(T),但是我是自定义的热流,所以采用了coded这种方式。我是按照patch("shell")索引的面,对所有属于这个patch的面进行了遍历,所以得到了一个长度为patch所含面数的数组。
    现在想把热流输出出来显示,用chtMultiRegionFoam -postProcess -func "grad (T)"后处理命令进行输出,报错如下:c132aa67-ee05-4b6d-8c90-55ed717f760e-image.png
    算例是可以正常计算的,已经计算结束了。这个错误也无从下手,所以还是尝试在coded里面直接写好把变量输出来。

  • 请教在求解器中编译函数调用的问题

    1
    1 帖子
    877 浏览
    chengan.wangC

    大家好、各位大佬好,咨询一个编译函数调用问题
    我在别的论坛找到了一个在求解器中定义函数编译调用的例子如下:

    //fun.H #ifndef fun_H #define fun_H #include "scalar.H" #include "vector.H" namespace Foam { vector TestMe(const vector& myvec); } #endif

    以及

    //fun.C #include "fun.H" Foam::vector Foam::TestMe(const vector& myvec) { Info<< "Hello from TestMe" << endl; Info<< "vector: " << myvec << endl; }

    然后在主函数里面调用 "fun.H"就行了。

    #include "fun.H" using namespace Foam; int main(int argc, char *argv[]) { vector A(1,2,3); Info<< "Hello from the main()" << endl; return 0; }

    我的思路就是直接在fun.H中写我想调用的代码。比如在“KinematicCloudI.H”文件中定义的

    template<class CloudType> inline const Foam::tmp<Foam::volVectorField> Foam::KinematicCloud<CloudType>::UPar() const { tmp<volVectorField> tUPar ( volVectorField::New ( this->name() + ":UPar", mesh_, dimensionedVector(dimVelocity, vector(0,0,0)), extrapolatedCalculatedFvPatchVectorField::typeName ) ); volVectorField& UPar = tUPar.ref(); forAllConstIter(typename KinematicCloud<CloudType>, *this, iter) { const parcelType& p = iter(); const label celli = p.cell(); UPar[celli] += (p.U())/(p.nParticle()); } UPar.primitiveFieldRef() /= 1;//mesh_.V(); UPar.correctBoundaryConditions(); return tUPar; }

    我直接移植过来,编译出现了如下错误:

    fun.H:16:41: error: no ‘const Foam::tmp<Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> > Foam::KinematicCloud<CloudType>::UPar() const’ member function declared in class ‘Foam::KinematicCloud<CloudType>’ Foam::KinematicCloud<CloudType>::UPar() const

    我尝试把“KinematicCloudI.H”里面的头文件粘上去也通过不了,应该怎么修改以及加一些什么内容呢?
    谢谢

  • 今天 win下的paraview打开出错怎么回事

    4
    4 帖子
    3k 浏览
    J

    @李东岳 是的李老师

  • pFinal 关键字的实现过程

    3
    3 帖子
    2k 浏览
    Y

    @李东岳 好的,李老师,原来是这样

  • 沸腾过程中气泡从出口出去的边界条件

    3
    3 帖子
    2k 浏览
    李东岳

    我也发现了这个情况,并且写在了ns方程笔记里面。不过没有详细研究。我觉得是因为多相流密度差异导致的压差计算不同导致的。有空是应该研究一下。

  • 虚拟机双版本OpenFOAM怎么切换使用

    3
    3 帖子
    2k 浏览
    U

    @李东岳谢谢东岳老师,已经搞定了在求解器窗口输入 of2106 就切换回V2106版本了:high:

  • 后处理输出Qdot

    2
    2 帖子
    2k 浏览
    C

    @adcpk 您好,请问这个问题后面解决了嘛,请教一下