14 #ifndef MODEL_VISCWALLINTERACTION_HPP
15 #define MODEL_VISCWALLINTERACTION_HPP
17 #include "Model/ViscWallInteraction.h"
44 Vec3 vp=this->m_p->getVel();
45 Vec3 vw=this->m_wall->getVel();
47 Vec3 force=m_nu*this->m_p->getMass()*(vp-vw);
48 Vec3 pos=this->m_p->getPos();
50 this->m_p->applyForce(force,pos);
51 if(this->m_inner_flag) this->m_wall->addForce(force);
60 Vec3 vp=this->m_p->getVel();
61 Vec3 vw=this->m_wall->getVel();
63 Vec3 force=m_nu*this->m_p->getMass()*(vw-vp);