{ "cells": [ { "cell_type": "markdown", "id": "8382e30e-2fac-41f1-ae50-baf0fc9f4f22", "metadata": {}, "source": [ "# Exploring the change of delta_x (spatial resolution) in diffusion accuracy.\n", "#### From the same initial setup, diffusion is carried out over a fixed time span,\n", "#### at different spatial resolutions - and then the respective results are compared" ] }, { "cell_type": "markdown", "id": "55e63560-e412-4463-9195-c8746a17add5", "metadata": {}, "source": [ "### TAGS : \"diffusion 1D\", \"under-the-hood\"" ] }, { "cell_type": "code", "execution_count": 1, "id": "7b9ba7bb-02d5-47e3-a6f4-01dad4a07de5", "metadata": {}, "outputs": [], "source": [ "LAST_REVISED = \"May 3, 2025\"\n", "LIFE123_VERSION = \"1.0.0rc3\" # Library version this experiment is based on" ] }, { "cell_type": "code", "execution_count": 2, "id": "b3ae84bc-72b9-4043-b5b9-70532eafdb44", "metadata": {}, "outputs": [], "source": [ "#import set_path # Using MyBinder? Uncomment this before running the next cell!" ] }, { "cell_type": "code", "execution_count": 3, "id": "6d524cb0-a0bf-4533-9803-a0e4985203ec", "metadata": {}, "outputs": [], "source": [ "#import sys\n", "#sys.path.append(\"C:/some_path/my_env_or_install\") # CHANGE to the folder containing your venv or libraries installation!\n", "# NOTE: If any of the imports below can't find a module, uncomment the lines above, or try: import set_path \n", "\n", "from life123 import BioSim1D, ChemData, Numerical, check_version" ] }, { "cell_type": "code", "execution_count": 4, "id": "73a03932-99be-49ac-8948-ca3d6f0f2e76", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "OK\n" ] } ], "source": [ "check_version(LIFE123_VERSION)" ] }, { "cell_type": "code", "execution_count": null, "id": "4c1e251e-3db1-4abf-8c6f-1828e5de2f9e", "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "markdown", "id": "28c89c3d-a7f7-41d9-a99b-b0185620b5c1", "metadata": {}, "source": [ "## Prepare the initial system" ] }, { "cell_type": "code", "execution_count": 5, "id": "01b3a969-5122-4c25-900b-ad6fba315553", "metadata": { "tags": [] }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "SYSTEM STATE at Time t = 0:\n", "60 bins and 1 chemical species:\n" ] }, { "data": { "text/html": [ "
| \n", " | Species | \n", "Diff rate | \n", "Bin 0 | \n", "Bin 1 | \n", "Bin 2 | \n", "Bin 3 | \n", "Bin 4 | \n", "Bin 5 | \n", "Bin 6 | \n", "Bin 7 | \n", "... | \n", "Bin 50 | \n", "Bin 51 | \n", "Bin 52 | \n", "Bin 53 | \n", "Bin 54 | \n", "Bin 55 | \n", "Bin 56 | \n", "Bin 57 | \n", "Bin 58 | \n", "Bin 59 | \n", "
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 0 | \n", "A | \n", "0.1 | \n", "10.0 | \n", "13.0 | \n", "17.0 | \n", "21.0 | \n", "25.0 | \n", "28.0 | \n", "30.0 | \n", "38.0 | \n", "... | \n", "37.0 | \n", "31.0 | \n", "25.0 | \n", "22.0 | \n", "20.0 | \n", "18.0 | \n", "15.0 | \n", "11.0 | \n", "9.0 | \n", "8.0 | \n", "
1 rows × 62 columns
\n", "| \n", " | Species | \n", "Diff rate | \n", "Bin 0 | \n", "Bin 1 | \n", "Bin 2 | \n", "Bin 3 | \n", "Bin 4 | \n", "Bin 5 | \n", "Bin 6 | \n", "Bin 7 | \n", "... | \n", "Bin 109 | \n", "Bin 110 | \n", "Bin 111 | \n", "Bin 112 | \n", "Bin 113 | \n", "Bin 114 | \n", "Bin 115 | \n", "Bin 116 | \n", "Bin 117 | \n", "Bin 118 | \n", "
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 0 | \n", "A | \n", "0.1 | \n", "10.0 | \n", "11.5 | \n", "13.0 | \n", "15.0 | \n", "17.0 | \n", "19.0 | \n", "21.0 | \n", "23.0 | \n", "... | \n", "19.0 | \n", "18.0 | \n", "16.5 | \n", "15.0 | \n", "13.0 | \n", "11.0 | \n", "10.0 | \n", "9.0 | \n", "8.5 | \n", "8.0 | \n", "
1 rows × 121 columns
\n", "| \n", " | Species | \n", "Diff rate | \n", "Bin 0 | \n", "Bin 1 | \n", "Bin 2 | \n", "Bin 3 | \n", "Bin 4 | \n", "Bin 5 | \n", "Bin 6 | \n", "Bin 7 | \n", "... | \n", "Bin 109 | \n", "Bin 110 | \n", "Bin 111 | \n", "Bin 112 | \n", "Bin 113 | \n", "Bin 114 | \n", "Bin 115 | \n", "Bin 116 | \n", "Bin 117 | \n", "Bin 118 | \n", "
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 0 | \n", "A | \n", "0.1 | \n", "10.0 | \n", "11.5 | \n", "13.0 | \n", "15.0 | \n", "17.0 | \n", "19.0 | \n", "21.0 | \n", "23.0 | \n", "... | \n", "19.0 | \n", "18.0 | \n", "16.5 | \n", "15.0 | \n", "13.0 | \n", "11.0 | \n", "10.0 | \n", "9.0 | \n", "8.5 | \n", "8.0 | \n", "
1 rows × 121 columns
\n", "| \n", " | Species | \n", "Diff rate | \n", "Bin 0 | \n", "Bin 1 | \n", "Bin 2 | \n", "Bin 3 | \n", "Bin 4 | \n", "Bin 5 | \n", "Bin 6 | \n", "Bin 7 | \n", "... | \n", "Bin 228 | \n", "Bin 229 | \n", "Bin 230 | \n", "Bin 231 | \n", "Bin 232 | \n", "Bin 233 | \n", "Bin 234 | \n", "Bin 235 | \n", "Bin 236 | \n", "Bin 237 | \n", "
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 0 | \n", "A | \n", "0.1 | \n", "10.0 | \n", "10.0 | \n", "11.5 | \n", "11.5 | \n", "13.0 | \n", "13.0 | \n", "15.0 | \n", "15.0 | \n", "... | \n", "11.0 | \n", "11.0 | \n", "10.0 | \n", "10.0 | \n", "9.0 | \n", "9.0 | \n", "8.5 | \n", "8.5 | \n", "8.0 | \n", "8.0 | \n", "
1 rows × 240 columns
\n", "| \n", " | Species | \n", "Diff rate | \n", "Bin 0 | \n", "Bin 1 | \n", "Bin 2 | \n", "Bin 3 | \n", "Bin 4 | \n", "Bin 5 | \n", "Bin 6 | \n", "Bin 7 | \n", "... | \n", "Bin 109 | \n", "Bin 110 | \n", "Bin 111 | \n", "Bin 112 | \n", "Bin 113 | \n", "Bin 114 | \n", "Bin 115 | \n", "Bin 116 | \n", "Bin 117 | \n", "Bin 118 | \n", "
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 0 | \n", "A | \n", "0.1 | \n", "10.0 | \n", "11.5 | \n", "13.0 | \n", "15.0 | \n", "17.0 | \n", "19.0 | \n", "21.0 | \n", "23.0 | \n", "... | \n", "19.0 | \n", "18.0 | \n", "16.5 | \n", "15.0 | \n", "13.0 | \n", "11.0 | \n", "10.0 | \n", "9.0 | \n", "8.5 | \n", "8.0 | \n", "
1 rows × 121 columns
\n", "| \n", " | Species | \n", "Diff rate | \n", "Bin 0 | \n", "Bin 1 | \n", "Bin 2 | \n", "Bin 3 | \n", "Bin 4 | \n", "Bin 5 | \n", "Bin 6 | \n", "Bin 7 | \n", "... | \n", "Bin 466 | \n", "Bin 467 | \n", "Bin 468 | \n", "Bin 469 | \n", "Bin 470 | \n", "Bin 471 | \n", "Bin 472 | \n", "Bin 473 | \n", "Bin 474 | \n", "Bin 475 | \n", "
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 0 | \n", "A | \n", "0.1 | \n", "10.0 | \n", "10.0 | \n", "10.0 | \n", "10.0 | \n", "11.5 | \n", "11.5 | \n", "11.5 | \n", "11.5 | \n", "... | \n", "9.0 | \n", "9.0 | \n", "8.5 | \n", "8.5 | \n", "8.5 | \n", "8.5 | \n", "8.0 | \n", "8.0 | \n", "8.0 | \n", "8.0 | \n", "
1 rows × 478 columns
\n", "