[slaners]

Hi all,

I need to modify simpleFoam to add energy equation. The problem is I am running openFoam on os x and it is encapsulated in an image used by docker.
I assume adding the equation and corresponding fields won't be too hard once I have access to the solvers folder.
Has anyone managed to do it before ?

Thank you in advance !

[clapointe]

There should already be a version of simpleFoam with the energy equation -- buoyantSimpleFoam. Here's a link to the solver files for OF dev : https://github.com/OpenFOAM/OpenFOAM...yantSimpleFoam.

Caelan

[slaners]

Hi Caelan,

Thanks for the answer. Sorry but I should have specified my case : circular pins between two plates placed in an incoming flow, this flow is not driven by buoyancy. I already tried with rhoSimpleFoam but it turns out to be very unstable when I use the k-omega SST turbulence model. However, I managed to get a stable solution with simpleFoam and therefore would like to add temperature to this solver.

[clapointe]

I'm not sure modifying simpleFoam is the answer then -- the product would essentially be a duplicate of rhoSimpleFoam. Have you tried other turbulence models with rhoSimpleFoam? Will it run without any turbulence modeling?

Caelan

[slaners]

Yes it worked with k-epsilon. I get messages errors involving compressible turbulence modeling with k-omega SST, that is why I tired to switch to an incompressible solver.

[clapointe]

Ok -- I misunderstood. What about trying buoyantBoussinesqSimpleFoam? Just set gravity to be the zero vector.

What errors are you getting when using kOmegaSST?

Caelan

[slaners]

I will try buoyantBoussinesqSimpleFoam.
The error message I get is :

#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::hePsiThermo<Foam:siThermo, Foam:ureMixture<Foam::sutherlandTransport<Foam:: species::thermo<Foam::hConstThermo<Foam:erfectGa s<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::calculate() at ??:?
#4 Foam::hePsiThermo<Foam:siThermo, Foam:ureMixture<Foam::sutherlandTransport<Foam:: species::thermo<Foam::hConstThermo<Foam:erfectGa s<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::correct() at ??:?
#5 ? at ??:?
#6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7 ? at ??:?
Floating point exception

[clapointe]

Looks like the temperature went out of the range the sutherland coefficients are good for -- I'd also make sure your case is set up correctly for using kOmegaSST. A quick scan of the tutorials revealed this as an example : https://github.com/OpenFOAM/OpenFOAM...rofoilNACA0012. You could also check your fvSchemes file against the one used in this tutorial.

Caelan

[slaners]

Thanks ! I used the tutorial you sent me to modify my case and it is now working ! I think that this is mainly due to the Temperature limiter in the fvOptions file.
Thanks again !!