Skip to content

OpenFOAM

OpenFOAM交流区

5.4k 主题 31.9k 帖子
  • 17 帖子
    13k 浏览
    O

    @vortex 是的 我课题组的代码也是您说的这种 是两部映射法里面中间步用f更新了速度再求压力泊松方程方程再更新最终速度 然而改到openfoam发现出了点小问题:136:

  • 关于openfoam中使用的双流体模型

    1
    1 帖子
    1k 浏览
    A

    在看《Computational Fluid Dynamics of Dispersed Two-Phase Flows at High Phase Fractions》的时候碰到一个不太理解的地方。
    一般的相间动量传递可以写成方程(3.10)
    1.png
    为了模拟phase inversion, 方程(3.11)加入了b相变成离散相时的曳力、升力等。
    但对于方程(3.11)的曳力相,如果仍然是a相为离散相的情况,$\alpha_a$→0 ,$ f_a$→0,
    那曳力相前面的部分不就变成$3/4\alpha_a\alpha_b(f_b\frac{C_{db}\rho_a}{d_b})$,和(3.10)里的不一样,b相变成分散相时候应该才是这个形式,虚拟质量力也是一样的情况。 升力也对不上。

    2.png

  • 网格变细后液面不正常

    17
    17 帖子
    12k 浏览
    J

    突然想到这个帖子,通过调整边界层厚度可以达到收敛了。

  • 场计算类型的问题

    9
    9 帖子
    6k 浏览
    X

    @fangyuanaza
    :146:

  • 关于interFoam或interMixingFoam空气密度的影响

    3
    3 帖子
    3k 浏览
    chengan.wangC

    @李东岳 李老师好,空气密度设为1000的情况流场数据在0.0489范围之内,虽然顶部看起来不好看。
    相分数初始场:

    #include "../include/initialConditions" defaultFieldValues ( volScalarFieldValue alpha.air 0 volScalarFieldValue alpha.saltwater 1 volScalarFieldValue alpha.water 0 ); regions ( boxToCell { box ($xyz_min 0.1 $xyz_min) ($xyz_max $xyz_max $xyz_max); fieldValues ( volScalarFieldValue alpha.air 1 volScalarFieldValue alpha.saltwater 0 ); } boxToCell { box ($xyz_min 0 $xyz_min) ($xyz_max 0.1 $xyz_max); fieldValues ( volScalarFieldValue alpha.saltwater 0 volScalarFieldValue alpha.water 1 ); } );

    速度场

    #include "../include/initialConditions" dimensions [0 1 -1 0 0 0 0]; internalField uniform ($velocity 0 0); boundaryField { inlet { type fixedValue; value uniform ($velocity 0 0); } outlet { type inletOutlet; inletValue uniform (0 0 0); value $internalField; } atmosphere { type pressureInletOutletVelocity; value uniform ($velocity 0 0); } bottom { type slip; } front { type symmetryPlane; } back { type symmetryPlane; } }

    压力场

    dimensions [1 -1 -2 0 0 0 0]; internalField uniform 0; boundaryField { inlet { type fixedFluxPressure; value uniform 0; } outlet { type zeroGradient; } atmosphere { type prghPressure; rho rho; p uniform 0; value uniform 0; // optional initial value } bottom { type fixedFluxPressure; value uniform 0; } front { type symmetryPlane; } back { type symmetryPlane; } }

    这个问题卡住好久了,但没有思路调整哪些参数,希望李老师能指点一下

  • 16 帖子
    12k 浏览
    R

    @李东岳 , 李老师您好,我已经在chtMultiRegionFoam中植入了第一张图所示的拉普拉斯及泊松方程,

    QQ截图20210619102650.jpg

    QQ截图20210619102636.jpg

    实现了对电势在流固区域的耦合求解,但电荷的求解似乎不对,然后电场力的求解结果也和参考文章里面差的太多,但我边界设置理论上是符合文章要求的,第二张图和第三张图是我和参考文章计算结果的对比。参考文章及我的算例和编译的求解器链接放在这里:https://pan.baidu.com/s/1U8NBBon8CFVvDoF9Fsi9UA 提取码:mppz
    麻烦李老师能抽空看一下,检查一下问题,不胜感激。一个头俩大:136: :136:

    QQ截图20210619095212.jpg QQ截图20210619095249.jpg

  • 气泡在静水中旋转上升

    2
    2 帖子
    2k 浏览

    @conghongchuan 提高一下热度,我可不可以理解成:把气泡替换成旋转机械,就是旋转速度了。

  • 有没有大佬做CFDEM Coupling的交流一下

    1
    1 帖子
    1k 浏览
    L

    关于网格尺度小了以后interpolation和mapping的问题

  • 后处理如何获得破碎液滴的粒径和速度分布?

    8
    8 帖子
    7k 浏览
    李东岳

    这个算例太酷炫了 小液滴都算出来了 还特别圆

  • 关于codedFixedValue使用

    3
    3 帖子
    3k 浏览
    李东岳

    @chengan-wang 在 关于codedFixedValue使用 中说:

    我定义了ymax,之后如果我使用if ( Cf[faceI].y() > ymax ) 设置相分数,时间步长越来越小。但如果改成if ( Cf[faceI].y() > 0.1 ) 一个具体数值,情况就好很多。我的使用方法是否有问题呢?

    理论上不应该有任何问题 :jingya:

  • 使用foamDictionary修改参数的问题

    3
    3 帖子
    3k 浏览
    W

    @浪迹天大 这样嘛。。。好的,谢谢!

  • 1 帖子
    1k 浏览
    五好青年

    如题,想用Foam-extend里的浸没边界法(Immersed Boundary Method)求解器,解决水与弹性杆的流固耦合数值模拟,有人有这方面经验吗?

  • 4 帖子
    5k 浏览

    @static 你好,我最近也在做三维飞机的相关数值模拟,可否给个联系方式,请教一些问题。谢谢,这是我的qq:3310792254

  • CFDEM虚拟机

    3
    3 帖子
    2k 浏览
    J

    @李东岳 很抱歉才看到老师的回复。感谢老师,我已经自己解决啦:xinxin:

  • 关于paraview显示分辨率的问题

    2
    2 帖子
    2k 浏览
    李东岳

    图挂了

  • OF计算动网格的问题

    3
    3 帖子
    3k 浏览
    W

    @李东岳 东岳老师,我在dynamicMeshDict里设置了constrints和restraints两个条目,控制圆柱仅在y方向移动,不转动,如下:

    constraints { yLine { sixDoFRigidBodyMotionConstraint line; centreOfRotation (0 0 0.5); direction (0 1 0); } noRotation { sixDoFRigidBodyMotionConstraint orientation; } } restraints { verticalSpring { sixDoFRigidBodyMotionRestraint linearSpring; anchor (0 0 0.5); refAttachmentPt (0 0 0.5); stiffness 2.1055; damping 0; restLength 0; } }

    不知道您说的是不是这个意思?

  • 压力泊松方程植入问题

    3
    3 帖子
    3k 浏览
    strive_DUTS

    @李东岳 感谢李老师百忙之中指点!这么来看correction(fvm::ddt(p))就代表前后时间步压力差值的时间变化率啦

  • 关于显示物体质量的问题

    1
    1 帖子
    1k 浏览

    各位大佬好,想请教一下,我想显示一个密度、体积不断变化的物体的质量,请问大家有什么好办法吗?

  • interPhaseChangeFoam求解器中的一些疑问?

    13
    13 帖子
    16k 浏览
    L

    屏幕截图 2021-06-05 161820.png openfoam自带案例interPhaseChangeFoam求解器下的cavitatingBullet,该图为我跑完0.05s后的水相分布,和0s时没啥变化,请问我哪里错了吗?

  • 边界条件设置

    3
    3 帖子
    2k 浏览
    S

    @李东岳 感谢李老师回复,我先试下