code-saturne / code-saturne/code_saturne

Incorrect wall distance using axisymmetric mesh

Open
#41 0 comments 0 reactions 0 assignees View on GitHub
Dominant language
C++
Stars
307
Forks
102
PR merge metrics
No merged PRs in 30d

Description

On "axysymmetric" cases (where the mesh uses a radially extruded slice, with symmetry on the sides), the wall distance approximated by a diffusion equation (in ``distpr``) can be imprecise, with a difference up to 40% having been observed.

This issue also occurs using the CDO operators, so seems due to the the choice of equation or boundary conditions rather than linear operator.

In any case, this can induce an error on turbulence variables.

More works needs to be done on this, but replacing the approximate wall distance with a
geometrically computed wall distance for the boundary cells might at least mitigate the effects.

Attached is a small case allowing to test this.
[test_pipe_axi.tar.gz](https://github.com/code-saturne/code_saturne/files/3135105/test_pipe_axi.tar.gz)

Contributor guide

Open the contributing guide

Research direction

Start with the attached test_pipe_axi.tar.gz case and the distpr wall-distance calculation; compare its axisymmetric result with the CDO-operator behavior. Investigate whether the equation or boundary conditions cause the reported 40% difference, and assess the proposed geometrically computed distance for boundary cells. Done means the test case no longer shows the problematic wall distance or its turbulence-variable error.

Written by the indexing model from the issue text.

Assessment

Tech stack
cpp
Domain
hpc
Issue type
Bug
Difficulty
5/5
Estimated time
Over a week
Activity status
Stale
Clarity
Needs clarification
Newbie friendliness
20/100

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.