DUNS里面错误多多，我快疯了

wangdingxi

calcwi.f里面的50行
         ca(kk) = max(0.,rcl*qr(k+2)/ri(kk))
我认为错了，应该为         ca(kk) = max(0.,rcl*qr(k+2)/wm(kk)*1000)
  ff = ener*dfa(nn)*tcl**dft(nn)*exp(freq*dfe(nn)/tcl)*tbc
         rr = ener*dba(nn)*tcl**dbt(nn)*exp(freq*dbe(nn)/tcl)*tbc
         do k=1,nst
            kk = kmap(k)
            ff = ff*ca(kk)**coff(kk,nn)
            bb = bb*ca(kk)**cofb(kk,nn)
上面的rr应该为bb
还有reactck.f里面的subroutine fun的变量gamdt没有赋值就使用了

origin-right

I had once checked DUNS with its test case, I found it never converge using
H2/O2 reaction scheme. It must have error. I coded my own code with H2/O2
reaction scheme, converge is easily obtained for hypsersonic flow past a
cylinder.

wangdingxi

Did you input the react mechenism of H2/O2 using the reactions mode in DUNS?Did you use a global one step mechanism? Could you change the corresponding codes as what I said and then test it?Thank you for your post.I hope we can have a detailed talk one day.My QQ:63537823

ivor23

留个脚印

菩提逍遥