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

diffusivity linear function of c

Browse source
This commit is contained in:
Yeongdo Park 2019-03-19 03:18:20 -04:00
parent df1b7fb30a
commit 363942207e
1 changed files with 2 additions and 4 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

6
code/ysolve.f90
View file

@ -108,8 +108,7 @@
nd=1 nd=1
CALL d2fnonp(nx,hx,ux(1,:),d2ux(1,:),nd,1) CALL d2fnonp(nx,hx,ux(1,:),d2ux(1,:),nd,1)
dm = diff dm = diff * (1.0 + (rvis - 1.0) * (1.0 - y1(:,1,1)))
where (y1(:,1,1)<(1-cstar)) dm = diff * rvis
DO i=1,nx DO i=1,nx
yr=ux(1,i) yr=ux(1,i)
@ -389,8 +388,7 @@
CALL dfnonp(nx,hx,ux(1:3,:),dux(1:3,:),3,1) CALL dfnonp(nx,hx,ux(1:3,:),dux(1:3,:),3,1)
CALL d2fnonp(nx,hx,ux(3:3,:),d2ux(1,:),1,1) CALL d2fnonp(nx,hx,ux(3:3,:),d2ux(1,:),1,1)
dm = diff dm = diff * (1.0 + (rvis - 1.0) * (1.0 - r1(:,1,1)))
where (r1(:,1,1)<(1.0-cstar)) dm = diff * rvis
DO i=1,nx DO i=1,nx
wrate=pre*ux(3,i)*exp(-ac/(1.+bc*(1.-ux(3,i)))) !wrate wrate=pre*ux(3,i)*exp(-ac/(1.+bc*(1.-ux(3,i)))) !wrate