CFD Online Discussion Forums - Block-Coupling of laplacian term in foam-extend

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)

- OpenFOAM Programming & Development (https://www.cfd-online.com/Forums/openfoam-programming-development/)

- - Block-Coupling of laplacian term in foam-extend (https://www.cfd-online.com/Forums/openfoam-programming-development/215712-block-coupling-laplacian-term-foam-extend.html)

Block-Coupling of laplacian term in foam-extend

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:

Code:

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.