Skip to content

GitLab

  • Menu
Projects Groups Snippets
    • Loading...
  • Help
    • Help
    • Support
    • Community forum
    • Submit feedback
    • Contribute to GitLab
  • Sign in / Register
  • openfoam openfoam
  • Project information
    • Project information
    • Activity
    • Labels
    • Planning hierarchy
    • Members
  • Repository
    • Repository
    • Files
    • Commits
    • Branches
    • Tags
    • Contributors
    • Graph
    • Compare
  • Issues 390
    • Issues 390
    • List
    • Boards
    • Service Desk
    • Milestones
  • Merge requests 13
    • Merge requests 13
  • Deployments
    • Deployments
    • Releases
  • Wiki
    • Wiki
  • Activity
  • Graph
  • Create a new issue
  • Commits
  • Issue Boards
Collapse sidebar
  • Development
  • openfoamopenfoam
  • Issues
  • #2085

Closed
Open
Created May 04, 2021 by Nikola Majksner@majksner

snappyHexMesh bumps (bulges) when using different refinementSurfaces levels

Summary

When using different refinementSurfaces levels, for example 5 6, mesh gets bumps (bulges) and while possible to get rid of them by changing minTetQuality to -1e+30 that change makes CFD very unstable and in most cases results in a crash.

  • Bumps mainly appear when reducing skewness from 20 4 to 4 2 and lower. Leaving skewness high again comes at the cost of stability during the run.
  • Bumps often occur at the transition from one refinement level to another
  • The bumps go away when using the same refinement level everywhere on the surface but that eliminates the possibility for local refinement

Steps to reproduce

Run attached case and open Paraview, and you'll see bumps the top of the wind shield.

Example case

See attached

What is the current bug behaviour?

Bumps appear at the transition from one refinement level to another

What is the expected correct behaviour?

Bumps should not be appearing on the surface when using different refinement levels

Relevant logs and/or images

This image shows bumps for maxBoundarySkewness 20 maxInternalSkewness 4 which is default value

aquilo_bumps

Using maxBoundarySkewness 2 maxInternalSkewness 2 bumps are even more pronounced

aquilo_bumps_2_2

Environment information

  • OpenFOAM version : 2012
  • Operating system : Ubuntu 20.04
  • Hardware info : AMD EPYC 7401P 24-Core Processor
  • Compiler : GCC

Possible fixes

  • Increase skewness (to the default values of 20 4) (but reduces solver stability)
  • set minTetquality to -1e+30 (but reduces solver stability significantly)

bumps.zip

Assignee
Assign to
Time tracking