 |
 |
 |
|
|
|
[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
|
|
|
|
|
|
|
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 |
Threaded Mode
