{"nbformat":4,"nbformat_minor":0,"metadata":{"colab":{"provenance":[],"authorship_tag":"ABX9TyOmzYQvATWJjCgD+dIZlx+z"},"kernelspec":{"name":"python3","display_name":"Python 3"},"language_info":{"name":"python"}},"cells":[{"cell_type":"markdown","source":["# Thermal Calibration Experiment\n","**Calculations and notes**\n","\n","Justin Tucker\n","\n","Last updated: April 1st 2024\n"],"metadata":{"id":"fS825N5WfnQx"}},{"cell_type":"code","source":["import numpy as np"],"metadata":{"id":"xczM7pl7fnmR","executionInfo":{"status":"ok","timestamp":1719962104456,"user_tz":420,"elapsed":317,"user":{"displayName":"Justin Tucker","userId":"10374969465321862369"}}},"execution_count":2,"outputs":[]},{"cell_type":"markdown","source":["### Test Tank volume"],"metadata":{"id":"xRsQVmImU3N0"}},{"cell_type":"code","source":["tank_volume_cubic_meters = 10 * 10 * 15\n","gallons_per_cubic_meter = 264.172\n","tank_volume_gallons = tank_volume_cubic_meters * gallons_per_cubic_meter\n","\n","print(\"Tank Volume: %d gallons\" % tank_volume_gallons)"],"metadata":{"colab":{"base_uri":"https://localhost:8080/"},"id":"G7tvW5ZbVytn","executionInfo":{"status":"ok","timestamp":1719962104764,"user_tz":420,"elapsed":3,"user":{"displayName":"Justin Tucker","userId":"10374969465321862369"}},"outputId":"5260af5a-a9cb-41a6-d374-55d49ed6cca4"},"execution_count":3,"outputs":[{"output_type":"stream","name":"stdout","text":["Tank Volume: 396258 gallons\n"]}]},{"cell_type":"markdown","source":["Concept A - heated plume"],"metadata":{"id":"zZpiIv6fpzX8"}},{"cell_type":"markdown","source":["### Power calcs for continuous hose delivery"],"metadata":{"id":"SH_VihrqVzBF"}},{"cell_type":"code","source":["cp = 4.186  # J/g°C water\n","rho_sea_water = 1024 * 10**3  # g/m^3 (1020-1029kg/m^3)\n","gpm_to_cubic_mps = 1 / 15850.3  # Gallons per min to m^3/s\n","\n","garden_hose_gpm = 15  # flow rate of hose 9 - 17gpm\n","garden_hose_cubic_mps = garden_hose_gpm * gpm_to_cubic_mps  # flow rate in m^3/s\n","\n","fire_hose_gpm = 150  # flow rate of fire hose 120-180gpm (from Gainesville)\n","fire_hose_cubic_mps = fire_hose_gpm * gpm_to_cubic_mps  # flow rate in m^3/s\n","\n","starting_temp_c = 15  # initial water temp in C\n","target_temp1_c = 50  # temp in C\n","delta_t1 = target_temp1_c - starting_temp_c\n","\n","power_garden_hose_temp1 = cp * rho_sea_water * garden_hose_cubic_mps * delta_t1  # Power in W\n","power_fire_hose_temp1 = cp * rho_sea_water * fire_hose_cubic_mps * delta_t1  # Power in W\n","\n","target_temp2_c = 90  # temp in C\n","delta_t2 = target_temp2_c - starting_temp_c\n","\n","power_garden_hose_temp2 = cp * rho_sea_water * garden_hose_cubic_mps * delta_t2  # Power in W\n","power_fire_hose_temp2 = cp * rho_sea_water * fire_hose_cubic_mps * delta_t2  # Power in W"],"metadata":{"id":"fUJXh7QlU4a7","executionInfo":{"status":"ok","timestamp":1719966481577,"user_tz":420,"elapsed":3,"user":{"displayName":"Justin Tucker","userId":"10374969465321862369"}}},"execution_count":16,"outputs":[]},{"cell_type":"code","source":["print(\"At a garden hose flow rate (%d gpm), the continuous power required to heat \\\n","water from %dC to %dC is %d W\" %(garden_hose_gpm,starting_temp_c,target_temp1_c,power_garden_hose_temp1))\n","\n","print(\"At a garden hose flow rate (%d gpm), the continuous power required to heat \\\n","water from %dC to %dC is %d W\" %(garden_hose_gpm,starting_temp_c,target_temp2_c,power_garden_hose_temp2))\n","\n","print(\"At a fire hose flow rate (%d gpm), the continuous power required to heat \\\n","water from %dC to %dC is %d W\" %(fire_hose_gpm,starting_temp_c,target_temp1_c,power_fire_hose_temp1))\n","\n","print(\"At a fire hose flow rate (%d gpm), the continuous power required to heat \\\n","water from %dC to %dC is %d W\" %(fire_hose_gpm,starting_temp_c,target_temp2_c,power_fire_hose_temp2))"],"metadata":{"colab":{"base_uri":"https://localhost:8080/"},"id":"VaChBTOhij4d","executionInfo":{"status":"ok","timestamp":1719962104764,"user_tz":420,"elapsed":2,"user":{"displayName":"Justin Tucker","userId":"10374969465321862369"}},"outputId":"f566429a-959c-41c0-b976-c2451bda91ea"},"execution_count":5,"outputs":[{"output_type":"stream","name":"stdout","text":["At a garden hose flow rate (15 gpm), the continuous power required to heat water from 15C to 38C is 93299 W\n","At a garden hose flow rate (15 gpm), the continuous power required to heat water from 15C to 90C is 304238 W\n","At a fire hose flow rate (150 gpm), the continuous power required to heat water from 15C to 38C is 932998 W\n","At a fire hose flow rate (150 gpm), the continuous power required to heat water from 15C to 90C is 3042385 W\n"]}]},{"cell_type":"markdown","source":["Power and time calcs for insulated storage tank"],"metadata":{"id":"Bi_2YS_ic_cd"}},{"cell_type":"code","source":["supplyPower=10000 # (W)\n","heatedTankVolumeGallons = 500 # (gallons)\n","heatedTankVolume = heatedTankVolumeGallons / gallons_per_cubic_meter # (m^3)\n","heatedTankMass = heatedTankVolume*rho_sea_water # (g)\n","heatedTankEnergy=cp*heatedTankMass*(delta_t1+0)\n","timeToHeatHr=heatedTankEnergy/supplyPower/(60*60)\n","print(\"Heating %d gallons, with %d W will take %d hrs assuming perfect efficency.\"\n"," %(heatedTankVolumeGallons, supplyPower,timeToHeatHr))\n"],"metadata":{"id":"sQtvC6JjdIvu","colab":{"base_uri":"https://localhost:8080/"},"executionInfo":{"status":"ok","timestamp":1719966746929,"user_tz":420,"elapsed":313,"user":{"displayName":"Justin Tucker","userId":"10374969465321862369"}},"outputId":"dc1c7d87-d1a6-4042-a51e-ff99a4c80908"},"execution_count":20,"outputs":[{"output_type":"stream","name":"stdout","text":["Heating 500 gallons, with 10000 W will take 7 hrs assuming perfect efficency.\n"]}]},{"cell_type":"markdown","source":["#### Voxel heat ratio\n","Assumes a perfect column of heated water extends from the hose of the same\n","diameter of the hose with no mixing.\n","\n","![gardenHoseDiamCyclinder.png](data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAmQAAAJLCAYAAABJ64pfAAAAAXNSR0IArs4c6QAAAARnQU1BAACxjwv8YQUAAAAJcEhZcwAADsMAAA7DAcdvqGQAACOJSURBVHhe7d1/kCRXYdjxvj1OlA4KGawzIoBdEbFFXJIhLqpMmUAAGYGIDdhEJDicDLJibFRyJITAMhDKNpVYqKKToBQhxC8jGYj9ByE/CidlJJcxlLBTAio4MT+klAUIEShACToZne8u/bpf78zOzu7OznbP6379+Whn+r03czqv7e396s1uz76TpQIAgGTW4hEAgEQEGQBAYoIMACAxQQYAkJggAwBITJABACQmyAAAEhNkAACJCTIAgMQEGQBAYoIMACAxQQYAkJggAwBITJABACQmyAAAEhNkAACJCTIAgMQEGQBAYoIMACAxQQYAkJggAwBITJABACQmyAAAEhNkAACJCTIAgMQEGQBAYoIMACAxQQYAkJggAwBITJABACQmyAAAEhNkAACJCTIAgMQEGQBAYoIMACAxQQYAkJggAwBITJABACQmyAAAEhNkAACJCTIAgMQEGQBAYoIMACAxQQYAkJggAwBITJABACQmyAAAEhNkAACJCTIAgMQEGQBAYoIMACAxQQYAkJggAwBITJABACQmyAAAEhNkAACJCTIAgMQEGQBAYoIMACAxQQYAkJggAwBITJABACQmyAAAEhNkAACJCTIAgMQEGQBAYoIMACCxfSdLcQxA5vbt2xdHO/PtAVZHkAEM3G4i687PfiGOiuKBB/5fHE1cfPFFxdW/+2+q8Ute/I+ro28T0D1BBtBTi4bWdGQF80JrUdNB1hBm0D1BBrBCbe5mdWFekDWaMAt864B2CTKAPVo2soJVhdaitguyaXbNoF2CDGALKV4yTG3RIGsIM2iHIANGJafdrC7sNsgaXs6EvRFkQBbGuJvVhWWDbJpdM9g9QQb0lt2s1WsjyBrCDBYnyICVs5vVX20GWUOYwc4EGdAKu1l56CLIGn7ODLYmyIBt2c0aly6DbJpdM9hIkMEI2c1iK6sKsoYwg5ogg4zYzWKvVh1kDS9nMnaCDHpu2d0skcUyUgXZNLtmjJEgg0TsZtFHfQiyhjBjTAQZtMhuFkPXpyBreDmTMRBksAM/AM+Y9DHIGmec8YTi6T/1lGrsWxe5EWSMlpcMYbM+B9k0L2eSG0FGVuxmwd4MJcgawoxcCDIGwW4WrMbQgqzh58wYOkFGMnazoH+GGmTT7JoxRIKM1tnNguHKIcgawowhEWQsxG4WjENOQdbwciZDIMhGzm4WMC3HIJtm14y+EmQZWnY3S2QBuQdZQ5jRN4JsQOxmAV0bS5A1hBl9IcgSs5sF9MnYgqzh58xITZB1xG4WMERjDbJpds1IQZDtgt0sIHeCbEKYsUqjD7JlIysQWkBuBNlmXs5kFbINMi8ZAuyeINueXTO6kl2QTYeYlw0BdkeQLUaY0ba1eMzKx/7otur2k089q7qJMQDa9B8++l+qW9gEaG6wF1kGWaMJs2c+42nVDQDa1ITZHZ/+nDBjT7IOssZsmIkzANp0331f3bRrBrsxiiBrNGFm1wyArggzljGqIJsmzADo0myYiTO2M9ogawgzALrUhJldM7Yz+iBrzIaZOAOgbcKMrQiyGU2Y2TUDoCuzYSbOEGTbEGYAdKkJM7tmCLIFzIaZOAOgbcJs3LJ866QQT107/wXPrY6f+OR/r47A7q3qP258nS7OWyf1h7dnGhdBtkfCDOZbJLZuu/2TxZYnoF2emSZP3/wHf+a5/zCO5vP1OyHI+qcJs0Cc5UuQtaQJs8DJnbGZF1+3/cmnFoiqk1s/ZZdnpsnT5/zBcmnu31Se/sLqeT/zrHo+Zaxfx4Ks3+ya5UuQdcCuGbmbDbDbbi/jK5j9sZcdzy79CLJ11aRee/7z/lGYrBvL17MgGwZhlh9B1iFhRm6aEFsPsFkZBdnG8cniBc9/Thhl//UsyIbFy5n5EGQrIMwYuibEPn7bJ6vo2repvKKMg6x+3sni/Ofn/fUsyIbLrtmwCbIV8nNmDFGIsT++7c8mETbyIKs+yvkLzz83y69jQTZ8wmyYXIdshUIoNrfwTW7253CA7X3846cXDz3ktAXbmb2eWbjRf85siQgz2L1zz/1WccopJ+IsPV+79FkTZnd8+nPCbAAEWWKzYeYET9+c9fd/YsfreI1NeLnyxnfdWv3vBvruvvu+umnXjP4RZD3RhJldM/roxptuKc597jOq25iFH+gPPwt6080fjCswLMKsvwRZDwkz+uidN/9+FSLnPucZxXOf89PVbQxecN6zq0tenF/ebn7vvy9ufs+H4iMwXLNhJs7S81uWA+CyGaQUfuvussuuLPatrVW/XVmdvNfCyXtf8S8u+qf1k6ZUV+hv7Hh2mfu7j7VdnpkmT5/zB8uluX/TnN+yPO95G6/af/N7P1x+zmuT31g7eaI4ceJkcaI8XnfkbcW73/3eej0TfstyvPx2ZlqCbEBcNoMUtguyoJqXw7899rfV0qsv/sVqfZ7NF5TdXZB9+UuPKH74Rx6c+4P9k6fP+YPl0vTf9LxznxlHm737fX9QrJWfa4i1INyH06QgYyyEWRqCbKDsmrEqiwZZmNcrQTkqP44fPx7ntXk7artzd3k7t7z972q2rPe+/w/jqKjia/okOAmv+hjuw9pkXZAxDt4FYLUE2cAJM7q2KcjKGAtfZ1VxhfttgqweTcaNmWlx/MTGHa/Zxxuvv/L84vWv/9Pi9EMPxJWtre3fH0elDWe5zae86ZVJeNXHcB/WJuuCjPGxa9Y9P9Q/cCE+w635BQC/BMAQ7V/bv+G21tzKqJq+/fiPf6t45KOOb1qfdwPa47czuyfIMtGE2XScQesSn4cvuvjO4uDBY3EGrJrfzuyOIMuQMAOgS02YeReA9giyjM2GmTgDoE3eBaA9gmwEmjCzawZAV4TZ3giykRFm0I4cf8MS2jAbZuJsMYJspIQZQ/TlLz0mjoC+a8LMrtliBNnIzYaZOKPPjvzbny6OHj0QZ8BQCLOdCTIqTZjZNYPtuSwmLG82zMTZhCBjE2EGQJeaMLNrNiHI2NJsmIkzAPYqvA3T9K0x9ijzXpbsivfOHJ8N72UZ/0u2PnHWJ896Xo2mLuRfjuKkWp08UJmZljavVGaWL/nVny2uufa/7v5q/RvOcptPedMrk/esrI/hPqw16ydPnihOnjhZXHvt72b5W5bey5JlTIfVTv7gDz8SRxu97IKfn3z9jZAgYynCbDwEWb3WrAsyxqKNyNoNQZbZZy/IVqsJs0Cc5UmQ1WvNuiBj6BYNrTYiazcEWWafvSBLx65ZnnIOsjAKf8X0w+unxCbAqqEgo99WvZvVBUGW2WcvyNITZnnpU5DddOPTiosu/kxx4MDxuLKgDWe5ySSMwl8x/fD6KbEJsGooyEijr7tZXRBkmX32gqw/hFke+hRkS9twlptMwij8FdMPr58SmwCrhoKM9uSwm9UFQZbZZy/I+sfPmQ2bIKvXmnVBxlbGtJvVBUGW2WcvyPrNrtnwCLJ6rVkXZONiN2t1BFlmn70gGwZhNhy5B9ms9VNiE2DVUJDlxm5W/wiyzD57QTYsXs7sP0FWrzXrgqy/7GYNmyDL7LMXZMNl16yf+hRkX/7SY4q/96PfjrNd2HCWm0zmnfzWT4lNgFVDQZaS3axxEGSZffaCbPiEWb/0Kci6uA7ZrPVTYhNg1VCQtc1uFrMEWWafvSDLhzDrB0FWrzXrgmx7drNYliDL7LMXZPnxc2Zp7S3I4trkgcrMtLR5pTKz3MWFYWetnxKbAKuG4w4yu1msgiDL7LMXZHmza7Z6fQqypW04y00m805+66fEJsCqYX5B9sxnPC2OdiayWAVBltlnL8jGQZitjiCr15r1vgfZoqF12+2fiqONDh9+eXHkyHVxBqsjyDL77AXZuHg5s3uCrF5r1lME2W52s7YKrUUJMlIRZJl99oJsvOyadUOQ1WvNeptBttfdrC4IMlIRZJl99oIMYdYuQVavNes7Bdkqd7O6IMhIRZBl9tkLMhpezmxHn4Is5YVhL37Vy6r5TvoYWbshyEhFkAkyRsCu2fL6FGRtX4fslRe+NI529o4b3lfenyxOVIV2vHjb1W8tbrnlQ9VjORFkpCLIBBkjIsx2b4hB9srDvxBH27vxpluqY30SLO/rj1o8NZ4M/4RhfSfIoCOCTJAxQsJscX0JsjqywsuVZ5a374alLd34rlvrQTi7hX/H3LPcZLEelff1R02QwUoJMkHGiPk5s511HWSLvmwYIuuKyy4o3vqvP1qcevChuDolnMlm/p71tblnucliPSrv64+aIIOVEmSCDCp2zeZbJsgueuUFcbSzG2+Ku1nzNP/CSJB1T5CRiiATZLCBMNuoCbJfe/XhuLKzG278wHocVYf1UKoH69PKxtkGMw8Jsu4JMlIRZIIM5hrLy5mLXjfryNtvLtbCVljYEQv/VPET5zGEJntkYX19tD5uBuvTysbZBjMPCbLuCTJSEWSCDHY0tF2zNi9OGr5BX3Hlb5ZfW2uCrBwLMuiGIBNksLDUYZbirXb2HGTNUqWebFiamW0w89D73vOM4hUX3lEcOHA8rkwJZ7LZf1WzNvcsN1msR+V9/VETZLBSgkyQwa61GWZ9f6udPgXZtsKZbPb5zdrcs9xksR6V9/VHbcsgC086IcigZYJMkMHStvs5sxS7WV0QZIIMVkGQCTJoxXScBUN/T8OGIBNksApjD7K1eARaECKsuQHAogQZAEBiggwYjLvvOhRHAHkRZMBgXHP1eXEEkBdBBgCQmCADBuPsc+6NI4C8CDJgT3ZzZYq9uuTS2+MIIC+CDGjfKisNIAOCDAAgMUEGAJCYIAMASEyQAYPhwrBArgQZMBguDAvkSpABACQmyAAAEhNkwEJcWgygO4IMACAxQQbsnu0ygFYJMgCAxAQZAEBiggwYjLPPuTeOAPIiyIDBuOTS2+MIIC+CDAAgMUEGAJCYIAMASEyQAXk5GY8AAyLIgMG4+65DcQSQF0EGDMY1V58XRwB5EWQAAIkJMmAwXBgWyJUgAwbDhWGBXAkyAIDEBBkAQGKCDAAgMUEGAJCYIAMGoy8XhvVmAEDbBBkwGC4MC+RKkAEAJCbIAAASE2QAAIkJMgCAxAQZkDe/EgkMgCADAEhMkAErtS8eAZgQZMBgnH3OvXEEkBdBBgzGJZfeHkcAeRFkAACJCTIAgMQEGQBAYoIMACAxQQYMxt13HYojgLwIMmAwrrn6vDgCyIsgAwBITJABg+HCsECuBBkwGC4MC+RKkAEAJCbIAAASE2QAAIkJMgCAxAQZMBguDAvkSpABg+HCsECuBBkAQGKCDAAgMUEGAJCYIAMASEyQAQAkJsgAABITZAAAiQkyYDDOPufeOALIiyADBuOSS2+PI4C8CDIAgMQEGQBAYoIMACAxQQYAkJggAwbj7rsOxRFAXgQZMBjXXH1eHAHkRZABACQmyIDBcGFYIFeCDBgMF4YFciXIAAASE2QAAIkJMgCAxAQZAEBiggwYDBeGBXIlyIBBePDoKcUN73h2nAHkRZABACQmyAAAEhNkAACJCTIAgMQEGbCt73//+3EEQFcEGQBAYoIMACAxQQYMwsMOHC/OPufeOAPIiyADBuFAGWSv+uVPxhlAXgQZAEBiggwAIDFBBgCQmCADAEhMkAGDcfddh+IIIC+CDBiEB4+eUtzwjmfHGUBeBBkAQGKCDBgEF4YFcibIgEFwYVggZ4IMACAxQQYAkJggAwBITJABACQmyIDBcGFYIFeCDBgEF4YFcibIAAASE2QAAIkJMmBb++IRgO4IMgCAxAQZAEBiggzojJc7ARYjyAAAEhNkwCA87MDx4uxz7o0zgLwIMmAQDpRB9qpf/mScAeRFkAHb84NgAJ0TZAAAiQkyAIDEBBkAQGKCDBiMu+86FEcAeRFkwCA8ePSU4oZ3PDvOAPIiyIBBeKAMsoMHH4ozgLwIMmBLhw+/vPiNq94SZ2k9ooyxs578jTgDyIsgAwbh1DLIXnHhHXEGkBdBBgCQmCADAEhMkAEAJCbIAAASE2TAIBw7tr/4yj2PjjOAvAgyYBDuv//U4l3vfFacAeRFkAE72BePAHRFkAEAJCbIAAASE2RAK7ywCbA8QQYAkJggAwBITJABACQmyIBBeMTBh4qznvyNOAPIiyADBuHUMsheceEdcQaQF0EGAJCYIAMASEyQAQAkJsgAABITZMAgHDu2v/jKPY+OM4C8CDJgEO6//9TiXe98VpwB5EWQAQAkJsiAQXBhWCBnggwYBBeGBXImyAAAEhNkAACJCTKA3ToZjwAtEWQAAIkJMmAQdnVhWDtYwMAIMmAQXBgWyJkgAwBITJABACQmyAAAEhNkAACJCTIAgMQEGQBAYoIMACAxQQYMxlP/wVfiCCAvggwYhNNP/17x0gvujDOAvAgyYEn74hGAvRJkAAvyFplAVwQZAEBiggwAIDFBBgzCsWP7i6/c8+g4A8iLIAMG4b6vP6q49QNPjzOAvAgyAIDEBBkwCKcePFb82FnfiDOAvAgyYBBcGBbImSADAEhMkAEAJCbIAAASE2QAAIkJMmAQXBgWyJkgAwZh4QvD7otHgAERZAAAiQkyAIDEBBkAQGKCDAAgMUEGAJCYIAMASEyQAQAkJsiAQTj14LHix876RpwB5EWQAYNw+unfK156wZ1xBpAXQQYswOXvAbokyAAAEhNkAACJCTIAgMQEGTAIx47tL75yz6PjDCAvggwYhPu+/qji1g88Pc4SOhmPAC0SZAAAiQkyYBBcGBbImSADWtbNNctcGBbImSADAEhMkAEAJCbIAAASE2QAAIkJMmAQXBgWyJkgAwahNxeGBeiAIAMG4ejRhxenHnwozgDyIsgAABITZAAAiQkyAIDEBBkAQGKCDAAgMUEGAJCYIAMG4bTTHizOPPNbcQaQF0EGDMIZj7u/eNFLPhdnAHkRZAAAiQkyAIDEBBkAQGKCDAAgMUEGDMKDR08p7vv6aXEGkBdBBgzCPfc8pvjwh54WZwB5EWQAAIkJMmAQXBgWyJkgAwbBhWGBnAkyAIDEBBkAQGKCDAAgMUEGAJCYIAMGwYVhgZwJMmAQXBgWyJkgAwBITJABACQmyAAAEhNkAACJCTIAgMQEGQBAYoIMACAxQQYMwmmnPViceea34gwgL4IMGIQzHnd/8aKXfC7OAPIiyAAAEhNkwMrsi0cANhJkAACJCTIAgMQEGTAIDx49pbjv66fFGUBeBBkwCPfc85jiwx96Wpx17GQ8AqyIIAMASEyQAYNw4MDx4qwnfyPOAPIiyIBBOPNJ3yzOf+Hn4wwgL4IMACAxQQaMkJ/aB/pFkAEAJCbIgMV43yOAzggyYBBcGBbImSADBuGv/uqxxX/86FPiDCAvggwAIDFBBgCQmCADAEhMkAEAJCbIAAASE2QAAIkJMgCAxAQZMAinnfY3xRN/+NtxBpAXQQYMwplP+mZx/gs/H2cAeRFkwPIyeH/Lk/FIURw+/PI4AlZt38lSHGdh3759xcf+6LY4g9U5/wXPLW67/VNxlofwDfq1V1xV7FtbK9bW9pX9VRZY+TUWvs7qcfhojkG4X5/UhzgOquduMm8t2uahTcKZrHl+M97y7DZ5YMOo/pi+K+/Lf8KwvKv+ORHmJ4pr3vbW4pZbPlQ9p+8WDa3/9J//W/FzP3tenBXFkSPXxRF072UX/Hz5tVV/3Y2RIIOWCLIg3K9P6kMcB9VzN5m3Fm3z0CbhTNY8vxlveXabPLBhVH9M35X3IcDCoBr1Jsh2s5sVQmsZTZwJM1ZBkAkyaIUgC8L9+qQ+xHFQPXeTeWvRNg9tEs5kzfOb8ZZnt8kDG0b1x/RdeR8CLAyqUedBtpvdrFURZqyCIBNk0ApBFoT79Ul9iOOgeu4m89aiqYcePHpKcf/9pxZnPO7+uDIjnMma5zfjLc9ukwc2jOqP6bvyPgRYGFSjpYJsFbtZq+DlTLokyAQZtEKQBeF+fVIf4jionrvJvLVo6qHP3PnE4i/+/O8Wv/KrfxpXZoQzWfP8Zrzl2W3ywIZR/TF9V96HAAuDarQhyBbV58hall0z2ibIBBm0QpAF4X59Uh/iOKieu8m8tWjqoVUF2c0331QeT1TfGE6ECAvjMsJOVMcTxfETx4vjx48Xf/PgA1mG1m4JM9oiyAQZtEKQBeF+fVIf4jionrvJvLVo6qG77zpUfOELj936WmThTNY8vxlPnd1+7/3vLaflQoiscl6Pp44xwl57xW+ESfVnwmPVKNxVf678J8bZm994hSCb4uVM9kqQCTJoxTiDrLqv4qc6Vvfrk/oQx0F89ox5a9E2DwW3fuD9dSuVt/DU+mw2OaXV0/K+/Pj1y66o1prH68eqYWXDn42nRUG2HLtmLEOQCTJohSCLg8mkPsRxEJ89Y+Parbf+Xnk/73nB1OmqHF5y6WVxshv1v6O6n/7XTS8IslYIM3ZDkAkyaIUgi4PJpD7E8Qd//9Z6MGXDySdMyudecsm/rOedqf/W6n7qfwBB1h1hxiIEmSCDVowxyD74wTKyYpSFtfK++nPV08pDOLk0x9e85tfL+ynxObXJqHv1Ka+6nzr7CbLu+TkztiPIBBm0IqcgCyHWOHToccX+hz2sDLK1OsSqWz3+tddcWj4jzOvnVuN4rAfxgWY6bcPapkc71MTW+rAiyFbLrhmzBJkgg1b0PcimI2snTWiEb5q/9dtvizEWbuViFWQxvOqF+lDpLsiOPnCg+M53H148/vHfiyvLamJrfVgRZGkIMxqCTJBBK1IF2aKhtUw89CnI/uLPzyj+7BN/p7j8ijvjyrKa2FofVgRZWl7ORJAJMmhFm0G2zG5WFwSZIEvBrtk4CTJBBq1YJMi63M3qgiATZCkJs3ERZIIM9uSCf/KS6vi97/3f4vFP+JFqvJWhfQMXZNNBdiIeBNmqCbNxEGSCDOZqQmsnf/zxT1THp//UU4o7P/uFapwLQSbI+qQJs0Cc5UeQCTJGZNHICprQWpQgCwQZq2HXLD+CTJCRgd3uZnVBkAWCjNUSZvkQZIKMnupyN6sLgiyYXq8G64f1pzQ2rG16dBNBxnaaMAvE2TAJMkHGCg0tsnZDkAXdBdkXv/jo4n9+/geLl/zCl+PKsupTXnU/dfYTZPmwazZMgkyQ0YI+vGSYmiALuguy9tSnvOp+6uwnyPIjzIZFkAkytpDzblYXBFkgyOifJswCcdZfgkyQjY7drG4IskCQ0W92zfpLkAmyLNjNSk+QBYKMYRBm/SPIBFmv2c0aDkEWCDKGRZj1hyATZCtnNytPgizoLsi++52HF8eO7S8O/dDRuLKs+pRX3U+d/QTZuDVhFoizNASZIGuN3axxE2RBd0H2J7c9sfjil36g+JVX/4+4sqz6lFfdT539BBkNu2ZpCDJBti27WSxKkAWCjHwIs9USZCMNMrtZtE2QBd0F2V9+/geLr33tkcV5z//ruLKs+pRX3U+d/QQZW2nCLBBn3RFkGQbZIx/5qDjbmtCibYIs6C7I2lOf8qr7qbOfIGMRds26I8gyDLLGHZ/+XBxB9wRZIMgYB2HWvrEH2Vo8ZiP8H7O5hW+Q4Qa0YJXNBD0XYjzcLr/8suoGe5VdkE0TZgB0aTbMxBnLyjrIGrNhJs4AaFMTZnbNWNYogqzRhJldMxiecGHYb/6fg3EG/SXMWMaogmyaMINh+exnfqj4yEeeFGd7NN6fG2aFZsNMnLGd0QZZQ5gB0KUmzOyasZ3RB1ljNszEGQBtE2ZsRZDNaMLMrhm0wWuDMM9smIkzBNk2hBkAXWrCzK4ZgmwBs2EmzgBomzAbt+zeOmlVmrdo8vZMNEKoZ/nWSb/ztmItvG1SePuksBjeNinc4jhcwn/9nZLCuDqsL5TKcf1R/M5vv7F+RvURVsKtvK8PU4P4yNT8O99+cXH06E8Uj3/CW6v5+iOTp5Ti3z+tXAirb7jqX8WF8pQXPsITp85+9ZkwLsTTYvmfYvVKuPPWSazY2N6eyXtZCrI9EWY0BFmwr/itt1xVzauvjep5zbH8xnLdkTgPH2El3KqHa+uD+MjU/GMfe2zxv/7yUcVrX/flar7+yOQppfrfukG5EFYvv/zymecGk4UrX/+m8j6eDgUZPdKEWZBznAkyQdYKYcYYg+wtb35DUT5Yjuu18Nh15TeMeljPq2dWx/Bnw0d1Vx+rW3lfH6YG8ZGp+Wc/8wPFPfccLF704q9X8/VHJk8pxb9nWrmwvjrvwagKtnX1afF1V/6mIKNXct41E2SCrFVNmAXibFzGEmRvfuOV5bhMnPL/16+/7vpNQVaP49dCnNfHUhjWd/WxupX39WFqEB+ZmtdmHg82PCX+PdPKhfXVeQ/OVZ8Wq0grh+E0+borrxJk9EaOYSbIBFlnmjgTZuOQe5C9qQyQtbX99Q5YDLKwnnOQVWKQXVbGWRVmr7tKkNEbTZgFQ48zQSbIOifMxiHnIHvTVVcU17/9+irIQtyMMcjCSjiG3347cfx4cf/93xZk9MrQd80EmSBbmSbMAnGWn1yD7KGHHirefv3bi7X9a4KsOl2erILsoot+SZDRS0MNs7EHmeuQrVB1Uo+38M073MhDjjHWOHLkSBwBQxD+QyHcXM9sWARZIsKMIQgn9Te84Q1xBgzJbJiJs37zkmVPeDlz2HLeIQvCSyDXXnutlyzL4yt/6Z9Xy+EbHQxNn1/O9DNkgqx3mjgTZsORe5AFzYn8+uvePsogu/DCX6ymQowc9DHMBJkg6y1hNhxjCLJGcyJv3PDvbqyCJ8cgO/yKl1fHQIiRo+mv59RxJsgEWe8Js/4bU5DNmg204Kabbq7//7a81QEWkygM67v6WN3K+/owNYiPTM1rM48HG54S/55p5cL66rwHo8Ov+GdxNCHCGJPUu2aCTJANRhNmgTjrlzEH2TzzIq3xnve8v+qgFEF2+PBkx2uW+IJaqjATZIJskOya9YsgW9x2sdY10QWLm/5aXUWcCTJBNmjCrB8EGZCzVeyaCTJBlgUvZ6YlyIAx6DLMBJkgy45ds9UTZMCYdBFmgkyQZUuYrY4gA8aozZ8zE2SCLHvCrHuCDBi7ve6aCTJBNhp+zqw7ggygtmyYCTJBNkp2zdolyAA22u3LmYJMkI2aMGuHIAPY2iK7ZoJMkFHycubeCDKAnW0XZoJMkDHDrtnuCTKAxc17OVOQCTK2IMwWJ8gAltPE2de++teCLI5hLi9n7kyQAWztJ596Vhxtbew5IsjYFbtm8wkyYGwWiayG1NiZIGMpwmwjQQbkYtHQkg/tEmTsiTCrCTKgz+xm9Z8goxVj/zkzQQasmsjKiyCjdWPcNRNkQFu8ZDhOgozOjCnMBBmwHbtZ7ESQ0bkxvJwpyGCc7GbRFkHGSuW6aybIIB92s0hBkJFEbmEmyKD/7GbRZ4KMpHJ5OVOQQRp2s8iFIKM3hrxrJsigXXazGBtBRu8MMcwEGezMbhZsTZDRW0MKM0HGmNnNgr0TZPTeEH7OTJCRG7tZsFqCjEHp666ZIGMo7GZBPwkyBqlvYSbISMluFgyfIGPQ+vJypiCjbSILxkWQkY2Uu2aCjEV5yRCYR5CRnRRhJsjGzW4WsFeCjGytMswEWZ7sZgGrIsjI3ip+zkyQDYfdLKCPBBmj0tWumSBLz24WMGSCjFFqO8wEWTfsZgFjIcgYtbZezhRku2M3C2AjQQbRXnbNBJndLIC9EGQwY5kwyznI7GYBdE+QwRZ283Lm0ILMbhZAvwgyWMBOu2Z9CTK7WQDDJMhgF7YKsy6DzG4WQP4EGSxhNsx2G2QiC4Bpggz2YPrnzEKQeckQgGUIMmhJiDNfTgAsYy0egT0SYwAsS5ABACQmyAAAEhNkAACJCTIAgMQEGQBAYoIMACAxQQYAkJggAwBITJABACQmyAAAEhNkAACJCTIAgMQEGQBAYoIMACAxQQYAkJggAwBITJABACQmyAAAEhNkAACJCTIAgMQEGQBAYoIMACAxQQYAkJggAwBITJABACQmyAAAEhNkAACJCTIAgMQEGQBAYoIMACAxQQYAkJggAwBITJABACQmyAAAEhNkAACJCTIAgMQEGQBAYoIMACAxQQYAkJggAwBITJABACQmyAAAEhNkAACJCTIAgMQEGQBAYoIMACAxQQYAkJggAwBITJABACQmyAAAEhNkAABJFcX/ByvzfazOX/K3AAAAAElFTkSuQmCC)\n","\n","\n","\n"],"metadata":{"id":"GwsUVLf4_EXT"}},{"cell_type":"code","source":["voxel_length = 0.1  # m\n","voxel_volume = voxel_length**3  # m^3\n","garden_hose_diameter_in = 0.75  # in\n","fire_hose_diameter_in = 2.5  # in\n","meter_per_in = 0.0254  # conversion ratio\n","\n","garden_hose_radius = garden_hose_diameter_in * meter_per_in * 0.5  # m\n","heated_column_volume = np.pi * garden_hose_radius**2 * voxel_length  # m^3\n","heated_ratio = heated_column_volume / voxel_volume  # m^3\n","\n","voxel_volume_temp_simple_average = (\n","    (target_temp2_c * heated_column_volume)\n","    + (starting_temp_c * (voxel_volume - heated_column_volume))\n",") / voxel_volume  # C\n","simple_average_temp_delta = voxel_volume_temp_simple_average - starting_temp_c  # C\n","\n","print(\"simple average across a %dcm voxel for a garden hose plume is %.1f C\" %(voxel_length*100,simple_average_temp_delta))"],"metadata":{"colab":{"base_uri":"https://localhost:8080/"},"id":"JNFJv5L4VAhz","executionInfo":{"status":"ok","timestamp":1719962104764,"user_tz":420,"elapsed":1,"user":{"displayName":"Justin Tucker","userId":"10374969465321862369"}},"outputId":"aec2c6f5-5428-47ca-99cf-03eb7a9abb67"},"execution_count":7,"outputs":[{"output_type":"stream","name":"stdout","text":["simple average across a 10cm voxel for a garden hose plume is 2.1 C\n"]}]},{"cell_type":"markdown","source":["## Concept B - heated blob\n","\n","Heat and deposite a volume of water into the tank such that 1 or more entire voxles see an elevated temperature of a significant temperature delta relative to the background water.\n","\n","Could be as simple as filling a cooler with hot water from the tap, sinking it to the bottom of the tank and opening the lid.\n","\n"],"metadata":{"id":"wSEVnGBdotH4"}},{"cell_type":"markdown","source":["### Power calcs for small Immersion cube"],"metadata":{"id":"Ajy5HUNm_AWL"}},{"cell_type":"code","source":["heated_volume = 4 / 1000  # m^3\n","energy_to_heat = heated_volume * rho_sea_water * delta_t2 * cp  # J\n","heater_power = 1000  # W\n","time_to_heat_min = energy_to_heat / heater_power / 60  # min\n","heated_voxels = heated_volume / voxel_volume"],"metadata":{"id":"-839Ov8EVjJs","executionInfo":{"status":"ok","timestamp":1719962111676,"user_tz":420,"elapsed":312,"user":{"displayName":"Justin Tucker","userId":"10374969465321862369"}}},"execution_count":9,"outputs":[]}]}