[SMallin]

Hi All!

I am struggling with getting my simulation to run completely through. I have refined the mesh, played with different fvSchemes and fvSolutions. I've also tried adding more faces to refine the mesh further in certain areas but wasn't very successful in running that simulation.

I get this error:
Courant Number mean: 2.88337e-05 max: 0.448034
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
DILUPBiCGStab: Solving for Ux, Initial residual = 0.00375051, Final residual = 1.05253e-14, No Iterations 1
DILUPBiCGStab: Solving for Uy, Initial residual = 0.00375002, Final residual = 1.08457e-14, No Iterations 1
DILUPBiCGStab: Solving for Uz, Initial residual = 0.000551018, Final residual = 4.80447e-09, No Iterations 1
DILUPBiCGStab: Solving for e, Initial residual = 0.000639895, Final residual = 1.63491e-13, No Iterations 1
DILUPBiCGStab: Solving for p, Initial residual = 0.000234893, Final residual = 2.43724e-13, No Iterations 1
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.35978e-13, global = -5.8114e-14, cumulative = 0.000207546
PIMPLE: iteration 2
DILUPBiCGStab: Solving for Ux, Initial residual = 0.0060651, Final residual = 5.04781e-09, No Iterations 1
DILUPBiCGStab: Solving for Uy, Initial residual = 0.00362514, Final residual = 3.78705e-09, No Iterations 1
DILUPBiCGStab: Solving for Uz, Initial residual = 0.000523088, Final residual = 4.79466e-10, No Iterations 1
DILUPBiCGStab: Solving for e, Initial residual = 0.000871031, Final residual = 2.6713e-09, No Iterations 1


--> FOAM FATAL ERROR:
Negative initial temperature T0: -5.55946

From function Foam::scalar Foam::species::thermo<Thermo, Type>::T(Foam::scalar, Foam::scalar, Foam::scalar, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar) const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar) const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar) const) const [with Thermo = Foam::hConstThermo<Foam:erfectGas<Foam::specie> >; Type = Foam::sensibleInternalEnergy; Foam::scalar = double; Foam::species::thermo<Thermo, Type> = Foam::species::thermo<Foam::hConstThermo<Foam:er fectGas<Foam::specie> >, Foam::sensibleInternalEnergy>]
in file /home/pawan/OpenFOAM/OpenFOAM-v1912/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 56.

FOAM aborting

#0 Foam::error:rintStack(Foam::Ostream&)sh: 1: addr2line: not found
addr2line failed
#1 Foam::error::abort()sh: 1: addr2line: not found
addr2line failed
#2 Foam::species::thermo<Foam::hConstThermo<Foam:er fectGas<Foam::specie> >, Foam::sensibleInternalEnergy>::TEs(double, double, double) constsh: 1: addr2line: not found
addr2line failed
#3 Foam::hePsiThermo<Foam:siThermo, Foam:ureMixture<Foam::constTransport<Foam::speci es::thermo<Foam::hConstThermo<Foam:erfectGas<Foa m::specie> >, Foam::sensibleInternalEnergy> > > >::calculate(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, bool)sh: 1: addr2line: not found
addr2line failed
#4 Foam::hePsiThermo<Foam:siThermo, Foam:ureMixture<Foam::constTransport<Foam::speci es::thermo<Foam::hConstThermo<Foam:erfectGas<Foa m::specie> >, Foam::sensibleInternalEnergy> > > >::correct()sh: 1: addr2line: not found
addr2line failed
#5 ?sh: 1: addr2line: not found
addr2line failed
#6 __libc_start_mainsh: 1: addr2line: not found
addr2line failed
#7 ?sh: 1: addr2line: not found
addr2line failed
Aborted (core dumped)

I'm not sure if I downloaded the STL files and if I should use binary rather than ASCII. Or maybe there is something still wrong with my mesh, but when I run checkMesh, it says OK.

I attached my system file to view the blockMesh, and snappyHexMesh. I'm not sure what else to try. I put in an fvOptions with a temperature limit and running that now, but I'm trying to pinpoint what my other error might be.

I appreciate any and all your help!

[SMallin]

Hello Everyone,

I've updated my simulation a bit to try and run further through the motor, but I still end up getting a negative initial temperature error.

I know this usually has something to do with the mesh. After running checkMesh -allTopography -allGeography, I find that I have some concave cells and faces. Those are the main errors and after reading through some other forums, I have found that those can mainly be ignored.

Is there another way to check if your mesh is okay?

Any advice or recommendations is greatly appreciated!! I'm still using sonicFoam and I'm just not sure which part is filled out incorrectly or why I can't get the simulation to finish. I've attached both my blockMesh dict and snappyHexMesh dict to see if I missed anything major.

[SMallin]

Please anyone that can help me, I would greatly appreciate it.

All I can tell with the mesh is concavity. Could that be the only reason that causes the negative initial temperature?

Any help or advice is greatly appreciated!