| ed_teller |
October 28, 2012 13: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 |
