Skip to content
  • 最新
  • 版块
  • 东岳流体
  • 随机看[请狂点我]
皮肤
  • Light
  • Cerulean
  • Cosmo
  • Flatly
  • Journal
  • Litera
  • Lumen
  • Lux
  • Materia
  • Minty
  • Morph
  • Pulse
  • Sandstone
  • Simplex
  • Sketchy
  • Spacelab
  • United
  • Yeti
  • Zephyr
  • Dark
  • Cyborg
  • Darkly
  • Quartz
  • Slate
  • Solar
  • Superhero
  • Vapor

  • 默认(不使用皮肤)
  • 不使用皮肤
折叠
CFD中文网

CFD中文网

火

火山口玩泥巴

@火山口玩泥巴
关于
帖子
53
主题
11
群组
0
粉丝
2
关注
0

帖子

最新

  • 大规模算例paraview看结果的一种方法
    火 火山口玩泥巴

    @xpqiu 好的谢谢:mianmo:


  • 大规模算例paraview看结果的一种方法
    火 火山口玩泥巴

    @火山口玩泥巴 因为准备申请个学校的云主机进行后处理,计算时用了128核,现在申请云主机的核数也需要至少128吗?


  • 大规模算例paraview看结果的一种方法
    火 火山口玩泥巴

    请问使用这个功能的时候,需要保证cpu核数大于等于并行时decompose的数吗?


  • 多相流计算中入口的速度过小导致计算结果不符
    火 火山口玩泥巴

    想尝试与一篇文献中实验结果作对比,但根据文中所提供的条件算出的入口速度很小(0.18m/s),以此作为入口的边界条件;同时入口边界的alpha.water也固定为1,计算一段时间后如图所示,发现入口处几乎没有水相的存在了。这是为什么?

    ae0a7d49-651d-4636-a27b-c8a0ff4b9334-image.png

    入口处的边界条件如下:

    U
    {
            type            fixedValue;
            value           uniform (0.18 0 0);
    }
    
    p_rgh
    {
            type            zeroGradient;
    }
    
    alpha.water
    {
            type            fixedValue;
            value           uniform 1;
    }
    
    alpha.air
    {
            type            zeroGradient;
    }
    

  • 三相空化计算明渠流动失败
    火 火山口玩泥巴

    @李东岳 LES和laminar都试过,只要把空化核数降一些就不会发散。


  • 三相空化计算明渠流动失败
    火 火山口玩泥巴

    @李东岳 是的老师,OF没提供三相空化的标准求解器,只能自己改着用,但又没啥经验,总怕自己改的这玩意稳定性太差。关键我对湍流模型了解的太浅了,k和omega越界这个东西我想改都无从下手。


  • 三相空化计算明渠流动失败
    火 火山口玩泥巴

    @李东岳 不是完全自己写的老师,是在interPhaseChangeFoam基础上改的,加了一个空气相进去。
    在粗网格下面,用interFoam来算,使用komegaSST,在一定的k和omega的边界条件下不会出现发散的情况;但用自己的求解器去算即使把相变系数调到0,k和omega的边条初条怎么调都还是会发散。


  • 三相空化计算明渠流动失败
    火 火山口玩泥巴

    直接先用upwind处理相方程了。然后尝试使用LES的话可以正常计算粗网格,前提是得把空化核密度降低两个数量级,否则出口处会出现奇怪的低压区,这里不懂为什么会这样。:134:


  • 三相空化计算明渠流动失败
    火 火山口玩泥巴

    同时我发现相分数存在越界情况,尤其是在细的那一套网格中,越界情况尤其明显,蒸气相甚至最小值达到了-1.4:136: ,相方程的求解使用的是显式MULES(图方便),这是因为显式MULES效果不佳嘛??我想在每次求解过后强行限定α的值:

    alpha1 = min(max(alpha1, scalar(0)), scalar(1));
    

    这样是否可行?


  • 摩擦雷诺数特征长度取值?
    火 火山口玩泥巴

    没听说过摩擦雷诺数啊,如果只是单纯的雷诺数的话特征长度管道流一般就取:截面面积/截面周长即水力直径吧。


  • 三相空化计算明渠流动失败
    火 火山口玩泥巴

    @李东岳 都试过了老师,都是会差不多的时间发散,即出现bounding k和bounding omega;
    刚刚我用粗网格的时候又出了一些别的问题:当空化模型里面的n0设的较大就会出现解压力方程的时候浮点溢出,具体是出现在GAMGSolver的Scale函数里面的缩放因子的分母的计算那里,我折腾半天没弄明白那里的溢出原因(我看了溢出单元的压力方程源项以及空化模型提供的主对角线项,都不怎么会导致溢出),因为我暂时还没弄懂那个函数的作用,就先给压力方程换成PCG求解了,这个问题暂时就没出现了;


  • 三相空化计算明渠流动失败
    火 火山口玩泥巴

    我现在尝试换了一套几何模型,重新画了更细的网格,网格质量除了几个面会有一点非正交性较大外其他没有问题,仍然会出现一样的问题,即在进行湍流模型参数更新的过程中浮点溢出。
    大佬们有啥思路解决嘛?是边界条件的问题吗?


  • 三相空化计算明渠流动失败
    火 火山口玩泥巴

    大佬们给点建议:mihu:


  • 三相空化计算明渠流动失败
    火 火山口玩泥巴

    顶一下


  • 三相空化计算明渠流动失败
    火 火山口玩泥巴

    在interPhaseChangeFoam的基础上改了改弄了个三相空化求解器来算明渠流动,第一步就会出现如下所示错误c4099a17-c7bf-4743-8bdf-530fb10e13c9-image.png
    Debug发现错误出现在湍流模型参数的更新这里,用的RANS模型:
    a23fc437-5049-4573-9d94-ffb245a7564f-image.png
    具体报错是FieldFunctions.C里的sumMag函数这里:
    21abe3c5-cca7-48ae-ab3f-588b535a35af-image.png
    请问各位大佬,是啥原因造成的呢?


  • 关于interPhaseChangeFoam里面的一些疑问
    火 火山口玩泥巴

    @火山口玩泥巴 似乎化简后变成了满足连续性方程限制的非守恒动量方程形式。为什么要这样建立方程呢?:141:


  • 关于interPhaseChangeFoam里面的一些疑问
    火 火山口玩泥巴

    借一下这个帖子,问一下interPhaseChangeFoam中,UEqn里的- fvm::Sp(fvc::ddt(rho) + fvc::div(rhoPhi), U)这一项是干啥用的嘞?似乎是连续性方程的左边和速度相乘?实际意义是啥呢?


  • 含沙水流计算问题
    火 火山口玩泥巴

    @Joe 我也是小白还在折腾这个三相的求解器😂,OF似乎没有现成的三相空化求解器,所以想着只能在两相模型的基础上推广到三相然后自己改一个求解器出来了。


  • 含沙水流计算问题
    火 火山口玩泥巴

    @李东岳 谢谢老师,看到了而且觉得效果可以的了,就准备试试这个方法了:146:

  • 登录

  • 登录或注册以进行搜索。
  • 第一个帖子
    最后一个帖子
0
  • 最新
  • 版块
  • 东岳流体
  • 随机看[请狂点我]