[wavefunction]

Hello.

I'm trying to simulate sessile droplet evaporation in a surrounding hotter atmosphere and I am using a custom solver that I found online called interThermalPhaseChangeFoam. I tried to modify the BubbleCondensation tutorial that they have given to make it a droplet evaporation case.

The velocity seems to be attaining ridiculously high values in this case as a result of which the time step is increasing by itself to huge orders.

blockMeshDict:

convertToMeters 1;

vertices
(
//Bot Face
(0 0 -1E-6) //00
(0 0 1E-6) //01
(150E-5 0 1E-6) //02
(150E-5 0 -1E-6) //03
//Top Face
(0 0.002 -1E-6) //04
(0 0.002 1E-6) //05
(150E-5 0.002 1E-6) //06
(150E-5 0.002 -1E-6) //07
);

blocks
(
hex ( 0 3 7 4 1 2 6 5) (210 344 1) simpleGrading (4 1 1) //00
);

edges
(
);

patches
(
patch Bottom
(
( 0 3 2 1 )
)

patch FarStream
(
( 3 7 6 2 )
)

patch Top
(
( 4 5 6 7 )
)

patch Axis
(
( 0 1 5 4 )
)

empty FrontNBack
(
( 1 5 6 2 )
( 0 4 7 3 )
)

);

mergePatchPairs
();

controlDict:

application ASRInterThermalPhaseFoam;

startFrom latestTime;

startTime 0;

stopAt endTime;

endTime 0.15;

deltaT 1E-3;

writeControl adjustableRunTime;
writeInterval 0.002;

purgeWrite 0;

writeFormat binary;

writePrecision 6;

writeCompression uncompressed;

timeFormat general;

timePrecision 6;

runTimeModifiable yes;

adjustTimeStep on;

maxAlphaCo 0.25;
maxCo 0.25;
maxFourier 100;
maxDeltaT 2E-5;;

I think the problem might lie in the boundary conditions I've used as well. I'm not very sure of the boundary conditions that I've given.

alpha1:
https://drive.google.com/open?id=148...ETj_0zNYcuY0KP

H:
https://drive.google.com/open?id=1ty...UfzuFcLEi85pc_

p_rgh:
https://drive.google.com/open?id=1pr...bMkKt-D5XIZpCU

pc:
https://drive.google.com/open?id=1PU...u583AEm8Qg-w3a

T:
https://drive.google.com/open?id=1dU...YxQAKIcKBEdsNW

U:
https://drive.google.com/open?id=1Pg...Ix1ELPC1Gy4I3O

Kindly give me tips on how to proceed solving the problem. Please let me know if there is any additional information I need to provide.