[grjmell]

Hi,
I'm trying to implement a logarithmic velocity profile using groovy BC. I've used this 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)";
    }

When I run it with pimpleFoam, I get the following error:

--> 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.

[Bernhard]

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.

[grjmell]

In case anyone is looking for this information... This works, to put a Karman-Prandtl velocity profile on your inlet:

    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";
    }

Greg

[Mojtaba.a]

Hi guys,
I am having some troubles creating my own velocity profile.
This is what I use as boundary condition:

type groovyBC;
variables "yp=pts().y;ph=(yp-2.712)/3.22;arctanph=ph-0.33*pow(ph,3)+0.2*pow(ph,5)-0.1428*pow(ph,7);vel=(arctanph+0.7)/0.19;";
valueExpression "vector(vel,0,0)";
// value uniform (25 0 0);

And here is what I get when i run the case:

--> FOAM FATAL ERROR:
Parser Error for driver PatchValueExpressionDriver at "1.5-9" :"syntax error, unexpected value"
"(yp-2.712)/3.22"
^^^^^
------|

Context of the error:


- From dictionary: /home/onixa/OpenFOAM/onixa-2.1.1/run/sehatSanj/0/U::boundaryField::velocity_inlet.2
Evaluating expression "(yp-2.712)/3.22"


From function parsingValue
in file lnInclude/CommonValueExpressionDriverI.H at line 1028.

FOAM exiting

[gschaider]

Originally Posted by Mojtaba.a

Hi guys,
I am having some troubles creating my own velocity profile.
This is what I use as boundary condition:

And here is what I get when i run the case:

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

[Mojtaba.a]

Originally Posted by gschaider

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

Thank you Gschaider,
I think I understand the problem. So I did something like this:

variables "yp=pts().y;ph=(yp-toPoint(2.712))/toPoint(3.22);arctanph=ph-toPoint(0.33)*pow(ph,3)+toPoint(0.2)*pow(ph,5)-toPoint(0.1428)*pow(ph,7);vel=((arctanph+toPoint(0 .7))/toPoint(0.19);";
valueExpression "vector(vel,0,0)";

and I got this error for valueExpression:

--> FOAM FATAL ERROR:
Parser Error for driver PatchValueExpressionDriver at "1.12" :"syntax error, unexpected value"
"vector(vel,0,0)"

What am I doing wrong? Is it again something about vertices and face centers?
Thank you so much for you help.

[gschaider]

Originally Posted by Mojtaba.a

Thank you Gschaider,
I think I understand the problem. So I did something like this:



and I got this error for valueExpression:




What am I doing wrong? Is it again something about vertices and face centers?
Thank you so much for you help.

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

[Mojtaba.a]

Originally Posted by gschaider

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

Yes it fixed the problem, Thank you Gschaider.
Now I have some convergency problems. I will play with it a little bit, maybe I can make it work.

[gschaider]

Originally Posted by Mojtaba.a

Yes it fixed the problem, Thank you Gschaider.
Now I have some convergency problems. I will play with it a little bit, maybe I can make it work.

Convergence is not my responsibility. Good luck

[LamiaOF2.1]

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.

[gschaider]

Originally Posted by LamiaOF2.1

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.

Something like

vel=min((u_f/K)*log(((pos().y+y0)/y0)+c),5)*normal();

should limit it (min(a,b) returns the minimum of a and b)

All available documentation is linked from the swak-page on the Wiki

[LamiaOF2.1]

thank's
I'll see if it gonna works

[LamiaOF2.1]

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

[gschaider]

Originally Posted by LamiaOF2.1

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

What have you tried and what doesn't work? Be a bit more specific. The expression should look like this: "exp(3.31453803168746*log(5)-5.62522091753)" (more elegant if you use variables)

[LamiaOF2.1]

so is that correct?
inlet
{
type groovyBC;
variables "vel=exp(3.31453803168746*log(5)-5.62522091753)*normal();";
valueExpression "vel";

}
thank's a lot

[gschaider]

Originally Posted by LamiaOF2.1

so is that correct?
inlet
{
type groovyBC;
variables "vel=exp(3.31453803168746*log(5)-5.62522091753)*normal();";
valueExpression "vel";

}
thank's a lot

Well Yeah. But with "elegant" I meant

inlet
    {
        type            groovyBC;
        variables  (
              "A=3.31453803168746;"
              "X=5;"
              "B=-5.62522091753;"
       );
        valueExpression "exp(A*log(X) + B)*normal()";
    }

BTW: that is a constant velocity for the whole patch ... which stretches the definition of "profile"

[LamiaOF2.1]

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

[gschaider]

Originally Posted by LamiaOF2.1

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
<snip>
0.00002, 0.19580
0.00001, 0.14580
0.00000, 0.00000

best regards

Yeah. You can either use timelines (for data that depends on time) or lookuptables (for a more general approach. This is what you'll want here). Have a look at the wobbler-case for a usage example.

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)

[LamiaOF2.1]

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!!

[gschaider]

Originally Posted by LamiaOF2.1

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!!

No you don't have to add anything to the controlDict.

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