Skip to content
  • 0 赞同
    3 帖子
    5k 浏览
    B

    尖的地方可以试试将多余的节点进行合并成一个节点映射到尖尖的点上,然后对锥块做一个Y剖,但是这样的话在尖尖位置的网格要么太密要么质量差

  • 0 赞同
    7 帖子
    10k 浏览
    chengan.wangC

    @黄昏晓为谁破 在 第三代涡识别方法liutex的Fortran程序运行出错原因,求助 中说:

    @Vortex学习 我没用过Fortran程序计算Rortex,不过我的方法可以借鉴下,就是使用Tecplot打开Fluent结果,并导出基于节点的data文件,里面只包含xyz坐标以及UVW速度,然后我是在matlab中进行处理获得的R值以及矢量值,在替换data文件里的数值,最后再使用Tecplot打开data文件。这样可以处理较大网格量的数据,不过这个网格量得取决于Tecplot软件,我之前有个2000万的网格,Tecplot无法打开,我就很无奈了。

    您好,这么多节点的UVW速度,如何在matlab里面计算速度梯度张量呢?

  • MFIX-DEM的并行性能验证

    已移动 OpenFOAM
    1
    0 赞同
    1 帖子
    4k 浏览
    散漫守望2016

    Liu et al., 2016 "A Comprehensive Benchmark Suite for Simulation of Particle Laden Flows Using the Discrete Element Method with Performance Profiles from the Multiphase Flow with Interface eXchanges (MFiX) Code" 一文中对MFIX-DEM的并行能力进行了验证,采用了Weak Scaling和Strong Scaling两种评估方法,结果证明在核数达到1000左右,MFIX-DEM code的并行能力还很不错。但是注意到,在做自由堆积等的时候,文章回避了负载均衡的问题,也就是如果流化床下半部分非常密集,而上半部分非常稀疏,这是对就需要一种新的分区方式,使得每个核都能计算差不多的颗粒数。现在的MFIX好像还达不到这个功能,但我知道有人开发出来了。文章还认为,当颗粒数达到10(8),即使用1000核,每个核也有10(5)个颗粒,还是非常有挑战性的。气固两相流的计算时间主要花费在了DEM颗粒彭碰撞计算,颗粒和流体之间的数据交换,曳力的求解。
    看到文章的一点心得,和大家分享一下。

  • 稳态求解器算出的湍流拟序结构

    Algorithm
    5
    0 赞同
    5 帖子
    8k 浏览
    C

    @evensun 求解器只负责求解数学式(控制方程),为什么要求解这个控制方程而不是别的控制方程是用户的选择,或者说取决于用户需要的是什么。
    如果你需要百分百精确的描述一个非定常的物理问题,那你用稳态求解器原则上来讲是不合理的,但这个不合理不是说解方程这一步不对,而是你没有求解与你的需求相匹配的方程。如果你只需要部分描述这个物理问题,比如你关注时间平均量,那这时候用RANS去解这个问题,这一套下来是合理的。因为你需要的只是时均量,只是这个物理问题的一部分。
    从这个角度去思考,你用稳态求解器计算了一个非定常的过程,求解器没有发散,甚至可能还体现出了一些有规律的震荡。这个现象本身其实是有意义的,但你在分析它的时候要非常小心。比如你认为“算得的结果呈现出了涡的发展变化过程”,这个表述是值得商榷的。

    再啰嗦一点,关于RANS我一直觉得有个很有意思的角度。我们都知道雷诺数是决定流动的重要因素,我们可以认为控制方程也有一个“数值雷诺数”,它和普通的雷诺数定义类似,只是粘度项应该替换为分子粘度+数值粘度。在RANS中,数值粘度一般高出分子粘度2个量级甚至3个量级,所以尽管物理问题可能雷诺数很大,但数值雷诺数可能只在10的3次方或者4次方量级。10的3次方已经和我们通常认为的临界雷诺数接近了,导致RANS结果看起来是稳定的,就好像层流一样。
    说白了,RANS模型对湍流的耗散非常强,把很多非定常的结构都给耗散掉了,留下了定常的部分。但这个耗散能力是有限度的,如果你求解的这个物理问题经过RANS模型后它的数值雷诺数正好处于临界雷诺数附近,那可能就是会有一些小的震荡。
    这个湍流模型是自己开发的吗?你可以通过尝试不同的湍流粘度来验证一下,比如量级直接调大10倍,你说的这些震荡还存在吗?

  • 求解压力方程中的nNonOrthCorr

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

    @搬运工不好当

    在对拉普拉斯项离散的时候,如果使用非正交网格会引入误差,因此需要非正交修正。所以nNonOrthCorr的作用就是消除非正交网格引入的误差。

  • 0 赞同
    2 帖子
    2k 浏览
    N

    使用fluent仿真同心圆管的融化过程,没有打开重力场,瞬态仿真,打开了能量方程、层流模型、融化凝固模型,密度、导热系数和比热容是分段线性,其余都是常数,仿真发现沿轴向有速度产生,但速度很小,数量级大概在e-7左右,固、液相区速度沿轴向指向相分界面,请问这是因为什么产生的速度场?并且由于速度场的产生,融化速度将近翻倍。

  • 0 赞同
    1 帖子
    3k 浏览
    袁宝强

    请问各位大佬,对于可压缩流大涡模拟的总湍动能如何定义呢?

    对于不可压缩,大涡模型采用的亚格子模型采用的雷诺时均的形式,因此可以用雷诺时均计算大涡湍动能和亚格子湍动能,两者的和就是总湍动能。
    $K=0.5*(\overline{uu}-\bar{u}\bar{u}+\overline{vv}-\bar{v}\bar{v}+\overline{ww}-\bar{w}\bar{w})+\bar{k}_{sgs}$ 但是对于可压缩流的大涡模拟,亚格子模型采用的为法夫尔时均形式,此时如何定义总湍动能呢?我的理解是
    $\widetilde{K}=0.5*(\overline{\rho uu}-\bar{\rho}\widetilde{u}\widetilde{u}+\overline{\rho vv}-\bar{\rho}\widetilde{v}\widetilde{v}+\overline{\rho ww}-\bar{\rho}\widetilde{w}\widetilde{w})+\widetilde{k}_{sgs}$
    请问各位大佬,我的理解对么
  • DEM如何考虑传热

    已移动 OpenFOAM
    2
    0 赞同
    2 帖子
    5k 浏览
    散漫守望2016

    在liggghts/src下面找到这个fix_heat_gran_conduction.cpp文件,里面讲了颗粒和壁面之间的静态导热的计算,看看代码应该可以知道的。

  • fluent meshing.怎样找到意外出现的baffles

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

    可以在join/intersect里面查找一下看看。

  • 0 赞同
    2 帖子
    4k 浏览
    W

    用surfaceCheck检查STL文件的结果呢?

  • 0 赞同
    2 帖子
    5k 浏览
    S

    如果我没记错的话,流体中文网上的ICEM版块里有人画过4119桨的结构网格,也传了文件在上面,你可以去上面找一下看还在不在。(耐心往后翻一下)

  • checkMesh failed,求助

    Meshy
    6
    0 赞同
    6 帖子
    9k 浏览
    Y

    @东岳 老师,我重新画网格后,decomposePar经常会报错,但是不影响后面计算,没有查到原因0_1537492983872_0403f9aa-ee80-4134-a228-ac6289ba58de-image.png

  • 计算流体力学软件研发工程师招聘

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

    央企中国电子信息集团所属企业长沙军民先进技术研究有限公司招聘
    数量:6人
    工作地点:长沙市
    岗位名称:流体计算软件研发工程师
    岗位要求:
    1、流体力学相关专业;
    2、熟悉GPU和CUDA编程;
    3、熟悉C++编程;
    以上三项条件满足二项即可,满足三项者优先录取,流体力学强相关专业者也可面谈。

    基本月薪待遇8K-15K,其他待遇面谈,特别优秀的硕士博士待遇还可以商谈。
    联系人:郭先生,18908480199
    bf037a90-2b9b-464c-83da-3b501a4e6431-image.png

  • 论坛搜索功能下线了?

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

    @李东岳 我发了帖子之后才发现你之前发过通知了,不是催更,抱歉抱歉~:chitang:

  • Fluent体积力驱动两相流

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

    你是不是开启了相间的传质? 还有一个 温度会引起流体密度变化,对应是不是体积也变化?

  • 小白问个问题。。。

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

    @搬运工不好当

    CFD界:相对于商业代码,OpenFOAM的特性是什么?

    Henry:OpenFOAM是完全免费并且开源的(CFD界注:开源的意思即为使用者可以随意更改软件的内核)!这提供了HUGE的方便,不仅仅在费用上,也表现在灵活性、以及从软件本身来学习CFD上。OpenFOAM不仅仅让CFD工程师可以模拟算例,它也允许CFD工程师来学习CFD方法、开发CFD算法并完善,甚至某一天创造一个可以替代OpenFOAM的软件。
    并且,许多CFD问题需要一些特别的功能、算法、方程、边界条件等。在商业软件中植入这些功能可能令人沮丧并且很可能失败。有些情况下可能有必要完全的重新写一个求解器、边界条件等。OpenFOAM提供了最大的灵活以及扩展性。这些对使用CFD做研究的人非常重要,因此我觉得OpenFOAM可能在某一天会替代闭源软件。

  • 0 赞同
    5 帖子
    8k 浏览
    F

    可以设计一个量纲的类,里面包含一个整形数组来记录量纲的分量,如[L, M, T]还要定义量纲的运算。然后设计一个标量的类,里面除了浮点数据之外,再加一个量纲,定义标量的动作时也要注意其中的量纲计算。同理还可以定义向量和张量。比如距离的量纲是[1, 0, 0],速度的量纲是[1, 0, -1]。

  • 0 赞同
    17 帖子
    20k 浏览
    G

    @chszkc 谢谢大佬!!学到啦!!!

  • 0 赞同
    3 帖子
    4k 浏览

    @gemini 不应该直接放进gx 这成了纵荡了, 纵摇的加速度要在三维进行分解才能放进去。

  • 0 赞同
    8 帖子
    11k 浏览
    hurricane007H

    @hxy970122 扫了一眼Fluent的组分输运方程,对于各个组分就没有重力的影响,所以如果组分分离就肯定是因为扩散率不同。但是在binary mixture里面,似乎两个组分的扩散系数是相同的?

    重力的影响是在质量方程和动量方程里面,所以就是所有组分受到的重力作用是平均重力……

    2e37529c-113a-4d19-a63f-e23bf85c3593-image.png