This is a branch to test the different methods to apply no slip boundary conditions on IBs, - The important part of this commit is applying the one-sided penalty terms only inside the IB, not outside; this is important to make sure that the boundary layer around the bodies are not modified; This finally makes sure that drag calculation is not affected by the source terms introduced at the fluid-solid interface. - Another thing that is tested in this commit is adding dV*(1.0 - H_s)*C_vol*(u-u_s)*rho term to the force terms; this is to make sure that the momentum loss inside the IBs are taken into account as well. However, this requiers that u=us is satisfied inside the IB up to a fairly tight tolerance, otherwise large alpha and beta constants will results in a non-zero contribution in this exrpression.
This is a branch to test the different methods to apply no slip boundary conditions on IBs,
- The important part of this commit is applying the one-sided penalty terms only inside the IB, not outside; this is important to make sure that the boundary layer around the bodies are not modified; This finally makes sure that drag calculation is not affected by the source terms introduced at the fluid-solid interface.
- Another thing that is tested in this commit is adding dV*(1.0 - H_s)*C_vol*(u-u_s)*rho term to the force terms; this is to make sure that the momentum loss inside the IBs are taken into account as well. However, this requiers that u=us is satisfied inside the IB up to a fairly tight tolerance, otherwise large alpha and beta constants will results in a non-zero contribution in this exrpression.

erdc/proteus,erdc/proteus,erdc/proteus,erdc/proteus