Algorithm

392 主题 2.0k 帖子
  • 一个小笔误让我搞了2小时

    4 帖子
    3k 浏览

    咦,这个工作我应该是2015年做的,怎么2017年才发帖..

  • 如何求整个流体域的y+

    2 帖子
    2k 浏览

    整个域中的摩擦速度用的是当地的还是壁面的呢?

    壁面的,我刚刚验证过,确定

    替代文字

  • 关于一个表达式的笛卡尔分解问题

    5 帖子
    2k 浏览

    @wwzhao 再次感谢:xinxin:

  • 12 帖子
    4k 浏览

    @金石为开 Any updates?

  • 关于变量分布假设的一个问题

    4 帖子
    2k 浏览

    看起来是的,找了下资料没看出来理由,看看大家什么看法?

    倒是发现一篇0阶精度的文章..一会发个CFD界

  • 升力、曳力、虚拟质量力等

    4 帖子
    4k 浏览

    @东岳 我也在做测吹方面的研究,但是E-E有点算不起(因为我的吹入速度很大)东岳先生的结晶器图片我看过,很漂亮。

  • 4 帖子
    2k 浏览

    是的,我觉得,本质上,还是应该走看书 - 手算 - 编程这个过程,如果看书直接去编程,可能理不清思路。

  • 新的Ryzen有用过的吗?

    已移动
    30 帖子
    21k 浏览

    @东岳 过两天我测试一下6142

  • 请教一下向量场的系数矩阵长什么样子?

    11 帖子
    3k 浏览

    又研究了一下明白了:如果手动添加各项异性的全隐式源项,
    需要拷贝一份A()并且构造成为张量,然后添加隐式源项,之后再构造HbyA,这样就能将各向异性的源线用全隐的方式添加到面通量上面了。

  • 2 帖子
    1k 浏览

    什么是粘性边界条件,粘度的边界条件?粘度不都一样嘛 :quwan:

  • 一个关于Fortran调试的问题

    3 帖子
    2k 浏览

    谢大婶答。
    一开始看到底下三个按钮直接崩溃了(原谅我两天两宿基本没睡觉,实在冷静不下来)。。后来发现了是给定的波浪参数有问题。。。
    PS:原来都说C语言亦或是Python啥的从入门到放弃,自从我看到了Fortran里面的各种变量定义之后还没有注释,要求硬生生理解看懂程序,于是。。。
    嗯,Fortran这种东西,哇,真是美妙。
    想退退不了呢,(^o^)/~

  • 5 帖子
    2k 浏览

    @东岳 哦哦,话说劳动人民筛麦糠好像用的就是这个原理:quwan:

  • 二阶张量的高斯公式如何理解?

    5 帖子
    2k 浏览

    @bestucan 3x3,不是2x2

  • S-A模型是如何处理湍动能k的?

    5 帖子
    2k 浏览

    @wwzhao 嗯嗯,明白了,非常感谢您的解答!:xiexie: :xiexie:

  • 1 帖子
    947 浏览

    雷诺数从物理意义上说是惯性力比上粘性力,想问为什么雷诺数的定义会涉及到惯性力这么一个虚拟力,是怎么考虑的,对于管道流动来说,如果是匀速运动,则流体的压力和粘性摩擦力平衡,对于所研究的流体单元来说是不存在惯性力的,如果管道流动存在加速度,则流体的压力和粘性力的合力导致流体的加速,可以采用一个虚拟的惯性力来等效流体的加速,则可以认为流体在压力、粘性力和惯性力三者作用下平衡。 所以,雷诺数为什么没有压力项,而用一个虚拟的惯性力来定义?

  • 12 帖子
    5k 浏览

    @lin_lia 三阶迎风格式的界面插值定义如何推导啊?也就是数值传热学第5章p184页习题5-9,对有限体积法截断误差的定义还是不太清楚?感觉有限差分法更直观啊,课本上p165页给出的格式,用泰勒展开法就可以证明,但界面插值定义该如何推导?

  • 4 帖子
    2k 浏览

    哦,这样啊,先看完这本书,再学新的。

  • 2 帖子
    1k 浏览

    直接相乘即可: $\mathbf S_b \cdot (\phi \mathbf u)_b$

  • Fluent中的VOF与组分运输模型

    1 帖子
    1k 浏览

    老师们,使用VOF能模拟气液两相流,使用组分运输能模拟同相态的质量传递。
    请问, 能用VOF与组分运输,模拟气液两相流的传质吗?

  • 19 帖子
    6k 浏览

    是我应该感谢cfd中文网的各位老师,对我的帮助很大.
    :xiexie: :xiexie: :xiexie: :xiexie: :xiexie: :xiexie: :xiexie: :xiexie:
    另外,我测试的这个构造函数和析构函数只承担了new和delete一个double[]的任务,如果是功能更复杂的构造函数,就不好说了。

  • 二维贴体坐标和常用坐标的转换问题

    3 帖子
    1k 浏览

    Fortran 代码,令人感动

  • front tracking 怎么翻译?

    4 帖子
    2k 浏览

    @东岳 不敢班门弄斧,我一直是你的粉丝:ok: :ok: :ok:

  • 多相流这面比较重要的研究方向

    1 帖子
    1k 浏览

    谁能帮忙填充一下没考虑到的内容么?

    捕获.PNG

  • 一个包含CFD各种研究方向的图

    3 帖子
    4k 浏览

    捕获.JPG

    重写了一个,以后用这个更新

    用的Adobe AI

  • 3 帖子
    2k 浏览

    @东岳 谢谢版主的解答。。这个问题暂时搁置了,以后有机会把UDF分享出来。

  • 求关于LBM方法的算例或者教程

    13 帖子
    7k 浏览

    @FluidGao 嗯嗯,感谢。

  • 欧拉方程按照声速传输的特征变量

    1 帖子
    1k 浏览

    前几天OKS的时候有个同学问欧拉方程是不是速度U和压力p按照声速传播。当时我给的答案是不是。在这里详细说一下这个问题。对于双曲系统
    \begin{equation}
    q_t+Aq_x=0
    \end{equation}
    其中$q$是求解的矢量,比如$[p,u]$;$_t$是对变量求时间导数,$_x$是对空间求导数,$A$是一个矩阵,但不是对角阵。双曲系统中的$A$是可以对角化的:
    \begin{equation}
    A=R \Lambda R^{-1}
    \end{equation}
    其中$R$是特征向量组建的矩阵,$\Lambda$是特征值组件的对角阵,如果定义$w=R^{-1}q$,最原始的方程就可以变成
    \begin{equation}
    w_t+\Lambda w_x=0
    \end{equation}
    这个方程式一系列解耦的双曲方程,可以所有的$w$都按照特征值的方向进行传播。因此对于欧拉方程,p和u并不是按照声速传播,按照声速传播的是$w$,$w$可通过$w=R^{-1}q$求出。

    $w$也就是特征变量哦。如果感觉困惑,或许可以看看这个:http://homepages.see.leeds.ac.uk/~amt6xw/Distance Learning/CFD5020COMP/node17.html

  • 2 帖子
    2k 浏览

    如果能实现顶盖驱动MAC算法,simple也没有问题啊~可以先试试MAC,速度边界条件都是第一类边界条件,其余都为0.只有顶盖有速度。

  • CPU在并行时候的指派

    5 帖子
    3k 浏览

    @程迪 请问linux如何并行计算,在ubuntu下没有问题,但是在linux下并行,显示这样的错误,希望您指教

    OpenFOAM-6(23) mpirun -np 2 interFoam parallel -------------------------------------------------------------------------- The value of the MCA parameter "plm_rsh_agent" was set to a path that could not be found: plm_rsh_agent: ssh : rsh Please either unset the parameter, or check that the path is correct
  • 有限面方法 Finite area method

    2 帖子
    1k 浏览

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

  • 实验室规模鼓泡床直接模拟

    8 帖子
    3k 浏览

    Examinination of bubble collisions and coalescence in bubbly flows
    D. Broder and M. Sommerfeld
    还有这一篇,这好像都是小的那个做的,大的具体尺寸我也不知道,不过什么时候能在做实验不好说,实验仪器快被转移了,大的那个。。

  • 25亿网格直接模拟

    1 帖子
    2k 浏览

    Direct numerical simulation of reactor two-phase flows enabled by high performance computing

    捕获.JPG

  • 开放合作系列

    3 帖子
    2k 浏览
    Simulations of gas-liquid dispersions by four-way coupled E-QBMM with aggregation and breakage processes Abstract Results Main References Job: Data processing Test cases debugging and running Writing paper in English if you are willing to Requirements: Benefits: You can get around 30000 yuan bonus from your host (if you write the manuscript). You can use this paper to get promoted faster (if you write the manuscript). Another 10000 yuan from my side If I write the paper, I can support you 40000 yuan as well. However, I dont have much time to write. So I prefer you writing the paper and to be the first author. Project duration
  • 滑移网格算法相关文献求助

    3 帖子
    2k 浏览

    @东岳 谢谢!

  • 颗粒表面化学反应的文献

    已移动
    3 帖子
    2k 浏览

    @ibelief颗粒表面化学反应的文献 中说:

    基于简单碰撞理论煤粉燃烧动力学模型的研究

    非常感谢

  • 7 帖子
    3k 浏览

    moment method对于nucleation和condensation的模拟简单吗?

    即为简单。但是溶解/蒸发数值上要复杂些。2017年我在HZDR的时候同事用fixed pivot method模拟增长,我用moment method,相比来看moment method要更简单,并且不需要处理相空间的drift:https://www.sciencedirect.com/science/article/pii/S1004954118301174/pdfft?md5=84216a8842ce3fd6d1a99d4f327c1ee1&pid=1-s2.0-S1004954118301174-main.pdf

  • 1 帖子
    895 浏览

    带摩擦的运输方程是 u_t = lambda * u_x = eps * u_xx, 通过变化, 我们得到了, 0_1545168405877_3bd01f52-6748-4ece-bac9-e55eeef96e15-image.png ,其中D是导数矩阵, D2 是二阶导数矩阵
    然后, 我们对M进行求特征值(EW), 得到了这种图片, 我想知道这是为什么呢?
    0_1545168538892_dba5eac1-7474-4ce6-b089-77d3d8903c03-image.png

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

    1 帖子
    1k 浏览

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

  • GPU计算: 多相流应用

    3 帖子
    2k 浏览

    DEM加速比很高
    CFD好像一般

  • 如何理解CFD问题只有三维

    3 帖子
    2k 浏览

    二维的vorticity只有一个方向 实际上从涡的形成来看 二维似乎并不适合一些涡量的模拟
    但是对于一些简单的流场 实际上2D具备非常大的优势 他的处理简单 很多算法很好自编程实现
    并且网格可以完全mapped 最大的优势 还是在于后处理的一目了然
    这就是为什么大多数理论研究还是集中2D 但是2D可能不是太适用Fluent Fluent主要接轨工程实际问题

  • 1 帖子
    1k 浏览

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

  • 欧拉双流体模型中应力张量的表示问题

    3 帖子
    2k 浏览

    @东岳 老师,没理解您的意思。后面的$\alpha_q\lambda_q\nabla·\vec{v_q}I$中单位二阶张量是表示体形变率仅出现在法相应力中,红框内等于$-\frac{2}{3}\lambda_q$时才能与法相应力对应啊

  • 二维轴对称动量方程

    3 帖子
    2k 浏览

    谢谢老师,我找出来看一下

  • 4 帖子
    2k 浏览

    0_1543372305704_捕获.JPG

    https://www.ansys.com/-/media/ansys/corporate/resourcelibrary/article/aa-v2-i2-coupling-momentum-and-continuity.pdf

    刚看了下,Fluent那面也只是速度和压力耦合求解。能量抛在了外面。

    但是另一个文档速度压力能量放在了一起 : http://www.bakker.org/dartmouth06/engs150/05-solv.pdf

    速度压力耦合求解在高超那面用的多,如果分析N-S方程的数学特征(双曲抛物),你需要把u, v, w, p当成矢向量,然后求特征值特征向量。接下来求解这个矩阵就是耦合算法。我目前只是了解到思路,没亲自植入过。可

  • 2 帖子
    2k 浏览

    多谢,已更新

  • 2 帖子
    1k 浏览

    5右边多写了个$\rho$吧,谢谢!回北京我更新一下

  • 超声速湍流平板计算问题

    3 帖子
    2k 浏览

    网格加密到Y+<1,计算结果变化不大,对流通量用二阶MUSCL, AUSM格式

  • 关于bubble flow中的通量计算问题

    5 帖子
    3k 浏览

    @东岳 如此说来,除了加密网格,没有好的解决方案啦:zoule: :zoule:

  • 2 帖子
    2k 浏览

    从特征来看,是双曲形方程。是步进形式的问题。

  • 一个小bug调试了一个月

    3 帖子
    2k 浏览

    对,然后不计算边界的话就是0