OF2006 Surface sampling problem : VTP to VTK format

alainislas · January 20, 2021, 12:34

Hello ,

I am using OF2006 release and I want to sample surface data for post-porcessing. Here is my file:

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v2006                                 |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
 
surfaces
{
    type        surfaces;
    libs        (sampling);
    writeControl    timeStep;
    writeInterval    10;

    surfaceFormat    vtk;
    fields       (p U T k epsilon);

    surfaces
    {

        zNormal
        {
            type    cuttingPlane;
            planeType    pointAndNormal;
            pointAndNormalDict
            
            {
                point    (0 0.126 0); // Overrides default basePoint (0 0 0)
                normal       (0 0 1);      // $z: macro for (0 0 1)
                
                interpolate        true;
            }
            
        }


    }
}

Somehow while running the case, the surface files are saved as .vtp and not as .vtk. When I open paraview, I can't read these files. (Also the time doesn't seem to be included). Does anyone know how to overcome this?

olesen · February 9, 2021, 15:39

Quote:

Originally Posted by alainislas

Hello ,

I am using OF2006 release and I want to sample surface data for post-porcessing. Here is my file:

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v2006                                 |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
 
surfaces
{
    type        surfaces;
    libs        (sampling);
    writeControl    timeStep;
    writeInterval    10;

    surfaceFormat    vtk;
    fields       (p U T k epsilon);

    surfaces
    {

        zNormal
        {
            type    cuttingPlane;
            planeType    pointAndNormal;
            pointAndNormalDict
            
            {
                point    (0 0.126 0); // Overrides default basePoint (0 0 0)
                normal       (0 0 1);      // $z: macro for (0 0 1)
                
                interpolate        true;
            }
            
        }


    }
}

Somehow while running the case, the surface files are saved as .vtp and not as .vtk. When I open paraview, I can't read these files. (Also the time doesn't seem to be included). Does anyone know how to overcome this?

This is a feature, not a bug. What you are looking for are the additional format options. Eg,

Code:

formatOptions
{
    vtk
    {
        legacy true;
        // format ascii;
    } 
}

olesen · February 9, 2021, 15:40

If they still cannot be read, it could actually be a bug. Would be useful to know if it is broken in v2012.

frantov · March 2, 2022, 02:27

You can also type in the terminal (or save it on an AllRun file)

foamToVTK -legacy

this will save it as vtk instead of vtp

March 2, 2022, 02:27		#4
frantov Member Francisco T Join Date: Nov 2011 Location: Melbourne, Australia Posts: 64 Blog Entries: 1 Rep Power: 15	You can also type in the terminal (or save it on an AllRun file) foamToVTK -legacy this will save it as vtk instead of vtp alainislas likes this.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
[Gmsh] gmshToFoam generates patches with 0 faces and 0 points	Simurgh	OpenFOAM Meshing & Mesh Conversion	4	August 25, 2023 08:58
OpenFOAM error	Vinay Kumar V	Main CFD Forum	0	February 20, 2020 10:17
[General] Paraview data format conversion from vtk to parallel prog data format.	odho	ParaView	0	September 20, 2016 08:01
Problem with surface monitor position	jf90cl	FLUENT	0	April 28, 2016 12:41
[ICEM] Automatic mesh generation script surface intersection problem	stuart23	ANSYS Meshing & Geometry	0	May 13, 2011 02:10

February 9, 2021, 15:40		#3
olesen Senior Member Mark Olesen Join Date: Mar 2009 Location: https://olesenm.github.io/ Posts: 1,714 Rep Power: 40	If they still cannot be read, it could actually be a bug. Would be useful to know if it is broken in v2012.