Summary
在这里，我们提出了一个协议，当水被引入到盐水的冷浴，作为二次制冷剂，在一定温度范围内大大低于水的凝固点以证明冰的生成。它可以用作制造用于工业冰的另一种方法。
Abstract
我们证明用于在过冷盐水环境中的热与质量传递的冷冻现象和研究的方法。我们的试验表明，在适当的条件下，可以在水被引入到冷盐水浴来制造冰。以制造冰的形式，除了具有盐水和水混合，传热的速率必须绕过传质。当在微小液滴的盐水表面的形式被引入的水，热和质量传递的方式是通过扩散。浮力从下方与盐水混合停止水，但随着冰成长厚，它减慢传热率，使得冰更难以成长的结果。当在流动的形式盐水内引入水，许多因素被发现影响多少冰能形成。盐水洗涤温度和浓度，它们的热与质量传递的驱动力，分别可影响水 - 冰转换拉提O;低浴温及盐水浓度鼓励更多的冰形成。流动流变性，它可以直接影响同时传热和传质系数，也是一个关键的因素。此外，流动流变学变化与散装流体的流动的接触面积。
Introduction
冰浆在工业中广泛使用，并且一种特别成功的应用是冰清管技术1，2。相比于常规的泡沫和固体猪，冰猪可以在很长的距离，通过复杂的拓扑行进的，因为液相的润滑效果和它的凝固点的，因为一些冰晶融化3，4，5的高度。即使猪被卡住，可以简单地等待冰浆熔化并在稍后继续清洁过程。管道清洗的这种方法是便宜的，并易于使用。
冰馏分起着冰猪的性能起关键作用。为了测量冰分数，可以使用一个壶煮（法国媒体），以确定是否冰浆是足够厚6，“> 7，一个高的壶煮冰馏分，通常为80％时，进行冰清管时需要。网上冰分数检测最近的研究表明，电磁和超声波适于作业8，9，10，11。
冰猪通常是由一个刮面制冰机从5％（重量）NaCl溶液（盐水）制成。它也是在工业制造冰浆的主要方式。这种类型的制冰机的冻结水或盐水在冷金属表面，通常是光滑的316不锈钢表面，然后循环剪冰粒关闭。液体-金属界面是非常复杂的，并通过广泛的那些对制冰12基本因素的影响。非金属和水之间的接口是非常不同的，一个特别有趣的例子是高岭石。该Kaolinite -水界面是特殊的，因为没有一个有利的冰结构邻近于固体表面，而是鼓励冰状氢键团簇形成在它的13顶，14两性基底流体的层。制作冰猪的另一种方式，而需要高浓度盐水同时加入粉碎预制冰块。对于此方法，所述制冷系统能以高得多的蒸发温度运行，因为没有降凝剂（FPD）之前被冰的形成增加;它因此认为更有效的，由于对于给定的冷却负荷15，16，17中的降低压缩率和减少功耗。
还有另外两个冰制作方法：从生产过冷的水冰，并把制冷剂与水直接接触
冰的形成需要的热与质量传递由于在过程中产生的熔化潜热。它首先是由Osborn的雷诺兹发现在1874年该热质在气体中的运输是强耦合，并且可以在类似的数学式20来表示。这项工作形成声势，热和质量的流体传递的主题开拓纸和被转载多次21,22。这个题目，然后通过一个研究他人的数目，从既分析和经验方法，对气体，液体，熔融金属23，24，25，26，27，28，29，30，31，32，33。除了热与质量传递，流体需要成核位点，其中树枝状冰增长可以开发。一个现代的洞察冰晶的生长采用构形法，由Adrian Bejan开发，解释为什么冰生长在这种方式34，35，36。
冰的形成在盐水不同于在纯水中非常不同，由于盐的存在。首先，盐变化的流体的热力学和按压其冰点。其次，盐可以不溶解在冰基质（除hydrohalite，当温度达到共晶点，这只能形成），并且当冰开始增长它被拒绝到散装流体。盐的排斥海冰和冰都在实验室中37，38研究发现。由于拒绝高浓度盐水的温度大大低于海水的冰点，因为它下降时，冰生长在流动盐水和静态散装流体之间的界面。这些冰的钟乳石，也叫brinicles，在麦克默多海峡，南极洲被首次发现并进行了实验39，40，41，42研究。 2011年，BBC摄制的冰冻星球系列brinicles形成“外部参照”> 43，44。
在我们的实验室中，人们发现，通过反转时水被引入到冷盐水浴中流动的静止流体，水可以转变成冰的正确条件45下。结果发现，其中，被引入的水的位置，流动流变学，和盐水洗涤温度和浓度是影响多少冰如何可以产生的所有的关键因素。本研究的总体目标是调查，如果一个制冰机可以通过这一机制被开发，以产生冰浆，考虑到提高的蒸发器的温度和液体到液体传热的高速率可以提高能源的使用效率。本文股份实验的关键方面。
Subscription Required. Please recommend JoVE to your librarian.
Protocol
注意：有两个有毒化学品，甲醇和乙二醇，在这些实验中使用。甲醇可以在人体内被代谢以产生甲醛，然后向甲酸或甲酸盐。这些物质是有毒的中枢神经系统，甚至可能导致死亡。乙二醇可被氧化成乙醇酸，然后可以变成草酸。这可能会导致肾功能衰竭而死亡。不要喝这些化学物质。立即咨询医生，如果发生事故。
1.冷却系统
注意：这是非常困难的，以保持盐水在-18℃左右，当环境温度是大约在室温下。重要的是存储乙二醇和盐水罐是公绝缘的合理大小的和，以避免多余的电力消耗，以确保最佳的系统性能。它建议在罐的大小不超过30 L。
- 准备二次冷却液
- 倒入1升乙二醇进入二次冷却罐，罐A（基数：400毫米×200毫米，高度：350毫米）。水约0.6-0.65 L（600-650 G）添加到油箱A.
- 重复步骤1.1.1几次，直到有坦克A（25L）足够的液体。
- 搅拌流体，使该流体是均匀的。
- 开关在箱体A两个泵的全容量设置（2,500升/小时）。确保被困在换热器所有的泡沫和管道被释放。
- 关掉泵，以观察是否所有气泡被释放。如果没有，重复步骤1.1.4。
- 盐水的制备
注意：在这个例子中，22％（重量）的盐水中制备。如果需要其他浓度，补充盐分的质量也应该相应改变。参考盐水的浓度和密度值可以在页D-257的第64 个版本（1983年）中找到化学与物理 46 电子CRC手册 。- 加入4公斤的水倒入5升的塑料烧杯中。
- 测量电子秤1公斤氯化钠盐，倒入这种盐与水的烧杯中。
- 搅拌混合物直至溶液澄清（ 即，不存在盐粒子或水中的气泡在溶液中可见）。
- 取一个样品，〜10毫升，用10毫升的注射器的溶液。
- 注入流体进入U型管密度计。
- 检查在管气泡。如果有任何，注入更多的流体推他们出去。
- 按“快速设置”，然后选择“密度温度。” Type在20°C，按“确定”。密度计现在将测量在该温度下的流体密度。
- 按下启动并等待结果。
- 密度读数与1164.00千克/ m 3的比较。
- 添加更多的盐，如果读的是下面的试比较氮密度。否则加水。
- 重复步骤1.2.3-1.2.10直到流体密度是正确的（1164.00公斤/米3）。
- 倾该溶液成一个较大的容器中，容器A.
- 使用步骤1.2.1-1.2.12，使35-40大号盐水，并把容器成冷柜在-40°C。保持盐水有48-72小时，直到温度达到-19.18°C（冻结这22重量％的盐水点）。
2.冰的注射和洗水的制备
- 准备冰为水注入
- 倒入1升水到一个小容器（200×200×50毫米）。
- 重复步骤2.1.1与另一容器，并放置在冷柜的两个容器在-40℃。
- 保持它们在冷冻10小时或更多，以确保所有的水被冻结。
- 准备洗涤水的冰壳
- 填有5升烧杯中以5升的水。
- 补的2L烧杯用2L的水。
- 放置两个烧杯放入卧式冷冻在-40℃下为8-10 h后，周围有未冻水冰包裹的厚壳。
- 在3-5米/秒的速度从水龙头使用高速水射流在冰壳的顶部打开一个3厘米直径的孔。
- 沥干水冰壳里面。
- 把两个烧杯放回冷冻。
- 如果冰壳的质量没有达到3公斤和1公斤的两个烧杯，分别重复步骤2.2.1-2.2.5，但在步骤2.2.3保持烧杯在冷冻更长。两个烧杯现在应该能够以包含2 L和1L水，分别
3.水简介地位和流变控制实验
- 在卤水表面水介绍
- 从容器中倒出2升22％（重量）的冷盐水进入冰淇淋机和开关冷却装置上的铝桶。
- 测量用温度计/热电偶盐水的温度（无论是K型或T型是合适的）。在实验进行，如果盐水是-15℃或更低。
- 在室温下填充用自来水100ml的玻璃注射器。附上一个2毫米内径，1毫米厚，和1米长的硅胶管与注射器的前端。
- 放置注射器在特定的位置，使得在注射器中的水和硅氧烷管的出口之间的头部。静水压力将挤压出来的水管。
- 淹没硅胶管的一定的长度，通常为70厘米，进入盐水洗涤。
- 调节注射器和管出口之间的相对位置，使得流体静压足够大以允许水离开注射器。如果管被堵塞，通过提高注射器向更高的垂直位置，UNT增加头部IL静水压力能够克服所述管内的剪切应力。
- 保持管出口大致为1cm以下的盐水表面上方。
- 调节浸没管的长度和注射器高度，以控制水的出口让温度，以便确定多少冰能够制造或者多少混合在盐水表面发生流速。冻结现象应现在盐水表面进行观察。请参见参考45进一步的方向。
- 通过盐水引入水
- 重复步骤3.1.1-3.1.6。
- 保持管出口处的盐水内，优选在容器的底部。
- 调节浸没管的长度和注射器的高度。
- 调整管出口的角度，以控制流的流变性。
- 重复步骤3.2.3-3.2.4找到最好耦合流动流变和流动，可以产生最冰速率。
- 制冰
- 如果有在管道中的气泡，打开箱A内的两个泵以释放气泡出乙二醇循环系统，然后关掉泵。
- 开关上的三个制冷机组，并让他们运行10-16小时，以降温乙二醇溶液。
- 测量用温度计/热电偶的乙二醇溶液。乙二醇的温度应在约-25℃。
- 测量在容器A中的盐水的温度，以确保它是在-19℃进到步骤4.1.5之前。
- 填写卤水罐，罐B，从容器中大约30升盐水和交换机上的坦克A.两个泵
- 测量罐A中的乙二醇如果温度低于-19℃的温度，关闭一个或多个冷却单元，以防止冰粒子的析出的热交换器外我Ñ罐B.如果温度高于预期盐水温度高，打开所有的三个冷却单元。开展在-17℃至-19℃的实验。
- 放置冰从步骤2.1在两个预制块到绝缘5升的烧杯中，容器B和倒大致3升的水进入烧杯中。
- 测量水的温度和由如果温度上升实验之间搅拌该混合物保持在2℃。
- 填用100mL的2℃水的玻璃注射器。
- 申请5-10毫升甲醇中，罐B的玻璃窗口以停止缩合和冰的形成。
- 通过调节注射器和管的出口之间的相对位置，以便有一个恒定的静水压力，因此以恒定流速注入水进入盐水洗涤。约70厘米硅胶管，应在盐水被淹没。调整喷射的角度为0°，使得在向上可怕的初始水流速度ction为0米/秒。
注：注射器可以是手持或夹持到一个支架上。当盐水温度是较冷的，因为它需要更多的时间进行调整的支架，和冰可能堵塞管子手持是比较合适的。保持流动流变通过确保以恒定的流速和喷射角（0°）和通过保持冻结前沿的管出口上方大致3厘米整个实验一致。不要让流进它开始转向动荡47的区域。请参见参考45进一步的方向。 - 如步骤4.2和4.3描述收集冰块。在不同温度下的盐水4.1.11 - 重复步骤4.1.8。
- 收集制作冰块和估计有多少冰生产（干集合）
- 通过按比例把一个容器（200×200×50毫米）和零读数的“打开”按钮。
- 使用筛舀出冰抖落盐水。
- 把这个冰在容器中。测量用的刻度冰的质量。
- 冰融化后，用10毫升的注射器取样品。注射液这个样品放入密度计。
- 执行步骤1.2.6-1.2.9。
- 记录密度读数。
- 从其密度计算净水的质量（ 即，水转化成冰的质量）使用下面的公式：
哪里
是由重量％的盐水的浓度和
和
分别是盐和水，群众。
- 收集制作冰块和估计有多少冰产生（湿集合）
- 填用冰壳（步骤2中的5升烧杯中。2）和室温的自来水。把它放回冷冻在-40°C。
- 倒出从5升烧杯中的冰壳的水倒入2升烧杯中，当它的温度在0℃。填补了5升烧杯中。保留两个烧杯在冷冻。
- 舀出在步骤4.1.8和4.1.9中产生的冰和倾200-500毫升从2升烧杯中的水到冰清洗。应用0℃的水之前，不要摇动筛。
- 抖落在筛中的流体。
- 重复步骤4.2.2-4.2.7。
是由重量％的盐水的浓度和
和
分别是盐和水，群众。 Subscription Required. Please recommend JoVE to your librarian.
Representative Results
图1比较的水在盐水表面水通过盐水注入引入的影响。在“冰盖”的情况下，形成的冰是固体，因为水不与散装流体混合得多。两种流体之间的温度和密度差异对水产生浮力，防止它们混合。两种流体是静态的（ 即，该传热比质量大得多;钪≈500，镨≈10，以及乐≈50），因此，冰能形成容易。既没有形成在这个实验糊状层也不盐排斥。一旦冰成长厚，它会阻碍热传递的速率，由于其低导热性，影响冰的形成速率。在这一点上，可以清楚地观察到，引入了“甜水”不再能够及时冻结成固体。此外，在没有对流，低热表面盐水本身l的电导率也阻碍了潜热从冷水槽的运输。的冰形成的速度直接与相关联，并且以盐水温度非常敏感。例如，在-15℃盐水洗涤水结冰比-13℃盐水快得多。在水喷射的情况下，冰的形状和大小是相关于流动流变性。在图1所示的冰的杆具有两个独特的部分：直链头后跟一个卷曲的尾巴。卷曲部形成更接近盐水表面，其中该流具有更多湍流给它。卷曲尾巴通常比因为紊流的发病，最大限度地减少热与质量传递速率之间的差的直头薄得多，特别是在流，其中所述热量和质量传输是相同的外层。因此，只有内芯可以冻结成冰块。如果管出口保持水平而不是垂直向上，冰WIL的片材生成湖冰的生成变得更加稳定，其结果是可重复的。最后，我们发现降低流速是不消除混合的有效方法。相反，它显著增加阻塞该管的可能性。
进行水到冰的转化率测量时，水喷射角被保持在0℃于水平轴。盐水洗涤温度和浓度的影响在图2中示出。转换率通常在0.4到0.9报考研究盐水温度和浓度。重要的是要保持冰的形成前沿恒定的流量流变性和位置在整个实验是重要的。大体积在箱体B盐水有助于减少对测量局部热梯度的影响。盐水温度和转化率之间的关系是对所研究的温度响第一阶即对于最佳拟合线的系数列于表1。如果使用不同的喷射角，水到冰转化率将不再按照下列关系，因为接触，因此面积，热和质量传递的速率是不同的。当收集冰中，保持施加到抖落盐水/洗涤水相一致，并尽量减少的水留在筛上的量的力是非常重要的。相似量的用于洗掉盐水应适用于避免不一致的结果。结果发现，如果超过500毫升水用于洗涤冰，任何进一步的盐度减少是不太可能发生。当体积低于200毫升，盐度可以高达4％（重量）。
由于蒸发器的温度高于一个刮面制冰机，它通常采用-40℃，如果这种方法用于生产冰高得多，较高的COP预计根据我们在图3中的计算。如果例如，蒸发器温度上升至-20℃，缔约方会议几乎可以使制冷剂R134A达到3。
图1：水介入位置。当在盐水表面引入的水的“冰冠”能形成。冰形式一根杆，当管出口保持直立。当水在盐水注入，冰的形状取决于流动流变性。 请点击此处查看该图的放大版本。
图2：在一个最佳拟合线不同浓度的盐水转化率比较。这两个布林È温度和浓度的影响多少水可冻成冰（转化率）时的流速和流变性被保持不变。转化率与盐水温度降低线性增加。在较低浴温盐水较低浓度产生更多的冰。洗涤方法收集较干收集方法更多的冰。 请点击此处查看该图的放大版本。
图3：在不同的蒸发器的温度为一系列冷却剂的性能系数。较高的蒸发器温度青睐的性能（COP）的冷却系统的系数。这两个过渡性制冷剂（R22和R134A）具有比已经禁止R502和混合物（R404A和R507A）的COP。请点击此处查看该图的放大版本。
| 盐浓度（重量％） | 收集干 | 湿集合 | ||
| P1 | P2 | P1 | P2 | |
| 23.3 | -0.09909 | -1.34 | -0.1196 | -1.439 |
| 22 | -0.1204 | -1.633 | -0.1439 | -1.839 |
| 21 | -0.1261 | -1.682 | -0.1545 | -1.98 |
表1： 系数FOr处的最佳拟合线的转换比与盐水温度图。转化率线性地根据下式与盐水温度相关： 
。无论干法和湿法收集方法在这里列出。
Subscription Required. Please recommend JoVE to your librarian.
Discussion
用盐水作为二次制冷剂冰生成的过程涉及的热与质量传递的组合。当热传递较大，则水之前冰形成具有与散装流体混合的机会。据观察，当存在引入的水和静止本体盐水之间的相对运动（即盐水内注入水）时，流程有助于传热，并鼓励冰以形成迅速。然而，当存在太多的湍流中的流动，可以产生没有冰。这种技术的最大限制是盐水的混合和稀释。盐水量将继续作为进程继续上升。因此，使得冰时这种方法，要意识到的上升盐水体积和滴盐水盐度的是重要的。此外，有人指出，如果没有收集所生成的冰，它会熔化。这可能是因为盐水是不在其熔化温度，从而允许热和形成的冰和散装流体之间的质量传递。的热与质量传递的方式是通过只扩散，和熔化的速度是缓慢的。然而，由于冰盐水表面上漂浮，来自周围环境的附加热量侵入提高冰融化的速率。由于这个原因，所产生的冰应迅速一旦被产生，以避免在盐水的体积进一步增加收集。
减少稀释或分离水和盐目前正在在我们的实验室研究。一项所述的许多思想是将注入的水重新引入到另一管是直径较大，使得水将只暴露于散装流体为很短的时间内，减少在二次制冷剂的量的变化。当水被暴露于盐水洗涤，随后在更大的管中的冰的生长的完成将发生冰核。通过加入这种固体表面，所产生的冰的堆积盐度是可控的。例如，如果需要在冰较低盐含量，可以向流体在二级管添加更多的“甜水”。该二次管的浸入长度可以容易地改变，这取决于产品所需的冰级分。
流动流变学对接触的表面积和在散装流体流动的面积 - 体积比的显著影响。我们的观察表明，接触的较大面积是鼓励更多的冰，以形成更有利的。接触的面积增加也应提高质量传递，但是还没有在所研究的盐水的温度和浓度范围内被观察到。看来，流进入过渡区，在动荡，流量分离开始发生之前，冰将永远被创建。如果流分离和大湍流存在，水分子的每个簇都需要自己的成核点，和冰可能无法在这些情况下形成。
“>盐水温度和水到冰的转化率之间的关系是线性的，而在恒定盐水浓度。该转换比与盐水温度最佳拟合线的偏移表明盐水浓度也起着重要的作用冰形成/水稀释过程。由于相变，该边界条件是在常规热和质量传递类比研究非常不同的，因此，那些类比不足以用于说明这种情况。这项研究还表明，由于该冷冻前沿可以被固定到从所述管的出口处的相对稳定的距离，则流程可以达到一个稳定状态。这表明，这种现象可以被用作冰生产行业一个可靠的新机制，因为高得多的蒸发器温度和COP预期相比于现有制冰技术。
Subscription Required. Please recommend JoVE to your librarian.
Disclosures
作者什么都没有透露。
Acknowledgments
作者没有确认。
Materials
| Name | Company | Catalog Number | Comments |
| DMA 4500 M | Anton Paar | 81546022 | Density Metre |
| GELATO Chef 2200 | magimix | 0036500504R13 | Ice Cream Maker |
| 280D | FREEZE MASTER | 241-1441 | Pipe Freezer |
| M17.5X2 | BLUE ICE MACHINES | GK924 | Slushy Puppy Machine |
| HH68K | OMEGA | 140045 | Thermometer |
| OHAUS | TS4KW | 1324 | Scale |
| ZFC321WA/BNI225 | ZANUSSI | 920672574-00 | Freezer |
| EIS Heater Matrix | Vauxhall | 214720041 | Heat Exchanger |
| 2500LPH | JBA | AP-2500 | Pump |
| Glass syringe | FORTUNA Optima | 100 mL | |
| OAT concentrated coolant | wilko | P30409014 | Ethylene Glycol |
| pure dried vacuum salt | INEOS Enterprise | 1433324 | NaCl Salt |
| Methylated Spirits | Barrettine | 1170 | Methanol |
References
- Cleaning and separation in conduits. UK patent. Quarini, G. L. , GB2358229, WO0151224 (2001).
- Quarini, J. Ice-pigging to reduce and remove fouling and to achieve clean-in-place. Appl. Therm. Eng. 22, 747-753 (2002).
- Evans, T. S., Quarini, G. L., Shire, G. S. F. Investigation into the transportation and melting of thick ice slurries in pipes. Int. J. Refrig. 31, 145-151 (2008).
- Shire, G. S. F., Quarini, G. L., Rhys, T. D. L., Evans, T. S. The anomalous pressure drop behaviour of ice slurries flowing through constrictions. Int. J. Multiph. Flow. 34, 510-515 (2008).
- Shire, G. S. F., Quarini, G. L., Evans, T. S. Pressure drop of flowing ice slurries in industrial heat exchangers. Appl. Therm. Eng. 29, 1500-1506 (2009).
- Evans, T. S. Technical Aspects of Pipeline Pigging with Flowing Ice Slurries [dissertation]. , Department of Mechanical Engineering, University of Bristol. (2007).
- Shire, G. S. F. The behaviour of ice pigging slurries [dissertation]. , Department of Mechanical Engineering, University of Bristol. (2006).
- Hales, A., et al. Ice fraction measurement of ice slurries through electromagnetic attenuation. Int. J. Refrig. 47, 98-104 (2014).
- Hales, A., et al. The effect of salinity and temperature on electromagnetic wave attenuation in brine. Int. J. Refrig. 51, 161-168 (2015).
- Hales, A. Ice slurry diagnostics through electromagnetic wave attenuation and other techniques [dissertation]. , Department of Mechanical Engineering, University of Bristol. (2015).
- Lucas, E. J. K., Hales, A., McBryde, D., Yun, X., Quarini, G. L. Noninvasive Ultrasonic Monitoring of Ice Pigging in Pipes Containing Liquid Food Materials. J. Food Process. Eng. 40, e12306 (2015).
- Carrasco, J., Hodgson, A., Michaelides, A. A molecular perspective of water at metal interfaces. Nat. Mater. 11, 667-674 (2012).
- Hu, X. L., Michaelides, A. Ice formation on kaolinite: Lattice match or amphoterism? Surf. Sci. 601, 5378-5381 (2007).
- Hu, X. L., Michaelides, A. The kaolinite (0 0 1) polar basal plane. Surf. Sci. 604, 111-117 (2010).
- Leiper, A. N., Ash, D. G., McBryde, D. J., Quarini, G. L. Improving the thermal efficiency of ice slurry production through comminution. Int. J. Refrig. 35, 1931-1939 (2012).
- Leiper, A. Carnot cycle optimisation of ice slurry production through comminution of bulk ice [dissertation]. , Department of Mechanical Engineering, University of Bristol. (2012).
- Leiper, A. N., Hammond, E. C., Ash, D. G., McBryde, D. J., Quarini, G. L. Energy conservation in ice slurry applications. Appl. Therm. Eng. 51, 1255-1262 (2013).
- Bédécarrats, J. -P., David, T., Castaing-Lasvignottes, J. Ice slurry production using supercooling phenomenon. Int. J. Refrig. 33, 196-204 (2010).
- Wijeysundera, N. E., Hawlader, M. N. A., Andy, C. W. B., Hossain, M. K. Ice-slurry production using direct contact heat transfer. Int. J. Refrig. 27, 511-519 (2004).
- Reynolds, O. On the extent and action of the heating surface of steam boilers. Proc. Lit. Philos. Soc. Manch. 14, 7-12 (1874).
- Reynolds, O. Papers on mechanical and physical subjects: reprinted from various transactions and journals. , Cambridge University Press. Vol. I, 1869-1882, Ch. 14. On the extent and action of the heating surface of steam boilers 81-85 (1900).
- Reynolds, O. Papers on mechanical and physical subjects. Int. J. Heat Mass Transfer. 12, 129-136 (1969).
- Prandtl, L. Eine Beziehung zwischen Wärmeaustausch und Strömungswiderstand der Flüssigkeiten (On the relation between heat exchange and stream resistance of fluid flow). Physik. Z. 11, 1072-1078 (1910).
- Prandtl, L. Bemerkung über den Wärmeübergang im Rohr (Note on heat transmission in pipes). Physik. Z. 29, 487-489 (1928).
- Taylor, G. I. Conditions at the surface of a hot body exposed to the wind. Rep. Memo. ACA. 272, (1916).
- Taylor, G. I. The Application of Osborne Reynolds' Theory of Heat Transfer to Flow through a Pipe. Proc. R. Soc. A. 129, 25-30 (1930).
- Kármán, T. v Proceedings of the Fourth International Congress for Applied Mechanics. , Cambridge, UK. 54-91 (1934).
- Kármán, T. v The analogy between fluid friction and heat transfer. Trans. Am. Soc. Mech. Eng. 61, 705-710 (1939).
- Martinelli, R. C. Heat transfer to molten metals. Trans. Am. Soc. Mech. Eng. 69, 947-959 (1947).
- Colburn, A. P. A method of correlating forced convection heat-transfer data and a comparison with fluid friction. Trans. Am. Inst. Chem. Eng. 29, 174-210 (1933).
- Colburn, A. P. A method of correlating forced convection heat-transfer data and a comparison with fluid friction. Int. J. Heat Mass Transfer. 7, 1359-1384 (1964).
- Chilton, T. H., Colburn, A. P. Mass Transfer (Absorption) Coefficients Prediction from Data on Heat Transfer and Fluid Friction. Ind. Eng. Chem. 26, 1183-1187 (1934).
- Friend, W. L., Metzner, A. B. Turbulent heat transfer inside tubes and the analogy among heat, mass, and momentum transfer. AIChE J. 4, 393-402 (1958).
- Bejan, A. Constructal-theory network of conducting paths for cooling a heat generating volume. Int. J. Heat Mass Transfer. 40, 799-816 (1997).
- Bejan, A., Lorente, S. Constructal theory of generation of configuration in nature and engineering. J. Appl. Phys. 100, 041301 (2006).
- Bejan, A., Lorente, S., Yilbas, B. S., Sahin, A. Z. Why solidification has an S-shaped history. Sci. Rep. 3, 1711 (2013).
- Lake, R. A., Lewis, E. L. Salt rejection by sea ice during growth. J. Geophys. Res. 75, 583-597 (1970).
- Wettlaufer, J. S., Worster, M. G., Huppert, H. E. Natural convection during solidification of an alloy from above with application to the evolution of sea ice. J. Fluid Mech. 344, 291-316 (1997).
- Paige, R. A. Stalactite Growth beneath Sea Ice. Science. 167, 171-172 (1970).
- Dayton, P. K., Martin, S. Observations of ice stalactites in McMurdo Sound, Antarctica. J. Geophys. Res. 76, 1595-1599 (1971).
- Eide, L. I., Martin, S. The formation of brine drainage features in young sea ice. J. Glaciol. 14, 137-154 (1975).
- Martin, S. Ice stalactites: comparison of a laminar flow theory with experiment. J. Fluid Mech. 63, 51-79 (1974).
- Jeffs, K., Attenborough, D. Frozen Planet: Episode 5 'Winter'. , BBC. (2011).
- Fothergill, A., Berlowitz, V., Attenborough, D. Ch. Winter: Life closes down. in Frozen Planet: A World Beyond Imagination. , BBC books. (2011).
- Yun, X., et al. Ice formation in the subcooled brine environment. Int. J. Heat Mass Transfer. 95, 198-205 (2016).
- Weast, R. C. CRC Handbook of Chemistry and Physics. 64, 64 edn, CRC Press. 257-258 (1983).
- Bejan, A., Lage, J. L. The Prandtl Number Effect on the Transition in Natural Convection Along a Vertical Surface. J. Heat Transfer. 112, 787-790 (1990).
