DreamBlooms/HansPub-15000-to-45000 huggingface dataset:
system: [SYSTEM]

user: 3. 结果与讨论 3.1. 水冷铜渣的综合特性分析 采用BT-1000型粉体综合特性测试仪对不同粒度铜渣的自然堆放料面与水平面之间的夹角(休止角)、振动后料面与水平面之间的夹角(崩溃角)、堆放在平板上振动前后料面和水平面之间夹角的平均值(平板角)、自然落下时飞散至空气中的粉尘百分比(分散度)等基本特性进行了检测。检测结果见表4，铜渣的颗粒表面光洁，流动性好，自然堆积的休止角和震动后的崩溃角都较小，粒度越大，堆角越小。铜渣的密度较大，粉尘含量少，不同粒度铜渣的分散度都在2.5%以下。 Table 4 Comprehensive properties of copper slag powders with different particle size 性能 <0.425 mm 0.425~2 mm 2~4 mm >4 mm 休止角θ r /˚ 21.3 20.1 19.5 19.4 崩溃角θ f /˚ 19.6 18.4 17.2 16.8 平板角/˚ 42.3 32.1 31.8 31.5 分散度D s /% 2.5 2.3 2.2 1.9 表4. 不同粒度铜渣的粉体综合性能数据表 3.2. 不同粒度铜渣的金相分析 为了研究不同粒度铜渣在常温条件下的金相组织，分别对4个粒级的铜渣进行金相试样镶嵌、磨制、抛光并用金相显微镜观察，结果见图1。 由图1(a)可知，当铜渣粒度在0.425 mm以下时，铁主要以铁橄榄石和一些氧化物的形式存在，如钙铁氧化物，而铁橄榄石分布在玻璃相中，其微观形貌呈树突状、板状或不规则形状晶体，少量的磁铁矿会被包裹在铁橄榄石相中，因而无法通过金相显微镜观察到。铜渣中的铜主要以铜锍的形式存在，黄铜矿、赤铜矿含量较少，主要分布在玻璃相中，斑铜矿在炉渣中主要嵌布在脉石矿物内，结晶良好，多于黄铜矿伴生，粒度较大。 由图1(b)可知，当铜渣粒度在0.425~2 mm之间时，渣中铁主要以铁橄榄石和铜铁氧化物的形式存在，铁橄榄石分布在玻璃相和脉石中，且含量较多，粒度较大。铜渣中的铜主要以铜锍的形式存在，其它含铜矿物较少，其中蓝铜矿结晶良好，粒度较粗，存在于玻璃相内。 由图1(c)可知，当铜渣粒度在2~4 mm时，渣中铁主要分布在铁橄榄石相和赤铁矿相中，铁橄榄石相分布在玻璃相中，含量较多，粒度较大。铜渣中的黄铜矿含量不高，与斑铜矿伴生，粒度较粗。 由图1(d)可知，含铁矿物细小，弥散分布在玻璃相中。铜渣中含铜矿物较多，主要以斑铜矿、黄铜矿、铜锍形式存在，分布在玻璃相和脉石中，单独分布或伴生存在，黄铜矿和斑铜矿粒度较粗。 图1. 不同粒度铜渣的金相显微组织：(a) ≤0.425 mm；(b) 0.425~2 mm；(c) 2~4 mm；(d) >4 mm 3.3. 不同粒度铜渣的物相分析 为了研究不同粒度铜渣在常温条件下的物相组成，选择其中3个粒度铜渣进行XRD分析，结果见图2。 由图2(a)可知，当铜渣粒度小于0.425 mm时，结晶不充分，未出现明显的结晶产物的衍射峰。衍射谱线中15˚~35˚之间呈现出典型的馒头峰，有非晶态的物质存在。通过分析可认为铜渣物相组成如下：结晶矿物主要有铁橄榄石、硅酸钙、钙铁氧化物及硫酸锌、SiO 2 等，存在少量的含铜矿物、磷酸盐矿物、MgO、TiO 2 等；非晶部分是无定形硅酸盐形成的玻璃体。 由图2(b)可知，当铜渣粒度在0.425~2 mm时，衍射谱线中15˚~40˚呈典型的馒头峰，有非晶态物质存在。通过分析可认为，铜渣物相组成中大多是无定形硅酸盐，有少量铁橄榄石和蓝铜矿存在。 由图2(c)可知，当铜渣粒度在4 mm以上时，结晶良好，出现明显的结晶产物的衍射峰。通过衍射图谱分析可知，铜渣主要由铁橄榄石和赤铁矿组成，少量MgO、CaO和Al 2 O 3 等主要以固溶体形式存在于铁橄榄石相中，无法检测到同时还存在少量的黄铜矿。 3.4. DTA-TG分析 为了研究不同粒度铜渣在升温过程中所产生的物理变化及化学反应，对其中三个粒度的铜渣进行热分析，分析结果见图3。 图2. 不同粒度铜渣XRD分析：(a) 0.425 mm；(b) 0.425~2 mm；(c) 2~4 mm 由图3(A)可知，在130℃以下，由于吸附水和层间水的失去，铜渣一直处于失重状态，失重约0.42%；200℃~600℃失重原因主要是硫酸盐矿物、碳酸盐矿物分解产生气体造成的；在590℃~615℃之间，对应的质量有所增加，增重约0.37%，此时的反应主要为金属铜等有机质的氧化反应；在652.6℃~679.9℃之间有一个大的放热峰，但重量并没有变化，可以推断此时发生的反应主要是γ-Fe 2 O 3 →α-Fe 2 O 3 的晶型转变过程；在997℃左右，DTA曲线上有一个吸热峰，对应的TG曲线上质量有所增加，增重约1.08%，此时的反应主要为橄榄石矿物的分解过程。 图3. 不同粒度铜渣的DTA-TG曲线：(A) <0.425 mm；(B) 0.425~2 mm；(C) 2~4 mm 由图3(B)可知，在400℃以下，铜渣一直处于失重状态，失重约2.03%，此阶段主要是水分的挥发以及一些硫酸盐矿物、碳酸盐矿物分解产生气体造成失重；在650℃左右，对应的TG曲线质量有所增加，此时发生的反应主要是Fe 3 O 4 被氧化为Fe 2 O 3 ；在1130℃左右，质量有所增加，此时发生的反应主要是橄榄石矿物的分解过程。 由图3(C)可知，在500℃以下，铜渣一直处于失重状态，失重约4.4%，此阶段主要是水分的挥发以及伴随一些盐类矿物的分解；在650℃左右，有一个明显的吸热峰，对应的TG曲线质量有所增加，增重约0.08%，此时的反应主要为金属铜等有机质的氧化反应；在700℃以后，重量一直在上升，发生的反应主要是橄榄石矿物的氧化过程。
system: [SYSTEM]

