[ELwardi]

It seems the current version of foam-extend only supports coupling of fvm::grad and fvm::div operators in equations.


I'm in dire need of an efficient fvm::laplacian descritization that works with existing block-coupling mechanisms.



The general form of target equations is:

fvScalarMatrix TEqn
(
     fvm::laplacian(phiT, T) + TSources == fvc::laplacian(phiTs, Ts) 
);


fvScalarMatrix TsEqn
(
     fvm::laplacian(phiTs, Ts) + TsSources 
);

I'm thinking I should add fvm::laplacian(phiTs, Ts) as a coupling term in a fvBlockMatrix. Any help/hints is greatly appreciated.