research/incomp-flame-1d
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

added an equation

Browse source
This commit is contained in:
ignis 2019-09-01 03:34:37 +09:00
parent 66adc81ee9
commit d5bb7d388a
1 changed files with 16 additions and 8 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

24
code/ysolve.f90
View file

@ -223,7 +223,7 @@
OPEN (305,FILE='sfield.dat')
DO i=1,nx
WRITE (305,'(e30.20)') y1(i,1)
WRITE (305,'(e30.20)') y1(i,2)
ENDDO
CLOSE (305)
@ -237,7 +237,7 @@
DO i=1,nx
yr=y1(i,1)
theta=y1(i,2)
theta=y1(i,3)
c=1.-yr
ux(1,i)=yr
@ -274,7 +274,8 @@
SUBROUTINE SET_BC
y1(1,1)=1.-ctmp
y1(1,2)=0.
y1(1,2)=1.
y1(1,3)=0.
END SUBROUTINE SET_BC
SUBROUTINE SET_IC
@ -304,8 +305,11 @@
y1(i,1)=(1.-xi) ! reactant mass fraction
ENDDO
! initialize intermediate species field
y1(:,2) = y1(:,1)
! initialize temperature field
y1(:,2) = 1.0 - y1(:,1)
y1(:,3) = 1.0 - y1(:,1)
! calculate flame position
pflame=0.
@ -409,7 +413,7 @@
! x-direction
DO i=1,nx
ux(1,i)=r1(i,1) ! Y
ux(2,i)=r1(i,2) ! T
ux(2,i)=r1(i,2) ! Y
ux(3,i)=r1(i,3) ! T
ENDDO
@ -418,18 +422,22 @@
CALL d2fnonp(nx,hx,ux(:,:),d2ux(:,:),3,1)
DO i=1,nx
wrate=rate_1step(ux(1,i), ux(2,i))
wrate=rate_1step(ux(1,i), ux(3,i))
! - u*dY/dx + D*d2Y/d2x
f(i,1) = - ( u(i)*dux(1,i) ) + diff * d2ux(1,i) - wrate
! - u*dY/dx + D*d2Y/d2x
f(i,2) = - ( u(i)*dux(2,i) ) + diff * d2ux(2,i) - wrate
! - u*dT/dx + D*d2T/d2x
f(i,2) = - ( u(i)*dux(2,i) ) + (lewis * diff) * d2ux(2,i) + wrate
f(i,3) = - ( u(i)*dux(3,i) ) + (lewis * diff) * d2ux(3,i) + wrate
! Boundary conditions
IF (i.eq.nx) THEN
f(nx,1) = -wrate - u(nx)*dux(1,nx)
f(nx,2) = wrate - u(nx)*dux(2,nx)
f(nx,2) = -wrate - u(nx)*dux(2,nx)
f(nx,3) = wrate - u(nx)*dux(3,nx)
ENDIF
ENDDO