CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > General Forums > Main CFD Forum

How to implement Robin Boundary Conditions?

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 8, 2013, 18:36
Default How to implement Robin Boundary Conditions?
  #1
New Member
 
Join Date: Feb 2013
Posts: 3
Rep Power: 13
giorgk is on a distinguished road
Hi, I'm currently working with a code that requires to assign the Boundary conditions using the coefficients on the general form of Robin Bc (a*u + b * du/dn = g)

From my understanding I can assign Dirichlet conditions by setting a=1, b=0, and g the desired specified value.

Similarly for Neuman bc I need to set a=0, b=1, and g the desired value.
is this correct so far?

Now lets considered a quadrilateral element
1--------2
| |
| |
3--------4

where there is Neuman boundary defined over the element and dirichlet on the nodes 1 and 3

Since I'm doing groundwater modeling, lets say I have specified hydraulic head for nodes 1 and 3 H=10 m and recharge flux equal to 1e-4 m/day. Lets say that the area of the element is 1e6, then 100 m^3/day are coming into that element. Using FE discretization the flux will split into 25 m^2/day per node.

I assume that for the nodes 2 and 4 the coefficients have to be a=0, b=1, g=25.
But what are the values for the nodes 1 and 3?

Thank you

George
giorgk is offline   Reply With Quote

Old   February 9, 2013, 04:52
Default
  #2
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,893
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by giorgk View Post
Hi, I'm currently working with a code that requires to assign the Boundary conditions using the coefficients on the general form of Robin Bc (a*u + b * du/dn = g)

From my understanding I can assign Dirichlet conditions by setting a=1, b=0, and g the desired specified value.

Similarly for Neuman bc I need to set a=0, b=1, and g the desired value.
is this correct so far?

Now lets considered a quadrilateral element
1--------2
| |
| |
3--------4

where there is Neuman boundary defined over the element and dirichlet on the nodes 1 and 3

Since I'm doing groundwater modeling, lets say I have specified hydraulic head for nodes 1 and 3 H=10 m and recharge flux equal to 1e-4 m/day. Lets say that the area of the element is 1e6, then 100 m^3/day are coming into that element. Using FE discretization the flux will split into 25 m^2/day per node.

I assume that for the nodes 2 and 4 the coefficients have to be a=0, b=1, g=25.
But what are the values for the nodes 1 and 3?

Thank you

George


Assume the point 3 is at x=y=0 coordinates. Assume a bi-linear approximation on 1-2-3-4:

u(x,y) = c0 + c1*x +c2*y +c3*xy

that over the boundary x=0 leads to

du/dx = c1+c3*y

Hence, your Robin Bc along edge 3-1 becomes

a*(c2*y) + b*(c1+c3*y)= g

where a,b, g must be known from your problem.
FMDenaro is offline   Reply With Quote

Old   June 12, 2017, 11:40
Default How to implement Robin Boundary Conditions in groovy?
  #3
New Member
 
Sachin
Join Date: Sep 2016
Location: Poitiers,France
Posts: 17
Rep Power: 10
smodh is on a distinguished road
Hello,

I also would like to implement robin boundary condition in my electrostatic program.
To implement my robin BC in code, I am comparing my problem with heat exchange robin BC.
Final robin BC in heat exchange problem is Tface = f*Tinf + (1-f)* Tcentre ;
f=(1+(k/alpha*delta))

So my question is how to find fraction value (f) in robin BC?

My robin BC is
(E*Np)-(grad*Np)=0
Here E is surfaceScalarField ; Np is volScalarField(both scalar value I am calculate).
So can anyone help me to find out fraction value for robin boundary condition and how to implement in groovy?

Thanks
smodh is offline   Reply With Quote

Reply

Tags
boundary condition, finite element method, robin conditions

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
symmetry boundary conditions in cfx lost.identity CFX 41 May 22, 2013 08:21
CFX13 Post Periodic interface EtaEta CFX 7 December 8, 2011 18:15
hexagonal boundary conditions s_h Main CFD Forum 0 August 9, 2009 23:15
Help with boundary conditions Dan CFX 0 April 3, 2006 12:32
A problem about setting boundary conditions lyang Main CFD Forum 0 September 19, 1999 19:29


All times are GMT -4. The time now is 22:52.