CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Pre-Processing

Problem with Rosin Rammler

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 12, 2013, 09:19
Default Problem with Rosin Rammler
  #1
Senior Member
 
Join Date: Jan 2010
Location: Stuttgart
Posts: 118
Rep Power: 6
Chrisi1984 is on a distinguished road
Hi all,

I want to use a RosinRammler pdf for my parcels.

In the the source code I found the following:

Quote:
Foam::scalar Foam::distributionModels::RosinRammler::sample() const
{
scalar K = 1.0 - exp(-pow((maxValue_ - minValue_)/d_, n_));
scalar y = rndGen_.sample01<scalar>();
scalar x = minValue_ + d_*:ow(-log(1.0 - y*K), 1.0/n_);
return x;
}
I tried to implement exactly that in excel to calculate the values for minValue, maxValue, d and n for my desired distribution.

I callculated K and with that later the x value. For y I take a list of 100 values from 0 to 1.

Then I thought for y= 0.5 I should get the value for dv50 of my destribution. The same I did with y= 0.9 for dv90.

But the distribution in FOAM later on does not correspond to my excel values.

What is the error that I do?

Thanks for you help in advance.

Kind regards

Chrisi
Chrisi1984 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
UDF compiling problem Wouter Fluent UDF and Scheme Programming 6 June 6, 2012 04:43
engineFoam new mesh problem ayhan515 OpenFOAM 2 May 1, 2012 04:29
Gambit - meshing over airfoil wrapping (?) problem JFDC FLUENT 1 July 11, 2011 05:59
natural convection problem for a CHT problem Se-Hee CFX 2 June 10, 2007 06:29
Adiabatic and Rotating wall (Convection problem) ParodDav CFX 5 April 29, 2007 19:13


All times are GMT -4. The time now is 22:40.