CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

manualInjection model in sprayFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   March 27, 2014, 06:42
Default manualInjection model in sprayFoam
  #1
New Member
 
Join Date: Dec 2012
Posts: 12
Rep Power: 4
Mentalo is on a distinguished road
Hi,

i'm trying to use the manualInjection model with the sprayFoam Solver. Unfortunately the simulation crashes prior the first iteration.

This is the corresponding error message:

Code:
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.2.x                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.2.x-0f3a889dd39b
Exec   : sprayFoam
Date   : Mar 27 2014
Time   : 11:35:58
Host   : "philipp-Home"
PID    : 5779
Case   : /home/philipp/cfd/OF_22x/EtF1/full_Spray_EtF1_8
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0


Reading g
Creating combustion model

Selecting combustion model PaSR<psiChemistryCombustion>
Selecting chemistry type 
{
    chemistrySolver ode;
    chemistryThermo psi;
}

Selecting thermodynamics package 
{
    type            hePsiThermo;
    mixture         reactingMixture;
    transport       sutherland;
    thermo          janaf;
    energy          sensibleEnthalpy;
    equationOfState perfectGas;
    specie          specie;
}

Selecting chemistryReader chemkinReader
Reading CHEMKIN thermo data in new file format
chemistryModel: Number of species = 5 and reactions = 1
Selecting ODE solver SIBS
Creating component thermo properties:
    multi-component carrier - 5 species
    liquids - 1 components
    solids - 0 components

Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting turbulence model type RASModel
Selecting RAS turbulence model kOmegaSST
kOmegaSSTCoeffs
{
    alphaK1         0.85034;
    alphaK2         1;
    alphaOmega1     0.5;
    alphaOmega2     0.85616;
    Prt             1;
    gamma1          0.5532;
    gamma2          0.4403;
    beta1           0.075;
    beta2           0.0828;
    betaStar        0.09;
    a1              0.31;
    b1              1;
    c1              10;
    F3              false;
}

Creating field dpdt

Creating field kinetic energy K

No finite volume options present


Constructing reacting cloud
Constructing particle forces
    Selecting particle force sphereDrag
Constructing cloud functions
    none
Constructing particle injection models
Creating injector: model1
Selecting injection model manualInjection
    Constructing 3-D injection
Selecting distribution model RosinRammler
Selecting dispersion model none
Selecting patch interaction model standardWallInteraction
Selecting surface film model none
Selecting U integration scheme Euler
Selecting heat transfer model RanzMarshall
Selecting T integration scheme analytical
Selecting composition model singlePhaseMixture
Selecting phase change model liquidEvaporationBoil
Participating liquid species:
    C2H5OH
Selecting AtomizationModel none
Selecting BreakupModel ReitzDiwakar
Selecting StochasticCollisionModel none
#0  Foam::error::printStack(Foam::Ostream&) in "/opt/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/opt/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2   in "/lib/x86_64-linux-gnu/libc.so.6"
#3  
 in "/opt/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/bin/sprayFoam"
#4  
 in "/opt/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/bin/sprayFoam"
#5  
 in "/opt/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/bin/sprayFoam"
#6  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7  
 in "/opt/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/bin/sprayFoam"
Gleitkomma-Ausnahme (Speicherabzug geschrieben)
My sprayCloudProperties file looks like this:
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.2.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      binary;
    class       dictionary;
    location    "constant";
    object      SprayCloudProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solution
{
    active          true;
    coupled         true;
    transient       yes;
    cellValueSourceCorrection on;

    sourceTerms
    {
        schemes
        {
            rho             explicit 1;
            U               explicit 1;
            Yi              explicit 1;
            h               explicit 1;
            radiation       explicit 1;
        }
    }

    interpolationSchemes
    {
        rho             cell;
        U               cellPoint;
        thermo:mu       cell;
        T               cell;
        Cp              cell;
        p               cell;
    kappa           cell;
    }

    integrationSchemes
    {
        U               Euler;
        T               analytical;
    }
}


