CFD Online Logo CFD Online URL
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

LinkBack Thread Tools Search this Thread Display Modes
Old   September 2, 2019, 06:34
Default Shear stress boundary condition at free surface
New Member
Join Date: Apr 2019
Posts: 12
Rep Power: 5
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:

velocity: fixed value
pressure: fixed gradient

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.

yummigummi is offline   Reply With Quote


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 08:44
stress free boundary condition Ari Main CFD Forum 2 September 23, 2015 16:32
How to implement shear stress boundary condition? Zubair959 Main CFD Forum 0 June 12, 2013 22:04
how to set up a wall boundary condition according to calculated wall shear stress? gameoverli OpenFOAM Pre-Processing 1 May 21, 2009 08:28
Shear Stress Boundary Condition Prabhas FLUENT 0 December 5, 2006 00:53

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