Friction in non-Newtonian DualSPHysics
Hello,
How can we apply friction between boundaries and fluid in non-Newtonian DualSPHysics?
Does ViscoBoundFactor works for NN_DSPH?
If it doesn't, Is there any solution to apply friction in NN_DSPH?
Thank you so much for sharing your knowledge.
Comments
@gfourtakas should know the answer!!!
Hello to the DualSPHysics team
I can see that ViscoBoundFactor is available in Interaction Fluid-Bound function.
However, I am not sure how ViscoBoundFactor works in non_Newtonian DualSPHysics.
//-ConsEq;
// Build strain rate tensor and compute eta_visco //what will happen here for boundary particles?
InteractionForcesFluid_NN_SPH_Visco_Stress_tensor<ftmode, tvisco>(t.npf, t.npb, Visco
, t.visco_eta, t.spstau, t.d_tensor, t.auxnn, t.pos, t.code, t.idp);
//-Compute viscous terms and add them to the acceleration
InteractionForcesFluid_NN_SPH_ConsEq<tker, ftmode, tvisco, tdensity, shift>(t.npf, t.npb, false, Visco
, t.divdata, t.dcell, t.visco_eta, t.spstau, t.auxnn, t.pos, t.velrhop, t.code, t.idp, t.ace, t.spsgradvel, t.d_tensor, t.pvf, t.apvf, t.sh_rate, t.viscosity, t.ar); //zxcv
//-Interaction Fluid-Bound.
InteractionForcesFluid_NN_SPH_ConsEq<tker, ftmode, tvisco, tdensity, shift>(t.npf, t.npb, true, Visco*ViscoBoundFactor
, t.divdata, t.dcell, t.visco_eta, t.spstau, t.auxnn, t.pos, t.velrhop, t.code, t.idp, t.ace, t.spsgradvel, t.d_tensor, t.pvf, t.apvf, t.sh_rate, t.viscosity, t.ar); //zxcv
}
Please advise me in this regard.
Thanks
@gfourtakas can help here
Hi, salehmeiabadi, I have the same question. Have you found an answer?
My understanding is that friction between fluid-boundary particles in non-Newtonian DualSPHysics is as much as friction between fluid-fluid particles.
I would be appreciated it if someone from non-Newtonian DualSPHysics confirm or reject that.
Regards,