adds output of IB state (posittion, velocities, angles)

Author: mrvandenboom

docs/documentation/case.md

@@ -615,7 +615,7 @@ To restart the simulation from $k$-th time step, see @ref running "Restarting Ca
 | `alpha_wrt(i)`          | Logical | Add the volume fraction of fluid $i$ to the database	|
 | `gamma_wrt`             | Logical | Add the specific heat ratio function to the database	|
 | `heat_ratio_wrt`        | Logical | Add the specific heat ratio to the database	|
-| `ib_force_wrt`          | Logical | Write IB force & torque datafile at each time step |
+| `ib_state_wrt`          | Logical | Write IB state and loads to a datafile at each time step |
 | `pi_inf_wrt`            | Logical | Add the liquid stiffness function to the database |
 | `pres_inf_wrt`          | Logical | Add the liquid stiffness to the formatted database	 |
 | `c_wrt`                 | Logical | Add the sound speed to the database	 |
@@ -676,7 +676,7 @@ If `file_per_process` is true, then pre_process, simulation, and post_process mu
 
 - `probe_wrt` activates the output of state variables at coordinates specified by `probe(i)%[x;y,z]`.
 
-- `ib_force_wrt` activates the output of data specified by patch_ib(i)%force(:) and patch_ib(i)%torque(:) into a single binary datafile for all IBs at all timesteps. During post_processing, this file is converted into separate force/torque time histories for each IB.
+- `ib_state_wrt` activates the output of data specified by patch_ib(i)%force(:) (and torque, vel, angular_vel, angles, [x,y,z]_centroid) into a single binary datafile for all IBs at all timesteps. During post_processing, this file is converted into separate time histories for each IB.
 
 - `output_partial_domain` activates the output of part of the domain specified by `[x,y,z]_output%%beg` and `[x,y,z]_output%%end`.
 This is useful for large domains where only a portion of the domain is of interest.

src/post_process/m_data_output.fpp

@@ -30,7 +30,7 @@ module m_data_output
  s_write_variable_to_formatted_database_file, &
  s_write_lag_bubbles_results_to_text, &
  s_write_lag_bubbles_to_formatted_database_file, &
- s_write_ib_force_files, &
+ s_write_ib_state_files, &
  s_write_intf_data_file, &
  s_write_energy_data_file, &
  s_close_formatted_database_file, &
@@ -1500,17 +1500,19 @@ contains
 
     end subroutine s_write_lag_variable_to_formatted_database_file
 
-    impure subroutine s_write_ib_force_files()
+    impure subroutine s_write_ib_state_files()
 
         character(len=len_trim(case_dir) + 4*name_len) :: in_file, out_file, file_loc
-        integer :: iu_in, ios, i, rec_step, rec_id
+        integer :: iu_in, ios, i, rec_id
         integer, allocatable, dimension(:) :: iu_out
         real(wp) :: rec_time
         real(wp), dimension(3) :: rec_force, rec_torque
+        real(wp), dimension(3) :: rec_vel, rec_angular_vel
+        real(wp), dimension(3) :: rec_angles, rec_centroid
 
         file_loc = trim(case_dir)//'/D'
 
-        in_file = trim(file_loc)//'/ib_force.dat'
+        in_file = trim(file_loc)//'/ib_state.dat'
         open (newunit=iu_in, file=trim(in_file), form='unformatted', access='stream', &
               status='old', action='read', iostat=ios)
         if (ios /= 0) then
@@ -1531,16 +1533,24 @@ contains
                 deallocate (iu_out)
                 return
             end if
-            write (iu_out(i), '(A)') 'mytime fx fy fz Tau_x Tau_y Tau_z'
+            write (iu_out(i), '(A)') &
+                'mytime fx fy fz Tau_x Tau_y Tau_z vx vy vz omega_x omega_y omega_z angle_x angle_y angle_z x_c y_c z_c'
         end do
 
         do
-            read (iu_in, iostat=ios) rec_step, rec_time, rec_id, rec_force, rec_torque
-            if (ios /= 0) exit  ! EOF (<0) or read error (>0)
+            read (iu_in, iostat=ios) rec_time, rec_id, &
+                rec_force, rec_torque, rec_vel, rec_angular_vel, rec_angles, &
+                rec_centroid(1), rec_centroid(2), rec_centroid(3)
+            if (ios /= 0) exit
+
             if (rec_id >= 1 .and. rec_id <= num_ibs) then
