|
[Sponsors] |
Free Surface Thermal Convection, wrong velocity in the corner |
![]() |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
![]() |
![]() |
#1 |
New Member
Join Date: Jan 2024
Posts: 8
Rep Power: 2 ![]() |
Hello everyone!
I am trying to simulate the freezing of water, using a customized version of rhoPimpleDyMFoam on OpenFOAM 5. Now I am running into a problem right from the start, when trying to simulate the density based heat convection. The cooling comes from the side wall on the right-hand side and the top patch is a surface open to the atmosphere. Now the usual boundary conditions for free surface: p - totalPressure 1e5 U - pressureInletOutletVelocity (0 0 0) result in a weird verlocity field, where the highest velocity is in the top corner on the cold side. I think the reason for this might be, that OpenFOAM treats this patch as an In-/Outlet where the sum of velocitys results in (0 0 0). Meaning incoming flow on one side directly results in outgoing flow on the other side, which does not seem right to me. The further the simulation run goes, the more this flow in the corner develops and disrupts my desired convection flow. I managed to achieve the convection I am aiming for, using a pressure reference cell somewhere in the top layer (1e5) and following BC´s: p: zero Gradient U: slip This runs well in the beginning (see screenshots), but eventually results in a rise in total pressure in the system, followed by velocity instability due to the unchanged pressure of the reference cell. Any idea what boundary conditions I could use to solve these problems? Or if there is any possibility of working with a reference pressure, while simultaneously limiting the system pressure? Thank you all P.S. BC´s for the other walls: p - zeroGradient U - noSlip frontAndBack are empty |
|
![]() |
![]() |
![]() |
![]() |
#2 |
New Member
Join Date: Jan 2024
Posts: 8
Rep Power: 2 ![]() |
Found the reason for this weird behaviour.
It was due to how Openfoam calculates the the gravitational effect of having different densities. Instead of increasing the pressure at the bottom, it lowered the pressure at the top so a gradient with the top boundary was created. Set my upper boundary on the y-axis and now it runs fine. |
|
![]() |
![]() |
![]() |
Tags |
atmosphere, boundary condition b, convection, free surface, totalpressure bc |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
SU2 7.0.7 Built on CentOS 7, parallel computation pyscript mpi exit error? | EternalSeekerX | SU2 | 3 | October 9, 2020 18:28 |
Shear stress boundary condition at free surface | yummigummi | OpenFOAM Programming & Development | 0 | September 2, 2019 06:34 |
Free Surface Convection Model | robo | CFX | 4 | September 30, 2015 08:01 |
UDF for inlet BC, Free Surface Open Channel Flow VOF | arshiya4 | Fluent UDF and Scheme Programming | 3 | March 6, 2012 18:13 |
free surface and dip velocity | nazari | FLUENT | 5 | January 14, 2009 04:09 |