strange problem with data

seav · April 11, 2014, 05:01

Hello FOMERS.

I dont know how to name my problem properly so it goes like this.

The problem is really strange for me. I want to make simple-test-checking simulation with setting and naming internal faces. I made simple blockMeshDict with vertices, blocks and patches. I created starting 0/ folder, system/ folder etc. and run my calculation without any problem. I was supprised when I was checking data collected.. 1st step at boundary conditions seems to be good, but then 2nd and whole steps after are red.

GEOMETERY : Its a small box in bigger box. attach 1, 2

1st STEP : Seems be ok. attach 3,4
PROBLEM : Red planes after calculation. attach 5

blockMeshDict:

Code:

// ------------------------- //
FoamFile
{
	version	2.0;
	format	ascii;
	class	dictionary;
	object	blockMeshDict;
}

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


convertToMeters 0.001;

vertices
(
	(0 0 0) // 0
	(0 0 1) // 1
	(0 0 2) // 2
	(0 0 3) // 3
	(1 0 0) // 4
	(1 0 1) // 5
	(1 0 2) // 6
	(1 0 3) // 7
	(2 0 0) // 8
	(2 0 1) // 9
	(2 0 2) // 10
	(2 0 3) // 11
	(3 0 0) // 12
	(3 0 1) // 13
	(3 0 2) // 14
	(3 0 3) // 15
	(0 1 0) // 16
	(0 1 1) // 17
	(0 1 2) // 18
	(0 1 3) // 19
	(1 1 0) // 20
	(1 1 1) // 21
	(1 1 2) // 22
	(1 1 3) // 23
	(2 1 0) // 24
	(2 1 1) // 25
	(2 1 2) // 26
	(2 1 3) // 27
	(3 1 0) // 28
	(3 1 1) // 29
	(3 1 2) // 30
	(3 1 3) // 31
	(0 2 0) // 32
	(0 2 1) // 33
	(0 2 2) // 34
	(0 2 3) // 35
	(1 2 0) // 36
	(1 2 1) // 37
	(1 2 2) // 38
	(1 2 3) // 39
	(2 2 0) // 40
	(2 2 1) // 41
	(2 2 2) // 42
	(2 2 3) // 43
	(3 2 0) // 44
	(3 2 1) // 45
	(3 2 2) // 46
	(3 2 3) // 47
	(0 3 0) // 48
	(0 3 1) // 49
	(0 3 2) // 50
	(0 3 3) // 51
	(1 3 0) // 52
	(1 3 1) // 53
	(1 3 2) // 54
	(1 3 3) // 55
	(2 3 0) // 56
	(2 3 1) // 57
	(2 3 2) // 58
	(2 3 3) // 59
	(3 3 0) // 60
	(3 3 1) // 61
	(3 3 2) // 62
	(3 3 3) // 63
	
	(1 1 1) // 64
	(2 1 1) // 65
	(2 2 1) // 66
	(1 2 1) // 67
	(1 1 2) // 68
	(2 1 2) // 69
	(2 2 2) // 70
	(1 2 2) // 71

);
blocks
(
	hex (0 4 20 16 1 5 21 17) (3 3 3) simpleGrading (1 1 1)
	hex (4 8 24 20 5 9 25 21) (3 3 3) simpleGrading (1 1 1)
	hex (8 12 28 24 9 13 29 25) (3 3 3) simpleGrading (1 1 1)
	hex (1 5 21 17 2 6 22 18) (3 3 3) simpleGrading (1 1 1)
	hex (5 9 25 21 6 10 26 22) (3 3 3) simpleGrading (1 1 1)
	hex (9 13 29 25 10 14 30 26) (3 3 3) simpleGrading (1 1 1)
	hex (2 6 22 18 3 7 23 19) (3 3 3) simpleGrading (1 1 1)
	hex (6 10 26 22 7 11 27 23) (3 3 3) simpleGrading (1 1 1)
	hex (10 14 30 26 11 15 31 27) (3 3 3) simpleGrading (1 1 1)
	hex (16 20 36 32 17 21 37 33) (3 3 3) simpleGrading (1 1 1)
	hex (20 24 40 36 21 25 41 37) (3 3 3) simpleGrading (1 1 1)
	hex (24 28 44 40 25 29 45 41) (3 3 3) simpleGrading (1 1 1)
	hex (17 21 37 33 18 22 38 34) (3 3 3) simpleGrading (1 1 1)
	hex (21 25 41 37 22 26 42 38) (3 3 3) simpleGrading (1 1 1)
	hex (25 29 45 41 26 30 46 42) (3 3 3) simpleGrading (1 1 1)
	hex (18 22 38 34 19 23 39 35) (3 3 3) simpleGrading (1 1 1)
	hex (22 26 42 38 23 27 43 39) (3 3 3) simpleGrading (1 1 1)
	hex (26 30 46 42 27 31 47 43) (3 3 3) simpleGrading (1 1 1)
	hex (32 36 52 48 33 37 53 49) (3 3 3) simpleGrading (1 1 1)
	hex (36 40 56 52 37 41 57 53) (3 3 3) simpleGrading (1 1 1)
	hex (40 44 60 56 41 45 61 57) (3 3 3) simpleGrading (1 1 1)
	hex (33 37 53 49 34 38 54 50) (3 3 3) simpleGrading (1 1 1)
	hex (37 41 57 53 38 42 58 54) (3 3 3) simpleGrading (1 1 1)
	hex (41 45 61 57 42 46 62 58) (3 3 3) simpleGrading (1 1 1)
	hex (34 38 54 50 35 39 55 51) (3 3 3) simpleGrading (1 1 1)
	hex (38 42 58 54 39 43 59 55) (3 3 3) simpleGrading (1 1 1)
	hex (42 46 62 58 43 47 63 59) (3 3 3) simpleGrading (1 1 1)
	
	hex (64 65 66 67 68 69 70 71) last_block (3 3 3) simpleGrading (1 1 1)
);
patches
(
	wall plate
		(
		(64 65 66 67)
		(68 69 70 71)
		(64 68 71 67)
		(69 65 66 70)
		(68 69 65 64)
		(71 70 66 67)
		)

	wall nate
		(
		(0 4 20 16)
		(4 8 24 20)
		(8 12 28 24)
		(16 20 36 32)
		(20 24 40 36)
		(24 28 44 40)
		(32 36 52 48)
		(36 40 56 52)
		(40 44 60 56)	// end front
		(3 7 23 19)	
		(7 11 27 23)
		(11 15 31 27)
		(19 23 39 35)
		(23 27 43 39)
		(27 31 47 43)
		(35 39 55 51)
		(39 43 59 55)
		(43 47 63 59)	// end back
		(0 1 17 16)
		(1 2 18 17)
		(2 3 19 18)
		(16 17 33 32)
		(18 19 35 34)
		(32 33 49 48)
		(33 34 50 49)
		(34 35 51 50)	// begin front - inlet
		(15 14 30 31)
		(14 13 29 30)
		(13 12 28 29)
		(31 30 46 47)	// begin back - outlet
		(29 28 44 45)
		(47 46 62 63)
		(46 45 61 62)
		(45 44 60 61)
		(51 55 54 50)
		(55 59 58 54)
		(59 63 62 58)
		(50 54 53 49)
		(54 58 57 53)
		(58 62 61 57)
		(49 53 52 48)
		(53 57 56 52)
		(57 61 60 56)	//end up plate
		(3 7 6 2)
		(7 11 10 6)
		(11 15 14 10)
		(2 6 5 1)
		(6 10 9 5)
		(10 14 13 9)
		(1 5 4 0)
		(5 9 8 4)
		(9 13 12 8)
		)
	
	patch inlet
		(
		(17 18 34 33)
		)
	
	patch outlet
		(
		(30 29 45 46)
		)
	
	/*defaultFaces
	{
		type empty;
		faces();
	}*/
);
mergePatchPairs
(
);