-                write (iu_out(rec_id), '(7(ES24.16E3,1X))') rec_time, &
+                write (iu_out(rec_id), '(19(ES24.16E3,1X))') rec_time, &
                     rec_force(1), rec_force(2), rec_force(3), &
-                    rec_torque(1), rec_torque(2), rec_torque(3)
+                    rec_torque(1), rec_torque(2), rec_torque(3), &
+                    rec_vel(1), rec_vel(2), rec_vel(3), &
+                    rec_angular_vel(1), rec_angular_vel(2), rec_angular_vel(3), &
+                    rec_angles(1), rec_angles(2), rec_angles(3), &
+                    rec_centroid(1), rec_centroid(2), rec_centroid(3)
             end if
         end do
 
@@ -1550,7 +1560,7 @@ contains
         end do
         deallocate (iu_out)
 
-    end subroutine s_write_ib_force_files
+    end subroutine s_write_ib_state_files
 
     !> @brief Extract the volume-fraction interface contour from primitive fields and write the coordinates to the interface data file.
     impure subroutine s_write_intf_data_file(q_prim_vf)

src/post_process/m_global_parameters.fpp

@@ -261,7 +261,7 @@ module m_global_parameters
     logical :: schlieren_wrt
     logical :: cf_wrt
     logical :: ib
-    logical :: ib_force_wrt
+    logical :: ib_state_wrt
     logical :: chem_wrt_Y(1:num_species)
     logical :: chem_wrt_T
     logical :: lag_header
@@ -495,7 +495,7 @@ contains
         sim_data = .false.
         cf_wrt = .false.
         ib = .false.
-        ib_force_wrt = .false.
+        ib_state_wrt = .false.
         lag_txt_wrt = .false.
         lag_header = .true.
         lag_db_wrt = .false.

src/post_process/m_mpi_proxy.fpp

@@ -101,7 +101,7 @@ contains
             & 'adv_n', 'ib', 'cfl_adap_dt', 'cfl_const_dt', 'cfl_dt',          &
             & 'surface_tension', 'hyperelasticity', 'bubbles_lagrange',        &
             & 'output_partial_domain', 'relativity', 'cont_damage', 'bc_io',   &
-            & 'down_sample','fft_wrt', 'hyper_cleaning', 'ib_force_wrt']
+            & 'down_sample','fft_wrt', 'hyper_cleaning', 'ib_state_wrt']
             call MPI_BCAST(${VAR}$, 1, MPI_LOGICAL, 0, MPI_COMM_WORLD, ierr)
         #:endfor
 

src/post_process/m_start_up.fpp

@@ -116,7 +116,7 @@ contains
             lag_rad_wrt, lag_rvel_wrt, lag_r0_wrt, lag_rmax_wrt, &
             lag_rmin_wrt, lag_dphidt_wrt, lag_pres_wrt, lag_mv_wrt, &
             lag_mg_wrt, lag_betaT_wrt, lag_betaC_wrt, &
-            alpha_rho_e_wrt, ib_force_wrt
+            alpha_rho_e_wrt, ib_state_wrt
 
         ! Inquiring the status of the post_process.inp file
         file_loc = 'post_process.inp'

src/post_process/p_main.fpp

@@ -93,8 +93,8 @@ program p_main
     end do
     ! END: Time-Marching Loop
 
-    if (proc_rank == 0 .and. ib_force_wrt) then
-        call s_write_ib_force_files()
+    if (proc_rank == 0 .and. ib_state_wrt) then
+        call s_write_ib_state_files()
     end if
 
     close (11)

src/simulation/m_data_output.fpp

@@ -37,19 +37,19 @@ module m_data_output
               s_open_run_time_information_file, &
               s_open_com_files, &
               s_open_probe_files, &
-              s_open_ib_force_file, &
+              s_open_ib_state_file, &
               s_write_run_time_information, &
               s_write_data_files, &
               s_write_serial_data_files, &
               s_write_parallel_data_files, &
               s_write_ib_data_file, &
               s_write_com_files, &
               s_write_probe_files, &
-              s_write_ib_force_file, &
+              s_write_ib_state_file, &
               s_close_run_time_information_file, &
               s_close_com_files, &
               s_close_probe_files, &
-              s_close_ib_force_file, &
+              s_close_ib_state_file, &
               s_finalize_data_output_module
 
     real(wp), allocatable, dimension(:, :, :) :: icfl_sf  !< ICFL stability criterion
