From 3baadc1d26d4b5eb2780d6c6813e420e4ed80e15 Mon Sep 17 00:00:00 2001
From: Lars Hubatsch <hubatsch@pks.mpg.de>
Date: Thu, 4 Feb 2021 20:10:10 +0100
Subject: [PATCH] Experiment with moving boundary entropy production. Not
finished.
---
prob_laplace.m | 13 +++++++++++--
1 file changed, 11 insertions(+), 2 deletions(-)
diff --git a/prob_laplace.m b/prob_laplace.m
index fc7e071..e4d074f 100644
--- a/prob_laplace.m
+++ b/prob_laplace.m
@@ -295,7 +295,7 @@ run_jump_lengths(-5, -0.00756985349, 7.7/3, 1, 4, '773_30_10', 10^-6, 4, 0.5);
%% %%%%%%%%%%%%%%%%%%% MOVING BOUNDARY TIME COURSE %%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% t = linspace(0, 10, 9000); % Entropy
-t = linspace(0, 10, 21); % Jump size distro
+t = linspace(0, 1, 1001); % Jump size distro
params_mov_bound = {-10, b(7/3, 10^-6), 0.5, e(7/3),...
0, 1, 20, 7, 0.5, 'Client', 0};
T_mov = Ternary_model(0, 'phi_tot', params_mov_bound, t, 0);
@@ -316,7 +316,8 @@ for i = 1:length(T_mov.t)
gra_a = Ternary_model.gradient_analytical(x_interp, T_mov.a, T_mov.b,...
T_mov.e, T_mov.v*T_mov.t(i));
dudx = diff(u)./diff(T_mov.x);
- f = -g0*(dudx+chi_phi.*u_interp.*gra_a);
+% g0*(dudx+chi_phi*u*gra_a)
+ f = g0*(dudx+chi_phi.*u_interp.*gra_a);
s_dot(i) = sum(diff(T_mov.x).*f.^2./(g0.*u_interp));
% figure(1); hold on;
% plot(0:10, zeros(1, 11), 'LineWidth', 2);
@@ -328,6 +329,14 @@ for i = 1:length(T_mov.t)
end
% csvwrite([s, 'Mov_Bou_Flux.csv'], [x_interp; f])
% csvwrite([s, 'Mov_Bou_Entr.csv'], [T_mov.t; s_dot])
+%%
+col = 'b';
+figure(2); hold on;
+% plot(u, col)
+% plot(u_interp, col)
+% plot(gra_a, col)
+plot(dudx, col)
+% plot(f, col)
%% MOVING BOUNDARY JUMP LENGTH DISTRIBUTION
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% load entropy_mov_bound_in_to_out.mat
--
GitLab