log velocity profile with groovyBC - syntax
Hi,
I'm trying to implement a logarithmic velocity profile using groovy BC. I've used this code: Code:
inlet Code:
--> FOAM FATAL ERROR: Can someone advise on what is wrong with the syntax? Thanks. |
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.
|
In case anyone is looking for this information... This works, to put a Karman-Prandtl velocity profile on your inlet:
Code:
inlet |
Hi guys,
I am having some troubles creating my own velocity profile. This is what I use as boundary condition: Quote:
Quote:
|
Quote:
|
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. |
Quote:
|
Quote:
Now I have some convergency problems. I will play with it a little bit, maybe I can make it work. |
Quote:
|
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. |
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 |
thank's
I'll see if it gonna works |
OpenFoam
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 |
Quote:
|
so is that correct?
inlet { type groovyBC; variables "vel=exp(3.31453803168746*log(5)-5.62522091753)*normal();"; valueExpression "vel"; } thank's a lot |
Quote:
Code:
inlet |
I need help
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 |
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) |
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!! |
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 |
All times are GMT -4. The time now is 16:59. |