user: 4. 结论 通过采用金相分析、X射线衍射分析和DTA-TG分析对不同粒度铜渣原矿进行基本特性分析，得出如下结论： 1) 铜渣原矿中含有较多Fe 2 O 3 、TiO 2 ，有害杂质S含量较高，流动性较好。粒度大多分布在0.425~4 mm之间。 2) 经分析确定铜渣原料存在铁橄榄石(Fe 2 SiO 4 )、赤铁矿(Fe 2 O 3 )、磁铁矿(Fe 3 O 4 )和无定型结构的脉石玻璃体，铁橄榄石是铜渣冷却过程中Fe析晶时与SiO 2 结合而成，是斜方晶系；随着煅烧温度的升高，首先是铜渣中的铁橄榄石氧化转变为赤铁矿和非晶态硅石(2FeO·SiO 2 + 0.5O 2 →α-Fe 2 O 3 + SiO 2 )，其次是磁铁矿的晶型发生转变(Fe 3 O 4 →γ-Fe 2 O 3 →α-Fe 2 O 3 )。 3) 通过DTA-TG分析可知，在650℃左右，发生的反应主要是铜渣中铜的氧化反应，在1100℃左右，发生的反应主要是铜渣中橄榄石矿物的分解反应。
system: [SYSTEM]

