Non-Newtonian Power-law fluid simulations in rectangular channel

Velocity distribution in rectangular channel for Power law fluid.

Using viscosity: eta(x)=K*(x/gamma0)^(n-1)/gamma0

Solved the PDE equation: https://doc.comsol.com/5.5/doc/com.comsol.help.comsol/comsol_ref_equationbased.23.008.htmlin rectangular channel in COMSOL withc=eta(sqrt(d(u,x)^2+d(u,y)^2), f=0.013333333 [Pa*s]e_a, d_a, alpha, beta, gamma=0gamma0=1[1/s]Boundary conditions on x=0, y=0, x=witdh, y=height is u=0

Output file contains:  Coordinates X, Y [mm] on triangular grid and Velocity u [m/s]

File specific parameters:"H2O_100x096.txt" - n=1, K=0.013333333 [Pa*s], height=1 [mm], width=0.96 [mm]"10pFBS_37C_100x115.txt" - n=0.5, K=0.00541540444218501 [Pa*s], height=1 [mm], width=1.15 [mm]"10pFBS_097x092.txt" - n=0.32, K=0.00541540444218501 [Pa*s], height=0.97 [mm], width=0.92 [mm]