Numerical Simulation of the Generalized Newtonian Free Surface Flows by a Density Reinitialization SPH Method

<div>The shape of a Newtonian droplet at impact angles of <svg height="11.0423pt" id="M242" style="vertical-align:-0.2456007pt" version="1.1" viewbox="-0.0498162 -10.7967 16.6412 11.0423" width="16.6412pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M384 0V27C293 34 287 42 287 114V635C232 613 172 594 109 583V559L157 557C201 555 205 550 205 499V114C205 42 199 34 109 27V0H384Z" id="g121-47"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="480" vert-adv-y="480"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,6.24,0)"><path d="M153 550H386L412 615L406 623H120L82 318C104 327 142 338 184 338C294 338 347 275 347 187C347 112 305 39 221 39C160 39 119 71 97 89C88 97 80 96 71 90C59 80 50 67 49 57C48 45 52 36 66 23C80 9 123 -12 169 -12C221 -11 288 15 342 59C403 109 431 165 431 225C431 308 366 395 238 395C212 395 165 379 127 364L153 550Z" id="g121-51"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="480" vert-adv-y="480"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,12.48,-5.741)"><path d="M339 262C339 340 273 403 193 403S48 340 48 262S113 121 193 121S339 184 339 262ZM283 262C283 210 242 168 193 168S104 210 104 262S144 356 193 356S283 314 283 262Z" id="g54-73"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="385" vert-adv-y="385"></glyph.data></g></svg> (first column), <svg height="11.0423pt" id="M243" style="vertical-align:-0.2456007pt" version="1.1" viewbox="-0.0498162 -10.7967 16.6412 11.0423" width="16.6412pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M285 378C315 398 338 416 353 432C373 451 384 474 384 503C384 579 325 635 236 635H235C182 635 136 610 108 579L65 516L85 496C110 533 150 575 205 575C258 575 300 543 300 481C300 407 232 369 141 339L147 310C163 315 188 321 211 321C268 321 338 284 338 192C338 94 288 40 217 40C160 40 119 68 93 91C85 98 77 97 69 91C60 84 47 71 46 58C44 46 48 35 62 22C75 10 116 -12 162 -12C234 -12 424 62 424 224C424 297 373 359 285 376V378Z" id="g121-49"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="480" vert-adv-y="480"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,6.24,0)"><path d="M241 635C89 635 35 457 35 312C35 153 89 -12 240 -12C390 -12 443 166 443 312C443 466 390 635 241 635ZM238 602C329 602 354 454 354 312C354 172 330 22 240 22C152 22 124 173 124 313S148 602 238 602Z" id="g121-46"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="480" vert-adv-y="480"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,12.48,-5.741)"><path d="M339 262C339 340 273 403 193 403S48 340 48 262S113 121 193 121S339 184 339 262ZM283 262C283 210 242 168 193 168S104 210 104 262S144 356 193 356S283 314 283 262Z" id="g54-73"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="385" vert-adv-y="385"></glyph.data></g></svg> (second column), and <svg height="11.0423pt" id="M244" style="vertical-align:-0.2456007pt" version="1.1" viewbox="-0.0498162 -10.7967 16.6412 11.0423" width="16.6412pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M456 178V225H360V632H320C217 496 115 347 20 206V178H280V106C280 40 276 34 189 27V0H445V27C364 34 360 39 360 106V178H456ZM280 225H82C149 335 214 431 278 520H280V225Z" id="g121-50"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="480" vert-adv-y="480"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,6.24,0)"><path d="M153 550H386L412 615L406 623H120L82 318C104 327 142 338 184 338C294 338 347 275 347 187C347 112 305 39 221 39C160 39 119 71 97 89C88 97 80 96 71 90C59 80 50 67 49 57C48 45 52 36 66 23C80 9 123 -12 169 -12C221 -11 288 15 342 59C403 109 431 165 431 225C431 308 366 395 238 395C212 395 165 379 127 364L153 550Z" id="g121-51"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="480" vert-adv-y="480"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,12.48,-5.741)"><path d="M339 262C339 340 273 403 193 403S48 340 48 262S113 121 193 121S339 184 339 262ZM283 262C283 210 242 168 193 168S104 210 104 262S144 356 193 356S283 314 283 262Z" id="g54-73"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="385" vert-adv-y="385"></glyph.data></g></svg> (third column) obtained using PDRI-SPH, varying with short dimensionless time.</div>

Mathematical Problems in Engineering

fig12

Figure 12

Figure 12: Numerical Simulation of the Generalized Newtonian Free Surface Flows by a Density Reinitialization SPH Method 