From 880c691b5ae6a2229fa6435154d68d26a495147d Mon Sep 17 00:00:00 2001
From: Lars Hubatsch <hubatsch@pks.mpg.de>
Date: Fri, 25 Sep 2020 16:18:30 +0200
Subject: [PATCH] Square root dependence between diffusion and partitioning.
---
FloryHugg_DiffUnbleached.ipynb | 20 ++++++++++++++++----
1 file changed, 16 insertions(+), 4 deletions(-)
diff --git a/FloryHugg_DiffUnbleached.ipynb b/FloryHugg_DiffUnbleached.ipynb
index 04e3b89..3d430d7 100644
--- a/FloryHugg_DiffUnbleached.ipynb
+++ b/FloryHugg_DiffUnbleached.ipynb
@@ -92,13 +92,18 @@
"c0_2 = create_func(F, 'x[0]<0.1 ? 0 :'+p_tot(0.9, 0.45), 1)\n",
"Ga0_2 = create_func(F,p_tot(1, 0.08), 1)\n",
"\n",
+ "c_tot_3 = create_func(F, p_tot(0.9, 0.3), 1)\n",
+ "c0_3 = create_func(F, 'x[0]<0.1 ? 0 :'+p_tot(0.9, 0.3), 1)\n",
+ "Ga0_3 = create_func(F,p_tot(1, 0.145), 1)\n",
+ "\n",
"c_tot_9 = create_func(F, p_tot(0.9, 0.1), 1)\n",
"c0_9 = create_func(F, 'x[0]<0.1 ? 0 :'+p_tot(0.9, 0.1), 1)\n",
"Ga0_9 = create_func(F,p_tot(1, 1), 1)\n",
"\n",
- "c0_1 = calc_sim(c0_1, c_tot_1, Ga0_1, 0.001, 50, 0)\n",
- "c0_2 = calc_sim(c0_2, c_tot_2, Ga0_2, 0.001, 50, 0)\n",
- "c0_9 = calc_sim(c0_9, c_tot_9, Ga0_9, 0.001, 50, 0)"
+ "c0_1 = calc_sim(c0_1, c_tot_1, Ga0_1, 0.001, 10, 0)\n",
+ "c0_2 = calc_sim(c0_2, c_tot_2, Ga0_2, 0.001, 10, 0)\n",
+ "c0_9 = calc_sim(c0_9, c_tot_9, Ga0_9, 0.001, 10, 0)\n",
+ "c0_3 = calc_sim(c0_3, c_tot_3, Ga0_3, 0.001, 10, 0)"
]
},
{
@@ -110,8 +115,15 @@
"x = np.linspace(0, 1, 10000)\n",
"plt.plot(x, eval_func(c0_1, x))\n",
"plt.plot(x, eval_func(c0_2, x))\n",
+ "plt.plot(x, eval_func(c0_3, x))\n",
"plt.plot(x, eval_func(c0_9, x))\n",
- "plt.xlim(0.0, 0.125)"
+ "plt.xlim(0.0, 0.125)\n",
+ "plt.show()\n",
+ "# Diffusion versus partitioning\n",
+ "plt.plot([Ga0_1(1)*(1-c_tot_1(1)), Ga0_2(1)*(1-c_tot_2(1)), Ga0_3(1)*(1-c_tot_3(1)), Ga0_9(1)*(1-c_tot_9(1))],\n",
+ " [c_tot_1(0)/c_tot_1(1), c_tot_2(0)/c_tot_2(1), c_tot_3(0)/c_tot_3(1), c_tot_9(0)/c_tot_9(1)])\n",
+ "plt.ylim(0, 9); plt.xlim(0, 1)\n",
+ "plt.plot(np.linspace(0, 1, 100), 9.5*np.sqrt(np.linspace(0, 1, 100)))"
]
},
{
--
GitLab