Hello Everybody...
I am simulating non-newtonian power-law fluid with Multi-Relaxation Time Lattice Boltzmann Method (MRT-LBM).
As you know, the relaxation time (τ) in the non-newtonian fluid is not constant. So when we consider the s8 and s7 array in the S diagonal matrix in the MRT-LBM as (1 / τ), it does not work and an error will occur in the Matlab software as following ...
How Can I fix the error?
Error using horzcat
Dimensions of matrices being concatenated are not consistent.