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.