Commit a61b0967 authored by Thomas Fetzer's avatar Thomas Fetzer
Browse files

[freeflow][rans] Set velocity gradient to zero, if it is a single element in on direction

parent a3cabb38
...@@ -265,6 +265,9 @@ public: ...@@ -265,6 +265,9 @@ public:
- velocity_[neighborIdx_[elementIdx][dimIdx][0]][velIdx]) - velocity_[neighborIdx_[elementIdx][dimIdx][0]][velIdx])
/ (cellCenter_[neighborIdx_[elementIdx][dimIdx][1]][dimIdx] / (cellCenter_[neighborIdx_[elementIdx][dimIdx][1]][dimIdx]
- cellCenter_[neighborIdx_[elementIdx][dimIdx][0]][dimIdx]); - cellCenter_[neighborIdx_[elementIdx][dimIdx][0]][dimIdx]);
if (abs(cellCenter_[neighborIdx_[elementIdx][dimIdx][1]][dimIdx]
- cellCenter_[neighborIdx_[elementIdx][dimIdx][0]][dimIdx]) < 1e-8)
velocityGradients_[elementIdx][velIdx][dimIdx] = 0.0;
} }
if (abs(velocity_[elementIdx][dimIdx]) > abs(velocityMaximum_[wallElementIdx][dimIdx])) if (abs(velocity_[elementIdx][dimIdx]) > abs(velocityMaximum_[wallElementIdx][dimIdx]))
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment