DES hybrid : error when using negative velocity input

Summary

When using DEShybrid scheme, if the input velocity is negative (according to co-ordinate system used), the case exit with error

Steps to reproduce

Use dummy negative x component of velocity for internalField (attached case, simply run Allrun)

Example case

annularThermalMixer.tgz

What is the current bug behaviour?

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in /lib64/libc.so.6
#3  Foam::sqrt(Foam::Field<double>&, Foam::UList<double> const&) at ??:?
#4  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> >
#5  Foam::DEShybrid<Foam::Vector<double> >::calcBlendingFactor(

What is the expected correct behavior?

Relevant logs and/or images

Environment information

OpenFOAM version : v1812 Operating system : centos Compiler : gcc

Possible fixes

magnitude or absolute value instead?