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

update diffusivity array

Browse source
This commit is contained in:
Yeongdo Park 2019-03-16 15:58:38 -04:00
parent fa1d5bcec8
commit b03880aa48
1 changed files with 6 additions and 0 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

6
code/ysolve.f90
View file

@ -108,6 +108,9 @@
nd=1
CALL d2fnonp(nx,hx,ux(1,:),d2ux(1,:),nd,1)
dm = diff
where (y1(:,1,1)<(1-cstar)) dm = diff * rvis
DO i=1,nx
yr=ux(1,i)
c=1.-yr
@ -386,6 +389,9 @@
CALL dfnonp(nx,hx,ux(1:3,:),dux(1:3,:),3,1)
CALL d2fnonp(nx,hx,ux(3:3,:),d2ux(1,:),1,1)
dm = diff
where (r1(:,1,1)<(1.0-cstar)) dm = diff * rvis
DO i=1,nx
wrate=pre*ux(3,i)*exp(-ac/(1.+bc*(1.-ux(3,i)))) !wrate
! IF ((1.-ux(3,i)).le.c_cut) THEN