Skip to content

Fix wallfunction variable computation on mpi rank edges #1791

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
bigfooted merged 3 commits into from
Oct 17, 2022
Merged

Fix wallfunction variable computation on mpi rank edges #1791

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
75cb6ea
fix yplus on mpi rank edge
bigfooted Oct 14, 2022
cce7ff4
fix for mg cycles
bigfooted Oct 14, 2022
19d6e5e
Apply suggestions from code review
pcarruscag Oct 16, 2022
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
7 changes: 2 additions & 5 deletions SU2_CFD/src/solvers/CIncNSSolver.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -647,11 +647,8 @@ void CIncNSSolver::SetTau_Wall_WF(CGeometry *geometry, CSolver **solver_containe

const auto iPoint = geometry->vertex[iMarker][iVertex]->GetNode();
const auto Point_Normal = geometry->vertex[iMarker][iVertex]->GetNormal_Neighbor();

/*--- Check if the node belongs to the domain (i.e, not a halo node)
* and the neighbor is not part of the physical boundary ---*/

if (!geometry->nodes->GetDomain(iPoint)) continue;
/*--- On the finest mesh compute also on halo nodes to avoid communication of tau wall. ---*/
if ((!geometry->nodes->GetDomain(iPoint)) && !(MGLevel==MESH_0)) continue;

/*--- Get coordinates of the current vertex and nearest normal point ---*/

Expand Down
1 change: 1 addition & 0 deletions SU2_CFD/src/solvers/CTurbSSTSolver.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -461,6 +461,7 @@ void CTurbSSTSolver::SetTurbVars_WF(CGeometry *geometry, CSolver **solver_contai

const auto iPoint = geometry->vertex[val_marker][iVertex]->GetNode();
const auto iPoint_Neighbor = geometry->vertex[val_marker][iVertex]->GetNormal_Neighbor();
if (!geometry->nodes->GetDomain(iPoint_Neighbor)) continue;

su2double Y_Plus = solver_container[FLOW_SOL]->GetYPlus(val_marker, iVertex);
su2double Lam_Visc_Wall = solver_container[FLOW_SOL]->GetNodes()->GetLaminarViscosity(iPoint);
Expand Down