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

Shear stress boundary condition at free surface

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   September 2, 2019, 07:34
Default Shear stress boundary condition at free surface
  #1
New Member
 
Join Date: Apr 2019
Posts: 12
Rep Power: 4
yummigummi is on a distinguished road
Hello everyone,


I am simulating a free rising bubble in a viscoelastic fluid with a moving mesh solver. I implemented the viscoelastic stress term in the momentum equation aswell as the equations to calculate velocity and pressure boundary conditions at the bubble interface. Since the newtonian stress can be directly expressed with viscosities and velocities, this leaves the just a pressure and a velocity field in the original solver. With boundary conditions at the free surface as follows:

Inside
velocity: fixed value
pressure: fixed gradient


Outside
velocity: fixed gradient
pressure: fixed value


While the viscoelastic stress can't be calculated directly from the velocity a PDE is solved returning a tau field for which I now have to find boundary conditions for.


What would be boundary conditions for my tau field at free surface?

Why were those boundary conditions for velocity and pressure chosen?
Is it done similar to inlet and outlet boundary conditions of a pipe for numerical stability? v_inlet=fixedValue, p_inlet=zeroGrad and v_outlet=zeroGrad, p_outlet=fixedValue.

Would I chose the boundary conditions similar to the velocity boundary conditions since the stress is a function of the velocity?
The fixed gradient boundary condition for the velocity is applied as n \bullet \nabla v with n being the normal vector on the free surface. If I would apply the same boundary condition for the shear stress this would be n \bullet \nabla \tauwhere I have to find the gradient of a tensorField, which is not defined in OpenFOAM. I was experimenting with snGrad() but im not sure, if it does what i request.


I am working with foam extend 4.0.



Many thanks for your help, I hope you understand my problem.


Paul
yummigummi is offline   Reply With Quote

Reply

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
Basic Nozzle-Expander Design karmavatar CFX 20 March 20, 2016 09:44
stress free boundary condition Ari Main CFD Forum 2 September 23, 2015 17:32
How to implement shear stress boundary condition? Zubair959 Main CFD Forum 0 June 12, 2013 23:04
how to set up a wall boundary condition according to calculated wall shear stress? gameoverli OpenFOAM Pre-Processing 1 May 21, 2009 09:28
Shear Stress Boundary Condition Prabhas FLUENT 0 December 5, 2006 01:53


All times are GMT -4. The time now is 15:58.