[Ant_Yang]

Hello all,

I am using the icoReactingMultiphaseInterFoam V2306 and a newer about the VOF. I try to simulate two different liquid droplets, i.e., liquid1 and liquid2, evaporating to the gas phase. When I set two liquid droplets, the vapor of the liquid, i.e., vapor1, occurs in the two locations, see the picture(the same as vapor2) in the attachment. Physically, vapor1 and vapour2 should occur at different locations.

The case I set is as follows:

type    massTransferMultiphaseSystem;

phases  (liquid water gas);

liquid
{
    type            pureMovingPhaseModel;
}

water
{
    type            pureMovingPhaseModel;
}

gas
{
    type            multiComponentMovingPhaseModel;
}

surfaceTension
(
    (gas and liquid)
    {
        type            constant;
        sigma           0.07;
    }
    
    (gas and water)
    {
        type            constant;
        sigma           0.07;
    }
    
    (liquid and water)
    {
        type            constant;
        sigma           0;
    }
);

massTransferModel
(
    (liquid to gas)
    {
        type            Lee;
        species         CH4.gas;
        C               1;
        Tactivate       366;
    }
    
    
    (water to gas)
    {
        type            Lee;
        species         C7H16.gas;
        C               1;
        Tactivate       366;
    }
);

I am not sure if it is a bug or if the solver can not simulate the evaporation of two different liquid droplets. I checked the solver and found the code to calculate the massSpeciesTransfer in multiphaseInter::MassTransferPhaseSystem.C

 forAllConstIters(massTransferModels_, iter)
     {
         if (iter()->transferSpecie() == speciesName)
         {
             // Explicit source
             Su +=
                   this->Su()[phase.name()]
                 + this->Sp()[phase.name()]*phase.oldTime();
         }
     }

No matter what species is transferred from the liquid phase (liquid1 or liquid2), the value of Su is the same.

Anyone can give some comments or suggestions about how to simulate it? Do I set wrong the cases or should I implement the code?

Thanks all,

Ant_