Add files via upload

Author: dindagustiayu

MB-interactive.ipynb

@@ -0,0 +1,103 @@
+{
+ "cells": [
+  {
+   "cell_type": "code",
+   "execution_count": 1,
+   "id": "c48badf8-e0a6-4d31-a3be-8391267f0c51",
+   "metadata": {},
+   "outputs": [
+    {
+     "data": {
+      "application/vnd.jupyter.widget-view+json": {
+       "model_id": "de4b4fa1f17446bbb1b4a5e4dcb1f408",
+       "version_major": 2,
+       "version_minor": 0
+      },
+      "text/plain": [
+       "interactive(children=(IntSlider(value=300, description='T', max=1000, min=100, step=50), Output()), _dom_class…"
+      ]
+     },
+     "metadata": {},
+     "output_type": "display_data"
+    },
+    {
+     "data": {
+      "text/plain": [
+       "<function __main__.plot_MB(T=300)>"
+      ]
+     },
+     "execution_count": 1,
+     "metadata": {},
+     "output_type": "execute_result"
+    }
+   ],
+   "source": [
+    "import numpy as np\n",
+    "import matplotlib.pyplot as plt\n",
+    "from ipywidgets import interact\n",
+    "\n",
+    "# Boltzman constant, J.K-1, atomis mass unit (kg)\n",
+    "kB, u = 1.381e-23, 1.661e-27\n",
+    "\n",
+    "# Gas particle masses, in kg\n",
+    "m_He, m_Ar = 4 * u, 40 * u\n",
+    "\n",
+    "# Temperature in K\n",
+    "T = 300\n",
+    "\n",
+    "# The value of the Maxwell-Boltzmann distribution\n",
+    "def fMB(v, T, m): # mass in kg, at temperature in K, at the speed v (in m.s-1)\n",
+    "    fac = m /2 /kB/T\n",
+    "    return (fac / np.pi)**1/5 * 4 * np.pi * v**2 * np.exp(-fac * v **2)\n",
+    "\n",
+    "def plot_MB(T=300):\n",
+    "    v= np.linspace(0, 3000, 3000)\n",
+    "    fv_He = fMB(v, T, m_He)\n",
+    "    fv_Ar = fMB(v, T, m_Ar)\n",
+    "\n",
+    "# Plot He and Ar distribution\n",
+    "    plt.plot(figsize=(8,6))\n",
+    "    plt.plot(v, fv_He, color='orange', label=f'He, T{T}K')\n",
+    "    plt.plot(v, fv_Ar, color='blue', ls='--', label=f'Ar, T{T}K')\n",
+    "    \n",
+    "    plt.xlabel('Speed v /m.s-1', fontsize=14)\n",
+    "    plt.ylabel('Density f(v) /s.m-1', fontsize=14)\n",
+    "    plt.legend()\n",
+    "    plt.legend(framealpha=0.7, fontsize=14, loc='upper right', bbox_to_anchor=(0.98, 1))\n",
+    "    plt.grid(True, color='gray', zorder=0, linestyle='--')\n",
+    "    plt.show()\n",
+    "\n",
+    "interact(plot_MB, T=(100, 1000, 50))"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "893b6a24-775b-48ec-ac44-66d188be717f",
+   "metadata": {},
+   "outputs": [],
+   "source": []
+  }
+ ],
+ "metadata": {
+  "kernelspec": {
+   "display_name": "myenv",
+   "language": "python",
+   "name": "myenv"
+  },
+  "language_info": {
+   "codemirror_mode": {
+    "name": "ipython",
+    "version": 3
+   },
+   "file_extension": ".py",
+   "mimetype": "text/x-python",
+   "name": "python",
+   "nbconvert_exporter": "python",
+   "pygments_lexer": "ipython3",
+   "version": "3.13.9"
+  }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 5
+}