Conditional Stochastic Simulations of Flow and  Transport with Karhunen-Loève Expansions, Stochastic  Collocation, and Sequential Gaussian Simulation

<div>Transport results with SGS-MC (a) and KL-MC (b) with <svg height="9.21094pt" id="M501" style="vertical-align:-0.2129908pt" version="1.1" viewbox="-0.0498162 -8.99795 55.1573 9.21094" width="55.1573pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.0135,0,0,-0.0135,0,0)"><path d="M962 650H795L470 145L347 650H176L170 622C268 613 275 606 244 503L190 322C150 188 132 126 118 91C102 50 80 33 18 28L12 0H245L251 28C175 35 170 48 174 93C177 128 191 180 220 284L292 542H294C331 392 383 150 409 4H432L774 555H776L714 137C700 40 694 34 612 28L606 0H868L874 28C793 34 784 37 797 137L849 533C859 612 863 616 956 622L962 650Z" id="g113-78"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="985" vert-adv-y="985"></glyph.data></g><g transform="matrix(.0135,0,0,-0.0135,17.132,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(.0135,0,0,-0.0135,28.87,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><g transform="matrix(.0135,0,0,-0.0135,35.373,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"></glyph.data></g><g transform="matrix(.0135,0,0,-0.0135,41.877,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"></glyph.data></g><g transform="matrix(.0135,0,0,-0.0135,48.38,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"></glyph.data></g></svg>. From left to right: moments of breakthrough curves <svg height="9.36053pt" id="M502" style="vertical-align:-0.2401896pt" version="1.1" viewbox="-0.0498162 -9.12034 25.2849 9.36053" width="25.2849pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.0135,0,0,-0.0135,0,0)"><path d="M372 352C427 368 507 407 507 502C507 548 486 592 451 616C420 638 376 650 301 650H44V622C124 616 128 610 128 527V123C128 40 122 34 36 28V0H257C336 0 402 12 453 40C512 72 548 124 548 191C548 291 471 335 372 350V352ZM213 362V558C213 591 216 602 226 608C235 614 258 618 282 618C377 618 415 558 415 490C415 406 369 362 260 362H213ZM213 328H258C380 328 450 283 450 181C450 76 377 35 296 34C231 34 213 53 213 125V328Z" id="g190-67"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="588" vert-adv-y="588"></glyph.data></g><g transform="matrix(.0135,0,0,-0.0135,7.83,0)"><path d="M592 498C586 556 582 629 581 675H560C545 656 537 650 510 650H115C87 650 74 652 60 675H40C38 620 33 555 28 495H57C69 539 77 568 89 585C102 607 122 616 205 616H267V123C267 44 259 34 168 28V0H454V28C360 34 352 44 352 123V616H422C498 616 514 608 531 584C542 568 552 543 562 495L592 498Z" id="g190-85"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="617" vert-adv-y="617"></glyph.data></g><g transform="matrix(.0135,0,0,-0.0135,16.081,0)"><path d="M614 175C564 76 510 21 408 21C256 21 146 149 146 336C146 488 235 629 402 629C510 629 570 586 597 480L626 488C620 541 614 582 606 638C578 643 510 665 429 665C206 665 44 527 44 316C44 157 153 -15 402 -15C474 -15 558 5 586 11C604 45 629 119 643 165L614 175Z" id="g190-68"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="665" vert-adv-y="665"></glyph.data></g></svg>  <svg height="12.0395pt" id="M503" style="vertical-align:-2.3616pt" version="1.1" viewbox="-0.0498162 -9.6779 31.3637 12.0395" width="31.3637pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.0135,0,0,-0.0135,0,0)"><path d="M300 -147C201 -63 143 98 143 270S200 602 300 686L282 710C136 610 70 450 70 271V270C70 89 136 -72 282 -170L300 -147Z" id="g113-41"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="346" vert-adv-y="346"></glyph.data></g><g transform="matrix(.0135,0,0,-0.0135,4.688,0)"><path d="M46 650V622C120 617 128 613 128 525V125C128 42 120 34 40 28V0H311V28C221 34 212 39 212 124V281L286 262C297 261 316 261 331 263C429 275 526 338 526 468C526 533 501 579 462 609C422 638 364 650 293 650H46ZM212 559C212 588 215 600 223 606C230 613 251 618 279 618C361 618 430 572 430 464C430 337 350 302 285 302C252 302 225 309 212 314V559Z" id="g190-81"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="563" vert-adv-y="563"></glyph.data></g><g transform="matrix(.0135,0,0,-0.0135,12.315,0)"><path d="M687 650H462V622C543 612 549 605 530 547C498 447 422 252 372 126H370C302 298 229 492 204 563C188 607 191 615 262 622V650H17V622C77 616 93 608 122 534C180 389 262 172 329 -11H360C436 196 541 450 568 516C606 605 619 614 687 622V650Z" id="g190-87"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="703" vert-adv-y="703"></glyph.data></g><g transform="matrix(.0135,0,0,-0.0135,21.84,0)"><path d="M303 0V28C221 34 213 39 213 125V525C213 610 221 616 303 622V650H38V622C120 616 128 610 128 525V125C128 40 120 34 38 28V0H303Z" id="g190-74"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="341" vert-adv-y="341"></glyph.data></g><g transform="matrix(.0135,0,0,-0.0135,26.541,0)"><path d="M275 270C275 450 212 609 64 710L45 686C145 604 203 442 203 270S147 -63 45 -147L64 -170C213 -68 275 89 275 270Z" id="g113-42"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="346" vert-adv-y="346"></glyph.data></g></svg>, and histograms of breakthrough times <svg height="9.17016pt" id="M504" style="vertical-align:-0.04981995pt" version="1.1" viewbox="-0.0498162 -9.12034 25.1081 9.17016" width="25.1081pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.0135,0,0,-0.0135,0,0)"><path d="M372 352C427 368 507 407 507 502C507 548 486 592 451 616C420 638 376 650 301 650H44V622C124 616 128 610 128 527V123C128 40 122 34 36 28V0H257C336 0 402 12 453 40C512 72 548 124 548 191C548 291 471 335 372 350V352ZM213 362V558C213 591 216 602 226 608C235 614 258 618 282 618C377 618 415 558 415 490C415 406 369 362 260 362H213ZM213 328H258C380 328 450 283 450 181C450 76 377 35 296 34C231 34 213 53 213 125V328Z" id="g190-67"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="588" vert-adv-y="588"></glyph.data></g><g transform="matrix(.0135,0,0,-0.0135,7.83,0)"><path d="M592 498C586 556 582 629 581 675H560C545 656 537 650 510 650H115C87 650 74 652 60 675H40C38 620 33 555 28 495H57C69 539 77 568 89 585C102 607 122 616 205 616H267V123C267 44 259 34 168 28V0H454V28C360 34 352 44 352 123V616H422C498 616 514 608 531 584C542 568 552 543 562 495L592 498Z" id="g190-85"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="617" vert-adv-y="617"></glyph.data></g><g transform="matrix(.0135,0,0,-0.0135,16.557,0)"><path d="M592 498C586 556 582 629 581 675H560C545 656 537 650 510 650H115C87 650 74 652 60 675H40C38 620 33 555 28 495H57C69 539 77 568 89 585C102 607 122 616 205 616H267V123C267 44 259 34 168 28V0H454V28C360 34 352 44 352 123V616H422C498 616 514 608 531 584C542 568 552 543 562 495L592 498Z" id="g190-85"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="617" vert-adv-y="617"></glyph.data></g></svg>.</div>

Journal of Applied Mathematics

fig8

Figure 8

Figure 8: Conditional Stochastic Simulations of Flow and  Transport with Karhunen-Loève Expansions, Stochastic  Collocation, and Sequential Gaussian Simulation 

Figure 8 | Conditional Stochastic Simulations of Flow and  Transport with Karhunen-Loève Expansions, Stochastic  Collocation, and Sequential Gaussian Simulation 