CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Pre-Processing

rSF: p divergence in combustor (wt negative value)

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 6, 2018, 04:03
Default rSF: p divergence in combustor (wt negative value)
  #1
New Member
 
Francesco Liuzzi
Join Date: Jul 2017
Posts: 15
Rep Power: 5
zonda is on a distinguished road
Good morning to all of you
Right now i'm working on a simulation involving the valuation of an air flux inside a combustion chamber (no combustion, only air flux) with rhoSimpleFoam.
You can see the geometry i'm working on in the picture below:
chamber.jpg

Where is the problem? It doesn't matter how i set the inputs (U, p, T and the relative k, epsilon and nut evaluated with the formulas found on this site), there will always be a divergence in pressure inside the swirlers. Openfoam returns a "floating point exception" error and when i visualize the last 2-3 timesteps before the crash on paraview it appears always like this:
press-divergence.jpg

If it can help, i don't have to work with a given air mass flux, i just need a stable situation after a large number of time step, but however, the higher the mass flux, the better.

Here i'm attaching the U, p, T, k, epsilon, nut and the controlDict files:

U
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  3.0.x                                 |
|   \\  /    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_swirl-raggera-1-intersecare2_3_
    {
        type            flowRateInletVelocity;
        massFlowRate    constant 0.0063021;
	rho		rho;
        rhoInlet        0.7;    // Guess for rho
    }
    inlet_swirl-raggera-1-intersecare2_1_
    {
        type            flowRateInletVelocity;
        massFlowRate    constant 0.0063021;
	rho		rho;
        rhoInlet        0.7;    // Guess for rho
    }
    inlet_swirl-raggera-1-intersecare2_8_
    {
        type            flowRateInletVelocity;
        massFlowRate    constant 0.0063021;
	rho		rho;
        rhoInlet        0.7;    // Guess for rho
    }
    inlet_swirl-raggera-1-intersecare2_6_
    {
        type            flowRateInletVelocity;
        massFlowRate    constant 0.0063021;
	rho		rho;
        rhoInlet        0.7;    // Guess for rho
    }
    inlet_swirl-raggera-1-intersecare2_4_
    {
        type            flowRateInletVelocity;
        massFlowRate    constant 0.0063021;
	rho		rho;
        rhoInlet        0.7;    // Guess for rho
    }
    inlet_swirl-raggera-1-intersecare2_2_
    {
        type            flowRateInletVelocity;
        massFlowRate    constant 0.0063021;
	rho		rho;
        rhoInlet        0.7;    // Guess for rho
    }
    inlet_swirl-raggera-1-intersecare2_7_
    {
        type            flowRateInletVelocity;
        massFlowRate    constant 0.0063021;
	rho		rho;
        rhoInlet        0.7;    // Guess for rho
    }
    inlet_swirl-raggera-1-intersecare2_5_
    {
        type            flowRateInletVelocity;
        massFlowRate    constant 0.0063021;
	rho		rho;
        rhoInlet        0.7;    // Guess for rho
    }
    inlet1
    {
        type            flowRateInletVelocity;
        massFlowRate    constant 0.0011458;
	rho		rho;
        rhoInlet        0.7;    // Guess for rho
    }
    inlet2
    {
        type            flowRateInletVelocity;
        massFlowRate    constant 0.0011458;
	rho		rho;
        rhoInlet        0.7;    // Guess for rho
    }
    inlet3
    {
        type            flowRateInletVelocity;
        massFlowRate    constant 0.0011458; 
	rho		rho;
        rhoInlet        0.7;    // Guess for rho
    }
    inlet4
    {
        type            flowRateInletVelocity;
        massFlowRate    constant 0.0011458; 
	rho		rho;
        rhoInlet        0.7;    // Guess for rho
    }
    walls_union
    {
        type            fixedValue;
	value		uniform (0 0 0);
    }
    Cavcilnolayer
    {
        type            fixedValue;
	value		uniform (0 0 0);
    }
    outlet_union
    {
        //type            zeroGradient;
	type		inletOutlet;
	inletValue	uniform (0 0 0);
	value 		uniform (0 0 0);
    }
}

// ************************************************************************* //
p

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

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

internalField   uniform 1.0e5;

