buoyantBoussinesqSimpleFoam heat flux boundary condition
 

Hi all,

I'm a complete beginner with openfoam. I'm trying to place a heat flux boundary condition with the buoyantBoussinesqSimpleFoam solver.
I have openfoam version 4.

There seems to be a compressible boundary condition, but I can't find one for incompressible solvers.

https://github.com/OpenFOAM/OpenFOAM...hScalarField.H

There seemed to have been one in the previous versions of openfoam for incompressible solvers:

https://github.com/OpenFOAM/OpenFOAM...hScalarField.H


Is there something I'm missing? I've searched online, but some of the explanations went way over my head.

Hello Sebastien,

by default, there is no heat flux BC for incompressible solvers. However, you can try to program it. Start form the compressible one, and then start to change it.

Hi agustinvo,

Thank you so much for your reply, and for confirming what I thought. Unfortunately, adding to the source code is too complex for me at this time. As I said, I'm just starting, and it is a lot to learn.

Thanks again.

Hi

do not worry, try to use swak4Foam, and use a groovyBC to set a gradient. That gradient you fix it since

q_w = k * gradT

you define your gradient as q_w/k, and you set the values for the heat flux and thermal conductivity. Please take care of the turbulent conductivity as well. If you have a refined mesh, you should not worry a lot about it, but with this it should be enough.

Hi,

I don't get it when it was available with the previous version what led to the removal of that boundary condition from the latest version OF-4.x.

I am also trying to use heatflux[W/m^2] BC for my problem,which I sucessfully used till OF2.4 but now it seems it's not there.
Any clue if it is avaliable in the development version? I am not sure how to proceed to code it myself hence asking.

Thank you

Hi Agustin,

I have the same problem as Sebastien, and I guess I will try to create this BC. However, I am not familiarized with swak4Foam either with groovyBC. After some research, what I found on the net is too vague for me (I am beginning to learn OpenFOAM). Could you please be more explicit about how to implement this incompressible heat flux boundary condition (For laminar flow it possible)?

Thanks in advance.