[neilk]

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:

    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.

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:

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