Yuuuuuuuuuuuuup

Author: JulesBlm

JulesD3D/new_sand_records_times.txt renamed to 7525_15.bcc

@@ -2293,6 +2293,10 @@
    "cell_type": "markdown",
    "metadata": {},
    "source": [
+    "# No longer developed because I don't know how to calculate it with the parameters in Delft3D\n",
+    "* How do I calculate 'channel depth'?\n",
+    "* Think it would still be interesting tho\n",
+    "\n",
     "# Densimetric Froude number\n",
     "\n",
     "$Fr = \\frac{u}{\\sqrt{g'\\cdot L}}=\\frac{\\text{Layer averaged flow velocity}}{\\text{Wave velocity in shallow water}}=\\frac{\\text{Inertial forces}}{\\text{Gravity forces}}$\n",
@@ -2308,6 +2312,7 @@
    "cell_type": "markdown",
    "metadata": {},
    "source": [
+    "\n",
     "The densimetric Froude number represents (1) the ratio of the current speed to the speed of a long wave on the dense flow interface, and (2) the ratio between the destabilizing effects of velocity shearing (inertial forces) over the stabilizing effects of vertical density stratification (gravity/ buoyancy forces). The densimetric Froude number Frd is defined as"
    ]
   },
@@ -3822,7 +3827,7 @@
    "name": "python",
    "nbconvert_exporter": "python",
    "pygments_lexer": "ipython3",
-   "version": "3.7.5"
+   "version": "3.7.6"
   }
  },
  "nbformat": 4,

FroudeNumber.ipynb

@@ -72,8 +72,8 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "sand_weight_vol = 39.74 # [kg/m3] # 19.875 # 25%\n",
-    "silt_weight_vol = 39.74 # [kg/m3] # 59.625 # 75%"
+    "sand_weight_vol = 19.875 # [kg/m3] #  # 25% 39.74\n",
+    "silt_weight_vol = 59.625 # [kg/m3] #  # 75% 39.74"
    ]
   },
   {
@@ -103,10 +103,9 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "timestep = 0.3 # ureg.minute [min]\n",
-    "start_discharge_time = 9 # ureg.minute [min]\n",
-    "end_discharge_time = 48 # ureg.minute [min]\n",
-    "# end_discharge_time = 123 # ureg.minute [min]"
+    "timestep = 0.3 # [min]\n",
+    "start_discharge_time = 9 # [min]\n",
+    "end_discharge_time = 48 # [min]"
    ]
   },
   {
@@ -143,7 +142,10 @@
     {
      "data": {
       "text/markdown": [
-       "Is multiple of timestep? It is! if True"
+       "\n",
+       "# Sanity check\n",
+       "Weird check for floating point errors\n",
+       "## Is multiple of timestep? True!\n"
       ],
       "text/plain": [
        "<IPython.core.display.Markdown object>"
@@ -156,7 +158,11 @@
    ],
    "source": [
     "from decimal import Decimal\n",
-    "md(f\"Is multiple of timestep? It is! if {Decimal(discharge_duration) % Decimal(timestep) < 10e-9}\") # weird check due to floatin point errors"
+    "md(f'''\n",
+    "# Sanity check\n",
+    "Weird check for floating point errors\n",
+    "## Is multiple of timestep? {Decimal(discharge_duration) % Decimal(timestep) < 10e-9}!\n",
+    "''') "
    ]
   },
   {
@@ -181,7 +187,7 @@
        " '9.0000000e+000',\n",
        " '4.8000000e+001',\n",
        " '4.8300000e+001',\n",
-       " '1.0200000e+003']"
+       " '1.2300000e+003']"
       ]
      },
      "execution_count": 11,
@@ -190,12 +196,13 @@
     }
    ],
    "source": [
+    "# TODO: hardcoded simultion start + edn time string\n",
     "times = [\"0.0000000e+000\", \n",
     "         f\"{formatSci(pre_start_discharge_time)}\", \n",
     "         f\"{formatSci(start_discharge_time)}\",\\\n",
     "         f\"{formatSci(end_discharge_time)}\", \n",
     "         f\"{formatSci(post_end_discharge_time)}\", \n",
-    "         \"1.0200000e+003\"]\n",
+    "         \"1.2300000e+003\"]\n",
     "times"
    ]
   },
@@ -273,7 +280,7 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "# TODO: should loop over sediments and then write each line\n",
+    "# TODO: should loop over arrray of sediments and then write each line but for only two like now it works fine\n",
     "sand_conc_records_discharge_layer = '  '.join(formatted_sand_conc for i in range(nr_discharge_layers))\n",
     "silt_conc_records_discharge_layer = '  '.join(formatted_silt_conc for i in range(nr_discharge_layers))\n",
     "\n",
@@ -596,7 +603,7 @@
      "data": {
       "text/markdown": [
        "\n",
-       "* **Total sediment volume flowing in = 3.1582052830188676 $\\times10^5 m^3$**\n"
+       "* **Total sediment volume flowing in = 3.159 $\\times10^5 m^3$**\n"
       ],
       "text/plain": [
        "<IPython.core.display.Markdown object>"
@@ -646,15 +653,15 @@
        "* Discharge duration = 39 min\n",
        "\n",
        "### Sediments\n",
-       "* Silt weight per volume = 39.74 $kg/m^3$\n",
-       "* Sand weight per volume = 39.74 $kg/m^3$\n",
+       "* Silt weight per volume = 59.625 $kg/m^3$\n",
+       "* Sand weight per volume = 19.875 $kg/m^3$\n",
        "\n",
-       "* Total silt volume 1.58 $\\times10^5 m^3$\n",
-       "* Total sand volume 1.58 $\\times10^5 m^3$\n",
-       "* **Total sediment volume = 3.158 $\\times10^5 m^3$**\n",
+       "* Total silt volume 0.79 $\\times10^5 m^3$\n",
+       "* Total sand volume 2.37 $\\times10^5 m^3$\n",
+       "* **Total sediment volume = 3.159 $\\times10^5 m^3$**\n",
        "\n",
-       "* Silt Concentration = 1.5%\n",
-       "* Sand Concentration =  1.5%\n",
+       "* Silt Concentration = 2.25%\n",
+       "* Sand Concentration =  0.75%\n",
        "* Total sediment concentration = 3.0%\n",
        "\n"
       ],