From 8c0d59c0b43bfe860594cda73ec9189fd20cfbc8 Mon Sep 17 00:00:00 2001
From: Lars Hubatsch <hubatsch@pks.mpg.de>
Date: Mon, 28 Sep 2020 14:59:54 +0200
Subject: [PATCH] Correcting weak form spherical diffusion equation.
---
FloryHugg_DiffUnbleached.ipynb | 9 +++++----
1 file changed, 5 insertions(+), 4 deletions(-)
diff --git a/FloryHugg_DiffUnbleached.ipynb b/FloryHugg_DiffUnbleached.ipynb
index 3d430d7..e3d4508 100644
--- a/FloryHugg_DiffUnbleached.ipynb
+++ b/FloryHugg_DiffUnbleached.ipynb
@@ -226,12 +226,13 @@
" return on_boundary\n",
"bc = df.DirichletBC(F, u_D, boundary)\n",
"# Weak form spherical symmetry\n",
- "form = (df.inner((c-c0)/dt, q*X[0]*X[0]) +\n",
- " df.inner(df.grad(c), df.grad(X[0]*X[0]*q))-\n",
+ "form = ((c-c0)/dt*q*X[0]*X[0] +\n",
+ " X[0]*X[0]*df.inner(df.grad(c), df.grad(q))-\n",
" c.dx(0)*2*X[0]*q) * df.dx\n",
"# Weak form 1D\n",
- "# form = (df.inner((c-c0)/dt, q) +\n",
- "# df.inner(df.grad(c), df.grad(q))) * df.dx\n",
+ "# form = ((c-c0)/dt* q + df.inner(df.grad(c), df.grad(q))) * df.dx\n",
+ "# Weak form 1D with .dx(0) notation for derivative in 1st direction.\n",
+ "# form = ((c-c0)/dt*q + c.dx(0)*q.dx(0)) * df.dx\n",
"t = 0\n",
"\n",
"# Solve in time\n",
--
GitLab