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. fvOptions中添加的源项需要考虑网格的体积吗?

fvOptions中添加的源项需要考虑网格的体积吗?

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

    大家好!

    最近我在OpenFOAM中使用fvOptions设置了一个动量方程的源项,定义为
    \begin{equation}
    S=-C_z \times |U| \times U
    \end{equation}
    $C_z$是一个常数,$|U|$是速度的大小,$U$是速度,我在fvOptions中是直接写了

    eqn += -Cz*mag(U)
    

    我看到源项的表达式应该是写成,是乘以网格体积了的
    \begin{equation}
    S=S_U \times V_P + S_P \times\ V_P \times T_P
    \end{equation}
    请问我在fvOptions中写的方式是需要乘以网格的体积吗

    李东岳李 1 条回复 最后回复
  • 李东岳李 在线
    李东岳李 在线
    李东岳 管理员
    在 中回复了 allanZhong 最后由 李东岳 编辑
    #2

    @allanZhong 2年过去了。。。是需要乘以体积。例如我刚写了个科氏力源项:

    nergySource
    {
        type            coded;
        selectionMode   all;
        field           U;
    
        codeInclude
        #{
        #};
    
        codeAddSup
        #{
            Pout<< "**codeAddSup**" << endl;
            const volVectorField& U = eqn.psi();
            vectorField& USource = eqn.source();
            vector omega(0, 0, 5e-5);
            const DimensionedField<scalar, volMesh>& V = mesh().V();
            USource += (2.0*omega)^U*V;
        #};
    
        codeAddRhoSup
        #{
            Pout<< "**codeAddRhoSup**" << endl;
        #};
    
        codeAddAlphaRhoSup
        #{
            Pout<< "**codeAddAlphaRhoSup**" << endl;
        #};
    }
    
    

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

    1 条回复 最后回复

  • 登录

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