汉代砖瓦窑的建造使用及数值模拟
2023-11-03张宜雯付小龙
温 睿,张宜雯,付小龙
[1. 中国-中亚人类与环境“一带一路”联合实验室(西北大学),陕西西安 710127;2. 文化遗产研究与保护技术教育部重点实验室(西北大学),陕西西安 710127;3. 西北大学文化遗产学院,陕西西安 710127; 4. 西安近代化学研究所,陕西西安 710065)]
0 引 言
新石器时代陶器的出现是中国陶瓷史的重要里程碑之一,从早期无窑烧制到裴李岗文化时期出现最早的横穴窑,陶器的质地得到了一定的提升;随着先民在生产生活上的需要,对于陶器的器型和使用范围需求不断增加,促使人们对陶窑的结构进行创新和改造。由于南北方的地理位置特征和原料选择的差异,南方自商代出现了依坡而建的小型龙窑,北方则以发展馒头窑为主。窑炉的改进促进了烧成温度的提高,为我国陶瓷发展提供了重要条件。汉代是制陶业发展的繁荣时期,砖瓦作为陶质建筑材料在汉代表现出卓越的工艺质量,这与汉代砖瓦窑结构的改进密切相关。半倒焰窑为汉代烧造陶器的主流窑炉,目前最早发现于商代[1],相较于升焰窑具有升温快、保温时间长、温度分布均匀等优点。经过升焰窑的结构变化和烧制技术改进,半倒焰窑在汉代发展趋于成熟并被广泛使用。半倒焰窑一般由前操作区、火膛、窑室、烟道等主要组成部分构成。根据汉代后北方黄河流域的窑址发掘显示,半倒焰窑的改进仅在规模和结构比例上得到体现,主要结构变化不大比较稳定。
目前已有学者对汉代陶窑的结构变化进行系统性的梳理研究[2-3],包括对汉代陶窑的建造形式、窑形结构、功能性质进行分类和总结概括,揭示了汉代陶窑结构的发展脉络,展现了汉代的窑业发展水平,但是未能深入探讨汉代陶窑性质与结构变化的科学逻辑。近年来,有学者运用数值模拟的方法复原陶窑内烧制过程中温度场变化的情况,建立了通过流体力学来研究窑炉结构科学性的方法,并对景德镇窑开展了研究。颜新华[4]通过对景德镇三种窑形变化分别进行模拟,结果表明从葫芦窑发展到蛋形窑温度分布更为均匀合理,并结合复原的葫芦窑测温实验验证了通过FLUENT软件模拟实验的可行性。本研究是在前人的基础上,采用数值模拟的方法对汉代数量最多的砖瓦窑烧成过程的温度场和流场分布进行模拟,进而讨论汉代砖瓦窑结构变化背后的科学道理,为从热能利用的角度研究中国陶窑的发展奠定基础。
1 汉代砖瓦窑发现概况和结构特点简述
目前为止见于公开发表资料的汉代陶窑主要发现集中于河南陕西省一带,黄河中下游地区汉代砖瓦窑目前已发现440余座,其中已发掘的窑为230座。陕西省最多为107座,河南省其次为85座,河北省发现22座,山西省6座,山东省10座。陶窑按照功能分类可分为日用陶器窑、砖瓦窑、烘范窑、陶俑窑、多功能窑等。其中,汉代砖瓦窑的基本结构包括窑门、火膛、窑室、烟道(图1),其他部分如操作间、挡火墙、排烟口存在差别。
图1 河南省焦作市苏蔺村汉代窑址Y6平、剖面图Fig.1 Plan and section of Kiln site Y6 of the Han Dynasties at Sulin Village, Jiaozuo, Henan Province
操作间因窑而异,根据发掘情况分为单窑操作间、双窑共用操作间和无操作间三种。单窑操作间是指各窑前的工作平台独立,各窑之间没有空间的联系,大部分砖瓦窑属于单窑操作间类型;双窑操作间是指相邻相对两座窑共用一个操作坑,双窑共用可见于陕西西安曲江春临村窑址Y2与Y6、Y4与Y5[5]、洛阳市东郊东汉砖瓦窑[6]等(图2),多窑共用操作间提高了工作效率。
汉代砖瓦窑大部分为单火膛结构,目前仅发现一座陶窑存在双火膛结构。河南省新乡市东干道发掘的一座双火膛对称的西汉陶窑[7],为半地穴式,窑床平面呈长方形,火膛位于窑室的东西两侧,排烟口和烟道在窑床南北方向各3个,中间的排烟口尺寸较大且烟道竖直向上,两侧的排烟口较小且烟道向中间靠拢(图3)。较单火膛而言,双火膛陶窑具有窑室温度高、分布均匀的优点。比较特殊的是在西安市未央区六村堡窑址发现的两座小规模窑址没有独立的火膛(图4),火门与窑室相接,窑室底部存在较多白色灰烬的扇形区域[8]。
图4 陕西西安汉长安城Y47Fig.4 Y47 at Han Chang’an City in Xi’an, Shaanxi Province
春秋早期的窑还存有升焰窑的结构,如洛阳东周王城内的古窑址Y5,窑床和火膛之间架设有窑箅[9]。战国中晚期至秦汉早期的砖瓦窑形制表现出逐渐统一的趋势,窑室的平面形状以圆形、梨形、卵圆形为主,还存有少量三角形、漏斗形的窑形,烟道则以单烟道为主,三个烟道数量的窑较少[2]。汉代砖瓦窑的火膛包括梯形、半圆形两种,窑床平面形状包括方形、长方形、椭圆形、梨形四种,火膛窑床面积比较秦汉早期大。根据汉代各窑址发掘报告和相关总结性研究来看汉代陶窑各结构发展规律如下,汉代砖瓦窑火膛以半圆形、梯形为主;西汉早期的窑室形状以圆形、椭圆形、梨形为主,西汉中晚期至东汉时期以近方形为主,表明烧造规模的扩大导致窑址结构的调整;西汉晚期陶窑的排烟系统则有汇合烟道和多烟道两种形式,其中砌筑方式由砖砌筑隔开形成烟道,早期则以生土掏挖夯实壁面作为烟道。同时烟道数量呈增多的趋势,除去单烟道的窑各个时期都有出现,大部分砖瓦窑烟道数量为3~5个。此外,根据李清临[10]的统计分析,汉代单烟道窑火膛、窑室、烟道的面积呈现正相关性,大部分多烟道和汇合烟道也呈正相关性。
为进一步研究汉代砖瓦窑结构变化的影响,针对单烟道窑和多烟道窑开展模拟研究,以期得到烟道数量对砖瓦窑的影响变化程度。由于大部分窑址在发掘或者使用时窑顶坍塌不存,甚至火膛、窑顶上部均不存,因此选择保存较好具有代表性的两个汉代陶窑进行建模模拟。单烟道以西汉前期的洛阳东周王城十一号窑(图5)为代表[9],多烟道以东汉时期的河南省焦作市苏蔺村汉代窑址Y6为代表[11],其中前者由于窑顶部分残存按照后者窑室高度比进行复原(表1)。根据窑内出土遗物多为建筑材料设定两窑均烧制尺寸为320 mm×140 mm×50 mm的汉砖,进口烟气物理性质、温度等设为定量参数,不同尺寸结构作为变量进行数值模拟。
表1 Y6和十一号窑的尺寸数据Table 1 Dimensional data of Kiln site Y6 and Kiln No.11
图5 多烟道、单烟道陶窑Fig.5 Multi-flue, single-flue pottery kilns
2 物理模型的建立及简化
为对汉代砖瓦窑进行复原模拟,首先需要对窑的建造方式、码窑方式、烧成方式进行总结研究,然后在此基础上得到初步的模型。
2.1 建造方式
多个窑址的简报提到了砖瓦窑的建筑方法[11-12],基本大致相同,为平地或斜坡上建造半地穴式窑,具体方法是先在生土中掏挖出窑室的平面形状如方形、椭圆形、梯形等,并在坑壁外留出一砖大小的边沿预留砌筑窑顶。然后在窑室前壁沿火膛平面形状如梯形、半圆弧形掏挖出一定深度的火膛,并在火膛前壁掏挖出窑门与操作坑相通。最后修整窑的基本形状夯实内壁,并涂抹草拌泥使其表面光滑,多次烘烤后内壁坚硬耐用。据已有的窑址显示窑门的砌筑方式主要是由土坯或砖券顶砌筑,如图6所示,先在生土中掏挖出拱形,再以双层土坯砌成土坯券顶横向平券,窑门下部一般有封火墙与火膛底齐平。根据已发掘的窑如河南省焦作武陟苗庄汉代陶窑等可判断,火膛一般是底部平整直壁掏挖而成,由火门向顶部缓慢延伸成穹窿顶形状至窑室上方地面,东西两壁为直壁。窑室后壁一般掏挖1~3个方形的排烟口,再由地面向下垂直掏挖烟道至与排烟口相通,此外还有多个烟道汇合中烟道的实例,掏挖此类烟道需要从中烟道对应地面的位置按照一定弧度向下掏挖,此类的烟道是由夯实的土壁隔开;西汉中晚期的汉墓陶窑中有以青砖隔开烟道的,此类排烟系统需要把烟道部分全部掏挖再以砖隔开砌成三条烟道。
图6 河南郑州马良寨遗址汉代陶窑Y3、Y4Fig.6 Han Dynasty Pottery kilns Y3 and Y4 at Malianzhai Ruins in Zhengzhou, Henan Province
窑室顶部保存下来的窑数量很少,河南省焦作市苏蔺村汉代窑址中发现一座保存较完整(图1),为四隅券进式穹窿顶结构(图7);从汉魏洛阳城东汉烧煤瓦窑Y1残存的窑壁顶砌砖来看,窑顶同样是四隅券进式结构[13],因此此种窑顶结构的窑址年代较晚均为东汉时期,而且目前的考古资料也能证明最早出现此种汉墓封顶方法在东汉时期[14]。此外西汉时期的窑也可据汉墓封顶方式推测封顶方式,推断为四面结顶式穹窿顶或者券顶结构。
图7 四隅券进式穹窿砌法示意图[15]Fig.7 Schematic diagram of quaripartite arched domes
2.2 码窑方式
汉代砖瓦窑窑床的平面较为规整,目前发掘的存有砖坯的窑址如新乡北站区前郭柳村汉代窑址Y3,窑床上撒有一层薄的草木灰,未经烧过的砖坯在窑内横立放垒起,窄面上皆沾有草木灰,推测是防止砖坯之间粘连所用,行与行之间摆放距离较密不足2 cm,砖坯与窑壁之间距离稍大约为5 cm[16]。汉魏洛阳城东汉烧煤瓦窑窑床上有一层排列整齐的东西向侧立砖坯(图8),文中指出这层砖坯是作为窑箅使用,行间距为5~10 cm,说明窑床底部铺设砖可能作为隔火道使用,为其他器物架设所用,河南郑州马良寨遗址Y1(图9)出土的窑床砖坯也出现了类似排列方式[17]。
图8 汉魏洛阳城东汉烧煤瓦窑Y1Fig.8 Coal-fired Tile kiln Y1 of the Eastern Han Dynasty at Han-Wei Luoyang City
图9 河南郑州马良寨遗址Y1Fig.9 Maliangzhai site Y1, Zhengzhou, Henan Province
以上结合历史文献记载、考古发掘资料及现代工业砖窑的多种码坯方式(图10)参考复原,即可对两窑模型进行建立,据东汉时期的河南省焦作市苏蔺村汉代窑址Y6考古报告显示,该窑窑顶为四隅券进式穹窿顶结构,而洛阳东周王城十一号窑为西汉早期的窑,因此窑顶可假设为四面结顶式穹窿顶,据此两窑建模完成后如图11所示。
图10 现代砖窑的部分码窑方式[18]Fig.10 Typical setting ways for modern brick kilns
图11 窑的物理模型Fig.11 Physical models of the kilns
2.3 烧成方式
明末宋应星的《天工开物》[19]中第十一章陶埏记载了砖瓦的制作过程和烧造过程,与汉代砖瓦窑最大的区别在于窑身均位于地表,而非汉代窑身大部分位于地下的结构;西汉中后期后砖窑形状结构逐渐稳定,至今我国部分地区还保存有类似的土窑结构,可见陶窑结构稳定后变化不大。《天工开物》中记载柴窑的烧造方法:“凡窑成坯之后,装入窑中。所装百钧则火力一昼夜,二百钧则倍时而足”“凡柴薪窑巅上侧凿三孔以出烟。火足止薪之候,泥固塞其孔,然后使水转釉。”关于浇水转釉过程的解释“凡转釉之法,窑巅作一平田样,四围稍弦起,灌水其上。砖瓦百钧用水四十石。水神透入土膜之下,与火意相感而成。水火既济,其质千秋矣。”即柴窑烧造后期,用泥封闭添柴入口,在窑顶开平面不断浇水冷却保持还原气氛,通过这种方法烧出来的砖更耐用,同时在窑顶采用浇水转釉的方法也证明了砖瓦窑从半地穴式转向地表式的这一过程。目前考古资料中最早出现这一结构的窑可追溯至战国[20],沿用至明清并在现代青砖生产中仍被使用。
2.4 模型的简化
本研究中数值模拟方法主要是通过使用ANSYS软件中FLUENT模拟窑室烟气的瞬态温度场和瞬态速度场,即烧窑后期封闭窑门窑内发生还原反应的阶段。由于仅对窑室进行研究,因此对操作间、窑门未进行建模,将燃料和火膛简化为带有一定进口速度和温度的烟气进入窑室加热器物制品,并忽略墙体壁厚,将窑壁、窑顶、窑底简化为具有一定热流密度的墙体,查取夯土的物性参数应用模拟;砖尺寸为320 mm×140 mm×50 mm,将其简化为具有一定吸热能力的墙体,查取青砖的物性参数应用模拟。根据火膛内残留的草木灰和时代背景判断燃料化学成分为木柴或秸秆,均属于传统生物质燃料,依据常见生物质燃料的主要元素组成[21],取木柴的主要化学成分均值作为此次模拟的燃料基本成分(表2)。本次使用ICEM CFD网格划分软件对窑进行非结构化网格划分,网格类型为四面体结构,数量为五百万左右。将生成后的网格文件导入FLUENT进行数值模拟,根据相关文献实地烧窑的资料,边界条件选择如下:烟气速度设为2 m/s,入口温度设为1 200 ℃(1 473 K),烟气的物理性质参数见表3[22];通过化学成分和经验公式计算可得入口烟气密度为1.56 kg/m3,出口边界压力出口,环境温度为300 K。由于烟气Re>105为湍流,湍流模拟方法设置为标准k-ε模型,算法设置为二阶迎风格式的SIMPLE算法。
表2 木柴的主要化学元素组成Table 2 Composition of the main chemical elements of firewood (%)
表3 大气压力下烟气热物理性质Table 3 Thermophysical properties of flue gas at the atmospheric pressure
3 数值模拟结果与分析
FLUENT仿真模拟完成后,通过后处理得到以下结果。
t=T-273.15 K
(1)
式中,t为摄氏温度,单位符号为℃;T为热力学温度,单位符号为K。
3.1 洛阳东周王城十一号窑模拟结果与分析
从图12中可以看出,洛阳东周王城十一号窑温度分布不均匀但差别不大,烧制砖的温度大致在1 000 K至1 400 K之间,出口烟气温度为1 250 K。从纵剖面的流线图来看,高温烟气由火膛升入窑室并穿过砖之间的空隙,最后汇聚于吸火孔通过烟道排出,符合半倒焰窑的特点;结合温度分布云图来看,高温烟气纵向传播范围由下至上(窑床至窑顶)依次减弱,横向传播范围由左至右(窑室前部至窑室后部)依次减弱,窑室前端受热温度最高热辐射最强,中端次之,后端最弱,但窑室前端烟气扰动最大,中端次之,后端最稳定;其中一部分高温烟气通过窑顶没有为烧成砖提供热量直接倒向吸火孔排出窑外,在单吸火孔部位产生局部高温区,总的来说,半倒焰窑的结构能够发挥很好的保温作用。
图12 洛阳东周王城十一号窑模拟结果Fig.12 Simulation results of Kiln No.11 at the Capital City site of the Eastern Zhou Dynasty in Luoyang
3.2 河南省焦作市苏蔺村汉代窑址Y6模拟结果与分析
从图13中看,河南省焦作市苏蔺村汉代窑址Y6的温度分布场特点与同样作为半倒焰窑的结构的洛阳东周王城十一号窑大体上差别不大,但是对比窑的纵向温度分布和同位置的平剖面温度云图,前者的高温烟气传播范围更广,窑室内的温度整体更高更均匀,在1 200 K至1 400 K之间,出口温度为1 250 K;由于Y6的吸火孔和烟道为三个,高温烟气在多个吸火孔的作用下可以更均匀地将高温烟气排出窑外,并在砖坯和窑顶之间反方向流动形成旋流,其中旋流能够增加烟气流动阻力、增强烟气扰流程度,从而达到增强烟气与砖坯之间对流换热能力的效果;同时旋流中心能形成低压区,有增加烟气流速的作用,因此Y6烟气流动速度比前者快,出口速度比前者慢,保证了砖坯能够充分地与高温烟气发生对流换热,使窑内温度分布更为均匀。
结合洛阳东周王城十一号窑和河南省焦作市苏蔺村汉代窑址Y6的模拟结果来看,汉代砖瓦窑的结构变化具有一定的合理性,特别是烟道由一个发展为多个的结构变化,改善了窑室内温度分布不均匀分布情况和保温效果,提高了砖坯的烧成温度范围和烧成质量,符合西汉中后期砖室墓和大型建筑兴起的需求,并为后续大规模烧制更高烧成温度条件的陶瓷器奠定了技术基础。
4 相关问题探讨
4.1 汉代砖瓦烧成温度的探究
砖瓦的生产技术源自陶器,它的坯体制作和烧成技术较陶器更为粗放。在实际生产过程中,由于砖瓦窑的窑室空间较陶窑更大,人们往往选择陶器和砖瓦建筑材料共烧的方式以提高产品生产效率,因此从制作技术和烧成技术层面看,砖瓦和陶器是密不可分的。我国古代陶器的烧成温度一般在800~1 000 ℃之间,据我国各时期部分砖瓦样品的烧成温度(表4)统计显示,砖瓦的烧成温度范围也在此之内。从各时期纵向对比来看,砖瓦的烧成温度在商周至东汉时期逐渐升高,在秦汉时期达到较高值,而后稳定在1 000 ℃左右。
砖瓦窑的窑形变化规律与其温度变化规律存在一定的相关性。单烟道窑作为商周至西汉早期的主要窑形代表,烧成温度较低,在850~1 000 ℃之间。西汉中晚期后,烟道数量增至三个及以上,温度也随之稳定在1 000 ℃左右,与模拟实验得到的多烟道窑内温度分布更均匀更稳定的结论相吻合。
4.2 挡火墙的初步研究
关于挡火墙的定义,南方窑址较多且与北方馒头窑有较大的区别,龙窑、阶级窑中的挡火墙设在窑室和烟道中间,墙上还开有排烟口,可以防止火焰流速过快使产品受热均匀,阶级窑中挡火墙还可用于分室,而半倒焰窑的挡火墙则设在火膛和窑室中间(图14)。早在秦陵周围发现的陶窑有此类似挡火墙的结构,包括下和村82下Y1的火膛有一道从底部筑起瓦砌的高出窑床约40 cm高的挡火墙,80西Y2的火膛后有细绳纹瓦片筑成的工字分火道,主要出土陶器以制陶工具、砖瓦、生活用具为主[25]。汉长安城西北部窑址中多个陶窑的火膛和窑室之间均有挡火墙,窑室中间纵向分布有分火道隔墙,其中Y22装满裸体陶俑坯体,隔火墙现存5层砖,宽0.2 m,残高0.46 m,通火孔约0.02~0.08 m,分火道将窑室南北隔开现存4层,残高0.38 m,通火孔约0.02~0.10 m,南北分别填满倒置的陶俑坯,文中据规模推测该窑群作为烧制西汉时期大型墓地的随葬陶俑所用,单窑内一次可装烧近四百个陶俑[26]。
图14 汉长安城西北部窑址Y22Fig.14 Kiln site Y22 in the northwest of Han Chang’an City
结合前文模拟实验的结果分析可知,这类结构的半倒焰窑在没有挡火墙的情况下,窑室前部最下端流经烟气温度最高,流速最快,对于烧制大批量的质量要求较高的器物较难满足要求,因此在下部设立一定高度具有通火孔的隔火墙一方面可以使窑室前部的制品不至于过热烧坏,一方面可以减缓底部烟气的流速不至于过快流出烟道,使更多的烟气往上流动,达到窑内温度更为均匀的效果。关于挡火墙对窑内温度场和流场的具体影响需要更进一步的模拟研究。
5 结 论
综上所述,根据前人研究及模拟实验结果可得到以下结论:
1) 汉代砖瓦窑的建造方式为半地穴式,掏挖出平面形状后夯实壁面并涂抹草拌泥,西汉时期窑顶的封顶方式为四面结顶式穹窿顶或者券顶,东汉时期窑顶的封顶方式为四隅券进式穹窿顶或者券顶,码窑方式为生坯侧立排列摆放,砖窑的烧成温度在1 000 ℃左右。
2) 汉代砖瓦窑的结构变化适应了当时社会对陶质材料和用具的扩大化需求,具有科学的合理性。通过模拟实验验证了半倒焰窑的烟气流动特点,特别是烟道由一个变为多个的结构变化使窑内温度分布更为均匀,提高了砖瓦的烧成温度和烧成质量。
3) 经过初步研究,汉代砖瓦窑在火膛和窑室之间设置挡火墙,起到了隔热和减缓高温烟气流速的作用,该结构对于烧制数量多质量较高的器物如陶俑具有较大的意义。