[mostafa kareem]

Hello foamers ,

I'm currently studying source code of coupledFvPatchField so i will use it in my future research extension


for Discretizing the Laplace operator , we need gradientInternalCoeffs and gradientBoundaryCoeffs

so we have

gamma*surface area*( gradientBoundaryCoeffs + gradientInternalCoeffs*internal field )

My question is how the those coefficients are evaluated for couple boundary condition so that we get the following results as stated in source code

template<class Type>
tmp<Field<Type> > coupledFvPatchField<Type>::gradientInternalCoeffs( ) const
{
return -pTraits<Type>:ne*this->patch().deltaCoeffs();
}


template<class Type>
tmp<Field<Type> > coupledFvPatchField<Type>::gradientBoundaryCoeffs( ) const
{
return -this->gradientInternalCoeffs();
}