|
[Sponsors] |
[swak4Foam] log velocity profile with groovyBC - syntax |
![]() |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
![]() |
![]() |
#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. |
|
![]() |
![]() |
![]() |
![]() |
#2 |
Senior Member
Bernhard
Join Date: Sep 2009
Location: Delft
Posts: 790
Rep Power: 21 ![]() |
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.
|
|
![]() |
![]() |
![]() |
![]() |
#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"; } |
|
![]() |
![]() |
![]() |
![]() |
#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:
|
|||
![]() |
![]() |
![]() |
![]() |
#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 |
|
![]() |
![]() |
![]() |
![]() |
#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. |
||||
![]() |
![]() |
![]() |
![]() |
#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 |
|
![]() |
![]() |
![]() |
![]() |
#8 | |
Senior Member
|
Quote:
Now I have some convergency problems. I will play with it a little bit, maybe I can make it work. |
||
![]() |
![]() |
![]() |
![]() |
#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 |
|
![]() |
![]() |
![]() |
![]() |
#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. |
|
![]() |
![]() |
![]() |
![]() |
#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 |
||
![]() |
![]() |
![]() |
![]() |
#12 |
New Member
Lamia
Join Date: Feb 2013
Posts: 18
Rep Power: 13 ![]() |
thank's
I'll see if it gonna works |
|
![]() |
![]() |
![]() |
![]() |
#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 |
|
![]() |
![]() |
![]() |
![]() |
#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 |
||
![]() |
![]() |
![]() |
![]() |
#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 |
|
![]() |
![]() |
![]() |
![]() |
#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 |
||
![]() |
![]() |
![]() |
![]() |
#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 |
|
![]() |
![]() |
![]() |
![]() |
#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 |
||
![]() |
![]() |
![]() |
![]() |
#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!! |
|
![]() |
![]() |
![]() |
![]() |
#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 |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
[UDF] Relative coordinates in logarithmic velocity profile | cfdworker | Fluent UDF and Scheme Programming | 25 | July 13, 2017 03:44 |
UDF error - parabolic velocity profile - 3D turbine | Zaqie | Fluent UDF and Scheme Programming | 9 | June 25, 2016 19:08 |
[swak4Foam] Problem with fully developed velocity profile groovyBC in 2D axi-sym mesh | harshad88 | OpenFOAM Community Contributions | 7 | July 24, 2014 14:28 |
[swak4Foam] groovyBC error: velocity profile (2D) >> what's wrong? | vitorspadeto | OpenFOAM Community Contributions | 4 | June 19, 2014 15:31 |
[boundary condition] logarithmic velocity profile | cfdworker | FLUENT | 2 | April 17, 2009 23:36 |