Flow and hydrodynamic shear stress inside a printing needle during biofabrication

• Published in: PloS one 2020-07, Vol.15 (7), p.e0236371-e0236371
• Main Authors: Müller, Sebastian J, Mirzahossein, Elham, Iftekhar, Emil N, Bächer, Christian, Schrüfer, Stefan, Schubert, Dirk W, Fabry, Ben, Gekle, Stephan
• Format: Article
• Language: English
• Subjects: Alginates; Alginates - chemistry; Alginic acid; Algorithms; Biology and life sciences; Chitosan; Chitosan - chemistry; Computational fluid dynamics; Computer simulation; Engineering and technology; Flow rates; Flow velocity; Fluid flow; Fluids; Hydrodynamics; Hydrogels; Hydrogels - chemistry; Incompressible flow; Laboratories; Mathematical models; Methods; Microfluidics; Navier-Stokes equations; Needles; Nozzles; Numerical analysis; Numerical simulations; Parameters; Physical Sciences; Polymer melts; Pressure data; Printing; Printing, Three-Dimensional - instrumentation; Rheological properties; Rheology; Shear flow; Shear rate; Shear Strength; Shear stress; Shear thinning (liquids); Simulation; Thinning; Velocity; Velocity distribution; Viscosity
• ISSN: 1932-6203; 1932-6203
• DOI: 10.1371/journal.pone.0236371

We present a simple but accurate algorithm to calculate the flow and shear rate profile of shear thinning fluids, as typically used in biofabrication applications, with an arbitrary viscosity-shear rate relationship in a cylindrical nozzle. By interpolating the viscosity with a set of power-law functions, we obtain a mathematically exact piecewise solution to the incompressible Navier-Stokes equation. The algorithm is validated with known solutions for a simplified Carreau-Yasuda fluid, full numerical simulations for a realistic chitosan hydrogel as well as experimental velocity profiles of alginate and chitosan solutions in a microfluidic channel. We implement the algorithm in an easy-to-use Python tool, included as Supplementary Material, to calculate the velocity and shear rate profile during the printing process, depending on the shear thinning behavior of the bioink and printing parameters such as pressure and nozzle size. We confirm that the shear stress varies in an exactly linear fashion, starting from zero at the nozzle center to the maximum shear stress at the wall, independent of the shear thinning properties of the bioink. Finally, we demonstrate how our method can be inverted to obtain rheological bioink parameters in-situ directly before or even during printing from experimentally measured flow rate versus pressure data.