OpenFOAM

4.5k 主题 26.2k 帖子
  • OF4.1这行代码释义,动网格问题

    4 帖子
    3k 浏览

    您好,您现在还有做SIXDOF得问题吗?有的话我们可以交流一下,谢谢啦

  • LPT-VOF coupling in OF 7

    3 帖子
    921 浏览

    @东岳 但是那个貌似是在v1906这种版本,我现在用的是OpenFoam 7, 貌似没有MPPICInterFoam

  • 关于interfoam出口压力设置的问题

    3 帖子
    1k 浏览

    @东岳 老师您好!是由于实际模型太大,为了减少计算量,截取了关键的部分作为计算域,但是由于是水下注汽,在实际模型中,该计算域最上方还有40mm的水深,我想设置一下该压力。

  • 进口通量等于出口通量的边界条件

    10 帖子
    2k 浏览

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

  • DNS求解VOF

    9 帖子
    3k 浏览

    @东岳 感谢!

  • 请问OF计算的结果可以用其他后处理软件么

    4 帖子
    2k 浏览

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

  • 11 帖子
    6k 浏览

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

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

    2 帖子
    805 浏览

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

  • reactingtwophaseeulerfoam?

    23 帖子
    3k 浏览

    alpha was solved in twoPhaseSystem.C

  • createPatch创建cyclic边界条件的问题

    7 帖子
    2k 浏览

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

  • Rhie-Chow interplation

    1 帖子
    529 浏览

    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?
    谢谢

  • 请教OpenFOAM不同网格之间场映射的原理

    2 帖子
    777 浏览

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

  • interFoam droplet

    4 帖子
    1k 浏览

    @东岳 :xinlei:

  • 请问OpenFOAM中如何在入口添加边界层厚度

    15 帖子
    3k 浏览

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

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

    3 帖子
    1k 浏览

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

  • TopoSet关于action的疑问

    2 帖子
    1k 浏览

    @D-Benjamin 求交集的意思

  • 湍流速度脉动

    7 帖子
    3k 浏览

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

  • 关于negSumDiag()的一个问题

    11 帖子
    5k 浏览

    @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 帖子
    806 浏览

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

  • 如何实现双重网格?

    6 帖子
    3k 浏览

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

  • paraView-v1906安装问题

    4 帖子
    1k 浏览

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

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

  • OpenFOAM22x编译错误

    3 帖子
    2k 浏览

    @赵一铭

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

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

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

    再次感谢!

  • interFoam解析——VOF模型推导问题(1)

    9 帖子
    2k 浏览

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

  • banana method是什么意思

    3 帖子
    897 浏览

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

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

    3 帖子
    879 浏览

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

  • -div(phi,T) 和div(-phi,T)结果为何不同?

    12 帖子
    6k 浏览

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

    微信图片_20200320103052.jpg

  • simpleFoam中的SIMPLE系列算法有区分吗?

    2 帖子
    760 浏览
  • 6 帖子
    4k 浏览

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

  • 湍流模型接口

    12 帖子
    3k 浏览

    @Samuel-Tu Lib那改到Foam_USER_Libbin

  • fvc::reconstruct 算法

    1 帖子
    483 浏览

    1.JPG
    2.JPG

  • 急求AdjointOptimisationFoam的使用说明

    3 帖子
    1k 浏览

    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 帖子
    2k 浏览

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

  • OF有网格变形或者移动的功能吗?

    2 帖子
    707 浏览

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

  • SnappyHexMesh为什么总是无法移除固体

    5 帖子
    1k 浏览

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

  • 用driftFluxFoam时出现报错

    3 帖子
    915 浏览

    @东岳 好的,谢谢老师!

  • 关于processor边界的问题

    5 帖子
    1k 浏览

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

  • 质量源项添加问题,参数rho无效

    13 帖子
    3k 浏览

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

  • 公式求助

    9 帖子
    2k 浏览

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

  • 影响parcel速度的因素

    1 帖子
    467 浏览

    如下图,再同样的位置上,我的模拟的速度再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 帖子
    1k 浏览

    @东岳 谢谢东岳老师!

  • autoPtr代码求教

    3 帖子
    978 浏览

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

  • 1 帖子
    514 浏览

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

  • 一个简单的速度压力非迭代求解器

    3 帖子
    884 浏览

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

  • 1 帖子
    463 浏览

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

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

    3 帖子
    2k 浏览

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

  • openfoam的sqr是什么意思?和sqrt有什么关系?

    3 帖子
    1k 浏览

    @东岳 谢谢东岳老师

  • salome

    1 帖子
    410 浏览

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

  • 模仿实验室的造流水槽,出口边界怎么设置

    4 帖子
    1k 浏览

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

    2.png

  • 绘制OF网格时如何设置Volume Conditions

    3 帖子
    823 浏览

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

  • kEpsilon和kOmegaSST仍未解决

    10 帖子
    2k 浏览

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