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

correction of KIVA Itape

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 2, 2006, 03:41
Default correction of KIVA Itape
  #1
hama
Guest
 
Posts: n/a
HELLO:

I'am a new KIVA user I have KIVA 2, and ITAPE also,I want to improve this my ITAPE.

I already made a simple simulation of an internal combustion engine which functions to the natural gas (i.e. I assume that this fuel is100% of methane CH4).

I use the KIVA-II code to calculate the outputs to this simulation and I animate it with TECPLOT software.

According to the (fig 1,2 and3) I note that the flame didn't propagate from the spark plug toward the walls of the cylinder; therefore, the methane/air mixture didn't burn , where I didn't observe any absence of methane and the oxygen after the simulation.

My question is why the spark plug didn't produce a flame that it should be propagates itself on the whole mixture that exists in the combustion chamber.

Or my question with another way is - how I can I get a good propagation of the flame through a methane/air mixture in the combustion chamber-.

If anyone can send me ITAPE or OUTPUT, that'll be very helpful to me.

If anyone want detailles just send and I will do.

thanks

Engine specifications:

Bore 10.47 cm

Stroke =9.53 cm

Clearance=1.26 cm

Compression ratio=8.56

Chamber shape is a disk

RPM =2000 rpm

Fuel : methane (100% CH4)

Equivalent ratio 0.8969

Sparking location : central

ITAPE file:

t3aaa k053190 20x1x22 2-d baseline, w/qsou,chem3190

irest 0

ipost 1

nx 45

ny 1

nz 30

lwall 1

nchop 5

lpr 0

jsectr 1

irez 2

ncfilm 99999

nctap8 99999

nclast 99999

cafilm 45.

cafin 91.0

cadump 9.99e+9

dcadmp 9.99e+9

angmom 1.0

cyl 1.0

dy 0.0

pgssw 0.0

sampl 0.0

dti 1.0e-6

dtmxca 9.99e+9

dtmax 5.0e-5

tlimd 0.0

twfilm 9.99e+9

twfin 9.99e+9

fchsp 0.25

stroke 9.53

squish 1.26

rpm 2000.0

atdc -181.0

conrod 16.34

offset 0.0

swirl 0.0

swipro 0.00

thsect 0.5

epsy 1.0e-2

epsv 1.0e-3

epsp 1.0e-4

epst 1.0e-3

epsk 1.0e-3

epse 1.0e-3

gx 0.0

gy 0.0

gz 0.0

tcylwl 653.0

thead 653.0

tpistn 653.0

tvalve 653.0

tempi 350.0

pardon 0.0

a0 0.0

b0 1.0

anc4 0.05

adia 0.0

anu0 0.0

visrat-.66666667

tcut 700.0

tcute 1200.0

epschm 0.02

omgchm 1.0

tkei 0.3

tkesw 1.0

sgsl 0.0

uniscal 0.0

airmu1 1.457e-5

airmu2 110.0

airla1 252.0

airla2 200.0

expdif 0.2

prl 0.74

rpr 1.11

rprq 1.0

rpre 0.769231

rsc 1.11

xignit 6.0e+3

t1ign -1.425e-2

tdign 1.5e-2

ca1ign-10.0e+0

cadign 10.0

iignl1 1

iignr1 6

jignf1 1

jignd1 1

kignb1 2

kignt1 4

iignl2 0

iignr2 0

jignf2 0

jignd2 0

kignb2 0

kignt2 0

kwikeq 1

modmh 1

amh 10.0

bmh 0.5

cebu 1.

numnoz 0

numvel 0

injdist 0

kolide 0

t1inj 2.5e-2

tdinj 1.40e-3

ca1inj 95.0

cadinj 12.672

tspmas 0.0000

pulse 2.0

tnparc 0.00

rhop 0.0000

tpi 300.0

turb 1.0

breakup 0.0

evapp 0.0

npo 46

nunif 10

1 1 0.000 0.0

2 1 0.116 0.0

