[albem]

Originally Posted by gschaider

Both vairants are equivalent for a sane person. But for some shells they are not. And the complain with "theFiles.sh not existing" although it exists.



Right. Just run them and see

Thank you Bernhard for your helpI am defining an extra force term in the N-S equations, for the moment I am defining this force using an analytical expression,and when I wrote in the terminal the command

funkySetFields -field A123 -expression 'vector(-exp(-(pow(pos().x,2)+pow(pos().y,2))),0, 0)' -time 0

it works!, but, if I use a /system/funkySetFieldDict file which contains the following

expressions
(
A123
{
field A123; //field initialise

create true;

expression “vector( -exp(-pow(pos().x,2)+pow(pos().y,2)),0, 0)”; //ring-shaped current

condition “sqrt(pow(pos().x,2)+pow(pos().y,2))<=1.0”;

}
);


and then if I write in the terminal

funkySetFields -time 0

the result is

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.7.x |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 1.7.x-a54150ce1198
Exec : funkySetFields -time 0
Date : Apr 20 2011
Time : 15:35:11
Host : acamapichtli
PID : 6554
Case : /home/albem/OpenFOAM/albem-1.7.x/run/Poiseuille2D_Fnonu2
nProcs : 1
SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Time = 0
Using funkySetFieldsDict



--> FOAM FATAL IO ERROR:
cannot open file

file: /home/albem/OpenFOAM/albem-1.7.x/run/Poiseuille2D_Fnonu2/system/funkySetFieldsDict at line 0.

From function regIOobject::readStream()
in file db/regIOobject/regIOobjectRead.C at line 61.

FOAM exiting


do you know what it could be wrong with the /system/funkySetFieldDict file?

I would like to just modified the file for the analytical expression instead of using the commands in the terminal

by the way I do not hoy to include the

condition “sqrt(pow(pos().x,2)+pow(pos().y,2))<=1.0”;

option in the command line terminal any idea?

thanks in advance

Alberto