Computation of Ice Shedding Trajectories Using Cartesian Grids, Penalization, and Level Sets

<table>Schematic representation of a computational domain <svg height="23.025" id="M4" style="vertical-align:-5.76988pt" version="1.1" viewbox="0 0 117.6 23.025" width="117.6" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,15.775)"><path d="M719 146l-18 -146h-243l-5 18q2 4 17 25t24 36l25 42t26.5 52.5t21 56.5t16 65t5.5 68q0 121 -55.5 194.5t-149.5 73.5q-95 0 -150 -72t-55 -199q0 -54 18.5 -114t36.5 -94t48.5 -82.5t32.5 -51.5l-7 -18h-244l-17 146l27 1q17 -49 31.5 -62t45.5 -13h89l-21 21t-26 26.5
l-26 29t-28 35.5l-23 37.5t-20.5 43t-12 45.5t-5.5 51q0 133 88 218.5t219 85.5q132 0 218 -83t86 -218q0 -40 -11.5 -78t-33 -71.5t-39 -56t-44 -51t-32.5 -35.5h84q35 0 48 12t31 63z" id="x3A9"></path></g><g transform="matrix(.017,-0,0,-.017,17.639,15.775)"><path d="M535 323h-483v50h483v-50zM535 138h-483v50h483v-50z" id="x3D"></path></g><g transform="matrix(.017,-0,0,-.017,32.343,15.775)"><path d="M719 146l-18 -146h-243l-5 18q2 4 17 25t24 36l25 42t26.5 52.5t21 56.5t16 65t5.5 68q0 121 -55.5 194.5t-149.5 73.5q-95 0 -150 -72t-55 -199q0 -54 18.5 -114t36.5 -94t48.5 -82.5t32.5 -51.5l-7 -18h-244l-17 146l27 1q17 -49 31.5 -62t45.5 -13h89l-21 21t-26 26.5
l-26 29t-28 35.5l-23 37.5t-20.5 43t-12 45.5t-5.5 51q0 133 88 218.5t219 85.5q132 0 218 -83t86 -218q0 -40 -11.5 -78t-33 -71.5t-39 -56t-44 -51t-32.5 -35.5h84q35 0 48 12t31 63z" id="x3A9"></path></g>
<g transform="matrix(.012,-0,0,-.012,45.225,19.862)"><path d="M619 670q0 -13 -9 -26t-18 -19q-13 -10 -25 2q-36 38 -66 38q-31 0 -54.5 -50t-45.5 -185h120l-20 -31l-107 -12q-23 -138 -57 -293q-27 -122 -55 -184.5t-75 -109.5q-60 -61 -114 -61q-25 0 -47.5 15t-22.5 31q0 17 31 44q11 8 20 -1q10 -11 31 -19t35 -8q26 0 47 19
q34 34 71 253l54 315h-90l-3 12l31 30h70q28 138 90 204q35 37 75 57.5t70 20.5q26 0 45 -14t19 -28z" id="x1D453"></path></g>
<g transform="matrix(.017,-0,0,-.017,57.262,15.775)"><path d="M535 194q0 -90 -71.5 -149.5t-169.5 -59.5t-169.5 59.5t-71.5 149.5v354h58v-354q0 -68 54.5 -113.5t128.5 -45.5t128.5 45.5t54.5 113.5v354h58v-354z" id="x222A"></path></g><g transform="matrix(.017,-0,0,-.017,71.014,15.775)"><path d="M457 488l-30 -3q-17 148 -131 148q-53 0 -84.5 -34.5t-31.5 -82.5q0 -42 25.5 -72t74.5 -62l33 -22q63 -42 95 -85t32 -102q0 -84 -67 -137t-163 -53q-58 0 -113 22t-70 43l-4 152l27 4q4 -32 15 -62.5t31 -59.5t53.5 -47t76.5 -18q56 0 92 35t36 96q0 39 -25 70t-78 68
l-31 22q-32 23 -53.5 41.5t-45 57t-23.5 77.5q0 82 58 132.5t156 50.5q46 0 101 -17l18.5 -6t17 -6t8.5 -3q-4 -55 0 -147z" id="x1D446"></path></g>
<g transform="matrix(.012,-0,0,-.012,79.175,7.613)"><path d="M384 0h-275v27q67 5 81.5 18.5t14.5 68.5v385q0 38 -7.5 47.5t-40.5 10.5l-48 2v24q85 15 178 52v-521q0 -55 14.5 -68.5t82.5 -18.5v-27z" id="x31"></path></g>
<g transform="matrix(.017,-0,0,-.017,89.288,15.775)"><path d="M535 194q0 -90 -71.5 -149.5t-169.5 -59.5t-169.5 59.5t-71.5 149.5v354h58v-354q0 -68 54.5 -113.5t128.5 -45.5t128.5 45.5t54.5 113.5v354h58v-354z" id="x222A"></path></g><g transform="matrix(.017,-0,0,-.017,103.039,15.775)"><path d="M457 488l-30 -3q-17 148 -131 148q-53 0 -84.5 -34.5t-31.5 -82.5q0 -42 25.5 -72t74.5 -62l33 -22q63 -42 95 -85t32 -102q0 -84 -67 -137t-163 -53q-58 0 -113 22t-70 43l-4 152l27 4q4 -32 15 -62.5t31 -59.5t53.5 -47t76.5 -18q56 0 92 35t36 96q0 39 -25 70t-78 68
l-31 22q-32 23 -53.5 41.5t-45 57t-23.5 77.5q0 82 58 132.5t156 50.5q46 0 101 -17l18.5 -6t17 -6t8.5 -3q-4 -55 0 -147z" id="x1D446"></path></g>
<g transform="matrix(.012,-0,0,-.012,111.2,7.613)"><path d="M412 140l28 -9q0 -2 -35 -131h-373v23q112 112 161 170q59 70 92 127t33 115q0 63 -31 98t-86 35q-75 0 -137 -93l-22 20l57 81q55 59 135 59q69 0 118.5 -46.5t49.5 -122.5q0 -62 -29.5 -114t-102.5 -130l-141 -149h186q42 0 58.5 10.5t38.5 56.5z" id="x32"></path></g>
</svg>. <svg height="18.612499" id="M5" style="vertical-align:-5.76988pt" version="1.1" viewbox="0 0 21.25 18.612499" width="21.25" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,11.363)"><path d="M719 146l-18 -146h-243l-5 18q2 4 17 25t24 36l25 42t26.5 52.5t21 56.5t16 65t5.5 68q0 121 -55.5 194.5t-149.5 73.5q-95 0 -150 -72t-55 -199q0 -54 18.5 -114t36.5 -94t48.5 -82.5t32.5 -51.5l-7 -18h-244l-17 146l27 1q17 -49 31.5 -62t45.5 -13h89l-21 21t-26 26.5
l-26 29t-28 35.5l-23 37.5t-20.5 43t-12 45.5t-5.5 51q0 133 88 218.5t219 85.5q132 0 218 -83t86 -218q0 -40 -11.5 -78t-33 -71.5t-39 -56t-44 -51t-32.5 -35.5h84q35 0 48 12t31 63z" id="x3A9"></path></g>
<g transform="matrix(.012,-0,0,-.012,12.925,15.45)"><path d="M619 670q0 -13 -9 -26t-18 -19q-13 -10 -25 2q-36 38 -66 38q-31 0 -54.5 -50t-45.5 -185h120l-20 -31l-107 -12q-23 -138 -57 -293q-27 -122 -55 -184.5t-75 -109.5q-60 -61 -114 -61q-25 0 -47.5 15t-22.5 31q0 17 31 44q11 8 20 -1q10 -11 31 -19t35 -8q26 0 47 19
q34 34 71 253l54 315h-90l-3 12l31 30h70q28 138 90 204q35 37 75 57.5t70 20.5q26 0 45 -14t19 -28z" id="x1D453"></path></g>
</svg>: fluid domain, <svg height="16.362499" id="M6" style="vertical-align:-0.23206pt" version="1.1" viewbox="0 0 12.15 16.362499" width="12.15" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,16.012)"><path d="M457 488l-30 -3q-17 148 -131 148q-53 0 -84.5 -34.5t-31.5 -82.5q0 -42 25.5 -72t74.5 -62l33 -22q63 -42 95 -85t32 -102q0 -84 -67 -137t-163 -53q-58 0 -113 22t-70 43l-4 152l27 4q4 -32 15 -62.5t31 -59.5t53.5 -47t76.5 -18q56 0 92 35t36 96q0 39 -25 70t-78 68
l-31 22q-32 23 -53.5 41.5t-45 57t-23.5 77.5q0 82 58 132.5t156 50.5q46 0 101 -17l18.5 -6t17 -6t8.5 -3q-4 -55 0 -147z" id="x1D446"></path></g>
<g transform="matrix(.012,-0,0,-.012,8.225,7.85)"><path d="M244 607q0 -25 -15.5 -43t-37.5 -18q-19 0 -32 13t-13 35q0 21 15 41t39 20q20 0 32 -14t12 -34zM222 91q-29 -33 -79 -68t-75 -35q-13 0 -19 7.5t-6 31t10 65.5l62 253q5 26 -1 26q-21 0 -72 -43l-13 24q43 40 91 68t71 28q30 0 10 -78l-71 -274q-8 -30 3 -30
q16 0 76 48z" id="x1D456"></path></g>
</svg>: solid domains.</table>

Modelling and Simulation in Engineering

fig1

Figure 1

Figure 1: Computation of Ice Shedding Trajectories Using Cartesian Grids, Penalization, and Level Sets 