Ready to test with a rotated square

Author: danieljvickers

src/simulation/m_ibm.fpp

@@ -962,6 +962,11 @@ contains
         ! Clears the existing immersed boundary indices
         ib_markers%sf = 0
 
+        ! recalulcate the rotation matrix based upon the new angles
+        do i = 1, num_ibs
+          if patch_ib(i)%moving_ibm call s_update_ib_rotation_matrix(i)
+        end if
+
         ! recompute the new ib_patch locations and broadcast them.
         call s_apply_ib_patches(ib_markers%sf(0:m, 0:n, 0:p), levelset, levelset_norm)
         call s_populate_ib_buffers() ! transmits the new IB markers via MPI

src/simulation/m_time_steppers.fpp

@@ -550,9 +550,14 @@ contains
                 ! start by using euler's method naiively, but eventually incorporate more sophistocation
                 if (patch_ib(i)%moving_ibm == 1) then
                     do j = 1, 3
-                        patch_ib(i)%vel(j) = patch_ib(i)%vel(j) + 0.0*dt ! TODO :: ADD EXTERNAL FORCES HERE
+                        patch_ib(i)%vel(j) = patch_ib(i)%vel(j) + 0.0*dt 
+                        patch_ib(i)%angular_vel(j) = patch_ib(i)%angular_vel(j) + 0.0*dt
+
+                        ! Update the angle of the IB
+                        patch_ib(i)%angles(j) = patch_ib(i)%angles(j) + patch_ib(i)%angular_vel(j)*dt
                     end do
 
+                    ! Update the position of the IB
                     patch_ib(i)%x_centroid = patch_ib(i)%x_centroid + patch_ib(i)%vel(1)*dt
                     patch_ib(i)%y_centroid = patch_ib(i)%y_centroid + patch_ib(i)%vel(2)*dt
                     patch_ib(i)%z_centroid = patch_ib(i)%z_centroid + patch_ib(i)%vel(3)*dt
@@ -709,8 +714,15 @@ contains
                 ! start by using euler's method naiively, but eventually incorporate more sophistocation
                 if (patch_ib(i)%moving_ibm == 1) then
                     patch_ib(i)%step_vel = patch_ib(i)%vel
+                    patch_ib(i)%step_angular_vel = patch_ib(i)%angular_vel
+                    patch_ib(i)%step_angles = patch_ib(i)%angles
+
                     do j = 1, 3
-                        patch_ib(i)%vel(j) = patch_ib(i)%step_vel(j) + 0.0*dt ! TODO :: ADD EXTERNAL FORCES HERE
+                        patch_ib(i)%vel(j) = patch_ib(i)%step_vel(j) + 0.0*dt 
+                        patch_ib(i)%angular_vel(j) = patch_ib(i)%step_angular_vel(j) + 0.0*dt
+
+                        ! Update the angle of the IB
+                        patch_ib(i)%angles(j) = patch_ib(i)%step_angles(j) + patch_ib(i)%angular_vel(j)*dt
                     end do
 
                     patch_ib(i)%step_x_centroid = patch_ib(i)%x_centroid
@@ -824,6 +836,9 @@ contains
                 if (patch_ib(i)%moving_ibm == 1) then
                     do j = 1, 3
                         patch_ib(i)%vel(j) = (patch_ib(i)%vel(j) + patch_ib(i)%step_vel(j) + 0.0*dt)/2._wp
+                        patch_ib(i)%angular_vel(j) = (patch_ib(i)%angular_vel(j) + patch_ib(i)%angular_step_vel(j) + 0.0*dt)/2._wp
+
+                        patch_ib(i)%angles(j) = (patch_ib(i)%angles(j) + patch_ib(i)%step_angles(j) + patch_ib(i)%step_angular_vel(j)*dt)/2._wp
                     end do
 
                     patch_ib(i)%x_centroid = (patch_ib(i)%x_centroid + patch_ib(i)%step_x_centroid + patch_ib(i)%step_vel(1)*dt)/2._wp
@@ -981,8 +996,15 @@ contains
               do i = 1, num_ibs
                 if (patch_ib(i)%moving_ibm == 1) then
                     patch_ib(i)%step_vel = patch_ib(i)%vel
+                    patch_ib(i)%step_angular_vel = patch_ib(i)%angular_vel
+                    patch_ib(i)%step_angles = patch_ib(i)%angles
+
                     do j = 1, 3
                         patch_ib(i)%vel(j) = patch_ib(i)%step_vel(j) + 0.0*dt
+                        patch_ib(i)%angular_vel(j) = patch_ib(i)%step_angular_vel(j) + 0.0*dt
+
+                        ! Update the angle of the IB
+                        patch_ib(i)%angles(j) = patch_ib(i)%step_angles(j) + patch_ib(i)%angular_vel(j)*dt
                     end do
 
                     patch_ib(i)%step_x_centroid = patch_ib(i)%x_centroid
@@ -1095,6 +1117,10 @@ contains
                 if (patch_ib(i)%moving_ibm == 1) then
                     do j = 1, 3
                         patch_ib(i)%vel(j) = (patch_ib(i)%vel(j) + 3._wp*patch_ib(i)%step_vel(j) + 0.0*dt)/4._wp
+                        patch_ib(i)%angular_vel(j) = (patch_ib(i)%angular_vel(j) + 3._wp*patch_ib(i)%step_angular_vel(j) + 0.0*dt)/4._wp
+
+                        patch_ib(i)%angles(j) = (patch_ib(i)%angles(j) + 3._wp*patch_ib(i)%step_angles(j) + patch_ib(i)%angular_vel(j)*dt)/4._wp
+
                     end do
 
                     patch_ib(i)%x_centroid = (patch_ib(i)%x_centroid + 3._wp*patch_ib(i)%step_x_centroid + patch_ib(i)%vel(1)*dt)/4._wp
@@ -1207,6 +1233,9 @@ contains
                 if (patch_ib(i)%moving_ibm == 1) then
                     do j = 1, 3
                         patch_ib(i)%vel(j) = (2._wp*patch_ib(i)%vel(j) + patch_ib(i)%step_vel(j) + 0.0*dt)/3._wp
+                        patch_ib(i)%angular_vel(j) = (2._wp*patch_ib(i)%angular_vel(j) + patch_ib(i)%step_angular_vel(j) + 0.0*dt)/3._wp
+
+                        patch_ib(i)%angles(j) = (2._wp*patch_ib(i)%angles(j) + patch_ib(i)%step_angles(j) + patch_ib(i)%angular_vel(j)*dt)/3._wp
                     end do
 
                     patch_ib(i)%x_centroid = (2._wp*patch_ib(i)%x_centroid + patch_ib(i)%step_x_centroid + 2._wp*patch_ib(i)%vel(1)*dt)/3._wp