[andy_andy10]

Hello,

I've been trying to get a sprayEngineFoam example or tutorial just to get started. I've been 90% successful because nothing gets past the following error.

sprayEngineFoam.exe
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.x |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
/* Windows port by CFD support (www.cfdsupport.com) [based on Symscape] *\
\*---------------------------------------------------------------------------*/
Build : 2.3.x-137bb2e4a64c
Exec : C:\cygwin64\opt\OpenFOAM\OpenFOAM-2.3.x\platforms\cygwin64mingw-w64DPOpt\bin\sprayEngineFoam.exe
Date : Sep 25 2014
Time : 16:04:20
Host : "BELL-PC"
PID : 7456
Case : C:/cygwin64/home/Bell/OpenFOAM/2.3.x/run/sprayenginefoam
nProcs : 1
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<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
using integrated reaction rate
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 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

No finite volume options present

Reading/calculating face velocity rhoUf


Constructing reacting cloud
Constructing particle forces
Selecting particle force sphereDrag
Constructing cloud functions
none
Constructing particle injection models
Selecting injection model none
Selecting dispersion model none
Selecting patch interaction model standardWallInteraction
Selecting stochastic collision model none
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
Average parcel mass: nan
Selecting radiationModel none
Courant Number mean: 0 max: 0
Total cylinder mass: 0.000739955

PIMPLE: Operating solver in PISO mode


Starting time loop

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

Solving 3-D 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
:
number of parcels added = 0
mass introduced = 0
Parcels absorbed into film = 0
New film detached parcels = 0
Parcel fate (number, mass)
- escape = 0, 0
- stick = 0, 0
Temperature min/max = 0, 0
Mass transfer phase change = 0
D10, D32, Dmax (mu) = 0, 0, 0
Liquid penetration 95% mass (m) = 0

diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 6.04675e-008, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 7.13203e-008, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 1.6471e-008, No Iterations 2
DILUPBiCG: Solving for C7H16, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for O2, Initial residual = 0.000145828, Final residual = 1.34218e-009, No Iterations 1
DILUPBiCG: Solving for CO2, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for H2O, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for h, Initial residual = 1, Final residual = 3.64821e-007, No Iterations 2
T gas min/max = 300, 433


--> FOAM FATAL ERROR:
LHS and RHS of + have different dimensions
dimensions : [0 3 -1 0 0 0 0] + [1 0 -1 0 0 0 0]


From function operator+(const dimensionSet&, const dimensionSet&)
in file dimensionSet/dimensionSet.C at line 478.

FOAM aborting

I have located my project at (http://expirebox.com/download/52eaac...976ed25df.html) but this dimensional issue suggests an internal solver issue I think. Can anybody help?

[Astrodan]

I can't look at your specific case right now, but looking at the dimension it appears that somewhere a scalar has been set without concentration (factor kg/m^3 is missing on LHS).

If you set some fields/properties yourself I'd suggest you to check those dimensions again. Otherwise you'd have to look at the code to see which fields it tries to evaluate at the moment it fails.

-Timm

[mturcios777]

To further help the debugging, I am running quite a few sprayFoam cases right now and can tell you that the next equation to be solved is that for p. Which means you likely have a kinematic pressure as opposed to a dynamic pressure units in your 0/p file.

[andy_andy10]

Originally Posted by mturcios777

To further help the debugging, I am running quite a few sprayFoam cases right now and can tell you that the next equation to be solved is that for p. Which means you likely have a kinematic pressure as opposed to a dynamic pressure units in your 0/p file.

My pressure units are
[ 1 -1 -2 0 0 0 0 ];

i.e. kg/m/s^2

I have gone through all my input files again based on these idea but no luck.

[andy_andy10]

Originally Posted by andy_andy10

My pressure units are
[ 1 -1 -2 0 0 0 0 ];

i.e. kg/m/s^2

I have gone through all my input files again based on these idea but no luck.

Still no luck, if anybody could take a look at my files and try and run them though that would be great. I will then post this as a diesel engine example case - which will solve this "no place to start" issue for a few people.

That is unless somebody can point to a sprayEngineFoam example which works out of the box.

[kmpang]

Hi,

Have you looked at / compared pEqn.H in both sprayFoam and sprayEngineFoam?

Regards,
Kar

[alexP]

Originally Posted by kmpang

Hi,

Have you looked at / compared pEqn.H in both sprayFoam and sprayEngineFoam?

Regards,
Kar

Hi
I'm having the same problem with sprayEngineFoam. I#m still quite new to CFD and Openfoam so sadly I don't understand what I can do to fix the problem. I have looked at the pEqn.H texts but although I see differences, I am not able to understand them fully and don't know what I should change.
Could you help me with that?