Pizza3 - WORKSHOP - PostTreatment - Part 2

INRAE\Han Chen, Olivier Vitrac - rev. 2023-08-28, 2023-09-06

Synopsis

The second part of the worshop will demonstrate how to interpret simulated data in space by projecting them on stuctured or unstructured grids, and to extract stresses.
The visualization techniques and use of the Verlet list are more sophisticated than in part 1. Computations are more expensive and require understanding of vectorial algebra.
The approaches are illustrated on more realistic data.

Description and dependencies

%% Workshop Main File - Part 2
% This script demonstrates the data preprocessing, analysis, and visualization techniques
% using 'lamdumpread2', a custom function for reading LAMMPS dump files.
 
% File Structure (change your local path to reflect the content)
% ├── example2.m (main developing file)
% ├── notebook
% │   ├── example2.mlx <-- this file
% └── ...
% └── data folder (dumps/pub1/)
 
%INRAE\Olivier Vitrac, Han Chen (rev. 2023-08-30,2023-08-31)
 
% Revision history
% 2023-09-04,05 implement grid interpolation
% 2023-09-06 deploy as a notebook
 

1. Initialization

1.1 General definitions

%% Initialization and Preprocessing
% -------------------------------------------------
% Turn preprocessing on or off based on `PREPROCESS_FLAG`.
% 'lamdumpread2' splits large 3D dump files if the flag is true.
 
statvec = @(f,before,after) dispf('%s: %s%d values | average = %0.5g %s', before, ...
cell2mat(cellfun(@(x) sprintf(' %0.1f%%> %10.4g | ', x, prctile(f, x)), {2.5, 25, 50, 75, 97.5}, 'UniformOutput', false)), ...
length(f), mean(f),after); % user function to display statistics on vectors (usage: statvec(f,'myvar','ok'))
 
PREPROCESS_FLAG = false;
addpath(rootdir(pwd),'-begin')
if PREPROCESS_FLAG
datafolder = '../dumps/pub1/';
lamdumpread2(fullfile(datafolder,'dump.*'),'split'); % for large 3D
end
 
%% File-specific Processing (same file as in Example1.m but larger)
% -----------------------------------------------------------------
% A specific dump file is selected and pre-processed to extract data and
% save to default folders for further operations.
datafolder = '../dumps/pub1/';
dumpfile = 'dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle';
datafolder = lamdumpread2(fullfile(datafolder,dumpfile),'search'); % fix datafolder based on initial guess1.2
Look for 'dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle'... (be patient) ...found in '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary' The dumpfile 'dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle' has been found with the method 'splitfolder' in the folder: ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary the original search started in ../dumps/pub1 The frame (split) folder is: ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle The first frame (split) is located in: ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle/TIMESTEP_000000000.mat The prefetch is split in several files. 1014 TIMESTEPS are availble: Column 01 Column 02 Column 03 Column 04 Column 05 Column 06 Column 07 Column 08 Column 09 Column 10 0 102000 204000 306000 408000 510000 612000 714000 816000 918000 1000 103000 205000 307000 409000 511000 613000 715000 817000 919000 2000 104000 206000 308000 410000 512000 614000 716000 818000 920000 3000 105000 207000 309000 411000 513000 615000 717000 819000 921000 4000 106000 208000 310000 412000 514000 616000 718000 820000 922000 5000 107000 209000 311000 413000 515000 617000 719000 821000 923000 6000 108000 210000 312000 414000 516000 618000 720000 822000 924000 7000 109000 211000 313000 415000 517000 619000 721000 823000 925000 8000 110000 212000 314000 416000 518000 620000 722000 824000 926000 9000 111000 213000 315000 417000 519000 621000 723000 825000 927000 10000 112000 214000 316000 418000 520000 622000 724000 826000 928000 11000 113000 215000 317000 419000 521000 623000 725000 827000 929000 12000 114000 216000 318000 420000 522000 624000 726000 828000 930000 13000 115000 217000 319000 421000 523000 625000 727000 829000 931000 14000 116000 218000 320000 422000 524000 626000 728000 830000 932000 15000 117000 219000 321000 423000 525000 627000 729000 831000 933000 16000 118000 220000 322000 424000 526000 628000 730000 832000 934000 17000 119000 221000 323000 425000 527000 629000 731000 833000 935000 18000 120000 222000 324000 426000 528000 630000 732000 834000 936000 19000 121000 223000 325000 427000 529000 631000 733000 835000 937000 20000 122000 224000 326000 428000 530000 632000 734000 836000 938000 21000 123000 225000 327000 429000 531000 633000 735000 837000 939000 22000 124000 226000 328000 430000 532000 634000 736000 838000 940000 23000 125000 227000 329000 431000 533000 635000 737000 839000 941000 24000 126000 228000 330000 432000 534000 636000 738000 840000 942000 25000 127000 229000 331000 433000 535000 637000 739000 841000 943000 26000 128000 230000 332000 434000 536000 638000 740000 842000 944000 27000 129000 231000 333000 435000 537000 639000 741000 843000 945000 28000 130000 232000 334000 436000 538000 640000 742000 844000 946000 29000 131000 233000 335000 437000 539000 641000 743000 845000 947000 30000 132000 234000 336000 438000 540000 642000 744000 846000 948000 31000 133000 235000 337000 439000 541000 643000 745000 847000 949000 32000 134000 236000 338000 440000 542000 644000 746000 848000 950000 33000 135000 237000 339000 441000 543000 645000 747000 849000 951000 34000 136000 238000 340000 442000 544000 646000 748000 850000 952000 35000 137000 239000 341000 443000 545000 647000 749000 851000 953000 36000 138000 240000 342000 444000 546000 648000 750000 852000 954000 37000 139000 241000 343000 445000 547000 649000 751000 853000 955000 38000 140000 242000 344000 446000 548000 650000 752000 854000 956000 39000 141000 243000 345000 447000 549000 651000 753000 855000 957000 40000 142000 244000 346000 448000 550000 652000 754000 856000 958000 41000 143000 245000 347000 449000 551000 653000 755000 857000 959000 42000 144000 246000 348000 450000 552000 654000 756000 858000 960000 43000 145000 247000 349000 451000 553000 655000 757000 859000 961000 44000 146000 248000 350000 452000 554000 656000 758000 860000 962000 45000 147000 249000 351000 453000 555000 657000 759000 861000 963000 46000 148000 250000 352000 454000 556000 658000 760000 862000 964000 47000 149000 251000 353000 455000 557000 659000 761000 863000 965000 48000 150000 252000 354000 456000 558000 660000 762000 864000 966000 49000 151000 253000 355000 457000 559000 661000 763000 865000 967000 50000 152000 254000 356000 458000 560000 662000 764000 866000 968000 51000 153000 255000 357000 459000 561000 663000 765000 867000 969000 52000 154000 256000 358000 460000 562000 664000 766000 868000 970000 53000 155000 257000 359000 461000 563000 665000 767000 869000 971000 54000 156000 258000 360000 462000 564000 666000 768000 870000 972000 55000 157000 259000 361000 463000 565000 667000 769000 871000 973000 56000 158000 260000 362000 464000 566000 668000 770000 872000 974000 57000 159000 261000 363000 465000 567000 669000 771000 873000 975000 58000 160000 262000 364000 466000 568000 670000 772000 874000 976000 59000 161000 263000 365000 467000 569000 671000 773000 875000 977000 60000 162000 264000 366000 468000 570000 672000 774000 876000 978000 61000 163000 265000 367000 469000 571000 673000 775000 877000 979000 62000 164000 266000 368000 470000 572000 674000 776000 878000 980000 63000 165000 267000 369000 471000 573000 675000 777000 879000 981000 64000 166000 268000 370000 472000 574000 676000 778000 880000 982000 65000 167000 269000 371000 473000 575000 677000 779000 881000 983000 66000 168000 270000 372000 474000 576000 678000 780000 882000 984000 67000 169000 271000 373000 475000 577000 679000 781000 883000 985000 68000 170000 272000 374000 476000 578000 680000 782000 884000 986000 69000 171000 273000 375000 477000 579000 681000 783000 885000 987000 70000 172000 274000 376000 478000 580000 682000 784000 886000 988000 71000 173000 275000 377000 479000 581000 683000 785000 887000 989000 72000 174000 276000 378000 480000 582000 684000 786000 888000 990000 73000 175000 277000 379000 481000 583000 685000 787000 889000 991000 74000 176000 278000 380000 482000 584000 686000 788000 890000 992000 75000 177000 279000 381000 483000 585000 687000 789000 891000 993000 76000 178000 280000 382000 484000 586000 688000 790000 892000 994000 77000 179000 281000 383000 485000 587000 689000 791000 893000 995000 78000 180000 282000 384000 486000 588000 690000 792000 894000 996000 79000 181000 283000 385000 487000 589000 691000 793000 895000 997000 80000 182000 284000 386000 488000 590000 692000 794000 896000 998000 81000 183000 285000 387000 489000 591000 693000 795000 897000 999000 82000 184000 286000 388000 490000 592000 694000 796000 898000 1000000 83000 185000 287000 389000 491000 593000 695000 797000 899000 1001000 84000 186000 288000 390000 492000 594000 696000 798000 900000 1002000 85000 187000 289000 391000 493000 595000 697000 799000 901000 1003000 86000 188000 290000 392000 494000 596000 698000 800000 902000 1004000 87000 189000 291000 393000 495000 597000 699000 801000 903000 1005000 88000 190000 292000 394000 496000 598000 700000 802000 904000 1006000 89000 191000 293000 395000 497000 599000 701000 803000 905000 1007000 90000 192000 294000 396000 498000 600000 702000 804000 906000 1008000 91000 193000 295000 397000 499000 601000 703000 805000 907000 1009000 92000 194000 296000 398000 500000 602000 704000 806000 908000 1010000 93000 195000 297000 399000 501000 603000 705000 807000 909000 1011000 94000 196000 298000 400000 502000 604000 706000 808000 910000 1012000 95000 197000 299000 401000 503000 605000 707000 809000 911000 1013000 96000 198000 300000 402000 504000 606000 708000 810000 912000 97000 199000 301000 403000 505000 607000 709000 811000 913000 98000 200000 302000 404000 506000 608000 710000 812000 914000 99000 201000 303000 405000 507000 609000 711000 813000 915000 100000 202000 304000 406000 508000 610000 712000 814000 916000 101000 203000 305000 407000 509000 611000 713000 815000 917000 Choose the time step you are interested in. Only the first one is returned for now. Use the prefetch (split: TIMESTEP 0) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000000000.mat 27-Jul-2023 20:52:07 1.1 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle ...loaded in 0.644 s

1.2 Typical frame

defaultfiles = lamdumpread2(fullfile(datafolder,dumpfile),'default'); % default folder (just for check)
 