3 1 0.233 0.0

4 1 0.349 0.0

5 1 0.465 0.0

6 1 0.582 0.0

7 1 0.698 0.0

8 1 0.814 0.0

9 1 0.931 0.0

10 1 1.047 0.0

11 1 1.163 0.0

12 1 1.280 0.0

13 1 1.396 0.0

14 1 1.512 0.0

15 1 1.629 0.0

16 1 1.745 0.0

17 1 1.861 0.0

18 1 1.978 0.0

19 1 2.094 0.0

20 1 2.210 0.0

21 1 2.327 0.0

22 1 2.443 0.0

23 1 2.559 0.0

24 1 2.676 0.0

25 1 2.792 0.0

26 1 2.908 0.0

27 1 3.025 0.0

28 1 3.141 0.0

29 1 3.257 0.0

30 1 3.374 0.0

31 1 3.490 0.0

32 1 3.606 0.0

33 1 3.723 0.0

34 1 3.839 0.0

35 1 3.955 0.0

36 1 4.072 0.0

37 1 4.188 0.0

38 1 4.304 0.0

39 1 4.421 0.0

40 1 4.537 0.0

41 1 4.653 0.0

42 1 4.770 0.0

43 1 4.886 0.0

44 1 5.002 0.0

45 1 5.119 0.0

46 1 5.235 0.0

nho 0

square 0.0

rcornr 0.0

nstrt 0

icont 11101111011110000000011000

mirror 0

nvzone 0

nvvvec 0

nvpvec 0

nvcont 0

nsp 5

ch4 rho1 4.9653E-5

o2 rho2 1.9861E-4 mw2 32.000 htf2 0.0

n2 rho3 6.5344E-4 mw3 28.016 htf3 0.0

co2 rho4 0.0 mw4 44.011 htf4 -93.965

h2o rho5 0.0 mw5 18.016 htf5 -57.103

rtout 0.0

topout 0.0

botin 0.0

nrk 1

cf1 8.3000e5 ef1 1.5780e+4 zf1 0.0

cb1 0.0 eb1 0.0 zb1 0.0

am1 1 2 0 0 0

bm1 0 0 0 1 2

ae1 1.000 1.000 0.000 0.000 0.000

be1 0.000 0.000 0.000 0.000 0.000

nre 0

if you need to see my chem.f subroutine is:

subroutine chem c c ================================================== ==================== c c calculates the change in species densities and internal energy c due to kinetic chemical reactions c c chem is called by: kiva c c chem calls the following subroutines and functions: (none) c c ================================================== ==================== c

implicit double precision (a-h,o-z)

include 'cvmg.inc'

include 'common.inc'

dimension domega(lnrk)

double precision kf,kb c c <><><><><><><><><><><><><><><><><><><><><><><><><> <><><><><><><><><><> c

modmh=1

cebu=1.0

amh=4.

bmh=0.5 c

tchem=1.0d-10

do 100 k=1,nz

i4b=(k-1)*nxpnyp

do 90 j=1,ny

i4=i4b+(j-1)*nxp+1

do 80 i=1,nx

if(f(i4).eq.0.) go to 80

tijk=temp(i4)

if(tijk.lt.tcut) go to 80

rtijk=1./tijk c++ EBU

roijk=ro(i4)

rtturi4=cebu*eps(i4)/tke(i4) c++

do 70 ir=1,nrk

rp=1.

pp=1.

ne=nelem(ir)

do 20 kk=1,ne

isp=cm(kk,ir)

rom=spd(i4,isp)*rmw(isp)

if(am(isp,ir).eq.0) go to 10

if(rom.le.0.) rp=0.

if(rom.gt.0.) rp=rp*rom**ae(isp,ir)

10 if(bm(isp,ir).eq.0) go to 20

if(rom.le.0.) pp=0.

if(rom.gt.0.) pp=pp*rom**be(isp,ir)

20 continue

kb=0.

kf=0.

teback=1.

