Skip to content
  • 最新
  • 版块
  • 东岳流体
  • 随机看[请狂点我]
皮肤
  • Light
  • Cerulean
  • Cosmo
  • Flatly
  • Journal
  • Litera
  • Lumen
  • Lux
  • Materia
  • Minty
  • Morph
  • Pulse
  • Sandstone
  • Simplex
  • Sketchy
  • Spacelab
  • United
  • Yeti
  • Zephyr
  • Dark
  • Cyborg
  • Darkly
  • Quartz
  • Slate
  • Solar
  • Superhero
  • Vapor

  • 默认(不使用皮肤)
  • 不使用皮肤
折叠
CFD中文网

CFD中文网

  1. CFD中文网
  2. Algorithm
  3. 简版NS方程用于CFD有完备性问题?这篇论文有看法

简版NS方程用于CFD有完备性问题?这篇论文有看法

已定时 已固定 已锁定 已移动 Algorithm
18 帖子 2 发布者 8.8k 浏览
  • 从旧到新
  • 从新到旧
  • 最多赞同
回复
  • 在新帖中回复
登录后回复
此主题已被删除。只有拥有主题管理权限的用户可以查看。
  • 数 离线
    数 离线
    数星少年
    在 中回复了 李东岳 最后由 编辑
    #5

    @李东岳 角动量与C成正比,它在时间上的导数就是受到的力矩,这里面写了粘性力矩; 压力应该没有力矩,因为是势能。另外套上了对流项,建立的方程

    1 条回复 最后回复
  • 李东岳李 在线
    李东岳李 在线
    李东岳 管理员
    写于 最后由 编辑
    #6

    我最近随机看论坛,才看到这个帖子。感觉挺有意思。

    文章里面的角动量变量C,是怎么封闭的?

    方程16进行了封闭。

    我在openfoam里面植入一下看看。我初步觉得这个影响太小了基本可以忽略。就类似二次流那种影响。

    http://dyfluid.com/index.html
    需要帮助debug算例的看这个 https://cfd-china.com/topic/8018

    1 条回复 最后回复
  • 李东岳李 在线
    李东岳李 在线
    李东岳 管理员
    写于 最后由 编辑
    #7

    我植入进去了,计算的时候说方程19量纲不一致

    http://dyfluid.com/index.html
    需要帮助debug算例的看这个 https://cfd-china.com/topic/8018

    1 条回复 最后回复
  • 李东岳李 在线
    李东岳李 在线
    李东岳 管理员
    写于 最后由 李东岳 编辑
    #8

    我详细看了一下:

    1. 首先原文的方程2的$\Lambda$和$\Lambda^T$与常规的$\nabla\bfU^T$以及$\nabla\bfU$不一致。这里面牵涉到分子布局与分母布局

    2. 正常来说虽然理论上$\nabla\cdot(\nu\nabla\bfU^T)=0$,但OpenFOAM都进行了求解,并没有直接省略掉

    3. 原文方程4我觉得他的符号写反了,应该是:$(-0.5)... + 0.5... +(-0.5) + 0.5...$,如果这样写的话方程4=0

    这个论文的根本是认为NS方程不满足角动量守恒。我把推导更新到《无痛苦ns方程笔记》第NS方程的角动量守恒一节了。

    http://dyfluid.com/index.html
    需要帮助debug算例的看这个 https://cfd-china.com/topic/8018

    1 条回复 最后回复
  • 数 离线
    数 离线
    数星少年
    写于 最后由 编辑
    #9

    如果流场稳定了,openfoam计算的那项,应该也不起作用了吧。

    1 条回复 最后回复
  • 数 离线
    数 离线
    数星少年
    写于 最后由 编辑
    #10

    原文方程4的力矩计算,是叉乘计算,力矩是有方向的。因此有另外的符号变化,原文是对的。

    1 条回复 最后回复
  • 李东岳李 在线
    李东岳李 在线
    李东岳 管理员
    写于 最后由 李东岳 编辑
    #11

    确实。

    但是那个量纲不一致。c的单位是$\frac{1}{ s}$,C的单位是$\frac{1}{m\cdot s}$,$\nabla\times C$的单位是$\frac{1}{m^2\cdot s}$。$\nabla^2\bfU$的单位是$\frac{1}{m\cdot s}$。

    $\frac{1}{m\cdot s}\neq\frac{1}{m^2\cdot s}$。OpenFOAM计算时候报错了。

    如果流场稳定了,openfoam计算的那项,应该也不起作用了吧。

    这个$\bfU$和$\bfU^T$的贡献一个在切向一个在法向。还约不下去。

    http://dyfluid.com/index.html
    需要帮助debug算例的看这个 https://cfd-china.com/topic/8018

    1 条回复 最后回复
  • 数 离线
    数 离线
    数星少年
    写于 最后由 编辑
    #12

    小c和大C量纲是一致的。它们都是旋转速度啊,量纲都是1/s。C只是c的一个扩展而已,包括了三个方向的旋转角速度分量

    1 条回复 最后回复
  • 李东岳李 在线
    李东岳李 在线
    李东岳 管理员
    写于 最后由 编辑
    #13

    我看他说$C=\nabla c$

    http://dyfluid.com/index.html
    需要帮助debug算例的看这个 https://cfd-china.com/topic/8018

    1 条回复 最后回复
  • 数 离线
    数 离线
    数星少年
    写于 最后由 编辑
    #14

    啊?大C不是公式19下面那行文字定义的吗?

    1 条回复 最后回复
  • 李东岳李 在线
    李东岳李 在线
    李东岳 管理员
    写于 最后由 编辑
    #15

    是的就是那个。原文是

    $C=[c_x, c_y, c_z]$

    那就是$C=\nabla c$。 还是说我理解错了

    http://dyfluid.com/index.html
    需要帮助debug算例的看这个 https://cfd-china.com/topic/8018

    1 条回复 最后回复
  • 数 离线
    数 离线
    数星少年
    写于 最后由 编辑
    #16

    你理解错了,xyz只是一个下标而已。如果是二维情况,就只有cz,表示c是垂直于xy平面(平行于z轴)的矢量。其余两个是三位的类似情况而已

    1 条回复 最后回复
  • 李东岳李 在线
    李东岳李 在线
    李东岳 管理员
    写于 最后由 编辑
    #17

    这样啊!那我试试

    http://dyfluid.com/index.html
    需要帮助debug算例的看这个 https://cfd-china.com/topic/8018

    1 条回复 最后回复
  • 数 离线
    数 离线
    数星少年
    写于 最后由 编辑
    #18

    三维的类似情况,我最后一个打错字了。

    1 条回复 最后回复

  • 登录

  • 登录或注册以进行搜索。
  • 第一个帖子
    最后一个帖子
0
  • 最新
  • 版块
  • 东岳流体
  • 随机看[请狂点我]