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中文网

C

Caijinjin

@Caijinjin
关于
帖子
49
主题
12
群组
0
粉丝
1
关注
0

帖子

最新

  • 多孔介质区域速度突变的现象
    C Caijinjin

    @han1112 好的,感谢大佬。按照这个意思来说的话,那么真实条件下的填料扰流是没办法用多孔介质源项来实现了。即使考虑了空隙率,速度也是不真实的。


  • 添加多孔介质源项对管式填充床中的液-液两相流模拟想到的一些问题
    C Caijinjin

    最近在用multiphaseEulerFoam模拟管式填充床中的液-液两相流。填充状态如下:
    95b7f9cd-a542-47e2-b536-70b2ac5d2ae6-image.png
    真实情况下管式填充床中呈现两相液滴流。因为填充床很长,所以打算利用多孔介质源项来模拟其中的流动,但是现在想到了一些问题还在探索:

    1. 使用多孔介质模型后,模拟使用的几何文件中既然没有填料结构,那么相同边界条件下,流型也会从真实情况下的湍流扰流变成层流(两相分层流动),这种情况是否有方法可以进行改变?
    2. openfoam中在多孔介质区域求解ns方程使用的表观流速,那么真实流速是在后处理的时候进行计算,还是我可以直接改变条件设置使得求解使用真实流速?

  • 多孔介质区域速度突变的现象
    C Caijinjin

    最近我也打算在openfoam中开始做多孔介质模拟和算法学习,目前是打算利用多孔介质模拟管式填充床中的液-液两相流,但是目前还有一些问题:1. 使用多孔介质源项之后,模拟的几何文件不是没有填料的三维空间占据了嘛,那么流动还会是湍流或者说还有填料的扰流作用嘛?2. 我之前使用fluent做过多孔介质模拟,可以选择物理速度和表观速度,那么openfoam中好像新版本的没有空隙率设置的内容了,所以of用的表观速度?3. 没有湍流了那液-液两相流是不是会形成分层流动?


  • 多孔介质多相流出口回流问题
    C Caijinjin

    请问您这个多孔介质源项是如何设置的呢?本人之前用fluent运算的时候,模拟条件实质是没有多孔介质的三维空间占据的,还想问一下您这个案例多孔介质区域如何设置的?


  • 球形颗粒接触点的处理方法求助
    C Caijinjin

    各位好,我最近用OpenFOAM也在做这方面的模拟,现在使用SHM划分网格遇到了一些问题,主要是网格质量、非正交性、凹度和负体积网格方面的,能否交流一下


  • 群体平衡方程中新增子液滴分布函数,关于韦伯数的程序编译问题
    C Caijinjin

    @李东岳
    李老师,这个函数里面,d,也就是xk,在这段程序里面是一个未知的自变量,只有在最后计算子液滴分布的时候才给定,我定义的分配函数里面也有这个定义:
    b10fda7f-6837-4b94-8d8b-125646f77767-image.png
    这里定义成dimensionedScalar就没问题,所以可能不是d的问题了


  • 群体平衡方程中新增子液滴分布函数,关于韦伯数的程序编译问题
    C Caijinjin

    李老师好,最终计算子液滴分布的结果是一个标量,那么在计算calcWe时返回场的话还能对应的使用吗?还有韦伯数是一个无量纲的数,如果是场的话这里应该选择是那种类型的场呢?


  • 群体平衡方程中新增子液滴分布函数,关于韦伯数的程序编译问题
    C Caijinjin

    各位大佬好,目前遇到一个问题。
    在群体平衡方程中新增子液滴分布函数,程序中需要用到韦伯数,韦伯数的计算公式是:
    6b9c9c38-8048-4cdb-8147-0b14bb860122-image.png
    公式中d是母液滴的直径,把它写到子液滴分布程序中换算成了关于母液滴体积的函数,如下:
    db288f95-d661-4e36-aea8-838eb0958da3-image.png
    遇到的问题是编译过程中出现下面的问题:
    00d17996-f245-4907-9b04-bf86a583c0d7-image.png
    请问有没有大佬了解这方面的问题,这是量纲或者场没对应上吗?


  • 集群上并行测试OpenFOAM,并行效率并没有比单节点提升
    C Caijinjin

    @李东岳 老师,是这样的,这种现象在单个节点跑案例的时候是没有的,CPU是接近满负荷在运算的,但是一到多个节点并行就出现。
    现在我们换了IB交换机之后多节点还没跑起来过。。。


  • 集群上并行测试OpenFOAM,并行效率并没有比单节点提升
    C Caijinjin

    @李东岳 李老师,我们组集群是新买的,应该不太可能是被黑了,除非是大厂他们给加了限制。现在我导给换了100G的IB交换机,就是运行案例的时候,它虽然显示是在run的,但是log文件里面没有实际的计算结果。我看了调度系统的日志说是:
    Open MPI accepted a TCP connection from what appears to be an another Open MPI process but cannot find a corresponding process entry for that peer.
    暂时还在寻找原因。


  • 气液搅拌釜顶部出口速度异常,导致epsilon等结果发散
    C Caijinjin

    @李东岳 好的,我先尝试着修改下再去跑跑看


  • 气液搅拌釜顶部出口速度异常,导致epsilon等结果发散
    C Caijinjin

    @李东岳李老师,最初使用过snappyHexMesh划分过结构化网格,模拟出错之后就改用icem非结构化网格来画稍微粗糙一点的网格去调试了


  • 气液搅拌釜顶部出口速度异常,导致epsilon等结果发散
    C Caijinjin

    搅拌釜中气液两相模拟,釜中搅拌桨底部设置了两个圆形气体入口,气速最大0.3 m/s,釜内为水,顶部出口。我用reactingTwoPhaseEulerFoam模拟的时候最终都会出现下图中的结果,顶部出口有个位置的速度非常大,超过10 m/s,最终导致epsilon发散,气泡粒径出现负值。53686cdd-aa5e-49a3-8066-b0bd42254420-image.png
    请问一下各位有没有遇到过类似的情况,出口的边界条件类型,离散格式都修改过还是会出现一样的情况,而且搅拌釜顶部流出一段空气之后,这块区域还会变大。


  • 集群上并行测试OpenFOAM,并行效率并没有比单节点提升
    C Caijinjin

    @李东岳 在 集群上并行测试OpenFOAM,并行效率并没有比单节点提升 中说:

    -genv FI_PROVIDER tcp你这个去掉不能跑的话。你如何确定走的是infiniband,而不是以太网模式。我们这面跑openfoam不需要这个参数。我们之前测试也出现过你这种情况。后来我们换交换机硬件了。但目前我还不确定现在我们这5节点能到什么样,得下周能出个测试结果。另外,openfoam离散设置差异(比如GAMG那个),我个人感觉不会引起特别大的差异。不过你可以实测看看,我也不100%确定

    @xpqiu 这位大佬之前好像测试过2048个核心 :mianmo:

    老师,想问一下,用openmpi在自己组里面的集群上用pbs跨节点并行(10g以太网交换机)命令用的mpirun --mca btl_tcp_if_include <ip地址> -np reactingTwoPhaseEulerFoam -parallel,运行的时候发现计算的节点上cpu的用户进程占比us只有50-60%,系统进程占比sy有40-50%,这个问题有没有什么好的解决办法?


  • 跨节点运行openfoam慢
    C Caijinjin

    你好,想问你一下,这个跨节点运算速度慢的问题有解决吗?我这边自己的集群也有这个问题,目前还在找原因中,集群centos7,openfoam运行前会有环境变量设置的脚本(gcc编译,openmpi并行),openfoam和openmpi等软件都安装在共享文件夹中。我这边的测试案例结果如下:
    08374518-59f3-4cde-9e4c-a771c6af1c9a-image.png


  • 管式填充床中两相流的模拟中,循环边界压力梯度、分散相含率运算出错的问题
    C Caijinjin

    @李东岳 好的,李老师,我去调试看看


  • 管式填充床中两相流的模拟中,循环边界压力梯度、分散相含率运算出错的问题
    C Caijinjin

    @李东岳
    单相流我还没尝试过,之前用过空的圆管做过一样的案例,问题都解决了。
    现在我这边入口fixedValue,mapped边界条件都能是正常运行,就cyclicAMI出现问题,所以我打算先看看mapped的结果。


  • 管式填充床中两相流的模拟中,循环边界压力梯度、分散相含率运算出错的问题
    C Caijinjin

    最近在做管式填充床中两相流动的模拟,采用循环边界cyclicAMI达到缩短管长,减少网格的目的。网格使用snappyHexMesh画完后没有问题,动力源项给的是连续项平均速度Ubar,为0.03m/s,填充床几何为
    07456e007b85fa4ee9f49c7e8a053a4.png
    运算时出现如下结果:
    bccffb96-29ae-41c1-9a37-5dd7b44ff9c0-image.png
    主要就是动力源计算的压力梯度一直在正负之间反复跳动,且数值非常大,压力P_rgh的迭代次数会超限。运算一段时间后会直接导致求解的相分数远远超出0-1的范围,如下图所示:
    4a88b6f4-03d2-4327-9c79-1a8102734d26-1666706031077.png
    想请教一下对这种循环边界熟悉的老师同学,能够解答一下我这些问题。
    另外附上我的fvSolution:

    solvers
    {
        "alpha.*"
        {
            nAlphaCorr      1;
            nAlphaSubCycles 3;
            implicitPhasePressure yes;
           
            solver          smoothSolver;
            smoother        symGaussSeidel;
            tolerance       1e-8;
            relTol          0;
            minIter         1;
        }
    
        p_rgh
        {
            solver          PCG;
            preconditioner  FDIC;
            tolerance       1e-6;
            relTol          0.01;
        }
    
        p_rghFinal
        {
            $p_rgh;
            tolerance       1e-6;
            relTol          0;
        }
        
        "U."
        {
            solver          smoothSolver;
            smoother        symGaussSeidel;
            tolerance       1e-7;
            relTol          0;
            minIter         1;
        }
    
        "(e|h).*"
        {
            solver          smoothSolver;
            smoother        symGaussSeidel;
            tolerance       1;
            relTol          0;
            minIter         0;
            maxIter         0;
        }
    
        "(k|epsilon|omega).*"
        {
            solver          smoothSolver;
            smoother        symGaussSeidel;
            tolerance       1e-7;
            relTol          0;
            minIter         1;
        }
    }
    
    PIMPLE
    {
        nOuterCorrectors    3;
        nCorrectors         2;
        nNonOrthogonalCorrectors 1;
        pRefPoint (0.10399 0 0);
        pRefValue 101325;
    }
    
    relaxationFactors
    {
        equations
        {
            ".*"            0.7;
        }
    }
    

    fvSchemes:

    ddtSchemes
    {
        default         Euler;
    }
    
    gradSchemes
    {
        default         Gauss linear;
        limited         cellLimited Gauss linear 1;
    }
    
    divSchemes
    {
        default                         none;
    
        "div\(phi,alpha.*\)"            Gauss vanLeer;
        "div\(phir,alpha.*\)"           Gauss vanLeer;
    
        "div\(alphaRhoPhi.*,U.*\)"      Gauss linearUpwind limited;
        "div\(phi.*,U.*\)"              Gauss linearUpwind limited;
    
        "div\(alphaRhoPhi.*,(h|e).*\)"        Gauss linearUpwind limited;
        "div\(alphaPhi.*,f.*\)"         Gauss linearUpwind limited;
        "div\(alphaRhoPhi.*,K.*\)"      Gauss linearUpwind limited;
        "div\(alphaPhi.*,p\)"           Gauss linearUpwind limited;
        "div\(alphaRhoPhi.*,\(p\|thermo:rho.*\)\)"     Gauss linearUpwind limited;
        "div\(alphaRhoPhi.*,(k|epsilon|omega).*\)"      Gauss linearUpwind limited;
    
        "div\(\(\(\(alpha.*\*thermo:rho.*\)\*nuEff.*\)\*dev2\(T\(grad\(U.*\)\)\)\)\)"  Gauss linear;
    }
    
    laplacianSchemes
    {
        default         Gauss linear corrected;
    }
    
    interpolationSchemes
    {
        default         linear;
    }
    
    snGradSchemes
    {
        default         corrected;
    }
    
    fluxRequired
    {
        default         no;
    }
    
    wallDist
    {
        method          meshWave;
        nRequired       yes;
    }
    

  • 计算出现负相率和负湍动能
    C Caijinjin

    @cyberk 抱歉,这个问题我没有遇到过,但我看了你之前用ICEM画的网格,这个六叶Ruston桨的搅拌槽的网格不太复杂的呀,这个你可以试试只画1/6个区域,然后对称过去,画完整体之后再去openfoam中用topoSet划分搅拌的区域


  • 计算出现负相率和负湍动能
    C Caijinjin

    @cyberk 我这边是因为缺少了createBaffles字典,如果案例的几何边界中存在零厚度面就需要对其进行设置,我这边的搅拌釜的搅拌桨就是零厚度面

  • 登录

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