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

Put "calculated" boundary condition

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 17, 2010, 15:07
Default Put "calculated" boundary condition
  #1
New Member
 
Marcelo Matos Martins
Join Date: Mar 2010
Location: Joinville, Santa Catarina, Brazil
Posts: 10
Rep Power: 16
demaram is on a distinguished road
Hi everyone,

I'm new user in OpenFoam and I need to solve a problem where I have to use a "calculated" boundary condition and I would like to know how I do this. Could somebody help me?

Thank your kind attention.
demaram is offline   Reply With Quote

Old   March 17, 2010, 16:15
Default
  #2
Senior Member
 
sega's Avatar
 
Sebastian Gatzka
Join Date: Mar 2009
Location: Frankfurt, Germany
Posts: 729
Rep Power: 20
sega is on a distinguished road
Quote:
Originally Posted by demaram View Post
Hi everyone,

I'm new user in OpenFoam and I need to solve a problem where I have to use a "calculated" boundary condition and I would like to know how I do this. Could somebody help me?

Thank your kind attention.
What is a "calculated" boundary condition to you?
A function depending on coordinates and/or time, or some kind of data (from a file for examle)?
__________________
Schrödingers wife: "What did you do to the cat? It's half dead!"
sega is offline   Reply With Quote

Old   March 17, 2010, 21:01
Default Put "calculated" boundary condition
  #3
New Member
 
Marcelo Matos Martins
Join Date: Mar 2010
Location: Joinville, Santa Catarina, Brazil
Posts: 10
Rep Power: 16
demaram is on a distinguished road
Hi Sega,

Thank your attention.

Let me explain what I want. I made a code to calculate the velocity field in an extrusion process with the following geometry:

[IMG]file:///C:/DOCUME~1/ADMINI~1/CONFIG~1/Temp/msohtml1/01/clip_image002.gif[/IMG]





























I want to put, on the indicate surfaces, a manner to calculate the tangential velocity. I already put the boundary conditions like “fixedValue”, “fixedGradient”, “zeroGradient”, “symmetryPlane”, “empty” and all the conditions failed because the profile solution gets different mainly on the face number 1 when I do a comparison with the literature solution. Then, I want to put a “calculated” boundary condition and I hope it is work. This boundary condition is a friction boundary and it will be calculated by equation that it is being described on the paper that I put a part of the paper where this equation is being described.


[IMG]file:///C:/DOCUME~1/ADMINI~1/CONFIG~1/Temp/msohtml1/01/clip_image002.gif[/IMG]
[IMG]file:///C:/DOCUME~1/ADMINI~1/CONFIG~1/Temp/msohtml1/01/clip_image002.gif[/IMG]


In the indicated surface, on the above figure, I want to calculate the tangential velocity where the used equation is inside the red rectangle described in the part of the paper above. This is, I have to implement this equation in my solver to calculate the tangential velocity on the boundary surfaces. I don’t know how I can make the implementation in the OpenFoam because I don’t know the commands and I don’t know which the place I have to put this comands.

Could you help me?

Thank so much!!!

















Last edited by demaram; March 17, 2010 at 21:11. Reason: Some figures didn't appear
demaram is offline   Reply With Quote

Old   March 18, 2010, 05:05
Default
  #4
Senior Member
 
sega's Avatar
 
Sebastian Gatzka
Join Date: Mar 2009
Location: Frankfurt, Germany
Posts: 729
Rep Power: 20
sega is on a distinguished road
Well, I can't see your images, but I have learned that you have some kind of algebraic function which describes your boundary condition?!

In this case there is a tool called groovyBC, which can apply functional boundary conditions.
__________________
Schrödingers wife: "What did you do to the cat? It's half dead!"
sega is offline   Reply With Quote

Old   March 18, 2010, 08:25
Default Put "calculated" boundary condition
  #5
New Member
 
