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?

changed milestone to %v1906

added bug label

changed the description

mentioned in commit f2b5b2527cee44ee60191751145d76c01dfb2ae5

By Andrew Heather on 2019-03-01T11:30:44 (imported from GitLab project)

Error is related to the (too aggressive) impulse start-up condition that can lead to negative densities for compressible solvers. Added a change in commit f2b5b252 that protects the sqrt calculation - case runs smoothly afterwards, but ideally a more appropriate initial condition should be imposed.

By Andrew Heather on 2019-03-01T11:33:47 (imported from GitLab project)

closed

By Andrew Heather on 2019-03-01T11:34:11 (imported from GitLab project)