[qi.yang@polimi.it]

Hi guys,

I am using funkySetFields to set inlet profiles and initial fields. The following two scripts are regarding the 0_org/Ua and funkySetFieldsDict files, respectively.

  fields
        (
        z Ua
        );
        lookuptables (
        {
        name inletprofilesedvel;
        outOfBounds clamp;
        file "$FOAM_CASE/1d_profil/Ua0.xy";
        }
        {
        name inletprofilesedvelv;
        outOfBounds clamp;
        file "$FOAM_CASE/1d_profil/Ua2.xy";

    inlet
    {
        type            groovyBC;
        refValue        uniform   (0 0 0);
        valueExpression "vector(inletprofileua(pos().z),0,inletprofilewa(pos().z))";
        gradientExpression "vector(0,0,0)";
        fractionExpression "1";
        evaluateDuringConstruction 1;
        variables       "";
        timelines       (
);
        lookuptables    (
{
        name            inletprofilewa;
        file        "$FOAM_CASE/1d_profil/Ua1.xy";
        outOfBounds     clamp;

Ua.xy is

(
(-0.09976763165984251 2.254282019e-06)
(-0.09930371004903199 6.748271668e-06)
(-0.0988414157182435 1.1219279619e-05)
(-0.0983807429595295 1.5667359493e-05)
(-0.0979216860849635 2.0092638377e-05)
(-0.09746423942657101 2.449525303e-05)
(-0.0970083973362585 2.8875340809e-05)
(-0.0965541541857445 3.3233038369e-05)
(-0.09610150436648901 3.7568481483e-05)
(-0.0956504422896245 4.1881805009e-05)
(-0.0952009623858885 4.6173142888e-05)
(-0.09475305910555301 5.0442628141e-05)
(-0.09430672691835701 5.4690392869e-05)
(-0.09386196031343699 5.8916568252e-05)
(-0.09341875379926 6.3121284546e-05)
(-0.092977101903556 6.7304671088e-05)
(-0.0925369991732495 7.1466856289e-05)
(-0.0920984401743925 7.5607967639e-05)
(-0.091661419492097 7.9728131703e-05)
(-0.0912259317304695 8.3827474123e-05)
(-0.09079197151254301 8.7906119615e-05)
(-0.090359533480211 9.1964191975e-05)
(-0.08992861229416149 9.600181407e-05)
(-0.089499202633811 0.000100019107847)
(-0.089071299197239 0.000104016194326)

after I executed ./Allrun, the following warning occurred.

 Modifying field Ua of type volVectorField

 Putting "vector(inletprofilesedvel(pos().z),inletprofilesedvelv(pos().z),inletprofilesedvelw(pos().z))" into field Ua at t = "0" if condition "true" is true
 Keeping patches unaltered

swak4Foam: Allocating new repository for sampledMeshes
--> FOAM Warning :
    From function groovyBCFvPatchField<Type>::groovyBCFvPatchField(const fvPatch& p,const DimensionedField<Type, volMesh>& iF,const dictionary& dict)
    in file groovyBCFvPatchField.C at line 136
    No value defined for Ua on inlet therefore using the internal field next to the patch
 Setting 5487152 of 5487152 cells
 Writing to "Ua"

and I checked 0/Ua, which is shown as follows.

FoamFile
{
    version     2.0;
    format      binary;
    class       volVectorField;
    arch        "LSB;label=32;scalar=64";
    location    "0";
    object      Ua;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 1 -1 0 0 0 0];


internalField   nonuniform List<vector>
5487152
(�/P�4l�?^@^@^@^@^@^@^@^@1^C�bk�e��/P�4l�?^@^@^@^@^@^@^@^@1^C�bk�e��/P�4l�?^@^@^@^@^@^@^@^@1^C�bk�e��/P�4l�?$
~�^M:�����Lyہˁ?^@^@^@^@^@^@^@^@^E~�^M:����сLyہˁ?^@^@^@^@^@^@^@^@^E~�^M:����сLyہˁ?^@^@^@^@^@^@^@^@^H~�^M:����$
~�^M:�����Lyہˁ?^@^@^@^@^@^@^@^@
~�^M:�����Lyہˁ?^@^@^@^@^@^@^@^@^H~�^M:����сLyہˁ?^@^@^@^@^@^@^@^@
6�Ζ?^@^@^@^@^@^@^@^@*ʖx5�[��c�
6��?^@^@^@^@^@^@^@^@*�x5�[��c��
6��?^@^@^@^@^@^@^@^@*�x5�[��c��
6��?^@^@^@^@^@^@^@^@*�x5�[��c��
6��?^@^@^@^@^@^@^@^@*�x5�[��c��
6��?^@^@^@^@^@^@^@^@*�x5�[��c��
6��?^@^@^@^@^@^@^@^@*�x5�[��c��
6��?^@^@^@^@^@^@^@^@*�x5�[��c��
6��?^@^@^@^@^@^@^@^@*�x5�[��c��
6��?^@^@^@^@^@^@^@^@*�x5�[��c��
6��?^@^@^@^@^@^@^@^@*�x5�[��c��
6��?^@^@^@^@^@^@^@^@*�x5�[��c��
6��?^@^@^@^@^@^@^@^@*�x5�[��c��
6��?^@^@^@^@^@^@^@^@*�x5�[��c��
6��?^@^@^@^@^@^@^@^@*�x5�[��c��
6��?^@^@^@^@^@^@^@^@*�x5�[��c��
6��?^@^@^@^@^@^@^@^@*�x5�[��c��
6��?^@^@^@^@^@^@^@^@*�x5�[��c��
6��?^@^@^@^@^@^@^@^@*�x5�[��c��
6��?^@^@^@^@^@^@^@^@*�x5�[��c��

Why the imported inlet values are not numbers? any help will be appreciated.

[qi.yang@polimi.it]

Does anybody have an idea?

[tomf]

Hi,

This is related to the setting writeFormat in your controlDict.

It looks like you have it as binary, so the computer can read it, but you probably can't. You could change it to ascii to check if the result is expected. However using binary actually saves some accuracy and disk space typically, so for the actual run I would advise to keep this as binary.

If you want to keep the original file readable I suggest to make a copy of your 0 folder as 0.org before running funkySetFields. This will also help in correcting mistakes as you can quickly return to the original file.

Have a look at some of the tutorials, they typically use this approach.

Best Regards,
Tom

[qi.yang@polimi.it]

Originally Posted by tomf

Hi,

This is related to the setting writeFormat in your controlDict.

It looks like you have it as binary, so the computer can read it, but you probably can't. You could change it to ascii to check if the result is expected. However using binary actually saves some accuracy and disk space typically, so for the actual run I would advise to keep this as binary.

If you want to keep the original file readable I suggest to make a copy of your 0 folder as 0.org before running funkySetFields. This will also help in correcting mistakes as you can quickly return to the original file.

Have a look at some of the tutorials, they typically use this approach.

Best Regards,
Tom

Thanks a lot, I finally realized this.