Friction in non-Newtonian DualSPHysics

salehmeiabadi

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.

Alex

@gfourtakas should know the answer!!!

salehmeiabadi

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

Alex

@gfourtakas can help here

LLL

Hi, salehmeiabadi, I have the same question. Have you found an answer?

salehmeiabadi

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,