user: 摘 要 奇异摄动问题一直需要高精度的数值方法，本文提出了重心拉格朗日插值配点法来求解一类奇异摄动抛物型反应扩散问题，与其他方法相比，数值实验表明该方法是解决这类问题的一种高精度方法。 关键词 :奇异摄动，重心拉格朗日插值配点法，抛物型反应扩散问题 Copyright © 2019 by author(s) and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 我们考虑下列的奇异摄动抛物型反应扩散问题： { ∂ u ∂ t − ε 1 ∂ 2 u ∂ x 2 + α 1 u + α 2 v + α 3 r = f 1 ( x ， t ) ， ∂ v ∂ t − ε 2 ∂ 2 v ∂ x 2 + β 1 u + β 2 v + β 3 r = f 2 ( x ， t ) ， ( x ， t ) ∈ [ 0 ， 1 ] × [ 0 ， 1 ] ， ∂ r ∂ t − ε 3 ∂ 2 r ∂ x 2 + γ 1 u + γ 2 v + γ 3 r = f 3 ( x ， t ) ， (1) 具有初始边界条件 u ( x ， 0 ) = f 0 ( x ) ， v ( x ， 0 ) = g 0 ( x ) ， r ( x ， 0 ) = s 0 ( x ) ， x ∈ [ 0 ， 1 ] u ( 0 ， t ) = f 1 ( t ) ， u ( 1 ， t ) = f 2 ( t ) ， v ( 0 ， t ) = g 1 ( t ) ， v ( 1 ， t ) = g 2 ( t ) ， t ∈ [ 0 ， 1 ] r ( 0 ， t ) = s 1 ( t ) ， r ( 1 ， t ) = s 2 ( t ) ， (2) 其中 ε 是任意小参数， α i ， β i ， γ i 为给定常数， f 1 ( x ， t ) ， f 2 ( x ， t ) ， f 3 ( x ， t ) ， f 1 ( t ) ， f 2 ( t ) ， g 1 ( t ) ， g 2 ( t ) ， s 1 ( t ) ， s 2 ( t ) 是已知函数， u ( x ， t ) ， v ( x ， t ) ， r ( x ， t ) 分别是未知函数。 奇异摄动抛物型反应扩散问题 [ 1 ] [ 2 ] [ 3 ] [ 4 ] [ 5 ] 在热科学和力学中有广泛的应用，有一些有价值的努力专注于找到解决问题的分析和数值方法 [ 6 ] - [ 11 ] 。在 [ 12 ] [ 13 ] [ 14 ] [ 15 ] [ 16 ] 中，J.P. Berrut引入了重心拉格朗日插值，并研究了它的数值稳定性和收敛性。重心拉格朗日插值在切比雪夫节点处无条件稳定。在 [ 17 ] [ 18 ] 中，Li和Wang给出了BLICM的算法和程序，用于普通的微分方程，积分方程和差分积分方程以及偏微分方程等。在 [ 19 ] - [ 24 ] 中，Wang运用了BLICM解决耦合粘性Burgers方程，奇异扰动问题和混沌系统。在 [ 25 ] 中，BLICM被提出用于求解第二类线性和非线性高维Fredholm积分方程。参考文献 [ 17 ] - [ 26 ] 的结果表明BLICM具有高精度，良好的稳定性和收敛性。在本文中，BLICM被提出用于求解一类奇异摄动抛物型反应扩散问题。
system: [SYSTEM]

user: 4. 结论 在本文中，我们首次使用重心插值配点法来解决奇异摄动问题。该方法的一些计算结果与其他方法的比较表明，该方法具有较高的精度和收敛性。从这篇文章中，我们可以发现我们的方法可以应用于解决这样的人口统计模型。因此，我们可以将此方法扩展到更广泛的领域。本文的所有数值计算结果都是通过数学软件MatlabR2007b给出。
system: [SYSTEM]

user: 摘 要 血清胱抑素C (serum cystatin C， SCysC)又名γ 2 痕迹碱性蛋白或后γ球蛋白，是一种低分子量的非糖基碱性分泌性蛋白，属半胱氨酸蛋白酶抑制剂家族。在病理状态下，尤其是在冠心病后慢性心衰(Chronic heart failure， CHF)后，SCysC表达水平会升高，提示SCysC在其病理生理过程中有重要作用。近年研究证实，冠心病后CHD患者SCysC水平高，且SCysC水平与冠心病后CHF的严重程度密切相关。本文主要对SCysC的分子结构与功能及其与冠心病后CHF的关系进行综述，为临床诊断和治疗CHF提供参考意义。 关键词 :冠心病后慢性心衰，血清胱抑素C，综述 Copyright © 2019 by authors and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 慢性心力衰竭(CHF)是大多数心血管疾病的最终归宿及最主要的死亡原因。近年来其发病率持续增长，有临床症状患者的5年存活率与恶性肿瘤相仿 [ 1 ] 。早期发现并处理CHF相当重要。《2018年中国心力衰竭诊断和治疗指南》流行病学显示 [ 2 ] CHF病因以冠心病居首。各年龄段心衰病死率均高于同期其他心血管病。有研究表明 [ 3 ] 血清CysC可抑制内源性半胱氨酸蛋白酶活性，参与细胞外基质产生和降解间的动态平衡，降低心肌顺应性，参与心室重构，导致心功能恶化。国外研究显示 [ 4 ] [ 5 ] [ 6 ] [ 7 ] [ 8 ] 血清CysC可通过对血管内皮细胞功能、血管平滑肌细胞、凝血及脂质过氧化的影响对心肌及内膜产生毒性作用，并参与炎症反应，抑制酶及激素前体的活性，与左室肥厚及舒张功能有关。近年国外研究表明SCysC参与心室重构，而国内少见报道 [ 9 ] 。本文主要对SCysC的分子结构与功能及其与冠心病后CHF的关系进行综述，为临床诊治CHF及其预后提供参考。
system: [SYSTEM]

user: 2. SCysC的分子结构与功能 SCysC由122个氨基酸组成，相对分子量为13 KD，在机体内广泛分布，生理作用为使细胞膜免受内源性或者外源性的蛋白酶水解，在动脉壁溶解平衡与抗蛋白溶解活性平衡中发挥重要作用 [ 10 ] [ 11 ] 。许多研究证实SCysC为评价肾功能的理想内源性标志物，与肌酐相比更敏感准确 [ 12 ] [ 13 ] 。近年研究发现，肾功正常的患者研究也存在SCysC水平升高 [ 14 ] 。国外相关研究表明其参与许多心管系统疾病的病理生理过程 [ 15 ] 。
system: [SYSTEM]

user: 3. SCysC与冠心病后慢性心衰的关系 冠心病后CHF的实质是心肌细胞损伤、心肌细胞凋亡及细胞外基质的改变致使心室重构，最后导致CHF [ 16 ] 。以往人们关注的焦点总是心肌细胞的改变，现人们逐渐认识到了心肌细胞外基质改变在CHF发生发展中的作用。细胞外基质重构的发生是心肌胶原合成与降解失衡的结果 [ 17 ] [ 18 ] 。有研究表明 [ 15 ] SCys-C可抑制内源性半胱氨酸蛋白酶活性，参与细胞外基质产生和降解间的动态平衡，降低心肌顺应性，参与心室重构，导致心功能恶化。国外研究显示 [ 4 ] [ 5 ] [ 6 ] [ 7 ] [ 8 ] SCys-C可通过对血管内皮细胞功能、血管平滑肌细胞、凝血及脂质过氧化的影响对心肌及内膜产生毒性作用，并参与炎症反应，抑制酶及激素前体的活性，与左室肥厚及舒张功能有关。相关研究表明SCysC水平与冠心病后心衰患者心功能NYHA分级呈正相关 [ 19 ] 。此外，有学者研究表明高水平SCysC可增加冠心病的全因死亡率、心血管事件及心衰的发生率 [ 20 ] 。SCysC可作为冠心病后慢性心衰的患者近期心血管事件及死亡的独立预测因子与近期的预后有密切的关系 [ 21 ] 。
system: [SYSTEM]

user: 4. 小结与展望 冠心病后慢性心力衰竭是心血管疾病的终末期表现及最主要的死因。近年来，CHF患病率逐年升高，严重威胁患者的生活质量及健康寿命，SCysC可作为临床诊断及判断冠心病后CHF严重程度的实验室指标 [ 20 ] 。早期根据SCysC水平早期进行危险分层，尽早给予合理的治疗措施可改善预后。但SCysC与冠心病后CHF相关性的更确切机制尚待进一步深入研究。
system: [SYSTEM]

user: 摘 要 目的：本文设计了一种自动控制的呼吸机呼出废气处理装置，该处理系统能在不影响患者吸气呼气的情况下，将病患呼出的气体排入医院中心吸引系统，从源头上减少感染。方法：利用STC12C5201AD型微处理器、呼气缓冲盒、压力传感器、真空泵和紧急呼气阀等组件进行系统电路和气路设计。结果：经测试，该装置可以在不影响病人接受机械通气治疗的前提下，将其呼出废气排入医院中心吸引系统。结论：该呼吸机呼出废气处理装置能有效防止病人呼出废气污染病房，降低交叉感染发生率 [ 1 ] ，系统具有操作便捷、实用性强、安全性高等特点。 关键词 :微处理器，机械通气治疗，废气处理，交叉感染 Copyright © 2019 by author(s) and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 呼吸机是挽救患者生命重要的生命支持与急救设备 [ 2 ] ，具有代替、控制或改变人的正常生理呼吸，增加肺通气量，改善呼吸功能，减轻呼吸功消耗，节约患者心脏储备能力等功能 [ 3 ] ，它在各层级医院的重症医学科、急诊科、呼吸内科等科室的医疗救治中起着不可代替的作用。与此同时，呼吸机也是风险为16分(依据美国佛蒙特大学的医疗设备风险评估系统，呼吸机风险值为16分，满分为17分，分数越高，风险越大)的高风险生命支持类设备 [ 4 ] ，需要医院医学工程部门投入大量的精力对其进行质量安全管理。一方面，大多数医院的医学工程部门目前已经意识到对呼吸机各项功能、参数、指标进行严格检测以保障其稳定运行的重要性 [ 5 ] ；但另一方面，却很少有人意识到，机械通气治疗的过程中，病人的呼出气体未经处理直接排放到急救病房而产生的各种危害。 因大多数进行机械通气治疗的病人均存在肺部感染，而大部分机械通气治疗患者需要长期使用呼吸机进行辅助通气治疗，呼吸机与其管道的温湿度都比较适合细菌生长 [ 6 ] ，长时间的通气会导致大量细菌滋生，引起呼吸机相关性肺炎(VAP)等并发症 [ 7 ] ；加之重症监护病房空气密闭，行机械通气的感染患者飞沫和接触传播风险明显 [ 8 ] ，如果病人呼出的气体不经处理直接排放到病房，长此以往，不仅会加大病患之间交叉感染的发生率 [ 9 ] ，危害病人生命安全，亦会增加医护人员和探视家属等人群的潜在感染风险。基于上述危害，本文设计了一种自动控制的呼吸机呼出废气处理装置，该处理系统能在不影响患者吸气呼气的情况下，将病患呼出的气体排入医院中心吸引系统，从源头上减少感染。
system: [SYSTEM]

user: 2. 呼出废气处理系统设计 2.1. 工作原理 本文所述呼吸机废气处理装置气路图如图1，箭头所示方向为气体流动方向，病人所呼出的气体首先进入呼气缓冲盒，然后由压力传感器对呼气缓冲盒内压力进行实时监测，若采集到呼气缓冲盒内的压力大于大气压，则真空泵开始工作，将缓冲盒内病人气体抽出并排入医院中心吸引系统，与此同时，始终保持呼气缓冲盒内压力与当地大气压力一致，保证病人呼气过程不受外力干扰。 2.2. 气路设计的要点 2.2.1. 呼气缓冲盒的使用 呼气缓冲盒在顾名思义在系统中起到呼吸缓冲的作用，如果不设计呼气缓冲盒，而将病人呼出气体端直接连接到真空泵输入端，则由于呼气管道腔体积的限制，呼出气路中阻力增大，病人在呼气过程中的做功增加，导致病人呼气时间延长，从而产生憋气的不适感，严重时，甚至导致二氧化碳潴留等严重后果。有了呼气缓冲盒，病人呼出的气体先经过缓冲盒，就不会在短时间加大呼出气路阻力，从而不会影响病人的呼吸治疗过程。 图1. 系统气路图 2.2.2. 真空泵的使用 本设计中真空泵一方面将呼气缓冲盒的气体及时排入医院中心吸引系统，另一方面将呼气缓冲盒与中心吸引系统进行隔离。若本系统中不设计该真空泵，而将呼气缓冲盒输出端与中心吸引系统输入端直接相连，虽然由于负压的作用也能及时将气体吸入中心吸引系统，但是在病人呼气过程中，亦由于负压的作用，病人呼气时间将会大大变短，使得病人呼气变得急促，干扰正常呼吸通气治疗过程。有了真空泵的隔离，病人的呼气过程便不会受中心吸引系统负压的影响。 2.2.3. 紧急呼气阀的使用 出于高级别安全考虑，医疗产品通常需要设计紧急处理装置。本设计中紧急呼气阀就是这样一个紧急装置，在未供电或掉电情况下该电磁阀在处于开启状态(即为常开)，系统上电后该电磁阀将会关闭，在紧急状况下，如停电或病人呼气回路压力异常(过高或过低)的情况下紧急开启，将呼气缓冲盒与大气直接连通，开放病人呼气通道，保证病人安全。 2.3. 电路设计的要点 2.3.1. 最小系统设计 本系统采用宏晶公司的STC12C5201AD型芯片作为本系统的微处理器，STC12C5201AD型处理器是STC生产的单时钟，机器周期为1T的单片机，是具有高速、低功耗，超强抗干扰的新一代8051单片机，指令代码完全兼容8051，但是速度却快8~12倍。内部集成2路PWM，8路高速8位A/D转换器，适用于针对电机控制等强干扰场合，能安全满足本系统设计需求 [ 10 ] 。图2中Sound Alarm、Led Alarm、Ex Valve、Rezero Valve、Vac pump和Pressure分别为声音报警、发光报警、紧急呼出阀控制、校零电磁阀、真空泵控制、压力采集的英文缩写网络标号。 图2. 最小系统电路图 2.3.2. 压力采集电路设计 本系统中压力采集电路采用MPXM2053GS芯片，该芯片为美国飞思卡尔半导体公司研制的硅材料压阻压力传感器，此传感器能提供一个高精度的线性差分输出电压，该输出电压与施加于其的压力值大小成正比，该传感器带有应变片和片上集成的薄膜电阻网络，可实现精确的偏移校准和温度补偿。具体采集电路如图3，压力传感器将压力信号转变为电信号后经INA333运算放大器进行放大，然后输出到微控制器进行AD采样，从而计算出精确的压力值。 图3. 压力采集电路图 2.3.3. 泵阀控制电路设计 因微处理器I/O口不足以直接驱动功耗较大的泵阀组件，故需要应用大功率驱动芯片来配合单片机进行控制。如图4本系统中采用摩托罗拉半导体公司的ULN2803芯片，该芯片中的八个NPN达林顿连接晶体管非常适用于低逻辑电平数字电路或其他类似电路的更高电流、电压要求之间的接口，适用于各种计算机，工业和消费类应用，非常适用于本系统中对泵阀组件的驱动。 图4. 泵阀驱动电路图 2.4. 系统运行流程图 呼吸机废气处理系统工作流程如图5，因各地区呼吸机使用地理环境不一样，系统上电后首先执行大气压力校准，获取当地大气压力值，然后通过压力传感器对呼气缓冲盒内压力进行监测，若缓冲盒内压力大于大气压则真空泵开始工作，将缓冲盒内气体抽出，直至其压力接近大气压则停止。若监测到缓冲盒内压力过高，则进入过压保护程序，打开紧急呼气阀，开放呼气通道，并进行声光报警，提示工作人员进行故障处理，如此循环往复。 图5. 系统运行流程图
system: [SYSTEM]

user: 3. 结语 目前，无论是临床科室还是医学工程部门都对呼吸机治疗过程中病患产生的气体污染关注较少，随着空气污染越来越严重，我国患有呼吸系统疾病的人越来越多，随之而来的需要接受机械通气治疗的病患也愈多。而呼吸机相关性肺炎已经成为机械通气患者常见的严重并发症 [ 11 ] ，我国VAP的发病率为4.7%~55.8%，病死率为19.4%~51.6% [ 12 ] ，可见医院需要进一步规范呼吸治疗过程，最大限度避免交叉感染 [ 13 ] ，医学工程技术人员应积极参与到此过程中去，利用自己的专业知识与临床科室一起寻找解决问题的方法，提高医院医疗设备的安全性。 因本文所设计呼吸机呼出废气处理系统与呼吸机主体在控制中互不影响，为两套独自运行的系统，故本系统适用于大多数呼吸机，适用性强，并不局限于某个品牌或某个型号的呼吸机。该呼吸机呼出废气处理装置在使用过程中无需对其进行任何的参数调节，具有操作便捷、实用性强、安全性能高的特点。
system: [SYSTEM]

user: 摘 要 本文在时变波动率Black-Scholes模型下，研究了到期区间理财产品的定价问题。首先，利用偏微分方程方法得到了金融理财产品的定价公式。接着，利用二次变差方法提取积分波动率，再选取四种理财产品进行实证分析，并对其选择给出建议。最后，为投资者选择金融理财产品提供了一系列的参考建议。 关键词 :金融理财产品，偏微分方程，定价公式，二次变差 Copyright © 2019 by author(s) and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 随着时代的发展，我国个人财富以及个人理财有了较大增长。为了抵补资本在时间上的损失，人们常常将流动资产用于投资。同时，银行为了吸纳闲散资金，推出了各种理财产品，到期区间理财产品就是其中一种。 理财产品的核心问题是定价，一个理财产品是否成功，关键在于定价是否合理，是否能让发行单位和投资人实现双赢。金融理财产品的定价源自于欧式期权定价。目前，相关期权定价方法的理论研究日趋成熟。文献 [ 1 ] 以零息票为计价单位，研究了利率服从HJM模型下的欧式期权定价问题。为了刻画欧式期权价格估值的不确定性和投资者的犹豫程度，文献 [ 2 ] 用三角直觉模糊数表示风险资产的变化因子，构建了三角直觉模糊数的二叉树模型，进而研究了欧式期权定价问题。文献 [ 3 ] [ 4 ] 采用了摄动方法研究了非线性Black-Scholes模型下的欧式障碍期权定价问题。在一定的假设条件下，用Green函数分析了近似结论的误差估计。之后文献 [ 5 ] 采用单参数方法分别研究了修正的欧式期权定价问题，并且采用了Feymann-Kac公式分析了相应结论的误差问题。 关于期权定价方法的研究国内外有很多，但具体到银行理财产品定价问题的文献却不多见。考虑理财产品的价值是投资人选择理财产品的重要指标之一，同时，理财产品的定价对理财产品的投资活动具有指导意义。基于此，本文研究了挂钩于沪深300指数的到期区间理财产品定价问题。
system: [SYSTEM]

user: 2. 金融市场数学模型 表1收集了四款挂钩于沪深300指数的到期区间理财产品。在进行理财产品价值分析之前，先对沪深300指数适合的随机模型进行识别和参数估计。本节选取了2018年5月16日至8月14日的沪深300指数数据，时间间隔为1天，共计63个数据(见表2)。 Table 1 Agricultural bank of China financial services product 序号 产品名称 发售期 理财期限/天 起购金额/ 万元 预期年化收益率 1 金钥匙·如意组合2016年第283期沪深300指数到期宽幅波动人民币金融理财产品 16.08.16~16.08.22 62 5 3%~3.95% 2 金钥匙·如意组合2016年第283期沪深300指数到期窄幅波动人民币金融理财产品 16.08.16~16.08.22 62 5 3%~3.95% 3 金钥匙·如意组合2016年第284期沪深300指数到期窄幅波动人民币金融理财产品 16.08.16~16.08.22 90 5 3%~4.2% 4 金钥匙·如意组合2016年第284期沪深300指数到期宽幅波动人民币金融理财产品 16.08.16~16.08.22 90 5 3%~4.2% 表1. 中国农业银行金融理财产品 Table 2 The csi 300 index dat 周次 日期 周一 周二 周三 周四 周五 1 5.16~5.22 3064.53 3095.52 3071.53 3060.34 3047.78 2 5.23~5.29 3078.52 3083.24 3079.75 3056.60 3059.73 3 5.30~6.05 3056.31 3068.60 3172.96 3158.03 3172.95 4 6.06~6.12 3192.78 3182.44 3171.81 - - 5 6.13~6.19 3134.05 3058.44 3043.96 3104.36 3096.09 6 6.20~6.26 3114.91 3124.90 3100.45 3129.72 3110.65 7 6.27~7.03 3056.13 3107.40 3142.48 3152.83 3156.93 8 7.04~7.10 3136.39 3199.16 3197.63 3206.55 3199.75 9 7.11~7.17 3199.04 3201.91 3274.02 3277.48 3278.84 10 7.18~7.24 3269.71 3260.44 3246.86 3238.34 3251.24 11 7.25~7.31 3220.17 3228.23 3270.08 3204.47 3217.19 12 8.01~8.07 3196.43 3173.76 3179.68 3190.55 3201.35 13 8.08~8.14 3199.57 3232.60 3255.18 3239.56 3231.81 表2. 沪深300指数数据 首先，将数据录入Eviews软件，取对数并一次差分后进行ADF检验，输出结果见表3，由结果可知t统计量为−8.333832，比三个临界值都要小，说明处理之后的数据不存在单位根，即通过单位根检验，故我们认定数据平稳。 Table 3 The csi 300 index opening price: the logarithmic difference the stationarity of dat t-Stalistlic Prob. * Augmented Dickey-Fuller test statistic −8.333832 0.0000 Test critical values: 1% level −3.542097 5% level −2.910019 10% level −2.592645 表3. 沪深300指数开盘价格：对数差分数据的平稳性 其次，由数据的自相关性(见图1)，数据的自相关系数全部位于虚线内并且对应的P值均大于0.05，这表明数据之间不相关，符合白噪声序列的条件之一。 图1. 沪深300指数开盘价格：对数差分数据的相关性 最后，由描述统计数据(见图2)，均值为−0.000204，中位数为−0.000641，最大值为0.034319，最小值为−0.038160，标准差为0.014282，偏度系数为0.035154，表明数据无偏，峰度系数为3.407955，表明数据的分布比标准正态分布略尖。计算得到的JB统计量的值为0.435568，对应的P值为0.804299，不能拒绝对数数据的一次差分是正态的。由表4也不能拒绝沪深300指数对数一次差分数据之间是独立的，从而我们可以认定沪深300指数对数一次差分为正态白噪声序列。 图2. 沪深300指数开盘价格：对数数据一次差分的统计描述 Table 4 The csi 300 index opening price: the independence of the log data of a differenc Dimension BDS Statistic Std.Error z-Statistic Prob. 2 0.047739 0.010722 4.452420 0.0000 3 0.061689 0.017262 3.573607 0.0004 4 0.059362 0.020827 2.850317 0.0044 5 0.047317 0.021996 2.151191 0.0315 6 0.037657 0.021497 1.751747 0.0798 Raw epsilon 0.020928 Pairs within epsilon 2635.000 V-Statistic 0.708143 Triples within epsilon 123249.0 V-Statistic 0.542993 Dimension C(m，n) c(m，n) C(1，n-(m-1)) c(1，n-(m-1) c(1，n(m-1)) 2 949.0000 0.536158 1237.000 0.698870 0.488419 3 687.0000 0.401520 1194.000 0.697838 0.339831 4 480.0000 0.290381 1146.000 0.693285 0.231019 5 326.0000 0.204261 1102.000 0.690476 0.156944 6 219.0000 0.142208 1057.000 0.686364 0.104550 表4. 沪深300指数开盘价格：对数数据一次差分的独立性 基于以上分析，本文采用如下修正的几何Brown运动模型来描述进行对数一次差分处理后的沪深300指数 d S t = μ ( t ) S t d t + σ ( t ) S t d ω t (1) 其中 ω t 表示定义在完备概率空间 ( Ω ， F ， Ρ ) 上的布朗运动，波动率 σ ( t ) ，无风险收益率 r ( t ) 和股票期望收益率 μ ( t ) 都是时间的函数。 在接下来的章节中，将用公式(1)刻画沪深300指数，并用其研究挂钩于沪深300指数的到期区间理财产品定价问题。
system: [SYSTEM]

user: 3. 区间金融理财产品的定价公式 根据到期区间理财产品的收益说明(见表1)，其到期之后的收益可以表示为 ( 1 + r 1 T ) I { S ( T ) ∈ [ K 1 ， K 2 ] } + ( 1 + r 2 T ) I { S ( T ) ∉ [ K 1 ， K 2 ] } (2) 其中 S ( T ) 表示在到期日T时刻黄金的价格， r 1 表示起始日的利率， r 2 表示到期日的利率，[K 1 ， K 2 ]表示约定的执行区间。根据金融市场的自融资策略 [ 6 ] ，投入一个单位的本金、并且损益公式为公式(2)的金融理财产品的价格满足抛物偏微分方程。 { L C = 0 ，     ( t ， S ) ∈ ( 0 ， T ] × R + C ( T ， S ) = ( 1 + r 1 T ) I { S ∈ [ K 1 ， K 2 ] } + ( 1 + r 2 T ) I ， { S ∉ [ K 1 ， K 2 ] }     S ∈ R + (3) 其中 L C = ∂ t C + 1 2 σ 2 ( t ) S 2 ( t ) S 2 ∂ s s C + r ( t ) S ∂ s C − r ( t ) C 下面我们考察公式(3)的解。令 x = ln S + ∫ t T r ( z ) d z − s ，     s = 1 2 ∫ t T σ 2 ( s ) d s (4) 并定义 C ( t ， S ) = exp { − ∫ t T r ( z ) d z } U ( s ， x ) (5) 则 ∂ C ∂ t = exp { − ∫ t T r ( z ) d z } [ r ( t ) U − 1 2 σ 2 ( t ) ∂ U ∂ s + ( − r ( t ) + 1 2 σ 2 ( t ) ) ∂ U ∂ x ] ∂ C ∂ S = 1 S exp { − ∫ t T r ( z ) d z } ∂ U ∂ x ∂ 2 C ∂ S 2 = 1 S 2 exp { − ∫ t T r ( z ) d z } [ ∂ 2 U ∂ x 2 − ∂ U ∂ x ] 将上述结果代入公式(3)得 ∂ U ∂ s = ∂ 2 U ∂ x 2 ，     ( s ， x ) ∈ [ 0 ， T ] × R 由公式(4)和公式(5)可知，当 t = T ， U ( 0 ， x ) = C ( T ， S ) 时 U ( 0 ， x ) = ( 1 + r 1 T ) I { e x ∈ [ K 1 ， K 2 ] } + ( 1 + r 2 T ) I { e x ∉ [ K 1 ， K 2 ] } (6) 根据热传导方程的经典解理论，其唯一强解可表示为 U ( s ， x ) = 1 2 π s ∫ R u 0 ( z ) exp { − ( z − x ) 2 4 s } d z = I 1 + I 2 (7) 其中 I 1 = ( 1 + r 1 T ) 1 2 π s ∫ R I { e x ∈ [ K 1 ， K 2 ] } exp { − ( z − x ) 2 4 s } d z I 2 = ( 1 + r 2 T ) 1 2 π s ∫ R I { e x ∉ [ K 1 ， K 2 ] } exp { − ( z − x ) 2 4 s } d z 利用公式(7)可得本文的主要结论。 定理1 [ 7 ] 损益为公式(2)的金融理财产品在t时刻的价格为 C ( t ， S ) = ( 1 + r 1 T ) exp { − ∫ t T r ( z ) d z } Φ ( d 1 ) − ( 1 + r 1 T ) exp { − ∫ t T r ( z ) d z } Φ ( d 2 )     + ( 1 + r 2 T ) exp { − ∫ t T r ( z ) d z } Φ ( − d 1 ) + ( 1 + r 2 T ) exp { − ∫ t T r ( z ) d z } Φ ( d 2 ) (8) 其中 d 1 = ln S − ln K 1 + ∫ t T r ( z ) d z − 1 2 ∫ t T σ 2 ( s ) d s ∫ t T σ 2 ( s ) d s d 2 = ln S − ln K 2 + ∫ t T r ( z ) d z − 1 2 ∫ t T σ 2 ( s ) d s ∫ t T σ 2 ( s ) d s 证明：*首先计算 I 1 ，利用积分换元，我们有 I 1 = ( 1 + r 1 T ) 1 2 π s ∫ ln K 1 ln K 2 exp { − ( z − x ) 2 4 s } d z = ( 1 + r 1 T ) 1 2 π ∫ ln K 1 − x 2 π ln K 2 − x 2 π exp { − y 2 2 } d y 从而 I 1 = ( 1 + r 1 T ) Φ ( x − ln K 1 2 s ) − ( 1 + r 1 T ) Φ ( x − ln K 2 2 s ) (9) 同理 I 2 = ( 1 + r 2 T ) Φ ( ln K 1 − x 2 s ) + ( 1 + r 2 T ) Φ ( x − ln K 1 2 s ) (10) 将公式(9)和公式(10)代入公式(7)，可得结论成立。* 证毕。
system: [SYSTEM]

user: 摘 要 采用ANSYS Fluent软件建立了方孔蜂窝陶瓷蓄热体内非稳态流固耦合传热的单元流道模型；以低热值燃料蓄热式燃烧为背景，对蓄热体内烟气、空气与骨架间复杂传热进行了三维数值模拟；预测了各种蓄热体长度下空气预热温度、蓄热体效能和压力损失等热流体性能，并据此进行结构优化；在此基础上研究蓄热体变流量条件下的传热和流动特性。此外，借助温度云图对蓄热式传热机理进行了分析。 关键词 :蜂窝蓄热体，三维CFD模拟，结构优化，变流量特性，蓄热式燃烧 Copyright © 2019 by author(s) and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 蓄热式燃烧是一项节能环保技术，在各种工业窑炉中得到了广泛应用 [ 1 ] [ 2 ] [ 3 ] [ 4 ] [ 5 ] 。蓄热式换热器是这种先进燃烧装置的关键部件之一，工作时高温烟气与常温空气以逆流形式周期性依次流过蓄热骨架，传热过程比较复杂 [ 5 ] 。与陶瓷小球相比，蜂窝陶瓷作为蓄热骨架能在保证传热性能基础上大幅度减小流动阻力 [ 6 ] 。除了实验以外 [ 7 ] [ 8 ] [ 9 ] [ 10 ] ，国内外学者对蜂窝蓄热体内部流动传热性能进行了数值模拟研究，比如，陈维汉和罗海兵从圆孔蜂窝中抽取单元结构，建立了二维轴对称的非稳态传热数值模型 [ 11 ] 。Rafidi和Blasiak采用钢玉和堇青石制作复合方孔蜂窝蓄热体，并在假定流动充分发展进而具有理论抛物线形速度场的基础上建立了二维数值模型 [ 12 ] 。Doh等借助商业软件CFX进行了整机数值模拟研究 [ 13 ] ，他们采用不到3万网格模拟1/4的蓄热体，计算结果与实验值相差比较大。值得指出，上述蓄热体研究都是针对高温空气燃烧开展的。在另一方面，工业生产产生大量低热值燃料，借助蓄热式技术，用烟气余热加热助燃空气，可以使这些燃料的燃烧变得稳定和充分。与高温空气燃烧相比，低热值燃料燃烧的空气过剩系数通常较大，其排烟温度较低且受到凝酸的限制，需要进行专门研究。当前商业CFD软件广泛用于各种换热器的数值模拟研究 [ 14 ] [ 15 ] [ 16 ] [ 17 ] [ 18 ] ，但蓄热式换热器由于存在冷热流体周期性切换，使用CFD软件模拟的研究报道不多。为此，本文拟在对ANSYS Fluent二次开发的基础上，建立方孔蜂窝蓄热体复杂传热的数值模型，并针对低热值燃料蓄热式燃烧应用开展三维数值模拟，分析蓄热体内烟气、空气与固体骨架间的非稳态耦合传热，研究蓄热体长度等参数对传热和流动性能的影响，并在结构优化的基础上研究蓄热体变流量条件的工作特性。
system: [SYSTEM]

user: 2. 数值模型 2.1. 物理模型 本次数值研究采用高效堇青石方孔蜂窝蓄热体，方孔边长2.07 mm，壁厚0.43 mm。鉴于低热值燃料燃烧温度较低，烟气和空气进口温度分别为673 K和285 K [ 12 ] ，两者切换时间为30 s或60 s；空气进口速度为1.95 m/s，对应烟气进口速度5.53 m/s。考虑到增加蓄热体长度，排烟温度下降，可能低于酸露点温度，为此作者在100 mm到200 mm间选择5种长度进行研究。 2.2. 控制方程 蓄热体内部为典型的非稳态流固耦合传热。由于方孔尺寸小，内部流动为层流；由于烟气温度不高，忽略辐射换热；由于流体速度不大，忽略动能和粘性耗散的影响；蓄热体长度适中，烟气凝酸影响不计。基于上述假设，蓄热体内工作流体和固体骨架的能量守恒方程分别表示为： ρ f C p f ∂ T f ∂ t + ρ f C p f ( u ⋅ ∇ ) T f = ∇ ⋅ ( k f ∇ T f ) (1) ∂ ( ρ s C p s T s ) ∂ t = ∇ ⋅ ( k s ∇ T s ) (2) 其中t和T分别代表时间和温度；ρ，k和Cp分别表示密度、热导率和定压比热；下标f和s分别代表流体和固体；u指流体速度，由Navier-Stokes方程组确定，即 ρ f ∂ u ∂ t + ρ f ( u ⋅ ∇ ) u = − ∇ p + ∇ ⋅ [ μ f ( ∇ u + ∇ u T − 2 3 ∇ ⋅ u I ) ] (3) ∂ ρ f ∂ t + ∇ ⋅ ( ρ f u ) = 0 (4) 其中，p和μ f 分别代表流体压力和动力粘度。 2.3. 数值方案和后处理 考虑到新型蜂窝蓄热体方孔通道均匀布置，本文采用一个截取的对称单元结构来代表整个蓄热体，以减少数值计算工作量。对称单元结构由方孔通道和周围骨架组成，如图1所示。 图1. 方孔蜂窝蓄热体单元通道与计算网格 单元结构的热端为质量流量边界，冷端则采用压力出口边界；流体和固体界面为耦合传热面。控制微分方程组采用ANSYS Fluent数值求解。烟气和空气均为不可压缩理想流体，通过用户自定义函数来描述其热导率、粘度与温度的函数关系，以及两种流体的切换。迭代计算收敛判据是质量和动量方程相对残差小于1e-04，能量方程小于1e-07。通过4套网格验证网格无关性，经过综合评估计算负荷和预测精度，最终模型约14万网格，时间步长为1 s。
system: [SYSTEM]

user: 3. 热流体性能 蜂窝蓄热体工作足够长时间后，其内部传热和流动将达到周期性稳态，空气预热温度就等于此阶段空气预热周期离开蓄热体的平均温度，即 T ¯ a o = ∫ τ 2 τ d t ∫ 0 A ( ρ u z T ) | z = 0 d S / ∫ τ 2 τ d t ∫ 0 A ( ρ u z ) | z = 0 d S (5) 其中t = τ表示空气预热周期初始时刻；A代表热端(z = 0)的横截面积；u z 为流体轴向速度。 蓄热体效能(ε)等于其实际传热量(Q)与最大可能传热量(Q max )之比 [ 12 ] ，即 ε = Q Q max = h a o − h a i h a ( T = T g i ) − h a i (6) 其中h表示流体比焓；下标g和a分别代表烟气和空气，i和o则分别代表进口和出口。 能量回收率(ERR)表示蓄热体实际回收能量值与最多可能利用能量值之比，即 E R R = ρ a q V a ( h a o − h a i ) ρ g q V g [ h g i − h g ( T = T a i ) ] (7) 其中q V 代表体积流量。 蓄热体流阻性能由烟气周期平均压力损失和空气周期平均压力损失两部分组成，即 Δ p g ¯ 和 Δ p a ¯ 。
system: [SYSTEM]

user: 4. 模型验证 为了评价上述CFD模型的预测精度，本文对Rafidi和Blasiak的蓄热体 [ 12 ] 进行三维数值模拟，并将典型工况性能的计算值与文献中的实验值和模拟值进行比较，见表1所示。从表1可以看出，切换时间(τ)等于30 s或60 s时，本文模型预测的温度效率( η T = ( T ¯ a o − T a i ) / ( T g i − T a i ) )均与实验值非常接近，相对偏差比文献 [ 12 ] 模型预测的小。另外，作者监测了蓄热体工作时冷端和热端的流体温度，发现它们均以2 τ为周期不断变化，具体数值与文献 [ 12 ] 一致。通过上述比较，可以认为本文CFD模型计算可靠，且有足够精度。 Table 1 The comparison between the predicted values and the literature value 蓄热体工况 温度效率 备注 τ = 30 s τ = 60 s 本文模型 空气和烟气进口温度380，1373 K； 进口速度2.6和10 m/s 0.925 0.913 三维 文献 [ 12 ] 模型 0.882 0.842* τ = 50 s 实验 [ 12 ] 0.913 0.902 二维 表1. 本文模型预测值与文献值的比较
system: [SYSTEM]

user: 5. 计算结果及分析 5.1. 蓄热式换热分析 为便于分析烟气通过蓄热骨架向空气传热的物理机制，我们提取蓄热体周期性稳定后不同时刻的温度云图来进行讨论。图2(a)~图2(f)分别描绘了蓄热体(含工作流体)在t = 5 s、15 s、25 s、35 s、45 s、55 s时刻的温度分布情况，其中，空气–烟气切换时间等于30 s，t = 0 s则指烟气周期初始时刻。 从图2(a)，图2(b)和图2(c)可以看到，在烟气周期，随着高温烟气不断进入，蓄热体内流体和骨架温度不断升高，且前者温度总高于后者。这表明烟气将部分热量转移到骨架。另外，随着骨架温度上升，烟气放热速率减慢，出口温度上升。空气预热周期温度场见图2(d)，图2(e)和图2(f)，其变化规律与烟气周期相反，即随着空气不断进入蓄热体，流体和骨架温度都在下降，且前者温度总低于后者，这表明骨架将烟气周期蓄存的热量转移到空气中。另外，随着骨架温度降低，空气得热速率减慢，出口温度降低。 5.2. 蓄热体长度对热流体性能的影响 本文对100 mm和200 mm间5个不同长度(L)的蓄热体在30 s和60 s两个切换时间(τ)下进行了三维 图2. 蓄热体对称面不同时刻温度云图 数值模拟，其中空气进口温度285 K，进口速度(u a ) 1.95 m/s；烟气进口温度673 K，进口速度5.53 m/s。数值计算得到了蓄热体的传热和流动阻力性能，其中不同长度的蓄热体的空气预热温度( T ¯ a o )，以及在经过一定长度的蓄热体后的烟气出口温度( T ¯ g o )见图3(a)所示。从图3(a)可以看出， T ¯ a o 随着L的增加而增加， T ¯ g o 的变化趋势相反；延长τ， T ¯ a o 下降， T ¯ g o 则上升，且τ的影响随L增加而变弱。图3(b)给出了上述5种蓄热体的效能(ε)和能量回收率(ERR)。对比图3(a)和图3(b)，可以看出，ε和ERR的变化规律与 T ¯ a o 相同，即随L的增加而增加，随着τ的增加而减小。当τ = 30 s时，随着L从100 mm增加到200 mm，ε和ERR分别从~76.7%，61.8%增加到~86.0%和69.3%；蓄热体长200 mm时，将τ延长到60 s，ε和ERR分别减少了~1.3%和1.04%，但在L = 100 mm时，两参数减少值分别为~16.15%和13.01%。 上图曲线变化趋势可用蓄热体非稳态传热来解释。随着长度L增加，相同时间内烟气传递给骨架的热量增多，空气也因此从骨架获得更多热量。另一方面，随着流体(烟气或空气)不断流过蓄热体，流体与骨架的温差逐渐减少，这导致热流量随切换时间(τ)的延长而减少。因此，增加L，或者减小τ， T ¯ a o 会升高， T ¯ g o 则降低，与图3(a)相一致。鉴于流体温度越高，比焓越大，由式(6)和(7)可知，L越长，或τ越小，蓄热体ε和ERR越大，与图3(b)中曲线一致。 图3(c)给出了烟气和空气压力损失( Δ p g ¯ 和 Δ p a ¯ )随蓄热体长度(L)的变化情况。可以看出，切换时间(τ) = 30 s时，随着L从100 mm增加到200 mm， Δ p g ¯ 和 Δ p a ¯ 分别近乎线性地从83.9 Pa和59.9 Pa增加到155.7 Pa和125.3 Pa，显然，压力损失的这种线性增长是符合粘性流体层流流动的一般规律 [ 19 ] 。另一方面，压力损失与τ的关系较复杂，将τ从30 s延长到60 s， Δ p g ¯ 增加了3.6%，而 Δ p a ¯ 减少了7.9%。这种变化规律可以解释如下： 图3. 蓄热体长度(L)对传热及流阻性能的影响 延长τ，相同时间切换次数越少，流体换向的局部损失越少，这将使压力损失减小。另一方面，根据图3(a)，τ变大后，烟气出口温度上升，流体速度变大，这可能导致 Δ p g ¯ 增大，所以 Δ p g ¯ 随τ的变化取决于上述两种因素的综合效果；空气预热温度随τ增加而降低，因而 Δ p a ¯ 随τ的增加单调减小。 5.3. 结构优化与变流量特性 增加蓄热体长度，有助于提高余热回收的效率，但烟气出口温度( T ¯ g o )降低。 T ¯ g o 下降到限制温度时，蓄热体出口容易出现凝酸腐蚀。金属管式空气预热器的 T ¯ g o 一般不低于140℃~150℃ [ 20 ] 。本文采用的是耐高温热腐蚀的陶瓷材料， T ¯ g o 不低于120℃。考虑到蓄热体工作流量偏离额定值后，实际 T ¯ g o 可能进一步降低，为此，作者优选长度150 mm的蜂窝蓄热体，并据此针对不同空气流速(烟气流速相应变化)开展变流量特性研究。 图4(a)示出切换时间30s条件下空气预热温度( T ¯ a o )和烟气出口温度( T ¯ g o )随空气进口速度(u a = 1.56~2.34 m/s，对应80%~120%额定流量)的变化情况。从图中曲线可以看出，随着u a 从额定值的80%增加到120%， T ¯ a o 从622.7 K降低到590.7 K， T ¯ g o 则从396.9上升到424.5 K，高于前面定义的 T ¯ g o 下限(120℃或393.15 K)。与 T ¯ a o 相似，蓄热体效能(ε)和能量回收率(ERR)随u a 的增加而减少，见图4(b)。就本文所研究的变工况条件，蓄热体ε和ERR的相对变化均小于~9.5%。蓄热体烟气和空气压力损失随工作流量的变化情况见图4(c)所示。从图中曲线可以看到，随着u a 从额定值的80%增加到120%，两流体压力损失分别增加了53%和46%。 图4. 蓄热体变流量(空气速度u a )条件下传热和流动性能
system: [SYSTEM]

user: 6. 总结和结论 本文以低热值燃料蓄热式燃烧应用为背景，通过对CFD软件ANSYS Fluent二次开发，建立方孔蜂窝陶瓷蓄热体复杂传热的单元流道模型，并开展相关三维数值模拟研究，结果表明： 1) 随着蓄热体长度(L)的增加，空气预热温度( T ¯ a o )、蓄热体效能(ε)和能量回收率(ERR)均增大，烟气和空气压力损失也相应增加；减小切换时间， T ¯ a o 、ε和ERR升高，L越小，变化越显著。 2) 烟气出口温度( T ¯ g o )随L的增加而减小。在防止出口结露酸腐的前提下，150 mm蓄热体传热和流动综合性能最佳，此时 T ¯ a o 、ε和ERR分别为606.1 K、82.7%和66.6%，而烟气 T ¯ g o 和压力损失分别为411.2 K和117.9 Pa。 3) 在偏离额定流量~20%条件下， T ¯ g o 较额定值降低~14.5 K，ε和ERR则分别增加了~4.2%和3.4%。
system: [SYSTEM]