Marcelo Matos Martins
Join Date: Mar 2010
Location: Joinville, Santa Catarina, Brazil
Posts: 10
Rep Power: 16
demaram is on a distinguished road
Hi Sega,

Thank your attention.

I’m sorry for the last post. Let me explain what I want. I made a code to calculate the velocity field in an extrusion process. I want to put, on the contact surface, a manner to calculate the tangential velocity in function of the internal velocity close the contact surface. How I said I tried all default boundary condition in OpenFoam but all failed. This boundary condition is a friction boundary and it will be calculated by following equation:

vt = vtp <Φ[1 – f(ε)]>

Where, vt is the tangential velocity on the contact surface, vtp is the internal velocity close on the contact surface, the terms McCauley’s < > brackets are “active” only for Φ>0. The function f(ε) is calculated by:

{ m · (pb/μb) · (hb/|vtp|) for Coulomb’s friction law
f(ε) ={
{ f · (ks/μb) · (hb/|vtp|) for constant friction law

Where, m and f are friction coefficient and friction factor, respectively, ks is the yield shear stress, and μb is the viscosity and pb is the pressure at the boundary node, hb is the distance from the internal point close of the contact surface to the boundary.



This is, I have to implement this equation in my solver to calculate the tangential velocity on the boundary surfaces. I don’t know how I can make the implementation in the OpenFoam because I don’t know the commands and I don’t know which the place I have to put this comands.

Could you help me?

I’m sorry again and thanks so much.
demaram is offline   Reply With Quote

Old   March 18, 2010, 08:50
Default
  #6
Senior Member
 
sega's Avatar
 
Sebastian Gatzka
Join Date: Mar 2009
Location: Frankfurt, Germany
Posts: 729
Rep Power: 20
sega is on a distinguished road
Quote:
Originally Posted by demaram View Post
Could you help me?

Sorry, I'm a little bit confused.
You are talking about contact and boundary surfaces.
What is that?
__________________
Schrödingers wife: "What did you do to the cat? It's half dead!"
sega is offline   Reply With Quote

Old   March 18, 2010, 10:12
Default Put "calculated" boundary condition
  #7
New Member
 
Marcelo Matos Martins
Join Date: Mar 2010
Location: Joinville, Santa Catarina, Brazil
Posts: 10
Rep Power: 16
demaram is on a distinguished road
Hi Sega,

My work has aim to calculate the velocity field in the extrusion process. My code solve the momentuum equation for all geometry, but I wanted that the my code solve the momentuum equation for the internal field and giving me the internal velocity and the tangential velocity, that is, the velocity on the surface should be calculated by the equation that I quoted in last post.

Do you understand? If not ask me again.

Thank you so much for your kind attention.
demaram is offline   Reply With Quote

Old   March 18, 2010, 10:41
Default
  #8
Senior Member
 
sega's Avatar
 
Sebastian Gatzka
Join Date: Mar 2009
Location: Frankfurt, Germany
Posts: 729
Rep Power: 20
sega is on a distinguished road
I'm still not getting it, sorry.
If you want to simulate an extrusion process you are dealing with multiphase simulation?! The multiphase solver will be capable of dealing with the whole velocity field (internal or external - whatever that is to you?? Can you explain?).

Would you be so kind to repost your images from above?
(You can't use pictures from your hard-drive)
Maybe this will help.
__________________
Schrödingers wife: "What did you do to the cat? It's half dead!"
sega is offline   Reply With Quote

Old   March 18, 2010, 14:09
Default Put "calculated" boundary condition
  #9
New Member
 
Marcelo Matos Martins
Join Date: Mar 2010
Location: Joinville, Santa Catarina, Brazil
Posts: 10
Rep Power: 16
demaram is on a distinguished road
Hi Sega

Let me try to explain again.

I made some modification in the icoFoam solver to work with extrusion process in aluminium. I am working with the geometry showed in the first figure on the attached.

I am getting to obtain the internal velocity field right when I compare with literature value. However, on the Face 1 in the first figure on the attached the velocity profile is giving different how you can see in the second figure on the attached that show the velocity contours for a central slice.


I would like to become the solutions more closely each other. How I tried all default boundary condition available in the OpeFoam I want to put a calculated boundary condition, that is, I want to calculate the tangential velocity from this before equation.

Thank your patient.
Attached Files
File Type: pdf CFD_Forum_2_Figura.pdf (50.0 KB, 81 views)
demaram is offline   Reply With Quote

Old   March 18, 2010, 14:14
Default
  #10
Senior Member
 
sega's Avatar
 
Sebastian Gatzka
Join Date: Mar 2009
Location: Frankfurt, Germany
Posts: 729
Rep Power: 20
sega is on a distinguished road
Quote:
Originally Posted by demaram View Post
I want to calculate the tangential velocity from this before equation.
Ok. I think I'm finally getting what you want.
But still sounds like groovyBC to me. This tool can apply boundary conditions from functions.
__________________
Schrödingers wife: "What did you do to the cat? It's half dead!"
sega is offline   Reply With Quote

Old   March 18, 2010, 14:26
Default Put "calculated" boundary condition
  #11
New Member
 
Marcelo Matos Martins
Join Date: Mar 2010
Location: Joinville, Santa Catarina, Brazil
Posts: 10
Rep Power: 16
demaram is on a distinguished road
Ok, Sega

I'm going to see how this groovybc works.

If I have doubt I can write you again.

Thank you so much.
demaram is offline   Reply With Quote

Old   November 26, 2014, 02:58
Default
  #12
New Member
 
zhanglimin
Join Date: Nov 2014
Posts: 1
Rep Power: 0
xiaoma is on a distinguished road
Quote:
Originally Posted by demaram View Post
Hi Sega,

Thank your attention.

Let me explain what I want. I made a code to calculate the velocity field in an extrusion process with the following geometry:

[IMG]file:///C:/DOCUME~1/ADMINI~1/CONFIG~1/Temp/msohtml1/01/clip_image002.gif[/IMG]





























I want to put, on the indicate surfaces, a manner to calculate the tangential velocity. I already put the boundary conditions like “fixedValue”, “fixedGradient”, “zeroGradient”, “symmetryPlane”, “empty” and all the conditions failed because the profile solution gets different mainly on the face number 1 when I do a comparison with the literature solution. Then, I want to put a “calculated” boundary condition and I hope it is work. This boundary condition is a friction boundary and it will be calculated by equation that it is being described on the paper that I put a part of the paper where this equation is being described.


[IMG]file:///C:/DOCUME~1/ADMINI~1/CONFIG~1/Temp/msohtml1/01/clip_image002.gif[/IMG]
[IMG]file:///C:/DOCUME~1/ADMINI~1/CONFIG~1/Temp/msohtml1/01/clip_image002.gif[/IMG]


In the indicated surface, on the above figure, I want to calculate the tangential velocity where the used equation is inside the red rectangle described in the part of the paper above. This is, I have to implement this equation in my solver to calculate the tangential velocity on the boundary surfaces. I don’t know how I can make the implementation in the OpenFoam because I don’t know the commands and I don’t know which the place I have to put this comands.

Could you help me?

Thank so much!!!
















hi,
recently ,I am simulating an extrusion process too .I meet the same problem .I do not know how to set up the friction boundary condition .Have you solved the problem now?
xiaoma 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
Boundary condition problem for open channel flow Andy CFX 9 June 11, 2016 08:20
inlet velocity boundary condition murali CFX 5 August 3, 2012 09:56
Transient outlet boundary condition problem jwillie2000 CFX 1 December 7, 2009 18:07
spanwiperiodic boundary condition for LES Sungho Yoon CFX 5 July 8, 2008 07:21
Help Urgent about changing boundary condition Anjum Naveed FLUENT 7 August 14, 2006 13:25


All times are GMT -4. The time now is 10:05.