OpenFOAM

4.5k Topics 26.2k Posts
  • 10 Posts
    2k Views

    @东岳 嗯嗯嗯,我去在这个方向上面试一试,这些李老师

  • DNS求解VOF

    9 Posts
    3k Views

    @东岳 感谢!

  • 4 Posts
    2k Views

    我转换ensight完成之后发现,只能读出几何模型,里面并没有数据,这是怎么回事

  • 11 Posts
    6k Views

    您好!我遇到类似问题,准备采用上述的第二种方法即topoSet + createPatch,但不太清楚具体应该如何操作,希望能指点一下,谢谢!@赵一铭 @小龙

  • 怎么把mapfields的code写进solver里?

    2 Posts
    804 Views

    您好,这里两种网格求解的方程形式一样吗?您是研究多尺度模拟的吗?我也想尝试OpenFoam里植入一种多尺度算法,是否方便交流一下

  • reactingtwophaseeulerfoam?

    23 Posts
    3k Views

    alpha was solved in twoPhaseSystem.C

  • 7 Posts
    2k Views

    @一颗橙子 您好,想问一下您这个问题解决了嘛?后续是如何操作的。。我似乎遇到了和您完全一样的问题

  • Rhie-Chow interplation

    1 Posts
    529 Views

    Original RC interplation
    No time dependency, No relaxation of the velocity field:

    // Rhie-Chow interplation phi = (U_avg_f & mesh.Sf()) - ( (DUf*( gradp_f - gradp_avg_f)) & mesh.Sf() );

    Time dependency, relaxation of the velocity field:

    // Rhie-Chow interplation phi = (U_avg_f & mesh.Sf()) - ( (DUf*( gradp_f - gradp_avg_f)) & mesh.Sf() ) + (scalar(1) - URFU)*(phi.prevIter() - (U_avg_prevIter_f & mesh.Sf())) + DTf*( phi_old - (U_old_f& mesh.Sf()));

    OpenFOAM
    No time dependency, with or without relaxation of the velocity field (check simpleFoam solver):

    surfaceScalarField phiHbyA("phiHbyA", fvc::interpolate(HbyA) & mesh.Sf());

    Questions: Where is the second term of the original RC? Where is the third term of the original RC if taking the relaxation of the velocity field into consideration?

    Time dependency, relaxation of the velocity field (check pisoFoam solver):

    surfaceScalarField phiHbyA ( "phiHbyA", (fvc::interpolate(HbyA) & mesh.Sf()) + fvc::ddtPhiCorr(rAU, U, phi) );

    Questions: ddtPhiCorr takes the function of (the third term+the fourth term of the original RC)? I guess NO. Let's recall the definition of ddtPhiCorr:

    tmp<fluxFieldType> ddtPhiCorr ( new fluxFieldType ( ddtIOobject, rDeltaT * this->fvcDdtPhiCoeff(U.oldTime(), phiAbs.oldTime()) * ( fvc::interpolate(rA*rho.oldTime())*phiAbs.oldTime() - ( fvc::interpolate(rA*rho.oldTime()*U.oldTime()) & mesh().Sf() ) ) ) );

    Look like ddtPhiCorr returns the following codes:

    fvcDdtPhiCoeff*DTf*( phi_old - (U_old_f& mesh.Sf()));

    Last question: How to take the effect of relaxation of the velocity field into ddtPhiCorr?
    谢谢

  • 2 Posts
    775 Views

    meshToMeshInterpolate.C 可以看看这个内容,对,就是插值

  • interFoam droplet

    4 Posts
    1k Views

    @东岳 :xinlei:

  • 15 Posts
    3k Views

    @东岳
    老师这是我目前算出来的一些结果:
    链接: https://pan.baidu.com/s/163fY1uXF_sCCcN9nt1MiCQ
    提取码: w3gp
    网格量很少的,我现在只是想通过调试把路给走通,方便以后的使用

  • 《无痛苦N-S方程笔记》求指导

    3 Posts
    1k Views

    @东岳 谢谢老师,这个问题当时思维定式了,把U_p提出来,立马就懂了。

  • TopoSet关于action的疑问

    2 Posts
    1k Views

    @D-Benjamin 求交集的意思

  • 湍流速度脉动

    7 Posts
    3k Views

    @小考拉 你好这个问题我和我师兄讨论过,你不可能输出瞬时的脉动速度场,因为湍流模型根本就没有求解脉动速度,只能输出二阶矩这种统计值,建模的话可以假设三个方向的速度大小相同,给出脉动速度的大小,但是不能给出速度的方向

  • 关于negSumDiag()的一个问题

    11 Posts
    5k Views

    @wwzhao

    是的,OpenFOAM中面的方向定义永远是从小编号指向大编号,Gauss定理中面的方向指向外。也许把这些面分为大面和小面区分一下更容易理解。

    对于编号大于$P$的单元$U$,他们之间面上的量为:

    $$\phi_ f = \varpi \phi_{P} + (1 - \varpi) \phi_{U}$$

    $$\dot { m } _ { f } = ( \rho \mathbf { u } ) _ { f } \cdot \mathbf { S } _ { f }$$

    对于编号小于$P$的单元$L$,他们之间面上的量为:

    $$\phi_ f = (1 - \varpi) \phi_{P} + \varpi\phi_{L}$$

    $$\dot { m } _ { f } =- ( \rho \mathbf { u } ) _ { f } \cdot \mathbf { S } _ { f }$$

    $$
    \sum _ { f \sim n b ( P ) } ( \rho \mathbf { u } \phi ) _ { f } \cdot \mathbf { S } _ { f }
    = \sum _ { N \in L ( P ) } \left( - \dot { m } _ { f } \phi _ { f } \right) + \sum _ { N \in U ( P ) } \dot { m } _ { f } \phi _ { f }
    $$
    $$
    = \sum _ { N \in L ( P ) } - \dot { m } _ { f } \left[ ( 1 - \varpi ) \phi _ { P } + \varpi \phi _ { N } \right] + \sum _ { N \in U ( P ) } \dot { m } _ { f } \left[ \varpi \phi _ { P } + ( 1 - \varpi ) \phi _ { N } \right]
    $$

    $$
    = \left( \sum _ { N \in L ( P) } - \dot { m } _ { f } ( 1 - \varpi ) + \sum _ { N \in U ( P ) } \dot { m } _ { f } \varpi \right) \phi _ { P } + \sum _ { N \in L ( P ) } - \dot { m } _ { f } \varpi \phi _ { N } + \sum _ { N \in U ( P ) } \dot { m } _ { f } ( 1 - \varpi ) \phi _ { N }
    $$

    这样可以看出lower存的是 $- \varpi _ { f } \dot { m } _ { f } $, upper存的是$\dot { m } _ { f } ( 1 - \varpi ) _ { f } $.并且diag存的是
    $$a _ { P } = -[\sum _ { N \in L ( P ) } \dot { m } _ { f } ( 1 - \varpi ) + \sum _ { N\in U ( P ) } (-\dot { m } _ { f } \varpi )]$$

  • 有人用过of里的fft吗

    2 Posts
    805 Views

    写个类调用下fftw的函数接口就好了

  • 如何实现双重网格?

    6 Posts
    3k Views

    @东岳 您好,东岳老师。我最近才看到您之前发的这个动图。这个是用mppicinterfoam做出来的吗?我尝试了很多次都不能堆积出您的这种沙床的结果。您能简单说一下您怎么做的吗?等您回复!

  • paraView-v1906安装问题

    4 Posts
    1k Views

    http://www.cfd-china.com/topic/3182

    有一些介绍,你试试二楼的

  • OpenFOAM22x编译错误

    3 Posts
    2k Views

    @赵一铭

    感谢分享!最近安装低版本的OpenFOAM做原有的程序二次开发,用到了这个语句~~,

    在此分享一个基于不同Linux系统的OpenFOAM不同版本安装总结,亲测OF211,OF231 再 Ubuntu 16.04 安装成功无Error!

    https://openfoamwiki.net/index.php/Installation/Linux

    再次感谢!

  • 9 Posts
    2k Views

    @东岳 谢谢老师,太久没看,给忘了

  • banana method是什么意思

    3 Posts
    895 Views

    @Cp_Zhao 明白了,感谢指点!!

  • 一个有意思的超音速算例

    3 Posts
    879 Views

    @gyzhangqm 是否有足够的数据支撑?有的话可以。

  • 12 Posts
    6k Views

    这是div(phi,T)与-div(phi,T)等的解释

    微信图片_20200320103052.jpg

  • 2 Posts
    760 Views
  • 6 Posts
    4k Views

    可以直接在of中sample提取数据

  • 湍流模型接口

    12 Posts
    3k Views

    @Samuel-Tu Lib那改到Foam_USER_Libbin

  • fvc::reconstruct 算法

    1 Posts
    483 Views

    1.JPG
    2.JPG

  • 3 Posts
    1k Views

    AdjointOptimisationFoam user Mannual 链接:
    https://www.openfoam.com/documentation/files/adjointOptimisationFoamManual.pdf
    建议使用之前看一下连续伴随方法中的相关公式推导、伴随变量的边界条件推导以及敏感度分析,目前仅有伴随S-A模型,适用于气动外形优化。目前OpenFOAM中实现连续伴随方法的主要难点在于动量伴随方程中ATC(adjoint transpose convection)项的处理,由于该项会导致矩阵刚度较大,如果处理不当会导致伴随方程求解发散。运用伴随方程可以方便的进行各种目标函数和约束的敏感度分析,包括表面敏度和体积敏度,分别对应形状优化和拓扑优化。欧洲相关学者已经研究得非常透彻了,未来应该会在openfoam新版本慢慢放出。已经商业化的是engys公司旗下的Helyx-adjoint产品。

  • 安装paraview时,./makeParaView出错

    6 Posts
    2k Views

    @队长别开枪 队长您好,您可否帮我解答一下这个帖子里问的问题,就是为什么snappyHexMesh画网格始终无法移除固体,感谢队长!
    http://www.cfd-china.com/topic/3408/snappyhexmesh为什么总是无法移除固体

  • 2 Posts
    707 Views

    就,就是普通的动网格,dynamicMesh

  • 5 Posts
    1k Views

    画网格过程中出现这个提示是什么意思
    1.png

  • 用driftFluxFoam时出现报错

    3 Posts
    914 Views

    @东岳 好的,谢谢老师!

  • 关于processor边界的问题

    5 Posts
    1k Views

    @东岳
    我现在是对所有网格做遍历,遍历每个网格上的面,判断是内部的面还是边界的面,如果是边界的面的话,通量值就会和另一块processor上对应的网格上的值相关,目前就是不知道怎么去找这个网格。

  • 13 Posts
    3k Views

    @dzw05
    不好意思, 刚仔细看了一遍看到了, 谢谢!

  • 公式求助

    9 Posts
    2k Views

    @东岳
    好的我试一下,谢谢李老师

  • 影响parcel速度的因素

    1 Posts
    466 Views

    如下图,再同样的位置上,我的模拟的速度再0.01以后就远小于实验数据。
    受力仅考虑drag和gravity。stochastic dispersion model。
    但是也有实验室同事写的collision model(已经在他们的毕业论文中验证过了)

    2-way coupling
    flow和实验数据匹配。
    injection 中的parcel 的size 和速度也和实验数据匹配。

    问题就是在径向一段距离以后,parcel的动量损失很多。
    不知道大佬们有没有什么相关经验?

    4ebe1f80-272c-4107-ab8d-ae747f0e738f-image.png sudu 的

  • fvSolution下的cache起什么作用?

    3 Posts
    1k Views

    @东岳 谢谢东岳老师!

  • autoPtr代码求教

    3 Posts
    977 Views

    @羽之下 这貌似是一个类型转换函数,搜到了类似帖子https://bbs.csdn.net/topics/390537437

  • 1 Posts
    513 Views

    turbulenceProperties中object为RASProperties或turbulenceProperties有什么区别?

  • 3 Posts
    884 Views

    参考Hirsch 2007: 12.4.1 Basic Approach of Pressure Correction Methods

  • 1 Posts
    462 Views

    请问OpneFoam复杂模型的多孔介质流动是否有成功的案例或者相关资料供参考。
    :xinlei: 心累,我啥也不懂,望各大神各显神通,感恩戴德。
    邮箱:2685306635@qq.com。:zoule:

  • OpenFOAM中暂停后如何继续计算

    3 Posts
    2k Views

    @东岳 谢谢东岳老师:ok3:

  • 3 Posts
    1k Views

    @东岳 谢谢东岳老师

  • salome

    1 Posts
    410 Views

    大家好,请问各位老哥有使用过salome前处理的吗,最近打算研究一下,可以交流交流:chouchou:

  • 4 Posts
    1k Views

    @东岳 李老师,我试了一下,把出口和入口速度设置的速度一样,在开始的时候水池的水就不稳定,在入口处水面上升,在出口处水面下降。得好久才能接近一种平衡,而且导致液面晃荡。
    1.png

    2.png

  • 3 Posts
    823 Views

    哦哦,那应该就对网格没什么影响,可能是Pointwise这个软件里面才有的设置,谢谢李老师:chitang:

  • kEpsilon和kOmegaSST仍未解决

    10 Posts
    2k Views

    插值从表格里面做热物性吧?那估计是超临界啥的。算例本身算起来可能是会有点不好收敛。你层流验证的时候好收敛么,还是也很难收敛?如果层流很好收敛但湍流不收敛,那大概率是湍流的边界问题;如果层流本身也很难收敛,那就可以去找别的原因了。其实你可以试试,同样的设置下把热物性换成普通的(理想气体之类的)看看能不能算,这样也能排除一些问题。
    单看log文件有两点我觉得值得考虑。1、你设置了pressureControl,限定了压力的最大最小值,你本身的工况是在这个压力范围之内么? 2、这个可能性比较低。你可以把SIMPLE的convergence criteria调小一些,看起来你的速度方程解了1步就收敛了,对线性求解器来说好像太容易了点。

  • fvc::grad(p)的计算问题

    2 Posts
    967 Views

    零梯度边界,边界上的值一直等于相邻网格单元的值。

  • Ubuntu 下 OpenFOAM代码的调试

    6 Posts
    2k Views

    @刘雄国 请教一下, OpenFOAM 的结果怎么一步一步的打印出来那?

  • unknown type这如何排除?

    2 Posts
    766 Views

    问题找到了,是newFunction1下的TrainingExp.H文件中TypeName(" ");的内容没对上。
    表现为:可用function1类型中包含相似名称的类型,边界条件与该类型存在字符偏差。
    相关链接:https://www.cfd-online.com/Forums/openfoam-programming-development/133466-unknown-patchfield-type.html