Skip to content

OpenFOAM

OpenFOAM交流区

5.0k 主题 29.3k 帖子
  • 合集:OpenFOAM安装准备工作贴

    已固定
    24
    24 帖子
    25k 浏览
    S

    @李东岳 好的,谢谢老师!

  • OpenFOAM教程 初级入门建议

    已固定
    55
    55 帖子
    133k 浏览
    李东岳

    新版的transportProperties改成了physicalProperties,你改个名

  • 有关DriftFluxFoam的相关设置

    30
    30 帖子
    625 浏览
    李东岳

    这个就是模型本身的缺陷了。目前无法处理。你要是写文章的话,我有个思路,我可以指导你来写。但是需要挂我的名字,位置随意,同时需要你导师同意。

  • OpenFOAM使用mpirun的时候报错

    8
    8 帖子
    870 浏览
    李东岳

    @zzkluck :146: :146: :146:

  • rhoCentralFoam求解shockTube算例出现数值震荡

    13
    13 帖子
    494 浏览

    @李东岳 老师我OF-4.1版本算例constant文件夹下得有turbulenceProperties字典,里面写的是laminar,不涉及湍流。好像OF-8以后都不需要turbulenceProperties字典了。

  • SHM为什么生成两个90度垂直面的边界层很烂?

    8
    8 帖子
    347 浏览
    zzkluckZ

    按照文档里的说法,有时换一下meshShrinker也许会有效果,感觉文档里那个例子和你的也有点像。

    addLayersControls { ... meshShrinker displacementMotionSolver; solver displacementLaplacian; displacementLaplacianCoeffs { diffusivity quadratic inverseDistance 1(wall); } ... }

    相应的还需要调整下fvSolution和fvScheme,具体文档里都有。

    话说新版文档虽然变好看了,但图居然是糊的,可能还需要参考下旧版文档。

    另外这里也请教下各位老师,按我的理解snappyHexMesh做完snap之后,在layer这步应该是从几何表面出发,向挤出来一块空间把边界层塞进去。怎么在楼主的这个例子中反而是向内把几何都挤变形了

  • precice tourorials中的受力问题

    6
    6 帖子
    113 浏览
    R

    @Mina-Lee 在 precice tourorials中的受力问题 中说:

    因为你编译的openfoam-adapter里,已经对流固耦合交界面的force变量引用了,OF在运行完成后,adapter会引用force变量,precice运行完毕后,会覆盖这个force变量。
    precice就是给openfoam施加了一个狄利克雷位移边界条件-理解了这句话就全明白了。

    其他的能理解,但是这个“precice运行完毕后,会覆盖这个force变量”指的是
    在controldict中定义了“forces”那么openfoam输出的forec是重新计算的,与流固耦合过程中交换的不是一个值的意思吗

  • 7 帖子
    69 浏览
    李东岳

    调试的时候用ascii,跑算例用binary。如果网格几十万以上,binary会小非常多

  • 15 帖子
    2k 浏览
    M

    你好

    PIMPLE: iteration 40 **time step continuity errors : sum local = 20238.39271, global = -1326.825454, cumulative = -1286.292262** **time step continuity errors : sum local = 1965.022984, global = -493.7266286, cumulative = -1780.018891** PIMPLE: iteration 22 GAMG: Solving for p, Initial residual = 0.1839153372, Final residual = 1.265686901e-06, No Iterations 20 GAMG: Solving for p, Initial residual = 0.1398796637, Final residual = 1.397575026e-06, No Iterations 17 GAMG: Solving for p, Initial residual = 0.1463017811, Final residual = 1.437462198e-06, No Iterations 16 GAMG: Solving for p, Initial residual = 0.121320625, Final residual = 6.97291136e-08, No Iterations 20 PIMPLE: iteration 23 .GAMG: Solving for p, Initial residual = 0.1600126604, Final residual = 1.256444652e-06, No Iterations 23 GAMG: Solving for p, Initial residual = 0.09022545056, Final residual = 7.317342264e-07, No Iterations 20 GAMG: Solving for p, Initial residual = 0.08084307981, Final residual = 7.297894683e-07, No Iterations 18 GAMG: Solving for p, Initial residual = 0.06296048891, Final residual = 7.500379861e-08, No Iterations 22

    你的问题在于:
    1、瞬态项上,你的瞬态项很难收敛,残差爆掉了,这个似乎与你的瞬态项离散方式有关,可以试一下二阶向后差分。

    ddtSchemes { default CrankNicolson 0.9; }

    2、你的泊松方程收敛性很差,外循环很难收敛,残差很高,你定义的外循环收敛条件根本达不到,泊松方程在非正交修正后,迭代次数反而变多,证明你的网格质量可能较差。
    建议:1.更换瞬态项离散方式。2.增加非正交修正次数,提高你的网格质量。
    注意:外循环定义50次合理,一般难收敛问题,最多10次就可以收敛了。流固耦合交界面前几个时间步不收敛很正常,但是你这个是流体不收敛,这就不对了,建议你先单独跑一下流体看看,你这个流体设置不太好。
    祝好~

  • 1 帖子
    22 浏览
    W

    最近想做空化的流固耦合问题,发现计算经常发散,表现为alpha.water大于1,然后逐渐发散掉,但是库朗数是正常的,大约为1.1。想问问大家有没有遇到过类似问题,有什么解决办法?

  • Openfoam计算泄露的问题

    8
    8 帖子
    2k 浏览
    1

    @李东岳 好的博主 但是为什么计算后在paraview里的质量分数大于1呢

  • 关于div的小问题

    4
    4 帖子
    50 浏览
    李东岳

    你这个方程,直接写成div(phi,e)

  • 6 帖子
    66 浏览
    U

    @李东岳 好的,谢谢老师!

  • 10 帖子
    356 浏览
    李东岳

    我想确认一下,7楼的代码真的在2206上运行成功了吗?谢谢

    可以,不过我没用的interfoam,我用的什么求解器我忘了。可能是simplefoam、icoFoam之类。你试试这俩能不能运行?

  • 储热模拟相关问题

    11
    11 帖子
    366 浏览
    Y

    @李东岳 东子老师,我把可压+稳态算法仔仔细细的看完了,大体上看懂了,但有一些疑问向您请教。您上次说我排版不利于阅读,这次我专门看了您提供的帮助里面的教程。

    问题1

    您在文章里说:上标n表示为当前迭代步(已知),上标表示下一个迭代步(待求),n+1也表示下一个迭代步。在稳态条件下,5b2efa88-9199-4b80-b4f0-c0c886fa0089-CodeCogsEqn.svg 应该等于658840c8-71e8-4642-9a77-9260eaa6b6f6-CodeCogsEqn (1).svg 吧。方程(18)-(14)等号左边不应该是0吗?为啥是d196288e-c3a2-4ce2-8580-b309954568ec-CodeCogsEqn (2).svg 呢?上标、**、n与n+1是啥关系呀?我猜测*、**是不是n到n+1之间的迭代步?

    问题2

    方程(23)没有c9eb4ec2-5110-41b3-abda-7fb63743d272-CodeCogsEqn (3).svg 呀?另外,为啥是15d26edd-0fbf-4db7-8143-5855742ca1e3-CodeCogsEqn (4).svg 而不是8e7d2219-1a52-48f2-8bb4-6a26b87935c8-CodeCogsEqn (5).svg 或者e88aae6e-7526-4d34-adc1-c5e70ac64eb3-CodeCogsEqn (6).svg 呢?
    e1cfe02f-da35-4150-9e43-c2a4356a47b5-image.png

    问题3

    由方程(3)可得670c6b75-18cf-4c6d-a39e-0a4f5d9bb832-CodeCogsEqn (7).svg ,又7854ecd4-8f21-40f4-8fd2-c5b73b97ff6b-CodeCogsEqn (8).svg ,所以 7fb1a991-7f26-4b14-a90e-6debc0afa22c-CodeCogsEqn (9).svg
    但由方程(27)得到的却是图中红框部分,如何得出方程(28)中蓝框部分的呢?方程(36)是我算的结果,方程(42)又是红框的结果了。
    9fc6e84b-b77d-45a3-af6b-94b05b4b0a23-image.png

    问题4

    方程(19)由(18)-(14)得来,1edcea56-096d-4193-84f5-b196f809f1b8-CodeCogsEqn (10).svg ,而图中这个推导过程:
    5cf4e3dd-7af1-4218-96d6-0d761db6ec75-image.png是 48040ba7-9ec7-4fcb-b8db-caad5cf8c058-CodeCogsEqn (11).svg ,是不是矛盾了呢?

    同样的,方程(43)中,0ef0f602-2a76-4be8-8342-0f7b005947b9-CodeCogsEqn (12).svg ;
    方程(25)中,dcb61c86-7441-4b63-8724-e1b4189113e2-CodeCogsEqn (13).svg
    我被这些上标n、'、*搞懵了:135:

    问题5

    方程(48)中,红框部分是不是少了一个c4cc7b73-f09a-4056-b44f-16484a4ee9c5-CodeCogsEqn (14).svg ?
    72a9548f-2d79-4c08-8db7-d50ec6ff9dbb-image.png

  • open foam10植入曳力模型

    3
    3 帖子
    47 浏览
    WeakFormW

    有这个文件
    image.png

  • compressibleInterFoam温度场没变化

    1
    1 帖子
    22 浏览
    1

    各位大佬,我将一部分的气体设置成了750K,作为喷管内部气体的温度,喷出去之后接触也为气体,但是计算之后温度场没发生变化。我采用的是compressibleInterFoam76b063eb-0cee-457b-b789-af3aaad2ef33-1732091843128.png
    以上图片是计算之后的结果,红色部分为我赋值的750K,蓝色为最初的初场300K,有大佬知道是什么原因吗

  • 求一个固液传热的算例

    2
    2 帖子
    2k 浏览
    1

    您好,你的问题解决了吗,我也遇到了这个问题,温度场不发生变化

  • 并行计算报错

    2
    2 帖子
    37 浏览
    李东岳

    我没有openfoam11,我测试了一下openfoam11,没有问题,不是bug。是不是你的机器只有一个核

  • 9 帖子
    163 浏览
    学流体的小明

    @yangom 无滑移边界条件,下面就是固壁面,无滑移就是壁面上速度为0,那下表面就需要从入口处开始发展发展出一个边界层。在这个边界层中,沿壁面法向向外,速度当然是从0逐渐变化成外流速度的
    你直接看也看不出来有没有速度,画几条曲线出来。
    1 靠近近壁面某个高度上,沿流动方向的一条线,画这条线上沿流动方向的速度
    2 选几个位置,沿垂直于壁面方向的一条线,画这条线上沿流动方向的速度