CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Particle Reynolds number calculation in Lagrangian tracking?

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 19, 2012, 02:21
Default Particle Reynolds number calculation in Lagrangian tracking?
  #1
Senior Member
 
Jie
Join Date: Jan 2010
Location: Australia
Posts: 119
Rep Power: 7
jiejie is on a distinguished road
Dear foamers

I am trying to implement the Srudhar and Katz bubble motion equation in the OpenFOAM. Hence, I looked at the exisiting Sphere Drag force model. I found there might be a possible mistake in calculating the particle Reynolds number. Particle Reynolds number is defined as Re_p = |U_rel|*d_p/nu_c, where U_rel is the relatvie velocity between the particle and carrying phase, d_p is the particle diameter and nu_c is the kinematic viscosity of the carrying phase. However, I found OpenFOMA uses the particle velocity (U_) instead of the relative velocity to calculate the particle Reynolds number (see line 110 in http://foam.sourceforge.net/docs/cpp/a05113_source.html).

I am not sure whether this makes sense to everyone. Did anyone have this problem before?

Thanks a lot.

Jie
jiejie is offline   Reply With Quote

Old   June 19, 2012, 02:31
Default
  #2
Super Moderator
 
niklas's Avatar
 
Niklas Nordin
Join Date: Mar 2009
Location: Stockholm, Sweden
Posts: 693
Rep Power: 19
niklas will become famous soon enough
You are refering to this line:
Code:
    // Reynolds number
    const scalar Re = this->Re(U_, d_, rhoc_, muc_);
but if you check KinematicParcelI.H you will find that
Code:

template<class ParcelType>
inline Foam::scalar Foam::KinematicParcel<ParcelType>::Re
(
    const vector& U,
    const scalar d,
    const scalar rhoc,
    const scalar muc
) const
{
    return rhoc*mag(U - Uc_)*d/(muc + ROOTVSMALL);
}
which is what you want
niklas is offline   Reply With Quote

Old   June 19, 2012, 02:50
Default
  #3
Senior Member
 
Jie
Join Date: Jan 2010
Location: Australia
Posts: 119
Rep Power: 7
jiejie is on a distinguished road
Quote:
Originally Posted by niklas View Post
but if you check KinematicParcelI.H you will find that
Code:
    return rhoc*mag(U - Uc_)*d/(muc + ROOTVSMALL);
which is what you want
Hi niklas

Thanks for the quick reply.This means the U_rel is defined as "U_p - U_c". The drag force is defined as

F_d = 0.5*Cd*rho_c*A_p*|U_c-U_p|*(U_c-U_p).

If OpenFOAM uses the U_rel = U_p - U_c.

Then, should F_d = -0.5*Cd*rho_c*A_p*|U_p-U_c|*(U_p-U_c) instead?

Jie
jiejie is offline   Reply With Quote

Old   June 19, 2012, 03:05
Default
  #4
Senior Member
 
Jie
Join Date: Jan 2010
Location: Australia
Posts: 119
Rep Power: 7
jiejie is on a distinguished road
I think I found the answer: F = Sp(U - Up) + Su

When they evaluate the force with the value.Sp, it is multiplied by (U_c-U_p), hence, there should not be a negative sign in the front.

Jie
jiejie is offline   Reply With Quote

Old   July 5, 2012, 10:13
Default
  #5
MMK
New Member
 
Join Date: Jan 2011
Posts: 9
Rep Power: 6
MMK is on a distinguished road
Quote:
Originally Posted by niklas View Post
You are refering to this line:
Code:
    // Reynolds number
    const scalar Re = this->Re(U_, d_, rhoc_, muc_);
but if you check KinematicParcelI.H you will find that
Code:

template<class ParcelType>
inline Foam::scalar Foam::KinematicParcel<ParcelType>::Re
(
    const vector& U,
    const scalar d,
    const scalar rhoc,
    const scalar muc
) const
{
    return rhoc*mag(U - Uc_)*d/(muc + ROOTVSMALL);
}
which is what you want
Hi nikolas,

I have a similar problem. In dieselFoam all non dimensional numbers are calculated in Parcel.H
but none of them gets written in the output data fields. Please can you elaborate how to get these non dimentional numbers at the output.

e.g. in parcel.H

scalar We
(
const vector& U,
const scalar rho,
const scalar sigma
) const;


And more over there is parcelFunctions.C which calculates theses non-dimensional numbers, is there any way to benefit from these.
MMK is offline   Reply With Quote

Old   July 6, 2012, 04:47
Default
  #6
MMK
New Member
 
Join Date: Jan 2011
Posts: 9
Rep Power: 6
MMK is on a distinguished road
Quote:
Originally Posted by MMK View Post
Hi nikolas,

I have a similar problem. In dieselFoam all non dimensional numbers are calculated in Parcel.H
but none of them gets written in the output data fields. Please can you elaborate how to get these non dimentional numbers at the output.

e.g. in parcel.H

scalar We
(
const vector& U,
const scalar rho,
const scalar sigma
) const;


And more over there is parcelFunctions.C which calculates theses non-dimensional numbers, is there any way to benefit from these.
Please Can any one answer,
MMK is offline   Reply With Quote

Reply

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
lagrangian solid particle tracking - OFv2.1 PelusDadidus OpenFOAM Programming & Development 0 April 3, 2012 07:08
Particle tracking in post (How to get number in a given plane?) DPD CFX 1 February 18, 2011 02:28
massless particle tracking problem Renold FLUENT 0 January 26, 2011 15:23
Particle Tracking Batis CFX 2 October 6, 2010 15:20
Warning 097- AB CD-adapco 6 November 15, 2004 05:41


All times are GMT -4. The time now is 19:48.