[RobertHB]

Dear Foamers,
i want to simulate the transient advective/diffusive transport of multiple species (Al/Na/Si) through a liquid phase. I'm just missing a suitable solver.

The same can be done with one species in scalarTransportFoam, but the equation does not allow for the use of multiple scalar quantaties.
Everything in multiphase/... can handle multiple liquid phases, but does not seem the simulate the advection/diffusion of dispersed species.


A quick idea that came to my mind, but has to be tested, is running scalarTransportFoam for 1 species, and add the other the via a function object during runtime. Similar what was suggested here. Or even better run icoFoam and add all three species via function objects?


Can you name a suitable solver for this kind of simulation?

With kind regards,
Robert

[RobertHB]

Results from the icoFoam - scalarTransport function object experiments.

Adding a species w/ a function object like it was explained here does nothing at first. Here is my example:

    Aluminium
    {
        type                scalarTransport;
        functionObjectLibs     ("libutilityFunctionObjects.so");
        field               Al;
        active              true;
        autoSchemes         false;
        D                   0.47e-8;
        nCorr               1;
        resetOnStartUp         false;
        fvOptions       
        {
            SourceNa
            {
                type scalarSemiImplicitSource;
                active true;

                scalarSemiImplicitSourceCoeffs
                {
                    volumeMode    absolute;         
                    selectionMode cellSet;
                    cellSet c0;
                 
                    injectionRateSuSp
                    {
                        Al (4.3e-4 0);
                    }
                }
            }
        }
    }

After adding #includeFunc scalarTransport to the controlDict icoFoam can start solving for a species called s. s can not be modified with the code presented above.
But, as soon as #includeFunc scalarTransport is inlcuded in the controlDict icoFoam starts also solving for the species defines by the code. At least the fvOptions source. It does not solve the transport.