[s.blaise]

Hi All,

I am trying to build a heat transfer case with chtMultiRegionSimpleFoam and use the codedMixed boundary condition below.

The code compiles fine and produces a .so. However, it seems that it is never called. The default refValue, refGradient and valueFraction are applied, and whatever I change in the code (provided that it compiles) does not affect the boundary value. None of my prints shows in the log.

Is there anything wrong with my boundary? The p boundary file has the following entry:

    Bottom
    {
    
    
        type codedMixed;
        refValue uniform 0;
        refGradient uniform 0;
        name hydrostaticBottomGradient;
        valueFraction uniform 1;
        code
        #{
            printMessage("CODED MIXED IN");
            Info << "CODED MIXED IN" <<  endl;    
            const fvPatch& boundaryPatch = patch();
            const fvBoundaryMesh& boundaryMesh = boundaryPatch.boundaryMesh();
            const fvMesh& mesh = boundaryMesh.mesh();
            const label patchid =  mesh.boundary().findPatchID(this->patch().name());
            const volScalarField& rho(mesh.lookupObject<volScalarField>("rho"));
            const auto& rhoPatchField = rho.boundaryField()[patchid];
            scalarField& val = this->refValue();
            scalarField& grad = this->refGrad();
            scalarField& frac = this->valueFraction();
            const vectorField& Cf = patch().Cf(); // Get centre coordinate (we don't use it actually);
            forAll(Cf, facei) // face loop
            {
                grad[facei] = 9.81 * rhoPatchField[facei];
                val[facei] = 200000.;
                frac[facei] = 1;
                Info << "Set val to: " << val[facei] << " for face " << facei <<  endl;

            }

        #};
        
    }

Thanks for any help!

[oswald]

Have you tried working with something like

this->refGrad() = 9.81 * rho.boundaryField()[patchid];

?


And at the moment you are prescribing a fixedValue only, meaning all values will be set to 200000. Is this intended? To use fixedGradient, you have to set valueFraction to 0.