boundaryField
{
    inlet_swirl-raggera-1-intersecare2_3_
    {
        type            zeroGradient;
    }
    inlet_swirl-raggera-1-intersecare2_1_
    {
        type            zeroGradient;
    }
    inlet_swirl-raggera-1-intersecare2_8_
    {
        type            zeroGradient;
    }
    inlet_swirl-raggera-1-intersecare2_6_
    {
        type            zeroGradient;
    }
    inlet_swirl-raggera-1-intersecare2_4_
    {
        type            zeroGradient;
    }
    inlet_swirl-raggera-1-intersecare2_2_
    {
        type            zeroGradient;
    }
    inlet_swirl-raggera-1-intersecare2_7_
    {
        type            zeroGradient;
    }
    inlet_swirl-raggera-1-intersecare2_5_
    {
        type            zeroGradient;
    }
    inlet1
    {
        type            zeroGradient;
    }
    inlet2
    {
        type            zeroGradient;
    }
    inlet3
    {
        type            zeroGradient;
    }
    inlet4
    {
        type            zeroGradient;
    }
    walls_union
    {
        type            zeroGradient;
    }
    Cavcilnolayer
    {
        type            zeroGradient;
    }
    outlet_union
    {
	type		fixedValue;
	value		uniform 1e5;
        /*type            totalPressure;
	U		U;
	phi		phi;
	rho		none;
	psi		thermo:psi;
        gamma		1.4;
	p0		uniform 1e5;*/
    }
}

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

dimensions      [0 0 0 1 0 0 0];

internalField   uniform 500;

boundaryField
{
    inlet_swirl-raggera-1-intersecare2_3_
    {
        type            fixedValue;
        value           uniform 500;
    }
    inlet_swirl-raggera-1-intersecare2_1_
    {
        type            fixedValue;
        value           uniform 500;
    }
    inlet_swirl-raggera-1-intersecare2_8_
    {
        type            fixedValue;
        value           uniform 500;
    }
    inlet_swirl-raggera-1-intersecare2_6_
    {
        type            fixedValue;
        value           uniform 500;
    }
    inlet_swirl-raggera-1-intersecare2_4_
    {
        type            fixedValue;
        value           uniform 500;
    }
    inlet_swirl-raggera-1-intersecare2_2_
    {
        type            fixedValue;
        value           uniform 500;
    }
    inlet_swirl-raggera-1-intersecare2_7_
    {
        type            fixedValue;
        value           uniform 500;
    }
    inlet_swirl-raggera-1-intersecare2_5_
    {
        type            fixedValue;
        value           uniform 500;
    }
    inlet1
    {
        type            fixedValue;
        value           uniform 500;
    }
    inlet2
    {
        type            fixedValue;
        value           uniform 500;
    }
    inlet3
    {
        type            fixedValue;
        value           uniform 500;
    }
    inlet4
    {
        type            fixedValue;
        value           uniform 500;
    }
    walls_union
    {
        type            zeroGradient;
    }
    Cavcilnolayer
    {
        type            zeroGradient;
    }
    outlet_union
    {
		type            zeroGradient;
    }
}

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

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

internalField   uniform 10.596;

boundaryField
{
    inlet_swirl-raggera-1-intersecare2_3_
    {
        type            fixedValue;
        value		uniform 10.596; //2.677;
    }
    inlet_swirl-raggera-1-intersecare2_1_
    {
        type            fixedValue;
        value		uniform 10.596; 
    }
    inlet_swirl-raggera-1-intersecare2_8_
    {
        /*type            turbulentIntensityKineticEnergyInlet;
        intensity       0.05;
        value           uniform 1;*/
        type            fixedValue;
        value		uniform 10.596;
    }
    inlet_swirl-raggera-1-intersecare2_6_
    {
        type            fixedValue;
        value		uniform 10.596;
    }
    inlet_swirl-raggera-1-intersecare2_4_
    {
        type            fixedValue;
        value		uniform 10.596;
    }
    inlet_swirl-raggera-1-intersecare2_2_
    {
        type            fixedValue;
        value		uniform 10.596;
    }
    inlet_swirl-raggera-1-intersecare2_7_
    {
        type            fixedValue;
        value		uniform 10.596;
    }
    inlet_swirl-raggera-1-intersecare2_5_
    {
        type            fixedValue;
        value		uniform 10.596;
    }
    inlet1
    {
        type            fixedValue;
        value		uniform 25.419; //6.401;
    }
    inlet2
    {
        type            fixedValue;
        value		uniform 25.419;
    }
    inlet3
    {
        type            fixedValue;
        value		uniform 25.419;
    }
    inlet4
    {
        type            fixedValue;
        value		uniform 25.419;
    }
    walls_union
    {
        type            kqRWallFunction;
        value           uniform 0;
    }
    Cavcilnolayer
    {
        type            kqRWallFunction;
        value           uniform 0;
    }
    outlet_union
    {
        type            zeroGradient;
    }
}

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

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

