The newly developed algorithm called the grid-by-grid inversion method is a very convenient method for converting an existing computer code for Newtonian flow simulations to that for viscoelastic flow simulations. In this method, the hyperbolic constitutive equation is split such that the term for the convective transport of stress tensor is treated as a source which is updated iteratively. This allows the stress tensors at each grid point to be expressed in terms of velocity gradient tensor at the same location, and the set of stress tensor components is found after inverting a small matrix at each grid point. To corroborate the robustness and accuracy of the grid-by-grid inversion method, we apply it to the 4 : 1 axisymmetric contraction problem. This algorithm is found to be robust and yields accurate results as compared with other finite volume methods. Any commercial CFD packages for Newtonian flow simulations can be easily converted to those for viscoelastic fluids exploiting the grid-by-grid inversion method.
H. M. Park. "Viscoelastic Flow through an Axisymmetric Contraction Using the Grid-by-Grid Inversion Method." J. Appl. Math. 2012 (SI07) 1 - 15, 2012. https://doi.org/10.1155/2012/420189