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

kivaTest tutorial-fault on temperature (engineFoam)

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

Like Tree5Likes
  • 2 Post By ayhan515
  • 1 Post By ayhan515
  • 1 Post By ayhan515
  • 1 Post By ayhan515

Reply
 
LinkBack Thread Tools Display Modes
Old   March 4, 2010, 06:51
Default kivaTest tutorial-fault on temperature (engineFoam)
  #1
Member
 
Ayhan Eses
Join Date: Mar 2009
Posts: 33
Rep Power: 8
ayhan515 is on a distinguished road
Hi all!

I'm an openfoam beginner. I used kivaTest tutorial (OpenFoam-1,5-dev on Slax). without any modification.

in combustionProperties file "ignite yes;"
Pressure is not increase but temperature is increasing continuously.


I think pressure should have increased to about ~60 -70 bar, temperature should have increased to ~2700-2900 K after TDC decrease again. So i cant find where is problem. Please inform me about what should i do for correct work.

Thanks for your help,
Best Regards.
ayhan515 is offline   Reply With Quote

Old   October 12, 2013, 13:23
Default Pressure Increase
  #2
New Member
 
Apoorv Talekar
Join Date: Oct 2013
Posts: 2
Rep Power: 0
apoorv121 is on a distinguished road
Rescale the pressure in paraFoam, display mode.
apoorv121 is offline   Reply With Quote

Old   August 21, 2015, 16:20
Default
  #3
Member
 
Join Date: Jul 2015
Posts: 31
Rep Power: 2
kaveh19 is on a distinguished road
hi
what's your T dictionary in boundary conditions?
kaveh19 is offline   Reply With Quote

Old   August 21, 2015, 16:25
Default kivaToFoam
  #4
Member
 
Join Date: Jul 2015
Posts: 31
Rep Power: 2
kaveh19 is on a distinguished road
hi everybody
1) I want to run the kivaTest by bore*stroke=76*66 in different compression ratio e.g.,12 and 20. How can i create this C.R ? I mean just by changing GoemetryDictionary I can't fix my favorite compression ratio?
kaveh19 is offline   Reply With Quote

Old   August 21, 2015, 17:34
Arrow
  #5
Member
 
Ayhan Eses
Join Date: Mar 2009
Posts: 33
Rep Power: 8
ayhan515 is on a distinguished road
hi kaveh19,
i suggest you should create new geometry and mesh for different CR. then check your mesh and engine compression ratio with checkMesh & engineCompRatio

wyldckat and kaveh19 like this.
ayhan515 is offline   Reply With Quote

Old   August 25, 2015, 16:08
Default
  #6
Member
 
Join Date: Jul 2015
Posts: 31
Rep Power: 2
kaveh19 is on a distinguished road
Dear ayhan515,
I have modified the below case by your guides:
bore=0.086
stroke=0.086
Fortunately I modified the compression ratio as same as I desired .I received ok for checkMesh .
Vmax = 0.00053768, Vmin = 3.83001e-05
Vmax/Vmin = 14.0386

************************************************** *
My geometry contains 2 cylinder ,one for Vmin and another for Vmax.
I have meshed the geometry in Gambite, but unfortunately I don’t know what boundary condition should I use fore jointed face between Vmin and Vmax volumes .
Some points:
1-wheater or not I consider it Gambite and openFoam will consider it as a wall so the wall appears it paraFoam as I expected .
2-In default case of kivaTest this wall condition doesn’t exist.
So when I want to run the test case .
What's your suggestion? Had you this problem?
Thanks in advance.
kaveh19 is offline   Reply With Quote

Old   August 25, 2015, 16:14
Default
  #7
Member
 
Join Date: Jul 2015
Posts: 31
Rep Power: 2
kaveh19 is on a distinguished road
hi ayhan515,
unfortunately I don't know how i can upload in cfd online .Could you please give me your email in order to I send my test case by details?
kaveh19 is offline   Reply With Quote

Old   August 26, 2015, 08:27
Arrow
  #8
Member
 
Ayhan Eses
Join Date: Mar 2009
Posts: 33
Rep Power: 8
ayhan515 is on a distinguished road
Dear kaveh19,