internalField   uniform 1227.4;

boundaryField
{
    inlet_swirl-raggera-1-intersecare2_3_
    {
        /*type            turbulentMixingLengthDissipationRateInlet;
        mixingLength    0.005;
        value           uniform 200;*/
	type		fixedValue;
	value		uniform 1227.4; //845.89;
    }
    inlet_swirl-raggera-1-intersecare2_1_
    {
	type		fixedValue;
	value		uniform 1227.4;
    }
    inlet_swirl-raggera-1-intersecare2_8_
    {
	type		fixedValue;
	value		uniform 1227.4;
    }
    inlet_swirl-raggera-1-intersecare2_6_
    {
	type		fixedValue;
	value		uniform 1227.4;
    }
    inlet_swirl-raggera-1-intersecare2_4_
    {
	type		fixedValue;
	value		uniform 1227.4;
    }
    inlet_swirl-raggera-1-intersecare2_2_
    {
	type		fixedValue;
	value		uniform 1227.4;
    }
    inlet_swirl-raggera-1-intersecare2_7_
    {
	type		fixedValue;
	value		uniform 1227.4;
    }
    inlet_swirl-raggera-1-intersecare2_5_
    {
	type		fixedValue;
	value		uniform 1227.4;
    }
    inlet1
    {
	type		fixedValue;
	value		uniform 11067; //7592.45;
    }
    inlet2
    {
	type		fixedValue;
	value		uniform 11067;
    }
    inlet3
    {
	type		fixedValue;
	value		uniform 11067;
    }
    inlet4
    {
	type		fixedValue;
	value		uniform 11067;
    }
    walls_union
    {
        type            epsilonWallFunction;
        value           $internalField;
    }
    Cavcilnolayer
    {
        type            epsilonWallFunction;
        value           $internalField;
    }
    outlet_union
    {
        type            zeroGradient;
    }
}

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

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

internalField   uniform 0.000823;

boundaryField
{
    inlet_swirl-raggera-1-intersecare2_3_
    {
        type            fixedValue;
	value		uniform 0.000823;
    }
    inlet_swirl-raggera-1-intersecare2_1_
    {
        /*type            calculated;
        value           uniform 0;*/
	type            fixedValue;
	value		uniform 0.000823;
    }
    inlet_swirl-raggera-1-intersecare2_8_
    {
	type            fixedValue;
	value		uniform 0.000823;
    }
    inlet_swirl-raggera-1-intersecare2_6_
    {
	type            fixedValue;
	value		uniform 0.000823;
    }
    inlet_swirl-raggera-1-intersecare2_4_
    {
	type            fixedValue;
	value		uniform 0.000823;
    }
    inlet_swirl-raggera-1-intersecare2_2_
    {
	type            fixedValue;
	value		uniform 0.000823;
    }
    inlet_swirl-raggera-1-intersecare2_7_
    {
	type            fixedValue;
	value		uniform 0.000823;
    }
    inlet_swirl-raggera-1-intersecare2_5_
    {
	type            fixedValue;
	value		uniform 0.000823;
    }
    inlet1
    {
	type            fixedValue;
	value		uniform 0.0005254;
    }
    inlet2
    {
	type            fixedValue;
	value		uniform 0.0005254;
    }
    inlet3
    {
	type            fixedValue;
	value		uniform 0.0005254;
    }
    inlet4
    {
	type            fixedValue;
	value		uniform 0.0005254;
    }
    walls_union
    {
        type            nutkWallFunction;
        value           $internalField;
    }
    Cavcilnolayer
    {
        type            nutkWallFunction;
        value           $internalField;
    }
    outlet_union
    {
        type            zeroGradient;
    }
}

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

application     rhoSimpleFoam;

startFrom       latestTime;

startTime       0;

stopAt          endTime;

endTime         5000; //10000

deltaT          1;