user: 摘 要 本文从僜人分布与族群由来、僜人村概况、居住地自然环境、民俗风情、近现代发展状况和今后发展思路等方面进行论述，主要是让大家进一步了解僜人、了解僜人的文化及发展情况，为今后僜人方面的系统研究提供可靠的科学依据。 关键词 :僜人，文化，发展情况 Copyright © 2019 by author(s) and Hans Publishers Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/
system: [SYSTEM]

user: 1. 引言 目前，我国学者尚未对僜人的发展现状开展过研究。只有个别学者对僜人的生存状况进行过少量报道 [ 1 ] 。通过此次调查，可以更好地促进僜人村的经济发展工作，为国家有关部门决策提供较为翔实的资料，为后来学者进一步开展人类学、民族学等研究提供最新的数据。
system: [SYSTEM]

user: 2 僜人分布与族群由来 僜人，又称僜巴，一个有待于国家正式承认为民族的部族。主要分布于西藏东南部林芝市察隅县，聚居于察隅河及其西支流贡日嘎布曲流域 [ 2 ] [ 3 ] ，总人数两万余人。但居住在我国实际控制线以北，察隅贡日嘎布曲到察隅河沙玛一带，上、下察隅两个镇9个自然村(沙琼、嘎腰、夏尼、新村、自更、京都、 体育、格拥、迟巴)的人总数仅2000余人 [ 4 ] 。 僜人有自己的语言，属汉藏语系藏缅语族，但无文字，通用藏语。分属“达让”和“格曼”两个语种。讲达让话的人自称“达让”，讲格曼话的人自称“格曼”。在察隅一些地区中，通用达让话。由于达让话里“美登”是穷人的意思，是与“莫北”即富人相对的称呼，因而汉语统称他们为“僜人”，又由于是族称，所以现在又把登写作“僜”。藏语中“巴”为“人”之意，故又有“僜巴”之称 [ 1 ] [ 4 ] 。
system: [SYSTEM]

