research/dns-hit3d-fdm
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

T^0.76 diffusivity

Browse source
This commit is contained in:
ignis 2019-03-24 13:59:40 +09:00
parent f38b995827
commit 2934fb0e44
1 changed files with 15 additions and 0 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

15
m_fdm_calc.f90
View file

@ -505,6 +505,21 @@ module m_fdm_calc
f_(:,:,:,1) = 0.0 ! continuity
! diffusivity
DO k=1,zz
DO j=1,yy
DO i=1,xx
y=(1.0 - r1_(i,j,k,2)/r1_(i,j,k,1)) * bc + 1.0
dm(i,j,k) = diff * (y ** 0.76)
ENDDO
ENDDO
ENDDO
if (myid.eq.0) write(*,*) 'min(dm)',minval(dm),'max(dm)',maxval(dm)
! reaction source term
DO k=1,zz
DO j=1,yy