1.1- create a geometry

just draw single cylinder (not 2 parts)
(for your CR=20 case) diameter: 76 mm height:69,473 mm

1.2-check z direction is height otherwise rotate your geometry.

1.3-rename the bottom surface as "piston"
rename the side surface as "liner"
rename the top surface as "cylinderHead"

2.create a CFD mesh from geometry and check named surfaces.

3.convert your mesh to Foam mesh.

https://yadi.sk/d/5OyB8MYJifqvu

sorry i've never used gambit i couldn't help.

i hope more professional helps come.

best regards.

Quote:
Originally Posted by kaveh19 View Post
Dear ayhan515,
I have modified the below case by your guides:
bore=0.086
stroke=0.086
Fortunately I modified the compression ratio as same as I desired .I received ok for checkMesh .
Vmax = 0.00053768, Vmin = 3.83001e-05
Vmax/Vmin = 14.0386

************************************************** *
My geometry contains 2 cylinder ,one for Vmin and another for Vmax.
I have meshed the geometry in Gambite, but unfortunately I don’t know what boundary condition should I use fore jointed face between Vmin and Vmax volumes .
Some points:
1-wheater or not I consider it Gambite and openFoam will consider it as a wall so the wall appears it paraFoam as I expected .
2-In default case of kivaTest this wall condition doesn’t exist.
So when I want to run the test case .
What's your suggestion? Had you this problem?
Thanks in advance.
kaveh19 likes this.
ayhan515 is offline   Reply With Quote

Old   August 28, 2015, 14:51
Default
  #9
Member
 
Join Date: Jul 2015
Posts: 31
Rep Power: 2
kaveh19 is on a distinguished road
Dear ayhan515,
I have modified the kivaTest for C.R=14. Fortunately I haven't problem with checkMesh,it is ok.In fact I have mentioned this post in other page ,too.
*************************
my engine geometry is:
engineMesh layered;

conRodLength conRodLength [ 0 1 0 0 0 0 0 ] 0.15449;

bore bore [ 0 1 0 0 0 0 0 ] 0.086;

stroke stroke [ 0 1 0 0 0 0 0 ] 0.086;

clearance clearance [ 0 1 0 0 0 0 0 ] 0.00115;

rpm rpm [ 0 0 -1 0 0 0 0 ] 1500;
*****************************
engineCompressioratio is:
engineMesh layered;

conRodLength conRodLength [ 0 1 0 0 0 0 0 ] 0.15449;

bore bore [ 0 1 0 0 0 0 0 ] 0.086;

stroke stroke [ 0 1 0 0 0 0 0 ] 0.086;

clearance clearance [ 0 1 0 0 0 0 0 ] 0.00115;

rpm rpm [ 0 0 -1 0 0 0 0 ] 1500;
****************************
engineCompressionratio is
Vmax = 0.00053768, Vmin = 3.83001e-05
Vmax/Vmin = 14.0386
*************************
the controlDict is:

application engineFoam;

startFrom latestTime;

startTime -180;

stopAt endTime;

endTime 30;

deltaT 0.05;

writeControl runTime;

writeInterval 5;

purgeWrite 0;

writeFormat ascii;

writePrecision 6;

writeCompression off;

timeFormat general;

timePrecision 6;

runTimeModifiable true;

adjustTimeStep true;

maxCo 0.2;

maxDeltaT 1;

functions
{
timeStep
{
type coded;
functionObjectLibs ("libutilityFunctionObjects.so");
redirectType setDeltaT;

code
#{
const Time& runTime = mesh().time();
if (runTime.timeToUserTime(runTime.value()) >= -15.0)
{
const_cast<Time&>(runTime).setDeltaT
(
runTime.userTimeToTime(0.025)
);
}
#};
}
}
**************************************
Other sittings are as same as the default kivaTest.
***************************************
I use mpirun -np 5 engineFoam -parallel >log
****************************************
but the run continues from -180 to -19.9968
*****************************************
my log file is:

Create engine time

Create mesh for time = -180

Selecting engineMesh layered
deckHeight: 0.086
piston position: 0
Reading combustion properties

