Cannot use codedFixedValue, Unknown patchField type

[Ananth-Narayan-IITM]

I want the inlet and outlet patch to have parabolic velocity profile. But upon debug I found that DAFoam doesn't support codedFixedValue. To confirm this, here are the debug steps

• I did run simpleFoam using OF2406 using the code (refer below), which ran cleanly.
• Ran DAFoam script with fixedValue which worked
• Replaced with codedFixedValue and crashes. Note that the profile given here is for demonstration, but actual profile varies along y direction

inlet
    {
        type            codedFixedValue;
        value           uniform (0 0 0);
        name            myDynamicInlet;
        code            #{
            const scalar t = this->db().time().value();
            vector Uin = vector(1.0 + 0.5*sin(2*M_PI*t), 0, 0); // Velocity changes based on defined sine function
            operator==(Uin);
        #};
    }

Error occured:

[7]PETSC ERROR: Caught signal number 11 SEGV: Segmentation Violation, probably memory access out of range
[7]PETSC ERROR: Try option -start_in_debugger or -on_error_attach_debugger
[7]PETSC ERROR: or see https://petsc.org/release/faq/#valgrind
[7]PETSC ERROR: or try http://valgrind.org on GNU/linux and Apple Mac OS X to find memory corruption errors
[7]PETSC ERROR: configure using --with-debugging=yes, recompile, link, and run
[7]PETSC ERROR: to get more information on the crash.
[2]PETSC ERROR: or see https://petsc.org/release/faq/#valgrind
[2]PETSC ERROR: or try http://valgrind.org on GNU/linux and Apple Mac OS X to find memory corruption errors
[2]PETSC ERROR: configure using --with-debugging=yes, recompile, link, and run
[2]PETSC ERROR: to get more information on the crash.
[2]PETSC ERROR: Run with -malloc_debug to check if memory corruption is causing the crash.
[5]PETSC ERROR: Caught signal number 11 SEGV: Segmentation Violation, probably memory access out of range
[5]PETSC ERROR: Try option -start_in_debugger or -on_error_attach_debugger
[5]PETSC ERROR: or see https://petsc.org/release/faq/#valgrind
[5]PETSC ERROR: or try http://valgrind.org on GNU/linux and Apple Mac OS X to find memory corruption errors
[5]PETSC ERROR: configure using --with-debugging=yes, recompile, link, and run
[5]PETSC ERROR: to get more information on the crash.
[5]PETSC ERROR: Run with -malloc_debug to check if memory corruption is causing the crash.
[6]
[6]
[6] --> FOAM FATAL IO ERROR: (openfoam-2506)
[6] Unknown patchField type myDynamicInlet for patch type patch

Valid patchField types :

103
(
.
.
codedFixedValue
.
.
)

Request to provide alternate where I can provide varying profile.

[Ananth-Narayan-IITM]

I found a workaround for this, OpenFOAM have setExprBoundaryFieldsDict which writes the values described by expression on certain patch, with this DAFoam can proceed further.

Initially I thought of modifying DAField.Cand got this working but not as expected but this setExprBoundaryFieldsDict is the best work around.

FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      setExprBoundaryFieldsDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
readFields
(
    U
);

pattern
{
    field U;

    expressions
    (
        // inletLeft
        {
            patch   inletLeft;
            target  value;

            expression
            #{
                vector
                (
                    1.0*(1.0 - pow((pos().y() - 3.5)/0.5,2)),
                    0,
                    0
                )
            #};
        }
    );
}