Files belonging to the open access publication:
A.S. Haase, J.A. Wood, L.M.J. Sprakel, R.G.H. Lammertink, Inelastic non-Newtonian flow over heterogeneously slippery surfaces, Phys. Rev. E 95, 23105 (2017).
This paper consists of two numerical parts:
1. Non-Newtonian flow (Carreau fluids) over a two-dimensional bubble mattress
Directory: NonNewtonianFlowBubbleMattress
2. Slip flow (power-law fluids) through a smooth tube with a heterogenously slippery wall
Directory: MaximumSlipTubeFlow
Concerning part 1.
BM_exp_nonNflow.mph
BM_exp_nonNflow_negativeProtrusion.mph
COMSOL 5.2 models (positive and negative protrusion angles) to simulate the flow of a Carreau fluid over a bubble mattress.
Excel-files
These files contain the simulation results: the (average) pressure gradients as function of the (average) liquid velocity and/or protrusion angle.
finitedifference_carreau.m (script)
finitedifference_carreau_dimless.m (function)
MATLAB codes to compute the velocity profile (dimensional and non-dimensional) for specified pressure gradient for a Carreau fluid
using the finite difference method.
finitedifference_carreau_evaluation.m
MATLAB code to evaluate the velocity profile and average velocity for a specified pressure gradient, channel height and slip length
by calling finitedifference_carreau_dimless.m. A mat-file is created, which functions as an interpolation map to determine the effective
slip length for the COMSOL results for the bubble mattress.
mat-files
Interpolation maps used by find_sliplength.m to evaluate the simulation results.
find_sliplength.m
MATLAB function to evaluate the main simulation results obtained using BM_exp_nonNflow.mph.
Concerning part 2.
Maximumslip_tubeflow.mph
COMSOL 5.2 model to compute the slip length through a periodic tube with different aspect ratios R/L.
R = 1, whereas L varies from 0.01 to 1000. Also the porosity e = Lg/L varies, as does the power-law index nn.
csv-files
These files contain the simulation results: the slip length as a function of the tube length L.
maximumslip_nn.m
MATLAB code to import the data from csv files to plot the slip length profiles for different power law fluid exponents nn.
Notes to the files:
The MATLAB codes run under MATLAB 8.1.0.604 (R2013a).
All MATLAB codes are annotated.
The MATLAB codes call the code matlab2tikz.m for creating TikZ figures/files. Please comments these lines if you do not have this code available.