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. 如何修改求解器计算速度分量的导数

如何修改求解器计算速度分量的导数

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

    大家好,我想在of中输出每个时间步下速度分量的导数,比如输出$\frac{\partial U(x)}{\partial x}$ ,$\frac{\partial U(y)}{\partial y}$。如何在求解器里面实现?我找到了标量的写法,速度矢量不知道如何处理。

            volVectorField gradT(fvc::grad(T));
    
            volScalarField gradTx
            (
                IOobject
                (
                    "gradTx",
                    runTime.timeName(),
                    mesh,
                    IOobject::NO_READ,
                    IOobject::AUTO_WRITE
                ),
                gradT.component(vector::X)
            );
    
            runTime.write();
        }
    
    1 条回复 最后回复
  • 李东岳李 在线
    李东岳李 在线
    李东岳 管理员
    写于 最后由 编辑
    #2

    fvc::grad(U.component(0))

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

    chengan.wangC 2 条回复 最后回复
  • chengan.wangC 离线
    chengan.wangC 离线
    chengan.wang
    在 中回复了 李东岳 最后由 编辑
    #3

    @东岳 谢谢东岳大神 :140:

    1 条回复 最后回复
  • chengan.wangC 离线
    chengan.wangC 离线
    chengan.wang
    在 中回复了 李东岳 最后由 编辑
    #4

    @东岳
    李老师,

        if (runTime.outputTime())
        {
            volVectorField gradUx(fvc::grad(U.component(0)));
            volVectorField gradUz(fvc::grad(U.component(2)));
    
            volScalarField dUxdx_add_dUzdz
            (
                IOobject
                (
                    "dUxdx_add_dUzdz",
                    runTime.timeName(),
                    mesh,
                    IOobject::NO_READ,
                    IOobject::AUTO_WRITE
                ),
                gradUx.component(vector::X)+gradUz.component(vector::Z)
            );
    
            runTime.write();
        }
    

    U.component(0)指的是vx,U.component(2)指的是vz
    我想问一下 gradUx.component(vector::X)中的X或者gradUz.component(vector::Z)中的Z分别指的是什么?是否与坐标轴一致?

    1 条回复 最后回复

  • 登录

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