teford=1.

ekback=1.

ekford=1.

if(cb(ir).le.0.) go to 30 c +++ c +++ backward reaction coefficient c +++

if(eb(ir).ne.0.) ekback=dexp(-eb(ir)*rtijk)

if(zetab(ir).ne.0.) teback=tijk**zetab(ir)

kb=cb(ir)*teback*ekback

30 if(cf(ir).le.0.) go to 40 c +++ c +++ forward reaction coefficient c +++

if(ef(ir).ne.0.) ekford=dexp(-ef(ir)*rtijk)

if(zetaf(ir).ne.0.) teford=tijk**zetaf(ir)

kf=cf(ir)*teford*ekford c +++ c +++ Magnussen ebu model

if(modmh.eq.1) then

tauchm=mw(2)/(roijk*kf)

if(rtturi4*tauchm.lt.1.)then

kf=amh*rtturi4

rp=dmin1(spd(i4,1)/(am(1,1)*mw(1)),spd(i4,2)/(am(2,1)*mw(2)),

1 bmh*(spd(i4,4)+spd(i4,5))/(bm(4,1)*mw(4)+bm(5,1)*mw(5)))

endif

endif c +++ c +++ if any rate coefficients cannot be put in standard c +++ form, code them by hand and put them here c +++ c +++ find the reference species (the one in greatest danger c +++ of being driven negative) c +++

40 omeg=kf*rp-kb*pp

rmin=0.

if(omeg.le.1.d-30) go to 70

do 50 kk=1,ne

isp=cm(kk,ir)

if(spd(i4,isp).le.0.) go to 50

rom=omeg*fbmam(isp,ir)*mw(isp)/spd(i4,isp)

if(rom.ge.0.) go to 50

if(rom.lt.rmin) iref=isp

rmin=dmin1(rmin,rom)

50 continue

rom=spd(i4,iref)*rmw(iref)

flam=fam(iref,ir)

flbm=fbm(iref,ir)

ctop=flam*kb*pp + flbm*kf*rp

cbot=flam*kf*rp + flbm*kb*pp

domega(ir)=rom*dt*(ctop-cbot)/((rom+dt*cbot)*(flbm-flam))

do 60 isp=1,nsp

spd(i4,isp)=spd(i4,isp)+mw(isp)*fbmam(isp,ir)*dome ga(ir)

60 continue c

dechem=qr(ir)*domega(ir)/ro(i4)

dechk=dabs(dechem/sie(i4))

sie(i4)=sie(i4)+dechem

tchem=dmax1(tchem,dechk) c

civaomega(i4)=domega(ir)/dt c

70 continue

80 i4=i4+1

90 continue 100 continue

write(*,*)ncyc,tchem c +++ c +++ ignition -- single point, dual point, or ring c +++

if(t.lt.t1ign .or. t.gt.t2ign) return

diignl=dfloat(iignl(2))

nspark=cvmgzi(1,2,diignl)

do 140 n=1,nspark

do 130 k=kignb(n),kignt(n)

i4k=(k-1)*nxpnyp

do 120 j=jignf(n),jignd(n)

i4jk=i4k+(j-1)*nxp

do 110 i=iignl(n),iignr(n)

i4=i4jk+i

if(temp(i4).lt.1600.) sie(i4)=sie(i4)*(1.+xignit*dt) 110 continue 120 continue 130 continue 140 continue

return

end

Your help will be very helpful to me.

  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
Floating point exception error Alan OpenFOAM Running, Solving & CFD 10 April 6, 2012 14:02
Correction of KIVA ITAPE hama Main CFD Forum 1 December 2, 2011 14:47
Pressure correction problem richard_larson OpenFOAM Running, Solving & CFD 8 October 30, 2008 09:48
PDF - conserved scalar approach in KIVA vega Main CFD Forum 0 July 5, 2005 02:49
mass flux correction at outflow boundaries Subhra Datta Main CFD Forum 2 November 24, 2003 14:11


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