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

[swak4Foam] groovyBC or parabolicVelocity + time

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 14, 2011, 10:12
Default groovyBC or parabolicVelocity + time
  #1
Senior Member
 
Rickard
Join Date: May 2010
Location: Lund, Skåne, Sweden
Posts: 143
Rep Power: 15
Rickard.Solsjo is on a distinguished road
Hi
I was trying to get the hang of groovyBc and the parabolic inlet velocity condition.
Im looking for an opinion on whether one or the other is better combined with a time restriction.
Putting it more clearly
At the moment I use timeVaryingUniformFixedValue with a ramp and outOfBounds clamp.
This is running okay, however I want to change the profile over the area so that it leaves my patch with a certain angle, 21.5degrees.
Can groovyBC and parabolic injection be combined with the time-extent? If so, how?
Appreciate the help!
R
Rickard.Solsjo is offline   Reply With Quote

Old   February 14, 2011, 10:34
Default
  #2
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 Rickard.Solsjo View Post
Hi
I was trying to get the hang of groovyBc and the parabolic inlet velocity condition.
Im looking for an opinion on whether one or the other is better combined with a time restriction.
Putting it more clearly
At the moment I use timeVaryingUniformFixedValue with a ramp and outOfBounds clamp.
This is running okay, however I want to change the profile over the area so that it leaves my patch with a certain angle, 21.5degrees.
Can groovyBC and parabolic injection be combined with the time-extent? If so, how?
Appreciate the help!
R
Yep. Search for "timeline" on http://openfoamwiki.net/index.php/Contrib_groovyBC
gschaider is offline   Reply With Quote

Old   February 14, 2011, 12:43
Default
  #3
Senior Member
 
Rickard
Join Date: May 2010
Location: Lund, Skåne, Sweden
Posts: 143
Rep Power: 15
Rickard.Solsjo is on a distinguished road
Thx
this timeline would simply be e.g.
--> "name of timeLine file"
(t(start)
t2
.
.
.
t(end)
);
?
Rickard.Solsjo is offline   Reply With Quote

Old   February 14, 2011, 13:22
Default
  #4
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 Rickard.Solsjo View Post
Thx
this timeline would simply be e.g.
--> "name of timeLine file"
(t(start)
t2
.
.
.
t(end)
);
?
It's the same format as for timeVaryingUniform. Have a look at the example case (groovyBC) that uses it

Bernhard
gschaider is offline   Reply With Quote

Old   February 15, 2011, 04:16
Default
  #5
Senior Member
 
Rickard
Join Date: May 2010
Location: Lund, Skåne, Sweden
Posts: 143
Rep Power: 15
Rickard.Solsjo is on a distinguished road
Ofc I looked at it. It says something about timevaryinguniform, however it only related to the "usual" conditions. I've been trying to find the original documentation of this without success. As I said above I've used timeVaryingUniformFixedValue before.
The only example of timelines at the contrib wiki are
timelines ();
or sth like
timelines (
file... impulse
...
);
without specifying the impulse input.
Thx
R
Rickard.Solsjo is offline   Reply With Quote

Old   February 15, 2011, 05:20
Default
  #6
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 Rickard.Solsjo View Post
Ofc I looked at it. It says something about timevaryinguniform, however it only related to the "usual" conditions. I've been trying to find the original documentation of this without success. As I said above I've used timeVaryingUniformFixedValue before.
The only example of timelines at the contrib wiki are
timelines ();
or sth like
timelines (
file... impulse
...
);
without specifying the impulse input.
Thx
R
That is in the example files. If you feel it should be on the Wiki page feel free to extend it

Bernhard
gschaider is offline   Reply With Quote

Old   February 15, 2011, 05:24
Default
  #7
Senior Member
 
Rickard
Join Date: May 2010
Location: Lund, Skåne, Sweden
Posts: 143
Rep Power: 15
Rickard.Solsjo is on a distinguished road
I didn't mean to be disrespectful, groovyBC is a wonderful tool.
My only concern was that the timelines input wasn't extended enough, so I dont know what the input file looks like
Regards
R
Rickard.Solsjo is offline   Reply With Quote

Old   February 15, 2011, 07:25
Default one more question..
  #8
Senior Member
 
