zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

Small changes to output during best guess determination.

Browse Source 
Making intersection test tolerance SMALL.
This commit is contained in:
graham
2009-10-09 11:01:26 +01:00
parent 3ef3638c44
commit 5488acac27
2 changed files with 2 additions and 29 deletions
Download Patch File Download Diff File Expand all files Collapse all files

2
src/lagrangian/basic/Cloud/Cloud.C
View File

@ -42,7 +42,7 @@ template<class ParticleType>
const Foam::scalar Foam::Cloud<ParticleType>::trackingRescueTolerance = 1e-6;
template<class ParticleType>
const Foam::scalar Foam::Cloud<ParticleType>::intersectionTolerance = 1e-3;
const Foam::scalar Foam::Cloud<ParticleType>::intersectionTolerance = SMALL;
template<class ParticleType>
const Foam::scalar Foam::Cloud<ParticleType>::planarCosAngle = (1 - 1e-6);

29
src/lagrangian/basic/Particle/Particle.C
View File

@ -281,6 +281,7 @@ void Foam::Particle<ParticleType>::trackToFaceConcave
Pout<< "Need a best guess " << nl
<< position_ << nl
<< stepFraction_ << nl
<< position_ + delta << nl
<< position_ - delta
<< endl;
@ -315,34 +316,6 @@ void Foam::Particle<ParticleType>::trackToFaceConcave
tmpLambdas.append(inter.distance());
tmpLambdaFaceIs.append(facei);
}
inter = mesh.faces()[facei].intersection
(
position_,
deltaPosition,
mesh.faceCentres()[facei],
mesh.points(),
intersection::HALF_RAY,
Cloud<ParticleType>::intersectionTolerance
);
Pout<< "Test repeat HALF "
<< facei << " " << inter.distance() << endl;
inter = mesh.faces()[facei].intersection
(
position_,
deltaPosition,
mesh.faceCentres()[facei],
mesh.points(),
intersection::FULL_RAY,
Cloud<ParticleType>::intersectionTolerance
);
Pout<< "Test repeat FULL "
<< facei << " " << inter.distance() << nl
<< (inter.distance()*deltaPosition) + position_
<< endl;
}
Pout<< tmpLambdaFaceIs << " " << tmpLambdas << endl;