I have already written a D2Q9 Lattice Boltzmann code which uses Immersed Boundary Method for complex geometries . Now, I want to elevate it by adding the ability to simulate the Non-Newtonian fluids.
I finally managed to add Power-Law model to my work.But, for some reason i can only simulate fluids when "n" is between 0.5 and 1.5. if i get out of this range my simulation become unstable. what can be the reason for this behavior?