Merge pull request #655 from ZhentongWang/Eulerian_codes_update

modify the velocity in wall boundary condition in a general form.
Xiangyu-Hu · Sep 12, 2024 · 1acaaa8 · 1acaaa8
2 parents 129179c + c75ceb2
commit 1acaaa8
Showing 1 changed file with 4 additions and 2 deletions.
diff --git a/...d/particle_dynamics/fluid_dynamics/eulerian_fluid_dynamics/eulerian_fluid_integration.hpp b/...d/particle_dynamics/fluid_dynamics/eulerian_fluid_dynamics/eulerian_fluid_integration.hpp
@@ -65,14 +65,15 @@ void EulerianIntegration1stHalf<Contact<Wall>, RiemannSolverType>::interaction(s
     {
         StdLargeVec<Vecd> &n_k = *(wall_n_[k]);
         StdLargeVec<Real> &Vol_k = *(wall_Vol_[k]);
+        StdLargeVec<Vecd> &vel_ave_k = *(wall_vel_ave_[k]); 
         Neighborhood &wall_neighborhood = (*contact_configuration_[k])[index_i];
         for (size_t n = 0; n != wall_neighborhood.current_size_; ++n)
         {
             size_t index_j = wall_neighborhood.j_[n];
             Vecd &e_ij = wall_neighborhood.e_ij_[n];
             Real dW_ijV_j = wall_neighborhood.dW_ij_[n] * Vol_k[index_j];
 
-            Vecd vel_in_wall = -state_i.vel_;
+            Vecd vel_in_wall = 2.0 * vel_ave_k[index_j] - state_i.vel_;
             Real p_in_wall = state_i.p_;
             Real rho_in_wall = state_i.rho_;
             FluidStateIn state_j(rho_in_wall, vel_in_wall, p_in_wall);
@@ -132,14 +133,15 @@ void EulerianIntegration2ndHalf<Contact<Wall>, RiemannSolverType>::interaction(s
     {
         StdLargeVec<Vecd> &n_k = *(this->wall_n_[k]);
         StdLargeVec<Real> &Vol_k = *(this->wall_Vol_[k]);
+        StdLargeVec<Vecd> &vel_ave_k = *(wall_vel_ave_[k]);
         Neighborhood &wall_neighborhood = (*contact_configuration_[k])[index_i];
         for (size_t n = 0; n != wall_neighborhood.current_size_; ++n)
         {
             size_t index_j = wall_neighborhood.j_[n];
             Vecd &e_ij = wall_neighborhood.e_ij_[n];
             Real dW_ijV_j = wall_neighborhood.dW_ij_[n] * Vol_k[index_j];
 
-            Vecd vel_in_wall = -state_i.vel_;
+            Vecd vel_in_wall = 2.0 * vel_ave_k[index_j] - state_i.vel_;
             Real p_in_wall = state_i.p_;
             Real rho_in_wall = state_i.rho_;