From ceca5c62d6f2635a59f4f4f361cdaea53add108c Mon Sep 17 00:00:00 2001
From: graham <g.macpherson@opencfd.co.uk>
Date: Thu, 24 Sep 2009 12:47:57 +0100
Subject: [PATCH] Revert "Adding constraint of tracking motion in reduced
dimension cases."
This reverts commit 9d3f0eaefb021db3ad0e765589d7e5f9ebb3c05f.
Need to do rigorously including implications for gathering data.
Bird's method for a 2D code involves all vector quantities being
actually 2D, and the density (for example) being particles per unit
area (flowfield) or length (surface). Needs some thought to
generalise.
---
src/lagrangian/dsmc/Make/options | 6 ++----
.../dsmc/parcels/Templates/DsmcParcel/DsmcParcel.C | 14 --------------
2 files changed, 2 insertions(+), 18 deletions(-)
diff --git a/src/lagrangian/dsmc/Make/options b/src/lagrangian/dsmc/Make/options
index 55865dfabcd..15874b7b55a 100644
--- a/src/lagrangian/dsmc/Make/options
+++ b/src/lagrangian/dsmc/Make/options
@@ -1,9 +1,7 @@
EXE_INC = \
-I$(LIB_SRC)/finiteVolume/lnInclude \
- -I$(LIB_SRC)/lagrangian/basic/lnInclude \
- -I$(LIB_SRC)/meshTools/lnInclude
+ -I$(LIB_SRC)/lagrangian/basic/lnInclude
LIB_LIBS = \
-llagrangian \
- -lfiniteVolume \
- -lmeshTools
+ -lfiniteVolume
diff --git a/src/lagrangian/dsmc/parcels/Templates/DsmcParcel/DsmcParcel.C b/src/lagrangian/dsmc/parcels/Templates/DsmcParcel/DsmcParcel.C
index 3742f3f8ea8..ffd7c67d428 100644
--- a/src/lagrangian/dsmc/parcels/Templates/DsmcParcel/DsmcParcel.C
+++ b/src/lagrangian/dsmc/parcels/Templates/DsmcParcel/DsmcParcel.C
@@ -25,7 +25,6 @@ License
\*---------------------------------------------------------------------------*/
#include "DsmcParcel.H"
-#include "meshTools.H"
// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
@@ -48,18 +47,8 @@ bool Foam::DsmcParcel<ParcelType>::move
scalar tEnd = (1.0 - p.stepFraction())*deltaT;
const scalar dtMax = tEnd;
- // Save the velocity to re-apply it after tracking
- vector U_save = U_;
-
- // Apply correction to velocity to constrain tracking for
- // reduced-D cases
- meshTools::constrainDirection(mesh, mesh.solutionD(), U_);
-
while (td.keepParticle && !td.switchProcessor && tEnd > ROOTVSMALL)
{
- // Apply correction to position for reduced-D cases
- meshTools::constrainToMeshCentre(mesh, p.position());
-
// Set the Lagrangian time-step
scalar dt = min(dtMax, tEnd);
@@ -78,9 +67,6 @@ bool Foam::DsmcParcel<ParcelType>::move
}
}
- // Restore the correct value of velocity
- U_ = U_save;
-
return td.keepParticle;
}
--
GitLab