CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Pre-Processing

Table defined BC

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 15, 2012, 14:32
Default Table defined BC
  #1
Senior Member
 
Daniele Vicario
Join Date: Mar 2009
Location: Novara, Italy
Posts: 142
Rep Power: 17
danvica is on a distinguished road
Is it possible to define a patch using its caracteristic data ?

For example, on an outlet we could know the p/U relation:

p=0.1bar ---> U=0.1l/s
p=0.5bar ---> U=0.4l/s
p=1.0bar ---> U=0.7l/s
....

It would be usefull in complex system. You can avoid simulating external part already characterized.

Can this be implemented ?

Daniele
danvica is offline   Reply With Quote

Old   March 15, 2012, 14:58
Default
  #2
Senior Member
 
Daniele Vicario
Join Date: Mar 2009
Location: Novara, Italy
Posts: 142
Rep Power: 17
danvica is on a distinguished road
Sorry to make questions before deep searching.

I think the best thing is to use groovyBC in following way.

Let's say we approximate the p/U relation with:

U=p/10+5 (note: just a random linear relation)

and the outlet flow is just in the x direction,
then:

For 0/U:

Code:
outlet          
{
     type            groovyBC;
     variables "Ux=p/10+5";
     valueExpression "vector(Ux,0,0)";
     value           uniform (0 0 0);
}
Is it correct ? Is there a better way ? Is there any stabilty problem ?

Daniele
danvica is offline   Reply With Quote

Old   March 15, 2012, 15:38
Default
  #3
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 4,225
Rep Power: 51
gschaider will become famous soon enoughgschaider will become famous soon enough
Quote:
Originally Posted by danvica View Post
Sorry to make questions before deep searching.

I think the best thing is to use groovyBC in following way.

Let's say we approximate the p/U relation with:

U=p/10+5 (note: just a random linear relation)

and the outlet flow is just in the x direction,
then:

For 0/U:

Code:
outlet          
{
     type            groovyBC;
     variables "Ux=p/10+5";
     valueExpression "vector(Ux,0,0)";
     value           uniform (0 0 0);
}
Is it correct ? Is there a better way ? Is there any stabilty problem ?

Daniele
Specifying the velocity on the outlet is always a bit problematic but if you leave it another boundary with a Neuman-condition you might be fine

In the variables-string a ; is missing

The valueExpression might be more elegant if you write it like "Ux*normal()". That way it works no matter how the boundary is oriented (not 100% sure. You might night a - there)

Concerning your original posting: expressions in swak4Foam can have lookup-tables where you specify a function by value pairs and values are interpolated between these. For an exmple on the usage see the wobbler-case in the examples
gschaider is offline   Reply With Quote

Old   March 15, 2012, 15:40
Default
  #4
Senior Member
 
calim_cfd's Avatar
 
mauricio
Join Date: Jun 2011
Posts: 172
Rep Power: 17
calim_cfd is on a distinguished road
it's hard, at least of me, to say.. but is your simulation ss or transient? i kinda of a need a very specific reason b4 implementing things.. so i can't tell..

as for stability.,, i guess you cannot fix pressure and flowrate, U, at the same patch and model.... so i guess you got pick one var to fix..
__________________
Best Regards
/calim

"Elune will grant us the strength"
calim_cfd is offline   Reply With Quote

Old   March 15, 2012, 17:17
Default
  #5
Senior Member
 
Daniele Vicario
Join Date: Mar 2009
Location: Novara, Italy
Posts: 142
Rep Power: 17
danvica is on a distinguished road
Many thanks ! Sure, I know fixing an outlet U is problematic. I'll make some simple bc tests and report here.

daniele
danvica 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
[ICEM] How can I define different zones in ICEM? llrr ANSYS Meshing & Geometry 14 February 12, 2017 14:44
Run time Selection Mechanism - Some help required to understand jaswi OpenFOAM Programming & Development 3 October 29, 2015 14:42
OpenFOAM static build on Cray XT5 asaijo OpenFOAM Installation 9 April 6, 2011 13:21
OpenFOAM13 for Mac OSX Darwin 104 hjasak OpenFOAM Installation 70 September 24, 2010 06:06
Table Manipulation furkan Siemens 0 April 4, 2006 11:23


All times are GMT -4. The time now is 13:17.