|
[Sponsors] |
Numerical artifacts within cells in OpenFOAM v2106 |
![]() |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
![]() |
![]() |
#1 |
New Member
Jacob H
Join Date: Aug 2021
Posts: 3
Rep Power: 5 ![]() |
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 |
|
![]() |
![]() |
![]() |
![]() |
#2 |
New Member
Jacob H
Join Date: Aug 2021
Posts: 3
Rep Power: 5 ![]() |
Code:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 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 |
|
![]() |
![]() |
![]() |
Tags |
cyclicami, les, openfoam, perturbu, pisofoam |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
Calculating source term total cell-zone volume parallel computation | RobV | Fluent UDF and Scheme Programming | 2 | May 30, 2022 03:46 |
[snappyHexMesh] snappyHexMesh does not detect highly skewed faces? | ptpacheco | OpenFOAM Meshing & Mesh Conversion | 1 | January 4, 2022 12:37 |
Map of the OpenFOAM Forum - Understanding where to post your questions! | wyldckat | OpenFOAM | 10 | September 2, 2021 05:29 |
OpenFOAM Training Jan-Jul 2017, Virtual, London, Houston, Berlin | CFDFoundation | OpenFOAM Announcements from Other Sources | 0 | January 4, 2017 06:15 |
Summer School on Numerical Modelling and OpenFOAM | hjasak | OpenFOAM | 5 | October 12, 2008 13:14 |