Files copyright Peter Stewart, University of Glasgow
code for the membrane tension model: semiimplicit_code.m produced in Matlab 2018
calls:
tubelaw.m: matlab function of the tube law
tubex.m: first derivative of the tube law
deriv.m compute derivative matrices using fourth order finite differences
MODEL SETUP: four region model
model parameters
DIMENSIONAL
DP: driving pressure in mmHg
TAP: application time of pressure in s
rho: density of the material
L: length of region 3
DIMENSIONLESS
tenref: reference membrane tension
ten1: tension in region 1
ten2: tension in region 3
dl2: dimensionless length of region 2
dl4: dimensionless length of region 4
pe: external pressure on region 1
pc: external pressure on region 3
OUTPUT: file 'bulge_vein_dp_50_ten_0_0001.dat'
13 columns of data
COL 1: dimensionless time
COL 2: dimensionless area of tube 1 at midpoint of region 1
COL 3: dimensionless area of tube 3 at midpoint of region 3
COL 4: dimensionless flux at junction between regions 1 and 2
COL 5: dimensionless flux at junction between regions 3 and 4
COL 6: dimensionless pressure at junction between regions 1 and 2
COL 7: dimensionless pressure at junction between regions 3 and 4
COL 8: minimum dimensionless area in tube 1
COL 9: minimum dimensionless area in tube 3
COL 10: maximum dimensionless area in tube 1
COL 11: maximum dimensionless area in tube 3
COL 12: area at second gridpoint adjancent to upstream boundary
COL 13: area at second gridpoint adjancent to downstream boundary