ITAPE5文件各参数说明.pdf

KIVA3V 说明文档之 1 ITAPE5文件各参数说明 deck epilog c c epilogue c c c a description of the quantities on the input data file itape5 c in the order in which they appear c c c ci name problem identification line up to 80 characters 课题名称 最多80个字符 ci irest restart flag where 0 no restart 0 dump number to c restart from read in in subr begin used in main driver c kiva and in newcyc which then resets it to zero 重新开始标志 当其值为0时 不重新开始 0时重新开始 作为求解器 需要中断时继续计算 在运行KIVA3V的时候重置为0 这是为了避免重复计算 如可在压缩过程之后输出避免再计算压缩过程 ci nohydro normally 0 set to 1 to turn off all hydro this is c intended for quick runs that have only piston and valve c motion grid rezoning and checks for cells for having c negative volumes or concavity it is a useful check to c ensure that the entire grid and moving surfaces will c behave in an expected fashion before investing a great c amount of time in a full hydro run 通常设为0 设为1时则关闭所有的Hydro 这是为了特意加快在只有一个 活塞和气阀运动时的运算速度 它会重新分区网格并检查是否有负的或 凹单元 这对于确定整个网格和运动表面按照预定规则计算是很有用的 它会在所有hydro运行之前花费大量的时间 ci lwall 1 law of the wall bdry layer 0 freeslip 1 noslip 1表示为壁函数 0 自由滑移壁面 气相速度 壁面速度 1 非滑移壁面 气相法向速度 壁面法向速度 详见KIVA2说明书 ci lpr long print control 1 write on fiche 0 omit entirely 输出控制 1 输出到胶片 0 完全忽略 通常为0 为1时多用来代码调试用 ci irez rezone flag 0 pure Eulerian 1 pure Lagrangian c 2 traditional Kiva engine rezone 3 is the same as 2 c except the squish region is never accordioned thus c preserving the original axial grid spacing irez 4 is c used for rezoning of grids with canted valves where c irez 4 is for quasi symmetric pentroof cylinders c irez 5 is for 2 valve OHV wedge geometries and c irez 6 is for geometries that have both pentroof and c wedge features note that the special rezone coding for c none of these three cases can be considered fully general 再分区标志 0 完全欧拉法 网格节点不动 所以网格形状不变 适用于稳态流 1 完全 拉各朗日法 表示网格节点以流体速度移动 2 传统KIVA法 允许网格改变形状做为活塞 移动计算 以活塞速度 3 除了压缩区域不可折叠外与2一样 因而保留了原有的轴向网格间 距 当其值 4时用来重分区斜向气阀的情况 4用于准平衡的封闭汽缸 5用于2气阀的OHV 锲型几何体 6用于同时有封闭汽缸与楔形体的几何体 注意 如果不是满足以上3种情况的 特殊分区法按一般方法处理 一般选2 ci ncfilm cycle interval between film outputs 图形输出的循环时间间隔 如100个循环 通常为9999表示以其他参数控制输出 如cafilm KIVA3V 说明文档之 2 ci nctap8 cycle interval between restart dumps on file otape8 从重新开始信息转储至otape8的时间间隔 ci nclast cycle number when calculation is terminated 计算终止时的循环数 当其值很大时 把其功能转移到其他参数控制 如cafin twfin 参考值 60000 ci ncmon cycle interval between summary prints of engine cylinder c data on special output files set ncmon 5000 to turn c off this option entirely 特殊输出文件中的汽缸数据进行摘要输出的时间间隔 设为 5000时 将会完全关闭这个选项 ci ncaspec 0 is the number of crank angles listed free format in c the following line s for which film output is to be made c this is useful when one wants output only at a number of c specific crank angles set ncaspec 0 to turn it off c note as dimensioned ncaspec must be 0表示的值为进行图形输出的曲轴转角的个数 这方便我们在特殊角度取出数据 设为0时关闭 注意 维数限定其值必须 atdc ci angmom 1 0 for angular momentum conservation in momflx cyl 1 0 1 0时表示在momflx里角动量守恒 cyl 1 0表示圆柱 0 0表示平面 angmom angmom cyl 见KIVA2 P152 ci pgssw 1 0 to employ pressure gradient scaling 0 0 to disable c note setup turns pgs off if pressure boundaries present 1 0为了使用压力缩放比例梯度 0 0则取消 在内燃机燃烧中 注意 在压力边界线出现时把pgs关闭 典型值0 0 ci dti initial time step dt supplied in input data file dt提供给输入数据文件的初始时间步长 若太大则计算开始时就失败 例如在不复杂的计 算中 如压缩过程 则大些 1 0e 5 在复杂计算中就选小些约为5 0e 7 此参数与引擎速 度和网格大小有关 范围 1 0e 7 5 0e 5 一般用于调试 为本人经验确定 ci dtmxca maximum number of crank angle degrees per time step use a c large value if you don t want it as an upper limit on dt 每时间步长的曲轴转角最大值 如果你不想超出dt极限就选个最大值 如果在KIVA计算失败 后 可调整此参数 更小些 可使计算恢复 在恢复后调大些 如果网格够好 发动机转 速够快时可小些 如果太大则表示用其他参数代替它 范围 0 25 1 0 典型值 0 5 ci dtmax maximum dt allowed during entire calculation 在整个循环中所允许的最大dt值 与前者相同意义 可相互代替 取大值 范围 1 0e 7 0 5e 5 典型值 5 0e 6 ci tlimd 1 0 forces a restart dump exit before job time limit KIVA3V 说明文档之 3 c 0 0 for no restart dump 1 0将在工作时间限度前使重新计算转储库强制退出 表示计算原过程 此值为典型值 0 0则没有此转储库 ci twfilm problem time interval between film outputs 图形输出的课题时间间隔 典型值 9 99e 9 ci twfin problem time when calculation is terminated 计算终止的课题时间 同上 ci fchsp time step factor for kinetic reactions and spray c evaporation if tchem q v fchsp dt is reduced by c the factor fchsp tchem similarly if tevap q v c fchsp dt is reduced by the factor fchsp tevap 动力反应和喷雾蒸发的时间步长因数 如果tchem q v fchsp dt将减少fchsp tchem 同样的 如果tevap q v fchsp dt也减少fchsp tevap 典型值 0 25 ci bore diameter of cylinder 汽缸直径 ci stroke piston travel from bottom dead center to top dead center 活塞冲程 ci squish clearance between cylinder head 0 0是其他情况 ci epsy allowed relative error in implicit solution of mass c diffusion typically 1 0e 3 在解决固有得质量扩散下允许相关错误 典型值 1 0e 3 ci epsv allowed relative error in implicit solution of momentum c diffusion typically 1 0e 3 在解决固有动力扩散下允许相关错误 典型值 1 0e 3 ci epsp allowed relative error in pressure iteration c typically between 1 0e 4 and 1 0e 8 在压力迭代下允许相关错误 典型值取1 0e 4和1 0e 8之间 ci epst allowed relative error in implicit solution of heat c diffusion typically 1 0e 3 在解决固有热传播下允许相关错误 典型值 1 0e 3 ci epskci epsk allowed relative error in implicit solution of tke c diffusion typically 1 0e 3 在解决tke扩散下允许相关错误 典型值 1 0e 3 ci epseci epse allowed relative error in implicit solution of epsilon c diffusion typically 1 0e 3 ci gx body acceleration in x direction or X方向的整体加速度 cm s 2 ci gy body acceleration in y direction or Y方向的整体加速度 ci gz body acceleration in z direction or Z方向的整体加速度 ci tcylwl cylinder wall or liner temperature kelvins this is also c the temp used for l r f d walls in 3 d cartesian meshes 汽缸壁线性温度 这也用来在3维笛卡尔坐标系下的l r f d壁网格温度 可由准维模型得 到 范围 400 600K ci thead cylinder head temperature top boundary of mesh 汽缸顶温度 顶边界网格 范围 450 600 ci tpistn piston face temperature and bottom of 3 d cartesian mesh 活塞表面温度 3维笛卡尔坐标系网格的底部 范围450 650 ci pardonci pardon 1 0 for partial donor cell fluxing see also a0 b0 or c 0 0 for quasi 2nd order upwind fluxing in general qsou c is preferable but use pdc when speed is more important c than accuracy or when cell reynolds numbers are 2 0 c see p 49 of la 11560 ms for more information 1 0是部分供给单元差分法 0 0是为了准2次迎风差分法 一般来说 qsou更可取 但当速度要求比精确性要求更高时或者当单元雷诺数0在上止点后点火 0则相反 ci cadign duration of ignition spark in crank angle degrees 以曲轴转角为单位的点火持续期 详见KIVA2 P81 ci xignl 1 xignr 1 yignf 1 yignd 1 zignb 1 zignt 1 c ignition region may be single cell or a block of cells c xignl1 etc describe the x y z physical coordinates in c cm of some volumetric region and are completely indep c endent of the actual zoning in the ignition region in c the ignition procedure near the end of subroutine chem c cell center coordinates are tested against xignl1 etc c if the cell center falls within the ignition volume the c cell is ignited ignition is very mesh dependent and in c many applications a block of cells that are around the c true spark location must be ignited to achieve a success c ful burn not just a single cell 指点火区域 可能指单一单元或连续的块单元 他们描述的是以厘米计算的在某些测定体积 和完全独立的点火实际区域的XYZ物理坐标 在点火程序中 快到子程序chem尾部 校验其相 反的匹配单元中心 如果单元中心落在点火体积里 则点火 点火完全由网格决定 在许多软 件中 围绕真实点火的块单元必须点火成功 不仅仅是一个单一单元 如果 KIVA3V 说明文档之 7 yignf1 1 yignd ny则表示环形点火 见KIVA2 P81 ci xignl 2 xignr 2 yignf 2 yignd 2 zignb 2 zignt 2 c 2nd dual ignition region optional 二重点火区域 可选 适用于双火化塞设置 ci kwikeq 1 to use fast equilibrium chemistry routine chmqgm c for 6 built in reactions 0 to use slower but general c routine chemeq instead note that chmqgm requires c the presence of carbon and hydrogen 1是用来对于6个嵌入式反应的子程序chmqgm来加快其平衡化学反应 0是用来减慢除了一般的程序chemeq的其他程序 注意 chmqgm需要C和H存在 ci numnoz number of injector nozzles 喷嘴数 如果ny 1 其必为0或1 见kiva2 P11 ci numvelci numvel number of entries in optional injection velocity table 进入可选的喷射速度表的数目 常为1 表示其它喷嘴也是一样的喷射速度 ci injdist 1 inject size distribution about smr i c 0 inject one size only smr i which is interpreted c as the nozzle radius 1 分给关于smr i 分配的喷射尺寸 0 仅一个喷射尺寸 可理解为喷嘴半径 指定油滴是否 有相同的初始半径或一个分配的半径 参考值 0 ci kolide 1 turns on droplet collision calculation 0 bypasses it 1则打开油滴碰撞计算 0则跳过 见KIVA2 P115 ci t1inj time to initiate droplet injection 开始油滴喷射的时间 常为 1 0 ci tdinj duration of injection pulse in seconds 喷射脉冲持续期 S 常为 1 0 ci ca1inj crank angle a t d c to initiate droplet injection c used only if t1inj is less than zero 开始油滴喷射的曲轴转角 仅用在t1inj小于0 0 points in the x c direction 0指向 X方向 is the nozzle radius 喷入cgs的油滴sauter平均半径 令1微米 0 0001厘米 10微米就输入1 e 3 如果injdist 0 则其值为喷嘴半径 Smr smr chf 1 2 ci amp0 initial amplitude of droplet oscillation at injector c based on weber number estimate applicable only if c breakup model is turned on 在喷嘴顶端的起始摆动振幅 基于韦伯数评估 仅适用于打开了破碎模型 ci velinj velocity at injector in cm sec for a continuous spray c or a simple pulse pulse 3 0 a single value of velinj c is supplied for the general case pulse 3 0 velinj c is a table of values that change with time in whatever c manner the user desires before using these values kiva c adjusts them if necessary to ensure that the correct KIVA3V 说明文档之 9 c total fuel mass tspmas will be injected 喷射速度 以cm sec为单位 适用于连续喷射或者单脉冲喷射 pulse 0 0 its value is used as a uniform length c scale when initializing eps in the setup if scli KIVA3V 说明文档之 10 c 0 0 initial eps will instead be proportional to the c distance to the nearest solid wall the code will c ensure that the initial length scale is le sgsl c when subgrid scale turbulence model is used 区域n的初始紊流长度尺寸 如果其值 0 0 当设置eps初始化时它的值被用作 统一的长度 如果其值 0 0 初始eps将代替最近的固壁面相对的距离 此编码 将在subgrid尺寸紊流模型运用时确定初始长度尺寸 0 valve c data follows 移动气阀的数量 0 5 如果其值 0 气阀按以下数据 c ci vliftmin mm minimum clearance of valve from its seat typically c this is 0 05 to 0 1 cm 从阀座离开的气阀的最小余隙 典型值 0 05 0 1cm ci skirtth mm thickness of valve mm at its edge or skirt in cm 气阀裙边的厚度 ci tmove mm temperature of valve mm in kelvins c dimensioned 1 nmovemx because it is referenced c as tmove idface i4 and idface i4 can 1 气阀温度 因为参考tmove idface i4 的值的尺寸 idface i4 1 ci vtiltxz mm tilt angle in degrees of valve mm in x z space c as the camshaft s are assumed parallel to the y c axis vtiltxz 0 is for valves whose stem is canted to the right c positive x direction vtiltxz 0 is for vertical c valves 在 XZ 空间的气阀倾斜角 假定凸轮轴平行于Y轴 其值0则相反 0则是垂直阀 ci nlift mm no of entries valve on file itape18 for each entry c the crank angle and lift in cm are specified only c non zero lift values need specification the crank c angles do not have to be uniformly spaced but their c values must increase from the beginning of the table c to the end the code will calculate the equivalent c crank angle if the current crank angle lies outside c the 0 to degrep 720 or 360 degree range 进入有关气阀的itape18文件次数 对于每次进入都要指定曲轴转角及升程 曲轴 转角没必要均匀分隔 但是他们的气阀则必须从开始到最后的增加 编码将计算其 相应的的曲轴转角 从0到720或者360度 c ci isootci isoot 1 to turn on the soot calculation 0 otherwise c if isoot 1 soot input data is read in following isoot KIVA3V 说明文档之 12 1则打开烟度计算 0为其他 如果其值 1 将从以下读入数据 c ci tcuts if cell temp is lt tcuts bypass soot calculation 如果单元温度小于此值 则跳过此计算 ci drni initial diameter of the growing fuel radical nuclei 初始增长的燃料活化核心直径 ci dnuci initial diameter of the growing soot particle nuclei 初始增长的烟粒子核心直径 ci drncr d r n cr critical diameter of the radical nuclei c chosen to obtain the correct induction time for the fuel 活化核心的临界直径 其用来获得正确的燃料感应时间 ci ra minimum radius of hc molecule e g 2e 8 2 angstroms hc的最小半径 2 e 8 2埃 cm 10 8 ci surf surovikin s coefficient f accounts for chemical c branching in nucleus formation observed experimentally surovikin s系数 为了说明核子形成的化学反应分支 由试验测得 ci rhornci rhorn rho of r n 2 gms cm 3 for carbon C的rho 2 gms cm 3 ci rhopn rho of particle nuclei has a similar value 粒子核的rho 有相似的值 ci xnc number of carbon atoms per fuel hydrocarbon molecule 燃料中的每炭氢分子的C原子的数量 ci ee1 e1 activation energy for the formation of radical nuclei e1 产生活化核心所需的活化能 ci ee2 e2 activation energy for the growth of radical nuclei e2 增长活化核心所需的活化能 ci ee3ci ee3 e3 activation energy for the growth of particle nuclei e3 增长粒子核心所需的活化能 ci ee4ci ee4 e4 activation energy for r n soot interaction e4 r n 烟交互作用所需的活化能 ci eka ekb ekt ekz are the activation energies a b t z and c fka fkb fkt fkz are the corresponding factors for the c empirical rate constants ka kb kt kz in the oxidation c of soot particles nagle c they are converted to ergs gm in subroutine rinput a b t z的活化能 其值是对应着烟粒子氧化时ka kb kt kz的经验率 常数的因素 nagle 在子程序rinput中转化为ergs gm ci fka fkb fkt fkz see eka ekb ekt ekz definition above 如上 c c suggested starting values for soot input data 烟度计算的建议初始值 c c tcuts 1200 0 drni 2 248e 8 dnuci drncr 12 5e 8 ra 2 0e 8 c surf 2 0e 4 rhorn rhopn 2 0 xnc 6 0 e g c6h6 c e1 116 0 e2 30 0 e3 8 0 e4 0 0 c eka 30 0 ekb 15 2 ekt 97 0 ekz 4 1 c fka 20 0 fkb 4 46e 3 fkt 1 51e 5 fkz 21 3 c KIVA3V 说明文档之 13 c ci distamb applies only to outflow boundaries if distamb gt 0 0 c the ambient pressure is imposed a distance distamb beyond c the boundary instead of right on it the purpose is to c enable the boundary to absorb acoustic waves reducing c any tendency they may have to be reflected back in if c used distamb is a typical mesh dimension e g the width 仅运用于出流边界 如果其值大于0 0 在距离边界d