Homec4science
Diffusion yearII_reports (master) Browse

rSTICAZZI/physnum/rap5/Approaching.m
5fb2fcbcde85master

Last Change
Raw File

Approaching.m
View Options

%% Setup %%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
repertoire = './'; % Chemin d'acces au code compile (NB: enlever le ./ sous Windows)
executable = 'Exercice5'; % Nom de l'executable (NB: ajouter .exe sous Windows)
input = 'configuration.in'; % Nom du fichier d'entree de base
output = 'm_outs/approaching';
N = 50;
L = 0.1;
kappa = 1e-3;
dt = 0.0001;
mode = 1;
% mode = 0, point
% mode = 1, final grid
% mode = 2, all steps grid
h = L / N;
X = zeros(N+1);
Y = zeros(N+1);
for i=1:N+1
X(i) = h * i;
Y(i) = h * i;
end
Xmid = X(1:N)+h/2;
Ymid = Y(1:N)+h/2;
% rectangles
xa = .02;
xb = .04;
xc = .075;
xd = .085;
ya = .03;
yb = .08;
Tc = 200;
Tf = -100;
rcstr = sprintf('{R, %f, {(%f, %f), (%f, %f)}}', Tc, xa, yb, xb, ya);
rfstr = sprintf('{R, %f, {(%f, %f), (%f, %f)}}', Tf, xc, yb, xd, ya);
%rc = [xa yb; xb ya];
%rf = [xc yb; xd ya];
%rboth = [xa yb; xd ya];
% triggered point
p = [0;0];
p(1) = (xb + xc) / 2;
p(2) = (ya + yb) / 2;
% transform point coordinates into indexes
p(1) = floor(p(1) / h) + 1;
p(2) = floor(p(2) / h) + 1;
%% Symulate and deduce stream %%
%%%%%%%%%%%%%%%%%%%%%%%%%
D = {};
%Jx = {};
%Jy = {};
Jnorm = {};
while xb < xc
% re-evaluate rectangles
rcstr = sprintf('{R, %f, {(%f, %f), (%f, %f)}}', Tc, xa, yb, xb, ya);
rfstr = sprintf('{R, %f, {(%f, %f), (%f, %f)}}', Tf, xc, yb, xd, ya);
% symulate
cmd = sprintf('%s%s %s dt=%.15g L=%f N=%d mode=%d \"chaud=%s\" \"froid=%s\" output=%s', repertoire, executable, input, dt, L, N, mode, rcstr, rfstr, output);
disp(cmd)
system(cmd);
% compute stream
T = load([output '_T.out']);
[jx, jy] = stream(T, p(1), p(2), h, kappa);
D = [D, xc - xb];
%Jx = [Jx, jx];
%Jy = [Jy, jy];
Jnorm = [Jnorm, sqrt(jx^2 + jy^2)];
% modify rectangles
xc = xc - h/2;
xd = xd - h/2;
xb = xb + h/2;
xa = xa + h/2;
end
D = cell2mat(D);
Jnorm = cell2mat(Jnorm);
%% Write data on output files %%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
out = fopen('outs/approaching.dat', 'w');
fprintf(out, 'D J\n');
for i = 1:length(D)
fprintf(out, '%.15g %.15g\n', D(i), Jnorm(i));
end
fclose(out);
%% Figures %%
%%%%%%%%%%%%%
% Stream :
figure
plot(D, Jnorm)
xlabel('Distance [s]')
ylabel('|J| [W/m]')
%legend('Jx', 'Jy', '|J|')
grid on
%% Utils %%
%%%%%%%%%%%%
%(i+1, j)%(i+1, j+1)
% %
% %
%(i, j)%%%(i, j+1)
function [jx, jy] = stream(T, i, j, h, kappa)
jx = -kappa * (T(i+1, j+1) + T(i, j+1) - T(i+1, j) - T(i, j)) /(2 * h);
jy = -kappa * (T(i+1, j) + T(i+1, j+1) - T(i, j) - T(i, j+1)) /(2 * h);
end
c4science · Help