user: 3. 下察隅镇僜人村概况 3.1. 下察隅镇概况 下察隅镇位于西藏自治区察隅县的南部，东与竹瓦根镇相邻，西北与上察隅镇相连，东南与缅甸相接，西南与印度接壤。全镇总面积5189.92平方千米，境内最高海拔5144米，最低海拔约600米，平均海拔约1600米。镇政府驻地洞冲位于北纬28˚30'，东经97˚1'，距察隅县城65千米，距上察隅镇45千米 [ 4 ] 。 2018年，下察隅镇辖扎巴、共同、自更、塔林、巴安通、新村、夏尼、嘎堆嘎美、沙琼、京都、塔玛、嘎腰、宗古、竹尼、日玛、拉丁、布巴、松古、沙玛和卡地等20个行政村，24个自然组，总户数1343户，人口5640人，主要居住藏、僜人、汉三个民族。经济以农业为主，是察隅县的粮食主产区，主要农作物有水稻、玉米、小麦、花生等。近年来，大力发展特色种植产业，主要有猕猴桃、石榴、茶叶等。 3.2. 僜人村概况 下察隅镇僜人村有6个，分别为沙琼村、夏尼村、嘎腰村、自更村、新村和京都村，有331户1340人；耕地面积2308.85亩(旱地1178.95亩，水田1129.9亩)。根据2018年下察隅镇国民经济统计数据，僜人村人均纯收入14，237.56元，人均现金收入11，674.8元。村民主要收入来源于种植业和国家政策性补贴。主要农作物有水稻、玉米、小麦、青稞、花生等。“十二五”至今，人均纯收入从2011年不足2300元提高到2018年的14，237.56元，6个僜人村人均纯收入增长超过500%。