下载文档原格式
溶胶的流变性质所谓流变性质，是指物质在外力作用下的变形和流动性质。
纯液体和大多数小分子溶液，在层流条件下其切应力和切变速率成正比τ＝ηＤ。
这就是著名的牛顿公式，式中的η是比例常数，称为黏度，国际单位是Pa·S。
η代表物质在流动时内摩擦力的大小。
液体流动时，为克服摩擦阻力而消耗一定能量（转变成热），所以从能量角度看，黏度也可定义为单位剪切速率下，单位体积和单位时间内所消耗的能量。
凡符合牛顿公式的流体称为牛顿流体，其特点是黏度只与温度有关，与切边速率无关。
不符合牛顿公式的流体称非牛顿流体，τ与Ｄ无正比关系，比值τ／Ｄ不是常数，而是Ｄ的表示，此时的τ／Ｄ，称为表观黏度。
事实上大多数的流体是非牛顿流体，物函数，用ηａ料随着剪切率或时间的变化会改变。
因此，在一定的条件下测量的粘度值不一样，所测得的粘度值是个曲线而不是一个恒定的常数层流和湍流体系的流动处于恒稳状态，体系中任一点的流速（大小和方向）不随时间而变化时，称为层流。
牛顿公式只有在层流状态下才成立。
当流速超过某一限度时，有不规则的或随时间而改变的旋涡生成，这时称为湍流。
流体的流动状态可以用雷诺数（Ｒｅ）来表示，它是一个无量纲量。
Ｒｅ超过某一临界值时，层流就变成湍流。
稀分散体系的粘度：液体流动时，为克服内摩擦阻力需消耗一定的能量，倘若液体中存在粒子，则流体的流线在粒子附近受到干扰，这就要消耗额外的能量，所以分散体系的黏度均高于纯溶剂的黏度。
通常溶胶黏度η与溶剂黏度η0之比称为相对黏度ηｒ。
1 分散相浓度的影响einstein黏度定律1906年，einstein根据流体力学理论推导出稀分散体系的黏度方程η＝η0（1＋2.5φ）式中，η、η0分别为溶液、溶剂的粘度，φ为体系中分散相的体积分数。
在推导此式时曾假设：①粒子是远大于介质分子的圆球；②粒子是刚体，完全为介质润湿；③分散体很稀，粒子间无相互作用；④无湍流。
对于较浓分散体系，粒子间相互干扰，Ｅｉｎｓｔｅｉｎ公式不再适用。
分类rcs曲线
RCS（雷达截面，Radar Cross Section）曲线是描述目标对雷达波的散射特性的图形。
分类RCS曲线通常是根据目标的形状和尺寸，以及其在不同方向上的散射特性而进行的。
以下是一些常见的RCS曲线分类：
1.球面散射体（Spherical Scattering）：
-球体是一种常见的散射目标，其RCS曲线通常在所有方向上均匀分布。
球体的RCS与入射波的方向无关，是一种各向同性散射。
2.角柱形散射体（Corner Reflector）：
-角柱形散射体是一种常见的RCS曲线，在某些方向上具有明显的峰值，而在其他方向上较低。
这是由于角柱形散射体对入射波的某些方向上的反射特性而产生的。
3.平板散射体（Flat Plate）：
-平板是一种具有明显后向散射特性的目标。
其RCS曲线在后向方向上通常会有一个峰值，而在其他方向上较小。
4.非均匀散射体（Non-uniform Scattering）：
-有些目标具有复杂的形状，使得其RCS曲线在不同方向上表现出不同的特征。
这些目标的RCS曲线被称为非均匀散射体，因为其散射特性在各个方向上不均匀分布。
5.各向同性散射体（Isotropic Scattering）：
-一些目标在所有方向上的散射特性相对均匀，其RCS曲线呈现出均匀分布的特征。
这类目标被称为各向同性散射体。
这些分类有助于理解不同形状和结构目标在雷达波照射下的散射特性，对于雷达系统性能评估和目标识别具有重要意义。
聚合物溶液的流变模式聚合物流变性是指其在流动过程中发生变形的性质，主要体现在有外力场作用时，溶液粘度与流速或压差之间的变化关系。
高分子的形态变化导致了聚合物溶液宏观性质的变化。
聚合物溶液是非牛顿流体，其流动行为可用Ostwald-Dewael幂律方程[56]来描述：根据流体力学对液体流型的分类，驱油用的部分水解聚丙烯酰胺溶液属于假塑性流型，即表观粘度呈现剪切稀化现象。
在低剪切速率下，溶液的流变曲线斜率n=1，符合牛顿流动定律，称为第一牛顿流动区，该区的粘度通常称为零切粘度η0，即γ&→0的粘度。
随着剪切速率的增大，流动曲线的斜率n<1，称假塑性区，该区的粘度为表观粘度ηa；剪切速率的增大，表观粘度ηa值变小，其表观粘度与剪切速率呈幂指数关系。
在假塑性区，剪切速率与表观粘度的关系可用Ostwald-Dewael幂律模型来描述。
HPAM水溶液的完整的流变曲线如图2-1所示。
流变曲线包括牛顿段、假塑段、极限牛顿段、粘弹段和降解段。
粘度随剪切速率的变化与高分子在溶液中的形态结构有关。
在很小的剪切速率下，大分子构象分布不改变，流动对结构没有影响，聚合物溶液的粘度不随剪切速率的变化而变化，此即牛顿段；当剪切速率较大时，在切应力的作用下高分子构象发生了变化，长链分子偏离平衡态构象，而沿流动方向取向，使聚合物解缠和分子链彼此分离，从而降低了相互运动阻力，这时表观粘度随剪切速率的增加而降低。
当剪切速率增加到一定程度以后，大分子取向达到极限状态，取向程度不再随剪切速率而变化，聚合物溶液遵守牛顿流动定律，表观粘度又成为常数，此即所谓的极限牛顿段。
当剪切速率再增加时，主链的相邻键偏离了正常的键角，从而产生了弹性恢复力，而表现出粘弹性，使表观粘度增加。
当剪切速率增加到足以使高分子链断裂时，发生能了聚合物降解，使聚合物粘度降低。
另外，随着剪切速率的增加，ηa下降，开始时ηa下降很快，随后变得缓慢，逐渐趋于平缓，进一步说明了HPAM属于假塑性流体。
化学纤维异形度指标的表征方法
化学纤维异形度指标（abbreviated as CFA)是衡量特定拉维分散体系中纤维异形度的常用
指标。
它是依据表面几何形态学原理，基于研究流变学特性，或者基于分散体系改良粉状溶液稳定性的基础上，来测量纤维的分布状况的指标。
化学纤维异形度指标是通过采用一系列测定参数来建立的，不同的参数会给出不同的纤维异形度指标值。
其主要参数包括：细度（用布尔系数和颗粒度指数来衡量），粒径（采用均匀和不均匀度指标），形貌（粒度统计、微结构、长纤维比例和极限纤维长度）。
化学纤维异形度指标的表征方法有两种：一种是以颗粒度指数增加，以粒度均匀度下降的模式来表示；另一种是以其余参数的增加或下降来表征。
以异形度从低到高为例，以前者表征则为：颗粒度指数增加，粒度均匀度下降；以后者表征则为：细度降低，粒径增大，形貌上极限纤维长度增加，数量分布发生剧变等。
因此，化学纤维异形度指标可以有效地表征不同拉维分散体系中纤维形貌和粒度结构特征，为拉维分散体系的配方优化和性能预测提供依据，从而更好地提高拉维分散体系的稳定性。
非牛顿型流体的分类Company Document number：WTUT-WT88Y-W8BBGB-BWYTT-199984． 非牛顿型流体的分类非牛顿型流体是一大类实际流体的统称。
一般地说，凡流动性能不能用方程（2-2）来描述的流体，统称为非牛顿型流体。
在高分子液体范畴内，可以粗略地把非牛顿型流体分为：纯粘性流体，但流动中粘度会发生变化，如某些涂料、油漆、食品等。
粘弹性流体，大多数高分子熔体、高分子溶液是典型的粘弹性流体，而且是非线性粘弹性流体。
一些生物材料，如细胞液，蛋清等也同属此类。
流动性质有时间依赖性的流体。
如触变性流体，震凝性流体。
4．1 Bingham 塑性体Bingham的可塑性质。
只有当外界施加的应力超过屈服应力y σ，物体才能流动。
流动方程为：⎩⎨⎧≥-<=yy y σσησσσσγ/)(0 （2-74）说明：有些Bingham 塑性体，在外应力超过y σ开始流动后，遵循Newton 粘度定律，流动方程为：γησσ p y += （2-75） 称为普通Bingham 流体，p η为塑性粘度。
有些Bingham 塑性体，开始流动后，并不遵循Newton 粘度定律，其剪切粘度随剪切速率发生变化，这类材料称为非线性Bingham 流体。
特殊地，若流动规律遵从幂律，方程为n y K γσσ += （2-76） 则称这类材料为Herschel-Bulkley 流体。
图2-16 Bingham 流体的流动曲线牙膏、油漆是典型Bingham 塑性体。
油漆在涂刷过程中，要求涂刷时粘度要小，停止涂刷时要“站得住”，不出现流挂。
因此要求其屈服应力大到足以克服重力对流动的影响。
润滑油、石油钻探用泥浆，某些高分子填充体系如碳黑混炼橡胶，碳酸钙填充聚乙烯、聚丙烯等也属于或近似属于Bingham 流体。
填充高分子体系出现屈服现象的原因可归结为，当填料份数足够高时，填料在体系内形成某种三维结构。
如CaCO 3形成堆砌结构，而碳黑则因与橡胶大分子链间有强烈物理交换作用，形成类交联网络结构。
板面为各向异性面的横观各向同性弯曲板的精化理论李荣荣;赵宝生【摘要】In this paper, based on the transversely isotropic theory, the refined theory of a transversely isotropic bending plate with the anisotropic plane is studied and derived. By using Elliott-Lodge solution and transverse-ly isotropic elastic theory, the displacement and stress components were obtained in terms of the mid-plane dis-placement and its derivatives without any assumptions. The refined equation for the plate with transverse sur-face loadings is obtained by using Lur'e method and boundary conditions and then omitting the high orders, the approximate governing differential equation is obtained. Making all physical quantities are irrelevant to x2 or x3 and the results are consistent with the refined theory of transversely isotropic beams and isotropic beams respectively.%在板面为各向同性面的横观各向同性板精化理论的基础上,对板面为各向异性面的横观各向同性板进行了研究,并推导出其精化理论.根据横观各向同性弹性理论和Elliott-Lodge通解,在不作任何预先假设的条件下,获得了由板中面上的位移和转角表示的位移场和应力场,根据Lur'e方法和边界条件获得了板面受横向载荷的精化方程,略去高阶项后获得可直接应用的近似控制微分方程.令所有物理量与x2或x3无关,得出的精化理论分别与横观各向同性梁和各向同性梁的精化理论一致.【期刊名称】《辽宁科技大学学报》【年(卷),期】2015(038)004【总页数】6页(P268-273)【关键词】弹性板;横观各向同性;精化理论;Elliott-Lodge通解;Lur'e方法【作者】李荣荣;赵宝生【作者单位】辽宁科技大学机械工程与自动化学院,辽宁鞍山 114051;辽宁科技大学机械工程与自动化学院,辽宁鞍山 114051【正文语种】中文【中图分类】O343对于三维板问题的分析，本文对沿板厚方向的位移等物理量做出适当的假设，从而简化成二维问题。
姓名：高墨尧学号：20150614专业：农业机械化非牛顿流体的分类根据非牛顿流体的粘度函数是否和剪切时间有关，可以把非牛顿流体分成两大类：非时变性非牛顿流体和时变性非牛顿流体。
1、非时变性非牛顿流体这类流体的切应力仅与剪切速率有关，即粘度函数仅与应变速率或（切应力）有关，而与时间（1（2与γ 的变体）（1（2τ处相（1B交。
τ时，宾汉姆流体并不产生流（2）当对流体施加的外力τ<τ时，体系才产生流动。
动，体积只产生有限的变形，只有当τ>Bτ是使体系产生流动所需要的最且流动后流体具有剪切稀释性。
B小剪切应力，即使流体产生大于0的剪切速率所需要的最小剪切应力，称之为屈服值。
屈服值的大小是体系所形成的空间网络结构的性质所决定的。
凡是具有屈服值的流体均称为塑性流体，外力克服其屈服值而产生的流动称为塑性流动。
2、时变性非牛顿流体这类流体的粘度函数不仅与应变速率有关,而且还与剪切持续时间有关。
大致可分为两类：触变性和流凝性流体：随着切应力作用时间的延长，表观粘度越来越小的流体叫做触变性流体随着切应力作用时间的延长，表观粘度越来越大的流体叫做流凝性流体，这种流体在实际中非常少见。
其特点：（1（2（3（4（1（2（3（4（5（6[1][2][3]刘海燕,庞明军,魏进家.非牛顿流体研究进展及发展趋势[J].应用化工,2010,05:740-746.[4]施庆珊,王计伟,欧阳友生,陈仪本.非牛顿流体粘度测定方法研究进展[J].发酵科技通讯,2011,02:42-45.[5]刘晓明,艾志久,黄俭波,陶云,熊昕.非牛顿流体与牛顿流体在旋流器内的流场分析[J].石油机械,2009,03:28-31.。
不饱和聚酯的流变性能如何测试不饱和聚酯是一种常见的聚合物材料，其在工业、建筑、汽车等领域中得到广泛应用。
不饱和聚酯的流变性能对其使用性能和耐久性有重要影响，因此需要进行相应的测试和研究。
一、不饱和聚酯的流变性能不饱和聚酯是指在树脂中加入不饱和单体（如乙烯基苯等）和交联剂（如双羰基化合物等）进行聚合反应得到的一种聚合物。
不饱和聚酯具有良好的成型性和涂装性能，因此广泛应用于制造复合材料、FRP、建筑材料、汽车外壳等领域。
不饱和聚酯的流变性能主要包括粘度、流动性、弹性模量、屈服强度等参数。
其中，粘度是指不饱和聚酯在一定条件下的黏度值，通常使用旋转黏度计测试。
流动性是指不饱和聚酯在一定条件下的流动性能，通常使用流变仪测试。
弹性模量是指不饱和聚酯在受力下的应变与应力的比值，属于材料的弹性机械性能之一。
屈服强度是指不饱和聚酯在受力下达到屈服的最大强度。
二、不饱和聚酯的流变性能测试方法1. 旋转黏度法旋转黏度法是一种常用的粘度测试方法，可用于测定不饱和聚酯的粘度值。
该方法使用旋转黏度计，将待测样品置于旋转锥或圆柱中，通过测量旋转锥或圆柱受到的扭矩和转速的比值得到粘度值。
旋转黏度法的优点是简单易行，精度高，能够测量多种不同粘度范围的样品。
但是，该方法需要待测样品具有一定的流动性，且需要较长时间的等温静置，所以不适用于测量高粘度和非牛顿流体。
2. 流变仪法流变仪法是一种流动性测试方法，可用于测定不饱和聚酯的流变性能。
该方法使用流变仪，通过施加旋转或振荡等外力，测量材料产生的应力响应和变形响应，从而得到流变性能参数。
流变仪法的优点是能够测量大范围的流动性，能够测量非牛顿流体和高粘度样品。
此外，流变仪法还能够测量不饱和聚酯的弹性模量、屈服强度等参数，从而全面了解材料力学性能。
3. 其他测试方法此外，还有其他测试方法可用于测定不饱和聚酯的流变性能，如拉伸、剪切、压缩等力学性能测试。
这些测试方法能够对不饱和聚酯的机械性能进行全面评估，从而为材料设计和工程应用提供重要参考。
催化剂负载量对三元催化转化器转化效率的影响孔祥华【摘要】According to the light-off behavior of three-way catalytic converter, the mathematical model has been developed. Mathematical model was solved by finite difference. It was found that the conversion efficiency was increased by coupling standard loading (30 000 m2/m3) when catalyst distribution was uniform. The conversion efficiencies of 4 times standard and 2 times standard were almost the same. For a fixed amount of metal loading, when catalyst distribution was nonuniform.the conversion efficiency was decreased with more catalyst near the outlet and increased with more catalyst near the inlet.%针对三元催化转化器的起燃特性,建立了转化器的数学模型,给出数值求解方法,模拟结果表明:以30 000 m2/m3为基准负载,催化剂均匀分布时,负载量增加2倍,转化率有所增加,但4倍基准负载与2倍基准负载的转化率变化不明显,在催化剂总量不变的情况下,沿轴向递增,转化率有所下降,沿轴向递减,起燃时间缩短,转化率增大.【期刊名称】《河南科学》【年(卷),期】2011(029)012【总页数】4页(P1495-1498)【关键词】催化转化器;数学模型;转化效率;催化剂负载量【作者】孔祥华【作者单位】宝鸡文理学院陕西省灾害监测与机理模拟重点实验室,陕西宝鸡721013【正文语种】中文【中图分类】TK411.5汽油发动机尾气中的有害气体通过安装三元催化转化器（Three－Way Catalytic Converter，TWC）被转化为无害的CO2、N2和H2O而得到有效控制.虽然现在汽车上的TWC在催化活性温度范围内对CO、HC和NOx的转化率已达95%以上，但在起燃温度点以前却不能有效地发挥作用[1].在装有TWC的汽车上，发动机冷启动后的前几分钟内CO和HC的排放率为60%～80%[2-3].这样的排放量难以满足日益严格的排放法规[4].为了降低冷启动排放量，就必须尽可能缩短起燃时间，提高转化率.但是在缺少足够理论分析与指导的情况下，转化器内复杂的物理化学过程和大量设计运行参数增加了实际模拟难度.已开发的模型大多集中在转化器内的瞬态行为分析[5-7]，传热传质过程分析[8-9]、催化过程动力学分析[10-11]、模型的建立都是在催化剂均匀负载模式下进行的.本文在考虑所有气体化学反应的情况下，假定进气温度和进气浓度恒定，建立数学模型，准确预测催化剂负载量对转化效率的影响，为催化器的优化设计提供参考.1．1 催化转化器数学模型催化转化器模型建立的基本假设为：①反应器表面的性质均匀；②与对流传热相比，气相的辐射传热与热传导可以忽略；③催化转化器具有良好的绝热性能；④所有小通道内的气体温度与浓度相等.通过以上的假设，催化转化器模型基本上就是转化器内气相与固相的质量与能量的平衡问题了.式中：δ为催化剂的空隙率，%；h为导热系数，W/（m·K）；S为催化剂几何表面积，m2/m3；Cgj和C sj分别表示第j种反应物在气相中的浓度与在固相表面上的浓度（j＝1～6，1 表示 CO；2 表示 C3H6；3 表示 CH4；4 表示 H2；5 表示NO；6表示O2）；hDj为气体物质j的传质系数；S cat为催化剂的活性表面积，m2/m3；Rj为第j种物质的反应速率，mol/（m2·s）；M为平均摩尔质量，mol/kg；ρs为催化剂表层密度，kg/m3；Tg，Ts分别表示气相和固相的温度；ks为催化剂表面的热导率，W/（m·K）；ΔHj为j物种的反应热，J/mol.模型的边界条件为：传热系数采用Baba[12]等的结果：式中：传质系数（hDj）为式中：所用的物质j的扩散系数Dj由Slattery－Bird[13]公式计算：其中：a＝2．745×10－4，b＝2．334.1．2 数值求解模型有偏微分方程、常微分方程，还有代数方程，因此找到一种合理的解法对求解数学模型很重要.对气体能量和质量方程，因为其对时间的微分可以忽略，所以用传热单元数法来求解比用传统方法更稳定.对于偏微分方程，用有限差分法对其进行离散求解，对于代数方程，则用牛顿－拉夫森法迭代求解.数值求解式：式中：通过对整个催化转化器的积分可以得到如下NTU形式的气体浓度关系：式中：根据质量守恒定律，质量平衡方程可写成.由边界条件式（7）与式（8），采用时间向前差分、空间中心差分的有限差分法可以从方程（4）求得催化转化器的温度Ts.方程（3）和（4）中的参数S cat定义为单位反应体积的催化剂表面积，本研究使用的整体催化剂的S cat值在整个催化器上都是一致的，其值为30 000 m2/m3.下面将要讨论的是起燃过程中贵金属负载量和负载模式对转化率的影响.为讨论方便，把均匀负载模式称为“基准模式”，把S cat＝30 000 m2/m3的负载量称为“基准负载”.研究贵金属负载量和负载模式参数的简单方法：如果负载量增加的同时，贵金属颗粒尺寸保持不变，催化剂表面积随负载量呈比例增加.因此如果贵金属负载量是基准负载的两倍，S cat也增加两倍，成为60 000 m2/m3，依此类推，给定的S cat分布就表示负载模式.2．1 均匀分布的负载量对转化率的影响图1中的3条曲线分别表示的是1倍、2倍和4倍于基准负载的负载模式对CO转化率的影响.从图中可以看出，随着负载量的增加，催化剂的起燃时间缩短，CO 的转化率增大.但是4倍于基准负载的负载模式与2倍于基准负载的负载模式相比，没有明显的变化.以上现象可以解释如下：因为反应速率常数的频率因子与催化剂负载量成比例，所以表面反应速率随着负载量的增加而增加.催化转化器中的总反应速率常数可用式（14）表示：在低温时，化学反应过程比传质过程要慢得多，所以总速率由化学反应速率决定，在高温时，化学反应速率比传质速率高，总反应速率就变成由传质速率控制.从式（14）可以看出，当负载量超过一定量而使表面反应速率足够高时，结果总反应速率受传质速率的限制，此时增加负载量对提高转化性能将不再起作用.2．2 沿轴向变化的负载量对转化率的影响保持贵金属总量不变，而使其负载量沿轴向呈线性变化，由此所得的3种转化率曲线如图2所示.从图2可以看出，以均匀负载为基准，负载量沿轴向递增时，其起燃时间与基准负载基本一致，而转化率较低，负载量沿轴向递减时，起燃时间大大缩短，转化率也比基准负载有所增加.这是因为从催化转化器内浓度分布曲线可以看出，高温气流进入反应器，首先加热反应器前端（距入口25 mm到30 mm），而使其先达到起燃温度，使大部分的反应物都在此处发生反应并放出大量的热，催化器前端也因此率先出现温度峰值.如图3所示，温度值也比同一时间的均匀负载量的催化器温度要高.因此可以通过在转化器前端加载高量贵金属来达到缩短起燃时间和改善转化率的目的.1）对于均匀负载模式，把负载量增加到基准负载的2倍时，转化率提高，当把负载量增加到4倍时，转化率增加不明显.这是由高负载量下，总反应速率由传质速率控制所造成的.2）在催化剂总量不变的情况下，使催化剂沿轴向呈线性分布，当负载量轴向增加，起燃时间与基准负载基本一致，但转化率下降；当负载量轴向递减时，起燃时间大大缩短，转化率提高.所以在对催化器优化设计时，可以选择负载量线性递减的模式，这样既可以缩短起燃时间，提高转化率，又可以节省催化剂的用量.【相关文献】［1］ Chan SH，Hoang DL．Chemical reactions in the exhaust system of a cold－start engine［J］．Chemical Engineering and Technology，2000，8（23）：727－730．［2］ Zhu J，Chan SH．An approach for rapid catalyst light off by high values of ignition retard［J］．Journal of the Energy Institute，1996，69：167－173．［3］ Chan SH，Hoang D L．Heat transfer and chemical reactions in exhaust system of a cold－start engine［J］．International Journal of Heat and Mass Transfer，1999，42（22）：4165－4183．［4］旻苏．机动车排放标准体系分析［J］．世界标准化与质量管理，2006（6）：29－31．［5］ Oh S H，Cavendish J C．Transients of monolithic catalytic converters：responsesto step changes temperature as related to controlling automobile emissions ［J］．Industrial and Engineering Chemical Production Research and Development，1982，21（1）：29－37．［6］ Keith JM，Chang N C，Leighton D T，et al．Designing a fast－igniting catalyticconverters［J］.AIChE journal，2001，47（3）：650－657．［7］ Ramanathan K，Balakltaiah V，West DH．Light－off criterion and transient analysis of catalytic monoliths［J］．Chemical Engineering Science，2003，58（8）：1381－1391．［8］ Hayes RE，Kolaczkowski ST．Mass and Heat Transfer Effects in Catalytic Monolith Reactors［J］．Chemical Engineering Science，1994，49（21）：3587－3599．［9］ Santos H，Costa M．Analysis of the mass transfer controlled regime in automotive catalytic converters［J］．International Journal of Heat and Mass Transfer，2008，51（1－2）：41－51．［10］ Kirchner T，Eigenberger G．On the dynamic behavior of automotive catalysts ［J］．Catalytic Today，1997，38（1）：3－20．［11］ Ma L P，Bart H J，Ning P，et al．Kinetic study of three－way catalyst of automotive exhaust gas：Modeling and application［J］．Chemical Engineering Journal，2009，155（1－2）：241－247．［12］ Baba N，Oshawa K，Sugiura S．Numerical approach for improving the conversion characteristics of exhaust catalysts under warming－up condition［J］．SAE paper 962076，1996．［13］ Bird R B，Stewart W E，Lightfoot E N．Transport Phenomena［M］．New York：Wiley，1960：217．。
