copyright Peter Stewart, University of Glasgow
05-07-2019
files for
Self-excited oscillations in a collapsible channel, with applications to retinal venous pulsation
for ANZIAM Journal
the files all run in MATLAB 2013a
globscript_channel_trace_dimensionless.m
is the global script to compute the growth rate of the mode. To do so, it first computes the static configuration using:
static_twocomp_channel.m
from which it then computes the growth rate using
globfun_channel.m
The global code uses a bisection method to compute the critical pc (isolating growth rate with zero real part) for a given pe and then takes a prescribed step in pe and restarts.
This file returns a file
eig_dimensionless_ten_1_rey_1.dat
with 13 columns
1) dimensionless wall tension
2) Reynolds number
3) oscillation frequency
4) external pressure on compartment 1
5) external pressure on compartment 3
6) minimum height of membrane in compartment 1
7) corresponding x position to above
8) maximum height of membrane in compartment 1
9) corresponding x position to above
10) minimum height of membrane in compartment 2
11) corresponding x position to above
12) maximum height of membrane in compartment 2
13) corresponding x position to above
and generates the corresponding static shape of the membrane in
static_channel_dimensionless_run.dat
with 3 columns
1) x positions
2) channel profile in compartment 1
3) channel profile in compartment 3