[CHARLES]

Hello,

For some reason I don't think that exp() is doing anything... can you please tell me what you think is wrong?

I have the following code:

argf2 = sqr(k_)/(6.0*nu()*epsilon_);
f2 = 1-2/9*Foam::exp(-1.0*sqr(argf2));

I can tell that argf2 is doing what is supposed to.
However, when I bring up f2 in paraview it just shows up as 1 everywhere.

If I define f2 in the following way:

f2 = 99-2/9*Foam::exp(-1.0*sqr(argf2));

f2 shows up as 99 everywhere.

If I define f2 in the following way:

f2 = -2/9*Foam::exp(-1.0*sqr(argf2));

Then the field f2 shows up as 0, even if it was initialized as a different value.

I have also tried to not use "Foam::" before exp() and it doesn't make a difference.

Additionally, when I look at the f2 that was output to the time directory it looks like this:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.2.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0.001";
    object      f2;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 0 0 0 0 0 0];

internalField   uniform 1;

boundaryField
{
    plate
    {
        type            calculated;
        value           uniform 1;
    }
    back
    {
        type            empty;
    }
    front
    {
        type            empty;
    }
    inlet
    {
        type            calculated;
        value           uniform 1;
    }
    outlet
    {
        type            calculated;
        value           uniform 1;
    }
    top
    {
        type            calculated;
        value           uniform 1;
    }
    symmetry
    {
        type            symmetryPlane;
    }
}


// ************************************************************************* //

when in reality it should be an internalField nonuniform List<scalar>


Thank you!