In order to demonstrate the capabilities of GENSMAC we present a number of calculations performed by the GENSMAC code as follows:
Two-dimensional Flows (Newtonian)
Axisymmetric Flows (Newtonian)
Two-dimensional Flows (Non-Newtonian)
where q is the shear rate defined by
q = [2(du/dx)^2 + 2(dv/dy)^2 + (dv/dx + du/dy)^2]^(1/2)
The scaling parameters were the same as for example 1 and the viscosity was defined by
nu_inf = 0.001 m^2/s , nu_0 = 0.01 m^2/s, K = 4, m = 1.0
Figure 11 displays several snapshots taken from this simulation.
nu(q) = Kq^(n-1)
where the power-law was defined by
K = 5.829, n = 0.4415, q_min = 0.1898 s^(-1), q_max = 40.0 s^(-1), rho = 1214 kg/m^3.
Due to the high viscosity, we employed the Crank-Nicolson solver for computing the tilde velocities. Figure 12 displays several snapshots taken from this simulation.