[4] Found ignition cells:
[4] 1(15267)

Ignition switched off
Reading thermophysical properties

Selecting thermodynamics package
{
type heheuPsiThermo;
mixture inhomogeneousMixture;
transport sutherland;
thermo janaf;
equationOfState perfectGas;
specie specie;
energy absoluteEnthalpy;
}

min(b) = 1

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;
}
--> FOAM Warning :
From function janafThermo<EquationOfState>::limit(const scalar T) const
in file /home/openfoam/OpenFOAM/OpenFOAM-2.4.0/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 108
attempt to use janafThermo<EquationOfState> out of temperature range 298.15 -> 5000; T = -5798.3
--> FOAM Warning :
From function janafThermo<EquationOfState>::limit(const scalar T) const
in file /home/openfoam/OpenFOAM/OpenFOAM-2.4.0/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 108
attempt to use janafThermo<EquationOfState> out of temperature range 298.15 -> 5000; T = -1720.57
.
and
.
.
.
--> FOAM Warning :
From function janafThermo<EquationOfState>::limit(const scalar T) const
in file /home/openfoam/OpenFOAM/OpenFOAM-2.4.0/src/thermophysicalModels/specie/lnInclude/janafThermoI.H at line 108
attempt to use janafThermo<EquationOfState> out of temperature range 298.15 -> 5000; T = -4512.94
DICPCG: Solving for p, Initial residual = 0.30828, Final residual = 4.00653e-18, No Iterations 1
DICPCG: Solving for p, Initial residual = 7.07815e-18, Final residual = 7.07815e-18, No Iterations 0
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 0.0996372, global = -0.0516891, cumulative = 3.35361e+06
DICPCG: Solving for p, Initial residual = 2.20533e-12, Final residual = 2.20533e-12, No Iterations 0
DICPCG: Solving for p, Initial residual = 2.20533e-12, Final residual = 2.20533e-12, No Iterations 0
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 31043.8, global = -0.0516891, cumulative = 3.35361e+06
DILUPBiCG: Solving for epsilon, Initial residual = 0.0415452, Final residual = 1.98579e-21, No Iterations 1
bounding epsilon, min: -1.78742e+63 max: 1.64625e+67 average: 4.90317e+63
DILUPBiCG: Solving for k, Initial residual = 0.996852, Final residual = 2851.55, No Iterations 1001
bounding k, min: -1.33881e+77 max: 6.24989e+75 average: -1.807e+73
Mean pressure:-1.21196e+16
Mean temperature:985.341
Mean u':8.03304e+34
ExecutionTime = 392.13 s ClockTime = 397 s

Courant Number mean: 0.727501 max: 8250.86
deltaT = 1.78964e-28
--> FOAM Warning :
From function Timeerator++()
in file db/Time/Time.C at line 1061
Increased the timePrecision from 11 to 12 to distinguish between timeNames at time -0.0019149
Crank angle = -17.2341 CA-deg
deltaZ = 0
clearance: 0.00245679
Piston speed = 0 m/s
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
*************************************************
Do you have any idea?
thanks in advanc
kaveh19 is offline   Reply With Quote

Old   August 28, 2015, 16:42
Arrow
  #10
Member
 
Ayhan Eses
Join Date: Mar 2009
Posts: 33
Rep Power: 8
ayhan515 is on a distinguished road
Dear kaveh19
my suggestion is decreasing of deltaT

for example

startTime -20;

stopAt endTime;

endTime 30;

deltaT 0.005; //instead 0.05

writeControl runTime;

writeInterval 1;//instead 5


best regards.
ayhan515 is offline   Reply With Quote

Old   Yesterday, 08:46
Default
  #11
Member
 
