OpenFOAM

4.5k 主题 26.2k 帖子
  • 2 帖子
    2k 浏览

    你用的什么进口条件?你的设置符合常规设置么?see Ferziger and Peric page 315

  • compressibleinterfoam能够计算激波吗?

    9 帖子
    6k 浏览

    @mohui 谢谢哈,

  • 最近国内有没有组织OpenFOam的培训?

    1 帖子
    2k 浏览

    感觉有人手把手带着学OPenfoam是每个openfoam初学者的梦想。。所以最近国内有没有相关机构或个人组织培训?:big_mouth:

  • OpenFOAM中utilities命令问题

    3 帖子
    2k 浏览

    好的,谢谢李老师!

  • OpenFOAM的并行性效率问题,以及CPU选择

    1 帖子
    1k 浏览

    请问有谁测试过OpenFOAM的并行效率?在使用酷睿或至强的工作站或集群上并行最多可以达到单核性能的多少倍?我用I7-4790和I7-3920XM测试的结果,用阿姆达尔定律推测最高是2.5倍不到,而且线程数超过物理核心之后就是负优化了,感觉不太可能,无法解释那些上亿网格是怎么处理的。现在正在鼓动老板至少搞一台高性能工作站,同等价位下是主频优先还是线程优先?优先考虑酷睿还是至强?

  • 平板降膜边界条件设定(interFoam求解器)

    1 帖子
    1k 浏览

    0_1512907622643_2.jpg
    fluent中设置是这样的:
    Boundary A, B and C, D
    are set as periodic boundary conditions. Boundary E is set as no slip
    wall. Boundary F is set as opening pressure boundary.
    不知道在OF怎样设置好

    FoamFile { version 2.0; format ascii; class volScalarField; object p_rgh; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [1 -1 -2 0 0 0 0]; internalField uniform 100000; boundaryField { inlet { type cyclic; } gasinlet { type cyclic; } outlet1 { type cyclic; } outlet2 { type cyclic; } left { type cyclic; } right { type cyclic; } TWall { type fixedFluxPressure; } atmosphere { type totalPressure; p0 uniform 100000; U U; phi phi; rho rho; psi none; gamma 1; value uniform 100000; } } dimensions [0 1 -1 0 0 0 0]; internalField uniform (0 0 0); boundaryField { inlet { type cyclic; } gasinlet { type cyclic; } outlet1 { type cyclic; } outlet2 { type cyclic; } left { type cyclic; } right { type cyclic; } TWall { type fixedValue; value uniform (0 0 0); } atmosphere { type pressureInletOutletVelocity; value uniform (0 0 0); } }

    体积项用setField,但这样不对,求大神帮忙:big_mouth:

  • 19 帖子
    9k 浏览

    @hurricane007 Hi,对您说的学习笔记特别感兴趣,方便的话可否发送一贴学习交流一下,我邮箱是zhhu5634@Colorado.EDU,谢谢!

  • openfoam implicit density-based solver

    12 帖子
    7k 浏览

    @nuaawqf 今天测试了dbnsTurbFoam自带的naca0012算例,最大CFL数在100依然能保持稳定,最大自适应步长7.14e-5。更大的CFL目前暂时没有测试。感觉使用隐式时间推进之后求解器的鲁棒性提高了不少。。

  • buoyantPimpleFoam的边界条件使用问题

    4 帖子
    3k 浏览

    果然是边界条件的问题。现在已经能够运行了,等待比较结果中。。。

  • 3 帖子
    2k 浏览

    是不是我的问题没有问到点子上?请各位赐教。

  • 子模型中添加中间变量(需求解方程)

    4 帖子
    2k 浏览

    @李东岳 问题已解决,直接在.H中添加了一段简单的二分法求解方程的程序,谢谢。

  • 有没有打算基于docker在windows下开发openfoam的

    5 帖子
    3k 浏览

    @sibo 其实大约都是基于虚拟化技术来的。我是给docker装的vncserver,虚拟一个主机,在的形式windows通过vnc访问。这样server和client

  • 6 帖子
    3k 浏览

    @青瓷碎了夏天 面积分可以分两部分,第一步先使用面两边的体变量插值出面上的物理量,比如一个面k,它左右两边为i,j两个控制体,则P_k=interpolate(P_i,P_j)。这个值相当于面上的平均值,与面矢量内积就能够得到面积分了。

  • compressiblemultiphaseinterfoam计算的bug

    8 帖子
    5k 浏览

    这个问题使用fluent算貌似同样会出现这个问题,其实这时候已经出现了部分空化,处理不当产生负压。

  • 可压缩求解器中的音速的定义和离散格式

    6 帖子
    4k 浏览

    @100yearsalone 恩,最近没过来逛逛咱们这个论坛。由于OpenFOAM官方版本一直以pressure-based为主导,所以会产生很多迷惑。

  • fluent 里的porous 边界如何在openfoam中实现

    1 帖子
    1k 浏览

    fluent 里的porous 边界如何在openfoam中实现

  • 请问有人做过源造波+动网格的算例吗?

    1 帖子
    1k 浏览

    大家好,之前我发了一个帖子,大意是基于interFoam,在连续性方程中添加源项实现质量的输入和输出,进而造波:
    http://www.cfd-china.com/topic/943/添加源项后相方程求解越界-请问这种情况有人遇到过吗

    经过多次调整,质量源造波的程序已经基本上完成了。我按照上贴子的公式,把源项加入相方程中并开启了MULES修正,同时在压力方程中也添加了源项,这样求解的时候,质量源区域内的alpha值(默认为1,源区位于水下)大概在1.0*(质量输入)-0.9*(质量输出)之间波动(因为我的源项是类似于正弦波动的)。

    然而,我采用相同的方法修改压力方程,使用动网格求解器interDyMFoam的时候,上述的方法出问题了,形成了类似于空化的现象(不懂空化,希望没说错),即相体积分数变得很小,最直观的描述是在水下质量源区域会形成一个气泡,后面相体积分数的求解救完全错误了。

    我仔细看了一下求解器的文件,发现interFoam.C和interDyMFoam.C求解时主要的区别在于动网格求解程序中多了一些通量计算以及压力方程中的部分不同。但是由于之前没看过动网格的求解器,因此还是看的不太懂,请问有做过这方面的大牛吗?可否分享一下经验,为什么网格不变形的时候完全正常,但是用了动网格求解器就不行了。:big_mouth:

    感激不尽:sunglasses:

  • CFD 行业报告_Aerosciences

    3 帖子
    2k 浏览

    是啊,相信CFD在国内也会有不错的发展!

  • 在OpenFOAM中实现新的Scheme

    1 帖子
    1k 浏览

    因OpenFOAM现有格式不满足需求,需要实现新的插值格式,有没有熟悉的人说明一下思路,此外我不清楚MeshObjects文件夹下的 upwindCFCCellToFaceStencilObject和upwindCPCCellToFaceStencilObject的作用和区别是什么?谢谢

  • 求解器编译求助

    5 帖子
    3k 浏览
    surfaceScalarField cPhi("cPhi", c & mesh.sf()); fvScalarMatrix c ( fvm::ddt(c) + fvc::div(cPhi) == source. );
  • paraview做slice问题

    5 帖子
    3k 浏览

    @wwzhao 谢谢

  • 12 帖子
    7k 浏览

    @小龙 filter 搜一下就可以找到

  • 群体平衡模型积分矩方法又一大作问世

    1 帖子
    989 浏览

    继Mcgraw于1997年提出QMOM之后,Marchisio et al. 于2005年提出DQMOM,Yuan and Fox于2011年提出CQMOM,Yuan et al. 于2012年提出EQMOM,Marchisio and Fox于2013年提出ECQMOM,Fox et al. 于17天前又提出HyQMOM和CHyQMOM并发布预印本。

    https://hal.archives-ouvertes.fr/hal-01632813

  • 自定义守恒方程时出问题

    3 帖子
    2k 浏览

    @李东岳
    级测试,如您所说改为NO_READ后不再出现警告,谢谢

  • 请问,如何改变输出的结果文件夹名称呢

    3 帖子
    2k 浏览

    @random_ran 非常感谢!问题已解决

  • OpenFoam张量取值运算

    3 帖子
    2k 浏览

    @李东岳 那我该怎么提取速度梯度中的量呢?

    gradU.component(tensor::ZX)

    这样对吗?

  • overPimpleDyMFoam计算时压力正负交替出现

    2 帖子
    2k 浏览

    我记得这个好像自带的mixer2DAMI也是这样的?

  • 方程离散!fvm:div到底output出什么?

    4 帖子
    3k 浏览

    是的你说的很正确。

  • compressibleinterfoam的压力突变问题。

    5 帖子
    3k 浏览

    @mohui
    考虑过,但是刚开始压力发生这种令人不解的变化,之后压力又能够恢复正常,气泡脉动在流场产生的压力又在合理的范围内。但是我对程序进行定点调试,就发现经过我图3的代码之后,压力就突然变了。

  • 有人在REHL的服务器上装过OpenFOAM吗

    1 帖子
    1k 浏览

    5.0版本的可以正常安装,但是低版本如2.4.0安装就会报错:

    -bash: /public/openfoam/OpenFOAM-2.4.0/etc/bashrc: line 140: syntax error near u nexpected token `;&' -bash: /public/openfoam/OpenFOAM-2.4.0/etc/bashrc: line 140: ` [ "$FOAM_V ERBOSE" -a "$PS1" ] && echo "Sourcing: $1" 1>&2'
  • 5 帖子
    3k 浏览

    四面体分解确定插值权重啊,用于相间耦合或者场平均(跟用多少耦合有关)。跟重心关联,还可以表示成局部坐标的形式。

  • 9 帖子
    7k 浏览
  • 1 帖子
    1k 浏览

    请教各位老师,有没有研究涉及OpenFOAM与PFC3D的耦合的。老板让用PFC3D研究DEM-CFD coupling,感觉很少有研究这个的,目前看的资料也就是github上一个itasca公司的老师编写的一系列代码和一篇今年6月份发表的文献,资料有点少,也不知道这是否是一个正确的方向,还望各位老师指点指点,谢谢各位老师!

  • 同位网格下的压力计算问题

    7 帖子
    4k 浏览

    @mohui 是的,interface处是有间断的,但是one-field形式的控制方程体现不出来间断
    压力分布不对,看一下理论解就明白了~

  • 求助!!!关于离散格式的选择

    4 帖子
    3k 浏览

    不是少部分,是一点点负值,比如说-0.00006。如果你的方程有一些要对浓度做开方,那一点点负值都不行。

  • OpenFoam中怎么处理ICEM生成的混合网格

    3 帖子
    2k 浏览

    @李东岳 CFDEMcoupling-PUBLIC-3.0.x用这个,耦合的是OpenFoam中的PisoFoam。createPatch这个是ICEM中的,还是OpenFoam中的?

  • 1 帖子
    1k 浏览

    我从OpenFOAMl里面导出了算例网格和数据,用Fluent进行后处理时,找不到浓度、含率参数,这是什么原因????!!!。其中C和alpha.water对应Fluent参数序号如下图。0_1510822319990_QQ图片20171116165029.jpg

  • OpenFoam能否将一个patch切分?

    3 帖子
    2k 浏览

    好的,我试试

  • 请问怎么实现一个圆管在风洞中的模拟?

    4 帖子
    3k 浏览

    0_1510557038187_1510484187969-捕获-resized.jpg

  • 二相流模型中流量入口怎么设置

    1 帖子
    1k 浏览

    想用二相欧拉模型做一个固液二相流问题,下层是固体(体积分数0.6),上层是水,入口处水给一个固定流量,入口的边界条件该怎么给呢?示意图:
    0_1510474822975_1.PNG

  • 怎么才能在代码里快速找到函数的定义?

    13 帖子
    8k 浏览

    @李东岳 感谢东岳老师,非常详细的步骤。我会好好看看这块

  • rhoCentralFoam求解器算湍流

    5 帖子
    3k 浏览

    @李东岳 您好!感谢您的提议。首先rhoCentralFoam应该是可以算湍流的,我是基于这个前提,想知道除了在constant文件下turbulenceModel选择RAS或LES以后,如何改动其system文件夹下的fvsolution和fvscheme字典以及在0文件夹下新加相应的变量字典来计算湍流。不知道我的这种想法是不是正确的。

  • OpenFOAM programming学习问题

    13 帖子
    6k 浏览

    东岳大神说了开发的问题。我来说说我的看法?清楚物理模型的物理意义貌似是最基本的也是最容易出错的部分。openfoam方面可以参考The Finite Volume Method in Computational Fluid Dynamics。还有就是最笨的办法:读代码

  • 下面的fvScheme方法计算起来有什么区别

    6 帖子
    4k 浏览

    @asininno 我没有详细对比过具体有多大影响。关于非正交修正,建议你仔细看一下公式,《The Finite Volume Method in Computational Fluid Dynamics An Advanced Introduction with OpenFOAM and Matlab》这本书有很详细的说明。说白了就是采用和正交网格上一样的方法去做离散,对于非正交网格会有误差,那么再加上一些措施进行修正。
    但是还是建议做好网格,控制网格的正交性和弯曲度;毕竟修正只是一种补救措施。

  • MRF源代码的问题

    6 帖子
    3k 浏览

    @dzw05
    哦了,有空整理整理。

  • 拉格朗日梯度项方程

    4 帖子
    2k 浏览

    好吧,你说的很对,看来就是:
    \begin{equation}
    \frac{\nabla p}{\rho}=\frac{\mathrm{D} \mathbf{U}}{\mathrm{D} t}
    \end{equation}

  • 能量方程中的剪切力做功

    1 帖子
    1k 浏览

    of中多相流这块不可压好像都没有能量方程,在模仿着可压的能量方程添加时发现,剪切力做功是用湍流粘性为形参的库函数表示的,现在我不调用库函数,怎么表示这一项呢?谢谢

  • 关于计算后的数值写入的问题

    3 帖子
    2k 浏览

    @wwzhao
    谢谢赵老师的回答,我尝试下写入不同的值看是否能正常写入文档

  • 三维Eo数计算公式

    1 帖子
    865 浏览

    有关Eo数的定义:http://47.93.231.95/doku.php?id=cfd百科#eoetvoes_number

    谁那见过三维Eo数的计算公式和文献出处?OpenFOAM植入的这个公式不知道从哪来的。

  • 5 帖子
    3k 浏览

    @李东岳 额,从controlDict里的时间来看,多相计算的物理时间要远大于单相,但还是没能得到合理结果。我再算下去看看。