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. 关于《icoFOAM解析》的一些疑问

关于《icoFOAM解析》的一些疑问

已定时 已固定 已锁定 已移动 Algorithm
4 帖子 2 发布者 3.5k 浏览
  • 从旧到新
  • 从新到旧
  • 最多赞同
回复
  • 在新帖中回复
登录后回复
此主题已被删除。只有拥有主题管理权限的用户可以查看。
  • 秦 离线
    秦 离线
    秦淮冷月
    写于 最后由 编辑
    #1

    以公式(19)作为理解的主体,操作过程为:(19)扣掉(18)->运用忽略临点影响假设,(19)再扣掉一项->(19)再次加上(18),具体过程如下:

    1.公式(19)扣掉(18)推得(20),

    \begin{equation}
    A_P \underbrace { \left( \vec U_P^{n+1} - \vec U_P^r \right) }_{U'_P}
    +\sum A_N \overbrace{\left( \vec U^{n+1}-\vec U^{n}\right) }^{U'_N}
    = -\sum\overbrace{\left(p_f^{n+1} - p_f^n \right)}^{p'_f}\vec S_f
    \end{equation}

    2.引入忽略临点假设后,LHS第二项被消除,有

    \begin{equation}
    A_P \overbrace { \left( \vec U_P^{n+1} - \vec U_P^r \right) }^{U'_P}
    =-\sum \overbrace{\left( p_f^{n+1} - p_f^n \right)}^{p'_f} \vec S_f
    \end{equation}

    3.再加上(18),下式中标绿的部分为(18)

    \begin{equation}
    A_P \overbrace { \left( \vec U_P^{n+1} - \vec U_P^r \right) }^{U'_P} + \color{green} {A_P \vec U_P^r}
    +\color{green}{\sum A_N \vec U_N^r}
    =-\sum \overbrace{\left( p_f^{n+1} - p_f^n \right)}^{p'_f} \vec S_f \color{green}{ -\sum p_f^n \vec S_f
    +S_P^n} \
    \rightarrow\
    A_P \overbrace { \vec U_P^{n+1} }^{\color{red}{U^*_P}}
    +\color{green}{\sum A_N \vec U_N^r}
    =-\sum \underbrace{ p_f^{n+1} }_{\color{red}{p^{*}_f}} \vec S_f
    +\color{green}{S_P^n}
    \end{equation}

    对比操作前后的变化,
    操作前,
    \begin{equation*}
    A_P \vec U_P^{n+1}
    +\color{green} {\sum A_N \vec U_N^{n+1}}
    =S_P^n
    -\sum p_f^{n+1} \vec S_f
    \end{equation*}
    操作后,
    \begin{equation*}
    A_P \underbrace { \vec U_P^{n+1} }_{\color{red}{U^{*}_P}}
    +\color{green}{\sum A_N \vec U_N^r}
    =S_P^n-\sum \underbrace{ p_f^{n+1} }_{\color{red}{p^*_f}} \vec S_f
    \end{equation*}

    结论:1.带星的符号和带撇的符号,是同一个量;
    2.上述操作主要是用预测量(上标$r$)替换了最终收敛量(上标$n+1$)

    上述理解是否正确?

    1 条回复 最后回复
  • 秦 离线
    秦 离线
    秦淮冷月
    写于 最后由 编辑
    #2

    尴尬....结论1应该去掉。但是不能编辑帖子

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

    你的方程是正确的。不过要注意这个是用来组建连续性方程的。

    我在想有空把$r$之类改成$t+\dt$

    你的排版挺好的

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

    秦 1 条回复 最后回复
  • 秦 离线
    秦 离线
    秦淮冷月
    在 中回复了 李东岳 最后由 编辑
    #4

    @东岳 感谢回复。建议原文适当增加配图,这将提高可阅读性。

    1 条回复 最后回复

  • 登录

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