controlDict:

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.2.2                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

application     icoFoam;

startFrom       startTime;

startTime       0;

stopAt          endTime;

endTime         0.5;

deltaT          0.005;

writeControl    timeStep;

writeInterval   20;

purgeWrite      0;

writeFormat     ascii;

writePrecision  6;

writeCompression off;

timeFormat      general;

timePrecision   6;

runTimeModifiable true;

fvSchemes:

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.2.2                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
    default         Euler;
}

gradSchemes
{
    default         Gauss linear;
    grad(p)         Gauss linear;
}

divSchemes
{
    default         none;
    div(phi,U)      Gauss linear;
}

laplacianSchemes
{
    default         none;
    laplacian(nu,U) Gauss linear orthogonal;
    laplacian((1|A(U)),p) Gauss linear orthogonal;
}

interpolationSchemes
{
    default         linear;
    interpolate(HbyA) linear;
}

snGradSchemes
{
    default         orthogonal;
}

fluxRequired
{
    default         no;
    p               ;
}


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

fvSolution:

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.2.2                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
    p
    {
        solver          PCG;
        preconditioner  DIC;
        tolerance       1e-06;
        relTol          0;
    }

    U
    {
        solver          PBiCG;
        preconditioner  DILU;
        tolerance       1e-05;
        relTol          0;
    }
}

PISO
{
    nCorrectors     2;
    nNonOrthogonalCorrectors 0;
    pRefCell        0;
    pRefValue       0;
}


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

0/p:

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.2.2                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 2 -2 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    	inlet 
    	{
        	type	fixedValue;
		value	uniform 3;
    	}
   	
	outlet
   	{
        	type 	zeroGradient;
   	}
    
	nate
	{
		type 	zeroGradient;
	}

	plate
	{
		type 	zeroGradient;
	}
}

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

0/U:

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.2.2                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volVectorField;
    object      U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (0 0 0);

boundaryField
{
	inlet
	{
		type            fixedValue;
		value		uniform (0.368 0 0);
	}

	outlet
	{
		type 		zeroGradient;
	}
	
	nate
	{
		type 		zeroGradient;
	}

	plate
	{
		type		zeroGradient;
	}
}

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

I am using OF 2.2.2 and ParaView-3.12

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Strange Flamelet Generation Problem	tar	FLUENT	11	March 22, 2014 10:52
Problem in Interpolating data command	shubham208011	FLUENT	6	July 15, 2013 08:18
[Commercial meshers] fluentMeshToFoam multidomain mesh conversion problem	Attesz	OpenFOAM Meshing & Mesh Conversion	12	May 2, 2013 10:52
Problem in running ICEM grid in Openfoam	Tarak	OpenFOAM	6	September 9, 2011 17:51
Anyone met this strange problem?	Lee	Siemens	1	June 7, 2005 03:15