Experiment with moving boundary entropy production. Not finished.

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