[Dommy]

Hello!

I currently writing my master thesis and therefore have some questions concerning an Euler- Lagrange simulation in OpenFoam 2.1.x.

A gas flow has to be dumped with a bunch of solid particles. I´m interested in the distribution of these particles at the end of the flow. Therefore a Lagrange- Solver has to be used but it is not necessary to take chemical reactions, collision within the particles, surface Films and atomization into account. Coupling between flow and particles is needed but the particles should not effect the flow.

This results in some questions:
Is there a proper way to simulate solid particles within an Euler-Lagrange solver?
I found solidParticleFoam but could not use it (http://openfoamwiki.net/index.php/Co...idParticleFoam).
Also I discovered a library in OpenFOAM-2.1.x/src/lagrangian/solidParticle but have no idea how to implement it.
All multiphase solvers are Euler-Euler and all Lagrangian solvers are made for liquid- solid systems - right?

I looked into all Lagrangian solvers and took a closer look at sprayFoam. Is my idea right that if I turn out all models (chemical reactions, collision within the particles, surface Films, atomization,particle Forces, ,...), set a constant volume and density OpenFoam will treat the particles as solid?

Any other advices or ideas?

Thanks for your time and help in advance!
Regards,
Thomas

[Dommy]

I made some achievements with my simulation.

First I changed to reactingParcelFoam where I was able to inject ash particles in Air - I also found out where and how I can add new materials to the list - until now only C, CaCo3 and ash are implemented.

But there are still some problems.
If I set coupled true; (which is my goal in order to have 2-way coupling and a comparison to experimental results) the solver crashes the first time a particle is injected.

Courant Number mean: 0.0048107 max: 0.0113791
deltaT = 6.79925e-05
Time = 0.000172535

[1]
[1]
[1] --> FOAM FATAL ERROR:
[1] Unable to determine global carrier id for phase 0 with local id 0
[1]
[1] From function Foam::label Foam::CompositionModel<CloudType>::localToGlobalCa rrierId(const label, const label) const

Solving cloud reactingCloud1
[1] in file /home/heinri8/OpenFOAM/OpenFOAM-2.1.x/src/lagrangian/intermediate/lnInclude/CompositionModel.C at line 262.
[1]
FOAM parallel run aborting
[1]
[1] #0 Foam::error:rintStack(Foam::Ostream&) in "/app2/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64Gcc46DPOpt/lib/libOpenFOAM.so"
[1] #1 Foam::error::abort() in "/app2/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64Gcc46DPOpt/lib/libOpenFOAM.so"
[1] #2
[1] in "/app2/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64Gcc46DPOpt/bin/reactingParcelFoam"
[1] #3
[1] in "/app2/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64Gcc46DPOpt/bin/reactingParcelFoam"
[1] #4
[1] in "/app2/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64Gcc46DPOpt/bin/reactingParcelFoam"
[1] #5
[1] in "/app2/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64Gcc46DPOpt/bin/reactingParcelFoam"
[1] #6
[1] in "/app2/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64Gcc46DPOpt/bin/reactingParcelFoam"
[1] #7
[1] in "/app2/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64Gcc46DPOpt/bin/reactingParcelFoam"
[1] #8 __libc_start_main in "/lib/libc.so.6"
[1] #9 Foam::regIOobject::writeObject(Foam::IOstream::str eamFormat, Foam::IOstream::versionNumber, Foam::IOstream::compressionType) const in "/app2/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64Gcc46DPOpt/bin/reactingParcelFoam"
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 1 in communicator MPI_COMM_WORLD
with errorcode 1.

NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.
--------------------------------------------------------------------------
--------------------------------------------------------------------------
mpirun has exited due to process rank 1 with PID 6970 on
node caepc110 exiting without calling "finalize". This may
have caused other processes in the application to be
terminated by signals sent by mpirun (as reported here).
--------------------------------------------------------------------------

Any ideas?

If I set coupled false; the simulation runs fine, but as soon as a particle touches a wall it stops without any errors (standardWallInteration rebound). It runs fine with escape or stick...



EDIT: This is the point where the montioned error apears (from CompostionModel.C)

template<class CloudType>
Foam::label Foam::CompositionModel<CloudType>::localToGlobalCarrierId
(
    const label phaseI,
    const label id
) const
{
    label gid = phaseProps_[phaseI].globalCarrierIds()[id];

    if (gid < 0)
    {
        FatalErrorIn
        (
            "Foam::label "
            "Foam::CompositionModel<CloudType>::localToGlobalCarrierId"
            "("
                "const label, "
                "const label"
            ") const"
        )   << "Unable to determine global carrier id for phase "
            << phaseI << " with local id " << id
            << abort(FatalError);
    }

    return gid;
}