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

<table class="table-group" width="800px"><tr class="fig-image" id="a"><td><img alt="274947.fig.008a" src="https://static.hindawi.com/articles/mse/volume-2011/274947/figures/274947.fig.008a.jpg"/></td></tr><tr class="fig-caption"><td align="center"><table><tr><td><b>(a) </b>Initial solution</td></tr></table></td></tr><tr class="fig-image" id="b"><td><img alt="274947.fig.008b" src="https://static.hindawi.com/articles/mse/volume-2011/274947/figures/274947.fig.008b.jpg"/></td></tr><tr class="fig-caption"><td align="center"><table><tr><td><b>(b) </b>Solution after 2800 time steps</td></tr></table></td></tr><tr class="fig-image" id="c"><td><img alt="274947.fig.008c" src="https://static.hindawi.com/articles/mse/volume-2011/274947/figures/274947.fig.008c.jpg"/></td></tr><tr class="fig-caption"><td align="center"><table><tr><td><b>(c) </b>Solution after 5600 time steps</td></tr></table></td></tr><tr class="fig-image" id="d"><td><img alt="274947.fig.008d" src="https://static.hindawi.com/articles/mse/volume-2011/274947/figures/274947.fig.008d.jpg"/></td></tr><tr class="fig-caption"><td align="center"><table><tr><td><b>(d) </b>Solution after 8400 time steps</td></tr></table></td></tr></table>

<table>Cylinder at <svg height="11.375" id="M154" style="vertical-align:-0.1638pt" version="1.1" viewbox="0 0 64.637497 11.375" width="64.637497" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,11.112)"><path d="M631 18l-4 -26q-24 0 -35 1q-57 3 -92.5 23t-69.5 68q-21 30 -101 160q-13 22 -30.5 31t-54.5 9h-32v-159q0 -63 14 -77.5t74 -19.5v-28h-262v28q62 5 76 19.5t14 77.5v402q0 62 -13 76t-73 19v28h248q109 0 159 -34q67 -42 67 -131q0 -115 -127 -172q27 -52 89 -144
q49 -71 81 -104q31 -36 72 -47zM212 316h48q73 0 108 29q55 42 55 124q0 76 -42 112t-104 36q-41 0 -54 -10q-11 -7 -11 -44v-247z" id="x52"></path></g><g transform="matrix(.017,-0,0,-.017,10.619,11.112)"><path d="M380 106l15 -22q-42 -52 -87 -74t-76 -22q-89 0 -142 62t-53 151q0 103 66 175q67 73 146 73h1q61 0 105.5 -46.5t44.5 -105.5q0 -16 -4.5 -22.5t-18.5 -9.5q-54 -10 -262 -17q1 -95 46 -143.5t104 -48.5q61 0 115 50zM225 412q-35 0 -64.5 -32.5t-40.5 -90.5q83 0 170 3
q14 0 19 4.5t5 15.5q1 41 -23.5 70.5t-65.5 29.5z" id="x65"></path></g><g transform="matrix(.017,-0,0,-.017,25.391,11.112)"><path d="M535 323h-483v50h483v-50zM535 138h-483v50h483v-50z" id="x3D"></path></g><g transform="matrix(.017,-0,0,-.017,40.095,11.112)"><path d="M153 550l-26 -186q79 31 111 31q90 0 141.5 -51t51.5 -119q0 -93 -89 -166q-85 -69 -173 -71q-32 0 -61.5 11.5t-41.5 23.5q-18 17 -17 34q2 16 22 33q14 9 26 -1q61 -50 124 -50q60 0 93 43.5t33 104.5q0 69 -41.5 110t-121.5 41q-53 0 -102 -20l38 305h286l6 -8
l-26 -65h-233z" id="x35"></path></g><g transform="matrix(.017,-0,0,-.017,48.254,11.112)"><path d="M153 550l-26 -186q79 31 111 31q90 0 141.5 -51t51.5 -119q0 -93 -89 -166q-85 -69 -173 -71q-32 0 -61.5 11.5t-41.5 23.5q-18 17 -17 34q2 16 22 33q14 9 26 -1q61 -50 124 -50q60 0 93 43.5t33 104.5q0 69 -41.5 110t-121.5 41q-53 0 -102 -20l38 305h286l6 -8
l-26 -65h-233z" id="x35"></path></g><g transform="matrix(.017,-0,0,-.017,56.413,11.112)"><path d="M241 635q53 0 94 -28.5t63.5 -76t33.5 -102.5t11 -116q0 -58 -11 -112.5t-34 -103.5t-63.5 -78.5t-94.5 -29.5t-95 28t-64.5 75t-34.5 102.5t-11 118.5q0 58 11.5 112.5t34.5 103t64.5 78t95.5 29.5zM238 602q-32 0 -55.5 -25t-35.5 -68t-17.5 -91t-5.5 -105
q0 -76 10 -138.5t37 -107.5t69 -45q32 0 55.5 25t35.5 68.5t17.5 91.5t5.5 105t-5.5 105.5t-18 92t-36 68t-56.5 24.5z" id="x30"></path></g>
</svg>. Vorticity contours corresponding to shed 13: green positive vorticity and blue negative vorticity.</table>

Modelling and Simulation in Engineering

fig8

Figure 8

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