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