From b2f016602cb312f8f43a555a10ee869729b83445 Mon Sep 17 00:00:00 2001
From: Lars Hubatsch <hubatsch@pks.mpg.de>
Date: Mon, 9 Nov 2020 22:19:39 +0100
Subject: [PATCH] Testing several D_i/D_o ratios, generally seems to work.
---
prob_laplace.m | 20 +++++++++++++++++++-
1 file changed, 19 insertions(+), 1 deletion(-)
diff --git a/prob_laplace.m b/prob_laplace.m
index d38125c..58e2774 100644
--- a/prob_laplace.m
+++ b/prob_laplace.m
@@ -186,11 +186,29 @@ for i = 1:100%length(T_mov.t)
shg; pause();
end
-
+%% Different ratios of Dm/Dp
+chi = 7/3; nu = 10^-6; u0 = 0.5; a = -5;
+b1 = b(chi, nu); P = (u0+e(chi))/(u0-e(chi));
+params = {a, b1, u0, e(chi), -0.15432893, 0, 10, 7, 0, 'Const_mob', 0};
+%% Stefano's way of having mobility with square. works, 0.5% diff
+chi = 7/3; nu = 10^-6; u0 = 0.5; a = -5;
+b1 = b(chi, nu); P = (u0+e(chi))/(u0-e(chi));
+params = {a, b1, u0, e(chi), -0.12314694, 0, 10, 7, 0, 'Const_mob', 0};
+%%
+D_i = 0.05672749; D_o = 0.85091231; chi = 7/3; nu = 10^-6; u0 = 0.5; a = -5;
+b1 = b(chi, nu); P = (u0+e(chi))/(u0-e(chi));
+[~, e_g, u_g] = calc_tanh_params(P, D_i, D_o, a, b1, u0);
+%%
+D_i = 0.03; D_o = 0.83; chi = 7/3; nu = 10^-6; u0 = 0.5; a = -5;
+b1 = b(chi, nu); P = (u0+e(chi))/(u0-e(chi));
+[~, e_g, u_g] = calc_tanh_params(P, D_i, D_o, a, b1, u0);
+%%
+params = {a, b1, u0, e(chi), e_g, u_g, 10, 7, 0, 'Constituent', 0};
%% %%%%%%%%%%%%%%%%% STEADY STATE JUMP LENGTH DISTRIBUTION %%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
params = {-5, b(7.7/3, 10^-6), 0.5, e(7.7/3), 0, 1, 10, 7, 0,...
'Constituent', 0};
+%%
t = [0, 0.05, 0.1, 1];
direc = 1;
x0 = sort(5-direc*(0:0.001:4.01));
--
GitLab