CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Running, Solving & CFD

Solid particals in gas flow

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Display Modes
Old   June 12, 2013, 07:18
Default Solid particals in gas flow
New Member
Thomas F.
Join Date: Jun 2013
Location: Graz, Austria
Posts: 14
Rep Power: 6
Dommy is on a distinguished road

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. Im 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 (
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!
Dommy is offline   Reply With Quote

Old   June 21, 2013, 04:06
New Member
Thomas F.
Join Date: Jun 2013
Location: Graz, Austria
Posts: 14
Rep Power: 6
Dommy is on a distinguished road
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] Unable to determine global carrier id for phase 0 with local id 0
[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.
FOAM parallel run aborting
[1] #0 Foam::error:rintStack(Foam::Ostream&) in "/app2/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64Gcc46DPOpt/lib/"
[1] #1 Foam::error::abort() in "/app2/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64Gcc46DPOpt/lib/"
[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/"
[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)
            "Foam::label "
                "const label, "
                "const label"
            ") const"
        )   << "Unable to determine global carrier id for phase "
            << phaseI << " with local id " << id
            << abort(FatalError);

    return gid;

Last edited by Dommy; June 21, 2013 at 05:19.
Dommy is offline   Reply With Quote


Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Review: Reversed flow CRT FLUENT 0 July 20, 2012 13:03
Converge problem for unsteady gas solid flow kongl1986 FLUENT 5 March 1, 2012 17:31
gas flow out of a C-D channel ljp FLUENT 0 March 7, 2011 15:32
Dispersed Solid particles in carrier gas (DPM) Frederik FLUENT 1 June 3, 2008 05:13
gas solid pipe flow tahar Phoenics 6 November 3, 2003 13:03

All times are GMT -4. The time now is 00:49.