A Comparative Study on Evaluation Methods of Fluid Forces on Cartesian Grids

<table class="figure-group"><tr class="fig-image" id="a"><td><object data="https://static.hindawi.com/articles/mpe/volume-2017/8314615/figures/8314615.fig.0014a.svgz" name="8314615.fig.0014a" type="image/svg+xml"></object></td></tr><tr class="fig-caption"><td><b>(a) </b>Pressure drag</td></tr><tr class="fig-image" id="b"><td><object data="https://static.hindawi.com/articles/mpe/volume-2017/8314615/figures/8314615.fig.0014b.svgz" name="8314615.fig.0014b" type="image/svg+xml"></object></td></tr><tr class="fig-caption"><td><b>(b) </b>Friction drag</td></tr></table>

<div>Estimated averaged drag and oscillation due to alignment of body and grid in Stokes flow around a moving sphere as functions of mesh spacing. Analytical estimates of pressure and friction drag are <svg height="13.9365pt" id="M195" style="vertical-align:-5.3006pt" version="1.1" viewbox="-0.0498162 -8.6359 76.1559 13.9365" width="76.1559pt" 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="M645 631C614 643 545 666 457 666C215 666 23 519 23 283C23 90 158 -16 337 -16C412 -16 489 2 522 10C543 39 590 127 606 167L580 181C519 89 459 18 348 18C201 18 122 136 122 287C122 464 244 632 435 632C544 632 602 595 608 472L639 475C643 526 645 581 645 631Z" id="g113-68"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="668" vert-adv-y="668" xpluscharwidth="8.684"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,8.619,3.132)"><path d="M43 650V617C122 612 130 606 130 522V128C130 44 122 37 34 32V0H275C400 0 497 25 573 82C652 142 698 229 698 350C698 445 662 518 613 565C549 625 456 650 327 650H43ZM217 544C217 580 220 591 229 599C240 607 264 613 311 613C380 613 437 599 482 570C561 523 598 436 598 336C598 177 511 40 324 40C248 40 217 61 217 135V544Z" id="g58-66"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="744" vert-adv-y="744" xpluscharwidth="15.389377679999999"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,15.435,3.132)"><path d="M173 382V461C125 442 68 425 25 418V390C88 382 92 382 92 317V-127C92 -191 84 -199 18 -204V-235H259V-204C180 -196 173 -192 173 -120V7C188 -1 212 -11 243 -12C375 12 495 109 495 261C495 361 427 451 316 451C303 451 283 446 266 435L173 382ZM173 345C200 365 243 388 274 388C346 388 409 329 409 222C409 111 353 40 266 40C231 40 196 56 173 79V345Z" id="g58-110"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="534" vert-adv-y="534" xpluscharwidth="20.29438398"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,20.294,3.132)"><path d="M98 134C72 134 46 117 46 90C46 73 55 65 60 64C95 55 124 32 124 -4C124 -42 95 -68 44 -89L57 -123C122 -104 194 -60 194 22C194 94 136 134 98 134Z" id="g50-45"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="237" vert-adv-y="237" xpluscharwidth="22.45069289"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,22.451,3.132)"><path d="M384 109C349 76 310 59 268 59C197 59 123 113 120 245C240 248 367 259 384 262C403 266 407 275 407 296C407 376 338 451 254 451H253C201 451 148 425 105 378C65 334 39 272 39 203C39 88 112 -12 235 -12C266 -12 339 5 402 84L384 109ZM230 411C285 411 319 365 319 313C319 298 313 293 294 293C236 291 180 290 125 290C141 370 185 411 230 411Z" id="g58-99"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="434" vert-adv-y="434" xpluscharwidth="26.40038698"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,26.4,3.132)"><path d="M482 0V30C420 38 407 48 370 103L273 248C308 299 331 332 352 357C387 397 401 404 463 409V439H277V409C318 405 324 400 307 369C290 337 270 307 251 279L194 367C175 396 180 402 220 409V439H18V409C74 401 87 395 117 349L204 213C173 168 147 130 119 95C78 46 67 38 5 30V0H195V30C144 38 143 47 162 79C181 112 202 149 225 181L297 70C314 43 307 35 264 30V0H482Z" id="g58-118"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="481" vert-adv-y="481" xpluscharwidth="30.777085569999997"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,30.787,3.132)"><path d="M442 39L431 68C420 63 406 57 394 57C377 57 359 71 359 118V301C359 353 349 393 314 424C291 443 261 451 230 451C171 438 107 401 79 381C59 366 46 355 46 340C46 316 72 295 91 295C106 295 117 303 122 319C131 350 139 371 152 384C163 397 178 404 197 404C246 404 281 365 281 293V278C260 262 186 232 124 213C69 194 41 161 41 112C41 47 92 -12 163 -12C195 -12 243 25 283 51C288 35 293 22 308 7C319 -4 339 -12 354 -12L442 39ZM281 87C262 68 226 51 199 51C169 51 129 74 129 126C129 162 151 181 190 198C211 208 259 229 281 240V87Z" id="g58-95"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="449" vert-adv-y="449" xpluscharwidth="34.87288653"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,34.845,3.132)"><path d="M396 114C349 71 317 59 274 59C216 59 123 104 123 241C123 342 180 400 246 400C283 400 317 386 348 359C355 352 361 348 366 348C378 348 401 371 402 394C402 405 398 413 385 424C369 437 337 451 294 451H293C254 451 193 434 141 394C73 342 39 277 39 198C39 90 114 -12 243 -12C303 -12 371 31 415 90L396 114Z" id="g58-97"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="432" vert-adv-y="432" xpluscharwidth="38.776187039999996"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,38.85,3.132)"><path d="M303 37L293 66C280 59 256 49 231 49S173 63 173 143V396H281C294 404 298 428 287 439H173V575L158 577L92 509V439H46L17 408L22 396H92V107C92 27 128 -12 192 -12C203 -12 218 -9 237 2L303 37Z" id="g58-114"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="312" vert-adv-y="312" xpluscharwidth="41.68919064"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,45.942,0)"><path d="M535 323V373H52V323H535ZM535 138V188H52V138H535Z" id="g117-34"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="587" vert-adv-y="587" xpluscharwidth="53.573"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,57.205,0)"><path d="M249 635C141 635 70 555 70 471C70 401 114 353 179 316C143 294 106 267 90 252C68 231 45 202 45 157C45 50 130 -12 237 -12C322 -12 435 52 435 169C435 256 372 304 303 343C349 374 375 398 383 407C401 429 411 458 411 487C411 569 344 635 249 635ZM238 603C285 603 337 567 337 482C337 422 310 385 276 358C205 393 145 426 145 500C145 552 179 603 238 603ZM248 20C183 20 125 70 125 163C125 218 158 268 206 300C284 261 355 217 355 143C355 66 308 20 248 20Z" id="g113-57"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="480" vert-adv-y="480" xpluscharwidth="63.445"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,63.445,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="g113-49"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="480" vert-adv-y="480" xpluscharwidth="69.685"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,69.685,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="g113-49"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="480" vert-adv-y="480" xpluscharwidth="75.925"></glyph.data></g></svg> and <svg height="13.9365pt" id="M196" style="vertical-align:-5.3006pt" version="1.1" viewbox="-0.0498162 -8.6359 82.4193 13.9365" width="82.4193pt" 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="M645 631C614 643 545 666 457 666C215 666 23 519 23 283C23 90 158 -16 337 -16C412 -16 489 2 522 10C543 39 590 127 606 167L580 181C519 89 459 18 348 18C201 18 122 136 122 287C122 464 244 632 435 632C544 632 602 595 608 472L639 475C643 526 645 581 645 631Z" id="g113-68"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="668" vert-adv-y="668" xpluscharwidth="8.684"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,8.619,3.132)"><path d="M43 650V617C122 612 130 606 130 522V128C130 44 122 37 34 32V0H275C400 0 497 25 573 82C652 142 698 229 698 350C698 445 662 518 613 565C549 625 456 650 327 650H43ZM217 544C217 580 220 591 229 599C240 607 264 613 311 613C380 613 437 599 482 570C561 523 598 436 598 336C598 177 511 40 324 40C248 40 217 61 217 135V544Z" id="g58-66"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="744" vert-adv-y="744" xpluscharwidth="15.389377679999999"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,15.435,3.132)"><path d="M173 382V461C125 442 68 425 25 418V390C88 382 92 382 92 317V-127C92 -191 84 -199 18 -204V-235H259V-204C180 -196 173 -192 173 -120V7C188 -1 212 -11 243 -12C375 12 495 109 495 261C495 361 427 451 316 451C303 451 283 446 266 435L173 382ZM173 345C200 365 243 388 274 388C346 388 409 329 409 222C409 111 353 40 266 40C231 40 196 56 173 79V345Z" id="g58-110"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="534" vert-adv-y="534" xpluscharwidth="20.29438398"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,20.294,3.132)"><path d="M98 134C72 134 46 117 46 90C46 73 55 65 60 64C95 55 124 32 124 -4C124 -42 95 -68 44 -89L57 -123C122 -104 194 -60 194 22C194 94 136 134 98 134Z" id="g50-45"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="237" vert-adv-y="237" xpluscharwidth="22.45069289"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,22.451,3.132)"><path d="M384 109C349 76 310 59 268 59C197 59 123 113 120 245C240 248 367 259 384 262C403 266 407 275 407 296C407 376 338 451 254 451H253C201 451 148 425 105 378C65 334 39 272 39 203C39 88 112 -12 235 -12C266 -12 339 5 402 84L384 109ZM230 411C285 411 319 365 319 313C319 298 313 293 294 293C236 291 180 290 125 290C141 370 185 411 230 411Z" id="g58-99"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="434" vert-adv-y="434" xpluscharwidth="26.40038698"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,26.4,3.132)"><path d="M482 0V30C420 38 407 48 370 103L273 248C308 299 331 332 352 357C387 397 401 404 463 409V439H277V409C318 405 324 400 307 369C290 337 270 307 251 279L194 367C175 396 180 402 220 409V439H18V409C74 401 87 395 117 349L204 213C173 168 147 130 119 95C78 46 67 38 5 30V0H195V30C144 38 143 47 162 79C181 112 202 149 225 181L297 70C314 43 307 35 264 30V0H482Z" id="g58-118"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="481" vert-adv-y="481" xpluscharwidth="30.777085569999997"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,30.787,3.132)"><path d="M442 39L431 68C420 63 406 57 394 57C377 57 359 71 359 118V301C359 353 349 393 314 424C291 443 261 451 230 451C171 438 107 401 79 381C59 366 46 355 46 340C46 316 72 295 91 295C106 295 117 303 122 319C131 350 139 371 152 384C163 397 178 404 197 404C246 404 281 365 281 293V278C260 262 186 232 124 213C69 194 41 161 41 112C41 47 92 -12 163 -12C195 -12 243 25 283 51C288 35 293 22 308 7C319 -4 339 -12 354 -12L442 39ZM281 87C262 68 226 51 199 51C169 51 129 74 129 126C129 162 151 181 190 198C211 208 259 229 281 240V87Z" id="g58-95"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="449" vert-adv-y="449" xpluscharwidth="34.87288653"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,34.845,3.132)"><path d="M396 114C349 71 317 59 274 59C216 59 123 104 123 241C123 342 180 400 246 400C283 400 317 386 348 359C355 352 361 348 366 348C378 348 401 371 402 394C402 405 398 413 385 424C369 437 337 451 294 451H293C254 451 193 434 141 394C73 342 39 277 39 198C39 90 114 -12 243 -12C303 -12 371 31 415 90L396 114Z" id="g58-97"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="432" vert-adv-y="432" xpluscharwidth="38.776187039999996"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,38.85,3.132)"><path d="M303 37L293 66C280 59 256 49 231 49S173 63 173 143V396H281C294 404 298 428 287 439H173V575L158 577L92 509V439H46L17 408L22 396H92V107C92 27 128 -12 192 -12C203 -12 218 -9 237 2L303 37Z" id="g58-114"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="312" vert-adv-y="312" xpluscharwidth="41.68919064"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,45.942,0)"><path d="M535 323V373H52V323H535ZM535 138V188H52V138H535Z" id="g117-34"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="587" vert-adv-y="587" xpluscharwidth="53.573"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,57.205,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" xpluscharwidth="63.445"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,63.445,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="g113-55"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="480" vert-adv-y="480" xpluscharwidth="69.685"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,69.685,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="g113-49"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="480" vert-adv-y="480" xpluscharwidth="75.925"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,75.925,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="g113-49"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="480" vert-adv-y="480" xpluscharwidth="82.16499999999999"></glyph.data></g></svg>, respectively.</div>

Mathematical Problems in Engineering

fig14

Figure 14

Figure 14: A Comparative Study on Evaluation Methods of Fluid Forces on Cartesian Grids 