Discussion Closed This discussion was created more than 6 months ago and has been closed. To start a new discussion with a link back to this one, click here.

Buoyancy Flow by Boussinesq Approximation

Michael Rembe Certified Consultant

Please login with a confirmed email address before reporting spam

Hello,
simulating buoyancy flow by aid of Boussinesq approximation I get problems by Fz. In my model there is a laminar navier stokes flow and a mass transfer of salt. density=f(c). I start simulation with a very high viscosity (for brine near water: 1e6 !!!) and a small effect of the Boussinesq approximation (gravfak=1e-4).

Fz=-9,81*gravfak*rhof(cc1,Tconst)

The axi-symetrical model has a fresh water inflow, a brine outflow and a wall with salt dissolution. You can see it in the attached file.

If I reduce viscosity I have to reduce the influence of buoyancy flow by reducing gravfak in the same manner. Cell-Reynold's number is <1, Cell-Peclet number is critical but not the problem.

I think the Boussinesq approximation causes the numerical problems. Is there anybody who has a good idea to modify the buoyancy flow ?

Working with -9,81*gravfak*(rhof(cc1,Tconst)-rho0) doesn't help.

Thank you very much!

Best regards
Michael R.


0 Replies Last Post 6 ott 2009, 04:05 GMT-4
COMSOL Moderator

Hello Michael Rembe

Your Discussion has gone 30 days without a reply. If you still need help with COMSOL and have an on-subscription license, please visit our Support Center for help.

If you do not hold an on-subscription license, you may find an answer in another Discussion or in the Knowledge Base.

Note that while COMSOL employees may participate in the discussion forum, COMSOL® software users who are on-subscription should submit their questions via the Support Center for a more comprehensive response from the Technical Support team.