CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Programming & Development (http://www.cfd-online.com/Forums/openfoam-programming-development/)
-   -   implementing dictionary and using values therefrom (http://www.cfd-online.com/Forums/openfoam-programming-development/93200-implementing-dictionary-using-values-therefrom.html)

Linse October 7, 2011 13:15

implementing dictionary and using values therefrom
 
Hello All,

I would like to implement an additional computation into some solvers. Basically, an additional action should be taken if a certain testvalue has been passed.
Basically it works, but only with the threshold value hardcoded within the code.

In order to avoid recompiling the specific solver every time the threshold value is changed, I would like to make a dictionary file such as the "transportProperties" residing in the constant-folder.

How can I get the values inside such files into my solver for comparison reason?

For example it works to write:

if (p[celli] > scalar(12))
{
action that is desired;
}

What do I have to write instead of "scalar(12)" for the routine to take over the value "testvalue" from the dictionary?

What I have working by now (in simplefied form) are following entries:

Within the solver I have:
Code:

forAll(mesh.cells(),celli)
{
if (p[celli] > scalar(12))
{action1;}
else
{action2;}
}

Within my createFields.H I have
Code:

IOdictionary dictionaryName
(
    IOobject
    (
        "dictionaryName",
        runTime.constant(),
        mesh,
        IOobject::MUST_READ,
        IOobject::NO_WRITE
    )
);

dimensionedScalar  testvalue(dictionaryName.lookup("testvalue"));


As I said: The first snippet works fine, the second complains if there is no file "dictionaryName" within the constant-folder and is fine if that is there. But all my tries in taking over a value from the dictionary-file failed so far.
Compilation is fine, though...

andyru January 8, 2012 06:08

Hi,

for example put in constant/myProperties
put in there:
myconst myconst [0 0 0 0 0 0 0] 1.0;

and read in with:
Info<< "\nReading myProperites" << endl;

IOdictionary myProperties
(
IOobject
(
"myProperties",
runTime.constant(),
mesh,
IOobject::MUST_READ,
IOobject::NO_WRITE
)
);
const dimensionedScalar myconst(myProperties.lookup("myconst"));

hope, this helps!

Linse January 8, 2012 18:09

Hi Andyru!

Thanks for the tip, but up to now that did not succeed...
What I had in the createFields.H is:

IOdictionary evaluationProperties
(
IOobject
(
"evaluationProperties",
runTime.constant(),
mesh,
IOobject::MUST_READ,
IOobject::NO_WRITE
)
);
const dimensionedScalar evalp(evaluationProperties.lookup("evalp"));

This far it compiles fine.

But as soon as I put in the evaluation-routine where it should check if the pressure is higher as evalp, it throws a compilation error as below:

forAll(mesh.cells(),celli)
{
if (p[celli] > evalp )
{action1;}
else
{action2;}
}

throws:
error: no match for 'operator>' in 'p.Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::<anonymous>.Foam::DimensionedField <double, Foam::volMesh>::anonymous>.Foam::Field<double>::<a nonymous>,Foam::List<double>::<anonymous>.Foam::UL ist<T>::operator[] [with T = double, Foam::label = int] (celli) > evalp'

I would have guesse it is a simple thing just to find my mistake, but I guess I am always again missing to take the wheel before pushing down the accelerator in front of the wall...

Thanks for any further suggestions!

marupio January 8, 2012 18:56

if (p[celli] > evalp.value() )

Linse January 9, 2012 05:32

Thank you, David!
That did the trick! :-)


All times are GMT -4. The time now is 02:04.