Tree-based grids bring the advantage of using fast Cartesian discretizations, such as finite differences, and the flexibility and accuracy of local mesh refinement. The main challenge is how to adapt the discretization stencil near the interfaces between grid elements of different sizes, which is usually solved by local high-order geometrical interpolations. Most methods usually avoid this by limiting the mesh configuration (usually to graded quadtree/octree grids), reducing the number of cases to be treated locally. In this work, we employ a moving least squares meshless interpolation technique, allowing for more complex mesh configurations, still keeping the overall order of accuracy. This technique was implemented in the HiG-Flow code to simulate Newtonian, generalized Newtonian and viscoelastic fluids flows. Numerical tests and application to viscoelastic fluid flow simulations were performed to illustrate the flexibility and robustness of this new approach.

In this work, we present a numerical study on the development length (the length from the channel inlet required for the velocity to reach $99\%$ of its fully-developed value) of a pressure-driven viscoelastic fluid flow (between parallel plates) modelled by the generalised Phan–Thien and Tanner (gPTT) constitutive equation. The governing equations are solved using the finite-difference method, and, a thorough analysis on the effect of the model parameters $\alpha$ and $\beta$ is presented. The numerical results showed that in the creeping flow limit ($Re=0$), the development length for the velocity exhibits a non-monotonic behaviour. The development length increases with $Wi$. For low values of $Wi$, the highest value of the development length is obtained for $\alpha =\beta =0.5$; for high values of $Wi$, the highest value of the development length is obtained for $\alpha =\beta =1.5$. This work also considers the influence of the elasticity number.

In this work, we implement models that are able to describe complex rheological behaviour (such as shear-banding and elastoviscoplasticity) in the HiGTree/HiGFlow system, which is a recently developed Computational Fluid Dynamics (CFD) software that can simulate Newtonian, Generalised-Newtonian and viscoelastic flows using finite differences in hierarchical grids. The system uses a moving least squares (MLS) meshless interpolation technique, allowing for more complex mesh configurations while still keeping the overall order of accuracy. The selected models are the Vasquez-Cook-McKinley (VCM) model for shear-banding micellar solutions and the Saramito model for viscoelastic fluids with yield stress. Development of solvers and numerical simulations of inertial flows of these models in 2D channels and planar-contraction 4:1 are carried out in the HiGTree/HiGFlow system. Our results are compared with those predicted by two other methodologies: the OpenFOAM-based software RheoTool that uses the Finite-Volume-Method and an in-house code that uses the Vorticity-Velocity-Formulation (VVF). We found an excellent agreement between the numerical results obtained by these three different methods. A mesh convergence analysis using uniform and refined meshes is also carried out, where we show that great convergence results in tree-based grids are obtained thanks to the finite difference method and the meshless interpolation scheme used by the HiGFlow software. More importantly, we show that our methodology implemented in the HiGTreee/HiGFlow system can successfully reproduce rheological behaviour of high interest by the rheology community, such as non-monotonic flow curves of micellar solutions and plug-flow velocity profiles of yield-stress viscoelastic fluids.

The phenomenon of electro-osmosis was studied by performing numerical simulations on the flow between parallel walls and at the nozzle microchannels. In this work, we propose a numerical approximation to perform simulations of vortex formation which occur after the passage of the fluid through an abrupt contraction at the microchannel. The motion of the charges in the solution is described by the Poisson–Nernst–Planck equations and used the generalized finite differences to solve the numerical problem. First, solutions for electro-osmotic flow were obtained for the Phan–Thien/Thanner model in a parallel walls channel. Later simulations for electro-osmotic flow were performed in a nozzle. The formation of vortices near the contraction within the nozzle was verified by taking into account a flow perturbation model.