Skip to content
  • 有关DriftFluxFoam的相关设置

    OpenFOAM
    29
    0 赞同
    29 帖子
    623 浏览

    @李东岳 李老师您好,我重新绘制了一下nut-y的关系,老师您看看,差异很大,我附上了原始数据屏幕截图 2024-11-21 140739.png 。不算时间平均的值,我就对速度求了一下时均值。屏幕截图 2024-11-21 140018.png

  • 0 赞同
    13 帖子
    493 浏览

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

  • 0 赞同
    8 帖子
    345 浏览
    zzkluckZ

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

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

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

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

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

  • precice tourorials中的受力问题

    OpenFOAM
    6
    0 赞同
    6 帖子
    111 浏览
    R

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

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

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

  • scalar transport如何考虑沉降

    Algorithm
    6
    0 赞同
    6 帖子
    4k 浏览
    wzhmiasanmiaW

    @李东岳
    请问老师,我最近刚开始用openfoam,主要是模拟车内CO2和PM2.5的分布然后进行空气质量分析,我希望加入这种沉降的作用以模拟出CO2比空气重然后在车舱内分层这样一种现象。想问问有没有具体的实现方式或者相关的研究的文章,谢谢!
    另外目前我想到的思路是比如在用simplefoam计算稳态速度场的时候加入源项考虑CO2的重力和浮力,得到稳态的速度场后在scalartransportfoam中作为初始速度场计算,但是用简单的case跑了一下感觉这样做对于标量传输基本就没什么影响,问题是出在哪里呢?

  • 0 赞同
    7 帖子
    69 浏览
    李东岳

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

  • 0 赞同
    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次就可以收敛了。流固耦合交界面前几个时间步不收敛很正常,但是你这个是流体不收敛,这就不对了,建议你先单独跑一下流体看看,你这个流体设置不太好。
    祝好~

  • 0 赞同
    1 帖子
    27 浏览

    我要用对流扩散方程模拟不同水动力下物质的扩散,因此我需要得到不同水动力条件下的物质的扩散系数D。目前有两种方案来研究:
    ① 通过实验来得到速度场与浓度场的变化来率定,得到了速度与浓度变化随时间t的变化数据,根据对流护散方程,反推出D。如果我通过实验得到了多个扩散系数,在模拟的时候要怎么带入到对流扩散方程中呢?是要平均吗?万一每个扩散系数量级差别很大 平均之后的D也不能很准确的有普适性。
    ②第二种方案是大多数文献采用的方法,用湍流涡粘性比上普朗特-施密特数:D=vt/σ
    ​现在已知我实验设备只能测出湍流动能与平均速度场。第二种方案不知道怎么进展下去。
    不知道这两种方案行不行的通,目前网上对于如何率定扩散系数并没有详细的说明。

  • 0 赞同
    1 帖子
    22 浏览
    W

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

  • Openfoam计算泄露的问题

    OpenFOAM
    8
    0 赞同
    8 帖子
    2k 浏览
    1

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

  • 关于div的小问题

    OpenFOAM
    4
    0 赞同
    4 帖子
    50 浏览
    李东岳

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

  • 0 赞同
    6 帖子
    66 浏览
    U

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

  • 0 赞同
    10 帖子
    356 浏览
    李东岳

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

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

  • 储热模拟相关问题

    OpenFOAM
    11
    0 赞同
    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植入曳力模型

    OpenFOAM
    3
    0 赞同
    3 帖子
    47 浏览
    WeakFormW

    有这个文件
    image.png

  • fluent Error: Surface creation failed!报错

    Fluent
    1
    0 赞同
    1 帖子
    26 浏览
    2

    第一个监测点没有超出域外,为什么创建面失败呢? 感谢各位解答疑惑。
    0406f157-633e-4582-bcb2-b9b35fb7470d-image.png

  • compressibleInterFoam温度场没变化

    OpenFOAM
    1
    0 赞同
    1 帖子
    19 浏览
    1

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

  • 求一个固液传热的算例

    OpenFOAM
    2
    0 赞同
    2 帖子
    2k 浏览
    1

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

  • CFD中文网 定位与规范

    已固定 公告
    101
    0 赞同
    101 帖子
    230k 浏览
    H

    我写 固我在,嘿嘿

  • 并行计算报错

    OpenFOAM
    2
    0 赞同
    2 帖子
    36 浏览
    李东岳

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