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. OpenFOAM
  3. buoyantPimpleFoam相关问题

buoyantPimpleFoam相关问题

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

    @李东岳 http://dyfluid.com/buoyantPimpleFoam.html
    在这个链接里

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

    那个是旧的链接,写的不完善。目前我得网页里面不提供公开访问权限,你找到的应该是搜到的,我先删掉了。你看看新版本的。

    http://dyfluid.com/compsimple.html

    http://dyfluid.com/compiso.html

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

    Y 1 条回复 最后回复
  • Y 离线
    Y 离线
    yjc95
    在 中回复了 李东岳 最后由 编辑
    #10

    @李东岳 谢谢东子老师

    1 条回复 最后回复
  • Y 离线
    Y 离线
    yjc95
    写于 最后由 编辑
    #11

    @李东岳 东子老师,我把可压+稳态算法仔仔细细的看完了,大体上看懂了,但有一些疑问向您请教。您上次说我排版不利于阅读,这次我专门看了您提供的帮助里面的教程。

    问题1

    您在文章里说:上标n表示为当前迭代步(已知),上标表示下一个迭代步(待求),n+1也表示下一个迭代步。在稳态条件下,5b2efa88-9199-4b80-b4f0-c0c886fa0089-CodeCogsEqn.svg 应该等于658840c8-71e8-4642-9a77-9260eaa6b6f6-CodeCogsEqn (1).svg 吧。方程(18)-(14)等号左边不应该是0吗?为啥是d196288e-c3a2-4ce2-8580-b309954568ec-CodeCogsEqn (2).svg 呢?上标、**、n与n+1是啥关系呀?我猜测*、**是不是n到n+1之间的迭代步?

    问题2

    方程(23)没有c9eb4ec2-5110-41b3-abda-7fb63743d272-CodeCogsEqn (3).svg 呀?另外,为啥是15d26edd-0fbf-4db7-8143-5855742ca1e3-CodeCogsEqn (4).svg 而不是8e7d2219-1a52-48f2-8bb4-6a26b87935c8-CodeCogsEqn (5).svg 或者e88aae6e-7526-4d34-adc1-c5e70ac64eb3-CodeCogsEqn (6).svg 呢?
    e1cfe02f-da35-4150-9e43-c2a4356a47b5-image.png

    问题3

    由方程(3)可得670c6b75-18cf-4c6d-a39e-0a4f5d9bb832-CodeCogsEqn (7).svg ,又7854ecd4-8f21-40f4-8fd2-c5b73b97ff6b-CodeCogsEqn (8).svg ,所以 7fb1a991-7f26-4b14-a90e-6debc0afa22c-CodeCogsEqn (9).svg
    但由方程(27)得到的却是图中红框部分,如何得出方程(28)中蓝框部分的呢?方程(36)是我算的结果,方程(42)又是红框的结果了。
    9fc6e84b-b77d-45a3-af6b-94b05b4b0a23-image.png

    问题4

    方程(19)由(18)-(14)得来,1edcea56-096d-4193-84f5-b196f809f1b8-CodeCogsEqn (10).svg ,而图中这个推导过程:
    5cf4e3dd-7af1-4218-96d6-0d761db6ec75-image.png是 48040ba7-9ec7-4fcb-b8db-caad5cf8c058-CodeCogsEqn (11).svg ,是不是矛盾了呢?

    同样的,方程(43)中,0ef0f602-2a76-4be8-8342-0f7b005947b9-CodeCogsEqn (12).svg ;
    方程(25)中,dcb61c86-7441-4b63-8724-e1b4189113e2-CodeCogsEqn (13).svg
    我被这些上标n、'、*搞懵了:135:

    问题5

    方程(48)中,红框部分是不是少了一个c4cc7b73-f09a-4056-b44f-16484a4ee9c5-CodeCogsEqn (14).svg ?
    72a9548f-2d79-4c08-8db7-d50ec6ff9dbb-image.png

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

    我昨天凑巧更新了一下,http://dyfluid.com/compsimple.html 你再看看

    打公式的时候,用美元符号最好 $\alpha$

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

    Y 1 条回复 最后回复
  • Y 离线
    Y 离线
    yjc95
    在 中回复了 李东岳 最后由 编辑
    #13

    @李东岳 谢谢东子老师

    1 条回复 最后回复
  • Y 离线
    Y 离线
    yjc95
    写于 最后由 编辑
    #14

    @李东岳 东子老师,我看完了您更新的文章,还有几个小问题。另外,有的时候我用 $\alpha$打不出来我想要的公式,可能前后公式会不统一,请见谅。

    问题1

    方程(4)为啥把速度处理成 $\textbf{U}^{n}\textbf{U}^{\ast}$呢?为啥不是$\textbf{U}^{n}\textbf{U}^{n}$或者$\textbf{U}^{\ast}\textbf{U}^{\ast}$呢?

    问题2

    方程(25)密度为啥写成$\rho_f^{n}+\rho_f^{'}$,为啥要修正呢?速度为啥用$\textbf{U}_f^{\ast\ast}$呢?用$\textbf{U}_f^{\ast}$不行吗?

    问题3

    方程(15)处,您说更新热物理模型库中的密度为$\rho^{\ast}=p^{n}/\psi^{n}$,这和方程(3)不一致呀?方程(43)处您又说是根据状态方程更新密度中间量。俩方程长的一样,一个是热物理模型库,一个是状态方程?

    问题4

    方程(42)中,您是用$\rho_{f}^{'}$来替代方程(41)中的$\bigtriangledown \rho$吗?

    问题5

    方程(48)写成(49),等号右边这一项为啥可以用924a8891-09b5-41fa-b8be-04dcf634a439-CodeCogsEqn.svg
    来代替$\bigtriangledown \rho^{*}_{rhg}$呢?

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

    问题1这是要做线性化。

    问题3,热物理模型库就是状态方程。$\psi=1/R/T$

    问题4,一个是离散形式,一个是连续形式。离散形式的是修正量。

    问题5,一个是离散形式,一个是连续形式。

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

    Y 1 条回复 最后回复
  • Y 离线
    Y 离线
    yjc95
    在 中回复了 李东岳 最后由 编辑
    #16

    @李东岳
    问题3,状态方程不是方程(3)$p=\rho RT$吗?$\psi=1/(RT)$那么不应该$\rho =p/(RT)=p\psi$吗?为什么是$\rho =p/\psi$呢?

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

    dyfluid上面应该写错了,明天我核实一下告诉你

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

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

    rho=p \psi 这个是对的,已更正

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

    Y 1 条回复 最后回复
  • Y 离线
    Y 离线
    yjc95
    在 中回复了 李东岳 最后由 编辑
    #19

    @李东岳 老师,好像没改过来:141: a5e72a24-fb56-445d-9530-cc02aa879351-image.png

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

    你刷新一下看看

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

    Y 1 条回复 最后回复
  • Y 离线
    Y 离线
    yjc95
    在 中回复了 李东岳 最后由 编辑
    #21

    @李东岳 刷新了,还是这样

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

    我想起了 ,我忘记上传了... 现在上传了

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

    1 条回复 最后回复
  • Y 离线
    Y 离线
    yjc95
    写于 最后由 编辑
    #23

    @李东岳 老师,前面算法理论部分我看懂了,但我在看关键代码这一部分的时候,发现和前面公式有点对应不上。比如高速流transonic中的压力方程通过下面的代码组建:
    9eb3bb56-48e1-4201-b258-a34bc21edf72-image.png
    fvc、fvm是啥意思啊?我查了下无痛苦N-S方程笔记
    50fbd12b-3ffe-400a-95d8-27acf94c2030-image.png
    为啥这里要除以$\bigtriangleup V\bigtriangleup t $ 呢?
    另外,$fvc::div(phiHbyA)$ 是对应方程中的哪一项呢?代码中每一项代表方程中的哪一项呢?这个我不是很理解。

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

    fvc是显性离散,fvm是隐性离散。

    除掉时间和体积就是因为代码本身就是这样做的

    divPhiHbyA是方程的第一项

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

    Y 1 条回复 最后回复
  • Y 离线
    Y 离线
    yjc95
    在 中回复了 李东岳 最后由 编辑
    #25

    @李东岳 代码等号右边那项是啥呀?

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

    可以忽略

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

    1 条回复 最后回复
  • Y 离线
    Y 离线
    yjc95
    写于 最后由 编辑
    #27

    @李东岳 东子老师,我看了可压+稳态的算法之后,紧接着看了可压+瞬态的算法,有几个问题想请教您。

    问题1

    方程(8)等号左边第一项离散后是$ \int (\int \frac{\partial \rho }{\partial t}dt)dV=\bigtriangleup V(\rho ^{t+\bigtriangleup t}-\rho ^{t})$
    第二项离散后是
    $\int (\int \bigtriangledown \cdot \rho \textbf{U}dV)dt=$ $\bigtriangleup t$ $\sum_{f} \rho _{f}^{t+\bigtriangleup t}$ $\textbf{U} _{f}^{t+\bigtriangleup t}$ $\cdot S$
    这俩公式相加化简后应该是48f2c6bc-5e7b-4a12-96f0-40a72b750163-CodeCogsEqn (1).svg
    为啥方程(9)左边第一项没有$\bigtriangleup V$呢?而方程(11)又有$V{p}$呢?

    问题2

    方程(19)左边第一项应该是bc3c7a4e-5360-4f16-a008-aecfd09dad9a-CodeCogsEqn (1).svg吧?
    另外,今晚回复帖子打字好卡:135:

    1 条回复 最后回复

  • 登录

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