writeControl    timeStep;

writeInterval   10; //100

purgeWrite      8; //2

writeFormat     ascii;

writePrecision  6;

writeCompression off;

timeFormat      general;

timePrecision   6;

graphFormat     raw;

runTimeModifiable true;


// ************************************************************************* //
Any suggestion on how to solve this problem?
Of course, let me know if you need more file or info and i will provide them.
Thanks in advance for the attention
zonda is offline   Reply With Quote

Old   April 6, 2018, 05:11
Default
  #2
Member
 
Hosein
Join Date: Nov 2011
Location: Germany
Posts: 81
Rep Power: 11
einstein_zee is on a distinguished road
Hi there,

It would be good if you could also provide the content of fvSolution and fvSchemes. I suppose based on the nature of your problem you need to perform unsteady simulations. But what I cannot see is a proper deltaT in your controlDict file. Also there is no control over the Courant number. Is your mesh big enough to compensate for a reasonable Courant number with your current time step?
One more important thing is that since your geometry is a bit complex, proper meshing also plays an important role here. Maybe it's also good that you provide us the result of either checkMesh or checkMesh -allTopology -allGeometry.
I'm afraid more support will be provided once your case information is given in detail.

Hope this helps.
einstein_zee is offline   Reply With Quote

Old   April 6, 2018, 05:51
Default
  #3
New Member
 
Francesco Liuzzi
Join Date: Jul 2017
Posts: 15
Rep Power: 5
zonda is on a distinguished road
Ok, first of all, thank you einstein_zee for your answer
Right now i'm running a steady simulation because, being an absolute beginner in openFoam (this is my first project), i received the instruction to run a steady one to start and then, with the obtained results, i will initialize the unsteady one (with rhoPimpeFoam).
Howevere, here are the files you asked for:

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

solvers
{
    p
    {
        solver          GAMG;
        tolerance       1e-08;
        relTol          0.1;
        smoother        GaussSeidel;

	directSolveCoarsest	false;
	agglomerator		faceAreaPair;
	mergeLevels		1;
	nPreSweeps		0;
	nPostSweeps		2;
	nFinestSweeps		2;
        nCellsInCoarsestLevel 8;
    }

    "(U|e|k|epsilon|nut)"
    {
        solver          GAMG;
        tolerance       1e-08;
        relTol          0.1;
        smoother        GaussSeidel;
        nCellsInCoarsestLevel 20;
    }

    Phi
    {
        $p;
        relTol          0;
    }

}

potentialFlow
{
	nNonOrthogonalCorrectors 10;
}

SIMPLE
{
    nNonOrthogonalCorrectors 3;
    rhoMin          0.2;
    rhoMax          1.4;
    //transonic	    yes;
    //consistent	    yes;

    residualControl
    {
        p               1e-3;
        U               1e-4;
        e               1e-3;

        // possibly check turbulence fields
        "(k|epsilon|omega|nut)" 1e-3;
    }
}

relaxationFactors
{

    fields
    {
        p               0.05;
        rho             0.05;
    }

    equations
    {
        p               0.1;
        U               0.3;
        e               0.3;
        k               0.3;
        epsilon         0.3;
	nut		0.3;
    }

    /*fields
    {
        p               0.05; //0.3;
        rho             0.05;
    }
    equations
    {
        U               	0.3; //0.7;
        "(k|epsilon|omega|nut)"	0.3; //0.7;
        e               	0.1; //0.5;
    }*/
}


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

ddtSchemes
{
    default         steadyState;
}

gradSchemes
{
    default         Gauss linear;
}

divSchemes
{
    default         none;

    div(phi,U)      				bounded Gauss GammaV 1;
    div(((rho*nuEff)*dev2(T(grad(U))))) 	Gauss linear;
    div(phi,e)      				bounded Gauss upwind;
    div(phi,omega) 				bounded Gauss upwind;
    div(phi,epsilon) 				bounded Gauss upwind;
    div(phi,k)      				bounded Gauss upwind;
    div(phi,Ekp)    				bounded Gauss upwind;
    div(phid,p)         			bounded Gauss upwind;
    div((phi|interpolate(rho)),p)  		bounded Gauss upwind;
}

laplacianSchemes
{
    default         Gauss linear corrected;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         corrected;
}

