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

    @Samuel-Tu 在 东岳老师CES论文中两相kEpsilon模型实现的疑问 中说:

    这说明a相的粘度是靠b相的粘度计算来的,但是我没有找到相关代码,希望大家能告诉我一下文件路径,让我学习一下如何实现。

    你可以参考下Gosman 1992年那篇文章

    对于公式(32)中C3项(绿框部分),是否少了一个epsilon。因为源码C3项中含有epsilon:

    应该是少了个epsilon

  • 这是什么类型的非结构网格

    已移动 Meshy
    3
    0 赞同
    3 帖子
    5k 浏览
    李东岳

    第一个是六面体的非结构网格
    第二个看不到内部,看起来像是四面体

  • 0 赞同
    1 帖子
    3k 浏览

    想找一些适合cfd入门学生学习的题目
    最好有把cfd的各个环节单列出来的题目(例如 贴体坐标系生成)
    以及这些问题用什么工具求解
    我是指学习使用哪些工具对于科研实际具有意义(例如使用哪些软件(fluent?turbo?) 还是用哪种语言自己编写代码(例如cpp)?有哪些部分已经出现类库,哪些需要自己编写呢?尤其是科研中是怎样实际操作的)

  • gambit最大网格数问题

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

    Solver 是Generic

  • 0 赞同
    4 帖子
    5k 浏览
    M

    @东岳 感谢东岳老师,意思就是 在 P 文件中写成
    inlet
    type codedfixedvalue
    value -ρ*(V-V'/δt)
    这种格式就好么
    东岳老师写的这个Δt是时间步长么,V-V(t-Δt)要怎么表述呢
    感谢东岳老师,谢谢你的指导

  • ICEMCFD中的remote有什么用?怎么设置?

    Meshy
    4
    0 赞同
    4 帖子
    5k 浏览
    bestucanB

    @东岳 ^_^倒不是嫌累,如果累的值,精神意志力会强化也要干。这种累有种被呼来喝去的疲劳感。这个remote在icemcfd的手册里解释的很清楚。

    如果回答,就类似于,人一张嘴问现在几点,我得屁颠屁颠去看看表然后汇报……

    他如果问如何连接,或者照着手册操作不通,还能共同探索一下。
    哪怕问“哪里能找到相关资料呢”

    甚至,在不知道有什么用的情况下,就问怎么设置。像不像导师说:“我也不知道为什么要往这个方向走,我让你走你就走吗”的感觉
    :135:

    当然我对题主没有任何恶意或建议或批评。只是谈下自己时常向大佬看齐,长期从事保姆行业,即将离职的疲劳感和惋惜感。和安慰一下题主,随问随答的人都被客观环境淘汰了,并不是世界冰冷
    :chouchou:

    以下是保姆的最后一保:
    remote就是远程的电脑,也叫计算节点(所以很多服务器默认host名称就是node)。为什么用计算节点,当然是自己用的电脑算力不足,需要外援。有两类目的。一类是图形渲染,一类是计算。远程计算节点,不只前处理(icem)可以用,计算(fluent)也可以用,后处理(paraview)也可以用。只要工作量够大的软件一般都带加远程节点的功能。

    这是个和“并行”很类似的概念。不过“并行”大部分情况下是同时进行的。这个更偏分布式计算,但分布式计算是不必同时进行的。比如这种公益的分布式计算: 即伯克利开放式网络计算平台 。发布任务、招募志愿者、安装客户端、领取计算任务、下载计算数据、返回结果。就像众筹,不过筹的是算力。

    如果说并行比串行高一个纬度,分布式计算(remote)就比并行高一个纬度。每个节点上都在并行跑,然后几个节点一起跑。

    不过有些软件上的remote都是假的分布式计算,看着是让节点来帮忙,大家一起干,其实是全给节点干……我只当个客户端。

    连接走网络协议,先开服务端的相同软件的服务器模式,然后在自己这边设置IP、端口、帐号密码验证(一般是ssh)

    我要开始无情了,无情地只回答认真提问的。
    :threaten:

  • 什么是averaging turbulence?

    Algorithm
    6
    0 赞同
    6 帖子
    10k 浏览
    .J..

    @showhand 多谢大神回复!

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

    没想到图越来越多。手机看的注意了。我这个图多。

    结构网格的拥护者简单说一下为什么我喜欢结构网格。以下只针对简单几何以及不是太复杂的几何。复杂到不行的网格我选择混合网格。我自己做的算例不用纯四面体网格。

    六面体网格具有以下优点:

    非矩形网格的歪斜(skewness)严重
    这需要做面法相梯度计算修正。稍微增加计算步骤。著名的CFD非矩形修正因此而来

    截面重构不光顺
    举例,一个正方盒子,一半是水,我们用矩形网格重构的气液截面是这样的:
    alt text
    黄色为一相,白色为一相。即使网格很糙,我们的截面也是平的。下图为三角形非结构网格:
    alt text
    可见糙网格下界面不平。这导致后处理需要耗费精力。

    高度非线性物理问题影响矩阵对角占优
    这也是非线性动力学,例如金属撞击、生物力学分析中结构网格肆虐的原因。(这个我没有测试,之前在网上看到的结论,来自链接文本

    非结构网格数量更多增加计算资源
    网格单元的面越多网格越少。六面体网格并不是数量最少的。数量最少的是多面体网格(star推荐的?)如下图
    alt text
    再看这个情况对比:只有第一个是结构网格,第二个网格多。第三个第四个和均为非正交且和流线不贴合
    alt text
    alt text
    因此目前自动生成的非结构网格一般数量是即为巨大的。计算的时间会超过你花一星期做结构网格的时间。另外,万一中途发散了怎么办。

    费结构网格高度不可控
    几何稍微一点点的变化会导致全场四面体网格做位移。结构网格(ICEMblock切出来的)具有即为即为高度可控性。可以对基本所有的网格单元进行控制。俩者区别天差地别。

    我了解的所有的离散格式对结构网格适应性更强
    也就是说结构网格对流场具有强的多的顺应性。例如旋转机械,结构网格和速度流线高度贴合。

    矩形网格对边界更好的捕获
    为何非结构网格要添加边界层网格?矩形网格对壁面贴合更好。

    复杂CFD模型在结构网格上收敛更强
    个人经验。比如双流体模型,LTS局部时间步框架下的模型,激波捕获等。当然那种不可压缩单相流不算

    太多就不写了...

    引用从业CFD27年OpenFOAM创始人Henry的一段话:

    hex and polyhedral cells are generally superior to tet cells in finite-volume and the difference in both accuracy and convergence behaviour can be large. However, it also depends on the quality of the hex and tet cells and many of the problems with tet meshes stem from the poor quality of some tets produced by many automatic tet mesh generators. In order to really appreciate the difference it might be a good idea for you to try both for your problems. For VOF calculations I would expect the solution on hex/poly cell meshes to be dramatically better.

    大意为:

    在有限体积中,六面体和多面体网格相对于四面体网格优势明显。主要在精度和收敛上。对于VOF(我上文提及的第2点),结构网格的优势是大大的。

    最后,需要耗费大量时间生成的六面体网格至今没有消失,必然有存在的道理。

  • 淘书了!

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

    淘书了!

    Macroscopic Transport Equations for Rarefied Gas Flows: Approximation Methods in Kinetic Theory

    Turbulence Modeling for CFD

  • 为什么LES具有随机性?

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

    这个问题也可以参考我之前的回答:

    http://120.25.245.130:4567/topic/94/les和rans中-离散的对象是相同的吗/3

  • 求助UDF并行计算的问题

    Fluent
    5
    0 赞同
    5 帖子
    7k 浏览
    waltenW

    @关伯兰在笑 逐行排查吧,所有分母都添加一个很小的量。

  • OpenFOAM网格序号重新命名?

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

    @东岳 谢谢老师,我已经找到网格生成规律了,从左至右,由下至上地规律依次生词。满足误之前需求工作,reNumberMesh有需要我再去琢磨吧,再谢!

  • Wind-US Flow Solver有没有大神了解这个

    Fluent
    5
    0 赞同
    5 帖子
    6k 浏览

    @东岳 因为时间成本问题吧,现阶段还没有这方面的打算,之后应该会学。。。其实只是想问下这个求解器,如果简单方便入手的话想要了解一下。。。

  • 并行运算,然后突然少了一个进程

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

    @星星星星晴 感谢

  • Moving reference frame

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

    好像foam_extend里面已经有了Rusche这个代码?你看看?
    这个算法我看过但没自己植入过。因此这些细节问题也不好说,没什么建议了

  • 0 赞同
    7 帖子
    8k 浏览

    @东岳 有关异相反应,我刚开始的想法是将煤等量换算成为气相,后来觉得比较困难,所以就开始考虑再加入一相离散颗粒,还能利用原本的一个煤燃烧案例,东岳老师,你觉得这个有可行性吗?

  • 在Euler-Euler类的求解器中估算界面速度

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

    关注点不太一样吧,我是做雾化模型的,主要是算的相对准,计算量要小,现在大佬们做雾化直接DNS,Quasi-DNS,SPH的@东岳

  • ICEMCFD14.0实例源文件

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

    @SnoWolf2018 @玉玺

    多谢!!

  • 形状优化

    Meshy
    6
    0 赞同
    6 帖子
    10k 浏览
    T

    这个是不是用到了基于物理量+动网格技术。。。

  • fluent中关于喷雾降尘的问题

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

    1.个人思路是在DPM中把雾滴设置成一个injection,粉尘颗粒设置成另一个injection,然后运算,想看到这个喷雾开启前后粉尘浓度的变化,但是冥冥之中感觉没有这么简单,希望各位老师能指点一下如何实现这个过程?
    2.DEM可以模拟两种不同颗粒的碰撞过程吗?比如水滴和粉尘碰撞。