[anthonygay1812]

Hello Everyone

I'm working on a project and I've been trying to get the coded fixed value BC running on the rhoCentralFoam solver however when I try and do this it produces a crash with no error output. This is my logFile

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

Create mesh for time = 0

Reading thermophysical properties

Selecting thermodynamics package 
{
    type            hePsiThermo;
    mixture         pureMixture;
    transport       const;
    thermo          hConst;
    equationOfState perfectGas;
    specie          specie;
    energy          sensibleInternalEnergy;
}

Using dynamicCode for patch sides on field p at line 1030 in "/home/anthonygay1812/OpenFOAM/tutorials/tutorials/compressible/rhoCentralFoam/shockTube/0/p.boundaryField.sides"

Then thats it. The compilation output is in a setfields log file and theres no error. The solver just stops. I've tried debugging with gdb and the code exits with this

[Inferior 1 (process 17977) exited with code 0220]

Which ive looked into but it didnt tell me a whole lot. I'm on WSL and running version 7.

Thank you!

p.orig file

/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  7
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform 0;

boundaryField
{
    sides
    {
        type            codedFixedValue;
        value           $internalField;
        redirectType    HPCFluid;
        code
  #{
                operator==(0.5*sin(3.14*this->db().time().value())/2);

  #};

codeOptions
#{
    -I$(LIB_SRC)/finiteVolume/lnInclude
#};

    //type            fixedValue;
    //value           uniform 1;
        
    }

    empty
    {
        type            empty;
    }
}

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

[HPE]

- Can you try to use the "fixedValue" to see if the same error occurs?
- Can you try to write a simpler "codedBC" to see if the same error occurs?
- Is it possible for you to run the same case in a Linux distro to isolate whether the issue is related to Windows?

[anthonygay1812]

Hi HPE

Thanks for the response

Originally Posted by HPE

- Can you try to use the "fixedValue" to see if the same error occurs?

I have tried fixed value and the program runs successfully, no crash.

Originally Posted by HPE

- Can you try to write a simpler "codedBC" to see if the same error occurs?

The crash does happen with less complicated/more complicated operators, I just now tried with operator==(min(10, 0.1)); and the same crash happened.

Originally Posted by HPE

- Is it possible for you to run the same case in a Linux distro to isolate whether the issue is related to Windows?

I'll get back to you on this, I normally have access to a computational cluster with linux/openFoam installed and was going to try running it on that but I've been having issues connecting remotely. I think its a good thing to try so I'll let you know if I can find a linux machine to test it on.

[anthonygay1812]

Originally Posted by HPE

- Is it possible for you to run the same case in a Linux distro to isolate whether the issue is related to Windows?

Bingo! It did work on a remote Linux machine so its a Windows related crash. Unfortunately I havent been able to compile my custom code related to my project on the cluster with Linux on it so it would be ideal to find a workaround on my windows machine. But, progress!

Thank you!

[HPE]

Good to hear "Bingo" .

I kindly suggest you that if you think this is a bug - please do report it to the Foundation issue tracker, for which you can find the link below.

Thank you.