(717b) Mathematical Aspects of Modeling the Rheology of Complex Material | AIChE

(717b) Mathematical Aspects of Modeling the Rheology of Complex Material


Armstrong, M. - Presenter, United States Military Academy
Bull, G., United States Military Academy
Horner, J. S., University of Delaware
Beris, A., University Of Delaware
Recent work using small angle neutron scattering under flow has identified the existence of a microstructure that is dependent on flow conditions in soft colloidal systems. Additionally, over the last twenty years, a large number of empirically based, thixotropic rheological models has been developed that involve, albeit phenomenologically, a single scalar structural parameter, lambda. These models involve the combination of a lambda-dependent decomposition of the shear stress to an elastic and a viscous contribution with a relaxation-based evolution equation for lambda. The rich behavior of the soft colloidal systems is primarily due to the interaction of the structure and the hydrodynamic force and this can be particularly enhanced applied through in transient flow, and/or large amplitude oscillatory shear flow. Additionally there is an analogous family of viscoelastic models, starting with the Oldroyd-8 family that has successfully modeled shear-thinning behavior. More recently there has been “unified” approaches in literature that attempt to combine the best features of both the viscoelastic and thixotropic modelling.

To successfully model thixo-elastic-visco-plastic (TEVP) complex material systems each modelling technique employs between 2 – 4 ODEs, and several algebraic equations, meaning that any optimization strategy attempting to fit models to rheological data must involve robust, and accurate ODE solvers. Due to the stiffness of the models, the ODE solvers must seamlessly be able to incorporate implicit ODE solution techniques. In addition we will discuss several other viscoelastic models like the recent work of Horner et al. (2018), de Souza Mendes and Thompson and the Blackwell and Ewoldt model using a modified Jeffrey’s model with a structural, thixotropic parameter, whose evolution is modeled with an ordinary differential equation in time.

[1] C.R. Lopez-Barron, A.K. Gurnon, A.P.R. Eberle, L. Porcar, and N.J. Wagner. Physical Review 89, 042301 (2014), 1-11.

[2] J. M. Kim. PhD thesis, University of Delaware, 2013.

[3] J. Mewis and N. J. Wagner, Colloidal Suspension Rheology, Cambridge Univ. Press (2012).

[4] P. de Souza Mendes and R. Thompson. Rheol. Acta (2013) 52:673-694.

[5] A. Mujumdar, A. N. Beris and A. B. Metzner. J. Non-Newtonian Fluid Mech. 102 (2003).

[6] C.J. Dimitriou, R.H. Ewoldt and G.H. McKinley. J. Rheol. 571(1), (2013).

[7] K. Dullaert and J. Mewis. Rheol. Acta (2005) 45: 23-32.

[8] K. Dullaert and J. Mewis. JNNFM 139 (2006) 21-30.

[9] M.J. Armstrong, PhD Thesis, University of Delaware (2015).

[10] M. J. Armstrong, A. N. Beris, S. Rogers, N. J. Wagner,J. Rheol. 60, 433 (2016).

[11] Ewoldt and McKinley Rheol. Acta 56, 195 (2017).

[12] Mwasame et al. J. Rheol. 60, 225 (2016).

[13] Rogers “Sequence of Physical Phenomena.” Rheol. Acata (2017).

[14] Blackwell and Ewoldt JNNFM 227 (2016).

[15] Horner, Armstrong, Wagner and Beris. JOR (2018).