linear momentum calculation in KinematicCloudI.H

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

 April 6, 2012, 12:22 linear momentum calculation in KinematicCloudI.H #1 New Member   Chi-Wei Tsang Join Date: Mar 2012 Posts: 2 Rep Power: 0 Sponsored Links Dear all, I am confused with the way OpenFOAM calculates linear momentum of system: template inline Foam::vector Foam::KinematicCloud::linearMomentumOfS ystem() const { vector linearMomentum(vector::zero); forAllConstIter(typename KinematicCloud, *this, iter) { const parcelType& p = iter(); linearMomentum += p.mass()*p.U(); } return linearMomentum; } I think that p.mass() returns the mass of a "particle" in a parcel. So the mass of a parcel should be p.mass() * p.nParticle(). Then the linearMomentum should be p.mass() * p.nParticle() * p.U() instead of p.mass() * p.U(). For your reference, below is the way how OpenFOAM calculates total mass of a system: template inline Foam::scalar Foam::KinematicCloud::massInSystem() const { scalar sysMass = 0.0; forAllConstIter(typename KinematicCloud, *this, iter) { const parcelType& p = iter(); sysMass += p.mass()*p.nParticle(); } return sysMass; } What do you think? Thank you! Regards

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Josh CFX 13 October 27, 2013 15:38 kavitagangisetti FLUENT 0 June 19, 2010 14:29 subash OpenFOAM 0 May 29, 2010 01:23 Demonwolf Main CFD Forum 2 October 29, 2009 20:53 AB Siemens 6 November 15, 2004 05:41