[Jacobhh]

Hi,
I am simulating turbulent flow within a cyclic pipe.
The flow perturbed with a modified version of Eugene De Villers perturbU, and when it the turbulent profile begins developing/ is developed certain cells within the mesh exhibits weird behavior. They have a significantly lower pressure than the surrounding cells, and nut and U is hence effected.
I have tried with three different meshes, and all of them behaves the same.
pisoFoam is used as the solver, cyclicAMI is used as boundary conditions for inlet/outlet and WALE is the turbulence model running on OF v2106

Have any of you had this problem before?

Please see attached pictures of the mean pressure and nut. Did not know the relevant code to post.

Best,
Jacob

[Jacobhh]

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Time = 0

Mesh stats
    points:           36818
    faces:            237648
    internal faces:   229522
    cells:            103937
    faces per cell:   4.49474
    boundary patches: 3
    point zones:      0
    face zones:       0
    cell zones:       0

Overall number of cells of each type:
    hexahedra:     0
    prisms:        51422
    wedges:        0
    pyramids:      0
    tet wedges:    0
    tetrahedra:    52515
    polyhedra:     0

Checking topology...
    Boundary definition OK.
    Cell to face addressing OK.
    Point usage OK.
    Upper triangular ordering OK.
    Face vertices OK.
    Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces...
    Patch               Faces    Points   Surface topology                  
    Outlet              390      320      ok (non-closed singly connected)  
    wall                7346     3704     ok (non-closed singly connected)  
    Inlet               390      320      ok (non-closed singly connected)  

Checking faceZone topology for multiply connected surfaces...
    No faceZones found.

Checking basic cellZone addressing...
    No cellZones found.

Checking geometry...
    Overall domain bounding box (-0.05 -0.05 -6.09182e-18) (0.05 0.05 1)
    Mesh has 3 geometric (non-empty/wedge) directions (1 1 1)
    Mesh has 3 solution (non-empty) directions (1 1 1)
    Boundary openness (4.67173e-16 -5.10727e-17 -1.49361e-17) OK.
    Max cell openness = 5.44332e-16 OK.
    Max aspect ratio = 7.55135 OK.
    Minimum face area = 6.0242e-06. Maximum face area = 0.000124292.  Face area magnitudes OK.
    Min volume = 1.88013e-08. Max volume = 4.37891e-07.  Total volume = 0.00781563.  Cell volumes OK.
    Mesh non-orthogonality Max: 49.1917 average: 11.8223
    Non-orthogonality check OK.
    Face pyramids OK.
    Max skewness = 0.462065 OK.
    Coupled point location match (average 0) OK.

Mesh OK.

End