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

sprayEngineFoam dimension inconsistency

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

Reply
 
LinkBack Thread Tools Display Modes
Old   September 25, 2014, 11:15
Default sprayEngineFoam dimension inconsistency
  #1
New Member
 
Join Date: Sep 2014
Posts: 3
Rep Power: 2
andy_andy10 is on a distinguished road
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?
andy_andy10 is offline   Reply With Quote

Old   September 25, 2014, 11:29
Default
  #2
Member
 
Timm Severin
Join Date: Mar 2014
Location: Munich
Posts: 34
Rep Power: 3
Astrodan is on a distinguished road
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
__________________
PhD Student at the Institute of Biochemical Engineering at TU München
Modelling of fluid dynamics in open photobioreactors.

System:
OpenFOAM 2.3.x, 64bit, 8 Core Xeon Workstation
Astrodan is offline   Reply With Quote

Old   September 25, 2014, 12:39
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
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.
mturcios777 is offline   Reply With Quote

Old   September 26, 2014, 06:51
Default
  #4
New Member
 
Join Date: Sep 2014
Posts: 3
Rep Power: 2
andy_andy10 is on a distinguished road
Quote:
Originally Posted by mturcios777 View Post
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 is offline   Reply With Quote

Old   September 26, 2014, 12:09
Default
  #5
New Member
 
Join Date: Sep 2014
Posts: 3
Rep Power: 2
andy_andy10 is on a distinguished road
Quote:
Originally Posted by andy_andy10 View Post
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.
andy_andy10 is offline   Reply With Quote

Old   October 15, 2014, 17:15
Default
  #6
New Member
 
Pang
Join Date: Mar 2011
Location: Denmark
Posts: 25
Rep Power: 6
kmpang is on a distinguished road
Hi,

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

Regards,
Kar
kmpang is offline   Reply With Quote

Old   February 16, 2015, 08:26
Default
  #7
New Member
 
Alex Perez
Join Date: Feb 2015
Posts: 1
Rep Power: 0
alexP is on a distinguished road
Quote:
Originally Posted by kmpang View Post
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?
alexP 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
dimension [0 0 -0 0 0 0 0]? sharonyue OpenFOAM Programming & Development 3 April 10, 2014 03:38
Dimension of pressure: M/LT2 or L2/T2? gcengineer OpenFOAM Running, Solving & CFD 6 August 7, 2013 16:41
ANSYS Fluent 14.0, dimension problem? SteveFinnan ANSYS 2 June 14, 2013 09:12
Tutorial for sprayEngineFoam ed_teller OpenFOAM Running, Solving & CFD 4 May 27, 2013 07:03
the problem about dimension bojiezhang OpenFOAM 2 October 16, 2011 00:22


All times are GMT -4. The time now is 16:26.