[viraj20feb]

Hello,


I am trying to simulate heat and mass transfer in two multicomponent stationary mixtures. I realize that in the current OF versions, we cannot define a multicomponent phase as stationary while using reactingEulerFoam solvers. Has anyone tried to implement this or can anyone direct me in the right direction so that I can proceed?


Thank you!

[viraj20feb]

I figured it out by adding a typedef "multiComponentStationaryPhaseModel" in the phaseModels.C file:

 
   typedef
        AnisothermalPhaseModel
        <
            MultiComponentPhaseModel
            <
                InertPhaseModel
                <
                    StationaryPhaseModel
                    <
                        ThermoPhaseModel<phaseModel, rhoReactionThermo>
                    >
                >
            >
        >
        multiComponentStationaryPhaseModel;

    addNamedToRunTimeSelectionTable
    (
        phaseModel,
        multiComponentStationaryPhaseModel,
        phaseSystem,
        multiComponentStationaryPhaseModel
    );

Now, when I try to use this PhaseModel, I get the following error:

reactingTwoPhaseEulerFoam | tee logreact
reactingTwoPhaseEulerFoam: symbol lookup error: /home/viraj/OpenFOAM-dev/platforms/linux64GccDPInt32Opt/lib/libreactingPhaseSystem.so: undefined symbol: _ZNK4Foam20StationaryPhaseModelINS_16ThermoPhaseModelINS_10phaseModelENS_9rhoThermoEEEE19continuityErrorFlowEv

Any help to resolve this error is greatly appreciated.

Thank you!