wallDist
{
    method meshWave;
}
// ************************************************************************* //
checkMesh
Code:
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  3.0.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 3.0.1-d8a290b55d28
Exec   : checkMesh
Date   : Apr 06 2018
Time   : 12:39:59
Host   : "an02"
PID    : 21005
Case   : /utenti/simulation/OpenFOAM/simulation-3.0.1/run/cases/AeroSpaceProp/GreenEngine/francesco/ge_rhoSimpleFoam_Q25
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create polyMesh for time = 0

Time = 0

Mesh stats
    points:           1404913
    faces:            4141591
    internal faces:   4066063
    cells:            1368486
    faces per cell:   5.99762
    boundary patches: 15
    point zones:      1
    face zones:       25
    cell zones:       15

Overall number of cells of each type:
    hexahedra:     1360217
    prisms:        514
    wedges:        0
    pyramids:      2882
    tet wedges:    0
    tetrahedra:    1581
    polyhedra:     3292
    Breakdown of polyhedra by number of faces:
        faces   number of cells
            6   102
            7   3084
            8   106

Checking topology...
    Boundary definition OK.
    Cell to face addressing OK.
    Point usage OK.
    Upper triangular ordering OK.
    Face vertices OK.
    Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces...
                   Patch    Faces   Points                  Surface topology
    inlet_swirl-raggera-1-intersecare2_3_      150      176  ok (non-closed singly connected)
    inlet_swirl-raggera-1-intersecare2_1_      150      176  ok (non-closed singly connected)
    inlet_swirl-raggera-1-intersecare2_8_      150      176  ok (non-closed singly connected)
    inlet_swirl-raggera-1-intersecare2_6_      150      176  ok (non-closed singly connected)
    inlet_swirl-raggera-1-intersecare2_4_      150      176  ok (non-closed singly connected)
    inlet_swirl-raggera-1-intersecare2_2_      150      176  ok (non-closed singly connected)
    inlet_swirl-raggera-1-intersecare2_7_      150      176  ok (non-closed singly connected)
    inlet_swirl-raggera-1-intersecare2_5_      150      176  ok (non-closed singly connected)
                  inlet1       47       58  ok (non-closed singly connected)
                  inlet2       47       58  ok (non-closed singly connected)
                  inlet3       47       58  ok (non-closed singly connected)
                  inlet4       47       58  ok (non-closed singly connected)
             walls_union    60332    60917  ok (non-closed singly connected)
            outlet_union     6158     6210  ok (non-closed singly connected)
           Cavcilnolayer     7650     7752  ok (non-closed singly connected)

Checking geometry...
    Overall domain bounding box (-0.0704992 -0.0704916 -5.91789e-17) (0.0704979 0.0704914 0.5175)
    Mesh has 3 geometric (non-empty/wedge) directions (1 1 1)
    Mesh has 3 solution (non-empty) directions (1 1 1)
    Boundary openness (3.94623e-16 5.57443e-16 -2.33326e-16) OK.
    Max cell openness = 3.15242e-15 OK.
    Max aspect ratio = 130.559 OK.
    Minimum face area = 9.18654e-09. Maximum face area = 1.02318e-05.  Face area magnitudes OK.
    Min volume = 1.16807e-12. Max volume = 1.57123e-08.  Total volume = 0.00480227.  Cell volumes OK.
    Mesh non-orthogonality Max: 85.5305 average: 10.6305
   *Number of severely non-orthogonal (> 70 degrees) faces: 13818.
    Non-orthogonality check OK.
  <<Writing 13818 non-orthogonal faces to set nonOrthoFaces
    Face pyramids OK.
    Max skewness = 3.47496 OK.
    Coupled point location match (average 0) OK.

Mesh OK.

End
checkMesh -allTopology -allGeometry
Code:
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  3.0.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 3.0.1-d8a290b55d28
Exec   : checkMesh -allGeometry -allTopology
Date   : Apr 06 2018
Time   : 12:36:52
Host   : "an02"
PID    : 20993
Case   : /utenti/simulation/OpenFOAM/simulation-3.0.1/run/cases/AeroSpaceProp/GreenEngine/francesco/ge_rhoSimpleFoam_Q25
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create polyMesh for time = 0

Enabling all (cell, face, edge, point) topology checks.

Enabling all geometry checks.

Time = 0

