CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   Tutorial for sprayEngineFoam (http://www.cfd-online.com/Forums/openfoam-solving/108612-tutorial-sprayenginefoam.html)

ed_teller October 28, 2012 14:52

Tutorial for sprayEngineFoam
 
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...?

Code:


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 November 16, 2012 00:16

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 November 16, 2012 13:28

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 December 26, 2012 00:51

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 May 27, 2013 07:03

Hello ,

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



Cheers,
Vijay.


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