Join Date: Jul 2015
Posts: 31
Rep Power: 2
kaveh19 is on a distinguished road
Dear ayhan515,
I really appreciate you for your mentions.
******************************************
I want to change the fuel of kivaTest to gasoline . Would you mind telling me how i should start?
In combustionProperties :
fuel IsoOctane; says that fuel is Isooctane , but
GuldersCoeffs
{
Methane
{
W 0.422;
eta 0.15;
xi 5.18;
alpha 2;
beta -0.5;
f 2.3;
}

Propane
{
W 0.446;
eta 0.12;
xi 4.95;
alpha 1.77;
beta -0.2;
f 2.3;
}

IsoOctane
{
W 0.4658;
eta -0.326;
xi 4.48;
alpha 1.56;
beta -0.22;
f 2.3;
}
} specify some weird coefficient ,of course for me , for Methane ,Propane,and Isooctane.
In thermophysicalProperties;
thermodynamics
{
Tlow 298.15;
Thigh 5000;
Tcommon 1000;
highCpCoeffs ( 19.3963 0.0488137 -1.55798e-05 2.17831e-09 -1.08391e-13 -37133.9 -78.8531 );
lowCpCoeffs ( 0.429357 0.0765769 8.41306e-06 -5.86603e-08 2.79414e-11 -30446.4 25.6735 );
}
transport
{
As 1.67212e-06;
Ts 170.672;
} are really unfamiliar for me.
Thanks in advance.
kaveh19 is offline   Reply With Quote

Old   Yesterday, 09:33
Arrow
  #12
Member
 
Ayhan Eses
Join Date: Mar 2009
Posts: 33
Rep Power: 8
ayhan515 is on a distinguished road
Dear kaveh19
i am not an expert and i am also asking questions to professionals to learn something.

if you want to use gasoline you can use isooctane there is no need to change fuel for kivaTest case.

if you want use to a different species or mixture
in my opinion you need
  • Gulders coefficients for laminar flame speed or implement different model for laminar flame speed and,
    coefficients for
  • thermophysical and
  • transport properties

if you have thermo files in chemkin format you can use chemkinToFoam command to get coefficients for thermophysical and transport properties

sorry i haven't used that case with different fuel so i couldn't help.
i hope more professional helps come.


if you use gasoline with 14.0386 compression ratio at SI engine it would be knock.

if you want to use at HCCI conditions then you should check sprayFoam & sprayEngineFoam.

best regards.
kaveh19 likes this.

Last edited by ayhan515; Yesterday at 09:48. Reason: typo
ayhan515 is offline   Reply With Quote

Old   Yesterday, 16:16
Default
  #13
Member
 
Join Date: Jul 2015
Posts: 31
Rep Power: 2
kaveh19 is on a distinguished road
hello ayhan515,
you are right Isooctane has a good harmony with gasoline but If I want to simulate a HCCI case, I can't use kivaTest and should change my solver to sprayEngineFoam?
sorry for bothering you.
bet regards.
kaveh19 is offline   Reply With Quote

Old   Yesterday, 18:41
Arrow
  #14
Member
 
Ayhan Eses
Join Date: Mar 2009
Posts: 33
Rep Power: 8
ayhan515 is on a distinguished road
Quote:
Originally Posted by kaveh19 View Post
hello ayhan515,
you are right Isooctane has a good harmony with gasoline but If I want to simulate a HCCI case, I can't use kivaTest and should change my solver to sprayEngineFoam?
sorry for bothering you.
bet regards.
For HCCI
Just use kivaTest and engineFoam is not sufficient because ignition starts with your control at that case.
HCCI ignition starts with high temperature like diesel engines.
and need to solve many chemical reactions for combustion.
firstly you need to combine
kivaTest (engineFoam) and aachenBomb (sprayFoam) cases and create a new case for sprayEngineFoam and if you define fuel and air species at initial condition you don't have to deal with solving spray that case.

and you can use more detailed reaction mechanisms for HCCI combustion instead of single reaction so your results would be more satisfying and accurate at least in this respect.

best regards.
kaveh19 likes this.
ayhan515 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
Calculation of the Governing Equations Mihail CFX 7 September 7, 2014 06:27
Weird problem with engineFoam tutorial chris1980 OpenFOAM Running, Solving & CFD 29 August 28, 2013 06:53
High temperature methane+air Peter FLUENT 5 January 26, 2009 19:04
Question about the kivaTest tutorial johan OpenFOAM Paraview & paraFoam 1 October 27, 2008 08:07
chemical reaction - decompostition La S. Hyuck CFX 1 May 23, 2001 00:07


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