Mesh stats
    points:           1404913
    faces:            4141591
    internal faces:   4066063
    cells:            1368486
    faces per cell:   5.99762
    boundary patches: 15
    point zones:      1
    face zones:       25
    cell zones:       15

Overall number of cells of each type:
    hexahedra:     1360217
    prisms:        514
    wedges:        0
    pyramids:      2882
    tet wedges:    0
    tetrahedra:    1581
    polyhedra:     3292
    Breakdown of polyhedra by number of faces:
        faces   number of cells
            6   102
            7   3084
            8   106

Checking topology...
    Boundary definition OK.
    Cell to face addressing OK.
    Point usage OK.
    Upper triangular ordering OK.
    Face vertices OK.
    Topological cell zip-up check OK.
    Face-face connectivity OK.
    Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces...
                   Patch    Faces   Points                  Surface topology Bounding box
    inlet_swirl-raggera-1-intersecare2_3_      150      176  ok (non-closed singly connected) (0.0114549 -0.0503399 0.1202) (0.0288562 -0.0404843 0.1317)
    inlet_swirl-raggera-1-intersecare2_1_      150      176  ok (non-closed singly connected) (-0.0503399 -0.0288562 0.1202) (-0.0404843 -0.0114549 0.1317)
    inlet_swirl-raggera-1-intersecare2_8_      150      176  ok (non-closed singly connected) (-0.0524698 0.011661 0.1202) (-0.0397846 0.0235849 0.1317)
    inlet_swirl-raggera-1-intersecare2_6_      150      176  ok (non-closed singly connected) (0.011661 0.0397846 0.1202) (0.0235849 0.0524698 0.1317)
    inlet_swirl-raggera-1-intersecare2_4_      150      176  ok (non-closed singly connected) (0.0397846 -0.0235849 0.1202) (0.0524698 -0.011661 0.1317)
    inlet_swirl-raggera-1-intersecare2_2_      150      176  ok (non-closed singly connected) (-0.0235849 -0.0524698 0.1202) (-0.011661 -0.0397846 0.1317)
    inlet_swirl-raggera-1-intersecare2_7_      150      176  ok (non-closed singly connected) (-0.0288562 0.0404843 0.1202) (-0.0114549 0.0503399 0.1317)
    inlet_swirl-raggera-1-intersecare2_5_      150      176  ok (non-closed singly connected) (0.0404843 0.0114549 0.1202) (0.0503399 0.0288562 0.1317)
                  inlet1       47       58  ok (non-closed singly connected) (-0.00247265 0.0147902 0.00439508) (0.00249994 0.015 0.00938048)
                  inlet2       47       58  ok (non-closed singly connected) (0.0147902 -0.00249994 0.00439508) (0.015 0.00247265 0.00938048)
                  inlet3       47       58  ok (non-closed singly connected) (-0.00249994 -0.015 0.00439508) (0.00247265 -0.0147902 0.00938048)
                  inlet4       47       58  ok (non-closed singly connected) (-0.015 -0.00247265 0.00439508) (-0.0147902 0.00249994 0.00938048)
             walls_union    60332    60917  ok (non-closed singly connected) (-0.0704992 -0.0704916 -5.91789e-17) (0.0704979 0.0704914 0.5175)
            outlet_union     6158     6210  ok (non-closed singly connected) (-0.070492 -0.0704916 0.5175) (0.0704918 0.0704914 0.5175)
           Cavcilnolayer     7650     7752  ok (non-closed singly connected) (-0.0704966 -0.0704913 0.33) (0.0704952 0.0704879 0.39)

