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

<div>The shape of a Newtonian droplet impact on horizontal (left) or inclined (<svg height="11.0423pt" id="M210" style="vertical-align:-0.2456007pt" version="1.1" viewbox="-0.0498162 -10.7967 38.9762 11.0423" width="38.9762pt" 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="M545 106L524 126C493 85 467 65 455 65C438 65 427 113 405 238C448 295 498 362 543 439L533 448L478 435C453 386 423 331 398 295H395C370 404 347 448 282 448C169 448 23 309 23 153C23 54 65 -12 128 -12C203 -12 283 70 339 155H341C360 29 380 -12 411 -12C444 -12 491 11 545 106ZM333 204C265 95 210 54 169 54C137 54 113 96 113 171C113 302 191 405 252 405C301 405 318 306 333 204Z" id="g113-223"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="568" vert-adv-y="568"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,11.017,0)"><path d="M535 323V373H52V323H535ZM535 138V188H52V138H535Z" id="g117-34"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="587" vert-adv-y="587"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,22.28,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,28.52,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,34.76,-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>, right) rigid plate obtained by PDRI-SPH method with <svg height="11.5303pt" id="M211" style="vertical-align:-3.18148pt" version="1.1" viewbox="-0.0498162 -8.34882 35.2273 11.5303" width="35.2273pt" 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="M545 106L524 126C493 85 467 65 455 65C438 65 427 113 405 238C448 295 498 362 543 439L533 448L478 435C453 386 423 331 398 295H395C370 404 347 448 282 448C169 448 23 309 23 153C23 54 65 -12 128 -12C203 -12 283 70 339 155H341C360 29 380 -12 411 -12C444 -12 491 11 545 106ZM333 204C265 95 210 54 169 54C137 54 113 96 113 171C113 302 191 405 252 405C301 405 318 306 333 204Z" id="g113-223"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="568" vert-adv-y="568"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,7.059,3.132)"><path d="M308 0V33C225 38 217 44 217 128V521C217 605 225 612 308 617V650H39V617C122 612 130 605 130 521V128C130 44 122 38 39 33V0H308Z" id="g58-71"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="348" vert-adv-y="348"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,10.226,3.132)"><path d="M308 0V33C225 38 217 44 217 128V521C217 605 225 612 308 617V650H39V617C122 612 130 605 130 521V128C130 44 122 38 39 33V0H308Z" id="g58-71"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="348" vert-adv-y="348"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,17.546,0)"><path d="M535 323V373H52V323H535ZM535 138V188H52V138H535Z" id="g117-34"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="587" vert-adv-y="587"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,28.809,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="g113-50"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="480" vert-adv-y="480"></glyph.data></g></svg>, <svg height="12.7178pt" id="M212" style="vertical-align:-3.42947pt" version="1.1" viewbox="-0.0498162 -9.28833 35.3969 12.7178" width="35.3969pt" 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="M558 587C558 666 497 712 432 712C379 712 330 691 284 650C212 586 178 508 148 348L71 -65C49 -185 35 -229 23 -235L27 -261C49 -259 101 -251 124 -224C131 -216 138 -178 159 -24C171 66 197 200 227 356C264 550 295 668 393 668C443 668 479 632 479 575C479 516 446 458 383 418C361 404 344 398 318 397L296 350C395 338 460 281 460 192C460 110 411 40 300 40C258 40 215 55 192 69L181 51C191 18 222 -16 266 -16C308 -16 351 1 397 26C471 67 545 142 545 221C545 315 486 365 401 395C469 437 558 498 558 587Z" id="g113-224"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="581" vert-adv-y="581"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,7.228,3.132)"><path d="M308 0V33C225 38 217 44 217 128V521C217 605 225 612 308 617V650H39V617C122 612 130 605 130 521V128C130 44 122 38 39 33V0H308Z" id="g58-71"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="348" vert-adv-y="348"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,10.395,3.132)"><path d="M308 0V33C225 38 217 44 217 128V521C217 605 225 612 308 617V650H39V617C122 612 130 605 130 521V128C130 44 122 38 39 33V0H308Z" id="g58-71"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="348" vert-adv-y="348"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,17.715,0)"><path d="M535 323V373H52V323H535ZM535 138V188H52V138H535Z" id="g117-34"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="587" vert-adv-y="587"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,28.978,0)"><path d="M412 140C382 77 369 73 315 73H129L270 222C362 320 402 379 402 466C402 571 322 635 234 635C177 635 130 609 99 576L42 495L64 475C90 514 133 568 201 568C274 568 318 519 318 435C318 349 255 267 193 193C144 135 87 78 32 23V0H405C417 45 427 89 440 131L412 140Z" id="g121-48"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="480" vert-adv-y="480"></glyph.data></g></svg> and different artificial stress parameter <svg height="9.49473pt" id="M213" style="vertical-align:-0.2063999pt" version="1.1" viewbox="-0.0498162 -9.28833 6.21949 9.49473" width="6.21949pt" 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="M452 333C452 398 425 448 375 448C329 448 235 404 140 292H138L229 683C233 701 234 712 225 712C208 712 149 686 66 677L64 652H97C135 652 140 651 129 598L38 167C23 96 29 57 44 33C60 7 98 -12 132 -12C169 -12 232 3 290 41C383 102 452 223 452 333ZM365 316C365 199 308 87 239 49C221 39 200 35 183 35C137 35 99 70 112 163C116 191 123 215 129 233C190 316 290 392 330 392C348 392 365 372 365 316Z" id="g113-99"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="469" vert-adv-y="469"></glyph.data></g></svg>, at dimensionless time <svg height="8.63352pt" id="M214" style="vertical-align:-0.2455397pt" version="1.1" viewbox="-0.0498162 -8.38798 34.9745 8.63352" width="34.9745pt" 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="M324 430H196L233 583L223 592L145 529L120 430H54L29 396L31 388H111L56 126C33 15 54 -12 77 -12C137 -12 214 57 250 95L233 119C208 92 155 59 138 59C126 59 120 70 131 125L186 390L298 394L324 430Z" id="g113-117"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="341" vert-adv-y="341"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,8.065,0)"><path d="M535 323V373H52V323H535ZM535 138V188H52V138H535Z" id="g117-34"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="587" vert-adv-y="587"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,19.328,0)"><path d="M137 343C167 482 260 545 321 574C357 591 397 603 429 609L423 641C382 634 335 622 295 608C189 570 37 457 37 238C37 84 125 -12 242 -12C362 -12 447 89 447 209C447 311 374 393 267 393C247 393 226 386 204 376L137 343ZM227 337C318 337 361 256 361 173C361 105 336 22 258 22C176 22 126 120 126 240C126 266 127 291 132 310C155 323 189 337 227 337Z" id="g121-52"></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,25.568,0)"><path d="M113 -12C146 -12 170 11 170 46C170 78 146 103 114 103S58 78 58 46C58 11 82 -12 113 -12Z" id="g113-47"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="228" vert-adv-y="228"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,28.532,0)"><path d="M412 140C382 77 369 73 315 73H129L270 222C362 320 402 379 402 466C402 571 322 635 234 635C177 635 130 609 99 576L42 495L64 475C90 514 133 568 201 568C274 568 318 519 318 435C318 349 255 267 193 193C144 135 87 78 32 23V0H405C417 45 427 89 440 131L412 140Z" id="g121-48"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="480" vert-adv-y="480"></glyph.data></g></svg>.</div>

Mathematical Problems in Engineering

fig7

Figure 7

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