Rickard
Join Date: May 2010
Location: Lund, Skåne, Sweden
Posts: 143
Rep Power: 15
Rickard.Solsjo is on a distinguished road
Mr GschaiderMay I ask yet another question.
At the moment I just have a mean velocity leaving the patch in the axial direction without any radial components.
I would like to create a boundary condition where the flow is leaving the patch at a certain angle relative to the center axis with a speed thats Ucos(theta) where U is the velocity where the angle theta=0.
My gridstructure is cartesian, so my question is if I have to calculate all of the x,y,z components or is multiplying a vector depending on the angle to the scalar velocity-value simply enough? Can groovyBC handle this?
I hope you understand my question
Regards
R
Rickard.Solsjo is offline   Reply With Quote

Old   February 15, 2011, 08:30
Default
  #9
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 Rickard.Solsjo View Post
I didn't mean to be disrespectful, groovyBC is a wonderful tool.
My only concern was that the timelines input wasn't extended enough, so I dont know what the input file looks like
Regards
R
Not a problem.

Don't know how good groovyBC actually is (I'm biased) but its a lot better than the documentation

As you mentioned above that you'd already used timeVaryingUniform I assumed that you were familiar with the file format for those files

If you think that that the file should be on the Wiki page feel free to add it
gschaider is offline   Reply With Quote

Old   February 15, 2011, 08:42
Default
  #10
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 Rickard.Solsjo View Post
Mr GschaiderMay I ask yet another question.
At the moment I just have a mean velocity leaving the patch in the axial direction without any radial components.
I would like to create a boundary condition where the flow is leaving the patch at a certain angle relative to the center axis with a speed thats Ucos(theta) where U is the velocity where the angle theta=0.
My gridstructure is cartesian, so my question is if I have to calculate all of the x,y,z components or is multiplying a vector depending on the angle to the scalar velocity-value simply enough? Can groovyBC handle this?
I hope you understand my question
Regards
R
I don't quite understand what you're trying to achive. But if it can be expressed with the usual vector operations (keep in mind that the operators & and ^ are inner and cross-product for vectors) then it can be done.

But I only help with problems with groovyBC. Not with vector arithmetics
gschaider is offline   Reply With Quote

Old   February 15, 2011, 10:33
Default
  #11
Senior Member
 
Rickard
Join Date: May 2010
Location: Lund, Skåne, Sweden
Posts: 143
Rep Power: 15
Rickard.Solsjo is on a distinguished road
so basically, it is something like this that i want to achieve

patch jet
{
type groovyBC;

variables: "Uz=150;rmax=0.0005;r=sqrt(pow(pos().x,2)+pow(pos( ).y,2));thetaMax=0.9825;theta=thetaMax*r/rmax";
"(r<=0.5) ? U=Uinlet:U=vector(0,0,0);" <-- is this included into the variables?

value: Uinlet="Uz*cos(theta)*vector(0,0,1)+Uz*sin(theta)* vector(0,1,0)*Uz*sin(theta)*vector(1,0,0):Uinlet=v ector(0,0,0)";

}

I appreciate the help
Rickard.Solsjo is offline   Reply With Quote

Old   March 27, 2018, 02:00
Default combining parabolic inlet velocity profile with time varying velocity vector values
  #12
New Member
 
Jawahar
Join Date: Nov 2016
Posts: 4
Rep Power: 9
jawapasu is on a distinguished road
Quote:
Originally Posted by Rickard.Solsjo View Post
so basically, it is something like this that i want to achieve

patch jet
{
type groovyBC;

variables: "Uz=150;rmax=0.0005;r=sqrt(pow(pos().x,2)+pow(pos( ).y,2));thetaMax=0.9825;theta=thetaMax*r/rmax";
"(r<=0.5) ? U=Uinlet:U=vector(0,0,0);" <-- is this included into the variables?

value: Uinlet="Uz*cos(theta)*vector(0,0,1)+Uz*sin(theta)* vector(0,1,0)*Uz*sin(theta)*vector(1,0,0):Uinlet=v ector(0,0,0)";

}

I appreciate the help
I am trying to do something similar to yours. can you help me out? pls refer to the title of this msg. I have posted a thread.
jawapasu 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
[Other] Contribution a new utility: refine wall layer mesh based on yPlus field lakeat OpenFOAM Community Contributions 58 December 23, 2021 02:36
How to export time series of variables for one point? mary mor OpenFOAM Post-Processing 8 July 19, 2017 10:54
pressure in incompressible solvers e.g. simpleFoam chrizzl OpenFOAM Running, Solving & CFD 13 March 28, 2017 05:49
dynamic Mesh is faster than MRF???? sharonyue OpenFOAM Running, Solving & CFD 14 August 26, 2013 07:47
Upgraded from Karmic Koala 9.10 to Lucid Lynx10.04.3 bookie56 OpenFOAM Installation 8 August 13, 2011 04:03


All times are GMT -4. The time now is 02:39.