timeVaryingMappedFixedValue error

mcfdma · May 5, 2020, 08:17

Hi. I am working on a 2D axi case where I want to create a parabolic velocity profile at the inlet. I am trying to use the tutorial case pitzDailyExptInlet as a reference for the timeVaryingMappedFixedValue aspect.

My U file

Code:

boundaryField
{
    axis
    {
        type            empty;
    }
    inlet
    {
        type            timeVaryingMappedFixedValue;
        offset          (0 0 0);
        setAverage      off;
    }
    wall
    {
        type            fixedValue;
	value		uniform (0 0 0);
    }
    outlet
    {
        type            zeroGradient;
    }
    front
    {
        type            wedge;
    }
    back
    {
        type            wedge;
    }
}

I am including the blockMeshDict to provide a clear understanding.

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.1.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version         2.0;
    format          ascii;
 
    root            "";
    case            "";
    instance        "";
    local           "";
 
    class           dictionary;
    object          blockMeshDict;
}
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
convertToMeters 0.1;

vertices
(
	(0 0 0)			//0
	(20 0 0)		//1			
	(0 1 -0.04366094291)	//2		
	(20 1 -0.04366094291)	//3
	(0 1 0.04366094291)	//4
	(20 1 0.04366094291)	//5

);
 
blocks
(
    hex (0 1 3 2 0 1 5 4) (100 100 1) simpleGrading (1.0 1.0 1.0)  //1

);

edges
(
);

boundary
(
    axis
    {
        type empty;
        faces
        (
            (0 1 1 0)
        );
    }

    inlet
    {
        type patch;
        faces
        (
            (0 0 4 2)
        );
    }
    
    wall
    {
        type wall;
        faces
        (
            (2 4 5 3)

        );
    }

    outlet
    {
        type patch;
        faces
        (
            (1 3 5 1)
        );
    }

    front
    {
        type wedge;
        faces
        (
            (0 1 5 4)

        );
    }

    back
    {
        type wedge;
        faces
        (
            (0 2 3 1)

        );
    }
);

mergePatchPairs
(
);

// ************************************************************************* //

I have generated a points file (constant/boundaryData/inlet) for the inlet based on the tutorial case.

Code:

