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 获取每个点的坐标向量

OpenFoam 获取每个点的坐标向量

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

    各位朋友:

    我想通过计算得到U在切向和轴向的速度分量,所以要用到如下的公式:

    Uradial = vector(position) & U / mag( vector(position))
    

    读取U场数据用如下代码:

    IOobject Uheader
    (
        "U"
        runTime.timeName(),
        mesh,
        IOobject::MUST_READ
    );
    

    可是在读取各个网格中心坐标的时候,遇到了问题,我的代码如下:

    void Foam::calc(const argList& args, const Time& runTime, const fvMesh& mesh)
    {
        //bool writeResults = !args.optionFound("noWrite");
    
        IOobject Uheader
        (
            "U",
            runTime.timeName(),
            mesh,
            IOobject::MUST_READ
        );
        // Check U exists
        if (Uheader.headerOk() )
        {
            Info<< "    Reading U" << endl;
            volVectorField U(Uheader, mesh);
          
           // 读取各个网格中心点的坐标
            volVectorField&  centres =  mesh.C().internalField();
          // 由于mesh.C()在网格内部返回的是中心点的坐标,在网格边界处返回的是面心坐标,所以用.internalField(),可是在这儿报错
    
    
            volScalarField URadial
            (
                IOobject
                (
                    "URadial",
                    runTime.timeName(),
                    mesh
                ),
                ((centres - centres.component(2) ) & U) //mag ((mesh.C() - (0., 0., mesh.C().component(2)))
                //这个地方我想得到一个平行于z轴的坐标,方法是用这一点的位置向量减去(0, 0, z)可是系统报错告诉我不能这样减。
            );
        }
        else
        {
            Info<< "    Missing p or T" << endl;
        }
        Info<< "\nEnd\n" << endl;
    }
    

    我想问一下,如果想实现mesh.C()- (0, 0, z)的话,需要怎么做?

    非常感谢!
    Janry

    Turbomachinery Enginee

    1 条回复 最后回复

  • 登录

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