%% Data Extraction
% -------------------------------------------------
% Extract the types of atoms and the list of available frames
X0 = lamdumpread2(fullfile(datafolder,dumpfile)); % default frame
The prefetch is split in several files. 1014 TIMESTEPS are availble: Column 01 Column 02 Column 03 Column 04 Column 05 Column 06 Column 07 Column 08 Column 09 Column 10 0 102000 204000 306000 408000 510000 612000 714000 816000 918000 1000 103000 205000 307000 409000 511000 613000 715000 817000 919000 2000 104000 206000 308000 410000 512000 614000 716000 818000 920000 3000 105000 207000 309000 411000 513000 615000 717000 819000 921000 4000 106000 208000 310000 412000 514000 616000 718000 820000 922000 5000 107000 209000 311000 413000 515000 617000 719000 821000 923000 6000 108000 210000 312000 414000 516000 618000 720000 822000 924000 7000 109000 211000 313000 415000 517000 619000 721000 823000 925000 8000 110000 212000 314000 416000 518000 620000 722000 824000 926000 9000 111000 213000 315000 417000 519000 621000 723000 825000 927000 10000 112000 214000 316000 418000 520000 622000 724000 826000 928000 11000 113000 215000 317000 419000 521000 623000 725000 827000 929000 12000 114000 216000 318000 420000 522000 624000 726000 828000 930000 13000 115000 217000 319000 421000 523000 625000 727000 829000 931000 14000 116000 218000 320000 422000 524000 626000 728000 830000 932000 15000 117000 219000 321000 423000 525000 627000 729000 831000 933000 16000 118000 220000 322000 424000 526000 628000 730000 832000 934000 17000 119000 221000 323000 425000 527000 629000 731000 833000 935000 18000 120000 222000 324000 426000 528000 630000 732000 834000 936000 19000 121000 223000 325000 427000 529000 631000 733000 835000 937000 20000 122000 224000 326000 428000 530000 632000 734000 836000 938000 21000 123000 225000 327000 429000 531000 633000 735000 837000 939000 22000 124000 226000 328000 430000 532000 634000 736000 838000 940000 23000 125000 227000 329000 431000 533000 635000 737000 839000 941000 24000 126000 228000 330000 432000 534000 636000 738000 840000 942000 25000 127000 229000 331000 433000 535000 637000 739000 841000 943000 26000 128000 230000 332000 434000 536000 638000 740000 842000 944000 27000 129000 231000 333000 435000 537000 639000 741000 843000 945000 28000 130000 232000 334000 436000 538000 640000 742000 844000 946000 29000 131000 233000 335000 437000 539000 641000 743000 845000 947000 30000 132000 234000 336000 438000 540000 642000 744000 846000 948000 31000 133000 235000 337000 439000 541000 643000 745000 847000 949000 32000 134000 236000 338000 440000 542000 644000 746000 848000 950000 33000 135000 237000 339000 441000 543000 645000 747000 849000 951000 34000 136000 238000 340000 442000 544000 646000 748000 850000 952000 35000 137000 239000 341000 443000 545000 647000 749000 851000 953000 36000 138000 240000 342000 444000 546000 648000 750000 852000 954000 37000 139000 241000 343000 445000 547000 649000 751000 853000 955000 38000 140000 242000 344000 446000 548000 650000 752000 854000 956000 39000 141000 243000 345000 447000 549000 651000 753000 855000 957000 40000 142000 244000 346000 448000 550000 652000 754000 856000 958000 41000 143000 245000 347000 449000 551000 653000 755000 857000 959000 42000 144000 246000 348000 450000 552000 654000 756000 858000 960000 43000 145000 247000 349000 451000 553000 655000 757000 859000 961000 44000 146000 248000 350000 452000 554000 656000 758000 860000 962000 45000 147000 249000 351000 453000 555000 657000 759000 861000 963000 46000 148000 250000 352000 454000 556000 658000 760000 862000 964000 47000 149000 251000 353000 455000 557000 659000 761000 863000 965000 48000 150000 252000 354000 456000 558000 660000 762000 864000 966000 49000 151000 253000 355000 457000 559000 661000 763000 865000 967000 50000 152000 254000 356000 458000 560000 662000 764000 866000 968000 51000 153000 255000 357000 459000 561000 663000 765000 867000 969000 52000 154000 256000 358000 460000 562000 664000 766000 868000 970000 53000 155000 257000 359000 461000 563000 665000 767000 869000 971000 54000 156000 258000 360000 462000 564000 666000 768000 870000 972000 55000 157000 259000 361000 463000 565000 667000 769000 871000 973000 56000 158000 260000 362000 464000 566000 668000 770000 872000 974000 57000 159000 261000 363000 465000 567000 669000 771000 873000 975000 58000 160000 262000 364000 466000 568000 670000 772000 874000 976000 59000 161000 263000 365000 467000 569000 671000 773000 875000 977000 60000 162000 264000 366000 468000 570000 672000 774000 876000 978000 61000 163000 265000 367000 469000 571000 673000 775000 877000 979000 62000 164000 266000 368000 470000 572000 674000 776000 878000 980000 63000 165000 267000 369000 471000 573000 675000 777000 879000 981000 64000 166000 268000 370000 472000 574000 676000 778000 880000 982000 65000 167000 269000 371000 473000 575000 677000 779000 881000 983000 66000 168000 270000 372000 474000 576000 678000 780000 882000 984000 67000 169000 271000 373000 475000 577000 679000 781000 883000 985000 68000 170000 272000 374000 476000 578000 680000 782000 884000 986000 69000 171000 273000 375000 477000 579000 681000 783000 885000 987000 70000 172000 274000 376000 478000 580000 682000 784000 886000 988000 71000 173000 275000 377000 479000 581000 683000 785000 887000 989000 72000 174000 276000 378000 480000 582000 684000 786000 888000 990000 73000 175000 277000 379000 481000 583000 685000 787000 889000 991000 74000 176000 278000 380000 482000 584000 686000 788000 890000 992000 75000 177000 279000 381000 483000 585000 687000 789000 891000 993000 76000 178000 280000 382000 484000 586000 688000 790000 892000 994000 77000 179000 281000 383000 485000 587000 689000 791000 893000 995000 78000 180000 282000 384000 486000 588000 690000 792000 894000 996000 79000 181000 283000 385000 487000 589000 691000 793000 895000 997000 80000 182000 284000 386000 488000 590000 692000 794000 896000 998000 81000 183000 285000 387000 489000 591000 693000 795000 897000 999000 82000 184000 286000 388000 490000 592000 694000 796000 898000 1000000 83000 185000 287000 389000 491000 593000 695000 797000 899000 1001000 84000 186000 288000 390000 492000 594000 696000 798000 900000 1002000 85000 187000 289000 391000 493000 595000 697000 799000 901000 1003000 86000 188000 290000 392000 494000 596000 698000 800000 902000 1004000 87000 189000 291000 393000 495000 597000 699000 801000 903000 1005000 88000 190000 292000 394000 496000 598000 700000 802000 904000 1006000 89000 191000 293000 395000 497000 599000 701000 803000 905000 1007000 90000 192000 294000 396000 498000 600000 702000 804000 906000 1008000 91000 193000 295000 397000 499000 601000 703000 805000 907000 1009000 92000 194000 296000 398000 500000 602000 704000 806000 908000 1010000 93000 195000 297000 399000 501000 603000 705000 807000 909000 1011000 94000 196000 298000 400000 502000 604000 706000 808000 910000 1012000 95000 197000 299000 401000 503000 605000 707000 809000 911000 1013000 96000 198000 300000 402000 504000 606000 708000 810000 912000 97000 199000 301000 403000 505000 607000 709000 811000 913000 98000 200000 302000 404000 506000 608000 710000 812000 914000 99000 201000 303000 405000 507000 609000 711000 813000 915000 100000 202000 304000 406000 508000 610000 712000 814000 916000 101000 203000 305000 407000 509000 611000 713000 815000 917000 Choose the time step you are interested in. Only the first one is returned for now. Use the prefetch (split: TIMESTEP 0) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000000000.mat 27-Jul-2023 20:52:07 1.1 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle ...loaded in 0.386 s
natoms = X0.NUMBER;
timesteps = X0.TIMESTEPS;
atomtypes = unique(X0.ATOMS.type);
ntimesteps = length(timesteps);
 