@@ -257,17 +257,17 @@ contains
 
     end subroutine s_open_probe_files
 
-    impure subroutine s_open_ib_force_file
+    impure subroutine s_open_ib_state_file
         character(LEN=path_len + 2*name_len) :: file_loc
 
-        write (file_loc, '(A)') 'ib_force.dat'
+        write (file_loc, '(A)') 'ib_state.dat'
         file_loc = trim(case_dir)//'/D/'//trim(file_loc)
         open (92, FILE=trim(file_loc), &
               FORM='unformatted', &
               ACCESS='stream', &
               STATUS='replace', &
               POSITION='append')
-    end subroutine s_open_ib_force_file
+    end subroutine s_open_ib_state_file
 
     !>  The goal of the procedure is to output to the run-time
         !!      information file the stability criteria extrema in the
@@ -1153,20 +1153,27 @@ contains
 
     end subroutine s_write_ib_data_file
 
-    !> @brief Writes IB force/torque records to D/ib_force.dat on rank 0.
-    impure subroutine s_write_ib_force_file(time_step)
+    !> @brief Writes IB state records to D/ib_state.dat on rank 0.
+    impure subroutine s_write_ib_state_file()
 
-        integer, intent(in) :: time_step
         if (proc_rank == 0) then
             block
                 integer :: i
                 do i = 1, num_ibs
-                    write (92) time_step, mytime, i, patch_ib(i)%force, patch_ib(i)%torque
+                    write (92) mytime, i, &
+                        patch_ib(i)%force, &
+                        patch_ib(i)%torque, &
+                        patch_ib(i)%vel, &
+                        patch_ib(i)%angular_vel, &
+                        patch_ib(i)%angles, &
+                        patch_ib(i)%x_centroid, &
+                        patch_ib(i)%y_centroid, &
+                        patch_ib(i)%z_centroid
                 end do
             end block
         end if
 
-    end subroutine s_write_ib_force_file
+    end subroutine s_write_ib_state_file
 
     !>  This writes a formatted data file where the root processor
     !!      can write out the CoM information
@@ -1937,11 +1944,11 @@ contains
 
     end subroutine s_close_probe_files
 
-    impure subroutine s_close_ib_force_file
+    impure subroutine s_close_ib_state_file
 
         close (92)
 
-    end subroutine s_close_ib_force_file
+    end subroutine s_close_ib_state_file
 
     !>  The computation of parameters, the allocation of memory,
         !!      the association of pointers and/or the execution of any

src/simulation/m_global_parameters.fpp

@@ -396,7 +396,7 @@ module m_global_parameters
     !> @{
     logical :: ib
     integer :: num_ibs
-    logical :: ib_force_wrt
+    logical :: ib_state_wrt
 
     type(ib_patch_parameters), dimension(num_patches_max) :: patch_ib
     type(vec3_dt), allocatable, dimension(:) :: airfoil_grid_u, airfoil_grid_l
@@ -716,7 +716,7 @@ contains
         ! Immersed Boundaries
         ib = .false.
         num_ibs = dflt_int
-        ib_force_wrt = .false.
+        ib_state_wrt = .false.
 
         ! Bubble modeling
         bubbles_euler = .false.

src/simulation/m_mpi_proxy.fpp

@@ -114,7 +114,7 @@ contains
             & 'cfl_adap_dt', 'cfl_const_dt', 'cfl_dt', 'surface_tension',       &
             & 'shear_stress', 'bulk_stress', 'bubbles_lagrange',                &
             & 'hyperelasticity', 'down_sample', 'int_comp','fft_wrt', &
-            & 'hyper_cleaning', 'ib_force_wrt']
+            & 'hyper_cleaning', 'ib_state_wrt']
             call MPI_BCAST(${VAR}$, 1, MPI_LOGICAL, 0, MPI_COMM_WORLD, ierr)
         #:endfor
 

src/simulation/m_start_up.fpp

@@ -144,7 +144,7 @@ contains
             x_domain, y_domain, z_domain, &
             hypoelasticity, &
             ib, num_ibs, patch_ib, &
-            ib_force_wrt, &
+            ib_state_wrt, &
             fluid_pp, bub_pp, probe_wrt, prim_vars_wrt, &
             fd_order, probe, num_probes, t_step_old, &
             alt_soundspeed, mixture_err, weno_Re_flux, &