[openfoammaofnepo]

Dear All,

Recently I am working on reactingTwoPhaseEulerFoam solver from OpenFOAM 3.0.1. I have to implement my own mass transfer model, and hence have to read a little bit about the mass transfer model in this solver. We all know that there is a source term in the RHS of the mass fraction governing equation:

  Y1iEqn==*massTransfer[Y1[i].name()]+ fvOptions(alpha1, rho1, Y1[i])

I know the meaning of the highlighted mass transfer term, for example, if we use Frossling model and spherical models, we will have the following lines of codes:

Foam::tmp<Foam::volScalarField>

Foam::massTransferModels::Frossling::K() const

{

    volScalarField Sh(scalar(2) + 0.552*sqrt(pair_.Re())*cbrt(Le_*pair_.Pr()));



    return 6.0*pair_.dispersed()*Sh/sqr(pair_.dispersed().d());

}

Foam::tmp<Foam::volScalarField>

Foam::massTransferModels::sphericalMassTransfer::K() const

{

    return 60.0*pair_.dispersed()/sqr(pair_.dispersed().d());

}

My question is: how *massTransfer[Y1[i].name()] is calculated based on the modelled K?

best regards,
OFFO

[openfoammaofnepo]

In the massTransfer() function, we have the following several lines:

              // Implicit transport through the phase
              *eqns[name] +=
                  phase.rho()*KD*Yf
                - fvm::Sp(phase.rho()*KD, eqns[name]->psi());
  
              // Sum the mass transfer rate
              dmdtExplicit += dmdtSign*phase.rho()*KD*Yf;
              dmdt -= dmdtSign*phase.rho()*KD*eqns[name]->psi();
  
              // Explicit transport out of the other phase
            if (eqns.found(otherName))
              {
                  *eqns[otherName] -=
                      otherPhase.rho()*KD*compositionModel.dY(member, Tf);
              }

So we can see that here dmdt is calculated. Here dmdt is related to the Yf which is the mass fraction at the phase interface. Anybody knows the references for this model?

[openfoammaofnepo]

Dear All,

About the BasePhaseSystem, it is widely used for the class template in the reactingTwoPhaseEulerFoam solver. For example, it is used for the heatTransferModel.H, and the piece of the code is as follows:

/*---------------------------------------------------------------------------*\
                    Class HeatTransferPhaseSystem Declaration
\*---------------------------------------------------------------------------*/

template<class BasePhaseSystem>
class HeatTransferPhaseSystem
:
    public BasePhaseSystem
{
protected:
         
    // Protected typedefs
         
        typedef HashTable
        <
            autoPtr<BlendedInterfacialModel<heatTransferModel> >,
            phasePairKey,
            phasePairKey::hash
        > heatTransferModelTable;

But I am not sure where this "BasePhaseSystem" is from. I searched in the \src and \reactingTwoPhaseEulerFoam, but did not find any relevant results. Do you where where it is defined in the OpenFOAM? Anybody also has this problem and do you have any comment for this question?

[LG Liang]

Hi,

The BasePhaseSystem can be found in reactingTwoPhaseEulerFoam/twoPhaseSystem/twoPhaseSystems.C

Say if the phaseSystem used is interfaceCompositionPhaseChangeTwoPhaseSystem, it can be seen in the code that

typedef
InterfaceCompositionPhaseChangePhaseSystem
<
PhaseTransferPhaseSystem
<
TwoResistanceHeatTransferPhaseSystem
<
MomentumTransferPhaseSystem<twoPhaseSystem>
>
>
>
interfaceCompositionPhaseChangeTwoPhaseSystem;

So the BasePhaseSystem of HeatTransferPhaseSystem would be MomentumTransferPhaseSystem<twoPhaseSystem>. By the way the code is from newer version though. Hope it helps.

Regards,
James