% Extract the middle frame (i.e. in the middle of the simulation duration)
% X0 is too far from steady state for advanced analysis
Xmiddle = lamdumpread2(fullfile(datafolder,dumpfile),'usesplit',[],timesteps(ceil(ntimesteps/2))); % middle frame
Use the prefetch (split: TIMESTEP 506000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000506000.mat 27-Jul-2023 21:17:15 8.5 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle ...loaded in 0.338 s
 
% Extract the number of beads for each type
% guess the bead type for the fluid (the most populated)
% same guess for the particle (the less populated)
T = X0.ATOMS.type;
natomspertype = arrayfun(@(t) length(find(T==t)),atomtypes);
[~,fluidtype] = max(natomspertype);
[~,solidtype] = min(natomspertype);
walltypes = setdiff(atomtypes,[fluidtype,solidtype]);
 
% Estimate the fluid bead size
% This step uses for more accuracy the buildVerletList()
fluidxyz = X0.ATOMS{T==fluidtype,{'x','y','z'}};
fluidid = X0.ATOMS{T==fluidtype,'id'};
nfluidatoms = length(fluidid);
nsolidatoms = natomspertype(solidtype);
% first estimate assuming that the bead is a cube
boxdims = X0.BOX(:,2) - X0.BOX(:,1);
Vbead_guess = prod(boxdims)/natoms;
rbead_guess = (3/(4*pi)*Vbead_guess)^(1/3);
cutoff = 3*rbead_guess;
[verletList,cutoff,dmin,config,dist] = buildVerletList(fluidxyz,cutoff);
Build Verlet list by searching in blocks... ... done in 0.2316 s with 729 search blocks | minimum distance 2.083e-05 Sort the Verlet list... ... done in 0.504 s buildVerletList: all done in 23.6 s for 215025 atoms
rbead = dmin/2;
 
% find the direction of the flow (largest dimension)
[~,iflow] = max(boxdims);
iothers = setdiff(1:size(X0.BOX,1),iflow);
 
% separate top and bottom walls
vel = {'vx','vy','vz'};
wall1vel = Xmiddle.ATOMS{Xmiddle.ATOMS.type==walltypes(1),vel{iflow}}; wall1vel = wall1vel(1);
wall2vel = Xmiddle.ATOMS{Xmiddle.ATOMS.type==walltypes(2),vel{iflow}}; wall2vel = wall2vel(1);
[wallvel,iwall] = sort([wall1vel,wall2vel],'descend'); % 1 is top (>0), 2 is bottom;
walltypes = walltypes(iwall);
 
% find the position of the particle (i.e., obstacle) placed in the flow (from the middle frame)
solidxyz = Xmiddle.ATOMS{T==solidtype,{'x','y','z'}};
solidid = Xmiddle.ATOMS{T==solidtype,'id'};
solidbox = [min(solidxyz);max(solidxyz)]';
 

2. Time-resolved interpretation: trajectories of selected particles

This section repeats the treatment performed in example1 (files are larger)

2.1 Select particles from one side of the channel

%% Data Preprocessing for Selected Particles
% -------------------------------------------------
% This part is a repetition of the content of example1. It generate the streamlines for particles
% at the left. It requires 5 min and it can be skipped as it not used for stresses.
% note: the duration of the simulation enables the particles in the center (slowest part) to
% cross only half of the cell.
 
% Pick n particles randomly from the left inlet and included in the mask of the solid (initial frame)
n = 300;
tol = 0.5; % add 40% particles around
selectionbox = NaN(3,2);
selectionbox(iflow,:) = [X0.BOX(iflow,1), X0.BOX(iflow,1)+2*rbead];
selectionbox(iothers,:) = (1+tol)*(solidbox(iothers,2)-solidbox(iothers,1))*[-1 1]/2 ...
+ (solidbox(iothers,1)+solidbox(iothers,2)) * [1 1]/2;
ok = true(nfluidatoms,1);
for c = 1:3
ok = ok & (fluidxyz(:,c)>=selectionbox(c,1)) & (fluidxyz(:,c)<=selectionbox(c,2));
end
icandidates = find(ok);
iselected = icandidates(MDunidrnd(length(icandidates),n));
selectedid = fluidid(iselected); % ID to be used
% plot selected particles and other ones
figure, hold on
plot3D(fluidxyz,'b.')
plot3D(fluidxyz(iselected,:),'ro','markerfacecolor','r')
plot3D(solidxyz,'ks','markerfacecolor','k')
view(3), axis equal
 
% Generate the trajectory for the selected particles (for all frames)
% All data are loaded with Lamdumpread2() with a single call using 'usesplit'
% each split is loaded individually and only atoms matching selectedid are stored
% all atoms are sorted as selectedid (no need to sort them)
Xselection = lamdumpread2(fullfile(datafolder,dumpfile),'usesplit',[],timesteps,selectedid);
Use the prefetch (split: TIMESTEP 0) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000000000.mat 27-Jul-2023 20:52:07 1.1 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 1000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000001000.mat 27-Jul-2023 20:52:09 6.3 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 2000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000002000.mat 27-Jul-2023 20:52:12 6.3 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 3000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000003000.mat 27-Jul-2023 20:52:15 6.0 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 4000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000004000.mat 27-Jul-2023 20:52:20 5.9 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 5000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000005000.mat 27-Jul-2023 20:52:23 5.9 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 6000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000006000.mat 27-Jul-2023 20:52:26 5.9 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 7000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000007000.mat 27-Jul-2023 20:52:30 5.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 8000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000008000.mat 27-Jul-2023 20:52:34 5.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 9000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000009000.mat 27-Jul-2023 20:52:38 5.7 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 10000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000010000.mat 27-Jul-2023 20:52:41 5.7 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 11000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000011000.mat 27-Jul-2023 20:52:46 5.7 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 12000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000012000.mat 27-Jul-2023 20:52:50 5.7 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 13000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000013000.mat 27-Jul-2023 20:52:55 5.7 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 14000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000014000.mat 27-Jul-2023 20:52:59 6.2 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 15000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000015000.mat 27-Jul-2023 20:53:03 6.5 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 16000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000016000.mat 27-Jul-2023 20:53:08 6.5 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 17000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000017000.mat 27-Jul-2023 20:53:13 6.5 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 18000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000018000.mat 27-Jul-2023 20:53:17 6.5 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 19000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000019000.mat 27-Jul-2023 20:53:21 6.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 20000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000020000.mat 27-Jul-2023 20:53:26 6.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 21000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000021000.mat 27-Jul-2023 20:53:30 6.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 22000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000022000.mat 27-Jul-2023 20:53:35 6.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 23000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000023000.mat 27-Jul-2023 20:53:39 6.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 24000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000024000.mat 27-Jul-2023 20:53:44 6.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 25000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000025000.mat 27-Jul-2023 20:53:49 6.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 26000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000026000.mat 27-Jul-2023 20:53:53 6.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 27000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000027000.mat 27-Jul-2023 20:53:57 6.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 28000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000028000.mat 27-Jul-2023 20:54:01 6.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 29000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000029000.mat 27-Jul-2023 20:54:04 6.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 30000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000030000.mat 27-Jul-2023 20:54:08 6.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 31000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000031000.mat 27-Jul-2023 20:54:10 6.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 32000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000032000.mat 27-Jul-2023 20:54:15 6.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 33000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000033000.mat 27-Jul-2023 20:54:19 6.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 34000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000034000.mat 27-Jul-2023 20:54:23 6.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 35000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000035000.mat 27-Jul-2023 20:54:26 6.7 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 36000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000036000.mat 27-Jul-2023 20:54:30 6.7 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 37000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000037000.mat 27-Jul-2023 20:54:33 6.7 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 38000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000038000.mat 27-Jul-2023 20:54:37 6.7 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 39000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000039000.mat 27-Jul-2023 20:54:40 6.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 40000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000040000.mat 27-Jul-2023 20:54:43 6.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 41000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000041000.mat 27-Jul-2023 20:54:47 6.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 42000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000042000.mat 27-Jul-2023 20:54:51 6.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 43000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000043000.mat 27-Jul-2023 20:54:54 6.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 44000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000044000.mat 27-Jul-2023 20:54:57 6.9 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 45000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000045000.mat 27-Jul-2023 20:55:01 6.9 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 46000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000046000.mat 27-Jul-2023 20:55:05 6.9 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 47000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000047000.mat 27-Jul-2023 20:55:09 6.9 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 48000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000048000.mat 27-Jul-2023 20:55:11 6.9 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 49000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000049000.mat 27-Jul-2023 20:55:15 7.0 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 50000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000050000.mat 27-Jul-2023 20:55:20 7.0 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 51000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000051000.mat 27-Jul-2023 20:55:24 7.0 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 52000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000052000.mat 27-Jul-2023 20:55:26 7.0 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 53000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000053000.mat 27-Jul-2023 20:55:30 7.1 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 54000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000054000.mat 27-Jul-2023 20:55:34 7.1 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 55000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000055000.mat 27-Jul-2023 20:55:38 7.1 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 56000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000056000.mat 27-Jul-2023 20:55:41 7.1 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 57000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000057000.mat 27-Jul-2023 20:55:45 7.1 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 58000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000058000.mat 27-Jul-2023 20:55:49 7.2 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 59000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000059000.mat 27-Jul-2023 20:55:53 7.2 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 60000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000060000.mat 27-Jul-2023 20:55:56 7.2 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 61000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000061000.mat 27-Jul-2023 20:56:00 7.2 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 62000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000062000.mat 27-Jul-2023 20:56:05 7.3 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 63000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000063000.mat 27-Jul-2023 20:56:10 7.3 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 64000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000064000.mat 27-Jul-2023 20:56:14 7.3 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 65000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000065000.mat 27-Jul-2023 20:56:19 7.3 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 66000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000066000.mat 27-Jul-2023 20:56:25 7.3 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 67000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000067000.mat 27-Jul-2023 20:56:30 7.4 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 68000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000068000.mat 27-Jul-2023 20:56:34 7.4 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 69000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000069000.mat 27-Jul-2023 20:56:38 7.4 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 70000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000070000.mat 27-Jul-2023 20:56:43 7.4 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 71000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000071000.mat 27-Jul-2023 20:56:48 7.4 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 72000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000072000.mat 27-Jul-2023 20:56:52 7.4 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 73000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000073000.mat 27-Jul-2023 20:56:57 7.5 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 74000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000074000.mat 27-Jul-2023 20:57:02 7.5 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 75000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000075000.mat 27-Jul-2023 20:57:06 7.5 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 76000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000076000.mat 27-Jul-2023 20:57:10 7.5 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 77000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000077000.mat 27-Jul-2023 20:57:15 7.5 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 78000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000078000.mat 27-Jul-2023 20:57:19 7.5 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 79000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000079000.mat 27-Jul-2023 20:57:24 7.5 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 80000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000080000.mat 27-Jul-2023 20:57:28 7.5 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 81000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000081000.mat 27-Jul-2023 20:57:32 7.5 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 82000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000082000.mat 27-Jul-2023 20:57:37 7.5 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 83000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000083000.mat 27-Jul-2023 20:57:42 7.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 84000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000084000.mat 27-Jul-2023 20:57:46 7.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 85000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000085000.mat 27-Jul-2023 20:57:49 7.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 86000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000086000.mat 27-Jul-2023 20:57:53 7.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 87000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000087000.mat 27-Jul-2023 20:57:56 7.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 88000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000088000.mat 27-Jul-2023 20:58:00 7.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 89000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000089000.mat 27-Jul-2023 20:58:02 7.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 90000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000090000.mat 27-Jul-2023 20:58:06 7.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 91000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000091000.mat 27-Jul-2023 20:58:10 7.7 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 92000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000092000.mat 27-Jul-2023 20:58:14 7.7 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 93000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000093000.mat 27-Jul-2023 20:58:18 7.7 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 94000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000094000.mat 27-Jul-2023 20:58:22 7.7 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 95000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000095000.mat 27-Jul-2023 20:58:26 7.7 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 96000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000096000.mat 27-Jul-2023 20:58:29 7.7 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 97000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000097000.mat 27-Jul-2023 20:58:32 7.7 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 98000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000098000.mat 27-Jul-2023 20:58:36 7.7 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 99000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000099000.mat 27-Jul-2023 20:58:40 7.7 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 100000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000100000.mat 27-Jul-2023 20:58:44 7.7 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 101000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000101000.mat 27-Jul-2023 20:58:47 7.7 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 102000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000102000.mat 27-Jul-2023 20:58:50 7.7 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 103000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000103000.mat 27-Jul-2023 20:58:54 7.7 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 104000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000104000.mat 27-Jul-2023 20:58:57 7.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 105000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000105000.mat 27-Jul-2023 20:59:00 7.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 106000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000106000.mat 27-Jul-2023 20:59:03 7.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 107000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000107000.mat 27-Jul-2023 20:59:07 7.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 108000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000108000.mat 27-Jul-2023 20:59:11 7.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 109000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000109000.mat 27-Jul-2023 20:59:13 7.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 110000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000110000.mat 27-Jul-2023 20:59:16 7.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 111000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000111000.mat 27-Jul-2023 20:59:20 7.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 112000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000112000.mat 27-Jul-2023 20:59:24 7.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 113000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000113000.mat 27-Jul-2023 20:59:27 7.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 114000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000114000.mat 27-Jul-2023 20:59:31 7.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 115000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000115000.mat 27-Jul-2023 20:59:34 7.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 116000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000116000.mat 27-Jul-2023 20:59:38 7.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 117000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000117000.mat 27-Jul-2023 20:59:42 7.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 118000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000118000.mat 27-Jul-2023 20:59:47 7.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 119000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000119000.mat 27-Jul-2023 20:59:51 7.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 120000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000120000.mat 27-Jul-2023 20:59:55 7.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 121000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000121000.mat 27-Jul-2023 20:59:59 7.9 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 122000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000122000.mat 27-Jul-2023 21:00:03 7.9 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 123000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000123000.mat 27-Jul-2023 21:00:07 7.9 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 124000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000124000.mat 27-Jul-2023 21:00:11 7.9 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 125000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000125000.mat 27-Jul-2023 21:00:15 7.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 126000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000126000.mat 27-Jul-2023 21:00:20 7.9 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 127000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000127000.mat 27-Jul-2023 21:00:23 7.9 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 128000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000128000.mat 27-Jul-2023 21:00:27 7.9 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 129000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000129000.mat 27-Jul-2023 21:00:30 7.9 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 130000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000130000.mat 27-Jul-2023 21:00:34 7.9 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 131000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000131000.mat 27-Jul-2023 21:00:38 7.9 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 132000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000132000.mat 27-Jul-2023 21:00:41 7.9 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 133000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000133000.mat 27-Jul-2023 21:00:44 7.9 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 134000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000134000.mat 27-Jul-2023 21:00:48 7.9 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 135000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000135000.mat 27-Jul-2023 21:00:50 7.9 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 136000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000136000.mat 27-Jul-2023 21:00:53 7.9 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 137000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000137000.mat 2...
 
% Collect the trajectory of the solid particle for all frames
Xsolid = lamdumpread2(fullfile(datafolder,dumpfile),'usesplit',[],timesteps,solidid);
Use the prefetch (split: TIMESTEP 0) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000000000.mat 27-Jul-2023 20:52:07 1.1 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 1000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000001000.mat 27-Jul-2023 20:52:09 6.3 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 2000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000002000.mat 27-Jul-2023 20:52:12 6.3 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 3000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000003000.mat 27-Jul-2023 20:52:15 6.0 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 4000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000004000.mat 27-Jul-2023 20:52:20 5.9 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 5000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000005000.mat 27-Jul-2023 20:52:23 5.9 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 6000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000006000.mat 27-Jul-2023 20:52:26 5.9 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 7000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000007000.mat 27-Jul-2023 20:52:30 5.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 8000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000008000.mat 27-Jul-2023 20:52:34 5.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 9000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000009000.mat 27-Jul-2023 20:52:38 5.7 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 10000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000010000.mat 27-Jul-2023 20:52:41 5.7 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 11000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000011000.mat 27-Jul-2023 20:52:46 5.7 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 12000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000012000.mat 27-Jul-2023 20:52:50 5.7 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 13000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000013000.mat 27-Jul-2023 20:52:55 5.7 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 14000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000014000.mat 27-Jul-2023 20:52:59 6.2 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 15000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000015000.mat 27-Jul-2023 20:53:03 6.5 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 16000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000016000.mat 27-Jul-2023 20:53:08 6.5 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 17000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000017000.mat 27-Jul-2023 20:53:13 6.5 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 18000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000018000.mat 27-Jul-2023 20:53:17 6.5 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 19000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000019000.mat 27-Jul-2023 20:53:21 6.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 20000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000020000.mat 27-Jul-2023 20:53:26 6.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 21000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000021000.mat 27-Jul-2023 20:53:30 6.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 22000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000022000.mat 27-Jul-2023 20:53:35 6.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 23000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000023000.mat 27-Jul-2023 20:53:39 6.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 24000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000024000.mat 27-Jul-2023 20:53:44 6.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 25000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000025000.mat 27-Jul-2023 20:53:49 6.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 26000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000026000.mat 27-Jul-2023 20:53:53 6.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 27000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000027000.mat 27-Jul-2023 20:53:57 6.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 28000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000028000.mat 27-Jul-2023 20:54:01 6.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 29000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000029000.mat 27-Jul-2023 20:54:04 6.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 30000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000030000.mat 27-Jul-2023 20:54:08 6.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 31000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000031000.mat 27-Jul-2023 20:54:10 6.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 32000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000032000.mat 27-Jul-2023 20:54:15 6.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 33000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000033000.mat 27-Jul-2023 20:54:19 6.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 34000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000034000.mat 27-Jul-2023 20:54:23 6.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 35000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000035000.mat 27-Jul-2023 20:54:26 6.7 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 36000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000036000.mat 27-Jul-2023 20:54:30 6.7 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 37000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000037000.mat 27-Jul-2023 20:54:33 6.7 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 38000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000038000.mat 27-Jul-2023 20:54:37 6.7 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 39000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000039000.mat 27-Jul-2023 20:54:40 6.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 40000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000040000.mat 27-Jul-2023 20:54:43 6.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 41000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000041000.mat 27-Jul-2023 20:54:47 6.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 42000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000042000.mat 27-Jul-2023 20:54:51 6.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 43000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000043000.mat 27-Jul-2023 20:54:54 6.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 44000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000044000.mat 27-Jul-2023 20:54:57 6.9 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 45000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000045000.mat 27-Jul-2023 20:55:01 6.9 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 46000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000046000.mat 27-Jul-2023 20:55:05 6.9 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 47000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000047000.mat 27-Jul-2023 20:55:09 6.9 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 48000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000048000.mat 27-Jul-2023 20:55:11 6.9 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 49000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000049000.mat 27-Jul-2023 20:55:15 7.0 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 50000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000050000.mat 27-Jul-2023 20:55:20 7.0 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 51000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000051000.mat 27-Jul-2023 20:55:24 7.0 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 52000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000052000.mat 27-Jul-2023 20:55:26 7.0 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 53000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000053000.mat 27-Jul-2023 20:55:30 7.1 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 54000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000054000.mat 27-Jul-2023 20:55:34 7.1 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 55000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000055000.mat 27-Jul-2023 20:55:38 7.1 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 56000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000056000.mat 27-Jul-2023 20:55:41 7.1 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 57000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000057000.mat 27-Jul-2023 20:55:45 7.1 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 58000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000058000.mat 27-Jul-2023 20:55:49 7.2 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 59000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000059000.mat 27-Jul-2023 20:55:53 7.2 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 60000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000060000.mat 27-Jul-2023 20:55:56 7.2 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 61000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000061000.mat 27-Jul-2023 20:56:00 7.2 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 62000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000062000.mat 27-Jul-2023 20:56:05 7.3 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 63000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000063000.mat 27-Jul-2023 20:56:10 7.3 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 64000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000064000.mat 27-Jul-2023 20:56:14 7.3 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 65000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000065000.mat 27-Jul-2023 20:56:19 7.3 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 66000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000066000.mat 27-Jul-2023 20:56:25 7.3 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 67000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000067000.mat 27-Jul-2023 20:56:30 7.4 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 68000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000068000.mat 27-Jul-2023 20:56:34 7.4 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 69000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000069000.mat 27-Jul-2023 20:56:38 7.4 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 70000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000070000.mat 27-Jul-2023 20:56:43 7.4 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 71000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000071000.mat 27-Jul-2023 20:56:48 7.4 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 72000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000072000.mat 27-Jul-2023 20:56:52 7.4 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 73000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000073000.mat 27-Jul-2023 20:56:57 7.5 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 74000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000074000.mat 27-Jul-2023 20:57:02 7.5 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 75000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000075000.mat 27-Jul-2023 20:57:06 7.5 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 76000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000076000.mat 27-Jul-2023 20:57:10 7.5 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 77000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000077000.mat 27-Jul-2023 20:57:15 7.5 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 78000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000078000.mat 27-Jul-2023 20:57:19 7.5 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 79000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000079000.mat 27-Jul-2023 20:57:24 7.5 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 80000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000080000.mat 27-Jul-2023 20:57:28 7.5 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 81000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000081000.mat 27-Jul-2023 20:57:32 7.5 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 82000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000082000.mat 27-Jul-2023 20:57:37 7.5 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 83000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000083000.mat 27-Jul-2023 20:57:42 7.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 84000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000084000.mat 27-Jul-2023 20:57:46 7.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 85000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000085000.mat 27-Jul-2023 20:57:49 7.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 86000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000086000.mat 27-Jul-2023 20:57:53 7.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 87000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000087000.mat 27-Jul-2023 20:57:56 7.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 88000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000088000.mat 27-Jul-2023 20:58:00 7.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 89000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000089000.mat 27-Jul-2023 20:58:02 7.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 90000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000090000.mat 27-Jul-2023 20:58:06 7.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 91000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000091000.mat 27-Jul-2023 20:58:10 7.7 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 92000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000092000.mat 27-Jul-2023 20:58:14 7.7 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 93000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000093000.mat 27-Jul-2023 20:58:18 7.7 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 94000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000094000.mat 27-Jul-2023 20:58:22 7.7 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 95000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000095000.mat 27-Jul-2023 20:58:26 7.7 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 96000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000096000.mat 27-Jul-2023 20:58:29 7.7 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 97000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000097000.mat 27-Jul-2023 20:58:32 7.7 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 98000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000098000.mat 27-Jul-2023 20:58:36 7.7 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 99000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000099000.mat 27-Jul-2023 20:58:40 7.7 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 100000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000100000.mat 27-Jul-2023 20:58:44 7.7 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 101000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000101000.mat 27-Jul-2023 20:58:47 7.7 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 102000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000102000.mat 27-Jul-2023 20:58:50 7.7 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 103000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000103000.mat 27-Jul-2023 20:58:54 7.7 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 104000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000104000.mat 27-Jul-2023 20:58:57 7.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 105000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000105000.mat 27-Jul-2023 20:59:00 7.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 106000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000106000.mat 27-Jul-2023 20:59:03 7.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 107000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000107000.mat 27-Jul-2023 20:59:07 7.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 108000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000108000.mat 27-Jul-2023 20:59:11 7.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 109000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000109000.mat 27-Jul-2023 20:59:13 7.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 110000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000110000.mat 27-Jul-2023 20:59:16 7.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 111000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000111000.mat 27-Jul-2023 20:59:20 7.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 112000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000112000.mat 27-Jul-2023 20:59:24 7.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 113000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000113000.mat 27-Jul-2023 20:59:27 7.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 114000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000114000.mat 27-Jul-2023 20:59:31 7.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 115000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000115000.mat 27-Jul-2023 20:59:34 7.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 116000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000116000.mat 27-Jul-2023 20:59:38 7.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 117000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000117000.mat 27-Jul-2023 20:59:42 7.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 118000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000118000.mat 27-Jul-2023 20:59:47 7.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 119000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000119000.mat 27-Jul-2023 20:59:51 7.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 120000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000120000.mat 27-Jul-2023 20:59:55 7.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 121000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000121000.mat 27-Jul-2023 20:59:59 7.9 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 122000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000122000.mat 27-Jul-2023 21:00:03 7.9 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 123000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000123000.mat 27-Jul-2023 21:00:07 7.9 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 124000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000124000.mat 27-Jul-2023 21:00:11 7.9 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 125000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000125000.mat 27-Jul-2023 21:00:15 7.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 126000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000126000.mat 27-Jul-2023 21:00:20 7.9 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 127000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000127000.mat 27-Jul-2023 21:00:23 7.9 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 128000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000128000.mat 27-Jul-2023 21:00:27 7.9 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 129000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000129000.mat 27-Jul-2023 21:00:30 7.9 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 130000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000130000.mat 27-Jul-2023 21:00:34 7.9 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 131000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000131000.mat 27-Jul-2023 21:00:38 7.9 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 132000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000132000.mat 27-Jul-2023 21:00:41 7.9 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 133000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000133000.mat 27-Jul-2023 21:00:44 7.9 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 134000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000134000.mat 27-Jul-2023 21:00:48 7.9 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 135000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000135000.mat 27-Jul-2023 21:00:50 7.9 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 136000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000136000.mat 27-Jul-2023 21:00:53 7.9 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle Use the prefetch (split: TIMESTEP 137000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000137000.mat 2...

2.3 Main loop to collect trajectories

 
% Store the trajectories stored in a ntimesteps x 3 x n matrix
% missing data will
[seltraj,selveloc] = deal(NaN(ntimesteps,3,n,'single'));
solidtraj = NaN(ntimesteps,3,nsolidatoms,'double');
goodframes = true(ntimesteps,1);
for it = 1:ntimesteps
% read fluid atoms
selframe = Xselection.ATOMS(Xselection.ATOMS{:,'TIMESTEP'} == timesteps(it),:);
nfoundatoms = size(selframe,1);
if nfoundatoms<n % incomplete dumped frame (it may happen in LAMMPS)
dispf('incomplete frame %d/%d (t=%0.4g): %d of %d fluid atoms missing',it,ntimesteps,timesteps(it),n-nfoundatoms,n)
[~,iatoms,jatoms] = intersect(selectedid,selframe.id,'stable');
seltraj(it,:,iatoms) = permute(selframe{jatoms,{'x','y','z'}},[3 2 1]);
goodframes(it) = false;
else
seltraj(it,:,:) = permute(selframe{:,{'x','y','z'}},[3 2 1]);
selveloc(it,:,:) = permute(selframe{:,{'vx','vy','vz'}},[3 2 1]);
end
% read solid atoms
solidframe = Xsolid.ATOMS(Xsolid.ATOMS{:,'TIMESTEP'} == timesteps(it),:);
nfoundatoms = size(solidframe,1);
if nfoundatoms<nsolidatoms % incomplete dumped frame (it may happen in LAMMPS)
dispf('incomplete frame %d/%d (t=%0.4g): %d of %d solid atoms missing',it,ntimesteps,timesteps(it),nsolidatoms-nfoundatoms,nsolidatoms)
[~,iatoms,jatoms] = intersect(solidid,solidframe.id,'stable');
solidtraj(it,:,iatoms) = permute(solidframe{jatoms,{'x','y','z'}},[3 2 1]);
goodframes(it) = false;
else
solidtraj(it,:,:) = permute(solidframe{:,{'x','y','z'}},[3 2 1]);
end
end
% velocity magnitude for the fluid particles
selveloc_magnitude = squeeze(sqrt(sum(selveloc.^2, 2)));
 

2.4 Trajectory visualization

% Trajectory Visualization
% -------------------------------------------------
% streamlines with color representing velocity magnitude
figure, hold on
col = parula(n);
for i=1:n
jumps = [1;ntimesteps+1];
for d=1:3
jumps = unique([jumps;find(abs(diff(seltraj(:,d,i)))>boxdims(d)/2)+1]);
end
for j=1:length(jumps)-1
u = jumps(j):jumps(j+1)-1;
streamline = seltraj(u, :, i);
color_line3(streamline(:, 1), streamline(:, 2), streamline(:, 3), selveloc_magnitude(u, i), 'linewidth', 3);
% faster method but without streamline
%plot3(traj(u,1,i),traj(u,2,i),traj(u,3,i),'-','linewidth',3,'color',col(i,:))
end
end
plot3D(solidxyz,'ko','markerfacecolor','k','markersize',5)
view(3), axis equal
title({'\rm{' strrep(dumpfile, '_', '\_') '}:' '\bf{Velocity of Selected Particles Along Streamlines}'}, 'Interpreter', 'tex')
xlabel('x (m)'), ylabel('Y (m)'), zlabel('z (m)')
hc = colorbar; hc.Label.String = 'velocity (m/s)';

3. Stress Analysis

This section is progressive and combine a broad variety of approaches with increasing refinements and sophistications
 
%% Stress Analysis in Reference Frame (PART 2)
% -------------------------------------------------
% This part focuses on evaluating the mechanical interactions and stress in a particular
% frame. Specifically, we employ a late-stage frame (closer to steady-state conditions) to
% derive insights into the types of contacts between fluid and solid atoms.
%
% * `timestepforstress`: Time step selected for stress analysis.
% * `Xstress`: Data structure containing all relevant atomic information at the chosen time step.
% * `verletListCross`: A specialized Verlet list that segregates contacts based on atom types.
% * `isincontact`, `isincontactwithsolid`, `isincontactwithfluid`: Boolean fields that signify
% whether an atom is in contact, and with which type of atom.
%
% Here we rely heavily on a short-cut-off Verlet list for efficiency, separating cross-terms
% and characterizing specific types of atomic interactions.
clear Xsolid Xmiddle X0 Xselection streamline
 

3.1 Initialization and possibly loop on several frames

% General container for storing data for different timesteps
postdata = repmat(struct('type','Landshoff|Hertz','timestep',NaN,'Xfluid',[],'Xsolid',[],'force',[]),0,0);
 
 
% Extract the Specific Frame for Stress Interpretation
% Selects a time step that is presumably closer to the steady-state to perform stress analysis.
list_timestepforstess = unique(timesteps(ceil((0.1:0.1:0.9)*ntimesteps)));
nlist_timestepforstess = length(list_timestepforstess);
consideredtimesteps = 1:nlist_timestepforstess;
 
% === Loop on all selected time steps ===
altxt = {'no plot','plot'};
 
 
for i_timestepforstess = consideredtimesteps
 
 
% current frame
timestepforstress = list_timestepforstess(i_timestepforstess);
Xstress = lamdumpread2(fullfile(datafolder,dumpfile),'usesplit',[],timestepforstress); % middle frame
doplot = ismember(timestepforstress,list_timestepforstess([1 end]));
dispf('\n%s-\n[%d/%d] TIMESTEP = %d (%s)',repmat('-*',1,40),i_timestepforstess,nlist_timestepforstess,timestepforstress,altxt{doplot+1})
 
 
% Categorization of Atom Types
% Identifies the most and least populous atom types to discern between fluid and solid types.
T = Xstress.ATOMS.type;
natomspertype = arrayfun(@(t) length(find(T==t)),atomtypes);
[~,fluidtype] = max(natomspertype);
[~,solidtype] = min(natomspertype);
walltypes = setdiff(atomtypes,[fluidtype,solidtype]);
 
% Boolean Classification of Atom Types
Xstress.ATOMS.isfluid = Xstress.ATOMS.type==fluidtype;
Xstress.ATOMS.issolid = Xstress.ATOMS.type==solidtype;
 
% Verlet List Construction with Short Cutoff
% Builds a Verlet list with a short cutoff distance, designed to identify only the closest neighbors.
[verletList,cutoff,dmin,config,dist] = buildVerletList(Xstress.ATOMS,3*rbead);
 
% Partition Verlet List Based on Atom Types
% This Verlet list is partitioned based on atom types, distinguishing between interactions
% that are exclusively fluid-fluid, solid-fluid, or solid-solid.
verletListCross = partitionVerletList(verletList,Xstress.ATOMS);
 
% Identify Contacting Atoms
% The atoms that are in contact with each other are identified.
Xstress.ATOMS.isincontact = ~cellfun(@isempty,verletListCross);
Xstress.ATOMS.contacttypes = cellfun(@(v) Xstress.ATOMS.type(v)',verletListCross,'UniformOutput',false);
 
% Identify Atoms in Contact with Solids and Fluids
% Further filters the atoms in contact to identify which are in contact with solids and which with fluids.
Xstress.ATOMS.isincontactwithsolid = cellfun(@(c) ismember(solidtype,c), Xstress.ATOMS.contacttypes);
Xstress.ATOMS.isincontactwithfluid = cellfun(@(c) ismember(fluidtype,c), Xstress.ATOMS.contacttypes);
 
% Flag Fluid Atoms in Contact with Solid and Vice Versa
% Identifies fluid atoms that are in contact with solid atoms and solid atoms in contact with fluid atoms.
Xstress.ATOMS.fluidincontactwithsolid = Xstress.ATOMS.isfluid & Xstress.ATOMS.isincontactwithsolid;
Xstress.ATOMS.solidincontactwithfluid = Xstress.ATOMS.issolid & Xstress.ATOMS.isincontactwithfluid;
 
% Identify Indices for Analysis
% Indices of fluid and solid atoms that are in contact with each other.
ifluidcontact = find(Xstress.ATOMS.fluidincontactwithsolid);
isolidcontact = find(Xstress.ATOMS.solidincontactwithfluid);
 
% Extend Fluid Atoms Set with Neighbors (Not Directly in Contact)
% Extends the set of fluid atoms under consideration to include their closest neighbors.
ifluidcontact_withneighbors = unique(union(ifluidcontact,cat(2,verletList{ifluidcontact}))); % include neighbors
ifluidcontact_withneighbors = ifluidcontact_withneighbors(T(ifluidcontact_withneighbors)==fluidtype); % filter, keep only fluids
[~,~,ind_withoutneighbors] = intersect(ifluidcontact,ifluidcontact_withneighbors,'stable');
notneighboringcontacts = setdiff(ifluidcontact_withneighbors,ifluidcontact_withneighbors(ind_withoutneighbors));
 
% Visualization for some frames
% filled red: solid "atoms" in contact with the fluid
% filled blue: fluid "atoms" in contact with the solid (crown/shell)
% empty blue: fluid "atoms" neighbor of previous ones but not contact
% empty blue are used to calculate Landshoff forces at the position of filled ones
if doplot
figure, hold on
plot3D(Xstress.ATOMS{ifluidcontact,{'x','y','z'}},'bo','markerfacecolor','b')
plot3D(Xstress.ATOMS{notneighboringcontacts,{'x','y','z'}},'bo','markerfacecolor','w')
plot3D(Xstress.ATOMS{isolidcontact,{'x','y','z'}},'ro','markersize',12,'markerfacecolor','r')
axis equal, view(3), drawnow
dispf('\tnumber of fluid atoms in contact with solid %d', length(ifluidcontact))
dispf('\tnumber of solid atoms in contact with fluid %d (neighbors %d)', length(isolidcontact),length(notneighboringcontacts))
end
end
Use the prefetch (split: TIMESTEP 101000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000101000.mat 27-Jul-2023 20:58:47 7.7 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle ...loaded in 0.29 s
-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- [1/9] TIMESTEP = 101000 (plot)
Build Verlet list by searching in blocks... ... done in 0.1869 s with 1584 search blocks | minimum distance 1.54e-05 Sort the Verlet list... ... done in 0.647 s buildVerletList: all done in 27.3 s for 313344 atoms
number of fluid atoms in contact with solid 853
number of solid atoms in contact with fluid 566 (neighbors 1094)
Use the prefetch (split: TIMESTEP 202000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000202000.mat 27-Jul-2023 21:03:44 8.1 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle ...loaded in 0.277 s
-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- [2/9] TIMESTEP = 202000 (no plot)
Build Verlet list by searching in blocks... ... done in 0.1811 s with 1584 search blocks | minimum distance 1.506e-05 Sort the Verlet list... ... done in 0.676 s buildVerletList: all done in 28.4 s for 313344 atoms
Use the prefetch (split: TIMESTEP 304000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000304000.mat 27-Jul-2023 21:08:00 8.3 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle ...loaded in 0.27 s
-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- [3/9] TIMESTEP = 304000 (no plot)
Build Verlet list by searching in blocks... ... done in 0.2024 s with 1584 search blocks | minimum distance 1.441e-05 Sort the Verlet list... ... done in 0.666 s buildVerletList: all done in 27.4 s for 313344 atoms
Use the prefetch (split: TIMESTEP 405000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000405000.mat 27-Jul-2023 21:12:42 8.4 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle ...loaded in 0.269 s
-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- [4/9] TIMESTEP = 405000 (no plot)
Build Verlet list by searching in blocks... ... done in 0.1682 s with 1584 search blocks | minimum distance 1.37e-05 Sort the Verlet list... ... done in 0.677 s buildVerletList: all done in 27.9 s for 313344 atoms
Use the prefetch (split: TIMESTEP 506000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000506000.mat 27-Jul-2023 21:17:15 8.5 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle ...loaded in 0.272 s
-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- [5/9] TIMESTEP = 506000 (no plot)
Build Verlet list by searching in blocks... ... done in 0.1637 s with 1584 search blocks | minimum distance 1.479e-05 Sort the Verlet list... ... done in 0.674 s buildVerletList: all done in 27.4 s for 313344 atoms
Use the prefetch (split: TIMESTEP 608000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000608000.mat 27-Jul-2023 21:21:52 8.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle ...loaded in 0.271 s
-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- [6/9] TIMESTEP = 608000 (no plot)
Build Verlet list by searching in blocks... ... done in 0.2936 s with 1584 search blocks | minimum distance 1.476e-05 Sort the Verlet list... ... done in 0.684 s buildVerletList: all done in 27.5 s for 313344 atoms
Use the prefetch (split: TIMESTEP 709000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000709000.mat 27-Jul-2023 21:26:29 8.6 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle ...loaded in 0.274 s
-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- [7/9] TIMESTEP = 709000 (no plot)
Build Verlet list by searching in blocks... ... done in 0.2175 s with 1584 search blocks | minimum distance 1.364e-05 Sort the Verlet list... ... done in 0.687 s buildVerletList: all done in 27.3 s for 313344 atoms
Use the prefetch (split: TIMESTEP 811000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000811000.mat 27-Jul-2023 21:31:05 8.7 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle ...loaded in 0.275 s
-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- [8/9] TIMESTEP = 811000 (no plot)
Build Verlet list by searching in blocks... ... done in 0.05926 s with 1584 search blocks | minimum distance 1.456e-05 Sort the Verlet list... ... done in 0.682 s buildVerletList: all done in 27.3 s for 313344 atoms
Use the prefetch (split: TIMESTEP 912000) folder (instead of '../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle')... TIMESTEP_000912000.mat 27-Jul-2023 21:35:16 8.8 MBytes ../dumps/pub1/numericalViscosimeter_reference_ulsphBulk_hertzBoundary/PREFETCH_dump.ulsphBulk_hertzBoundary_referenceParameterExponent+1_with1SuspendedParticle ...loaded in 0.278 s
-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*- [9/9] TIMESTEP = 912000 (plot)
Build Verlet list by searching in blocks... ... done in 0.2728 s with 1584 search blocks | minimum distance 1.461e-05 Sort the Verlet list... ... done in 0.689 s buildVerletList: all done in 27.4 s for 313344 atoms
number of fluid atoms in contact with solid 886
number of solid atoms in contact with fluid 580 (neighbors 1035)
 
 

3.2. Projection of velocity field on a regular grid

The interpolation process of large systems requires to be performed with interp3SPHVerlet() with the help of a Verletlist. The search for neighbors can be done for beads but also for arbitrary nodes from a mesh. The interpolation can be applied to scalar and vectorial fields. Once velocities are propertly calculated, it is possible to derive streamlines.
With large grids, the calculations are very intensive and may require a lot of memory. Reduce the grid size if your experience memory errors.
%% === 3D view of the fluid velocity around the solid ===
DOINTENSIVECALC = true; % set it to false, to prevent intensive calculations
% Definition of the Viewbox (for all grid-based visualizations)
coords = {'x','y','z'}; % Cartesian coordinates
vcoords = cellfun(@(c) ['v',c],coords,'UniformOutput',false); % vx, vy, vz
% Extract the box around solid atoms and include fluid ones
fluidbox = [ min(Xstress.ATOMS{Xstress.ATOMS.isfluid,coords})
max(Xstress.ATOMS{Xstress.ATOMS.isfluid,coords}) ]';
solidbox = [ min(Xstress.ATOMS{Xstress.ATOMS.issolid,coords})
max(Xstress.ATOMS{Xstress.ATOMS.issolid,coords}) ]';
[~,iflow] = max(max(Xstress.ATOMS{Xstress.ATOMS.isfluid,vcoords})); % direction of the flow
viewbox = fluidbox; viewbox(iflow,:) = solidbox(iflow,:);
if DOINTENSIVECALC % THIS SECTION IS VERY INTENSIVE (more than 5 mins)
% 3D Cartesian Grid from the widow xw, yw, zw
nresolution = [300 300 30]/2; % /2 for the notebook (removed in a real Matlab code)
resolutionmax = max(nresolution);
for icoord = 1:3
viewbox(icoord,:) = mean(viewbox(icoord,:)) + [-1.2 1.2]*diff(viewbox(icoord,:))/2*nresolution(icoord)/resolutionmax;
viewbox(icoord,1) = max(viewbox(icoord,1),fluidbox(icoord,1));
viewbox(icoord,2) = min(viewbox(icoord,2),fluidbox(icoord,2));
end
hLandshoff = 5*rbead; %1.25e-5; % m
xw = linspace(viewbox(1,1),viewbox(1,2),nresolution(1));
yw = linspace(viewbox(2,1),viewbox(2,2),nresolution(2));
zw = linspace(viewbox(3,1),viewbox(3,2),nresolution(3));
insidewindow = Xstress.ATOMS.isfluid;
for icoord = 1:3
insidewindow = insidewindow ...
& Xstress.ATOMS{:,coords{icoord}}>=viewbox(icoord,1) ...
& Xstress.ATOMS{:,coords{icoord}}<=viewbox(icoord,2);
end
dispf('%d fluid atoms have be found in the cross section around the solid',length(find(insidewindow)));
[Xw,Yw,Zw] = meshgrid(xw,yw,zw);
XYZgrid = [Xw(:),Yw(:),Zw(:)];
 
% Verlet Grid List: lists the beads close to a grid node/vertex
XYZ = Xstress.ATOMS{insidewindow,coords}; % kernel centers
VXYZ = buildVerletList({XYZgrid XYZ},1.2*hLandshoff); % special grid syntax
 
% Interpolation of all velocity components (vx, vy, vz)
% v3XYZgrid is the interpolated vectorial field
% vxXYZgrid, vyXYZgrid, vzXYZgrid are the components
% vXYZgrid is the magnitude
% note: For interpolating the density only: interpolates ones(size(XYZ,1),1,'single')
W = kernelSPH(hLandshoff,'lucy',3); % kernel expression
vXYZ = Xstress.ATOMS{insidewindow,vcoords}; % kernel centers
vXYZmag = sqrt(sum(vXYZ.^2,2)); % velocity magnitude
mbead = 9.04e-12;
Vbead = mbead/1000;
v3XYZgrid = interp3SPHVerlet(XYZ,vXYZ,XYZgrid,VXYZ,W,Vbead);
vxXYZgrid = reshape(v3XYZgrid(:,1),size(Xw)); vxXYZgrid(isnan(vxXYZgrid)) = 0;
vyXYZgrid = reshape(v3XYZgrid(:,2),size(Xw)); vyXYZgrid(isnan(vyXYZgrid)) = 0;
vzXYZgrid = reshape(v3XYZgrid(:,3),size(Xw)); vzXYZgrid(isnan(vzXYZgrid)) = 0;
vXYZgrid = reshape(sqrt(sum(v3XYZgrid.^2,2)),size(Xw));
 
% 3D Figure (for control only, we plot one iso-velocity surface)
% we add the solid phase as a meshed solid
figure, hold on
visocontour = max(vXYZmag)/10;
isosurface(Xw,Yw,Zw,vXYZgrid,visocontour)
DT = delaunayTriangulation(double(Xstress.ATOMS{isolidcontact,{'x','y','z'}}));
K = convexHull(DT);
trisurf(K, DT.Points(:,1), DT.Points(:,2), DT.Points(:,3), 'FaceColor', 'w','Edgecolor','k','FaceAlpha',0.6);
lighting gouraud, camlight left, axis equal, view(3) % shading interp
 
%% cuts and slices (caps are put on the smallest contour)
figure, hold on
vXYZcut = vXYZgrid;
boxcenter = (fluidbox(:,1)+fluidbox(:,2))/2; % boxcenter = (solidbox(:,1)+solidbox(:,2))/2;
[Xwcut,Ywcut,Zwcut] = deal(Xw,Yw,Zw);
vXYZcut(:,yw>boxcenter(2),:) = [];
Xwcut(:,yw>boxcenter(2),:) = [];
Ywcut(:,yw>boxcenter(2),:) = [];
Zwcut(:,yw>boxcenter(2),:) = [];
% vXYZcut(:,:,zw>boxcenter(3),:) = [];
% Xwcut(:,:,zw>boxcenter(3),:) = [];
% Ywcut(:,:,zw>boxcenter(3),:) = [];
% Zwcut(:,:,zw>boxcenter(3),:) = [];
% iso-surface val
colors = parula(256);
visomax = max(vXYZmag);
visocontour = visomax * [1/20,1/10,1/5, 1/4, 1/3, 1/2];
visocolors = interp1(linspace(0,visomax,size(colors,1)),colors,visocontour);
for icontour = 1:length(visocontour)
patch(isosurface(Xwcut,Ywcut,Zwcut,vXYZcut, visocontour(icontour)),...
'FaceColor',visocolors(icontour,:),'EdgeColor','none','FaceAlpha',0.6);
end
p2 = patch(isocaps(Xwcut,Ywcut,Zwcut,vXYZcut, visocontour(1)),'FaceColor','interp','EdgeColor','none','FaceAlpha',0.6);
colormap(colors), camlight('right'), lighting gouraud, view(3), axis equal
trisurf(K, DT.Points(:,1), DT.Points(:,2), DT.Points(:,3), 'FaceColor', 'w','Edgecolor','k','FaceAlpha',0.6);
hs= slice(Xw,Yw,Zw,vXYZgrid,boxcenter(1), ...
[solidbox(2,1) boxcenter(2) solidbox(2,2)], ...
[viewbox(3,1) boxcenter(3)]);
set(hs,'edgecolor','none','facealpha',0.5)
view(115,34)
% add quiver plot
quiver3(Xw(1:5:end,1:5:end,1:1:end),Yw(1:5:end,1:5:end,1:1:end),Zw(1:5:end,1:5:end,1:1:end), ...
vxXYZgrid(1:5:end,1:5:end,1:1:end),vyXYZgrid(1:5:end,1:5:end,1:1:end),vzXYZgrid(1:5:end,1:5:end,1:1:end), ...
'color','k','LineWidth',1)
% add streamlines (note that all arguments required to be double)
[startX,startY,startZ] = meshgrid(double(xw(1)),double(yw(5:5:end-5)),double(zw(1:1:end)));
vstart = interp3(Xw,Yw,Zw,vxXYZgrid,startX,startY,startZ); startX(vstart<0) = double(xw(end));
hsl = streamline(double(Xw),double(Yw),double(Zw),vxXYZgrid,vyXYZgrid,vzXYZgrid,startX,startY,startZ);
set(hsl,'linewidth',2,'color',[0.4375 0.5000 0.5625])
plot3(startX(:),startY(:),startZ(:),'ro','markerfacecolor',[0.4375 0.5000 0.5625])
end
3935 fluid atoms have be found in the cross section around the solid
Calculate the 1 x 341435 pair distances... ... done in 1.09 s Find atoms within cuttoff... ... done in 2.32 s Build the Verlet list... Building the Verlet list [completed: 94.1% | elapsed: 16 s | remaining: 1.01 s] ... done in 17.3 s [min=0, median=66, max=122] Sort the Verlet list... ... done in 1.61 s buildVerletList: all done in 23.5 s for 341435 atoms
INTERP3SPHVERLET interpolates 337500 x 3 grid points with a Verlet list including from 0 to 122 neighbors... ...done in 38.96 s. INTERP3SPHVerlet completed the interpolation of 337500 points with 3935 kernels
 

3.3 Landshoff forces

%% Landshoff forces in the fluid
%
% This section is dedicated to the calculation of Landshoff forces in the fluid medium.
% The computation is confined to a shell or "crown" around the solid particles.
% Here, we build a new Verlet list specific to the fluid particles around solid structures.
%
% * `Xfluidcontact_withneighbors`: Atomic information specifically for the fluid atoms in the vicinity of the solid.
% * `Vfluidcontact_withneighbors`: Verlet list constructed for these neighboring fluid atoms.
% * `Flandshoff`: Landshoff forces calculated based on the new Verlet list.
% * `f`: Magnitude of the Landshoff forces.
%
% Finally, visualization aids are employed to graphically represent these forces and their directions.
 
 
% Extract Information for Fluid Atoms in Contact with Neighbors
% Obtains the atomic information for fluid atoms in the vicinity of solid atoms, as identified earlier.
Xfluidcontact_withneighbors = Xstress.ATOMS(ifluidcontact_withneighbors,:);
 
% Verlet List Construction for Fluid Atoms
% Constructs a Verlet list for these neighboring fluid atoms with a larger cutoff (4*rbead).
Vfluidcontact_withneighbors = buildVerletList(Xfluidcontact_withneighbors,4*rbead);
Calculate the 1921 x 1921 pair distances... ... done in 0.00487 s Find atoms within cuttoff... ... done in 0.00582 s Build the Verlet list... ... done in 0.0473 s [min=9, median=21, max=29] Sort the Verlet list... ... done in 0.00591 s buildVerletList: all done in 0.0767 s for 1921 atoms
 
%%% Compute Landshoff Forces
% Computes Landshoff forces based on the new Verlet list.
% The config structure (configLandshoff) should match the input values of the simulation or equivalent
hLandshoff = 5*rbead; %1.25e-5; % m
configLandshoff = struct( ...
'gradkernel', kernelSPH(hLandshoff,'lucyder',3),...kernel gradient
'h', hLandshoff,...smoothing length (m)
'c0',0.32,...speed of the sound (m/s)
'q1',30,... viscosity coefficient (-)
'rho', 1000, ...density
'm', 9.04e-12 ...
);
Flandshoff = forceLandshoff(Xfluidcontact_withneighbors,[],Vfluidcontact_withneighbors,configLandshoff);
Calculate Landshoff forces between [1921 x 3] atoms... ... done in 0.0539 s
Flandshoff_all = Flandshoff;
 
% Restrict to the Crown/Shell Around the Solid
% Filters out the Landshoff forces, focusing only on the "crown" or shell around the solid atoms.
Flandshoff = double(Flandshoff(ind_withoutneighbors,:)); % we restrict only to the crown/shell
flandshoff = sqrt(sum(Flandshoff.^2,2));
statvec(flandshoff,'Landshoff',sprintf('<-- TIMESTEP: %d',timestepforstress))
Landshoff: 2.5%> 4.96 | 25.0%> 12.42 | 50.0%> 17.58 | 75.0%> 23.4 | 97.5%> 36.87 | 886 values | average = 18.463 <-- TIMESTEP: 912000
 
%%% Store the current result
postdata(end+1).type = 'Landshoff'; % end+1 triggers a new object
postdata(end).timestep = timestepforstress; % we fill it with end
postdata(end).Xfluid = Xstress.ATOMS(ifluidcontact,:);
postdata(end).Xsolid = Xstress.ATOMS(isolidcontact,:);
postdata(end).force = Flandshoff;
 
% loop here to accumulate more postdata

3.4 Hertz Contacts

%%% Compute Hertz contacts
% Extract Information for Fluid Atoms in Contact
Xcontactregion = Xstress.ATOMS(union(ifluidcontact,isolidcontact),:);
% Verlet List Construction for Solid Atoms (only) and considering only Fluid Atoms as neighbors
[Vcontactregion,~,dmincontact] = buildVerletList(Xcontactregion,3*rbead,[],[],[],Xcontactregion.isfluid,Xcontactregion.issolid);
Calculate the 1 x 1466 pair distances... ... done in 0.00145 s Find atoms within cuttoff... ... done in 0.00146 s Build the Verlet list... ... done in 0.00554 s [min=0, median=0, max=11] Sort the Verlet list... ... done in 0.00477 s buildVerletList: all done in 0.0268 s for 1466 atoms
% Compute Hertz Forces
Rfluid = 1.04e-5; % m
Rsolid = 1.56e-5; % m
Rfluid = Rsolid;
configHertz = struct('R',{Rsolid Rfluid},'E',2000);
FHertz = forceHertz(Xcontactregion,Vcontactregion,configHertz);
fhertz = sqrt(sum(FHertz.^2,2));
fhertzcontacts = fhertz(fhertz>0);
statvec(fhertzcontacts,'Hertz',sprintf('<-- TIMESTEP: %d\n\tsubjected to Rsolid=[%0.4g %0.4g] dmin/2=%0.4g',timestepforstress,configHertz(1).R,configHertz(2).R,dmincontact/2))
Hertz: 2.5%> 1.667e-08 | 25.0%> 4.027e-08 | 50.0%> 6.772e-08 | 75.0%> 8.849e-08 | 97.5%> 1.05e-07 | 580 values | average = 6.378e-08 <-- TIMESTEP: 912000 subjected to Rsolid=[1.56e-05 1.56e-05] dmin/2=1.376e-05
 
%%% Store current results
% Store the current result
postdata(end+1).type = 'Hertz'; % end+1 triggers a new object
postdata(end).timestep = timestepforstress; % we fill it with end
postdata(end).Xfluid = [];
postdata(end).Xsolid = Xcontactregion(fhertz>0,:);
postdata(end).force = FHertz(fhertz>0);
 
% loop here to accumulate more postdata
 
 

4. Force Visualization

4.1 Solid visualization using a convex-hull

%% Visualization (bead-based) of Landshoff forces
% Creates a 3D visualization featuring both the solid structure and the fluid atoms under consideration.
% The Landshoff forces are represented as arrows originating from the fluid atoms.
 
% Tesselated Solid Visualization
% Utilizes Delaunay triangulation to represent the solid structure.
DT = delaunayTriangulation(double(Xstress.ATOMS{isolidcontact,{'x','y','z'}}));
K = convexHull(DT);
figure, hold on
trisurf(K, DT.Points(:,1), DT.Points(:,2), DT.Points(:,3), 'FaceColor', 'w','Edgecolor','k','FaceAlpha',0.6);
axis equal; view(3)
 

4.2. Basic visualization of Landshoff forces

%% Basic Visualization of Fluid Atoms and Landshoff Forces
% Landshoff forces are scaled and then visualized as arrows.
% fluid around
% scatter3D(Xfluidcontact_withneighbors{ind_withoutneighbors,{'x','y','z'}},f);
fmedian = median(flandshoff);
fmin = fmedian/50;
fscale = 4*rbead/fmedian;
start = Xfluidcontact_withneighbors{ind_withoutneighbors,{'x','y','z'}};
stop = start + Flandshoff * fscale;
start(flandshoff<fmin,:) = []; % non-significant forces are removed
stop(flandshoff<fmin,:) = [];
ha = arrow(start,stop,'length',4,'BaseAngle',60,'color','r');
Warning: ARROW changed the axis limits on several axes when adding the arrows.
Call ARROW FIXLIMITS to reset them now.
axis equal, view(3)
title('Landshoff Forces around the Tessellated Surface');
xlabel('X'); ylabel('Y'); zlabel('Z')
 

4.3. Advanced grid visualization of Lanshoff forces

The followed approach generalize the quiver plot used for the velocity field
%% Advanced LandShoff Visualization based on Grid Interpolation
coords = {'x','y','z'};
landshofbox = [ min(Xfluidcontact_withneighbors{:,coords})
max(Xfluidcontact_withneighbors{:,coords}) ]';
xlw = linspace(viewbox(1,1),viewbox(1,2),50);
ylw = linspace(viewbox(2,1),viewbox(2,2),50);
zlw = linspace(viewbox(3,1),viewbox(3,2),50);
[Xlw,Ylw,Zlw] = meshgrid(xlw,ylw,zlw);
XYZlgrid = [Xlw(:),Ylw(:),Zlw(:)];
XYZl = Xfluidcontact_withneighbors{:,coords}; % kernel centers
% Interpolation at Grid Points of Flandshoff_all
VXYZl = buildVerletList({XYZlgrid XYZl},1.2*hLandshoff); % special grid syntax
Calculate the 1 x 126921 pair distances... ... done in 0.222 s Find atoms within cuttoff... ... done in 0.396 s Build the Verlet list... ... done in 1.44 s [min=0, median=0, max=68] Sort the Verlet list... ... done in 0.286 s buildVerletList: all done in 2.67 s for 126921 atoms
W = kernelSPH(hLandshoff,'lucy',3); % kernel expression
mbead = 9.04e-12;
Vbead = mbead/1000;
FXYZgrid = interp3SPHVerlet(XYZl,Flandshoff_all,XYZlgrid,VXYZl,W,Vbead);
INTERP3SPHVERLET interpolates 125000 x 3 grid points with a Verlet list including from 0 to 68 neighbors... ...done in 4.137 s. INTERP3SPHVerlet completed the interpolation of 125000 points with 1921 kernels
FXYZgridx = reshape(FXYZgrid(:,1),size(Xlw));
FXYZgridy = reshape(FXYZgrid(:,2),size(Ylw));
FXYZgridz = reshape(FXYZgrid(:,3),size(Zlw));
% Quiver plot to show the forces with an adjusted step
quiver3(...
Xlw(1:2:end,1:2:end,1:2:end), ...
Ylw(1:2:end,1:2:end,1:2:end), ...
Zlw(1:2:end,1:2:end,1:2:end), ...
FXYZgridx(1:2:end,1:2:end,1:2:end), ...
FXYZgridy(1:2:end,1:2:end,1:2:end), ...
FXYZgridz(1:2:end,1:2:end,1:2:end), ...
'color','k','LineWidth',1)

4.4 Basic visualization of Hertz contacts

% Similar Visualization for Hertz Contacts
warning off
figure, hold on
fmedian = median(fhertzcontacts);
fmin = fmedian/50;
fscale = 4*rbead/fmedian;
start = Xcontactregion{fhertz>0,{'x','y','z'}};
stop = start - FHertz(fhertz>0,:) * fscale; % outwards forces (preferred)
start(fhertzcontacts<fmin,:) = []; % non-significant forces are removed
stop(fhertzcontacts<fmin,:) = [];
D = pdist2(DT.Points, start); % Compute pairwise distance
[~, idx] = min(D, [], 2); % Find closest particles on the mesh
trisurf(K, DT.Points(:,1), DT.Points(:,2), DT.Points(:,3), fhertzcontacts(idx),'Edgecolor','k','FaceAlpha',0.6);
axis equal; view(3), camlight('headlight'); camlight('left'); lighting phong; colorbar;
shading interp
ha = arrow(start,stop,'length',4,'BaseAngle',60,'color','r');
axis equal, view(3)
title('Hertz Forces onto the Tessellated Surface');
xlabel('X'); ylabel('Y'); zlabel('Z')
 

4.5 Advanced visualization of Hert contacts

The forces are calculated at the surface of the solid particle defined by a refined and smoothed mesh.
%% Grid Visualization of Hertz Contacts
% === STEP 1/5 === Original Delaunay triangulation and the convex hull
DT = delaunayTriangulation(double(Xstress.ATOMS{isolidcontact, {'x', 'y', 'z'}}));
K = convexHull(DT);
% === STEP 2/5 === refine the initial mesh by adding midpoints
% Extract the convex hull points and faces
hullPoints = DT.Points;
hullFaces = DT.ConnectivityList(K, :);
% Initialize a set to keep track of midpoints to ensure they are unique
midpointSet = zeros(0, 3);
% Calculate midpoints for each edge in each triangle and add to the point list
for faceIdx = 1:size(hullFaces, 1)
face = hullFaces(faceIdx, :);
for i = 1:3
for j = i+1:3
p1 = hullPoints(face(i), :);
p2 = hullPoints(face(j), :);
midpoint = (p1 + p2) / 2;
% Store the midpoint if unique
if isempty(midpointSet) || ~ismember(midpoint, midpointSet, 'rows')
midpointSet = [midpointSet; midpoint]; %#ok<AGROW>
end
end
end
end
% Merge the original points and the new midpoints
newPoints = [hullPoints; midpointSet];
% Re-calculate the Delaunay triangulation and convex hull)
newDT = delaunayTriangulation(newPoints);
newK = convexHull(newDT);
% === STEP 3/5 === Laplacian Smoothing
points = newDT.Points; % === Extract points and faces
faces = newDT.ConnectivityList(newK, :);
n = size(points, 1); % === Initialize new points
newPoints = zeros(size(points));
neighbors = cell(n, 1); % List of neighbors
for faceIdx = 1:size(faces, 1) % === Find the neighbors of each vertex
face = faces(faceIdx, :);
for i = 1:3
vertex = face(i);
vertex_neighbors = face(face ~= vertex);
neighbors{vertex} = unique([neighbors{vertex}; vertex_neighbors(:)]);
end
end
for i = 1:n % === Laplacian smoothing
neighbor_indices = neighbors{i};
if isempty(neighbor_indices) % Keep the point as is if it has no neighbors
newPoints(i, :) = points(i, :);
else % Move the point to the centroid of its neighbors
newPoints(i, :) = mean(points(neighbor_indices, :), 1);
end
end
% Update the Delaunay triangulation with the smoothed points
newDT = delaunayTriangulation(newPoints);
newK = convexHull(newDT);
% === STEP 4/5 === Interpolate the Hertz forces on the triangular mesh
coords = {'x','y','z'};
XYZhtri = newDT.Points;
XYZh = Xcontactregion{:,coords}; % kernel centers
VXYZh = buildVerletList({XYZhtri XYZh},1.2*hLandshoff); % special grid syntax
Calculate the 1 x 2547 pair distances... ... done in 0.00249 s Find atoms within cuttoff... ... done in 0.00441 s Build the Verlet list... ... done in 0.0479 s [min=0, median=0, max=70] Sort the Verlet list... ... done in 0.00787 s buildVerletList: all done in 0.0738 s for 2547 atoms
W = kernelSPH(hLandshoff,'lucy',3); % kernel expression
mbead = 9.04e-12;
Vbead = mbead/1000;
FXYZtri = interp3SPHVerlet(XYZh,FHertz,XYZhtri,VXYZh,W,Vbead);
INTERP3SPHVERLET interpolates 1081 x 3 grid points with a Verlet list including from 0 to 70 neighbors... ...done in 0.1381 s. INTERP3SPHVerlet completed the interpolation of 1081 points with 1466 kernels
FXYZtrix = reshape(FXYZgrid(:,1),size(Xlw));
FXYZtriy = reshape(FXYZgrid(:,2),size(Ylw));
FXYZtriz = reshape(FXYZgrid(:,3),size(Zlw));
% === STEP 5/5 === Extract tagential forces
% Calculate face normals and centroids
points = newDT.Points;
faces = newK;
v1 = points(faces(:, 1), :) - points(faces(:, 2), :);
v2 = points(faces(:, 1), :) - points(faces(:, 3), :);
faceNormals = cross(v1, v2, 2);
faceNormals = faceNormals ./ sqrt(sum(faceNormals.^2, 2));
centroids = mean(reshape(points(faces, :), size(faces, 1), 3, 3), 3);
% Interpolate force at each centroid using scatteredInterpolant for each component
FInterp_x = scatteredInterpolant(XYZhtri, double(FXYZtri(:,1)), 'linear', 'nearest');
FInterp_y = scatteredInterpolant(XYZhtri, double(FXYZtri(:,2)), 'linear', 'nearest');
FInterp_z = scatteredInterpolant(XYZhtri, double(FXYZtri(:,3)), 'linear', 'nearest');
FXYZtri_at_centroids = [FInterp_x(centroids), FInterp_y(centroids), FInterp_z(centroids)];
% Calculate normal and tangential components of the force at each face centroid
normalComponent = dot(FXYZtri_at_centroids, faceNormals, 2);
normalForce = repmat(normalComponent, 1, 3) .* faceNormals;
tangentialForce = FXYZtri_at_centroids - normalForce;
% Calculate the magnitude of the tangential force
tangentialMagnitude = sqrt(sum(tangentialForce.^2, 2));
 
 
% === FINAL PLOTS ===
figure, hold on
trisurfHandle = trisurf(newK, newDT.Points(:, 1), newDT.Points(:, 2), newDT.Points(:, 3), 'Edgecolor', 'k', 'FaceAlpha', 0.6);
set(trisurfHandle, 'FaceVertexCData', tangentialMagnitude, 'FaceColor', 'flat');
colorbar;
% Quiver plot to show the forces with an adjusted step
quiver3(...
XYZhtri(:,1), ...
XYZhtri(:,2), ...
XYZhtri(:,3), ...
-FXYZtri(:,1), ...
-FXYZtri(:,2), ...
-FXYZtri(:,3), ...
2, ...scale
'color','k','LineWidth',3)
axis equal; view(3), camlight('headlight'); camlight('left'); lighting phong; colorbar;
title('Hertz Forces onto the Tessellated Surface');
xlabel('X'); ylabel('Y'); zlabel('Z')

5. Stress calculation and visualization

The difference between force and stress visualization is that you need to divide force by the surface area, where the force applies.
The current version use of example2 demonstrate the strategy without combining the method with the most advanced interpolation of forces on structured or unstructured mesh. For production, it is recommended to use forces which have be smoothed (interpolated) over a grid.

5.1 Template strategy

In the absence of grid, the results are not as smooth.
%% Project Landshoff Forces onto the Tessellated Surface
% Allocate space for the projected forces.
FprojectedN = zeros(size(Flandshoff));
FprojectedT = zeros(size(Flandshoff));
% Extract vertex coordinates from Delaunay Triangulation object.
vertices = DT.Points;
% Loop through all fluid atoms in contact with the solid.
for i = 1:length(ind_withoutneighbors)
% Get the position of the current fluid atom.
fluidPos = Xfluidcontact_withneighbors{ind_withoutneighbors(i),{'x','y','z'}};
% Find the closest vertex on the tessellated surface.
dists = sum((vertices - fluidPos).^2, 2);
[~, closestVertexIdx] = min(dists);
closestVertex = vertices(closestVertexIdx, :);
% Calculate the vector from the fluid atom to the closest surface vertex.
surfaceToFluidVec = fluidPos - closestVertex;
% Normalize the vector.
surfaceToFluidVec = surfaceToFluidVec / norm(surfaceToFluidVec);
% Get the Landshoff force acting on the fluid atom.
landshoffForce = Flandshoff(i, :);
% Project the Landshoff force onto the normal.
% Here, surfaceToFluidVec serves as the approximation of the outward normal at the closest vertex.
projectedForceN = dot(landshoffForce, surfaceToFluidVec) * surfaceToFluidVec;
% Store the projected force.
FprojectedN(i, :) = projectedForceN;
FprojectedT(i, :) = landshoffForce-projectedForceN;
end
 
% At this point, Fprojected contains the Landshoff forces projected onto the surface normal.
% Interpolation and Stress Magnitude Visualization
% Get faces and vertices from the tessellation.
faces = K; % Faces are given by the convex hull indices
vertices = DT.Points; % Vertex coordinates
 
% Loop over all faces to interpolate force and stress.
% Compute Barycentric Interpolation of Forces
% Initialize an array to store the interpolated forces and stresses
FNinterp_face = zeros(size(faces, 1), 3);
FTinterp_face = zeros(size(faces, 1), 3);
SNinterp_face = zeros(size(faces, 1), 3); % Normal stress per unit area
STinterp_face = zeros(size(faces, 1), 3); % Tangential stress per unit area
 

5.1.1 Stress at faces

% Loop over all faces.
for i = 1:size(faces, 1)
% Vertices of the triangle.
vertex_indices = faces(i, :);
A = vertices(vertex_indices(1), :);
B = vertices(vertex_indices(2), :);
C = vertices(vertex_indices(3), :);
% Surface area of the triangle using Heron's formula.
AB = norm(A - B);
AC = norm(A - C);
BC = norm(B - C);
s = (AB + AC + BC) / 2;
area_triangle = sqrt(s * (s - AB) * (s - AC) * (s - BC));
% Normal to faces
N = cross(B - A, C - A);
N = N / norm(N); % Normalizing the normal vector
% Forces at the vertices.
F_A = FprojectedN(vertex_indices(1), :) + FprojectedT(vertex_indices(1), :);
F_B = FprojectedN(vertex_indices(2), :) + FprojectedT(vertex_indices(2), :);
F_C = FprojectedN(vertex_indices(3), :) + FprojectedT(vertex_indices(3), :);
% The centroid (P) of the triangle for interpolation.
P = (A + B + C) / 3;
% Barycentric weights.
ABC_inv = inv([A; B; C]);
w = ABC_inv * P';
% Interpolate the force at the centroid.
Finterp = w(1) * F_A + w(2) * F_B + w(3) * F_C;
% Decompose the force into normal and tangential components
F_N = dot(Finterp, N) * N; % only the normal force is trivial to compute
F_T = Finterp - F_N; % by subtraction we get the tangential
% Store interpolated force.
Finterp_face(i, :) = Finterp;
% Compute normal and tangential stress per unit area (Force/Area)
SNinterp_face(i, :) = F_N / area_triangle;
STinterp_face(i, :) = F_T / area_triangle;
end
% Compute the magnitude of interpolated stress.
face_force_magnitude = sqrt(sum(Finterp_face.^2, 2));
face_stressN_magnitude = sqrt(sum(SNinterp_face.^2, 2));
face_stressT_magnitude = sqrt(sum(STinterp_face.^2, 2));
 
 

5.1.2. STress at vertices

% Add the smoothed force at vertices
% In computational geometry, the "valence" of a vertex refers to the number of edges
% (or equivalently, faces for a triangular mesh) incident to it. Valence can serve as
% an important metric in mesh quality and adaptivity considerations.
% In our case, the valence is expected high (up to 11,12 faces per vertex), then poor mesh
uniq_vertex_indices = unique(faces(:)); % Find unique vertices used in faces
uniq_vertices = vertices(uniq_vertex_indices, :); % Create a new vertices array based on these unique indices
mapping = NaN(size(vertices, 1), 1); % Create a mapping from old vertex indices to new ones
mapping(uniq_vertex_indices) = 1:length(uniq_vertex_indices);
uniq_faces = mapping(faces); % Update the faces array to reflect new vertex indices
% Initialize interpolated forces at vertices
Finterp_vertex = zeros(size(uniq_vertices, 1), 3);
SNinterp_vertex = zeros(size(uniq_vertices, 1), 3);
STinterp_vertex = zeros(size(uniq_vertices, 1), 3);
% Create a mapping from vertices to faces
vertex_to_faces = cell(size(uniq_vertices, 1), 1);
for i = 1:size(uniq_faces, 1)
for j = 1:3 % for each Vertex
vertex_to_faces{uniq_faces(i, j)} = [vertex_to_faces{uniq_faces(i, j)}; i]; % we add the ith face
end
end
% Interpolate forces at each vertex
for i = 1:size(uniq_vertices, 1)
adjacent_faces = vertex_to_faces{i};
if ~isempty(adjacent_faces)
Finterp_vertex(i, :) = mean(Finterp_face(adjacent_faces, :), 1);
SNinterp_vertex(i, :) = mean(SNinterp_face(adjacent_faces, :), 1);
STinterp_vertex(i, :) = mean(STinterp_face(adjacent_faces, :), 1);
end
end
vertex_force_magnitude = sqrt(sum(Finterp_vertex.^2, 2));
vertex_stressN_magnitude = sqrt(sum(SNinterp_vertex.^2, 2));
vertex_stressT_magnitude = sqrt(sum(STinterp_vertex.^2, 2));
 

5.2. Visualization

5.2.1. Face-based visualization

% Visualization (face-based, less smooth)
figure;
trisurf(faces, vertices(:, 1), vertices(:, 2), vertices(:, 3), face_force_magnitude, 'EdgeColor', 'none');
axis equal; view(3);
colorbar;
title('Interpolated Force Magnitude on Tessellated Surface (face-based)');
xlabel('X'); ylabel('Y'); zlabel('Z')
 

5.2.2. Vertex-based visualization

% Visualization (vertex-based, more smooth)
figure;
trisurf(uniq_faces, uniq_vertices(:, 1), uniq_vertices(:, 2), uniq_vertices(:, 3), vertex_force_magnitude, 'EdgeColor', 'none');
axis equal; view(3);
colorbar;
title('Interpolated Force Magnitude on Tessellated Surface (vertex-based)');
xlabel('X'); ylabel('Y'); zlabel('Z')
 

5.2.3 Slighly optimized visualization

% Optimized visualization using 'patch' for stresses
% only vertex_stressT_magnitude is meaningful
figure;
patch_data.Vertices = uniq_vertices; % Prepare data for 'patch' function
patch_data.Faces = uniq_faces;
patch_data.FaceVertexCData = vertex_stressT_magnitude;
patch_data.FaceColor = 'interp';
patch_data.EdgeColor = 'none';
% Calculate vertex normals for smooth shading
patch_data.VertexNormals = -vertexNormal(triangulation(uniq_faces, uniq_vertices));
% Create the surface plot
p = patch(patch_data);
% Set view and lighting
axis equal;
view(3); camlight('headlight'); camlight('left'); lighting phong; colorbar;
xlabel('X'); ylabel('Y'); zlabel('Z');
title('Tangential Landshoff Stress Magnitude onto Solid Surface');
% Option to remove grid for a cleaner look
% grid off;