Files belonging to the open access publication:
A.S. Haase, R.G.H. Lammertink, Heat and mass transfer over slippery, superhydrophobic surfaces. Phys. Fluids. 28, 42002 (2016).
Description of COMSOL Multiphysics 5.0 models:
Graetz-Nusselt_het_N.mph
Graetz-Nusselt_het_N_v2_caseT.mph
Graetz-Nusselt_het_N_v2_caseM.mph
COMSOL models for the Graetz-Nusselt problem for heterogeneously slippery surfaces. The v2 models are based on Graetz-Nusselt_het_N.mph.
There are no essential differences between these models, except for a different set of boundary conditions.
Graetz-Nusselt_het_N_infinitelylong.m
Graetz-Nusselt_het_N_infinitelylong.mph
The MATLAB code is used to model the specific Greatz-Nusselt problem for infinitely long models using COMSOL 5.0 with MATLAB.
The temperature/concentration profile at the outlet of a domain is used as the inlet boundary condition for a next simulation.
This is possible, as the flow profile is fully periodic and does not depend on the temperature/concentration.
Description of files in the directory GRAETZhet:
Each csv-file in this directory contains a solution generated by COMSOL Multiphysics 5.0,
containing the (x,r)-coordinates with the velocity and temperature field at those coordinates.
It also includes the Graetz and Peclet number and the effective slip length for that solution.
The name of the file consists of the case (C, T or M; see paper), the length of the domain
(for instance L10R means L = 10R), the Peclet number and, if axial diffusion is included in
the simulations, the letter D.
Description of the MATLAB codes:
GRAETZhet.m
Main MATLAB code. Requires GRAETZhet_comsol.m and GRAETZ_hom_b.m to run.
This script can to used to compare the solution(s) of the COMSOL model
and the semi-analytical MATLAB model (see Haase et al., J. Fluid Mech. R3 (2015))
and to combine the results of both models in various plots.
GRAETZhet_comsol.m
MATLAB code to import the data from a csv-file in the directory GRAETZhet.
GRAETZ_hom_b.m
MATLAB code to solve the partial differential equation describing the Graetz-Nusselt problem
for finite homogeneous wall slip (see Haase et al., J. Fluid Mech. R3 (2015)).
GRAETZ_Nuinf.m
Stand-alone MATLAB code for plotting the developed Nusselt number Nu_inf as a function of
the Peclet number. The values for Nu_inf can be obtained using GRAETZhet.m and the Tdata#.csv
files. These Tdata#.csv have been overwritten for each case and are only used to compute Nu_inf
for that case. The Tdata#.csv in the directory GRAETZhet are for Pe = 1000 for either case T
or case M (not known anymore).
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.