// Points
102
(

// min z
(0.00   0.0000 -0.000000000000) 
(0.00   0.0020 -0.000087321886) 
(0.00   0.0040 -0.000174643772) 
(0.00   0.0060 -0.000261965657) 
(0.00   0.0080 -0.000349287543) 
(0.00   0.0100 -0.000436609429) 
(0.00   0.0120 -0.000523931315) 
(0.00   0.0140 -0.000611253201) 
(0.00   0.0160 -0.000698575087) 
(0.00   0.0180 -0.000785896972) 
(0.00   0.0200 -0.000873218858) 
(0.00   0.0220 -0.000960540744) 
(0.00   0.0240 -0.001047862630) 
(0.00   0.0260 -0.001135184516) 
(0.00   0.0280 -0.001222506401) 
(0.00   0.0300 -0.001309828287) 
(0.00   0.0320 -0.001397150173) 
(0.00   0.0340 -0.001484472059) 
(0.00   0.0360 -0.001571793945) 
(0.00   0.0380 -0.001659115831) 
(0.00   0.0400 -0.001746437716) 
(0.00   0.0420 -0.001833759602) 
(0.00   0.0440 -0.001921081488) 
(0.00   0.0460 -0.002008403374) 
(0.00   0.0480 -0.002095725260) 
(0.00   0.0500 -0.002183047145) 
(0.00   0.0520 -0.002270369031) 
(0.00   0.0540 -0.002357690917) 
(0.00   0.0560 -0.002445012803) 
(0.00   0.0580 -0.002532334689) 
(0.00   0.0600 -0.002619656575) 
(0.00   0.0620 -0.002706978460) 
(0.00   0.0640 -0.002794300346) 
(0.00   0.0660 -0.002881622232) 
(0.00   0.0680 -0.002968944118) 
(0.00   0.0700 -0.003056266004) 
(0.00   0.0720 -0.003143587889) 
(0.00   0.0740 -0.003230909775) 
(0.00   0.0760 -0.003318231661) 
(0.00   0.0780 -0.003405553547) 
(0.00   0.0800 -0.003492875433) 
(0.00   0.0820 -0.003580197318) 
(0.00   0.0840 -0.003667519204) 
(0.00   0.0860 -0.003754841090) 
(0.00   0.0880 -0.003842162976) 
(0.00   0.0900 -0.003929484862) 
(0.00   0.0920 -0.004016806748) 
(0.00   0.0940 -0.004104128633) 
(0.00   0.0960 -0.004191450519) 
(0.00   0.0980 -0.004278772405) 
(0.00   0.1000 -0.004366094291) 

// max z
(0.00   0.0000 0.000000000000) 
(0.00   0.0020 0.000087321886) 
(0.00   0.0040 0.000174643772) 
(0.00   0.0060 0.000261965657) 
(0.00   0.0080 0.000349287543) 
(0.00   0.0100 0.000436609429) 
(0.00   0.0120 0.000523931315) 
(0.00   0.0140 0.000611253201) 
(0.00   0.0160 0.000698575087) 
(0.00   0.0180 0.000785896972) 
(0.00   0.0200 0.000873218858) 
(0.00   0.0220 0.000960540744) 
(0.00   0.0240 0.001047862630) 
(0.00   0.0260 0.001135184516) 
(0.00   0.0280 0.001222506401) 
(0.00   0.0300 0.001309828287) 
(0.00   0.0320 0.001397150173) 
(0.00   0.0340 0.001484472059) 
(0.00   0.0360 0.001571793945) 
(0.00   0.0380 0.001659115831) 
(0.00   0.0400 0.001746437716) 
(0.00   0.0420 0.001833759602) 
(0.00   0.0440 0.001921081488) 
(0.00   0.0460 0.002008403374) 
(0.00   0.0480 0.002095725260) 
(0.00   0.0500 0.002183047145) 
(0.00   0.0520 0.002270369031) 
(0.00   0.0540 0.002357690917) 
(0.00   0.0560 0.002445012803) 
(0.00   0.0580 0.002532334689) 
(0.00   0.0600 0.002619656575) 
(0.00   0.0620 0.002706978460) 
(0.00   0.0640 0.002794300346) 
(0.00   0.0660 0.002881622232) 
(0.00   0.0680 0.002968944118) 
(0.00   0.0700 0.003056266004) 
(0.00   0.0720 0.003143587889) 
(0.00   0.0740 0.003230909775) 
(0.00   0.0760 0.003318231661) 
(0.00   0.0780 0.003405553547) 
(0.00   0.0800 0.003492875433) 
(0.00   0.0820 0.003580197318) 
(0.00   0.0840 0.003667519204) 
(0.00   0.0860 0.003754841090) 
(0.00   0.0880 0.003842162976) 
(0.00   0.0900 0.003929484862) 
(0.00   0.0920 0.004016806748) 
(0.00   0.0940 0.004104128633) 
(0.00   0.0960 0.004191450519) 
(0.00   0.0980 0.004278772405) 
(0.00   0.1000 0.004366094291) 

)

// ************************************************************************* //

The U file (constant/boundaryData/inlet/0)

Code:

