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. 均匀各项同性湍流的初场设置问题

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

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

    =========
    最近在学《可压缩湍流直接数值模拟》一书的过程中,学生我针对7.2.2节各向同性湍流初场设置的内容有几点疑惑,想请教一下各位老师/前辈。
    !==============================分割线====================================
    背景介绍:
    根据书中7.2.2节第一部分的描述,公式(7.2.3)给出了谱空间中的速度谱${\hat v_i}$为:

    \begin{equation}
    {\hat v_i}\left( {{k_1},{k_2},{k_3}} \right) = \sqrt {{\textstyle{2 \over 3}}E\left( {{k_1},{k_2},{k_3}} \right)} {e^{{\rm{\hat i}}\theta }}
    \end{equation}

    将三维能谱${E(k1,k2,k3)}$代入上式中,并对${\hat v_i}$进行投影,得到满足散度为0的速度谱${\hat u_i}$:

    \begin{equation}
    {\hat u_i} = \left( {\frac{{{k_i}{k_j}}}{{{k^2}}} - {\delta _{ij}}} \right){\hat v_j}
    \end{equation}

    此时,由于(1)中包含有随机相位${\theta}$导致速度谱${\hat u_i}$中包含有虚数部分。我针对这样的速度谱${\hat u_i}$直接进行Fourier逆变换:

    \begin{equation}
    {u_i}\left( {{x_1},{x_2},{x_3}} \right) = {\left( {{\textstyle{1 \over {2\pi }}}} \right)^3}\int_{ - \infty }^\infty {\int_{ - \infty }^\infty {\int_{ - \infty }^\infty {{{\hat u}_i}\left( {{k_1},{k_2},{k_3}} \right){\rm{d}}{k_1}{\rm{d}}{k_2}{\rm{d}}{k_3}} } }
    \end{equation}

    所得到的物理空间域内的速度${u_i}$存在虚部。

    !==============================分割线====================================
    疑惑点:
    所以我的疑惑是:

    1. 虚部在实际物理问题中意义不明,我利用上述方法所得到的物理空间域内速度${u_i}$的虚部不为0,这是否意味着我针对速度谱${\hat u_i}$直接进行Fourier逆变换(3)是一种错误的操作?
    2. 如果不能直接针对速度谱${\hat u_i}$直接进行Fourier逆变换,那么我应该采用怎样的操作才能得到物理空间域内速度${u_i}$呢?
    3. 速度谱${\hat v_i}$通过公式(2)的投影产生新的速度谱 ,两者模之间的大小并不相同。这会不会导致${\hat u_i}$不满足原始的能谱${ E(k1,k2,k3) }$?
    1 条回复 最后回复
  • 李东岳李 离线
    李东岳李 离线
    李东岳 管理员
    写于 最后由 编辑
    #2

    很遗憾,完全不懂,目前还没研究过有关虚部的东西.. :big_mouth: :big_mouth:

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

    S 1 条回复 最后回复
  • S 离线
    S 离线
    ShiFC
    在 中回复了 李东岳 最后由 编辑
    #3

    @李东岳 哈哈哈,我也是来咨询一下的:tongue_out:

    1 条回复 最后回复
  • B 离线
    B 离线
    Briareos
    写于 最后由 编辑
    #4

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

    1 条回复 最后回复

  • 登录

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