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

Tutorial for sprayEngineFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   October 28, 2012, 14:52
Default Tutorial for sprayEngineFoam
  #1
New Member
 
Rui Alexandre Trigo Ribeiro Pereira
Join Date: Mar 2009
Location: Coimbra, Coimbra, Portugal
Posts: 23
Rep Power: 8
ed_teller is on a distinguished road
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
ed_teller is offline   Reply With Quote

Old   November 16, 2012, 00:16
Default
  #2
u22
Member
 
Anthony Nitski
Join Date: Aug 2009
Location: Earth
Posts: 34
Rep Power: 8
u22 is on a distinguished road
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.
u22 is offline   Reply With Quote

Old   November 16, 2012, 13:28
Default
  #3
Senior Member
 
mturcios777's Avatar
 
Marco A. Turcios
Join Date: Mar 2009
Location: Vancouver, BC, Canada
Posts: 725
Rep Power: 18
mturcios777 will become famous soon enough
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.
mturcios777 is offline   Reply With Quote

Old   December 26, 2012, 00:51
Default
  #4
New Member
 
pushkaraj
Join Date: Nov 2012
Posts: 8
Rep Power: 4
pushkaraj is on a distinguished road
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
pushkaraj is offline   Reply With Quote

Old   May 27, 2013, 07:03
Default
  #5
New Member
 
Vijay Bhaskar Devarapalli
Join Date: May 2013
Posts: 4
Rep Power: 4
vijay480 is on a distinguished road
Hello ,

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



Cheers,
Vijay.
vijay480 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
[ANSYS Meshing] I need tutorial files colorful ANSYS Meshing & Geometry 5 September 19, 2014 01:21
Tutorial for subcooled nucleate boiling Asghari FLUENT 34 September 27, 2013 04:58
Problem on Fluent Tutorial: Horizontal Film Boilig Feng FLUENT 2 April 13, 2013 05:34
[Virtualization] OpenFOAM oriented tutorial on using VMware Player - support thread wyldckat OpenFOAM Installation 2 July 11, 2012 16:01
Rotor/stator tutorial, and how to... gilberto CFX 5 January 21, 2002 10:41


All times are GMT -4. The time now is 06:44.