[CoolHand]

I'm posting this to Solving because I think it's a problem with Solver/Mesh interaction. Please let me know if its not the right place for it.

I have a 2D simulation of flow over a half cylinder. It works perfectly until I use refineMesh on it. After using refineMesh, checkMesh reports:

 
***Number of edges not aligned with or perpendicular to non-empty directions: 1374
  <<Writing 2061 points on non-aligned edges to set nonAlignedEdges

and if I try to run the simulation it gives the error:

Starting time loop

Time = 1

#0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2   in "/lib/x86_64-linux-gnu/libc.so.6"
#3   at symmTensorField.C:0
#4  Foam::inv(Foam::UList<Foam::SymmTensor<double> > const&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#5  Foam::leastSquaresVectors::makeLeastSquaresVectors() const in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#6  Foam::leastSquaresVectors::pVectors() const in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#7  Foam::fv::leastSquaresGrad<Foam::Vector<double> >::calcGrad(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::word const&) const in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#8  Foam::fv::cellLimitedGrad<Foam::Vector<double> >::calcGrad(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::word const&) const in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#9  Foam::fv::gradScheme<Foam::Vector<double> >::grad(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::word const&) const in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libincompressibleTurbulenceModel.so"
#10  Foam::tmp<Foam::GeometricField<Foam::outerProduct<Foam::Vector<double>, Foam::Vector<double> >::type, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::grad<Foam::Vector<double> >(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libincompressibleTurbulenceModel.so"
#11  Foam::incompressible::RASModels::kEpsilon::divDevReff(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>&) const in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
#12  
 in "/opt/openfoam220/platforms/linux64GccDPOpt/bin/simpleFoam"
#13  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#14  
 in "/opt/openfoam220/platforms/linux64GccDPOpt/bin/simpleFoam"
Floating point exception (core dumped)

The refinement goes across the bottom of the domain from the inlet to the outlet. If I reduce the refinement to be only 1 radius upstream and downstream of the cylinder I don't get the errors.

Here is my case file:

http://www.filedropper.com/matlabmesh

The mesh was created in the GMSH format by a MatLab code that I wrote. The refined mesh is in the 1 folder.

Any ideas?