CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Programming & Development

liquidEvaporation in simpleReactingFoam

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 16, 2019, 20:37
Default liquidEvaporation in simpleReactingFoam
  #1
New Member
 
Nilay Kulkarni
Join Date: May 2018
Posts: 24
Rep Power: 6
neilk is on a distinguished road
Hello,

I am trying to model mass transfer phase change in simpleReactingFoam. When setting the properties in the reactingCloudProperties in solution for either transient - yes or no:
Code:
    active          true;
    coupled         no;
    transient      no; //yess
   calcFrequency   1;
    cellValueSourceCorrection no;
    //maxCo           1;
     maxTrackTime 		1;
      resetSourcesOnStartup	false;
For Steady, the diameter get reset to the original injected diameter but the dt used for the dMassPC calculation remains 1.
Code:
Solving 3-D cloud reactingCloud1

Cloud: reactingCloud1 injector: model1
    Added 1 new parcels

Ni = 1.80455e-11 d= 4e-06 dt = 1
Cloud: reactingCloud1
    Current number of parcels       = 1
    Current mass in system          = 4.646091741e-14
    Linear momentum                 = (8.545341647e-20 3.137059054e-22 1.413548384e-22)
   |Linear momentum|                = 8.54541092e-20
    Linear kinetic energy           = 7.858653236e-26
    Injector model1:
      - parcels added               = 1
      - mass introduced             = 4.691445029e-14
    Parcel fate: system (number, mass)
      - escape                      = 0, 0
    Parcel fate: patch INLET (number, mass)
      - escape                      = 0, 0
      - stick                       = 0, 0
    Parcel fate: patch TOP (number, mass)
      - escape                      = 0, 0
      - stick                       = 0, 0
    Parcel fate: patch BOTTOM (number, mass)
      - escape                      = 0, 0
      - stick                       = 0, 0
    Parcel fate: patch OUTLET (number, mass)
      - escape                      = 0, 0
      - stick                       = 0, 0
    Parcel fate: patch WALL (number, mass)
      - escape                      = 0, 0
      - stick                       = 0, 0
    Temperature min/max             = 300, 300
    Mass0 = 0
 dMass = 4.535328818e-16
    Mass transfer phase change      = 4.535328818e-16
    Mass transfer devolatilisation  = 0
    Mass transfer surface reaction  = 0

smoothSolver:  Solving for Ux, Initial residual = 0.03017707843, Final residual = 0.002718135185, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.9966951568, Final residual = 0.06951720496, No Iterations 4
smoothSolver:  Solving for Uz, Initial residual = 0.2575426795, Final residual = 0.01208645105, No Iterations 2


surfaceFieldValue surfaceFieldValue1 write:
    sum(BOTTOM) of H2O = 0
    sum(BOTTOM) of air = 19238

Time = 2


Solving 3-D cloud reactingCloud1

Cloud: reactingCloud1 injector: model1
    Added 1 new parcels

Ni = 1.804549861e-11 d= 4e-06 dt = 1
Cloud: reactingCloud1
    Current number of parcels       = 1
    Current mass in system          = 4.646091745e-14
    Linear momentum                 = (2.536591448e-20 -6.572352722e-19 -1.02175801e-19)
   |Linear momentum|                = 6.656136469e-19
    Linear kinetic energy           = 4.76789473e-24
    Injector model1:
      - parcels added               = 1
      - mass introduced             = 4.691445029e-14
Where as when I put solution method as transient, it gives out multiple iterations for every timestep and it solves for mass transfer multiple in the 1st time step:
Code:
Time = 1

 Entered the transient time loop here 

Solving 3-D cloud reactingCloud1
Ni = 1.80455e-11 0 4e-06dt = 1

Cloud: reactingCloud1 injector: model1
    Added 1 new parcels

Ni = 1.804549861e-11 4.535328818e-16 3.987068597e-06dt = 1
Cloud: reactingCloud1
    Current number of parcels       = 1
    Current mass in system          = 4.600885077e-14
    Linear momentum                 = (8.198791097e-20 -1.215618056e-22 1.352201464e-22)
   |Linear momentum|                = 8.19881126e-20
    Linear kinetic energy           = 7.305171174e-26
    Injector model1:
      - parcels added               = 1
      - mass introduced             = 4.691445029e-14
    Parcel fate: system (number, mass)
      - escape                      = 0, 0
    Parcel fate: patch INLET (number, mass)
      - escape                      = 0, 0
      - stick                       = 0, 0
    Parcel fate: patch TOP (number, mass)
      - escape                      = 0, 0
      - stick                       = 0, 0
    Parcel fate: patch BOTTOM (number, mass)
      - escape                      = 0, 0
      - stick                       = 0, 0
    Parcel fate: patch OUTLET (number, mass)
      - escape                      = 0, 0
      - stick                       = 0, 0
    Parcel fate: patch WALL (number, mass)
      - escape                      = 0, 0
      - stick                       = 0, 0
    Temperature min/max             = 300, 300
    Mass0 = 0
 dMass = 9.055995246e-16
    Mass transfer phase change      = 9.055995246e-16
    Mass transfer devolatilisation  = 0
    Mass transfer surface reaction  = 0

smoothSolver:  Solving for Ux, Initial residual = 0.03017707843, Final residual = 0.002718135185, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 0.9966951568, Final residual = 0.06951720496, No Iterations 4
smoothSolver:  Solving for Uz, Initial residual = 0.2575426795, Final residual = 0.01208645105, No Iterations 2
smoothSolver:  Solving for H2O, Initial residual = 0, Final residual = 0, No Iterations 20
smoothSolver:  Solving for h, Initial residual = 0.9999999996, Final residual = 0.02940408779, No Iterations 2


Time = 2

 Entered the transient time loop here 

Solving 3-D cloud reactingCloud1
Ni = 1.804549722e-11 9.055995246e-16 3.97409498e-06dt = 0.7257109082
Ni = 1.804549621e-11 1.23260168e-15 3.964657458e-06dt = 0.2742890918
Cloud: reactingCloud1
    Current number of parcels       = 1
    Current mass in system          = 4.555854866e-14
    Linear momentum                 = (1.200638905e-19 -1.249845965e-18 -2.45281995e-21)
   |Linear momentum|                = 1.255601963e-18
    Linear kinetic energy           = 1.730231028e-23
    Injector model1:
      - parcels added               = 1
Is there a way to make sure that the diameter of the particle gets updated after every time step according to the dMassPC and make sure that it calculates dMassPC only once every time step?

Thank you
neilk is offline   Reply With Quote

Reply

Tags
liquidevaporation, phasechange, reactingcloud, simplereactingparcelfoam

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
LiquidEvaporationBoil & LiquidEvaporation amirhoshangm OpenFOAM Programming & Development 2 July 31, 2019 06:03
cyclone-case with lower outlet using simpleReactingFoam (steady-state) jango OpenFOAM Running, Solving & CFD 3 April 17, 2018 06:06
LiquidEvaporation and Particle&Wall Interaction liliu OpenFOAM 1 March 6, 2017 10:02


All times are GMT -4. The time now is 09:16.