CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Free Surface Thermal Convection, wrong velocity in the corner

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 22, 2024, 06:46
Post Free Surface Thermal Convection, wrong velocity in the corner
  #1
New Member
 
Join Date: Jan 2024
Posts: 8
Rep Power: 2
danylo is on a distinguished road
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
Attached Images
File Type: png free Surface Convection.png (180.8 KB, 7 views)
File Type: png pRef - Convection.png (168.7 KB, 8 views)
danylo is offline   Reply With Quote

Old   April 23, 2024, 09:21
Default
  #2
New Member
 
Join Date: Jan 2024
Posts: 8
Rep Power: 2
danylo is on a distinguished road
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.
danylo is offline   Reply With Quote

Reply

Tags
atmosphere, boundary condition b, convection, free surface, totalpressure bc

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
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 19:28
Shear stress boundary condition at free surface yummigummi OpenFOAM Programming & Development 0 September 2, 2019 07:34
Free Surface Convection Model robo CFX 4 September 30, 2015 09:01
UDF for inlet BC, Free Surface Open Channel Flow VOF arshiya4 Fluent UDF and Scheme Programming 3 March 6, 2012 19:13
free surface and dip velocity nazari FLUENT 5 January 14, 2009 05:09


All times are GMT -4. The time now is 01:54.