CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Visualization & Post-Processing Software > ParaView

[OpenFOAM] PlotOverLine gives NaN values on surface sample

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 23, 2024, 09:47
Default PlotOverLine gives NaN values on surface sample
  #1
New Member
 
Ege Batmaz
Join Date: Dec 2018
Location: Germany
Posts: 13
Rep Power: 7
egebat7 is on a distinguished road
Hi,



I am facing a small issue when I plot my 2D data over a line. I used this function to create my VTK format surface data:


Code:
    zSlice1
    {
        type    surfaces;
        libs ( "libsampling.so" );
        writeControl    runTime;
        writeInterval   0.01;
        format          binary;
        fields          ( U T p_vapor k epsilon nut alphat rhok );
        interpolationScheme cellPoint;
        surfaceFormat   vtk;
        
        surfaces
        (
            zC1
            {
                type    cuttingPlane;
                planeType   pointAndNormal;
                pointAndNormalDict
                {
                    point   ( 0 0 0 );
                    normal  ( 0 0 1 );
                }
                interpolate true;
            }
        );
    }
When I try to visualize this surface it looks well. However, when I try to plot over a line using "PlotOverLine" feature in Paraview, some values don't appear (the line is discontinues in some parts) and give NaN.

Is this a bug or a numerical situation where some data is not exactly at z = 0 plane? You can see how the plot looks like from the attachment. For every field (T, U, p ...), same locations are missing. However, missing locations change when I plot a different line, but always with holes in it.


Best regards,
Ege


Edit: The missing points are not related with the local mesh density, since it is a very uniform mesh and contains only block cells.
Attached Images
File Type: jpg Paraview_screenshot.jpg (77.2 KB, 4 views)
egebat7 is offline   Reply With Quote

Old   March 24, 2024, 06:56
Default
  #2
New Member
 
Ege Batmaz
Join Date: Dec 2018
Location: Germany
Posts: 13
Rep Power: 7
egebat7 is on a distinguished road
Okay I solved the issue by changing my cutting plane z-axis point from 0 to 0.001:


Code:
            zC1
            {
                type    cuttingPlane;
                planeType   pointAndNormal;
                pointAndNormalDict
                {
                    point   ( 0 0 0.001 );
                    normal  ( 0 0 1 );
                }
                interpolate true;
            }
I am working with hexahedral cells and the plane passing through z = 0 is always intersects with the cell edges. Maybe because of this, sampling function was not working properly. I am still not sure but this could be the issue.
egebat7 is offline   Reply With Quote

Reply


Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
icoFoam, pisoFoam, pimpleFoam unsuccessful transient solution of laminar flow channel Andrew Gow OpenFOAM Running, Solving & CFD 0 December 24, 2022 11:13
problem during mpi in server: expected Scalar, found on line 0 the word 'nan' muth OpenFOAM Running, Solving & CFD 3 August 27, 2018 04:18
thermoPhysical: simulation crashes with perfectGas-rhoCoeffs piu58 OpenFOAM Running, Solving & CFD 1 December 24, 2017 10:42
[CFD-Post] velocity values from user surface iltis CFX 1 January 11, 2017 13:16
sample: free surface & more Phicau OpenFOAM Post-Processing 2 October 11, 2011 10:19


All times are GMT -4. The time now is 09:21.