 |
 |
 |
|
|
|
[Sponsors] | ||||
rSF: p divergence in combustor (wt negative value) |
 |
|
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
April 6, 2018, 04:03
|
rSF: p divergence in combustor (wt negative value)
|
#1 |
|
New Member
Francesco Liuzzi
Join Date: Jul 2017
Posts: 15
Rep Power: 8  |
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);
}
}
// ************************************************************************* //
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;*/
}
}
// ************************************************************************* //
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;
}
}
// ************************************************************************* //
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;
}
}
// ************************************************************************* //
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;
}
}
// ************************************************************************* //
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;
}
}
// ************************************************************************* //
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;
// ************************************************************************* //
Of course, let me know if you need more file or info and i will provide them. Thanks in advance for the attention
|
|
|
|
|
|
April 6, 2018, 05:11
|
|
#2 |
|
Member
Hosein
Join Date: Nov 2011
Location: Germany
Posts: 93
Rep Power: 14 |
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. |
|
|
|
|
|
|
April 6, 2018, 05:51
|
|
#3 |
|
New Member
Francesco Liuzzi
Join Date: Jul 2017
Posts: 15
Rep Power: 8 |
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;
}*/
}
// ************************************************************************* //
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;
}
// ************************************************************************* //
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
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
|
|
|
|
|
|
|
April 7, 2018, 05:54
|
|
#4 |
|
Member
Hosein
Join Date: Nov 2011
Location: Germany
Posts: 93
Rep Power: 14 |
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... |
|
|
|
|
|
|
April 10, 2018, 06:59
|
|
#5 |
|
New Member
Francesco Liuzzi
Join Date: Jul 2017
Posts: 15
Rep Power: 8 |
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
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
|
|
|
|
|
|
|
|
|
Similar Threads
|
||||
| Thread | Thread Starter | Forum | Replies | Last Post |
| [Other] mesh airfoil NACA0012 | anand_30 | OpenFOAM Meshing & Mesh Conversion | 13 | March 7, 2022 17:22 |
| 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 |
| [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 |
Linear Mode
