  Despite the recent interest in the discontinuous shear-thickening (DST)
behaviour, few computational works tackle the rich hydrodynamics of these
fluids. In this work, we present the first implementation of a microstructural
DST model in Smoothed Particle Hydrodynamic (SPH) simulation. The scalar model
was implemented in an SPH scheme and tested in two flow geometries. Three
distinct ratios of local to non-local microstructural effects were probed:
weak, moderate, and strong non-locality. Strong and moderate cases yielded
excellent agreement with flow curves constructed via the Wyart--Cates (WC)
model, with the moderate case exhibiting banding patterns. Weak non-locality
produced stress-splitting instability, resulting in discontinuous stress fields
and poor agreement with the WC model. The mechanism of the stress-splitting has
been explored and contextualised by the interaction of local microstructure
evolution and the stress-control scheme. Velocity profiles obtained in body
force-driven channel flow were found to be in excellent agreement with the
analytical solution, yielding an upward inflection corresponding to the typical
S-curve. Simulations carried out at increasing driving forces exhibited a
decrease in flow. We showed that even the simple scalar model can capture some
of the key properties of DST materials, laying the foundation for further SPH
study of instabilities and pattern formation.
