CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > Siemens

Dropro.f subroutine for multi-component droplets

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

LinkBack Thread Tools Search this Thread Display Modes
Old   February 20, 2012, 08:26
Default Dropro.f subroutine for multi-component droplets
Join Date: Oct 2011
Location: Thessaloniki, Greece
Posts: 75
Rep Power: 13
crevoise is on a distinguished road
Helo all.

I am trying to compute the evaporation process of AdBlue (UWS) droplets along a pipe.
I am defining the droplets as a mixture of two components (urea and water), with one evaporating to water(g) and the other one to urea(g).

I wish to know few detailed in order to write my subroutine in a proper way. For the moment, the computation is running, but the results are not satisfactory. The urea should evaporate only once the water totally evaporated; however, I am getting some vapor urea at the same as water vapor.

About the bulk properties, the subroutine should be ok, as I tested these datae from a single component droplet.
The concern is about the datae definition of the components, especially the vapor pressure, as this one is of main importance in my case:
1/ While defining the data, I am using this way:
if ico 1 (water), then pvH2O= .... , and finish with pv=pvH20
if ico 2 (urea), then pvurea= .... , and finish with pv=pvurea

Is it done in the right way? It seems that Star will then manage my two components in independent way, so will compute their evaporation in parallel? Can't it be return to Star a pv1 and pv2?

2/ About the same vapor pressure, I was wondering if I need to write in dropro.f the Raoult's law - thus to write pv=0.6664 * pvH20 ? (with 0.6664 the water mass fraction)

3/ Together with the fact evaporation of urea is happening too early, the droplets remains at 320K (which is the evaporation temperature of the problem one using only pure water)
Trying solving this point, I have wrote in the subroutine to start computing the vapor pressure of urea, only once the water mass fraction was 0 (if DRMFR(1) .le. 0), then...)
However, the droplet temperature is then increasing without stopping to the urea 'boiling' point

So the question is, do I need to define for every component a boiling temperature, or the bulk boiling temperature is enough?
If so, how to do it?

4/ I was wondering if using the mixture model was appropriated to this problem. Does anybody think I need to use a emulsified mixture?

5/Does anybody knows where to find some advice or supplementary note about the definition of user subroutine in StarCD, to get deeper knowledge on them?

Thank you for your advices and any help.
crevoise is offline   Reply With Quote


droplet, dropro.f, emulsified, mixture, multi component

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Momentum source coefficient, cylindrical coordinates, circumferential component? polakse CFX 15 January 18, 2016 02:40
Fortran Subroutine Error sergiorbrett CFX 12 September 3, 2012 07:40
Errorin subroutine appeared when applying cavitation model pitisrisuk CFX 1 July 2, 2012 03:36
Mass & Heat Transfer : Problems with the droplets sheintz STAR-CD 0 June 7, 2010 08:48
UDF component in Cyliderical Velocity? Hongchan Kim FLUENT 3 April 24, 2001 12:20

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