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

    恭喜恭喜

  • 0 赞同
    5 帖子
    6k 浏览
    T

    @bestucan 万分感谢!

  • 如何在把下面这个东西画成图

    CFD彩虹条
    8
    0 赞同
    8 帖子
    9k 浏览
    bestucanB

    @李东岳 :140: ,效果很好。还能把点调大点,调成正方形,然后就可以没有空隙了。pointtype, pointsize这些

  • 0 赞同
    2 帖子
    3k 浏览
    D

    @小令 你可以在plot—>blanking—>value Blanking中,选择把右侧的数据隐藏掉。

  • 自适应网格优化 求助!!!

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

    基于2011年发表在Journal of Computational Physics的文章Adaptive mesh refinement based on high order finite difference WENO scheme for multi-scale simulations,并根据您文章的网格自适应算法编写了一套计算程序,在计算双马赫反射算例过程中,遇到了一些问题,受困许久,希望大神指点迷津~ 在此诚挚感谢
    遇到的问题具体问题如下:
    在计算双马赫反射问题时,在Kelvin-Helmholtz rolls两侧流场出现较大的震荡,局部密度云图及网格自适应加密情况如图1所示,震荡区域存在粗细网格过渡,故怀疑在粗细网格过渡中存在算法错误或者代码错误(我已多次检查代码,未发现代码错误,但不能完全排除仍有代码错误存在)。
    0_1536922532180_9e8d3f5c-b1c2-4828-bec3-8e472986a4d2-image.png
    (a)局部密度云图(红框处为存在流场震荡处)
    0_1536922552890_10e26979-d84c-4d89-956b-e81584a6adb5-image.png
    (b)局部网格生成
    0_1536922574620_04fd03c7-eca8-4851-a7e4-3c342cae3939-image.png
    (c)总体密度云图
    0_1536922591916_530a5e07-8d83-41c6-b1e3-338c5acb0b3a-image.png
    (d)总体网格生成(黑色区域表示网格加密)
    图1 局部密度云图及网格自适应加密情况
    该程序采用C++语言编写,流场计算使用有限差分法,五阶WENO格式,Steger-Warming通量分裂,2D欧拉方程,时间推进采用3阶Runge-Kutta方法;网格自适应加密方法和加密准则与您的论文有所不同,采用基于原始根网格进行加密,每个网格单元生成3×3个子网格单元,加密准则使用密度梯度(∂ρ⁄(∂x+∂ρ⁄∂y))。一共使用三层网格,在计算域[0,3]×[0,1]中,原始根网格单元数量为300×100。Layer 0表示原始根网格,Layer 2表示最密的网格。Layer 2比Layer 1略大,覆盖Layer 1,两层网格的间隙为3个根网格单元的尺度,如图2所示。以此希望精准捕捉激波,并适当减少网格生成数量。子网格层的时间步长为父网格层的1/3。
    父网格层每计算完一步,对父网格节守恒量(ρ,ρu,ρv,E)和守恒量对时间的导数[ρ’,(ρu)’,(ρv)’,E’]进行五阶加权(类似WENO加权)空间插值,求解出子网格层边界的虚网格节点的相应守恒量及其导数。并将前一时间步和本时间步插值出的通量及通量导数使用Hermite插值,以获得三个子时间步的相应变量,作为加密网格层的边界条件,并使用式(1)作为Runge-Kutta时间推进的中间步的边界条件。
    0_1536922724975_5a02e4bb-a1eb-4fce-884b-9d73584cc27c-image.png
    在子网格层三个时间步计算完成后,父网格层被子网格层覆盖的节点数值被子网格对应同位置节点更新。每一或两个根网格时间步重新生成一次加密网格。
    0_1536922743488_5f7e9dc2-9e8c-44f9-88fb-73083e652ec1-image.png
    (a)粗网格和细网格
    0_1536922779864_ecb19b47-59e2-4210-97fb-4d0f9de2d195-image.png
    b)不同层次网格
    图2 网格加密示意图
    其他相关的情况介绍:
    在一次测试中使网格只加密,不减少,结果显示流场基本没有前述震荡,结果如图3所示。此算例应该可以证明该程序计算格式是正确的。然而,由于生成网格实在太多,不具实际应用意义。
    0_1536922795688_29716dff-e9d7-4afb-a475-b7a750cebe66-image.png
    (a)密度云图
    0_1536922809856_e6bd76e2-af3b-4966-9f15-c3e59d18ea3b-image.png
    (b)网格生成情况(黑色为加密区域)
    图3 网格只加密不减少的计算结果

    不知大家在调试程序过程中,是否遇到过类似的情况。
    希望大家能够指出
    (1)我对算法的理解是否存在错误,特别是粗细网格过渡部分的算法;
    (2)解决该问题有何建议;
    (3)能够告知粗细网格过渡区域,造成的流场震荡的原因可能有哪些;

    再次感谢大家!

  • 上下振动混合器的流场模拟问题

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

    东岳大神、各路大神好:
    最近在模拟一个上下振动混合器中的单相三维流场,如图。 考虑到计算时间和成本问题,想采用moving frame of reference的方法来代替动网格进行振动的模拟。目前设置inlet速度为振动速度,壁面和出口都设置的是wall,速度也是振动速度,而板和轴的速度为0.但得到的流场似乎不太对,想请教下边界条件的设置存在哪些问题?WeChat Image_20190929112939.jpg

  • CFD青年成长支持计划(2022)

    公告
    87
    0 赞同
    87 帖子
    216k 浏览

    我是南京航空航天大学航空学院2022级硕士研究生吴祥清,跟随张老师链接文本从事CFD研究,本人的主要研究方向为高速出入水,初步学习商软入门CFD现在准备编程写算法。查找学习资料了解到CFD中文网以及老师,所以想申请2022CFD青年成长支持计划,我目前正在使用LS-DYNA软件进行垂直破冰仿真,开始着手发相关论文。在本科期间接触流体力学便想深入了解这个专业,也申请了这个方向的研究生,如果论文写作顺利,便会在研究生二年级申请硕博连读,为国家军事建设贡献自己的一份力量。目前在上课之余也开始着手C++入门。但是CFD学习之路道阻且长,希望可以得到老师资助来更好的学习CFD。

  • turbulence->validate有何意义

    OpenFOAM
    5
    0 赞同
    5 帖子
    7k 浏览

    其实也不难,你用grep大概能猜出个所以然来。
    grep -rn "validate"
    其他都是空函数只有一个函数有点意思

    template<class BasicTurbulenceModel> void Foam::ReynoldsStress<BasicTurbulencemModel>::validate() { correctNut(); }

    然后我又grep 了correctNut()
    常用的是kOmega 我看了下

    template<class BasicTurbulenceModel> void kOmega<BasicTurbulenceModel>::correctNut() { this->nut_ = k_/omega_; correctBC; BasicTurbulenceModel::correctNut(): }

    我虽然没有一行行debug, 大概猜出来就是湍流初始化。最开始给定k omega 初始值计算出第一步的nut 初始值给k omega方程,即previous time step,欢迎指正。

  • twoPhaseEulerFoam中的IATE模型系数

    OpenFOAM
    34
    0 赞同
    34 帖子
    35k 浏览
    U

    @李东岳 对的,大概有十几公里,那谢谢东岳老师,我再使劲努努力钻研一下。个人觉得根据文献植入算法太难了,虽然能看懂离散和求解过程,但是想在openfoam里面实现觉得难度挺大的,当时老师您带我们植入的单相可压缩的,目前这是两相可压缩,我努力努力,再次感谢东岳老师的持续关心:xinxin:

  • 如何在计算中开启使用多重网格

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

    可以啊 这个跟几何没关系啊

  • 0 赞同
    1 帖子
    2k 浏览

    Inked捕获_LI.jpg 替代文字

  • 颗粒流动量方程中固相压力项的疑惑

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

    在fluent手册介绍固相压力处写到“The solids pressure is composed of a kinetic term and a second term due to particle collisions”,即固相压力包含动力及碰撞两个部分,不包含摩擦压力。但是在介绍摩擦压力时,又写到“The frictional stress is added to the stress predicted by the kinetic theory when the solids volume fraction exceeds a critical value.”所以,当固相体积分数超过临界摩擦体积分数时,固相动量方程中的固相压力项是否包括摩擦压力呢?

  • V锥流量计网格划分求助

    Meshy
    13
    0 赞同
    13 帖子
    22k 浏览

    @simthere 好久没有画结构网格了,回归下,练练手,暂时还能画出来,哈哈0_1503487106341_l1.png 0_1503487116733_l2.png

  • paraview作图求助

    OpenFOAM
    3
    0 赞同
    3 帖子
    3k 浏览
    Y

    謝謝老哥,後來用的歪門邪道搞定的,你説的這個方法不錯 @bestucan

  • 关于cht的网格划分问题

    OpenFOAM
    5
    0 赞同
    5 帖子
    7k 浏览
    R

    @李东岳 在 关于cht的网格划分问题 中说:

    writeZones

    万分感谢李老师的解答,转化格式(写入命令语句)为
    fluentMeshToFoam xxx.msh -writeZones,
    本人OpenFoam小白,望老师多多包涵:chouchou:

  • 压力方程和解

    Algorithm
    2
    0 赞同
    2 帖子
    4k 浏览

    这个是理论大坑啊。

    如果压力基的求解器,求解的是速度和压力场,那么得用连续性方程来求压力。
    如果是密度基的,连续性可以求密度,然后用状态方程求压力。

  • 流体区域内部可以设置进出口吗

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

    不可以,但是可以设置为风扇边界

  • 计算流体力学开源软件推荐

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

    PeleC(AMReX)

  • 如何提高混合网格的质量

    Meshy
    7
    0 赞同
    7 帖子
    13k 浏览

    @simthere 这个非结构网格可以在结构网格上生成吧,不需要强制合并节点,感觉这样会好点,强制合并出来的太差了

  • ANSYS Fluent中QMOM的植入缺陷

    Algorithm
    4
    0 赞同
    4 帖子
    7k 浏览
    S

    @东岳 谢谢:baobao: