IMRv2/s_run_fdmodel.m at main · InertialMicrocavitationRheometry/IMRv2 · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
clc;
clear;
close;

addpath(genpath('src'));

% equation options
R0 = 100e-6;
Req = R0/10;
tfin = 1.5*R0; %160E-6;
kappa = 1.4;
Lheat = 2.378193575129533e+04;
T8 = 298.15;
rho8 = 998.2;
mu = 1e-2;
Gelastic = 0;
tvector = linspace(0,tfin,256);
radial = 2;
vapor = 1;
collapse = 0;
bubtherm = 0;
medtherm = 0;
masstrans = 0;
stress = 1;
varin = {'progdisplay',0,...
         'radial',radial,...
         'bubtherm',bubtherm,...
         'tvector',tvector,...
         'vapor',vapor,...
         'medtherm',medtherm,...
         'masstrans',masstrans,...
         'method',23,...
         'stress',stress,...
         'collapse',collapse,...
         'mu',mu,...
         'g',Gelastic,...
         'lambda1',1e-7,...
         'lambda2',0,...
         'alphax',1e-3,...
         'r0',R0,...
         'req',Req,...
         'kappa',kappa,...
         't8',T8,...
         'rho8',rho8};

[tfd,Rfd,Rfddot,Pfd,Tfd,Tmfd,kvfd] = f_imr_fd(varin{:},'Nt',70,'Mt',70);
% [tsp,Rsp,Rspdot,Psp,Tsp,Tmsp,kvsp] = f_imr_spectral(varin{:},'Nt',10,'Mt',10);

figure(1)
hold on;
plot(tfd,Rfd,'b-');
% plot(tsp,Rsp,'r^');
ylim([0 1]);