I got an physical problem in paste extrusion that I am unable to solve:
A pasty material flows through a conical nozzle. The material follows Hurschel-Bulkley flow model (H-B model): $\tau=\tau_y+K(\dot\gamma)^n$, where $\tau_y$ is yield strength, $K$ is flow consistency and $n$ is flow index (they are all known values). The nozzle geometry is also known (inlet and outlet radius $R_0$ and $R_1$, half cone angle $\alpha$).
Find the analytic relation between inlet pressure $P$ and outlet flow rate $Q$.
Assumptions: (1) material is homogeneous, isotropic, incompressible; (2) laminar, steady state flow, no inertia; (3) No boundary slippery; (4) isothermal, the effects of heat generated from viscous flow can be neglected; (5) Von-Mises yield criterion.
! This is a 2D flow problem due to the symmetry of nozzle. Through literature I am aware of, people have been working on paste extrusion problems mainly via two ways:
One is originally from Snelling's paper dates back to 1960 (Snelling, G.R. and J.F. Lontz, 1960). In their work they used spherical coordinates and assumed the material in the conical nozzle all flows towards the apex, and that the material on a virtual spherical cap with same distance to the apex is flowing at the same speed (so velocity is only dependent on $r$ in spherical coordinates). The assumption was also used by A.P. Metzner (1970). And later, M.J. Adams and R.A. Basterfield (2005) used H-B model and Von-Mises yield criterion to derive the extrusion equation in conical nozzle. However, they seemed to neglect shear flow, and only considered elongational flow. And the assumption basically implies the material at the wall of nozzle flows at the same speed with that at centerline.
Another approach is by F.N. Cogswell (1972) who derived the conical flow using cylindrical coordinates where he considered both shear and elongational flow. But they neglected yield strength of the material which made the problem way easier to solve...
I am by no means an expert in mathematics. I tried to introduce H-B model into Cauchy momentum equation, while ended up with equations I am unable to solve. I am wondering if there is a way to get the analytic expression of flow rate $Q$ in this case?