// Data on points
102
(

//minz
(  1.0000   0.0000   0.0000) 
(  0.9996   0.0000   0.0000) 
(  0.9984   0.0000   0.0000) 
(  0.9964   0.0000   0.0000) 
(  0.9936   0.0000   0.0000) 
(  0.9900   0.0000   0.0000) 
(  0.9856   0.0000   0.0000) 
(  0.9804   0.0000   0.0000) 
(  0.9744   0.0000   0.0000) 
(  0.9676   0.0000   0.0000) 
(  0.9600   0.0000   0.0000) 
(  0.9516   0.0000   0.0000) 
(  0.9424   0.0000   0.0000) 
(  0.9324   0.0000   0.0000) 
(  0.9216   0.0000   0.0000) 
(  0.9100   0.0000   0.0000) 
(  0.8976   0.0000   0.0000) 
(  0.8844   0.0000   0.0000) 
(  0.8704   0.0000   0.0000) 
(  0.8556   0.0000   0.0000) 
(  0.8400   0.0000   0.0000) 
(  0.8236   0.0000   0.0000) 
(  0.8064   0.0000   0.0000) 
(  0.7884   0.0000   0.0000) 
(  0.7696   0.0000   0.0000) 
(  0.7500   0.0000   0.0000) 
(  0.7296   0.0000   0.0000) 
(  0.7084   0.0000   0.0000) 
(  0.6864   0.0000   0.0000) 
(  0.6636   0.0000   0.0000) 
(  0.6400   0.0000   0.0000) 
(  0.6156   0.0000   0.0000) 
(  0.5904   0.0000   0.0000) 
(  0.5644   0.0000   0.0000) 
(  0.5376   0.0000   0.0000) 
(  0.5100   0.0000   0.0000) 
(  0.4816   0.0000   0.0000) 
(  0.4524   0.0000   0.0000) 
(  0.4224   0.0000   0.0000) 
(  0.3916   0.0000   0.0000) 
(  0.3600   0.0000   0.0000) 
(  0.3276   0.0000   0.0000) 
(  0.2944   0.0000   0.0000) 
(  0.2604   0.0000   0.0000) 
(  0.2256   0.0000   0.0000) 
(  0.1900   0.0000   0.0000) 
(  0.1536   0.0000   0.0000) 
(  0.1164   0.0000   0.0000) 
(  0.0784   0.0000   0.0000) 
(  0.0396   0.0000   0.0000) 
(  0.0000   0.0000   0.0000) 

// maxz
(  1.0000   0.0000   0.0000) 
(  0.9996   0.0000   0.0000) 
(  0.9984   0.0000   0.0000) 
(  0.9964   0.0000   0.0000) 
(  0.9936   0.0000   0.0000) 
(  0.9900   0.0000   0.0000) 
(  0.9856   0.0000   0.0000) 
(  0.9804   0.0000   0.0000) 
(  0.9744   0.0000   0.0000) 
(  0.9676   0.0000   0.0000) 
(  0.9600   0.0000   0.0000) 
(  0.9516   0.0000   0.0000) 
(  0.9424   0.0000   0.0000) 
(  0.9324   0.0000   0.0000) 
(  0.9216   0.0000   0.0000) 
(  0.9100   0.0000   0.0000) 
(  0.8976   0.0000   0.0000) 
(  0.8844   0.0000   0.0000) 
(  0.8704   0.0000   0.0000) 
(  0.8556   0.0000   0.0000) 
(  0.8400   0.0000   0.0000) 
(  0.8236   0.0000   0.0000) 
(  0.8064   0.0000   0.0000) 
(  0.7884   0.0000   0.0000) 
(  0.7696   0.0000   0.0000) 
(  0.7500   0.0000   0.0000) 
(  0.7296   0.0000   0.0000) 
(  0.7084   0.0000   0.0000) 
(  0.6864   0.0000   0.0000) 
(  0.6636   0.0000   0.0000) 
(  0.6400   0.0000   0.0000) 
(  0.6156   0.0000   0.0000) 
(  0.5904   0.0000   0.0000) 
(  0.5644   0.0000   0.0000) 
(  0.5376   0.0000   0.0000) 
(  0.5100   0.0000   0.0000) 
(  0.4816   0.0000   0.0000) 
(  0.4524   0.0000   0.0000) 
(  0.4224   0.0000   0.0000) 
(  0.3916   0.0000   0.0000) 
(  0.3600   0.0000   0.0000) 
(  0.3276   0.0000   0.0000) 
(  0.2944   0.0000   0.0000) 
(  0.2604   0.0000   0.0000) 
(  0.2256   0.0000   0.0000) 
(  0.1900   0.0000   0.0000) 
(  0.1536   0.0000   0.0000) 
(  0.1164   0.0000   0.0000) 
(  0.0784   0.0000   0.0000) 
(  0.0396   0.0000   0.0000) 
(  0.0000   0.0000   0.0000) 
)

// ************************************************************************* //

