[AleQ95]

Hello everyone,

I'm working on my final bachelor degree thesis and I could use some help. I'm running a simulation using rhoCentralFoam for a mesh I did in pointwise. Here are my BCs for p,T,U

p
boundaryField
{
Inlet
{
type totalPressure;
gamma 1.4;
p0 uniform 8.25;

}



Boxin
{

type zeroGradient;
}




Walls
{
type zeroGradient;
}

Box
{
type fixedValue;
value uniform 1;
}
Boxout
{
type fixedValue;
value uniform 1;
}

}
T
internalField uniform 2.75;

boundaryField
{
Inlet
{
type totalTemperature;
gamma 1.4;
T0 uniform 2.75;
}



Walls
{
type zeroGradient;
}
Boxin
{
type zeroGradient;
}
Box
{
type zeroGradient;
}
Boxout
{
type zeroGradient;
}
}
U
internalField uniform (0 0 0);

boundaryField
{
Inlet
{
type zeroGradient;
}



Walls
{
type noSlip;
}

Boxin
{
type slip;

}
Box
{
type zeroGradient;
}
Boxout
{
type zeroGradient;
}
}
This are the thermophysical properties and turbulence prop.

thermoType
{
type hePsiThermo;
mixture pureMixture;
transport const;
thermo hConst;
equationOfState perfectGas;
specie specie;
energy sensibleInternalEnergy;
}


mixture
{
specie
{
nMoles 1;
molWeight 8314.472;
}
thermodynamics
{
Cp 29.1;
Hf 0;
}
transport
{
mu 0;
Pr 1;
}
}
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

simulationType laminar;


// ************************************************** *********************** //
And now controldict,fv schemes and solution
application rhoCentralFoam;

startFrom startTime;

startTime 0;

stopAt endTime;

endTime 30;

deltaT 1e-6;

writeControl adjustableRunTime;

writeInterval 0.1;

purgeWrite 0;

writeFormat ascii;

writePrecision 6;

writeCompression off;

timeFormat general;

timePrecision 6;

runTimeModifiable true;

adjustTimeStep yes;

maxCo 0.6;

maxDeltaT 1;

fluxScheme Kurganov;

ddtSchemes
{
default Euler;
}

gradSchemes
{
default Gauss linear;
}

divSchemes
{
default none;
div(tauMC) Gauss linear;
}

laplacianSchemes
{
default Gauss linear corrected;
}

interpolationSchemes
{
default linear;
reconstruct(rho) vanAlbada;
reconstruct(U) vanAlbadaV;
reconstruct(T) vanAlbada;
}

snGradSchemes
{
default corrected;
}
solvers
{
"(rho|rhoU|rhoE)"
{
solver diagonal;
}

U
{
solver smoothSolver;
smoother GaussSeidel;
nSweeps 2;
tolerance 0;
relTol 0;
}

h
{
$U;
tolerance 0;
relTol 0;
}
}


By establishing the pressure difference between inlet and outlet (exterior) I aim to set several nozzle regimes to study the flow but when the pressure difference it's too big I always get this error.


[0] #0 Foam::error:rintStack(Foam::Ostream&) at ??:?
[0] #1 Foam::sigFpe::sigHandler(int) at ??:?
[0] #2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #3 Foam::sqrt(Foam::Field<double>&, Foam::UList<double> const&) at ??:?
[0] #4 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::sqrt<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<doub le, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
[0] #5 ? at ??:?
[0] #6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #7 ? at ??:?
[Aspire-F5-573G:18398] *** Process received signal ***
[Aspire-F5-573G:18398] Signal: Floating point exception (8)
[Aspire-F5-573G:18398] Signal code: (-6)
[Aspire-F5-573G:18398] Failing at address: 0x3e8000047de

Please can someone help me? I'm so stuck and I need to fix my problem. I attach my main folders. Thank you so much. You can see the mesh in the attached image.


Alex

[AleQ95]

Please I really need help with this error.

[0] #0 Foam::error:rintStack(Foam::Ostream&) at ??:?
[0] #1 Foam::sigFpe::sigHandler(int) at ??:?
[0] #2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #3 Foam::sqrt(Foam::Field<double>&, Foam::UList<double> const&) at ??:?
[0] #4 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::sqrt<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<doub le, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
[0] #5 ? at ??:?
[0] #6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #7 ? at ??:?
[Aspire-F5-573G:18398] *** Process received signal ***
[Aspire-F5-573G:18398] Signal: Floating point exception (8)
[Aspire-F5-573G:18398] Signal code: (-6)
[Aspire-F5-573G:18398] Failing at address: 0x3e8000047de