Projection-based solver for viscoelastic Stokes flow using Fast Fourier Transforms

Understanding the flow of complex media is relevant for a wide range of research fields and industrial applications. Several numerical approaches exist by which approximate solutions can be determined for the Stokes equations that describe microhydrodynamic flows at the continuum level. However, achieving efficiency and accuracy for an incompressible fluid remains challenging. Here, we present an algorithm for solving the Stokes equations for an Oldroyd-B fluid using Fourier transforms. We gain efficiency by leveraging the 'Fastest Fourier Transform in the West' (FFTW). We validate our approach for the well-characterized four-roll mill, which exhibits nearly singular points of stress at the extensional points of the flow. We capture this divergence and showcase the potential of our method without making the usual diffusive renormalization. We also focus on characterizing the power-law behavior and numerically assess the divergence criterion. Future work will concentrate on active systems, the introduction of moving boundaries, and application to microfluidic devices.