Skip to content

Algorithm

基本理论区

408 主题 2.1k 帖子
  • 欧拉双流体模型中应力张量的表示问题

    3
    3 帖子
    3k 浏览
    AlexA

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

  • 二维轴对称动量方程

    3
    3 帖子
    3k 浏览
    AlexA

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

  • 4 帖子
    4k 浏览
    李东岳

    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 帖子
    2k 浏览
    李东岳

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

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

    3
    3 帖子
    3k 浏览

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

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

    5
    5 帖子
    4k 浏览
    yhdthuY

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

  • 连续性方程是双曲型方程,还是椭圆型方程

    2
    2 帖子
    2k 浏览
    李东岳

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

  • 一个小bug调试了一个月

    3
    3 帖子
    3k 浏览
    李东岳

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

  • 求助贴!关于无稳态流场的CFD计算

    5
    5 帖子
    4k 浏览

    @him 谢谢解答,另外一般商业软件的定常求解器是用的时间推进来求解的吗?如果是时间推进的定常求解器,这种不稳定的振荡可否理解为流场随时间的变化?

  • LES和RANS中,离散的对象是相同的吗?

    31
    31 帖子
    37k 浏览
    H

    如果使用standard KE,计算二维流场,在某个高雷诺数下把Cd, Cl, St都算对,那么说明碰对了,意思就是好运。这并不能说明standard KE可以算对漩涡脱离。

  • 21 帖子
    20k 浏览
    H

    但是RC插值没有离散,因为它是显式计算源项 -- “速度插值+高阶压力导数项”。

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

    1
    1 帖子
    2k 浏览

    基于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)能够告知粗细网格过渡区域,造成的流场震荡的原因可能有哪些;

    再次感谢大家!

  • OpenFOAM里面有热通量heat flux的边界条件吗?

    4
    4 帖子
    5k 浏览

    @东岳 这个边界条件已经整合到exteranlHeatFluxTemperature里面了,算的也是温度梯度,也就是$q=-\kappa \frac{\partial T}{\partial n}$,这种边界按照FDS技术手册里面应该是DNS的计算方法。
    LES计算时给定的方法是$q=h \Delta T$,$h$是对流换热系数。openfoam里面的边界条件要么是固定值要么是梯度,这种差值不知道在of里面能不能实现

  • 关于气液分离的采出分率问题

    2
    2 帖子
    2k 浏览
    C

    补充一下,不知道这些图中的这些数据点是如何提取的,谢谢大家

  • 均匀各项同性湍流的初场设置问题

    4
    4 帖子
    4k 浏览
    B

    @shifc 不知道你现在解决了没有,我最近做格栅湍流也碰到了同样的问题,我的解决方法是谱空间的复速度对三个坐标轴必须都是共轭对称的。

  • Hill vortex是怎么产生的?

    2
    2 帖子
    2k 浏览
    S

    你好👋

  • 数值传热学P94,关于内节点法的离散过程。

    4
    4 帖子
    4k 浏览
    李东岳

    额,自问自答?

  • 计算对流扩散方程没有扩散效果?

    8
    8 帖子
    7k 浏览
    李东岳

    不好意思,那个算例里面求解器我这面编译不了(openfoam-4.x and 5.x).. 估计后来问题已经解决了吧

  • 压力方程和解

    2
    2 帖子
    3k 浏览

    这个是理论大坑啊。

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