[deepblue17]

Hi,
I try to setup a simple case using the new MaxwellStefanFourier ThermophysicalTransportModel in multiphaseEulerFoam.

My first step in using the MaxwellStefanFourier model was to use chtMultiRegionFoam. Therefore I added in the thermophysicalTransport file the binary and thermal diffusion coefficients. That worked flawlessly in
chtMultiRegionFoam, and should do also the job in reactingFoam.

My next step was to use the MaxwellStefanFourier model with multiphaseEulerFoam and the bubbleColumnEvaporatingReacting tutorial. But multiphaseEulerFoam doesn't recognise it out of the box. So I figured out, that the models are included in

fluidReactionThermophysicalTransportModels.C

but are missing in:

phaseFluidReactionThermophysicalTransportModels.C

So I added the following lines of code to the latter and recompiled the libso:

#include "FickianFourier.H"
makeLaminarThermophysicalTransportModel(FickianFourier);

#include "MaxwellStefanFourier.H"
makeLaminarThermophysicalTransportModel(MaxwellStefanFourier);

Then at least multiphaseEulerFoam was recognising the MaxwellStefanFourier model, but kicked me out soon after, when it tried to solve for the first component in the gas phase with the following error:

 --> FOAM FATAL ERROR:  
CO.gas not found in table.  Valid entries:  
5 
( 
CO2 
AIR 
H2 
CO 
H2O 
) 


    From function const T& Foam::HashTable<T, Key, Hash>::operator[](const Key&) const [with T = int; Key 
= Foam::word; Hash = Foam::string::hash] 
    in file /home/ubuntu/OpenFOAM/OpenFOAM-9/src/OpenFOAM/lnInclude/HashTableI.H at line 126.

I understand that CO.gas is not a species, but contains only the CO field of the gas phase. So contrary to single phase reactions, the species are stored in CO.gas and CO.liquid files, if both are multicomponent fluids as e.g. multiComponentPhaseModel or reactingPhaseModel.

Any idea how to get the MaxwellStefanFourier also running in multiphases?

[SvenBent]

Hi.

I have just had the same issue with the FickianFourier model and found the following workaround.

The problem is in the lookup of the species as you point out, and it can be traced to the composition.

composition.index(Yi)

here Yi has the phase group added to its name e.g. CO2.gas and the composition only has the specie name CO2. So basically all index calls has to be replaced.

To get the species name without the group the and then the index of the specie in the specieTable find the label as :

const label idx = composition.species()[Foam::IOobject::member(Yi.name())];

Before compiling its important to add the OpenFoam libaray to the Make/options file to be able to use IOobject.

LIB_LIBS = \
     -lOpenFOAM \

Best Regard
Mathias

[bertigo]

Greetings,

I would like to ask you something, since I could not found any information on user guide. Is it possible to use the multi-component transport tool with functions, not with scalar values? If there is any examples of that could you point it for me?

Thank you