CFD Online Discussion Forums - RegionCoupling

Dear All

Im trying to solve my previous case which was based on chtmultiregionfoam with conjucateheatfoam.
my boundary at the interface of solid and fluid looks like:

-K1grad(T1)=Ksgrad(Ts)+sigma

I dont use coorectboundarycondition for K . so I dont have harmonic interpolation of K. and for T I am doing the interpolation similar to energyRegionCoupledFvPatchScalarField in openfoam 2.2

Code:

 Field<Type> fOwn = this->patchInternalField();

    Field<Type> fNei = this->patchNeighbourField();

    Field<Type> fOwnsigma = this->patchsigmaField();

    Field<Type> fNeisigma = this->patchsigmaNeighbourField();







    Field<Type> Alpha=this->patchDTField();

    scalarField fAlpha=mag(Alpha);

deltas(patch.nf() & patch.delta());     const fvPatch& patch = regionCoupledPatch_.patch();

    const scalarField alphaDelta(fAlpha/deltasii);

    Field<Type> nbrAlpha=this->patchDTNeighbourField();

    scalarField fnbrAlpha=mag(nbrAlpha);

    const scalarField nbrAlphaDelta(fnbrAlpha/nbrDeltasii);

    scalarField weightsd(fOwn.size());  

    scalarField gweights(fOwn.size());  



    forAll (weightsd, faceI) 

    {

        scalar di = alphaDelta[faceI];

        scalar dni = nbrAlphaDelta[faceI];

        weightsd[faceI] = di/(di + dni);

        gweights[faceI] = 1.0/(di + dni);

    }

if (sigmaflag_==1)

{

    // Do interpolation

   Field<Type>::operator=(weightsd*fOwn + (1.0 - weightsd)*fNei-fOwnsigma*sigmaflag_*gweights-fNeisigma*sigmaflag_*gweights); 

}

which sigma is just defined on surface of the interface. in the correct result I should get some sort of generation due to sigma and the internal field should be effected by that.

the boundary values that I get after modifying regioncoupling is the same as the values at the boundary that I get by chtmultiregion :D.
but its not puting any effect on the internal field !!!!!:confused::(.

I would appreciate any comment

best
mahdi