CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Bugs

possible bug in KinematicParcel

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Display Modes
Old   April 19, 2011, 04:20
Default possible bug in KinematicParcel
New Member
Guanghao Wu
Join Date: Mar 2009
Location: Tokyo, Japan
Posts: 15
Rep Power: 8
guanghaowu is on a distinguished road

Suppose a particle is in a cell next to rebound type wall, while the fluid velocity of the cell center is toward to the wall.
If the particle hits the wall,
U_ = rebound of U_,
U_ = Uc_;
where Uc_ is toward to the wall, so in the next WHILE iteration, dt becomes very small 1E-19 order, but still
and the WHILE iteration does not stop...

In KinematicParcel.C

template<class ParcelType>
template<class TrackData>
bool Foam::KinematicParcel<ParcelType>::move(TrackData& td)
while (td.keepParticle && !td.switchProcessor && tEnd > ROOTVSMALL)
if (
dt *= p.trackToFace(p.position() + dt*U_, td);

tEnd -= dt;
p.stepFraction() = 1.0 - tEnd/deltaT;

// Avoid problems with extremely small timesteps
if (dt > ROOTVSMALL) //=> if (dt > ROOTVSMALL && !p.onBoundary())
// Update cell based properties
p.setCellValues(td, dt, cellI);

if (
p.cellValueSourceCorrection(td, dt, cellI);

p.calc(td, dt, cellI);

if (p.onBoundary() && td.keepParticle)
if (isA<processorPolyPatch>(pbMesh[p.patch(p.face())]))
td.switchProcessor = true;

Best regards,
guanghaowu is offline   Reply With Quote


Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Serious bug in LES interface fs82 OpenFOAM Bugs 21 November 16, 2009 09:15
bug in Rcomponents utility cedric_duprat OpenFOAM Bugs 1 May 7, 2009 02:56
Maybe a bug in CFX12.0 BETA? frank CFX 4 November 26, 2008 19:08
Bug reports Mattijs Janssens (Mattijs) OpenFOAM 0 January 10, 2005 11:05
Forum y2k Bug Jonas Larsson Main CFD Forum 1 January 5, 2000 11:22

All times are GMT -4. The time now is 04:50.