Dem volume insertion initial velocities (chaos-polymtl#1080)

Follow-up to PR chaos-polymtl#1078 ... Former-commit-id: 13e9dcb
M-Badri · Mar 29, 2024 · b017a0f · b017a0f
1 parent 5ec8f67
commit b017a0f
Show file tree

Hide file tree

Showing 3 changed files with 5 additions and 5 deletions.
diff --git a/include/core/parameters_lagrangian.h b/include/core/parameters_lagrangian.h
@@ -187,7 +187,7 @@ namespace Parameters
       Point<3> insertion_box_point_1, insertion_box_point_2;
 
       // Insertion initial velocity conditions
-      Tensor<1, 3> initial_vel, initial_ang_vel;
+      Tensor<1, 3> initial_vel, initial_omega;
 
       // Insertion distance threshold
       double distance_threshold;

diff --git a/source/core/parameters_lagrangian.cc b/source/core/parameters_lagrangian.cc
@@ -444,7 +444,7 @@ namespace Parameters
         seed_for_insertion       = prm.get_integer("insertion prn seed");
 
         initial_vel = entry_string_to_tensor3(prm, "initial velocity");
-        initial_ang_vel =
+        initial_omega =
           entry_string_to_tensor3(prm, "initial angular velocity");
 
         // Read x, y and z lists

diff --git a/source/dem/insertion.cc b/source/dem/insertion.cc
@@ -82,9 +82,9 @@ Insertion<dim>::assign_particle_properties(
       double vel_x        = dem_parameters.insertion_info.initial_vel[0];
       double vel_y        = dem_parameters.insertion_info.initial_vel[1];
       double vel_z        = dem_parameters.insertion_info.initial_vel[2];
-      double omega_x      = dem_parameters.insertion_info.initial_ang_vel[0];
-      double omega_y      = dem_parameters.insertion_info.initial_ang_vel[1];
-      double omega_z      = dem_parameters.insertion_info.initial_ang_vel[2];
+      double omega_x      = dem_parameters.insertion_info.initial_omega[0];
+      double omega_y      = dem_parameters.insertion_info.initial_omega[1];
+      double omega_z      = dem_parameters.insertion_info.initial_omega[2];
       double fem_force_x  = 0.;
       double fem_force_y  = 0.;
       double fem_force_z  = 0.;