Skip to content
  • 关于不可压缩求解器流体密度的设置

    OpenFOAM
    5
    0 赞同
    5 帖子
    8k 浏览
    OYYOO

    @李东岳 在 关于不可压缩求解器流体密度的设置 中说:

    对于不可压缩流动,密度/压力的大小是无关紧要的。重要的是压力差。因此你不需要保证求解的密度是你要的密度。

    \begin{equation}
    \frac{\partial \rho \mathbf{U} }{\partial t}+\nabla \cdot \left(\rho \mathbf{U} \mathbf{U} \right)-\nabla \cdot \left(\mu \nabla \mathbf{U} \right)=-\nabla p
    \label{NS2}
    \end{equation}

    比如这个方程,如果密度不变,你直接除掉密度就可以。因此在不可压缩求解器中你不需要设置密度的大小

    李老师,但是如果计算壁面剪切应力wall shear stress的话,$\mu \frac{dU}{dy}$,式中 $\mu = \nu * \rho$的是与密度有关的。这时候算出的Wall shear stress是不是与真实值存在差距呀,不知道有啥办法调整后处理输出的wallshearstress没,不然只能在paraView中自己乘以真实密度了。

  • 圆柱相切的划分策略

    Meshy
    4
    0 赞同
    4 帖子
    7k 浏览
    I

    谢谢老师~~

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

    厉害厉害,这个图非常好

  • 一个积分

    CFD彩虹条
    5
    0 赞同
    5 帖子
    8k 浏览
    C

    find . -type f | xargs -n 5 touch

  • 0 赞同
    2 帖子
    4k 浏览

    @cfd-china 直接用WIN自带的远程连接就可以了

  • 有限面方法 Finite area method

    Algorithm
    2
    0 赞同
    2 帖子
    3k 浏览
    W

    以前简单看过,好像是研究表面活性剂的方法。大概是在三维曲面上用类似FVM的方法离散控制方程并求解。

  • 到底静压值是怎么得来的?

    Algorithm
    1
    0 赞同
    1 帖子
    3k 浏览
    L

    压差是分析很多工程设备的关键参数,在之前,主要都是针对多孔介质的Ergun模式来处理这一块。
    也跟试验对比很好。但是,最近新开发的一个模型,虽然与试验的结果对比比较好,但感觉总是有问题。
    于是不知道CFD中静压值到底是怎么得出来的?
    对于一块多孔介质好像比较好处理,但是如果这个时候多孔介质表面沉积了灰,这个时候这块的整个压降会升高,然后灰越堆越多整个表面的压降都会升高。然而,这个局部压降和表面灰的分布十分相关,造成整个表面的这个灰分的压降相差一千多Pa。而最后一般的检测压降点,都是压差变松器的两个接口,一般都是设置在容器外壁的。但是如果我模拟检测点和试验一样的时候,这个时候的压差比我的表面计算压差平均值要小得多。(压差是1200多Pa,但是我的多孔介质有800多Pa,加上平均粉尘层的压降有640多Pa,合起来应该有1400多Pa,而不是1200多Pa.)
    这就产生了一个问题,这个通过N-S方程求解的压力项到底是怎么处理的?
    当他通过一层多孔介质的时候,多孔介质内部有压力源项,而在外部没有,那么当他气体从多孔介质外向内流通的时候,他的这层压力和速度是怎么计算的?如果这个时候我再添加一个源项在多孔介质层,而这个源项在多孔介质外是没有的,那么这层额外的源项是不是也会通过增加流体的静压来表达呢?(通过伯努利方程似乎很好理解,但是N-S是怎么实现的,他的这个\Delta P 是怎么获得的?)

  • 过饱和度计算

    Fluent
    1
    0 赞同
    1 帖子
    1k 浏览
    V

    fluent自定义过饱和函数时应该用那个压力数据来计算?我的静压只有不到100pa,太小了

  • refineMesh加密网格出现错误

    OpenFOAM
    1
    0 赞同
    1 帖子
    945 浏览

    我尝试用refineMesh进行多层加密,第一层加密完成后,在其内部的一个子区域进行第二次加密时,出现useHexTopology specified but cell 5018638 on face 7394 of patch bottom is not a hex的问题,我检查发现5018638单元并不在我的第二次加密区域内,为什么仍然要检查他是否为hex?
    以下为refineMeshDict:

    set expand7; coordinateSystem patchLocal; patchLocalCoeffs { patch bottom; tan1 (0 0 1); } directions ( tan1 tan2 ); useHexTopology yes; geometricCut no; writeMesh no;

    如果改用geometricCut会不会导致网格质量下降?

  • guns and roses复合

    CFD彩虹条
    1
    0 赞同
    1 帖子
    3k 浏览
    李东岳

    也不能说是复合,2016年Slash 51岁生日,Axl和Slash同台表演,时隔20年的分道扬镳。
    http://www.gunsnroses.com/tour
    0_1485104660529_20057521-mmmain.jpg

  • 0 赞同
    9 帖子
    14k 浏览
    李东岳

    充值100送100

  • fluent计算微通道收敛不了

    Fluent
    1
    0 赞同
    1 帖子
    1k 浏览
    D

    针对一个高度为3mm,宽度为1mm的矩形微通道,用fluent mesh 画网格后,在入口流速为0.3 0.4 m/s时,计算都可以收敛。但是当入口流速为0.5m/s时,计算收敛不了,这可能是什么原因导致的,尝试过增加网格数量也无效

  • 搅拌反应器边界层网格划分

    Meshy
    8
    0 赞同
    8 帖子
    12k 浏览

    @alex 看到你的网格了,感谢分享,感觉划分的非常不错,好厉害!学习了!

  • 滑移网格算法相关文献求助

    Algorithm
    3
    0 赞同
    3 帖子
    4k 浏览
    肖艺

    @东岳 谢谢!

  • 如何在动网格工况中如何设置移动测点

    Fluent
    1
    0 赞同
    1 帖子
    1k 浏览
    L

    请教各位:
    目前有一个6自由度运动的sloshing两相流模型,通过udf文件给定运动速度。

    需要测出壁面位置某些点的压强时程变化,以及某些区域的水位(波高)时程变化。请问该如何设置这些测点,使其与动网格区域一起运动。

  • 0 赞同
    1 帖子
    3k 浏览
    L

    各位专家好:
    在多孔介质模拟的过程中,我发现了一个问题,虽然在发表文章中省略了。但是确实是一个直观的问题,我的模型,RANS方法中,不管采用哪一种,在多孔介质的近壁面处的拟合的速度与实测值总是具备一定的差距。这种差距我分析认为是在较高阻力降的多孔介质环境下,气流在流经多孔介质表面时,一部分流体会因为高阻降而产生流向变化,形成比较复杂的近壁面流场。而在这个近壁面流场中,因为Fluent默认是采用的多孔介质域或者是多孔跳跃面,应该是对这层近壁面没有处理的,依然按照充分发展层处理,造成这层近壁面的流速(甚至是流向)与实测值都有较大区别。
    虽然能够较好的拟合实时压降(Endo based equations),但是如果在近壁面流场不符合实际的话,那么很显然对于最后的表面覆尘的预算难以符合实际情况。不知道大家对于这种多孔介质的近壁面流是否有过研究,或者能够有比较好的解决思路?欢迎大家提出宝贵的想法和建议,感谢,再次感谢!
    下图为提出新模型对于滤管外壁的模拟和实测结果对比:
    0_1544421734755_compare.png

  • 函数调用出现未知类型

    OpenFOAM
    3
    0 赞同
    3 帖子
    2k 浏览
    S

    @李东岳 在 函数调用出现未知类型 中说:

    const Thermo& thermo是什么变量?

    是个Thermo变量。Thermo是一个类型。

    谢谢李老师!该如何定义Thermo类型的变量呢?我想调用Tvenewdefinition这个函数,需要输入该类型的变量,但我输出thermo发现是个关于组分信息的矩阵,不知道该如何定义。

  • 0 赞同
    5 帖子
    9k 浏览

    是用八叉树算法吗?设置part上的最大尺寸也需要是2的幂吗?对这方面有了解吗?我也有这方面的一个疑问呢,之前在论坛发帖子还没人回。帮助文档里面也没写清楚,网上查了一些也没查到

  • 关于划分网格的无量纲数的定义

    Meshy
    2
    0 赞同
    2 帖子
    5k 浏览
    D

    @hhj0219 和y+定义一样。其实就是:流向网格间距除以壁面第一层距离,然后再乘以y+。

  • 颗粒剪切应力

    OpenFOAM
    3
    0 赞同
    3 帖子
    2k 浏览
    F

    @李东岳 感谢老师回复。是拉格朗日颗粒,我用的MPPICFoam