Checking geometry...
    Overall domain bounding box (-0.0704992 -0.0704916 -5.91789e-17) (0.0704979 0.0704914 0.5175)
    Mesh has 3 geometric (non-empty/wedge) directions (1 1 1)
    Mesh has 3 solution (non-empty) directions (1 1 1)
    Boundary openness (3.94623e-16 5.57443e-16 -2.33326e-16) OK.
    Max cell openness = 3.15242e-15 OK.
    Max aspect ratio = 130.559 OK.
    Minimum face area = 9.18654e-09. Maximum face area = 1.02318e-05.  Face area magnitudes OK.
    Min volume = 1.16807e-12. Max volume = 1.57123e-08.  Total volume = 0.00480227.  Cell volumes OK.
    Mesh non-orthogonality Max: 85.5305 average: 10.6305
   *Number of severely non-orthogonal (> 70 degrees) faces: 13818.
    Non-orthogonality check OK.
  <<Writing 13818 non-orthogonal faces to set nonOrthoFaces
    Face pyramids OK.
    Max skewness = 3.47496 OK.
    Coupled point location match (average 0) OK.
    Face tets OK.
    Min/max edge length = 3.1e-05 0.00452139 OK.
    All angles in faces OK.
    Face flatness (1 = flat, 0 = butterfly) : min = 0.960373  average = 0.999985
    All face flatness OK.
    Cell determinant (wellposedness) : minimum: 8.6638e-06 average: 5.02002
 ***Cells with small determinant (< 0.001) found, number of cells: 20768
  <<Writing 20768 under-determined cells to set underdeterminedCells
 ***Concave cells (using face planes) found, number of cells: 3144
  <<Writing 3144 concave cells to set concaveCells
    Face interpolation weight : minimum: 0.00225489 average: 0.493188
 ***Faces with small interpolation weight (< 0.05) found, number of faces: 3569
  <<Writing 3569 faces with low interpolation weights to set lowWeightFaces
    Face volume ratio : minimum: 0.00549366 average: 0.968486
 ***Faces with small volume ratio (< 0.01) found, number of faces: 1
  <<Writing 1 faces with low volume ratio cells to set lowVolRatioFaces

Failed 4 mesh checks.

End
zonda is offline   Reply With Quote

Old   April 7, 2018, 05:54
Default
  #4
Member
 
Hosein
Join Date: Nov 2011
Location: Germany
Posts: 81
Rep Power: 11
einstein_zee is on a distinguished road
okay, one thing I noticed from your pics is that the geometries are different isn't it? Since in the picture where you made a cut there is a smaller diameter channel in the middle.

But anyways, one thing to mention is that you may change the BC for "outlet_union" in U and p and use pressureInletOutletVelocity and totalPressure respectively. And try to remove any complexity! you don't have to model turbulence in your first simulation. So instead use laminar in your turbulenceProperties file. Also you can initialize your internal temperature field with 500 and for all boundaries define zeroGradient.
You may try to test PCG for pressure equation and see its performance also. Also play with your tolerances! you may reduce it from e-08 to lower orders of magnitude (say e-04).

These are general guidelines at least for now I can suggest. But you can also compress everything and put your case here so that somebody may take a look and solve the issue for you . Finally to make your query more complete, you may also provide a figure with your residuals.

hope this helps...
einstein_zee is offline   Reply With Quote

Old   April 10, 2018, 06:59
Default
  #5
New Member
 
Francesco Liuzzi
Join Date: Jul 2017
Posts: 15
Rep Power: 5
zonda is on a distinguished road
Ok, according to your suggestions i have modified my case, but now all i get is this error:

Code:
--> FOAM FATAL IO ERROR:
Attempt to return dictionary entry as a primitive

file: /utenti/simulation/OpenFOAM/simulation-3.0.1/run/cases/AeroSpaceProp/GreenEngine/francesco/ge_rhoSimpleFoam_Q25/constant/thermophysicalProperties::thermoType from line 20 to line 26.

    From function ITstream& primitiveEntry::stream() const
    in file db/dictionary/dictionaryEntry/dictionaryEntry.C at line 82.

FOAM aborting
For what i have understood, it's related to an error in the sintax of some file (in this case probably the thermoPhisicalProperties that i didn't even touched) but i'm not completely sure.

However here is the case :
https://app.box.com/s/0iyi5kb0ckeleq1cu2hpylv0vn4fsoeg
Hope is the correct way to post a case, otherwise let me know
Thank you for your help
zonda is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

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
fluent divergence for no reason sufjanst FLUENT 2 March 23, 2016 16:08
[blockMesh] error message with modeling a cube with a hold at the center hsingtzu OpenFOAM Meshing & Mesh Conversion 2 March 14, 2012 09:56
[Other] mesh airfoil NACA0012 anand_30 OpenFOAM Meshing & Mesh Conversion 12 December 12, 2011 04:16
[blockMesh] BlockMesh FOAM warning gaottino OpenFOAM Meshing & Mesh Conversion 7 July 19, 2010 14:11
[blockMesh] Axisymmetrical mesh Rasmus Gjesing (Gjesing) OpenFOAM Meshing & Mesh Conversion 10 April 2, 2007 14:00


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