We introduce a nonlinear generalized tensorial Maxwell-type constitutive equation to describe shear-thinning glass-forming fluids, motivated by a recent microscopic approach to the nonlinear rheology of colloidal suspensions. The model captures a nonvanishing dynamical yield stress at the glass transition and incorporates normal-stress differences. A modified lattice-Boltzmann (LB) simulation scheme is presented that includes non-Newtonian contributions to the stress tensor and deals with flow-induced pressure differences. We test this scheme in pressure-driven 2D Poiseuille flow of the nonlinear generalized Maxwell fluid. In the steady state, comparison with an analytical solution shows good agreement. The transient dynamics after startup and cessation of the pressure gradient are studied; the simulation reproduces a finite stopping time for the cessation flow of the yield-stress fluid in agreement with previous analytical estimates.

We present a hybrid lattice Boltzmann algorithm for the simulation of flow glass-forming fluids, characterized by slow structural relaxation, at the level of the Navier-Stokes equation. The fluid is described in terms of a nonlinear integral constitutive equation, relating the stress tensor locally to the history of flow. As an application, we present results for an integral nonlinear Maxwell model that combines the effects of (linear) viscoelasticity and (nonlinear) shear thinning. We discuss the transient dynamics of velocities, shear stresses, and normal stress differences in planar pressure-driven channel flow, after switching on (startup) and off (cessation) of the driving pressure. This transient dynamics depends nontrivially on the channel width due to an interplay between hydrodynamic momentum diffusion and slow structural relaxation.

We review recent progress on a microscopic theoretical approach to describe the nonlinear response of glass-forming colloidal dispersions under strong external forcing leading to homogeneous and inhomogeneous flow. Using mode-coupling theory (MCT), constitutive equations for the rheology of viscoelastic shear-thinning fluids are obtained. These are, in suitably simplified form, employed in continuum fluid dynamics, solved by a hybrid-Lattice Boltzmann (LB) algorithm that was developed to deal with long-lasting memory effects. The combined microscopic theoretical and mesoscopic numerical approach captures a number of phenomena far from equilibrium, including the yielding of metastable states, process-dependent mechanical properties, and inhomogeneous pressure-driven channel flow.

We present a hybrid-lattice Boltzmann (LB) algorithm for calculating the flow of glass-forming fluids that are governed by integral constitutive equations with pronounced nonlinear, non-Markovian dependence of the stresses on the flow history. The LB simulation for the macroscopic flow fields is combined with the mode-coupling theory (MCT) of the glass transition as a microscopic theory, in the framework of the integration-through transients formalism. Using the combined LB-MCT algorithm, pressure-driven planar channel flow is studied for a schematic MCT model neglecting spatial correlations in the microscopic dynamics. The cessation dynamics after removal of the driving pressure gradient shows strong signatures of oscillatory flow both in the macroscopic fields and the microscopic correlation functions.