[ed_teller]

Hi

I have tried to make a tutorial for sprayEngineFoam, using sprayFoam files for combustion, chemistry and spray definition, and engineFoam for cylinder movement.

I also updated the injecto position description so as to match the cylinde description geometry in engineFoam... everything seems allright, but whan i try to launch sprayEngineFoam, some weird error occurs... can someo ne pls help me here...?

alex@iskandhar:~/OpenFOAM/alex-2.1.x/run/tutorials/lagrangian/sprayFoam/dieselEngine$ sprayEngineFoam
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.1.x                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.1.x-c1ffc2074f88
Exec   : sprayEngineFoam
Date   : Oct 28 2012
Time   : 20:46:59
Host   : "iskandhar"
PID    : 4955
Case   : /home/alex/OpenFOAM/alex-2.1.x/run/tutorials/lagrangian/sprayFoam/dieselEngine
nProcs : 1
sigFpe : Floating point exception trapping - not supported on this platform
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

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

Create mesh for time = -180

Selecting engineMesh layered
deckHeight: 0.085639
piston position: 0

Reading g
Creating combustion model

Selecting combustion model PaSR<psiChemistryCombustionModel>
Selecting psiChemistryModel ODEChemistryModel<gasThermoPhysics>
Selecting thermodynamics package hsPsiMixtureThermo<reactingMixture<gasThermoPhysics>>
Selecting chemistryReader chemkinReader                                                                                                                                                                                          
ODEChemistryModel: Number of species = 15 and reactions = 39                                                                                                                                                                     
Selecting ODE solver SIBS                                                                                                                                                                                                        
Creating component thermo properties:                                                                                                                                                                                            
    multi-component carrier - 15 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 kEpsilon                                                                                                                                                                                          
kEpsilonCoeffs
{
    Cmu             0.09;
    C1              1.44;
    C2              1.92;
    C3              -0.33;
    sigmak          1;
    sigmaEps        1.3;
    Prt             1;
}

Creating field dpdt

Creating field kinetic energy K


Constructing sources
No field sources present


Constructing reacting cloud
Constructing particle forces
    Selecting particle force sphereDrag
Constructing cloud functions
    none
Selecting dispersion model none
Selecting injection model coneNozzleInjection
    Constructing 3-D injection
Selecting distribution model RosinRammler
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:
    C7H16
Selecting AtomizationModel none
Selecting BreakupModel ReitzDiwakar
Selecting StochasticCollisionModel none
Average parcel mass: 3e-06
Selecting radiationModel none
Courant Number mean: nan max: nan
Total cylinder mass: -nan

PIMPLE: Operating solver in PISO mode


Starting time loop

Courant Number mean: nan max: nan
Crank angle = -179.75 CA-deg
deltaZ = 2.86047e-07
clearance: 0.0856387
Piston speed = 0.0102977 m/s

Solving cloud sprayCloud
Cloud: sprayCloud
    Current number of parcels       = 0
    Current mass in system          = 0
    Linear momentum                 = (0 0 0)
   |Linear momentum|                = 0
    Linear kinetic energy           = 0
    Rotational kinetic energy       = 0
    Total number of parcels added   = 0
    Total mass introduced           = 0
    Parcel fate (number, mass)
      - escape                      = 0, 0
      - stick                       = 0, 0
    Mass transfer phase change      = 0
    D10, D32, Dmax (mu)             = 0, 0, -1e+21
    Liquid penetration 95% mass (m) = 0

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG:  Solving for Ux:  solution singularity
DILUPBiCG:  Solving for Uy:  solution singularity
DILUPBiCG:  Solving for Uz:  solution singularity


--> FOAM FATAL ERROR: 
Singular matrix

    From function LUdecompose(scalarSquareMatrix& matrix, labelList& rowIndices)
    in file matrices/scalarMatrices/scalarMatrices.C at line 60.

FOAM exiting

alex@iskandhar:~/OpenFOAM/alex-2.1.x/run/tutorials/lagrangian/sprayFoam/dieselEngine$

Best Regards

Alex

[u22]

Alex, didyou solved a problem? I also try to get sprayEngineFoam working... Can you share your case? My mail is nitskiy.a (at) gmail.com.

[mturcios777]

I would check both the velocity and pressure fields and the timestep. Because you are running multiphysics (fluid + spray), the timestep will need to be smaller. Try setting your starting dT to 0.01, and your maxDeltaT to 0.1. A max Courant number of 0.1 is advisable.

[pushkaraj]

hi ed_teller

can you share your spryengineFoam tutorial with me..... im new here n i have no idea how exactly spray simulation works......i tried sprayFoam its not working out.....can you give me easy algorithm to carry out simulations for spray

thanks
pushkaraj

[vijay480]

Hello ,

Has anyone solved the problems associated with sprayEngineFoam ?
Please reply !



Cheers,
Vijay.