Skip to content
Snippets Groups Projects

CONFIG: add settings for Cray compiler and cray mpich

Merged CONFIG: add settings for Cray compiler and cray mpich
feature-config-cray into develop
Merged Mark OLESEN requested to merge feature-config-cray into develop
Viewing commit f1780abe
Show latest version
4 files
+ 19
− 26
Preferences
Compare changes
Files
4
  • f1780abe 
    The 4.x tracking enforces reduced dimensionality on the parcels by
moving them to the centre of the mesh at the start of each track,
without considering the topology. This can leave the parcel outside it's
associated tetrahedron.

The barycentric algorithm isn't tolerant to incorrect topology, so
instead of changing position, it was written to track to the mesh
centre. This worked, but effectively doubled the number of tracking
calls. This additional cost has now been removed by absorbing the
constraint displacement into the existing motion track, so that the same
number of tracking steps are performed as before.

Partially resolves bug report https://bugs.openfoam.org/view.php?id=2688
src/lagrangian/DSMC/parcels/Templates/DSMCParcel/DSMCParcel.C
+ 4
− 4
@@ -53,17 +53,17 @@ bool Foam::DSMCParcel<ParcelType>::move
while (td.keepParticle && !td.switchProcessor && p.stepFraction() < 1)
{
// Apply correction to position for reduced-D cases
p.constrainToMeshCentre();
Utracking = U_;
// Apply correction to velocity to constrain tracking for
// reduced-D cases
meshTools::constrainDirection(mesh, mesh.solutionD(), Utracking);
// Deviation from the mesh centre for reduced-D cases
const vector d = p.deviationFromMeshCentre();
const scalar f = 1 - p.stepFraction();
p.trackToAndHitFace(f*trackTime*Utracking, f, cloud, td);
p.trackToAndHitFace(f*trackTime*Utracking - d, f, cloud, td);
}
return td.keepParticle;