Improved convergence of Fraggle by reducing minimum velocity used to …

…compute available kinetic energy to below the escape velocity on each successive try
MintonGroup · Feb 10, 2023 · 6be045f · 6be045f
1 parent ae92f15
commit 6be045f
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 3 deletions.
diff --git a/examples/Fragmentation/Fragmentation_Movie.py b/examples/Fragmentation/Fragmentation_Movie.py
@@ -59,9 +59,9 @@
                }
 
 vel_vectors = {"disruption_headon"         : [np.array([ 0.00,  6.280005, 0.0]),
-                                              np.array([ 0.00,  4.28,     0.0])],
+                                              np.array([ 0.00,  3.28,     0.0])],
                "disruption_off_axis"       : [np.array([ 0.00,  6.280005, 0.0]),
-                                              np.array([ 0.05,  4.28,     0.0])],
+                                              np.array([ 0.05,  3.28,     0.0])],
                "supercatastrophic_headon":   [np.array([ 0.00,  6.28,     0.0]),
                                               np.array([ 0.00,  4.28,     0.0])],
                "supercatastrophic_off_axis": [np.array([ 0.00,  6.28,     0.0]),

diff --git a/src/fraggle/fraggle_generate.f90 b/src/fraggle/fraggle_generate.f90
@@ -717,7 +717,7 @@ module subroutine fraggle_generate_vel_vec(collider, nbody_system, param, lfailu
                   ! Remove a constant amount of velocity from the bodies so we don't shift the center of mass and screw up the momentum 
                   ke_avail = 0.0_DP
                   do i = fragments%nbody, 1, -1
-                     ke_avail = ke_avail + 0.5_DP * fragments%mass(i) * max(fragments%vmag(i) - vesc,0.0_DP)**2
+                     ke_avail = ke_avail + 0.5_DP * fragments%mass(i) * max(fragments%vmag(i) - vesc/try,0.0_DP)**2
                   end do
 
                   ke_remove = min(E_residual, ke_avail)