基于CFD模拟探究自然通风对日光温室小气候的影响

基于CFD模拟探究自然通风对日光温室小气候的影响一、引言1.1研究背景与目的1.1.1研究背景随着全球人口的增长和人们对农产品需求的不断增加，高效、可持续的农业生产变得愈发重要。日光温室作为一种重要的农业生产设施，能够利用自然光照和太阳能，为农作物提供适宜的生长环境，在现代农业中占据着举足轻重的地位。它不仅可以延长农作物的生长周期，增加农作物的产量，还能提高农作物的品质，满足市场对多样化农产品的需求。例如，在北方寒冷地区，日光温室使得冬季种植蔬菜成为可能，丰富了人们的餐桌。在日光温室的运行过程中，小气候环境对作物的生长发育起着关键作用。适宜的温度、湿度、气流速度和二氧化碳浓度等条件，能够促进作物的光合作用、呼吸作用和蒸腾作用，有利于作物的健康生长和高产优质。而自然通风作为调节日光温室小气候的一种重要且经济环保的手段，具有不可替代的作用。通过自然通风，可以引入室外新鲜空气，排出室内污浊空气，调节温室内的温度和湿度，补充二氧化碳，改善室内空气品质，为作物生长创造良好的环境。同时，自然通风无需消耗大量能源，符合可持续发展的理念，对于降低农业生产成本、减少环境污染具有重要意义。近年来，随着计算机技术的飞速发展，计算流体动力学（CFD）模拟技术在农业工程领域的应用越来越广泛。CFD模拟技术能够通过数值计算的方法，对流体的流动、传热和传质等物理现象进行精确模拟，为研究日光温室自然通风提供了有力的工具。与传统的实验研究方法相比，CFD模拟具有成本低、周期短、可重复性强等优点，可以在虚拟环境中对不同的通风方案进行快速评估和优化，大大提高了研究效率和准确性。通过CFD模拟，可以深入了解自然通风条件下日光温室内气流场、温度场和湿度场的分布规律，为温室的设计、改造和管理提供科学依据。1.1.2研究目的本研究旨在利用CFD模拟技术，深入探究自然通风对日光温室小气候的影响，具体目标如下：揭示小气候参数变化规律：通过CFD模拟，详细分析自然通风过程中日光温室内温度、湿度、气流速度等小气候参数的时空变化规律，明确不同通风条件下这些参数的分布特征，为深入理解自然通风对温室小气候的影响机制提供理论支持。分析通风口位置与面积的影响：研究通风口的位置和面积对自然通风效果及小气候环境的影响，找出最优的通风口设置方案，以提高自然通风效率，改善温室内小气候的均匀性，为日光温室的通风设计提供科学指导。为温室设计和管理提供依据：综合考虑自然通风对日光温室小气候的影响，结合作物生长需求，提出合理的温室设计和管理建议，包括温室结构优化、通风系统调控等方面，以实现日光温室的高效、节能和可持续生产，提高农作物的产量和品质。1.2国内外研究现状1.2.1自然通风的研究现状自然通风作为一种传统且环保的通风方式，在建筑领域有着广泛的研究。其原理主要基于风压和热压作用。风压是由于大气中的压力差，当风吹过建筑物时，在迎风面和背风面形成压力差，促使空气从迎风面的窗缝和其他空隙流入室内，而室内空气则从背风面孔口排出，从而形成全面换气的风压自然通风，其效果受到建筑的几何形状、相对于风向的方位、风速以及建筑周围自然地形等因素影响。热压则是由室内外空气的温度差引起，当室内温度高于室外时，建筑物上部压力较高，下部压力较低，空气通过下部开口进入，从上部流出；反之，气流方向相反，热压大小取决于两个开口处的高度差和室内外的空气温度差。在实际建筑中，自然通风通常是风压和热压共同作用的结果。在绿色建筑领域，自然通风被视为一项关键技术。它不仅能够有效地实现室内环境的降温，还能节约常规能源、减少环境污染，同时将新鲜空气引入室内并及时排出污染物，极大地改善室内环境品质。有调查显示，人们对自然风存在偏爱，即使在有点热的环境中，多数人还是选择自然环境而非空调环境。因此，自然通风有着空调所无法取代的优点。然而，自然通风技术在实际应用中仍面临一些问题，例如在不同气候条件下如何优化通风设计，以确保室内环境的舒适性和通风效率；以及如何解决自然通风受外界环境因素影响较大的问题，如风向、风速的不稳定等。1.2.2日光温室小气候的研究现状日光温室小气候涵盖了温室内的温度、湿度、气流速度和二氧化碳浓度等多个关键因素，这些因素相互关联、相互影响，共同对作物的生长发育起着决定性作用。适宜的小气候条件能够为作物创造良好的生存环境，促进其光合作用、呼吸作用和蒸腾作用的顺利进行，进而实现作物的健康生长和高产优质。国内外众多学者围绕日光温室小气候展开了丰富的研究。研究内容涉及不同地区日光温室小气候的特点分析，以及温室结构、覆盖材料、通风方式等因素对小气候的影响。例如，在温室结构方面，墙体材料、温室跨度和高度等的变化都会对温室内的温度分布产生影响；不同的覆盖材料，其透光率、保温性能等有所差异，进而影响温室内的光照和温度条件。通风方式也是影响小气候的重要因素，顶通风与底通风相结合的方式可以促进温室内空气流动，有效降低温室内温度和湿度，并及时补充室内二氧化碳浓度，满足植物的生长需求。然而，目前的研究在小气候的精准调控方面仍存在不足，对于如何根据不同作物的生长阶段和需求，实现小气候参数的动态优化调控，还需要进一步深入研究。1.2.3CFD模拟技术在日光温室自然通风研究中的应用现状随着计算机技术的飞速发展，CFD模拟技术在农业工程领域，尤其是日光温室自然通风研究中的应用日益广泛。CFD技术通过计算机模拟流体的运动、热量传递等物理现象，能够对日光温室内的气流场、温度场和湿度场等进行精确模拟，为研究日光温室自然通风提供了有力的工具。在相关研究中，学者们利用CFD技术对不同通风口位置和面积的日光温室进行模拟，分析通风效果及小气候环境的变化。通过模拟可以直观地看到温室内空气的流动路径、速度分布以及温度和湿度的分布情况，从而为优化通风设计提供科学依据。有研究通过CFD模拟发现，合理设置通风口可以使温室内空气流动更加均匀，降低温度梯度，提高通风效率。同时，CFD模拟还可以与实际测试相结合，验证模拟结果的准确性，进一步完善对日光温室自然通风的研究。不过，目前CFD模拟在日光温室自然通风研究中，对于复杂的温室结构和作物冠层的模拟还存在一定局限性，如何更准确地考虑作物与环境之间的相互作用，以及如何提高模拟的计算效率和精度，仍是需要解决的问题。1.3研究意义与创新点1.3.1理论意义本研究运用CFD模拟技术，对自然通风条件下日光温室小气候进行深入研究，在理论层面具有重要意义。通过模拟分析，能够更加全面、细致地揭示日光温室内温度、湿度、气流速度等小气候参数在自然通风过程中的时空变化规律。这不仅有助于完善温室环境调控理论，填补目前在自然通风与温室小气候相互作用机制研究方面的部分空白，还能为进一步深入研究温室环境系统提供坚实的理论基础。在过往的研究中，对于日光温室小气候的研究虽然涉及多个方面，但在自然通风对小气候影响的系统性和深入性研究上仍有不足。本研究通过CFD模拟，能够精确地量化不同通风条件下小气候参数的变化，为构建更加完善的温室环境调控理论体系提供关键的数据支持和理论依据。此外，本研究中CFD模拟技术的应用，丰富了CFD模拟在农业工程领域尤其是日光温室研究中的应用案例。通过探索CFD模拟在日光温室自然通风研究中的具体应用方法和技巧，为后续相关研究提供了可借鉴的研究思路和方法，推动CFD模拟技术在该领域的进一步发展和应用，促进多学科交叉融合，为解决农业工程领域的复杂问题提供新的视角和方法。1.3.2实践意义从实践角度来看，本研究成果具有广泛的应用价值。准确掌握自然通风对日光温室小气候的影响规律，能够为日光温室的设计提供科学指导。通过CFD模拟不同通风口位置和面积对小气候的影响，可以优化通风系统设计，提高自然通风效率，改善温室内小气候的均匀性，减少温度和湿度的梯度差异，为作物生长创造更加适宜的环境。合理的通风设计还能降低温室建设成本，提高资源利用效率，促进日光温室产业的可持续发展。对于温室管理者而言，研究结果有助于制定更加科学合理的管理策略。根据模拟结果，可以根据不同的天气条件、作物生长阶段和小气候参数的变化，精准地调控通风系统，实现温室内环境的动态优化。在高温天气下，通过合理开启通风口，增加通风量，降低室内温度，避免作物遭受热害；在高湿天气下，及时通风排湿，减少病虫害的发生几率。这样可以提高作物的产量和质量，增加农民的经济收入。自然通风作为一种经济环保的通风方式，利用本研究成果优化自然通风系统，能够减少对机械通风等耗能设备的依赖，降低能源消耗和生产成本，符合可持续农业发展的要求，对于推动农业绿色发展、实现节能减排目标具有积极意义。1.3.3创新点本研究在多个方面具有创新之处。在研究方法上，采用CFD模拟与理论分析相结合的方式，突破了传统研究方法的局限性。CFD模拟能够直观地展示日光温室内气流、温度和湿度等参数的分布情况，而理论分析则为模拟结果提供了深入的解释和理论支撑，两者相辅相成，提高了研究结果的准确性和可靠性。在模型构建方面，综合考虑了多种因素对自然通风和日光温室小气候的影响，建立了更加全面和准确的CFD模型。不仅考虑了温室的结构、通风口的位置和面积等常规因素，还将作物冠层对气流和热量传递的影响纳入模型中，更加真实地反映了实际的温室环境。这种多因素耦合的模型构建方法，能够更准确地预测自然通风条件下日光温室小气候的变化，为温室设计和管理提供更具针对性的建议。在影响因素分析上，本研究创新性地对不同通风工况下的多因素耦合作用进行了深入研究。通过改变通风口的开启方式、室外风速和风向等参数，分析这些因素之间的相互作用对日光温室小气候的综合影响，为全面理解自然通风机制提供了新的视角。这种对多因素耦合作用的系统研究，有助于发现自然通风过程中的关键影响因素和规律，为优化通风设计和调控策略提供了科学依据。二、自然通风与日光温室小气候相关理论2.1自然通风原理2.1.1热压通风原理热压通风，又称“烟囱效应”，其形成机制基于热力学基本原理。当日光温室内部空气因太阳辐射、作物呼吸作用以及土壤蓄热等因素被加热时，空气分子获得更多能量，运动加剧，导致空气密度减小。根据阿基米德原理，密度较小的热空气会在浮力作用下上升，趋向温室顶部。此时，若温室顶部设有通风口，热空气便会通过该通风口排出室外。与此同时，室外相对低温且密度较大的冷空气，在压力差的驱动下，会从温室底部的通风口或其他开口处流入室内。如此循环往复，便形成了由热压驱动的空气自然循环流动，不断更新室内空气，带走热量与污浊气体。热压的大小与多个因素密切相关。其中，通风口之间的高度差起着关键作用，高度差越大，热空气上升的动力就越强，热压通风效果也就越显著。室内外的空气温度差同样至关重要，温度差越大，空气密度差异越明显，热压作用就越强烈。以常见的日光温室为例，若在夏季晴朗的中午，温室内温度可达35℃，而室外温度为25℃，此时室内外较大的温度差会促使热压通风高效运行，大量热空气迅速从顶部通风口排出，新鲜冷空气不断从底部进入，有效降低室内温度。2.1.2风压通风原理风压通风的原理基于空气动力学。当风吹向日光温室时，在温室的迎风面，气流受到阻挡，速度降低，动能转化为静压能，使得迎风面的气压升高，形成正压区。而在温室的背风面，气流由于绕过温室而形成局部涡流，导致空气流速增加，静压降低，形成负压区。这种迎风面与背风面之间的压力差，即为风压通风提供了动力。在风压的作用下，室外空气会从迎风面的通风口或缝隙进入温室内，室内空气则会从背风面的通风口排出，从而实现空气的交换和通风换气。通风强度与正压侧和负压侧的开口面积以及风力大小密切相关。开口面积越大，空气流通的通道越顺畅，通风量就越大；风力越大，风压也就越大，通风效果越明显。例如，在有风的天气里，当风速为5m/s时，若温室迎风面和背风面的通风口总面积为5平方米，室内外的空气交换量会明显大于通风口总面积为2平方米时的情况，能够更有效地改善室内空气品质。2.1.3热压与风压共同作用机制在实际的日光温室自然通风过程中，热压通风和风压通风通常并非孤立存在，而是相互影响、共同作用，推动温室内外的空气交换。这种共同作用的机制较为复杂，受到多种因素的综合影响。当室外有风且温室内外存在温度差时，热压和风压会相互叠加。在一些情况下，热压和风压的作用方向可能一致，此时它们会相互增强，显著提高通风效率。在春季的晴天，室外微风吹拂，温室内因太阳辐射温度升高，热压和风压都促使空气从温室底部进入，从顶部排出，两者的协同作用使得室内空气能够快速更新，温度和湿度得到有效调节。然而，在另一些情况下，热压和风压的作用方向可能相反，此时它们会相互削弱，降低通风效果。在夏季的傍晚，室外风力较大，但温室内外温度差较小，风压可能促使空气从温室一侧进入，从另一侧排出，而热压的作用相对较弱，甚至可能与风压方向相反，导致通风效率下降。此外，热压和风压的共同作用还受到温室结构、通风口位置和面积等因素的影响。合理设计温室的结构和通风口布局，能够充分利用热压和风压的协同作用，优化通风效果。将顶部通风口和底部通风口合理设置，使得热压和风压在不同情况下都能有效地促进空气流通，避免两者相互抵消，从而为作物生长创造良好的室内环境。2.2日光温室小气候构成要素2.2.1温度温度作为影响作物生长发育的关键环境因子之一，对作物的生理过程起着至关重要的作用。不同作物在各个生长阶段，如种子萌发、幼苗生长、开花结果等，都对温度有着特定的要求。适宜的温度能够确保作物各项生理活动的正常进行，例如，在种子萌发阶段，合适的温度能促进种子内酶的活性，加快种子的吸水膨胀和营养物质的转化，从而顺利发芽。对于黄瓜来说，其种子萌发的适宜温度一般在25-30℃之间，若温度过低，种子萌发会受到抑制，甚至可能导致烂种；温度过高，则可能使种子失去活力。在幼苗生长阶段，适宜的温度有助于根系的生长和对养分的吸收，促进茎叶的健壮发育。番茄幼苗生长的适宜温度白天一般为20-25℃，夜间为10-15℃，在此温度范围内，番茄幼苗能够茁壮成长，叶片翠绿，茎干粗壮。当温度不适宜时，作物的生长发育会受到严重影响，可能导致生长缓慢、发育不良、病虫害易感性增加等问题。在低温环境下，作物的新陈代谢减缓，根系对养分和水分的吸收能力下降，容易出现叶片发黄、枯萎等现象；高温则可能引发作物的热害，导致光合作用受阻，呼吸作用增强，消耗过多的光合产物，影响作物的产量和品质。自然通风在调节日光温室内温度方面发挥着关键作用。通过热压通风和风压通风，自然通风能够实现室内外空气的交换，从而有效调节室内温度。在白天，随着太阳辐射的增强，温室内温度迅速升高。此时，热压通风作用下，温室内热空气上升并通过顶部通风口排出，室外相对低温的空气从底部通风口进入，带走温室内的热量，降低室内温度。若室外有风，风压通风也会协同作用，进一步增强通风效果，加速室内热量的散失。在夏季晴朗的中午，当温室内温度高达35℃以上时，合理开启通风口，利用自然通风可使室内温度降低3-5℃，为作物生长创造较为适宜的温度环境。然而，自然通风对温室内温度分布的影响较为复杂，受到通风口位置、面积以及温室结构等多种因素的制约。通风口位置的不同会导致室内空气流动路径和速度的差异，进而影响温度分布。若通风口设置不合理，可能会出现通风死角，导致局部区域温度过高或过低。通风口面积过小，通风量不足，无法有效降低室内温度；通风口面积过大，则可能导致室内温度波动过大，不利于作物生长。温室结构，如温室的跨度、高度、墙体材料等，也会对自然通风的效果产生影响，进而影响温室内的温度分布。跨度较大的温室，空气流通相对困难，容易出现温度不均匀的情况；而墙体保温性能好的温室，在自然通风时，热量散失相对较慢，室内温度下降幅度较小。2.2.2湿度湿度在作物生长过程中扮演着不可或缺的角色，它对作物的蒸腾作用、光合作用以及病虫害的发生发展都有着显著的影响。作物通过蒸腾作用，将水分从根部吸收并运输到叶片，再通过气孔散发到大气中，这个过程不仅能够帮助作物吸收和运输养分，还能调节作物的体温。适宜的空气湿度能够保证作物蒸腾作用的正常进行，维持作物的水分平衡。一般来说，大多数作物适宜的空气相对湿度在60%-80%之间。当空气湿度过低时，作物蒸腾作用过强，会导致水分散失过快，使作物出现萎蔫现象，影响光合作用和生长发育；空气湿度过高，则会使作物叶片表面凝结水珠，为病原菌的滋生和传播提供有利条件，增加病虫害的发生几率，如黄瓜在高湿环境下容易发生霜霉病、白粉病等病害。自然通风是调节日光温室内湿度的重要手段之一。其主要通过引入相对干燥的室外空气，与室内潮湿空气进行混合和交换，从而降低室内空气湿度。在通风过程中，室内潮湿空气被排出室外，同时室外干燥空气进入室内，使室内空气的水汽含量减少，湿度降低。在夏季高温高湿的天气条件下，温室内相对湿度常常会超过90%，此时开启自然通风，能够有效地将室内湿度降低到适宜的范围内。通过自然通风，还能促进室内空气的流动，加速水分的蒸发和扩散，进一步降低湿度。在温室内，水分主要来源于土壤蒸发和作物蒸腾，自然通风可以加快这些水分的排出速度，保持室内湿度的稳定。自然通风对湿度的调节效果同样受到多种因素的影响。通风口的大小和开启时间直接决定了通风量的大小，进而影响湿度调节效果。通风口开启时间过短或通风量不足，无法及时排出室内潮湿空气，导致湿度下降缓慢；通风口开启时间过长或通风量过大，则可能使室内湿度过低，对作物生长产生不利影响。室外空气的湿度和温度也会对自然通风的湿度调节效果产生重要影响。当室外空气湿度较高时，自然通风降低室内湿度的效果会受到一定限制；而室外温度较低时，通风可能会导致室内温度下降过快，影响作物生长。此外，温室的密封性、作物的种植密度以及灌溉方式等因素也会间接影响自然通风对湿度的调节作用。密封性好的温室，通风时湿度变化相对较小；种植密度大的温室，作物蒸腾作用强，产生的水汽多，需要更大的通风量来调节湿度；不合理的灌溉方式，如过量灌溉，会增加室内湿度，加大自然通风调节湿度的难度。2.2.3气流速度气流速度对日光温室内的空气流通、热量传递和气体交换过程有着深远的影响，进而对作物的生长发育产生重要作用。适宜的气流速度能够促进温室内空气的均匀混合，避免局部区域出现温度过高或过低、湿度过大或过小的情况，为作物创造一个相对稳定且适宜的生长环境。气流能够将作物周围的二氧化碳及时补充，满足作物光合作用的需求，同时带走作物呼吸作用产生的二氧化碳和水汽，维持良好的气体环境。当气流速度过小时，温室内空气流通不畅，容易导致二氧化碳分布不均，影响作物的光合作用效率；而气流速度过大，则可能对作物造成机械损伤，影响作物的正常生长。自然通风是改变日光温室内气流速度的主要方式之一。通过热压通风和风压通风，室外空气进入温室内，形成一定的气流。通风口的位置和面积是影响气流速度和分布的关键因素。不同位置的通风口会引导气流在温室内形成不同的流动路径，从而导致气流速度在温室内的分布不均匀。顶部通风口和底部通风口同时开启时，会形成上下对流的气流，使温室内不同高度的空气得到较好的混合；而仅开启顶部通风口或底部通风口时，气流分布则会相对单一。通风口面积越大，通风量越大，气流速度也会相应增加。但过大的通风口面积可能会导致气流速度过快，对作物生长不利，因此需要根据实际情况合理选择通风口面积。温室内的作物冠层也会对气流速度产生显著影响。作物冠层具有一定的阻力，会使气流在通过时速度降低。作物的种植密度、高度和叶片形态等都会影响作物冠层对气流的阻力。种植密度大、植株高、叶片宽大的作物，其冠层对气流的阻力较大，气流速度在作物冠层内下降明显；而种植密度小、植株矮、叶片狭长的作物，冠层对气流的阻力相对较小，气流速度下降幅度较小。在自然通风条件下，气流速度在作物冠层内的变化会影响冠层内的气体交换和热量传递，进而影响作物的光合作用和蒸腾作用。合理调整作物种植密度和布局，能够优化温室内的气流分布，提高自然通风效果，促进作物生长。2.2.4二氧化碳浓度二氧化碳作为作物进行光合作用的重要原料，在作物的生长发育过程中具有不可替代的作用。光合作用是作物将光能转化为化学能，合成有机物质的关键过程，而二氧化碳的充足供应是保证光合作用正常进行的前提条件。在适宜的光照、温度和水分等条件下，增加二氧化碳浓度能够显著提高作物的光合速率，促进作物的生长和发育，增加作物的产量和品质。研究表明，在一定范围内，二氧化碳浓度的增加可使作物的光合产物积累增加，从而提高作物的产量。对于黄瓜、番茄等蔬菜作物，当二氧化碳浓度从大气中的300-400μmol/mol增加到800-1000μmol/mol时，其产量可提高20%-50%。自然通风在调节日光温室内二氧化碳浓度方面起着重要作用。在自然通风过程中，室外含有较高浓度二氧化碳的新鲜空气进入温室内，补充室内因作物光合作用而消耗的二氧化碳，维持室内二氧化碳浓度的相对稳定。在白天，随着作物光合作用的进行，温室内二氧化碳浓度会迅速下降，若不及时补充，会限制作物的光合作用。通过自然通风，能够及时将室外的二氧化碳引入室内，满足作物的生长需求。当室外风速为3m/s，通风口面积为5平方米时，在通风1小时后，温室内二氧化碳浓度可从200μmol/mol迅速上升到350μmol/mol左右，接近大气中的二氧化碳浓度水平，为作物光合作用提供了充足的原料。然而，自然通风对温室内二氧化碳浓度的调节效果并非一成不变，它受到多种因素的制约。通风量的大小直接影响二氧化碳的补充速度和浓度调节效果。通风量越大，单位时间内进入温室内的新鲜空气越多，二氧化碳的补充量也越大，能够更快地提高室内二氧化碳浓度；反之，通风量过小，则二氧化碳补充不足，难以满足作物光合作用的需求。室外空气中二氧化碳浓度的变化也会对自然通风的调节效果产生影响。在不同的时间和地点，室外空气中二氧化碳浓度可能存在差异，如在城市郊区，由于人类活动和工业排放等因素，室外二氧化碳浓度可能会相对较高；而在偏远山区，室外二氧化碳浓度则更接近自然水平。此外，温室内作物的种类、种植密度和生长阶段等因素也会影响二氧化碳的消耗速度和对通风调节的需求。种植密度大、生长旺盛的作物，对二氧化碳的需求量较大，需要更大的通风量来补充二氧化碳。2.3CFD模拟技术基础2.3.1CFD模拟的基本原理CFD模拟的核心是基于数值求解流体力学的基本方程，主要包括连续性方程、动量方程（纳维-斯托克斯方程，N-S方程）和能量方程。这些方程描述了流体流动过程中的质量守恒、动量守恒和能量守恒定律。连续性方程从质量守恒的角度出发，反映了在流体运动过程中，单位时间内流入和流出控制体的质量差与控制体内质量变化率之间的关系。其数学表达式在笛卡尔坐标系下为：\frac{\partial\rho}{\partialt}+\frac{\partial(\rhou_i)}{\partialx_i}=0，其中\rho表示流体密度，t为时间，u_i是速度矢量在x_i方向的分量。动量方程，即N-S方程，是CFD模拟的关键方程，它体现了动量守恒定律。该方程描述了作用在流体微元上的力（包括压力、粘性力和惯性力等）与流体微元动量变化之间的关系。在笛卡尔坐标系下，其一般形式为：\rho(\frac{\partialu_i}{\partialt}+u_j\frac{\partialu_i}{\partialx_j})=-\frac{\partialp}{\partialx_i}+\frac{\partial\tau_{ij}}{\partialx_j}+\rhog_i，其中p是压力，\tau_{ij}为粘性应力张量，g_i是重力加速度在x_i方向的分量。能量方程则基于能量守恒原理，它考虑了流体的内能、动能以及由于热传导、对流和辐射等方式引起的能量传递。在考虑热传导和粘性耗散的情况下，能量方程在笛卡尔坐标系下的表达式为：\rho(\frac{\partialh}{\partialt}+u_j\frac{\partialh}{\partialx_j})=\frac{\partial}{\partialx_j}(\lambda\frac{\partialT}{\partialx_j})+\tau_{ij}\frac{\partialu_i}{\partialx_j}+S_h，其中h是比焓，\lambda为热导率，T是温度，S_h表示其他热源项。在实际的CFD模拟过程中，由于这些方程通常是非线性且复杂的偏微分方程，难以获得解析解。因此，需要采用数值方法将连续的计算区域离散化，把偏微分方程转化为代数方程组，通过迭代求解这些代数方程组来获得流场的数值解。常见的离散化方法包括有限差分法、有限体积法和有限元法等。有限差分法是将连续的偏微分方程通过差分近似转化为代数方程组，利用离散节点上的函数值来近似导数；有限体积法基于守恒定律，将计算域划分为有限数量的小体积单元，对每个单元应用守恒定律，以积分形式表达物理定律；有限元法则将计算域划分为由节点和单元组成的网格，并通过最小化一个能量泛函来求解控制方程。通过这些离散化方法，CFD模拟能够实现对流体流动和传热过程的数值模拟，得到流场中各物理量（如速度、压力、温度等）在空间和时间上的分布情况，从而为研究日光温室自然通风对小气候的影响提供了有力的工具。2.3.2CFD模拟在农业领域的应用现状在农业领域，CFD模拟技术近年来得到了广泛的应用，为解决农业工程中的诸多问题提供了新的思路和方法。在温室环境研究方面，CFD模拟已成为分析温室内气流分布、温度场和湿度场变化的重要手段。通过建立温室的CFD模型，可以模拟不同通风条件下温室内的空气流动情况，预测温度和湿度的分布规律。学者们利用CFD模拟研究了不同通风口设置对温室内气流和温度分布的影响，发现合理的通风口布局可以有效改善温室内的通风效果，使温度分布更加均匀，为温室的通风设计和环境调控提供了科学依据。CFD模拟还可以考虑温室内部作物冠层的影响，研究作物与空气之间的热量和质量交换过程，进一步优化温室环境，满足作物生长的需求。在畜禽舍通风设计中，CFD模拟同样发挥着重要作用。畜禽舍内的空气质量直接影响畜禽的健康和生产性能，良好的通风是保证畜禽舍空气质量的关键。利用CFD模拟可以对畜禽舍的通风系统进行优化设计，分析不同通风方式和通风参数下舍内的气流速度、温度和有害气体浓度分布，找出最佳的通风方案，降低舍内氨气、硫化氢等有害气体的浓度，提高畜禽的生活环境质量，减少疾病的发生，提高养殖效益。在灌溉系统设计中，CFD模拟可以用于研究喷头的喷雾特性和灌溉水在田间的分布情况。通过模拟不同喷头类型、工作压力和布置方式下的喷雾过程，可以优化喷头的选型和布置，提高灌溉水的利用效率，实现精准灌溉，减少水资源的浪费。然而，CFD模拟在农业领域的应用也存在一些局限性。一方面，农业环境往往较为复杂，涉及到多种物理过程和因素的相互作用，如温室中的作物生长、土壤水分蒸发、畜禽舍内的生物热产生等，准确地将这些因素纳入CFD模型中仍具有一定的难度，可能导致模拟结果与实际情况存在一定偏差。另一方面，CFD模拟对计算资源的要求较高，尤其是对于复杂的农业系统模型，计算时间较长，限制了其在一些实时性要求较高的应用场景中的应用。2.3.3常用CFD软件介绍目前，市场上有多种CFD软件可供选择，它们各自具有不同的特点和适用场景。以下对几款常用的CFD软件进行简要介绍：Fluent：是一款功能强大、应用广泛的CFD软件，由美国ANSYS公司开发。它具有丰富的物理模型库，能够模拟各种复杂的流体流动现象，包括层流、湍流、多相流、传热传质等。Fluent支持多种网格划分方式，能够适应不同形状和复杂程度的计算区域，并且提供了友好的用户界面和后处理功能，方便用户进行模型设置、计算求解和结果分析。在农业领域，Fluent被广泛应用于温室、畜禽舍等环境的模拟研究，能够准确地预测气流、温度和湿度等参数的分布，为农业设施的设计和优化提供了有力支持。CFX：也是一款知名的CFD软件，由英国AEATechnology公司开发，后被ANSYS公司收购。CFX采用了先进的数值算法，在处理复杂几何形状和多物理场耦合问题时具有较高的精度和计算效率。它具备强大的并行计算能力，能够大大缩短计算时间，适用于大规模的CFD模拟。CFX还提供了丰富的边界条件和求解器选项，用户可以根据具体问题进行灵活设置。在农业工程中，CFX常用于研究灌溉系统的水流特性、温室通风的多因素耦合作用等复杂问题。STAR-CCM+：是CD-adapco公司开发的一款综合性CFD软件，以其高度的集成性和强大的功能而受到用户的青睐。STAR-CCM+采用了独特的多面体网格技术，能够在复杂几何模型上生成高质量的网格，减少网格数量，提高计算效率。它还具备丰富的物理模型和求解器，支持多相流、化学反应、热管理等多种物理过程的模拟。在农业领域，STAR-CCM+可用于模拟农业机械的空气动力学性能、农产品干燥过程中的传热传质等问题。本研究选择Fluent软件进行日光温室自然通风对小气候影响的CFD模拟，主要基于以下依据：Fluent软件具有广泛的应用案例和成熟的技术支持，在农业领域的温室环境模拟方面已经取得了大量的研究成果，其可靠性和准确性得到了充分验证。Fluent软件的物理模型库能够较好地满足本研究中对自然通风、传热传质以及作物冠层影响等多物理过程的模拟需求。Fluent提供了丰富的网格划分工具和边界条件设置选项，能够方便地对日光温室的复杂结构进行建模和模拟。其友好的用户界面和强大的后处理功能，也有助于提高研究效率，便于对模拟结果进行直观的分析和展示。三、CFD模拟模型的建立与验证3.1日光温室模型构建3.1.1温室结构参数确定本研究选取某位于[具体地理位置]的实际日光温室作为案例进行研究。该日光温室主要用于蔬菜种植，为了确保CFD模拟的准确性，对其关键结构参数进行了详细测量与确定。通过实地测量，该日光温室的长度为60m，采用钢卷尺沿温室的长轴方向进行测量，测量精度精确到0.1m。宽度为10m，同样使用钢卷尺在温室的短轴方向进行测量。脊高为3.5m，利用激光测距仪从地面垂直测量至温室屋脊的最高点，以保证测量的准确性。墙体厚度方面，后墙和山墙均采用砖石结构，后墙厚度为0.5m，山墙厚度为0.4m，使用游标卡尺在墙体的不同位置进行多点测量，取平均值作为墙体厚度。屋面坡度通过测量屋面的高度和水平投影长度计算得出。屋面的高度为从温室屋脊到前底角的垂直距离，测量值为2.5m；水平投影长度为温室跨度的一半，即5m。根据三角函数关系，屋面坡度的正切值为屋面高度与水平投影长度的比值，即tan\theta=\frac{2.5}{5}=0.5，通过反正切函数计算得出屋面坡度\theta\approx26.6^{\circ}。此外，温室的前屋面采用塑料薄膜覆盖，透光率约为85%，通过专业的透光率测量仪在不同位置对塑料薄膜进行测量，取平均值得到该透光率数据。后屋面采用保温材料覆盖，保温材料为聚苯乙烯泡沫板，厚度为0.1m，导热系数为0.04W/(m・K)，这些参数通过查阅保温材料的产品说明书以及相关标准确定。通风口设置方面，温室顶部设有屋脊通风口，通风口长度与温室长度相同，为60m，宽度为0.5m；前屋面底部设有底部通风口，通风口长度为10m，宽度为0.3m。通风口的面积和位置对自然通风效果有着重要影响，因此在测量和记录时确保数据的准确性和完整性。3.1.2模型简化与假设在利用CFD软件进行建模的过程中，为了提高计算效率并简化模型，对温室结构和物理过程进行了一系列合理的简化处理和假设。考虑到温室内的栽培架、灌溉管道等设备虽然在实际中存在，但它们对自然通风和小气候参数分布的影响相对较小，且详细模拟这些设备会大大增加模型的复杂性和计算量。因此，在模型中忽略了这些设备的细节，将温室内部视为一个相对空旷的空间，主要关注空气的流动和热量、质量的传递过程。假设温室内的空气为理想气体，遵循理想气体状态方程pV=nRT，其中p为压强，V为体积，n为物质的量，R为理想气体常数，T为温度。这一假设简化了对空气物理性质的描述，忽略了空气的压缩性和粘性等复杂因素对计算结果的影响。在自然通风过程中，空气的流动速度相对较低，这种假设在一定程度上是合理的，能够满足本研究对自然通风效果和小气候参数分布的分析需求。假设温室墙体和覆盖材料为均质、各向同性的介质，且其物理性质（如导热系数、比热容等）在模拟过程中保持不变。实际的温室墙体和覆盖材料可能存在一定的不均匀性和温度依赖性，但在本研究中为了简化计算，忽略了这些因素。假设温室的密封性良好，忽略了温室缝隙和孔洞等可能导致的空气泄漏现象。虽然在实际情况中，温室存在一定的空气泄漏，但在本研究中，主要关注自然通风口引起的空气流动和小气候变化，因此对这些次要因素进行了简化处理。在模拟过程中，不考虑作物的生长过程和生理活动对温室内环境的动态影响。作物在生长过程中会进行光合作用、呼吸作用和蒸腾作用等，这些生理活动会消耗或释放热量、水分和二氧化碳等物质，从而影响温室内的小气候环境。然而，为了简化模型，本研究将作物视为固定的障碍物，仅考虑其对气流的阻挡作用，而不考虑其生理活动对环境的影响。3.1.3几何模型创建本研究利用三维建模软件SolidWorks进行日光温室几何模型的创建，以确保模型的准确性和完整性，为后续的CFD模拟分析提供可靠的基础。首先，打开SolidWorks软件，在新建的零件文件中，选择合适的坐标系，以保证模型在空间中的位置和方向准确。利用软件的拉伸、旋转等基本建模工具，根据测量得到的温室结构参数，逐步构建温室的各个部分。使用拉伸工具创建温室的墙体，根据测量的墙体厚度、高度和长度参数，绘制墙体的截面轮廓，然后沿相应方向拉伸得到墙体的三维模型。对于后墙，在草图绘制界面中绘制一个长60m、高3.5m、厚0.5m的矩形轮廓，通过拉伸操作生成后墙模型；山墙则根据其厚度0.4m和相应的高度、宽度参数，同样利用拉伸工具创建。接着，创建温室的屋面结构。屋面包括前屋面和后屋面，前屋面为弧形结构，利用样条曲线工具绘制屋面的弧形轮廓，然后通过旋转工具将轮廓绕温室的长轴旋转，生成前屋面的三维模型。后屋面为倾斜平面，根据测量的后屋面坡度和尺寸参数，在草图中绘制后屋面的平面轮廓，再通过拉伸操作创建后屋面模型。在创建过程中，确保各个部分的尺寸和形状与实际测量值一致，以提高模型的准确性。完成温室主体结构的创建后，添加通风口模型。对于屋脊通风口，在温室屋脊位置绘制一个长60m、宽0.5m的矩形通风口轮廓，通过拉伸切除操作在屋面模型上创建屋脊通风口；底部通风口则在前屋面底部，按照测量的长度10m和宽度0.3m，同样采用拉伸切除的方式创建通风口模型。在建模过程中，对模型的细节进行了仔细处理，确保模型的表面光滑、连接部位准确无误。对模型进行了多次检查和修正，以避免出现几何错误和不连续的情况。在完成日光温室几何模型的创建后，将模型保存为STL格式文件，以便后续导入到CFD软件Fluent中进行网格划分和模拟计算。3.2网格划分3.2.1网格划分方法选择在CFD模拟中，网格划分是至关重要的环节，其质量直接影响模拟结果的准确性和计算效率。常见的网格划分方法包括结构化网格和非结构化网格，两者各有优缺点。结构化网格是指网格区域内所有的内部点都具有相同的毗邻单元，其网格节点排列规则，网格线之间相互平行或正交。这种网格的优点显著，它能够便捷地实现区域的边界拟合，对于流体和表面应力集中等方面的计算表现出色。在处理简单几何形状的模型时，结构化网格生成速度快，质量高，数据结构简单，便于存储和计算。对一些规则形状的物体进行流体力学模拟时，结构化网格可以快速生成高质量的网格，大大提高计算效率。然而，结构化网格的适用范围相对较窄，仅适用于形状规则的图形。当面对复杂几何形状，如具有不规则边界或内部结构复杂的日光温室模型时，结构化网格的划分难度较大，可能需要对模型进行过多的简化和分割，导致模型失真，影响模拟结果的准确性。非结构化网格则是指网格区域内的内部点不具有相同的毗邻单元，网格的形状和大小可以根据模型的几何形状和计算要求灵活调整。非结构化网格的优点在于其适应性强，能够很好地处理复杂几何形状的模型。在划分日光温室模型的网格时，非结构化网格可以根据温室的复杂结构，如弯曲的屋面、不规则的通风口等，自动生成贴合模型的网格，无需对模型进行过多的简化。非结构化网格在生成过程中采用一定的准则进行优化判断，能生成高质量的网格，并且很容易控制网格大小和节点密度，采用随机的数据结构有利于进行网格自适应。一旦在边界指定网格的分布，在边界之间可以自动生成网格无需分块或者用户的干预，而且不需要在子域之间传递信息。不过，非结构化网格也存在一些缺点，它不能很好地处理粘性问题，在附面层内只采用三角形或四面体网络，其网格数量将极其巨大。对于相同的物理空间，网格填充效率不高，在满足同样流场计算条件的情况下，它产生的网格数量要比结构网格大得多，这会导致计算量增加，计算时间延长。考虑到日光温室模型的结构较为复杂，包含弯曲的屋面、不同形状和位置的通风口等不规则几何特征，若采用结构化网格划分，不仅划分过程困难重重，而且可能无法准确地描述模型的几何形状，导致模拟结果出现偏差。因此，本研究选择非结构化网格划分方法，以更好地适应日光温室模型的复杂几何形状，提高网格划分的质量和效率，为后续的CFD模拟提供可靠的网格基础。3.2.2网格质量控制为了确保模拟结果的准确性，在进行非结构化网格划分后，需要对网格质量进行严格控制。通过一系列的方法和措施，提高网格的质量，减少网格相关误差对模拟结果的影响。网格加密是提高网格质量的重要手段之一。在日光温室模型中，一些关键区域，如通风口附近、作物冠层区域以及墙体边界层等，气流和温度变化较为剧烈，对这些区域进行网格加密，可以更精确地捕捉物理量的变化。在通风口附近，气流速度和压力梯度较大，加密网格能够更准确地模拟空气的进出和流动情况；在作物冠层区域，由于作物对气流的阻挡和扰动作用，加密网格有助于更好地描述气流在作物间的流动特性。通过局部加密网格，可以在不显著增加整体计算量的前提下，提高关键区域的模拟精度。平滑处理也是改善网格质量的常用方法。在网格划分过程中，可能会出现一些形状不规则、质量较差的网格单元，这些单元会影响计算的稳定性和准确性。通过平滑处理，可以调整网格节点的位置，使网格单元的形状更加规则，边长更加均匀，从而提高网格质量。常用的平滑处理算法包括拉普拉斯平滑算法、算术平均平滑算法等。拉普拉斯平滑算法通过迭代调整节点位置，使节点周围的网格单元更加均匀；算术平均平滑算法则是将节点的位置调整为其相邻节点位置的算术平均值，以达到平滑网格的目的。在进行网格质量控制时，还需要关注网格的纵横比、雅克比行列式等指标。纵横比是指网格单元最长边与最短边的比值，较小的纵横比表示网格单元形状更接近正方形或正四面体，有利于提高计算精度。雅克比行列式用于衡量网格单元的扭曲程度，其值越接近1，表示网格单元的扭曲程度越小，质量越好。通过检查和调整这些指标，可以确保网格质量满足模拟要求。在实际操作中，利用Fluent软件自带的网格质量检查工具，对划分好的网格进行全面检查，对于不满足质量要求的网格单元，进行针对性的调整和优化，直至网格质量达到可接受的范围。通过这些网格质量控制措施，可以有效提高网格的质量，为准确模拟自然通风对日光温室小气候的影响奠定坚实的基础。3.2.3网格无关性验证网格无关性验证是CFD模拟中确保模拟结果准确性和可靠性的关键步骤。其目的是通过改变网格数量进行模拟计算，对比分析不同网格数量下的模拟结果，确定模拟结果与网格数量无关的最优网格方案，以避免因网格数量不足或过多而导致的模拟误差。在本研究中，为了进行网格无关性验证，首先基于已划分好的非结构化网格，通过调整网格生成参数，如网格尺寸、增长率等，生成了一系列不同网格数量的网格模型。分别生成了包含100万、200万、300万、400万和500万个网格单元的模型。针对每个网格模型，在相同的边界条件和模拟参数下，进行自然通风对日光温室小气候影响的CFD模拟计算。模拟计算过程中，重点关注温室内气流速度、温度和湿度等关键小气候参数在不同网格数量下的变化情况。模拟计算完成后，对不同网格数量下的模拟结果进行详细的对比分析。以温室内某一特征点的气流速度为例，绘制了该点气流速度随网格数量变化的曲线。当网格数量从100万增加到200万时，该点气流速度的模拟值变化较为明显；继续增加网格数量，从200万增加到300万时，气流速度模拟值的变化逐渐减小；当网格数量达到400万和500万时，气流速度模拟值基本保持稳定，变化幅度小于1%。对于温室内的温度和湿度分布，也进行了类似的对比分析，发现随着网格数量的增加，温度和湿度的模拟结果逐渐趋于稳定。综合考虑模拟结果的准确性和计算效率，当网格数量达到300万时，模拟结果已经基本不受网格数量的影响，且计算时间和计算资源消耗在可接受范围内。因此，确定包含300万个网格单元的网格方案为最优网格方案，在后续的模拟研究中采用该网格方案进行计算，以确保模拟结果的准确性和可靠性。通过网格无关性验证，不仅保证了本研究模拟结果的有效性，也为类似的CFD模拟研究提供了网格划分和验证的参考方法。3.3边界条件设置3.3.1进口边界条件进口边界条件的设置需紧密依据实际的自然通风情况。对于温室进风口，风速是关键参数之一。在模拟过程中，通过查阅当地气象资料，获取不同季节和时段的平均风速数据。考虑到自然通风的复杂性，将风速设置为一个范围值，以更真实地反映实际情况。在春季，当地平均风速范围为2-4m/s，因此在模拟时，将进风口风速设定在这个范围内，具体取值根据不同的模拟工况进行调整。温度方面，同样参考当地气象数据以及温室内外温度的实际测量值。室外空气温度随季节和时间变化明显，在夏季，室外平均气温较高，可达30-35℃；而在冬季，室外平均气温较低，一般在5-10℃。在设置进口温度边界条件时，根据不同季节和模拟时间，将进风口空气温度设定在相应的范围内。湿度也是重要的边界条件参数。通过湿度传感器对温室内外空气湿度进行长期监测，得到不同季节和天气条件下的湿度数据。室外空气相对湿度在不同季节有所差异，夏季相对湿度较高，可达70%-80%；冬季相对湿度较低，一般在40%-50%。在模拟中，根据实际情况，将进风口空气的相对湿度设定在对应的范围内。在Fluent软件中，将进风口定义为速度入口边界条件（Velocity-Inlet）。在设置速度入口边界条件时，除了输入风速、温度和湿度等参数外，还需要指定湍流参数。选择合适的湍流模型是准确模拟自然通风的关键，本研究采用标准k-ε湍流模型，该模型在处理自然通风等复杂湍流问题时具有较好的准确性和稳定性。对于湍流参数，根据风速和进风口尺寸，计算出湍流强度和水力直径。湍流强度I可通过公式I=0.16(Re_D)^{-1/8}计算，其中Re_D为基于水力直径的雷诺数；水力直径D_h可根据进风口的几何形状和尺寸计算得出。将计算得到的湍流强度和水力直径输入到Fluent软件的速度入口边界条件设置中，以确保模拟的准确性。3.3.2出口边界条件对于温室出风口，采用自由出流边界条件（Outflow），这是因为自由出流边界条件能够较好地模拟空气在不受外界约束的情况下自由流出的情况，符合温室出风口的实际物理过程。自由出流边界条件的设置相对简单，在Fluent软件中，只需将出风口区域定义为Outflow边界类型即可。在这种边界条件下，Fluent软件会自动根据计算域内的流动情况，确定出口处的速度、压力等物理量，无需用户手动输入具体数值。这是因为自由出流边界条件假设出口处的流动充分发展，即出口处的物理量仅取决于上游的流动情况，而与下游的外部环境无关。在实际的自然通风过程中，温室出风口的空气在排出室外时，受到的外部阻力相对较小，可以近似认为是自由出流。通过采用自由出流边界条件，能够有效地确保空气能够顺利地从温室出风口流出，避免因边界条件设置不当而导致的模拟误差。在设置自由出流边界条件时，还需要注意与其他边界条件和模拟参数的协调。确保计算域的完整性和连续性，避免出现边界不匹配或物理量突变等问题，以保证模拟结果的准确性和可靠性。3.3.3壁面边界条件温室墙体、屋面等壁面的边界条件设置对于准确模拟自然通风对小气候的影响至关重要。在本研究中，采用无滑移边界条件（No-SlipBoundaryCondition）来描述壁面与空气之间的相互作用。无滑移边界条件假设壁面与流体之间没有相对滑动，即壁面上的流体速度为零。这一假设基于实际物理现象，在壁面附近，由于流体分子与壁面之间的粘性力作用，流体的速度会逐渐减小，直至在壁面上速度降为零。在Fluent软件中，设置无滑移边界条件时，只需在壁面区域选择相应的边界类型，并确保速度分量在壁面法向和切向均为零。除了考虑壁面的无滑移特性，还需要考虑壁面的热传导和对流换热。在热传导方面，根据温室墙体和屋面材料的物理性质，确定其导热系数。墙体采用砖石材料，其导热系数为0.8W/(m・K)；屋面覆盖材料为塑料薄膜，导热系数为0.1W/(m・K)。在Fluent软件的壁面边界条件设置中，输入相应的导热系数，以准确模拟热量通过壁面的传导过程。对于壁面的对流换热，考虑自然对流和强制对流两种情况。在自然通风条件下，自然对流起着重要作用。根据壁面与空气之间的温度差和壁面的几何形状，利用自然对流换热关联式计算壁面的对流换热系数。对于垂直壁面，自然对流换热系数可通过努塞尔数关联式Nu=C(GrPr)^{n}计算，其中Nu为努塞尔数，Gr为格拉晓夫数，Pr为普朗特数，C和n为与壁面形状和流动状态相关的常数。在Fluent软件中，将计算得到的对流换热系数输入到壁面边界条件中，以考虑壁面与空气之间的对流换热过程。在存在强制通风或其他因素导致的强制对流情况下，也需要综合考虑强制对流的影响，通过相应的计算方法确定对流换热系数，并在边界条件中进行合理设置。3.4物理模型选择3.4.1湍流模型在CFD模拟中，选择合适的湍流模型对于准确预测自然通风条件下日光温室内的气流特性至关重要。常见的湍流模型包括k-ε模型、k-ω模型等，它们各自具有独特的特点和适用范围。k-ε模型是一种基于涡粘性假设的双方程湍流模型，它通过求解湍动能k和湍动能耗散率ε的输运方程来封闭雷诺应力项。标准k-ε模型在工程领域应用广泛，具有计算效率高、对高雷诺数湍流流动模拟效果较好的优点。在处理简单几何形状和均匀流场的问题时，标准k-ε模型能够快速准确地给出结果。然而，该模型在处理复杂流动，如强逆压梯度、分离流和旋转流等情况时，存在一定的局限性。在模拟日光温室自然通风时，若通风口附近存在较强的气流分离现象，标准k-ε模型可能无法准确捕捉气流的变化。为了改进标准k-ε模型的不足，出现了RNGk-ε模型和Realizablek-ε模型。RNGk-ε模型通过重整化群理论对湍流进行修正，在处理高应变率和流线弯曲较大的流动时表现更好；Realizablek-ε模型则引入了新的湍动能耗散率方程和涡粘性公式，能够更准确地模拟具有复杂流动特征的问题。k-ω模型也是一种双方程湍流模型，它求解湍动能k和比耗散率ω的输运方程。该模型对近壁面流动的模拟具有较高的精度，在处理低雷诺数流动和边界层流动时表现出色。SSTk-ω模型是k-ω模型的一种改进形式，它结合了k-ε模型和k-ω模型的优点，在近壁面区域采用k-ω模型，在远离壁面区域采用k-ε模型，通过混合函数实现两者的平滑过渡。这种模型能够更好地处理复杂的流动问题，如分离流、激波与边界层相互作用等。考虑到日光温室自然通风过程中，气流在通风口、温室内部以及作物冠层等区域存在复杂的流动现象，包括气流的加速、减速、分离和再附着等。通风口附近气流速度变化剧烈，容易出现分离流；作物冠层对气流有阻挡和扰动作用，导致气流在冠层内形成复杂的流动结构。因此，需要选择能够准确模拟这些复杂流动的湍流模型。综合比较各湍流模型的特点和适用范围，本研究选择SSTk-ω模型作为模拟日光温室自然通风的湍流模型。该模型能够较好地处理近壁面流动和复杂流动现象，对于日光温室自然通风过程中气流在不同区域的变化能够进行准确模拟，从而为后续分析自然通风对日光温室小气候的影响提供可靠的基础。3.4.2辐射模型日光温室内的辐射传热过程是一个复杂的物理过程，涉及太阳辐射的吸收、反射和透射，以及温室内各物体之间的长波辐射交换。准确模拟辐射传热对于研究日光温室内的温度分布和能量平衡至关重要。离散坐标（DO）辐射模型是一种常用的辐射模型，它基于辐射传递方程（RTE），将空间离散为有限个立体角，通过求解每个立体角方向上的辐射强度来计算辐射热流。该模型能够准确地模拟辐射在参与性介质（如空气、作物冠层等）中的传播和散射过程，适用于处理复杂几何形状和多介质的辐射问题。在日光温室中，DO辐射模型可以考虑温室覆盖材料对太阳辐射的吸收、反射和透射特性，以及温室内空气、作物和墙体等物体之间的长波辐射交换。通过准确模拟辐射传热，能够更真实地反映温室内的能量分布，进而准确预测温度场的变化。与其他辐射模型相比，如P-1模型、Rosseland模型等，DO辐射模型具有较高的精度。P-1模型是一种简化的辐射模型，它将辐射强度近似为线性分布，在处理简单几何形状和均匀介质的辐射问题时具有一定的优势，但在处理复杂的日光温室环境时，其精度相对较低。Rosseland模型主要适用于高温、光学厚度较大的介质，对于日光温室这种相对低温、光学厚度较小的环境，其适用性有限。考虑到日光温室内辐射传热的复杂性，包括太阳辐射在不同方向上的传播、温室内物体的非均匀分布以及长波辐射的多次反射和散射等。本研究选择离散坐标（DO）辐射模型来模拟日光温室内的辐射传热过程。该模型能够充分考虑温室内的各种辐射现象，准确计算辐射热流，为准确模拟自然通风条件下日光温室内的温度场提供有力支持。3.4.3传热模型温室内的传热过程涉及空气与壁面、作物之间的热量交换，是一个复杂的多物理场耦合过程。选择合适的传热模型对于准确模拟温室内的温度场至关重要。能量方程是描述传热过程的基本方程之一，它基于能量守恒原理，考虑了流体的内能、动能以及由于热传导、对流和辐射等方式引起的能量传递。在CFD模拟中，通过求解能量方程，可以得到温室内各点的温度分布。在日光温室中，能量方程能够考虑空气的对流换热、与壁面的导热以及辐射传热等多种传热方式。空气与壁面之间的导热通过壁面的导热系数和温度梯度来计算；辐射传热则通过前面选择的离散坐标（DO）辐射模型进行计算，并将其结果耦合到能量方程中。在考虑空气与作物之间的传热时，需要进一步细化能量方程。作物冠层具有一定的热容量和导热性能，会与周围空气进行热量交换。可以通过建立作物冠层的传热模型，将作物冠层对热量的吸收、储存和释放过程考虑到能量方程中。一种常用的方法是将作物冠层视为多孔介质，利用多孔介质传热模型来描述空气与作物之间的传热过程。在多孔介质模型中，定义作物冠层的孔隙率、渗透率以及传热系数等参数，通过这些参数来反映作物冠层对空气流动和传热的影响。通过这种方式，能够更准确地模拟温室内空气与作物之间的传热过程，从而提高温度场模拟的准确性。本研究采用基于能量方程的传热模型，并结合多孔介质传热模型来考虑空气与作物之间的传热，以实现对日光温室温度场的准确模拟。这种传热模型能够全面考虑温室内的各种传热因素，为深入研究自然通风对日光温室小气候的影响提供了可靠的工具。3.5模拟结果验证3.5.1实验设计与数据采集为了验证CFD模拟结果的准确性和可靠性，在实际的日光温室中开展了自然通风实验。选择了与CFD模拟中相同结构参数的日光温室，该温室位于[具体地理位置]，主要用于蔬菜种植，具备良好的代表性。在温室内布置了多个温湿度传感器和风速仪，以全面采集不同位置和时间的小气候数据。在温室的纵向和横向方向上，分别均匀设置了5个和3个测量点，共计15个测量点，以确保能够准确捕捉温室内小气候参数的空间变化。测量点的高度分别设置为0.5m（作物冠层高度）、1.5m（温室中部高度）和2.5m（接近屋面高度），以获取不同高度处的小气候信息。温湿度传感器选用了高精度的SHT30型传感器，该传感器具有测量精度高、响应速度快等优点，温度测量精度可达±0.3℃，相对湿度测量精度可达±2%RH。风速仪采用了热线式风速仪，型号为Testo405i，其测量范围为0-30m/s，测量精度为±(0.05m/s+3%测量值)，能够满足本实验对风速测量的需求。实验时间选择在[具体实验时间段]，该时间段内天气晴朗，风力稳定，便于获取稳定的实验数据。在实验过程中，记录了不同时刻各个测量点的温度、湿度和气流速度数据。每隔10分钟采集一次数据，每次采集持续5分钟，以确保数据的稳定性和可靠性。同时，记录了实验期间的室外气象数据，包括温度、湿度、风速和风向等，这些数据通过位于温室附近的气象站获取。通过对这些数据的采集和分析，为后续与CFD模拟结果的对比提供了丰富的实验依据。3.5.2模拟结果与实验数据对比分析将CFD模拟结果与实验数据进行了详细的对比分析，通过绘制温度、湿度、气流速度等参数的对比曲线，直观地展示模拟结果与实验数据的差异，从而深入分析模拟结果的准确性和可靠性。以温度为例，选取了温室内具有代表性的测量点，绘制了该点在不同时刻的模拟温度与实验温度对比曲线。在上午10点时，实验测得该点温度为28.5℃，而CFD模拟结果为28.8℃，两者相差0.3℃，相对误差约为1.05%；在中午12点，实验温度为32.0℃，模拟温度为32.5℃，相差0.5℃，相对误差为1.56%。从整体温度变化趋势来看，模拟曲线与实验曲线基本吻合，能够较好地反映温室内温度随时间的变化规律。对于湿度，同样选取典型测量点进行对比。在下午2点，实验测得该点相对湿度为70%，模拟结果为68%，相差2%；在下午4点，实验相对湿度为65%，模拟结果为66%，相差1%。湿度的模拟结果与实验数据也较为接近，能够准确地反映温室内湿度的变化情况。在气流速度方面，选取了通风口附近和温室中部的测量点进行对比。通风口附近，由于气流速度变化较为剧烈，模拟结果与实验数据存在一定差异，但整体趋势一致。在某一时刻，实验测得通风口附近气流速度为2.5m/s，模拟结果为2.3m/s，相差0.2m/s；在温室中部，气流速度相对稳定，模拟结果与实验数据更为接近，某时刻实验测得气流速度为0.8m/s，模拟结果为0.78m/s，相差0.02m/s。通过对多个测量点不同时刻的温度、湿度和气流速度的对比分析，可以看出CFD模拟结果与实验数据在整体趋势上基本一致，大部分测量点的模拟值与实验值的相对误差在可接受范围内。这表明本研究建立的CFD模型能够较为准确地模拟自然通风条件下日光温室内的小气候参数分布和变化情况，为后续研究提供了可靠的基础。然而，在一些局部区域，如通风口附近等气流复杂的区域，模拟结果与实验数据仍存在一定偏差，需要进一步分析原因并对模型进行优化。3.5.3模型验证与修正根据模拟结果与实验数据的对比分析，对CFD模型进行了验证和修正，以进一步提高模型的模拟精度。通过对比发现，在通风口附近和作物冠层区域，模拟结果与实验数据的偏差相对较大。通风口附近，由于气流受到通风口形状、尺寸以及周围环境的影响，流动情况较为复杂，导致模拟结果与实验存在一定差异。在作物冠层区域，作物对气流的阻挡和扰动作用较为明显，而模型在考虑作物冠层对气流的影响时，可能存在一定的简化，导致模拟精度受到影响。针对这些问题，对CFD模型进行了以下修正。在边界条件设置方面，进一步优化了通风口的边界条件，考虑了通风口的局部阻力系数和流量系数，以更准确地描述通风口处的气流进出情况。通过查阅相关文献和实验数据，确定了通风口的局部阻力系数和流量系数，并将其输入到Fluent软件中，对通风口边界条件进行了重新设置。在物理模型参数优化方面，对作物冠层模型进行了改进。考虑了作物的高度、种植密度以及叶片形态等因素对气流的影响，通过建立更准确的作物冠层阻力模型，提高了模型对作物冠层区域气流模拟的准确性。在模型中，根据实际作物的参数，调整了作物冠层的孔隙率、渗透率以及传热系数等参数，使模型能够更真实地反映作物冠层对气流和热量传递的影响。经过修正后，再次进行CFD模拟，并将模拟结果与实验数据进行对比。结果表明，修正后的模型模拟精度得到了显著提高，在通风口附近和作物冠层区域，模拟结果与实验数据的偏差明显减小。在通风口附近，某时刻模拟气流速度与实验值的相对误差从修正前的8%降低到了3%；在作物冠层区域，某测量点的模拟温度与实验值的相对误差从修正前的5%降低到了2%。这表明通过对CFD模型的验证和修正，有效提高了模型的模拟精度，使其能够更准确地模拟自然通风对日光温室小气候的影响。四、自然通风对日光温室小气候影响的模拟结果分析4.1温度场分布特征4.1.1不同通风工况下的温度分布在不同通风工况下，日光温室内温度分布呈现出明显的差异。当通风口大小发生变化时，对温室内温度分布有着显著影响。在仅开启屋脊通风口且通风口宽度为0.3m时，温室内温度分布存在较大的不均匀性。靠近屋脊通风口的区域，由于热空气能够迅速排出，温度相对较低，约为28℃；而在温室底部靠近后墙的区域，空气流通不畅，温度较高，可达32℃。这是因为通风口较小，通风量不足，导致室内空气无法充分混合，热空气难以扩散到通风口排出。当通风口宽度增大到0.5m时，温室内温度分布的均匀性得到明显改善。此时，通风量增加，空气流通速度加快，热空气能够更迅速地从通风口排出，室内各区域温度差异减小。靠近屋脊通风口区域的温度降至27℃左右，温室底部靠近后墙区域的温度也降低到30℃，整个温室内的温度分布更加均匀，有利于作物的生长。不同通风时间也会导致温室内温度分布的变化。在上午9点开启通风口，通风1小时后，温室内温度开始逐渐下降。此时，温室前部靠近通风口的区域温度下降较为明显，从初始的30℃降至28℃；而温室后部由于空气流通相对较慢，温度下降幅度较小，仍保持在30.5℃左右。随着通风时间延长至3小时，温室内各区域温度均有显著下降，且分布更加均匀。温室前部温度降至26℃，后部温度降至27℃，整个温室内的温度差缩小到1℃以内，表明通风时间的延长有助于改善温室内温度分布的均匀性。外界风速的变化同样对温室内温度分布产生影响。当外界风速为2m/s时，风压通风作用较弱，温室内主要以热压通风为主。此时，温室内温度分布呈现出上低下高的趋势，靠近顶部通风口的区域温度较低，为27℃，而底部温度较高，为31℃。当外界风速增大到4m/s时，风压通风作用增强，与热压通风共同作用，使温室内空气流动更加剧烈。此时，温室内温度分布更加均匀，各区域温度均在28-29℃之间，温度差明显减小。这说明外界风速的增加能够增强通风效果，改善温室内温度分布。4.1.2温度随时间的变化规律通过模拟不同时间段内温室内温度的变化，深入分析自然通风对温室温度日变化和季节变化的影响，找出温度变化的峰值和谷值。在典型的晴天，温室内温度日变化呈现出明显的规律。在清晨6点，温室内温度为20℃，此时太阳辐射较弱，温室内热量主要来自于夜间的蓄热。随着太阳辐射的增强，温室内温度逐渐升高，在上午10点左右，温度达到26℃。此时，自然通风开始发挥作用，通过热压通风和风压通风，室内热空气排出，冷空气进入，温度上升速度减缓。到中午12点，太阳辐射最强，温室内温度达到峰值32℃。之后，随着通风的持续进行，以及太阳辐射强度的逐渐减弱，温室内温度开始下降。在下午4点，温度降至28℃，此时通风仍在继续，室内外空气不断交换，热量持续散失。到傍晚6点，太阳辐射消失，温室内温度继续下降，降至24℃。夜间，由于没有太阳辐射，温室内热量逐渐散失，温度持续降低，在清晨6点达到谷值20℃。在不同季节，温室内温度变化也有所不同。在春季，室外温度相对较低，但随着太阳辐射的增强，温室内温度上升较快。在晴朗的春季，温室内最高温度可达30℃左右，最低温度在18℃左右。自然通风能够有效地调节温室内温度，避免温度过高对作物生长造成影响。在夏季，室外温度较高，温室内温度上升更为迅速。在没有有效通风的情况下，温室内最高温度可达35℃以上，对作物生长极为不利。通过自然通风，温室内最高温度可控制在32℃左右，能够满足大多数作物的生长需求。在秋季，室外温度逐渐降低，温室内温度变化相对较为平缓。自然通风的作用主要是调节室内湿度和补充二氧化碳，对温度的调节作用相对较小。在冬季，室外温度较低，温室内需要依靠保温措施来保持温度。自然通风在冬季的主要作用是在白天适当通风，排出室内的湿气和污浊空气，同时补充二氧化碳，但通风时间和通风量需要严格控制，以避免室内温度过低对作物造成冻害。4.1.3温度分布不均匀性分析利用温度梯度等指标，对温室内温度分布的不均匀性进行分析，探讨不均匀性产生的原因，如通风方式、温室结构等。通过计算温室内不同位置的温度梯度，可以直观地反映温度分布的不均匀程度。在温室的纵向方向上，从温室前部到后部，温度梯度在不同通风工况下有所不同。在通风不良的情况下，如通风口面积较小或通风时间较短时，温度梯度较大，可达1-2℃/m。这是因为通风不足导致热空气在温室后部积聚，无法及时排出，从而使后部温度明显高于前部。而在通风良好的情况下，温度梯度较小，一般在0.2-0.5℃/m之间。此时，通风量充足，空气能够充分混合，热空气能够均匀地分布在温室内，减少了温度差异。在温室的横向方向上，温度梯度也存在一定的变化。靠近通风口的一侧，温度相对较低，而远离通风口的一侧，温度相对较高。在仅开启顶部通风口时，温室靠近顶部通风口一侧的温度比远离通风口一侧低1-3℃，这是由于通风口附近的空气流动较快，热量能够迅速带走，而远离通风口的区域空气流动较慢，热量积聚，导致温度升高。通风方式是导致温度分布不均匀的重要原因之一。顶通风和底通风相结合的方式，能够形成较好的空气对流，使温室内温度分布更加均匀。当顶部通风口和底部通风口同时开启时，热空气从顶部通风口排出，冷空气从底部通风口进入，形成上下对流，能够有效地混合室内空气，降低温度梯度。而仅采用顶通风或底通风时，空气流动路径相对单一，容易出现通风死角，导致温度分布不均匀。温室结构对温度分布也有显著影响。温室的跨度和高度会影响空气的流通和热量的分布。跨度较大的温室，空气流通相对困难，容易出现温度不均匀的情况。跨度为10m的温室比跨度为8m的温室更容易出现温度梯度较大的问题，因为跨度大使得空气在室内的流动距离增加，热量难以均匀扩散。温室的墙体和屋面材料的保温性能也会影响温度分布。保温性能好的材料能够减少热量的散失，使温室内温度更加稳定，但如果通风设计不合理，也容易导致温度分布不均匀。4.2湿度场分布特征4.2.1湿度分布情况在自然通风条件下，日光温室内湿度场呈现出复杂的分布情况。从垂直方向来看，靠近地面区域的湿度相对较高，随着高度的增加，湿度逐渐降低。在温室底部靠近地面0.2m高度处，由于土壤水分的蒸发以及作物蒸腾作用产生的水汽聚集，相对湿度可达80%-85%。而在靠近屋面2.5m高度处，相对湿度降低至65%-70%。这是因为热空气上升过程中，携带的水汽逐渐扩散，同时与通风口进入的相对干燥的室外空气混合，导致湿度下降。在水平方向上，湿度