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

comments to the input parameters and added new parameter d_turb

Browse source
This commit is contained in:
ignis 2020-02-25 09:07:10 +09:00
parent a2ca7c6ba0
commit 1c78049f75
2 changed files with 37 additions and 27 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

63
code/m_parameters.f90
View file

@ -11,6 +11,7 @@ module m_parameters
! Transport Properties
real :: vis,sc,diff
real :: d_turb
real :: scp,prp,lep,vis0p,rod
real :: lewis, le_a, le_x
@ -86,33 +87,41 @@ contains
END IF
OPEN(otape,FILE='otape')
CALL READ_PARAMETER (nx)
CALL READ_PARAMETER (l_0)
CALL READ_PARAMETER (int_pr)
CALL READ_PARAMETER (tar_lo)
CALL READ_PARAMETER (dt)
CALL READ_PARAMETER (sc)
CALL READ_PARAMETER (vis)
CALL READ_PARAMETER (pre)
CALL READ_PARAMETER (ac)
CALL READ_PARAMETER (bc)
CALL READ_PARAMETER (u0)
CALL READ_PARAMETER (tf)
CALL READ_PARAMETER (dt_uf)
CALL READ_PARAMETER (ctmp)
CALL READ_PARAMETER (c_cut)
CALL READ_PARAMETER (c_ref)
CALL READ_PARAMETER (min_wr)
CALL READ_PARAMETER (prof_wr)
CALL READ_PARAMETER (lewis)
CALL READ_PARAMETER (lambda1)
CALL READ_PARAMETER (lambda2)
CALL READ_PARAMETER (beta1)
CALL READ_PARAMETER (hrp)
CALL READ_PARAMETER (le_a)
CALL READ_PARAMETER (le_x)
CALL READ_PARAMETER (nsp)
CALL READ_PARAMETER (nrxn)
CALL READ_PARAMETER (nx) ! n grid points
CALL READ_PARAMETER (l_0) ! domain length
CALL READ_PARAMETER (int_pr) ! print interval
CALL READ_PARAMETER (tar_lo) ! target flame location
CALL READ_PARAMETER (dt) ! time step size
CALL READ_PARAMETER (sc) ! Schmidt Number
CALL READ_PARAMETER (vis) ! viscosity
CALL READ_PARAMETER (pre) ! 1-step chemistry
CALL READ_PARAMETER (ac) ! 1-step chemistry
CALL READ_PARAMETER (bc) ! 1-step chemistry
CALL READ_PARAMETER (u0) ! inlet velocity
CALL READ_PARAMETER (tf) ! end time
CALL READ_PARAMETER (dt_uf) ! flame control time step
CALL READ_PARAMETER (ctmp) ! unburned c value
CALL READ_PARAMETER (c_cut) ! reaction rate profile parameter
CALL READ_PARAMETER (c_ref) ! reaction rate profile parameter
CALL READ_PARAMETER (min_wr) ! reaction rate profile parameter
CALL READ_PARAMETER (prof_wr) ! reaction rate profile parameter
CALL READ_PARAMETER (lewis) ! Lewis number
CALL READ_PARAMETER (lambda1) ! 2-step chemistry
CALL READ_PARAMETER (lambda2) ! 2-step chemistry
CALL READ_PARAMETER (beta1) ! 2-step chemistry
CALL READ_PARAMETER (hrp) ! 2-step chemistry
CALL READ_PARAMETER (le_a) ! 2-step chemistry
CALL READ_PARAMETER (le_x) ! 2-step chemistry
CALL READ_PARAMETER (nsp) ! number of species
CALL READ_PARAMETER (nrxn) ! number of reactions
CALL READ_PARAMETER (d_turb) ! turbulent diffusivity
! READ(itape,*) cdum,table_size
! WRITE(otape,*) cdum,table_size

1
code/sample/itape.py
View file

@ -39,6 +39,7 @@ le_a 1.
le_x 0.20
nsp 2
nrxn 1
d_turb 0
'''
).split()