constantProperties
{
    T0              294;

    // place holders for rho0 and Cp0 - reset from liquid props using T0
    rho0            1000;
    Cp0             4187;

    youngsModulus   1e9;
    poissonsRatio   0.35;

    epsilon0        1;
    f0              0.5;
    Pr              0.7;
    Tvap            273;
    Tbp             351;

    constantVolume  false;
}


subModels
{
    particleForces
    {
        sphereDrag;
    }

    injectionModels
    {
        model1
        {
        type            manualInjection;  //coneNozzleInjection;
            massTotal       1.25e-4;
            parcelBasisType mass;
            SOI             0;
            positionsFile   "reactingCloud1Positions";
            U0              (0 0 0);
            //injectionMethod disc;
            //flowType        constantVelocity;
        //UMag            0;
            //outerDiameter   2.0e-4;
            //innerDiameter   0;
            //duration        0.1;
            //position        ( 0 0 -0.215 );
            //direction       ( 0 0 1 );
            //parcelsPerSecond 900000;
        //flowRateProfile constant 1;
        //velocityProfile constant 0;
          
  

            //Cd              constant 0.9;

            //thetaInner      constant 0.0;
            //thetaOuter      constant 10.0;

            sizeDistribution
            {
                type        RosinRammler;

                RosinRammlerDistribution
                {
                    minValue        1e-06;
                    maxValue        0.00015;
                    d               0.00015;
                    n               3;
                }
            }
        }
    }

    injectionModel manualInjection;

    dispersionModel none;

    patchInteractionModel standardWallInteraction;

    heatTransferModel RanzMarshall;

    compositionModel singlePhaseMixture;

    phaseChangeModel liquidEvaporationBoil;

    surfaceFilmModel none;

    atomizationModel none;

    breakupModel    ReitzDiwakar; // ReitzKHRT;

    stochasticCollisionModel none;

    radiation       off;

    standardWallInteractionCoeffs
    {
        type            rebound;
    }

    RanzMarshallCoeffs
    {
        BirdCorrection  true;
    }

    singlePhaseMixtureCoeffs
    {
        phases
        (
            liquid
            {
                C2H5OH               1;
            }
        );
    }

    liquidEvaporationBoilCoeffs
    {
        enthalpyTransfer enthalpyDifference;

        activeLiquids    ( C2H5OH );
    }

    ReitzDiwakarCoeffs
    {
        solveOscillationEq yes;
        Cbag            6;
        Cb              0.785;
        Cstrip          0.5;
        Cs              10;
    }

/*
    ReitzKHRTCoeffs
    {
        solveOscillationEq yes;
        B0              0.61;
        B1              40;
        Ctau            1;
        CRT             0.1;
        msLimit         0.2;
        WeberLimit      6;
    }

    TABCoeffs
    {
        y0              0;
        yDot0           0;
        Cmu             10;
        Comega          8;
        WeCrit          12;
    }
*/
}


cloudFunctions
{}


// ************************************************************************* //
Does someone have an idea whats this about? Help is appreciated.

Regards
Mentalo is offline   Reply With Quote

Old   April 2, 2014, 09:29
Default
  #2
Member
 
Ananda Kannan
Join Date: Feb 2014
Location: Horsholm, Denmark
Posts: 44
Rep Power: 3
ansubru is on a distinguished road
Hi!!

If I were you, I would first run the case using the debug build of openfoam. You would get a much better insight into the origin of the error.

Regards

ansubru
ansubru is offline   Reply With Quote

Reply

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
Mixture model + VOF? blacksoil2012 Fluent Multiphase 2 February 9, 2014 00:59
eddy dissipation model: combustion doesn't occur roukaia FLUENT 0 December 25, 2011 14:50
eddy dissipation model: combustion doesn't occur roukaia FLUENT 0 December 24, 2011 10:10
Problems bout CFD model of biomass gasification, Downdraft gasifier wanglong FLUENT 2 November 26, 2009 00:27
Reynolds Stress model in CFX vs Fluent Tim CFX 1 October 7, 2009 06:19


All times are GMT -4. The time now is 08:19.