Skip to content
  • 0 赞同
    4 帖子
    6k 浏览
    李东岳

    如果不考虑时间项
    \begin{equation}
    \int^{f+1}_ {f} \frac{\partial}{\partial x}\left(\alpha \frac{\partial \phi}{\partial x}\right)\mathrm{d}x=\left(\alpha \frac{\partial \phi}{\partial x}\right)_ {f+1}-\left(\alpha \frac{\partial \phi}{\partial x}\right)_ {f}
    \end{equation}
    把后面两项$\left(\alpha \frac{\partial \phi}{\partial x}\right)_ {f+1}-\left(\alpha \frac{\partial \phi}{\partial x}\right)_{f}$把所有的面连起来,看起来是守恒的,剩下的只有边界条件

  • Laplace测试,对表面张力模型进行测试

    Fluent
    3
    0 赞同
    3 帖子
    4k 浏览

    @guanghuilineu 请问最后解决这个问题吗?你这个液滴是三维的吧?我最近也在测试VOF计算表面张力的效果,误差也不小,但没你这么大

  • 0 赞同
    3 帖子
    3k 浏览
    C

    您好,我在使用comsol网格自适应算两相流会出现网格自适应速度跟不上两相运动速度,导致,这方面您是如何解决的呢

  • 大深宽比微通道结构网格绘制

    Meshy
    2
    0 赞同
    2 帖子
    2k 浏览
    李东岳

    这种没法搞,注定就是深宽比很大

  • 0 赞同
    6 帖子
    8k 浏览
    H

    ef44221d-062f-4315-b9c1-50e96688d88e-image.png

    一个公式真是折腾死了。

    我看了好几个人讲关于von Neumann稳定性分析。

    12步入门Python这个女老师讲的。 MIT qiqi wang讲的PDE课程

    还是有些地方很模糊。

  • 0 赞同
    2 帖子
    3k 浏览
    李东岳

    这个要是openfaom的网格,就好处理了。2亿网格确实太多了,icem我也只画过最多几百万的网格。卡的都不行。你这2亿网格太疯狂了,机器内存应该很大。

  • 0 赞同
    2 帖子
    4k 浏览
    李东岳

    @流动魅力 我比较好奇你的相间作用力的直径怎么给的

  • 0 赞同
    1 帖子
    2k 浏览

    想开个贴讨论下多相流界面追踪类方法目前的前沿水平,例如能算到的最高参数,各类方法都有什么限制等等。目前多相流算法里,计算界面追踪的多相动力学,比如droplet collision, droplet splashing, droplet impact on solid wall,这几个典型案例,目前有用VOF,level-set,phase field,LBM, SPH, 欧拉拉格朗日耦合其它模型的。这几个案例在不同Re和We下会产生许多种演化形态,实验都做出来了,但感觉计算还没跟上。

    调研过不下百篇的文献,基本上关键词相关的都翻了一遍,在文献中看到droplet collision大概最好的也就是Re=6500,We=1500,密度比DR=1000左右的,大致是这个级别,像droplet impact有篇08年的用TransAT算到了RE=11000,而最近16年左右的用levelset这些也算到Re=8000,We=1500左右。大多数都是采用了一些改进的数值技巧来稳定模拟,而且网格也都比较密,计算资源消耗也都比较大,估计都得上超算,当然和实验也都对的上,基本都是到了prompt splattering这样的形态。

    想讨论下大家都做过或者见过的这类模拟最高参数是多少? 或者往上算各类方法的限制在哪? 因为多相流在算到高RE和We的时候会数值发散,稳定性应该是主要限制,还有例如VOF和level set会在薄液层处由于数值格式产生人工破碎的非物理情况,需要单独使用人工识别薄液层再稳定的方式,对不规则碰撞并不实用。

  • fluent 模拟油滴的对流扩散

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

    各位老师w71h11131532_1601213641_696.png ,我想模拟一下油滴在超临界水里下落过程中的对流和扩散问题:
    初步打算用VOF模型来捕捉相界面,但有一个问题VOF并不能解决扩散的问题。请问有没有什么比较好的解决方案?
    如果用UDF/UDS来解决扩散项的问题,请问能不能分享下类似的代码?

  • 0 赞同
    7 帖子
    12k 浏览
    X

    @李东岳
    参考文献:
    Knacke, T. (2013).
    Potential effects of Rhie & Chow type interpolations in airframe
    noise simulations. In: Schram, C., Dénos, R., Lecomte E. (ed):
    Accurate and efficient aeroacoustic prediction approaches for
    airframe noise, VKI LS 2013-03.

  • 0 赞同
    5 帖子
    7k 浏览
    zzkluckZ

    说一点我的看法,不一定对,仅供参考:

    non-orthogonality

    dca809a8-9cd5-4f74-9329-6d8459966264-image.png

    non-orthogonality由snappyHexMeshDict.meshQualityControls.maxNonOrtho直接控制,你的snappyHexMeshDict里这个值是45,这个值一般不需要控制的这么小,65是更常见的配置值。

    另外,你的snappyHexMeshDict里,relaxed.maxNonOrtho给的是75,这个值通常是和maxNonOrtho = 65相配合的。如果你真的要把maxNonOrtho设为45,按逻辑来说relaxed.maxNonOrtho也应该相应的减少一点。

    网格纵横比

    这里截取你的背景网格配置:

    vertices ( (-340 -302 -3) // Slightly smaller than the STL model bounds (332 -302 -3) (332 341 -3) (-340 341 -3) (-340 -302 95) // Slightly larger than the STL model bounds (332 -302 95) (332 341 95) (-340 341 95) ); blocks ( hex (0 1 2 3 4 5 6 7) (100 100 100) simpleGrading (1 1 1) // Adjust mesh density as needed );

    如上面B老师所说,你网格的大纵横比就是这里导致的。如果不是有意要在Z方向上做加密的话没必要这么分块,我口算一下,大概(70 65 10)就可以。

    castellate

    从上面可以看到,你背景网格的量级是100 * 100 * 100 = 1,000,000;而你SHM的maxGlobalCells给了2,000,000。这看起来不太对,这几乎没给SHM进一步细分的空间,从结果来看也是这样,网格划分的不够细。

    另外,截取一段level的配置:

    features ( { file "building.eMesh"; level 3; } ...... ); refinementSurfaces { building { level (3 4); } ...... }

    我理解features的level应该不低于Surface的。

    胡乱写了一些个人看法,希望能有所帮助。

  • STAR-CCM+网格转换后,checkMesh报错

    Meshy
    2
    0 赞同
    2 帖子
    3k 浏览
    I

    1.我猜是symmetry这个边界条件,试一下改成wall,然后在0文件夹里自己定义上下的slip边界条件。openfoam会对生成的网格编号进行检查,外部软件的网格节点编号和blockMesh的逻辑不一样。
    2.可能是前后的empty边界问题,你可以在starccm里生成3D网格后,重新在openfoam里extrude来生成前后empty边界
    类似这样
    constructFrom patch;
    sourceCase "../w3-d_hc1-3";
    sourcePatches (symFront);

    // If construct from patch: patch to use for back (can be same as sourcePatch)
    exposedPatchName symBack;

    // Flip surface normals before usage. Valid only for extrude from surface or
    // patch.
    flipNormals false;

    //- Linear extrusion in point-normal direction
    extrudeModel linearNormal;

    nLayers 1;

    expansionRatio 1.0;

    linearNormalCoeffs
    {
    thickness 1;
    }

    // Do front and back need to be merged? Usually only makes sense for 360
    // degree wedges.
    mergeFaces false; //true;

    // Merge small edges. Fraction of bounding box.
    mergeTol 0;

  • 北师大学生性行为调查报告2016版

    CFD彩虹条
    3
    0 赞同
    3 帖子
    6k 浏览
    C

    @Wayne

    我已联系你:sunglasses:

  • 崩溃啊,这个符号什么意思?

    Algorithm
    5
    0 赞同
    5 帖子
    7k 浏览
    李东岳

    确实和概率论有关,谢谢,记录一下

    \begin{equation}
    \left( \begin{matrix}
    k\\
    x
    \end{matrix}\right )=\frac{k!}{(k-x)!x!}
    \end{equation}

  • 关于fluent19.0

    Fluent
    7
    0 赞同
    7 帖子
    9k 浏览
    L

    @发芽的土豆 :
    感谢赐教,多谢了

  • 0 赞同
    3 帖子
    5k 浏览

    @haaDoo 物性参数都可以变,温度不一样物性参数就不一样,或者换不同的流体介质

  • OpenFOAM大涡模拟的滤波尺度

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

    @zhangxy 是一个网格的尺度

  • 0 赞同
    2 帖子
    5k 浏览

    退Fortran用c++保平安,虽然没用过Fortran,但看知乎上专门论述Fortran的缺陷,发现用c++一下就好了。但是c++的学习周期长,东西多,也是问题。

  • 帝国理工

    CFD彩虹条
    1
    0 赞同
    1 帖子
    4k 浏览

    0_1460680809776_捕获.JPG

  • Fluent 如何设置降雪粒子

    Fluent
    2
    0 赞同
    2 帖子
    3k 浏览
    I

    @caizhao https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1002/2016JD025316
    可以参考一下这篇文章
    上面链接打不开的话,可以用下面这个
    https://agupubs.onlinelibrary.wiley.com/doi/full/10.1002/2016JD025316