GitLab

I was testing the lagrangian/sprayFoam/aachemBomb tutorial to check if it is possible to spray solid particles instead of fluid particles.

Is the following error a bug or the reason is that the solver accepts only fluid particles sprays?

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v3.0+                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : v3.0+-e941ee6c15e9
Exec   : sprayFoam


...

Starting time loop

Courant Number mean: 0 max: 0
deltaT = 2.94118e-06
Time = 2.94118e-06


Solving 3-D cloud sprayCloud
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigSegv::sigHandler(int) at ??:?
#2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::liquidMixtureProperties::X(Foam::Field<double> const&) const at ??:?
#4  ? at ??:?
#5  ? at ??:?
#6  ? at ??:?
#7  ? at ??:?
#8  ? at ??:?
#9  ? at ??:?
#10  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#11  ? at ??:?
Segmentation fault (core dumped)

Only two modifications done in the tutorial:

1. in the constant/thermophysicalProperties file, comment the C7H16 in the liquid subdictionary, and add C in solid:

liquids
{
/*
    C7H16
    {
        defaultCoeffs   yes;
   
    }
*/
}

solids
{
   C
    {
        defaultCoeffs   no;
        // if  defaultCoeffs no properties should be :
        CCoeffs
        {
            rho             9000;  // Density [kg/m3]
            Cp              462;   // Specific heat capacity [J/(kg.K)]
            K               0.04;  // Thermal conductivity [W/(m.K)]
            Hf              0.312e+06;     // Heat of formation [J/kg]
            emissivity      0.28;  // 1;  // Emissivity  [based on http://www.monarchinstrument.com/pdfs/TableofEmissivity.pdf]
        }
    }
}

2. in the sprayCloudProperties file:

    phaseChangeModel none; //liquidEvaporationBoil;

    singlePhaseMixtureCoeffs
    {
        phases
        (
            solid
            {
                C               1;
            }
        );
    }