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
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?