search

luchete80 / WeldForm

Weakly-Compressible Smoothed Particle Hydrodynamics Parallel Solver for Elasto Plastic and thermal coupled Solid Mechanics
GNU General Public License v3.0
10 stars 5 forks source link

Contact Neighbours and Surface Calculation times are excessive #177

Open luchete80 opened 2 years ago

luchete80 commented 2 years ago

Output No. 50 at 0.00480005 has been generated Current Time Step = 8.64393e-08 Total CPU time: 8415.86 , Nb: 1272.06, StAcc: 0.007, PrAcc: 44.138Ls Acc: 1801.86Contact: 4828.14Msh: 4.715, BC: 9.35, mv: 111.913 Max plastic strain: 0.0607853in particle20936 Steps count in this interval: 1153Total Step count44197 Total Nb search time in this interval: 36.635Average Nb search time in this interval: 0.0317736 Avg Neighbour Count45 Max Contact Force: 148128

luchete80 commented 2 years ago

Save Nb Data is called twice 

            if (contact) {
                SaveNeighbourData();                //Necesary to calulate surface! Using Particle->Nb (count), could be included in search
                CalculateSurface(1);                //After Nb search           
                ContactNbSearch();
                SaveNeighbourData();    //Again Save Nb data
            }//contact
luchete80 commented 2 years ago

SaveNeighbourData is not parallelized

luchete80 commented 2 years ago

Original (dx = 0.015)

Output No. 22 at 0.00200011 has been generated Current Time Step = 1.10891e-07 Total CPU time: 671.437 , Nb: 52.151, StAcc: 0, PrAcc: 8.172Ls Acc: 254.489Contact: 286.16Msh: 0.801, BC: 1.936, mv: 16.377 Max plastic strain: 0.00633061in particle11681 Steps count in this interval: 904Total Step count11651 Total Nb search time in this interval: 10.248Average Nb search time in this interval: 0.0113363 Avg Neighbour Count39 Max Contact Force: 159616

New (with parallel contact and not saving Contact Neighbours)

Output No. 22 at 0.00200004 has been generated Current Time Step = 1.11872e-07 Total CPU time: 675.29 , Nb: 51.245, StAcc: 0.004, PrAcc: 7.721Ls Acc: 259.951Contact: 286.048Msh: 0.751, BC: 1.65, mv: 15.02 Max plastic strain: 0.00626818in particle11653 Steps count in this interval: 898Total Step count11649 Total Nb search time in this interval: 9.929Average Nb search time in this interval: 0.0110568 Avg Neighbour Count38 Max Contact Force: 160121