A numerical scheme for non-Newtonian fluids and plastic solids under the GPR model

Abstract A method for modeling non-Newtonian fluids (dilatants and pseudoplastics) by a power law under the Godunov-Peshkov-Romenski model is presented, along with a new numerical scheme for solving this system. The scheme is also modified to solve the corresponding system for power-law elastoplastic solids. The scheme is based on a temporal operator splitting, with the homogeneous system solved using a finite volume method based on a WENO reconstruction, and the temporal ODEs solved using an analytical approximate solution. The method is found to perform favorably against problems with known exact solutions, and numerical solutions published in the open literature. It is simple to implement, and to the best of the authors' knowledge it is currently the only method for solving this modified version of the GPR model.

[1]  A. Almgren,et al.  Highly parallelisable simulations of time-dependent viscoplastic fluid flow with structured adaptive mesh refinement , 2018, Physics of Fluids.

[2]  Vladimir A. Titarev,et al.  Exact and approximate solutions of Riemann problems in non-linear elasticity , 2009, J. Comput. Phys..

[3]  H. S. Green,et al.  A Kinetic Theory of Liquids , 1947, Nature.

[4]  Dimitris Drikakis,et al.  An Eulerian method for multi-component problems in non-linear elasticity with sliding interfaces , 2010, J. Comput. Phys..

[5]  Michael Dumbser,et al.  A unified hyperbolic formulation for viscous fluids and elastoplastic solids , 2016, 1705.02151.

[6]  Michael Dumbser,et al.  Finite volume schemes of very high order of accuracy for stiff hyperbolic balance laws , 2008, J. Comput. Phys..

[7]  B. C. Bell,et al.  p-version least squares finite element formulation for two-dimensional, incompressible, non-Newtonian isothermal and non-isothermal fluid flow , 1994 .

[8]  Randall J. LeVeque,et al.  A study of numerical methods for hyperbolic conservation laws with stiff source terms , 1990 .

[9]  Chi-Wang Shu,et al.  Efficient Implementation of Weighted ENO Schemes , 1995 .

[10]  Eleuterio F. Toro,et al.  Conservative Models and Numerical Methods for Compressible Two-Phase Flow , 2010, J. Sci. Comput..

[11]  E. Toro,et al.  CONSERVATIVE HYPERBOLIC MODEL FOR COMPRESSIBLE TWO-PHASE FLOW WITH DIFFERENT PHASE PRESSURES AND TEMPERATURES , 2004 .

[12]  M. Giles.  An extended collection of matrix derivative results for forward and reverse mode algorithmic dieren tiation , 2008 .

[13]  Ilya Peshkov,et al.  On a pure hyperbolic alternative to the Navier-Stokes equations , 2014 .

[14]  Haran Jackson,et al.  A fast numerical scheme for the Godunov-Peshkov-Romenski model of continuum mechanics , 2017, J. Comput. Phys..

[15]  Rémi Abgrall,et al.  A nominally second-order cell-centered Lagrangian scheme for simulating elastic-plastic flows on two-dimensional unstructured grids , 2013, J. Comput. Phys..

[16]  Michael Dumbser,et al.  A Simple Extension of the Osher Riemann Solver to Non-conservative Hyperbolic Systems , 2011, J. Sci. Comput..

[17]  J. M. Nóbrega,et al.  Analytical solutions for Newtonian and inelastic non-Newtonian flows with wall slip , 2012 .

[18]  Michael Dumbser,et al.  High order ADER schemes for a unified first order hyperbolic formulation of Newtonian continuum mechanics coupled with electro-dynamics , 2016, J. Comput. Phys..

[19]  Eftychios Sifakis,et al.  Computing the Singular Value Decomposition of 3x3 matrices with minimal branching and elementary floating point operations , 2011 .

[20]  Eleuterio F. Toro,et al.  Conservative hyperbolic formulation for compressible two-phase flow with different phase pressures and temperatures , 2007 .

[21]  E. Toro.  Riemann Solvers and Numerical Methods for Fluid Dynamics , 1997 .

[22]  Michael Dumbser,et al.  On Universal Osher-Type Schemes for General Nonlinear Hyperbolic Conservation Laws , 2011 .

[23]  J. Frankel.  Kinetic theory of liquids , 1946 .

[24]  Michael Dumbser,et al.  High order ADER schemes for a unified first order hyperbolic formulation of continuum mechanics: Viscous heat-conducting fluids and elastic solids , 2015, J. Comput. Phys..

[25]  Michael Dumbser,et al.  ADER-WENO finite volume schemes with space-time adaptive mesh refinement , 2012, J. Comput. Phys..

[26]  Michael Dumbser,et al.  Efficient, high accuracy ADER-WENO schemes for hydrodynamics and divergence-free magnetohydrodynamics , 2008, Journal of Computational Physics.

[27]  A. Malyshev,et al.  Hyperbolic equations for heat transfer. Global solvability of the Cauchy problem , 1986 .

[28]  Panagiotis Neofytou,et al.  A 3rd order upwind finite volume method for generalised Newtonian fluid flows , 2005, Adv. Eng. Softw..

[29]  Michael Dumbser,et al.  Theoretical and numerical comparison of hyperelastic and hypoelastic formulations for Eulerian non-linear elastoplasticity , 2018, J. Comput. Phys..

[30]  Dimitris Drikakis,et al.  An Eulerian finite‐volume scheme for large elastoplastic deformations in solids , 2010 .

[31]  E. I. Romenskii.  Hyperbolic equations of Maxwell's nonlinear model of elastoplastic heat-conducting media , 1989 .

[32]  Michael Dumbser,et al.  Cell centered direct Arbitrary-Lagrangian-Eulerian ADER-WENO finite volume schemes for nonlinear hyperelasticity , 2016 .