|
[Sponsors] |
linear momentum calculation in KinematicCloudI.H |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
April 6, 2012, 13:22 |
linear momentum calculation in KinematicCloudI.H
|
#1 |
New Member
Chi-Wei Tsang
Join Date: Mar 2012
Posts: 2
Rep Power: 0 |
Dear all,
I am confused with the way OpenFOAM calculates linear momentum of system: template<class CloudType> inline Foam::vector Foam::KinematicCloud<CloudType>::linearMomentumOfS ystem() const { vector linearMomentum(vector::zero); forAllConstIter(typename KinematicCloud<CloudType>, *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<class CloudType> inline Foam::scalar Foam::KinematicCloud<CloudType>::massInSystem() const { scalar sysMass = 0.0; forAllConstIter(typename KinematicCloud<CloudType>, *this, iter) { const parcelType& p = iter(); sysMass += p.mass()*p.nParticle(); } return sysMass; } What do you think? Thank you! Regards |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Momentum and displacement thickness: du/dy calculation | Josh | CFX | 13 | October 27, 2013 15:38 |
Momentum source calculation | kavitagangisetti | FLUENT | 0 | June 19, 2010 15:29 |
solution diverges when linear upwind interpolation scheme is used | subash | OpenFOAM | 0 | May 29, 2010 02:23 |
Derivation of Momentum Equation in Integral Form | Demonwolf | Main CFD Forum | 2 | October 29, 2009 20:53 |
Warning 097- | AB | Siemens | 6 | November 15, 2004 05:41 |