Create Matrix for pressure correction
 

Hey Everybody,

i try to create a Matrix UEqn for the correction Loop.
The Equatin is:
U-Ra*T*k=-grad(p)
U and k are vectors and T and Ra are scalars.

Is there a way to create a Matrix like
fvVectorMatrix UEqn
(
U
- Ra*T*k
);
but this doesnt work:mad:

Thanks for Help!

You need to use fvm::Sp(U) (or fvm::S for the linear source term, so something like:

fvVectorMatrix UEqn
(
fvm::Sp(U) - Ra*T*k
);

should work. Note that this is just a diagonal matrix.