|
[Sponsors] |
[swak4Foam] log velocity profile with groovyBC - syntax |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
March 9, 2011, 06:36 |
log velocity profile with groovyBC - syntax
|
#1 |
Member
Join Date: Nov 2010
Posts: 41
Rep Power: 15 |
Hi,
I'm trying to implement a logarithmic velocity profile using groovy BC. I've used this code: Code:
inlet { type groovyBC; variables "u_f=0.05;ks=0.0132;z0=ks/30;K=0.41;vel=2.303*(u_f/K)*log(pos().z/z0)*normal();"; valueExpression "vector(vel,0,0)"; } Code:
--> FOAM FATAL ERROR: Parser Error at "1.11" :"syntax error, unexpected ','" "vector(vel,0,0)" " ^ " From function parsingValue in file PatchValueExpressionDriver.C at line 192. FOAM exiting Can someone advise on what is wrong with the syntax? Thanks. |
|
March 9, 2011, 07:48 |
|
#2 |
Senior Member
Bernhard
Join Date: Sep 2009
Location: Delft
Posts: 790
Rep Power: 22 |
In the expression for vel you have '*normal()', which means that vel is a vector. You can just put "vel" in het valueexpression as far as I know.
|
|
March 18, 2011, 05:19 |
|
#3 |
Member
Join Date: Nov 2010
Posts: 41
Rep Power: 15 |
In case anyone is looking for this information... This works, to put a Karman-Prandtl velocity profile on your inlet:
Code:
inlet { type groovyBC; variables "u_f=0.0239;z0=0.0000333;K=0.41;vel=(u_f/K)*log(pos().z/z0)*normal();"; valueExpression "-vel"; } |
|
March 11, 2013, 15:53 |
|
#4 | ||
Senior Member
|
Hi guys,
I am having some troubles creating my own velocity profile. This is what I use as boundary condition: Quote:
Quote:
|
|||
March 11, 2013, 18:16 |
|
#5 |
Assistant Moderator
Bernhard Gschaider
Join Date: Mar 2009
Posts: 4,225
Rep Power: 51 |
yp is defined on the vertices. Constants are defined on the face centers. Something like "toPoint(2.712)" would fix that. There are examples for that in .... the Examples
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request |
|
March 12, 2013, 04:04 |
|
#6 | |||
Senior Member
|
Quote:
I think I understand the problem. So I did something like this: Quote:
Quote:
What am I doing wrong? Is it again something about vertices and face centers? Thank you so much for you help. |
||||
March 12, 2013, 11:26 |
|
#7 |
Assistant Moderator
Bernhard Gschaider
Join Date: Mar 2009
Posts: 4,225
Rep Power: 51 |
Same problem in reverse: vel is a value defined on the points but valueExpression needs it on the faces. A well-set "toFace" should fix that
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request |
|
March 12, 2013, 14:42 |
|
#8 | |
Senior Member
|
Quote:
Now I have some convergency problems. I will play with it a little bit, maybe I can make it work. |
||
March 12, 2013, 14:47 |
|
#9 |
Assistant Moderator
Bernhard Gschaider
Join Date: Mar 2009
Posts: 4,225
Rep Power: 51 |
Convergence is not my responsibility. Good luck
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request |
|
April 15, 2013, 10:32 |
|
#10 |
New Member
Lamia
Join Date: Feb 2013
Posts: 18
Rep Power: 13 |
Hi Bernhard
first I would thank you for all your help well, I have a velocity profile on my inlet as condition, I used groovy Bc and it works perfectly, but in the experimental the maximum velocity should'nt exceed 5 m/s, how can I respect it and add the value of U max = 5 m/s if I am using the following: inlet { type groovyBC; variables "u_f=0.43;y0=0.0007;y=0.78;c=5.75;K=0.41;vel=( u_f/K)*log(((pos().y+y0)/y0)+c)*normal();"; valueExpression "-vel"; }would you please help me with best regards another thing, I asked you for a tutorial about the groovy bc, can you send me a link of your workshop or of any documentation concerning groovy bc, I need it so much to learn more about it. |
|
April 15, 2013, 14:53 |
|
#11 | |
Assistant Moderator
Bernhard Gschaider
Join Date: Mar 2009
Posts: 4,225
Rep Power: 51 |
Quote:
Code:
vel=min((u_f/K)*log(((pos().y+y0)/y0)+c),5)*normal(); All available documentation is linked from the swak-page on the Wiki
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request |
||
April 16, 2013, 11:16 |
|
#12 |
New Member
Lamia
Join Date: Feb 2013
Posts: 18
Rep Power: 13 |
thank's
I'll see if it gonna works |
|
April 28, 2013, 11:30 |
OpenFoam
|
#13 |
New Member
Lamia
Join Date: Feb 2013
Posts: 18
Rep Power: 13 |
Hi
from tecPlot I've got the equation of my experimental velocity profile , I need to program it with groovy BC, I tried many codes but it dose'nt work with me this is the equation from tecplot Y = e^(A*log(X) + B), A = 3.31453803168746710000E+000 B = -5.62522091753798930000E+000 X= 5 m/s where Y means velocity can anyone help me I need to have the same profile than my experimental case to generate my modelisation with OpenFoam best regards |
|
April 28, 2013, 18:21 |
|
#14 | |
Assistant Moderator
Bernhard Gschaider
Join Date: Mar 2009
Posts: 4,225
Rep Power: 51 |
Quote:
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request |
||
April 29, 2013, 12:45 |
|
#15 |
New Member
Lamia
Join Date: Feb 2013
Posts: 18
Rep Power: 13 |
so is that correct?
inlet { type groovyBC; variables "vel=exp(3.31453803168746*log(5)-5.62522091753)*normal();"; valueExpression "vel"; } thank's a lot |
|
April 29, 2013, 15:07 |
|
#16 | |
Assistant Moderator
Bernhard Gschaider
Join Date: Mar 2009
Posts: 4,225
Rep Power: 51 |
Quote:
Code:
inlet { type groovyBC; variables ( "A=3.31453803168746;" "X=5;" "B=-5.62522091753;" ); valueExpression "exp(A*log(X) + B)*normal()"; }
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request |
||
May 12, 2013, 12:30 |
I need help
|
#17 |
New Member
Lamia
Join Date: Feb 2013
Posts: 18
Rep Power: 13 |
In OpenFoam I need to enter a tabular data to program my velocity profile
is it possible to do it with groovy bc here is the data of both altitude and velocity values please help me for I can pass to the next level of my work y (altitude)[m], U Velocity [m s^-1] 0.87500, 4.15580 0.85000, 4.15580 0.80000, 4.15580 0.75000, 4.15580 0.70000, 4.15580 0.65000, 4.15580 0.60000, 4.15580 0.55000, 4.15580 0.50000, 4.15580 0.45000, 4.15580 0.40187, 4.15580 0.39866, 4.15580 0.39548, 4.15580 0.39231, 4.15580 0.38916, 4.15580 0.38603, 4.15580 0.38291, 4.15580 0.37980, 4.15580 0.37673, 4.15580 0.37367, 4.05580 0.37063, 4.04580 0.36760, 4.03580 0.36459, 4.02580 0.36160, 4.01580 0.35862, 4.00580 0.35566, 3.99580 0.35272, 3.98580 0.34979, 3.97580 0.34689, 3.96580 0.34400, 3.95580 0.34112, 3.94580 0.33827, 3.93580 0.33542, 3.92580 0.33260, 3.91580 0.32979, 3.90580 0.32700, 3.89580 0.32423, 3.88580 0.32147, 3.87580 0.31873, 3.86580 0.31601, 3.85580 0.31330, 3.84580 0.31061, 3.83580 0.30793, 3.82580 0.30527, 3.81580 0.30263, 3.80580 0.30000, 3.79580 0.29739, 3.78580 0.29479, 3.77580 0.29221, 3.76580 0.28965, 3.75580 0.28710, 3.74580 0.28457, 3.73580 0.28205, 3.72580 0.27955, 3.71580 0.27706, 3.70580 0.27459, 3.69580 0.27214, 3.68580 0.26970, 3.67580 0.26727, 3.66580 0.26486, 3.65580 0.26247, 3.64580 0.26009, 3.63580 0.25773, 3.62580 0.25538, 3.61580 0.25305, 3.60580 0.25073, 3.59580 0.24842, 3.58580 0.24614, 3.57580 0.24386, 3.56580 0.24160, 3.55580 0.23936, 3.54580 0.23713, 3.53580 0.23491, 3.52580 0.23271, 3.51580 0.23052, 3.50580 0.22835, 3.49580 0.22619, 3.48580 0.22405, 3.47580 0.22192, 3.46580 0.21980, 3.45580 0.21770, 3.44580 0.21562, 3.43580 0.21354, 3.42580 0.21148, 3.41580 0.20944, 3.40580 0.20741, 3.39580 0.20539, 3.38580 0.20339, 3.37580 0.20140, 3.36580 0.19942, 3.35580 0.19746, 3.34580 0.19551, 3.33580 0.19357, 3.32580 0.19165, 3.31580 0.18974, 3.30580 0.18784, 3.29580 0.18596, 3.28580 0.18409, 3.27580 0.18224, 3.26580 0.18039, 3.25580 0.17856, 3.24580 0.17675, 3.23580 0.17494, 3.22580 0.17315, 3.21580 0.17137, 3.20580 0.16961, 3.19580 0.16786, 3.18580 0.16612, 3.17580 0.16439, 3.16580 0.16267, 3.15580 0.16097, 3.14580 0.15928, 3.13580 0.15760, 3.12580 0.15594, 3.11580 0.15429, 3.10580 0.15265, 3.09580 0.15102, 3.08580 0.14940, 3.07580 0.14780, 3.06580 0.14621, 3.05580 0.14463, 3.04580 0.14306, 3.03580 0.14150, 3.02580 0.13996, 3.01580 0.13842, 3.00580 0.13690, 2.99580 0.13540, 2.98580 0.13390, 2.97580 0.13241, 2.96580 0.13094, 2.95580 0.12948, 2.94580 0.12803, 2.93580 0.12659, 2.92580 0.12516, 2.91580 0.12374, 2.90580 0.12233, 2.89580 0.12094, 2.88580 0.11956, 2.87580 0.11818, 2.86580 0.11682, 2.85580 0.11547, 2.84580 0.11413, 2.83580 0.11280, 2.82580 0.11149, 2.81580 0.11018, 2.80580 0.10888, 2.79580 0.10760, 2.78580 0.10632, 2.77580 0.10506, 2.76580 0.10380, 2.75580 0.10256, 2.74580 0.10133, 2.73580 0.10011, 2.72580 0.09889, 2.71580 0.09769, 2.70580 0.09650, 2.69580 0.09532, 2.68580 0.09415, 2.67580 0.09299, 2.66580 0.09184, 2.65580 0.09069, 2.64580 0.08956, 2.63580 0.08844, 2.62580 0.08733, 2.61580 0.08623, 2.60580 0.08514, 2.59580 0.08405, 2.58580 0.08298, 2.57580 0.08192, 2.56580 0.08087, 2.55580 0.07982, 2.54580 0.07879, 2.53580 0.07776, 2.52580 0.07675, 2.51580 0.07574, 2.50580 0.07474, 2.49580 0.07375, 2.48580 0.07278, 2.47580 0.07181, 2.46580 0.07084, 2.45580 0.06989, 2.44580 0.06895, 2.43580 0.06802, 2.42580 0.06709, 2.41580 0.06618, 2.40580 0.06527, 2.39580 0.06437, 2.38580 0.06348, 2.37580 0.06260, 2.36580 0.06173, 2.35580 0.06086, 2.34580 0.06001, 2.33580 0.05916, 2.32580 0.05832, 2.31580 0.05749, 2.30580 0.05667, 2.29580 0.05585, 2.28580 0.05505, 2.27580 0.05425, 2.26580 0.05346, 2.25580 0.05268, 2.24580 0.05190, 2.23580 0.05114, 2.22580 0.05038, 2.21580 0.04963, 2.20580 0.04889, 2.19580 0.04816, 2.18580 0.04743, 2.17580 0.04671, 2.16580 0.04600, 2.15580 0.04530, 2.14580 0.04460, 2.13580 0.04391, 2.12580 0.04323, 2.11580 0.04256, 2.10580 0.04189, 2.09580 0.04123, 2.08580 0.04058, 2.07580 0.03994, 2.06580 0.03930, 2.05580 0.03867, 2.04580 0.03805, 2.03580 0.03743, 2.02580 0.03682, 2.01580 0.03620, 2.00580 0.03562, 1.99580 0.03504, 1.98580 0.03445, 1.97580 0.03388, 1.96580 0.03331, 1.95580 0.03275, 1.94580 0.03220, 1.93580 0.03165, 1.92580 0.03111, 1.91580 0.03057, 1.90580 0.03004, 1.89580 0.02952, 1.88580 0.02901, 1.87580 0.02850, 1.86580 0.02799, 1.85580 0.02750, 1.84580 0.02701, 1.83580 0.02652, 1.82580 0.02604, 1.81580 0.02557, 1.80580 0.02510, 1.79580 0.02464, 1.78580 0.02419, 1.77580 0.02374, 1.76580 0.02330, 1.75580 0.02286, 1.74580 0.02243, 1.73580 0.02200, 1.72580 0.02158, 1.71580 0.02117, 1.70580 0.02076, 1.69580 0.02036, 1.68580 0.01996, 1.67580 0.01957, 1.66580 0.01918, 1.65580 0.01880, 1.64580 0.01843, 1.63580 0.01805, 1.62580 0.01769, 1.61580 0.01733, 1.60580 0.01697, 1.59580 0.01662, 1.58580 0.01628, 1.57580 0.01594, 1.56580 0.01560, 1.55580 0.01527, 1.54580 0.01495, 1.53580 0.01463, 1.52580 0.01431, 1.51580 0.01400, 1.50580 0.01370, 1.49580 0.01340, 1.48580 0.01310, 1.47580 0.01281, 1.46580 0.01252, 1.45580 0.01224, 1.44580 0.01196, 1.43580 0.01169, 1.42580 0.01142, 1.41580 0.01115, 1.40580 0.01089, 1.39580 0.01063, 1.38580 0.01038, 1.37580 0.01013, 1.36580 0.00989, 1.35580 0.00965, 1.34580 0.00941, 1.33580 0.00918, 1.32580 0.00895, 1.31580 0.00873, 1.30580 0.00851, 1.29580 0.00830, 1.28580 0.00808, 1.27580 0.00788, 1.26580 0.00767, 1.25580 0.00747, 1.24580 0.00727, 1.23580 0.00708, 1.22580 0.00689, 1.21580 0.00670, 1.20580 0.00652, 1.19580 0.00634, 1.18580 0.00617, 1.17580 0.00600, 1.16580 0.00583, 1.15580 0.00566, 1.14580 0.00550, 1.13580 0.00534, 1.12580 0.00518, 1.11580 0.00503, 1.10580 0.00488, 1.09580 0.00474, 1.08580 0.00459, 1.07580 0.00445, 1.06580 0.00432, 1.05580 0.00418, 1.04580 0.00405, 1.03580 0.00392, 1.02580 0.00380, 1.01580 0.00368, 1.00580 0.00356, 0.99580 0.00344, 0.98580 0.00332, 0.97580 0.00321, 0.96580 0.00310, 0.95580 0.00300, 0.94580 0.00289, 0.93580 0.00279, 0.92580 0.00269, 0.91580 0.00260, 0.90580 0.00250, 0.89580 0.00241, 0.88580 0.00232, 0.87580 0.00224, 0.86580 0.00215, 0.85580 0.00207, 0.84580 0.00199, 0.83580 0.00191, 0.82580 0.00184, 0.81580 0.00176, 0.80580 0.00169, 0.79580 0.00162, 0.78580 0.00155, 0.77580 0.00149, 0.76580 0.00143, 0.75580 0.00136, 0.74580 0.00130, 0.73580 0.00125, 0.72580 0.00119, 0.71580 0.00114, 0.70580 0.00108, 0.69580 0.00103, 0.68580 0.00098, 0.67580 0.00094, 0.66580 0.00089, 0.65580 0.00085, 0.64580 0.00080, 0.63580 0.00076, 0.62580 0.00072, 0.61580 0.00068, 0.60580 0.00065, 0.59580 0.00061, 0.58580 0.00058, 0.57580 0.00055, 0.56580 0.00051, 0.55580 0.00048, 0.54580 0.00046, 0.53580 0.00043, 0.52580 0.00040, 0.51580 0.00038, 0.50580 0.00035, 0.49580 0.00033, 0.48580 0.00033, 0.47580 0.00029, 0.46580 0.00027, 0.45580 0.00025, 0.44580 0.00023, 0.43580 0.00021, 0.42580 0.00020, 0.41580 0.00018, 0.40580 0.00017, 0.39580 0.00015, 0.38580 0.00014, 0.37580 0.00013, 0.36580 0.00012, 0.35580 0.00011, 0.34580 0.00010, 0.33580 0.00009, 0.32580 0.00008, 0.31580 0.00007, 0.30580 0.00006, 0.29580 0.00005, 0.27580 0.00004, 0.25580 0.00003, 0.22580 0.00002, 0.19580 0.00001, 0.14580 0.00000, 0.00000 best regards |
|
May 12, 2013, 18:24 |
|
#18 | |
Assistant Moderator
Bernhard Gschaider
Join Date: Mar 2009
Posts: 4,225
Rep Power: 51 |
Quote:
This builds on OF-routines for such files so available file formats is the OF-specific format or CSV (a short documentation is found in the Documentation-folder of the swak-sources)
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request |
||
May 14, 2013, 10:03 |
|
#19 |
New Member
Lamia
Join Date: Feb 2013
Posts: 18
Rep Power: 13 |
Hi again
for introducing a tabular data, I did that: boundaryField { inlet { type groovyBC; fields ( U Yn ); lookuptables ( { name profile; outOfBounds clamp; fileName "$FOAM_CASE/profile.data"; } valueExpression "profile(U)"; ); } is that correct?? and should I add something in the other files as the controlDict or another!! |
|
May 14, 2013, 10:53 |
|
#20 | |
Assistant Moderator
Bernhard Gschaider
Join Date: Mar 2009
Posts: 4,225
Rep Power: 51 |
Quote:
I was just going through my own documentation "what does the fields-entry mean". Anyway. Lookups only work for scalars. This is a boundary condition for U where you want the velocity on the inlet to be a function of the height, right (my interpretation from the dat file you posted). In that case the valueExpression could either be "-normal()*profile(pos().y)" or "vector(profile(pos().y),0,0)" (depends on your case what is more appropriate). If the format of your data file is like in your posting then you'll want to add to the dictionary of profile the entry "readerType csv;". This will require some additional options which are described in swak4FoamReference.org in the Documentation folder
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request |
||
Tags |
groovybc, logarithmic, syntax, velocity profile |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[UDF] Relative coordinates in logarithmic velocity profile | cfdworker | Fluent UDF and Scheme Programming | 25 | July 13, 2017 04:44 |
UDF error - parabolic velocity profile - 3D turbine | Zaqie | Fluent UDF and Scheme Programming | 9 | June 25, 2016 20:08 |
[swak4Foam] Problem with fully developed velocity profile groovyBC in 2D axi-sym mesh | harshad88 | OpenFOAM Community Contributions | 7 | July 24, 2014 15:28 |
[swak4Foam] groovyBC error: velocity profile (2D) >> what's wrong? | vitorspadeto | OpenFOAM Community Contributions | 4 | June 19, 2014 16:31 |
[boundary condition] logarithmic velocity profile | cfdworker | FLUENT | 2 | April 18, 2009 00:36 |