[vbchris]

Hi all,

I am looking for help with the following code problem.
I'd like the ability to adjust the boundary fields in code without using swak4foam.

I've added the following snippet of code to bubbleFoam (OF v2.0):

label outletPatchID = mesh.boundaryMesh().findPatchID("outletBottom");   
Ua.boundaryField()[outletPatchID] = Uabar/
                    gSum(
                    mesh.magSf().boundaryField()[outletPatchID] *
                    alpha.boundaryField()[outletPatchID]) / 
                    gSum(
                    mesh.magSf().boundaryField()[outletPatchID]
                    );


Info<< nl << "Corrected Boundry Velocity = " << 
                    Uabar/
                    gSum(
                    mesh.magSf().boundaryField()[outletPatchID] *
                    alpha.boundaryField()[outletPatchID]) / 
                    gSum(
                    mesh.magSf().boundaryField()[outletPatchID]
                    )
                    << endl;

Info<< "Current Boundry Field = " << endl << Ua.boundaryField()[outletPatchID] << nl << endl;

However it runs with the following console output:

Courant Number mean: 0.00492361038 max: 0.025444501
deltaT = 0.01

Corrected Boundry Velocity = (0 -0.144231961 0)
Current Boundry Field = 
type            fixedValue;
value           uniform (0 -0.05 0);

UaBar is a vector(0 -0.01 0)

I would have expected both the outputs for current/corrected boundary to be the same value. However, as the debug info shows Ua.boundryField[] doesn't change after I assign it a new value.

Any help would be appreciated.

[vbchris]

Just as a follow-up

I was mimicking this snippet of code also found in bubbleFoam, which allows explicit declaration of boundry field values using the method .boundaryField()[patchi].

// PEqn.H snippit

    forAll(p.boundaryField(), patchi)
    {
        if (isA<zeroGradientFvPatchScalarField>(p.boundaryField()[patchi]))
        {
            phiDraga.boundaryField()[patchi] = 0.0;
            phiDragb.boundaryField()[patchi] = 0.0;
        }
    }