OpenFOAM

4.5k Topics 26.2k Posts
  • 时均值,面均值?

    1 Posts
    537 Views

    各位大佬,of里面除了时均值有面均值吗?就是假如有个1m长的水管,网格dx=0.001m,能不能在某个时间步下输出1000个值代表着一共1000个剖面上的某个物理量的面平均值?如果要自己编程解决,那下面这25284个网格怎么和实际坐标联系起来?谢谢

    dimensions [0 0 0 0 0 0 0];//整个计算域所有网格的空气相分数 internalField nonuniform List<scalar> 25284 ( 0.878566 0.854193 0.876398 0.860644 0.879329 0.859877 0.878446 0.860776
  • 求通量问题

    4 Posts
    1k Views

    @东岳 谢谢,改完您建议的就没有那个警告了

  • 3 Posts
    1k Views

    把原来流场的小眼睛点开

  • 判断壁面类型

    1 Posts
    432 Views
    #include "fixedValueFvPatchFields.H" #include "zeroGradientFvPatchFields.H" #include "wallPolyPatch.H" if isA<wallPolyPatch> ( mom[0][0].boundaryField()[P] ) else if ( isA<fixedValueFvPatchScalarField> ( mom[0][0].boundaryField()[P] ) )
  • 请教?如何控制网格变形细节

    3 Posts
    871 Views

    @东岳 嗯嗯 好的 谢谢啦:shangxue:

  • rhoCentralFoam计算发散问题

    2 Posts
    1k Views

    你好,我也想在rhocenteralFoam中添加燃烧,做超音速和跨音速燃烧,不止能够请假该如何将燃烧模块添加到rhocenteralFoam中?期待您的答复。

  • 如何直接利用OF求解小矩阵?

    7 Posts
    2k Views

    @Samuel-Tu 自己写一个也不难,OF的数学库都是针对GeometricField的,估计是没有

  • OpenFOAM识别层流过渡流和湍流

    9 Posts
    2k Views

    @Samuel-Tu 哦哦 有点启发了 多谢了

  • 1 Posts
    355 Views

    大家好,我最近在使用interfoam求解器通过openmpi算一个算例,但是每次算一会儿就会出现以下错误提示:
    e86f3c41-75e6-4103-a136-858021b7b253-image.png
    有人能解答一下么,万分感谢。

  • v1812 decmopose

    1 Posts
    494 Views

    在v1812里,使用了 cyclic boundary

    preservePatches (LEFT RIGHT) (TOP DOWN);

    decompose时报错:

    --> FOAM Warning : Using 'preservePatches' constraint specification. This constraint keyword is deemed to be 54 months old. --> FOAM FATAL IO ERROR: Entry 'preservePatches' has 4 excess tokens in stream 8(( LEFT RIGHT ) ( TOP DOWN ))

    怎么回事:chouchou:

  • mapFields如何实现并行映射呢

    20 Posts
    8k Views

    请问先用粗网格跑至充分发展湍流再用mapfield映射到细网格上,那粗网格可以有多粗?
    我目前DNS网格是7500W,那我可否就用五六百万的网格先跑至充分发展?
    另外既然是用于DNS的初始场,那粗网格也是设成laminar跑是吗?谢谢。

  • 5 Posts
    1k Views
  • 求助一篇文献关于LEMOS-2.4.x的

    1 Posts
    480 Views

    Kornev N, Shchukin E, Taranov E, et al. Development and implementation of inflow generator
    for LES and DNS applications in OpenFOAM. Proceedings of the Open Source CFD International
    Conference; 2009 Nov 12–13; Barcelona, Spain; 2009

  • 未知非对称矩阵求解器

    6 Posts
    2k Views

    @人来人往 relaxationFactors实在PIMPLE下面的,稍稍调整一下,另外删除一些没必要的东西,GAMG不需要preconditioner。
    具体上网找找每个求解器的定义吧,都需要什么,什么是default的,到底哪个参数需要修改

    GAMG generalised geometric-algebraic multi-grid. { solver GAMG; tolerance 1e-07; relTol 0; cacheAgglomeration true; switch specifying caching of the agglomeration strategy (default true). nCellsInCoarsestLevel 2; approximate mesh size at the most coarse level in terms of the number of cells (default 10). directSolveCoarset false use a direct solver at the coarsest level (default false). agglomerator faceAreaPair; mergeLevels 2; keyword controls the speed at which coarsening or refinement is performed. simple mesh use 2 preSweepsLevelMultiplier multiplier for the the number of sweeps between each coarsening level (default 1). maxPreSweeps maximum number of sweeps as the algorithm is coarsening (default 4). nPostSweeps number of sweeps as the algorithm is refining (default 2). postSweepsLevelMultiplier multiplier for the the number of sweeps between each refinement level (default 1). maxPostSweeps maximum number of sweeps as the algorithm is refining (default 4). nFinestSweeps number of sweeps at finest level (default 2). smoother DILUGaussSeidel; nPreSweeps 0; number of sweeps as the algorithm is coarsening (default 0). nPostSweeps 2; minIter 1; }
  • 6 Posts
    2k Views

    @yfclarkundefined symbol: _ZN4Foam23ControllerStateProvider8positionE错误求助 中说:

    ic或sta

    我是用的静态数据变量static 然后也是public成员

  • 7 Posts
    2k Views
  • 12 Posts
    5k Views

    @linhan-ge 请问您考虑了什么作用,现在颗粒润湿是不是还没人实现

  • OpenFOAM求解声学问题?

    1 Posts
    587 Views

    想要用OpenFoam计算螺旋桨或者风扇的声传播问题。所以想问一下各位在OpenFoam中有没有求解FW——H方程的求解器?

  • pyFoam实时监测残差

    5 Posts
    3k Views

    网上的有些命令可能不对,我照着安装了,没法运行

    pyFoamPlotWatch.py cavity.log,而正确的命令应该为:

    pyFoamPlotWatcher.py cavity.log

  • 7 Posts
    2k Views

    @wwzhao
    真的非常感谢。
    抱歉, 因offline没及时发现您的回复。

  • 6 Posts
    1k Views

    @东岳 老师好,如何在指定时间范围内对新变量求平均,我尝试在controldict添加fieldaverage,但是计算时没有Mean和prime2mean文件输出,应该如何解决?

  • 1 Posts
    549 Views

    TIM图片20191111100118.png

  • UList使用

    3 Posts
    831 Views

    @东岳 对,我后来是用List<type>解决的。最后用List来初始化UList

    List<label> IBMarks(50,(label)1); labelUList affectedList(IBMarks); Info << affectedList<<endl; Info << affectedList[25]<<endl;

    结果如下

    50{1} 1

    v_感觉是类似数组头的指针,运算符[]可能被重载了,可以直接取值,就像我上面的affectedList[25],输出为1。
    另外我原来

  • openfoam并行出现错误。

    1 Posts
    400 Views

    大家好,我最近在用openfoam中的interfoam进行并行计算时出现了以下的错误。

    068f8746-20e6-4aab-8566-00ba4f1862d0-image.png

    我看相关论坛之后说是初值设置有问题或者求解器有问题,但修改之后还是出错误,希望大家能够解答一下,万分感谢。

  • 通量重构算法

    1 Posts
    389 Views

    看同济硕士开题有人开通量重构算法的课题。。这个大家有了解吗。说是高阶的算法,实现简单。想了解一下。另外OpenFOAM可能实现吗。。

  • openfoam并行错误求助

    2 Posts
    465 Views

    不好意思大家,图片未能上传成功,浪费大家时间了,非常抱歉。

  • 5 Posts
    1k Views

    @yfclark 方便私聊么?我也想咨询一下你关于平面 测parcel的直径问题,
    我邮箱在资料

  • 求助。

    12 Posts
    2k Views

    @win 如果涉及到parcel到边界的话,注意cell size 和parcel size的关系~

  • DPMFoam计算粒子分布

    7 Posts
    2k Views

    @ZYj

    interpolationSchemes//添加离散格式 { rho.air cell; U.air cellPoint; mu.air cell; DUcDt cellPoint; } pressureGradient//添加虚拟质量力 { U U.air; } virtualMass { Cvm 0.5; }
  • 13 Posts
    2k Views

    @yfclark 真是太不好意思,又来打扰您,今天我又尝试了bluecfd的方法,按照他给的步骤,wmake完成后把库也链接到算例下的controlDict了,按照他的描述应该在每个时间步下/lagrangian/kinematicCloud`内出现一个position.orig但是无论我怎么尝试,也都出现不了,终端显示也是读取了这个命令了

    writeCloudOldStyle writeCloudOldStyle1: applying to clouds: kinematicCloud

    这中间哪里出了问题呢???

  • 关于算例问题

    9 Posts
    1k Views

    @东岳 解决了,确实是内存的问题,把并行编译核数减少就可以了

  • 4 Posts
    1k Views

    求解器编译成功了,,,
    将Options中动态链接库中的:-lincompressibleRASModel 去掉就编译通过了。

  • 2 Posts
    743 Views

    这个看不出什么问题啊

  • 4 Posts
    1k Views

    主要是D的先进性起了决定性作用

  • 边界条件设定问题?

    8 Posts
    5k Views

    你好,我也遇到了你一样斗问题,请问你最后是怎么设置成功的呢,是给定压力进出口边界条件,压力固定值进口,固定值出口,速度全部零法向梯度吗

  • 速度、压力边界条件

    1 Posts
    622 Views

    请问入口处有压力阀控制气压,已知压力值,入口边界条件该怎么设置?压力边界可以用fixValue,速度边界用哪种啊?
    出口跟外界连通,出口压力和速度边界又该如何设置呢?

  • AMI 和WaveDyMFOAM

    1 Posts
    586 Views

    用icem合并生成网格,之后在openfoam里定义ami边界条件(水槽中间有个封闭的矩形边界,里面是动网格)后,用waveDyMfoam求解器求解时,原本设置好的探针波高检测器(cellpointface)无法产生实际效果,运行后没有生成对应的surfaceElevation文件

  • 10 Posts
    9k Views

    前辈,您好,在这个设置区域的浓度方面这个fvoptions与setfields的用法上,fvoptions能实现随时间的变化,其他还要什么呢。

  • DPMFoam后处理

    1 Posts
    547 Views

    麻烦各位大佬看一下,小白问的可能比较蠢的问题:mihu:
    Q1:最近在用DPMFoam做颗粒流,想在paraview中显示出粒子的轨迹有没有相关的例子可以参考下;
    Q2:在计算管道流动过程中,流体还未流出,为什么在outLet处就可以看到相应的速度,这个地方的速度不应该是0吗

  • OpenFOAM7卸载

    4 Posts
    3k Views

    @东岳 按上诉方法已经成功删除了,还成功重装了of1806。删除主要是想腾出空间。

  • SnappyHexMesh中的使用问题

    10 Posts
    4k Views

    @秋山结月
    1.对的,最好另外再建立新的几何体,用这个新的几何体包裹要refine的局部区域。
    2.这个问题倒是没有遇到过,是不是stl文件的问题?
    3.snap结束后不应该直接没掉啊,即便是没有需要施加边界层的区域,也会有“No layers to generate ...”,这样的信息。看看你后面有没有这样的提示信息?

  • 3 Posts
    826 Views

    感谢楼主分享

  • 7 Posts
    2k Views

    @星星星星晴
    谢谢你,我已经把这个错误解决了,是因为我在files文档的第一行kOmega231/kOmega231.C应该写成kOmega231.C

  • 动态数组的使用

    2 Posts
    551 Views

    new一次,delete一次,完美

  • cyclic和cyclicAMI的区别

    2 Posts
    1k Views

    cyclic和cyclicAMI的value的值没有意义,这是类的成员,必须赋值不管有没有用,随便飞翔吧老铁!

  • 5 Posts
    1k Views

    现在明白了,谢谢东岳老师讲解

  • 6 Posts
    1k Views

    @东岳 李老师,关掉考虑自转之后出现倒流现象更严重了

  • PIMPLE算法的使用

    8 Posts
    3k Views

    @星星星星晴 好的谢谢,我目前网格有7500W,质量0.72以上,结构是简单的圆管,在算DNS,ncorrector设置的2,您觉得可以吗?

  • 1 Posts
    495 Views

    大家好,我参考openfoam自带的求解器,在求解器中添加自适应网格的功能。我的算例网格总数是50000,运行算例的时候,遇到下面的错误
    1.png
    以下是reactMinMax.H 17-30行

    17 const label tautableDim = 3; 18 19 forAll(p,i) 20 { 21 List<scalar> lookupCoordinates(tautableDim,0.0); 22 lookupCoordinates[0]=fH[i]; 23 lookupCoordinates[1]=pHI[i]; 24 lookupCoordinates[2]=THI[i]; 25 26 const List<scalar> &LOG_t_ign_temp = tignTable.lookUpBounded(lookupCoordinates,false); 27 tIgnHI[i] = min(max(1e-14,pow(10,LOG_t_ign_temp[0])),1e9); 28 29 } 30 tIgnHI.correctBoundaryConditions();

    请教大家问题在哪呢

  • 8 Posts
    6k Views

    请问PBiCGStab在更稳定,收敛速度更快的情况下,会不会牺牲精度呢?还是两者都会给出正确结果?我目前正在算DNS,对于计算准确性比较关心。谢谢。