My blockMesh runs fine but when I run my solver, I am getting this error:

Code:

Create time

Create mesh for time = 0


PIMPLE: Operating solver in PISO mode

Reading electricalProperties

Reading field p_rgh

Reading field U


DTRIS2 - Fatal error!


--> FOAM FATAL ERROR: 
Failed dtris2 with vertices:102

    From function static Foam::triSurface Foam::triSurfaceTools::delaunay2D(const Foam::List<Foam::Vector2D<double> >&)
    in file triSurface/triSurfaceTools/triSurfaceTools.C at line 2568.

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::error::abort() at ??:?
#2  Foam::triSurfaceTools::delaunay2D(Foam::List<Foam::Vector2D<double> > const&) at ??:?
#3  Foam::pointToPointPlanarInterpolation::calcWeights(Foam::Field<Foam::Vector<double> > const&, Foam::Field<Foam::Vector<double> > const&) at ??:?
#4  Foam::pointToPointPlanarInterpolation::pointToPointPlanarInterpolation(Foam::Field<Foam::Vector<double> > const&, Foam::Field<Foam::Vector<double> > const&, double, bool) at ??:?
#5  Foam::PatchFunction1Types::MappedFile<Foam::Vector<double> >::checkTable(double) const at ??:?
#6  Foam::PatchFunction1Types::MappedFile<Foam::Vector<double> >::value(double) const at ??:?
#7  Foam::timeVaryingMappedFixedValueFvPatchField<Foam::Vector<double> >::updateCoeffs() at ??:?
#8  Foam::timeVaryingMappedFixedValueFvPatchField<Foam::Vector<double> >::timeVaryingMappedFixedValueFvPatchField(Foam::fvPatch const&, Foam::DimensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::dictionary const&) at ??:?
#9  Foam::fvPatchField<Foam::Vector<double> >::adddictionaryConstructorToTable<Foam::timeVaryingMappedFixedValueFvPatchField<Foam::Vector<double> > >::New(Foam::fvPatch const&, Foam::DimensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::dictionary const&) at ??:?
#10  Foam::fvPatchField<Foam::Vector<double> >::New(Foam::fvPatch const&, Foam::DimensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::dictionary const&) at ??:?
#11  Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::Boundary::readField(Foam::DimensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::dictionary const&) at ??:?
#12  Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::readFields(Foam::dictionary const&) at ??:?
#13  Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::readFields() at ??:?
#14  Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::GeometricField(Foam::IOobject const&, Foam::fvMesh const&, bool) at ??:?
#15  ? at ??:?
#16  __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#17  ? at ??:?
Aborted

hk19 · September 8, 2020, 14:21

Hey mcfdma.

I am too facing the same problem, have you solve it ?

petros · September 8, 2020, 16:47

Welcome hk19,

I am quoting a possible solution to the problem:

Quote:

Originally Posted by mattijs

timeVaryingMapped does a triangulation of the supplied points. You'll get this error if you don't have enough points for a triangulation or they are not in a plane. It works out the plane from the first three points. Switch on its debug flag to see which points it uses.

P.S. Always remember that your problem has been most probably answered elsewhere.

Best,
Petros

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
[OpenFOAM.org] compile error in dynamicMesh and thermophysicalModels libraries	NickG	OpenFOAM Installation	3	December 30, 2019 01:21
[blockMesh] blockMesh with double grading.	spwater	OpenFOAM Meshing & Mesh Conversion	92	January 12, 2019 10:00
[swak4Foam] GroovyBC the dynamic cousin of funkySetFields that lives on the suburb of the mesh	gschaider	OpenFOAM Community Contributions	300	October 29, 2014 19:00
OpenFOAM without MPI	kokizzu	OpenFOAM Installation	4	May 26, 2014 10:17
DecomposePar links against liblamso0 with OpenMPI	jens_klostermann	OpenFOAM Bugs	11	June 28, 2007 18:51

September 8, 2020, 14:21		#2
hk19 New Member Hrishi Join Date: Sep 2020 Posts: 1 Rep Power: 0	Hey mcfdma. I am too facing the same problem, have you solve it ?