[tehache]

to whom it may interest:

Figured out that the problems where due to wrong fluxes calculated in peqn.flux() from the gradients at my pressure B.C.
I changed gradientInternalCoeffs() and gradientBoundaryCoeffs() in my B.C. to return 0, now everything is fine.

[nsghaisas]

Thomas,

Can you upload your boussinesq solver? I am facing the same problem which you had mentioned. I can try the things you have mentioned in the subsequent post in order to fix the problems, but with not much background in C++ and OpenFOAM, it will not be that easy for me. Your boussinesq solver will sureley be a great help.

Many Thanks

Niranjan

[tehache]

Niranjan,

I will have a look at the solver, clean up a little, and can upload it hopefully next week.

No warranty, though

[nsghaisas]

Hi Thomas

Any luck with getting that solver? In the meanwhile, I have gained some familiarity with C++ and OpenFoam and can now try to do the things you have mentioned. It would be of great help if you could answer a few of my queries.

I cant find wallBuoyantPressure BC. Is it that the BC was there in OpenFOAM 1.3, but has been removed from future versions? I am using 1.4 and there is a fixedFluxBuoyantPressure BC in 1.4.

Could you give me some pointers regarding how to go about building that new boundary condition? Also, what does the existing fixedFluxBuoyantPressure BC do?

Thanks for your time.

Niranjan

[tehache]

Niranjan,

sorry for promising too much. I am too busy with other stuff currently, hope to find time to get back to OpenFOAM in December, and then will need to learn and adapt my stuff to OpenFoam 1.5 first. I cant help you, for the moment.