[anthonygay1812]

Hello

I'm trying to implement the following BC

        type               codedFixedValue; 
        value              $internalField;
        redirectType       hybridParticleContinuim; 

        code 
        #{ 
            const fvMesh& fluidMeshRef = patch().boundaryMesh().mesh();
            const volScalarField& pFluid = db().lookupObject<volScalarField>("p"); 
            const volScalarField& pParticle = db().lookupObject<volScalarField>("p_Mean");
            const volScalarField& pFluidOld=pFluid.oldTime();
            Info<< "code accessed";

            //const scalar t = this->db().time().value();
            const label fluidIdx=fluidMeshRef.nCells();
            const label particleIdx=1;
            const scalar w=0.75; 

            operator==(1-w)*pFluidOld[fluidIdx]+w*(0.5*pFluid[fluidIdx]+0.5*pParticle[particleIdx]); 
        #};

Unfortunately the

pFluidOld[fluidIdx]

is triggering the following error

/p.boundaryField.fluidBound:43:28: error: invalid operands of types ‘void’ and ‘const double’ to binary ‘operator*’
/opt/openfoam7/wmake/rules/General/transform:25: recipe for target 'Make/linux64GccDPInt32Opt/fixedValueFvPatchFieldTemplate.o' failed
make: *** [Make/linux64GccDPInt32Opt/fixedValueFvPatchFieldTemplate.o] Error 1


--> FOAM FATAL IO ERROR

I can't quite figure out how to approach this. I have figured out that its pFluidOld[fluidIdx] thats returning a void value. Is there a better way to get the pressure from the previous timestep?

Any help is appreciated thank you!

[anthonygay1812]

So it turns out my initial guess was way off and I simply had a syntax error in the operator statement, the correct line is

            operator==((1-w)*pFluidOld[fluidIdx]+w*(0.5*pFluid[fluidIdx]+0.5*pParticle[particleIdx]));

Notice the parenthesis. Unfortunately now my code is crashing when it tries to access the BC

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

Create mesh for time = 0

Create mesh for time = 0

Meshes created
Reading thermophysical properties

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

Using dynamicCode for patch fluidBound on field p at line 30 in "/home/anthonygay1812/OpenFOAM/Working/finalShockTube/0/fluid/p.boundaryField.fluidBound"

I have run into this sort of issue before and it happened when I had not properly initialized with the setFields command. So my guess is that something within the BC code may not accessible.