Three test cases are selected to show the performances and the limitations of the parallel code. In the first example, the flow of a Giesekus fluid in a pressure-driven square-shaped channel is examined. In the second test, the unsteady drag on a sphere is computed and, in the third example, a solid particle in a sheared viscoelastic fluid is considered, thus involving a moving boundary.
The parallel software gives high performances for problems such that the finite element mesh is independent of time. In this case, the coefficients of the Stokes-like system are constant in time and the matrix factorization is computed at the first time step only and kept during the whole simulation, drastically reducing the total computational time. On the other hand, whether the factorization needs to be computed at each time step, the solution time increases by more than one order of magnitude, suggesting to seek for alternative iterative solvers.
Finally, for all the test cases, scalability of parallel software is evaluated using the speed-up normalized with respect to the software execution time on a number of processors greater than 1. This was necessary due to the high memory requirements. Experimental results validate the performance gain of the parallel code as both the problem size and the processor number increase.