forked from jagan-crl/SRADIF
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathInput.py
71 lines (51 loc) · 2.42 KB
/
Input.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
"""
How to Cite:
J.K. Tavares, V. Gururajan, J. Jayachandran, Effects of radiation heat loss on planar and spherical
hydrofluorocarbon/air flames, Combust. Flame 258 (2023) 113067.
"""
import cantera as ct
import numpy as np
import SRADIF as OTL
import time
######################################################################
# General Controls
######################################################################
# Polynomial fit for exp data
POLY = False
######################################################################
# Shell Generation Controls
######################################################################
# Number of layers, default was 1000
NLAY = 1000 #So that polyfit uses odd number
# Chamber volume in cm^3:
Chamber_radius = 3.0e0
VOLU = 4.0e0/3.0e0*np.pi*Chamber_radius**3.0e0
# Ignition shell volume (LTORC data only, set to 0 to ignore)
ign = 0.0e0 #cm
IGN = (4.0e0/3.0e0*np.pi*ign**3.0e0)/7.0e0
######################################################################
# Tolerance Controls
######################################################################
# Volume tolerance in cm^3
ATOL = 1.0E-6
# Gama tolerance
GTOL = 1.0E-6
# Cooling time tolerance in s
TTOL = 1.0E-6
######################################################################
# Gas Initial Conditions
######################################################################
# Initial temperature
TEMP = 300.0e0
# Initial pressure in atm
PRES = 1.0e0
# Reactant composition (moles or mole fraction)
reactants = 'CH2F2:1.2, O2:1.0, N2:3.76'
# Chemical mechanism
scheme = 'hfcmech_R32_2021Aug20.cti'
#######################################################################
startAll = time.time()# program timer
OTL.spherical(reactants,scheme,PRES,TEMP,NLAY,VOLU,\
ATOL,GTOL,TTOL,POLY,IGN)
endAll = time.time() # end time
print(' Elapsed Time {0:4.2f} seconds'.format(endAll-startAll))