多定日镜风场数值模拟：理论、方法与应用的深度剖析

多定日镜风场数值模拟：理论、方法与应用的深度剖析一、绪论1.1研究背景与意义随着全球对清洁能源需求的不断增长，太阳能作为一种丰富、清洁且可再生的能源，受到了广泛关注。塔式太阳能热发电技术作为太阳能利用的重要方式之一，具有发电效率高、可储能等优势，在未来能源结构中占据着重要地位。定日镜是塔式太阳能热发电系统的关键组成部分，其作用是将太阳光线反射并汇聚到位于高塔顶部的接收器上，实现太阳能到热能的转换。定日镜通常数量众多，组成规模庞大的定日镜场，如美国SolarTwo电站拥有1926台定日镜，西班牙PS10电站则包含624台大型定日镜。这些定日镜在运行过程中，长期暴露在自然环境中，不可避免地受到风载荷的作用。风载荷是定日镜设计和运行中需要考虑的重要因素之一，其对定日镜的影响是多方面的。强风可能导致定日镜结构变形、损坏，甚至倒塌，影响整个发电系统的正常运行和安全性。据相关研究表明，在一些风力较大的地区，定日镜因风致破坏的事故时有发生，造成了严重的经济损失和发电中断。风还会引起定日镜的振动，影响其跟踪精度和聚光效率，进而降低发电系统的整体性能。由于风的随机性和复杂性，准确评估风载荷对定日镜的影响具有一定难度。多定日镜风场数值模拟对于塔式太阳能热发电技术的发展具有重要意义。通过数值模拟，可以深入了解风场特性以及风载荷在定日镜上的分布规律，为定日镜的结构设计和优化提供科学依据，从而降低结构成本。精确的数值模拟能够提高定日镜的抗风能力和稳定性，保障发电系统的安全可靠运行，减少因风灾导致的停机时间和维修成本，提高发电效率和经济效益。数值模拟还可以为定日镜场的布局设计提供参考，优化定日镜的排列方式和间距，减小风干扰，提高整个镜场的性能。1.2定日镜结构简介定日镜主要由平面反射镜、跟踪机构、镜架和控制系统等部分组成。平面反射镜是定日镜的核心部件，其作用是反射太阳光线。目前，常用的反射镜材料为玻璃，表面镀银以提高反射率，反射率可达97%。为了使阳光经定日镜反射后能更好地聚集到集热器内，部分定日镜采用了镜表面具有微小弧度（16′）的平凹面镜。跟踪机构负责调整定日镜的角度，使其能够实时跟踪太阳的位置，确保太阳光准确反射到集热器上。跟踪机构通常包括方位角调整机构和高度角调整机构，可实现定日镜在两个方向上的转动。方位角调整机构一般由电机、减速机和齿轮等组成，通过电机驱动齿轮转动，实现定日镜在水平方向上的旋转；高度角调整机构则通过电机带动丝杆或链条，使定日镜在垂直方向上进行俯仰运动。镜架用于支撑平面反射镜和跟踪机构，要求具有足够的强度和刚度，以保证定日镜在各种工况下的稳定性。镜架一般采用金属材料制作，如钢结构或铝合金结构，其结构形式有多种，常见的有矩形独立柱支架式、三角形支架式等。控制系统是定日镜的“大脑”，它根据太阳的位置信息，控制跟踪机构的动作，实现定日镜的精确跟踪。控制系统通常采用计算机控制技术，结合传感器实时监测太阳的位置、定日镜的角度等参数，通过算法计算出跟踪机构的控制指令，确保定日镜始终对准太阳。在太阳能塔式发电系统中，大量的定日镜按照一定的布局组成定日镜场。定日镜场的布局设计需要考虑多种因素，如太阳的运行轨迹、地形条件、集热器的位置等，以确保每个定日镜都能充分接收阳光，并将反射光准确地汇聚到集热器上。常见的定日镜场布局有圆形、方形、扇形等。在圆形布局中，集热器位于圆心位置，定日镜呈环形分布，这种布局能够充分利用太阳的辐射能，且便于管理和维护；方形布局则适用于地形较为规则的场地，可提高土地利用率；扇形布局则根据太阳的主要照射方向进行设计，能在一定程度上减少定日镜之间的遮挡。1.3结构风工程研究方法综述1.3.1现场实测方法现场实测是获取真实风场数据和结构风效应的直接手段。通过在实际工程现场布置风速仪、风压传感器等设备，可以实时测量风的速度、方向、风压等参数，以及结构在风作用下的响应，如位移、加速度、应力等。现场实测能够反映真实环境中风的特性和结构的实际受力情况，数据具有较高的可靠性和真实性。在一些大型建筑或桥梁的风工程研究中，现场实测可以获取建筑物或桥梁在自然风作用下的动态响应，为结构的抗风设计提供重要依据。对于定日镜风场的研究，现场实测可以直接测量定日镜在实际风场中的受力和变形情况，验证数值模拟和理论分析的结果。现场实测也存在一些局限性。现场实测受到环境条件的限制，如地形、气候等因素会影响风场的特性，使得测量结果的代表性受到一定影响。在山区等地形复杂的地区，风的流动会受到地形的干扰，导致风速和风向的分布不均匀，增加了测量的难度和不确定性。现场实测需要投入大量的人力、物力和时间，成本较高。安装和维护测量设备需要专业人员，且测量周期较长，才能获取足够的有效数据。测量数据的处理和分析也较为复杂，需要专业的技术和软件支持。1.3.2风洞试验方法风洞试验是结构风工程研究中常用的方法之一。它通过在风洞中模拟真实的风场环境，将结构模型放置在风洞中，测量模型在风作用下的受力和响应情况。风洞试验可以精确控制风速、风向、湍流度等风场参数，能够较为准确地模拟不同工况下的风场，研究结构的风效应。在风洞试验中，首先需要根据实际结构制作缩尺模型，模型的制作要求尽可能准确地反映实际结构的几何形状、材料特性和边界条件。然后，将模型安装在风洞的试验段中，通过调节风洞的风机和相关设备，产生不同风速和风向的风场。利用压力传感器、应变片、位移传感器等测量设备，测量模型表面的风压分布、结构的应力和变形等参数。通过对测量数据的分析，可以得到结构在风作用下的受力特性和响应规律。风洞试验在研究高层建筑、桥梁、大跨度结构等的风效应方面取得了丰硕的成果，为这些结构的抗风设计提供了重要的技术支持。对于定日镜风场的研究，风洞试验可以研究单个定日镜和多定日镜之间的风干扰效应，分析不同布局和间距下定日镜的风载荷分布情况。风洞试验也存在一些缺点。风洞试验的成本较高，建设和维护风洞需要大量的资金投入，且试验过程中需要消耗一定的能源。制作高精度的结构模型和安装测量设备也需要较高的成本。风洞试验的周期较长，从模型制作、试验准备到数据采集和分析，需要花费较多的时间。风洞试验中模型与实际结构之间存在一定的相似性误差，虽然可以通过相似理论进行模型设计和参数换算，但在某些情况下，模型试验结果与实际结构的响应仍可能存在差异。1.3.3数值模拟方法数值模拟方法是随着计算机技术的发展而兴起的一种结构风工程研究方法。它基于计算流体力学（CFD）、计算结构力学（CSM）等理论，通过求解控制方程，对风场和结构在风作用下的响应进行数值计算。数值模拟方法可以模拟复杂的风场环境和结构形式，具有高效、灵活、可重复性强等优点。在数值模拟中，首先需要建立风场和结构的数学模型，将实际问题转化为数学问题。对于风场的模拟，常用的方法有基于雷诺平均N-S方程（RANS）的方法、大涡模拟（LES）方法和直接数值模拟（DNS）方法等。基于RANS方程的方法通过对N-S方程进行时间平均，引入湍流模型来封闭方程，求解得到平均流场和湍流特性，计算效率较高，适用于工程实际问题的计算；LES方法则直接模拟大尺度湍流结构，对小尺度湍流采用亚格子模型进行模拟，计算精度较高，但计算量较大；DNS方法不引入任何湍流模型，直接求解N-S方程，能够精确模拟湍流的所有尺度，但计算量巨大，目前仅适用于简单流动和低雷诺数情况。对于结构在风作用下的响应计算，常用的方法有有限元法（FEM）、有限差分法（FDM）等。有限元法将结构离散为有限个单元，通过求解单元的平衡方程，得到结构的应力、应变和位移等响应；有限差分法则将控制方程在空间和时间上进行离散化，通过差分格式求解离散后的代数方程，得到结构的响应。数值模拟方法可以方便地改变风场参数、结构形式和边界条件，进行多种工况的模拟分析，为结构的抗风设计提供全面的参考。对于多定日镜风场的研究，数值模拟可以考虑定日镜之间的相互干扰、地形的影响等复杂因素，深入分析风场特性和定日镜的风载荷分布。数值模拟还可以与现场实测和风洞试验相结合，相互验证和补充，提高研究结果的可靠性和准确性。1.4数值模拟研究风工程的成果及现状随着计算机技术的飞速发展，数值模拟在风工程研究中的应用越来越广泛，取得了众多成果。在网格生成技术方面，从早期简单的结构化网格发展到现在的非结构化网格、混合网格等，能够更好地适应复杂的几何形状和边界条件。多块分区网格技术、自适应网格技术的出现，进一步提高了网格生成的效率和质量，使得数值模拟能够更加准确地模拟复杂风场。在多定日镜风场数值模拟研究方面，国内外学者已经开展了大量工作。一些研究通过数值模拟分析了定日镜的单体风载荷特性，研究了不同形状、尺寸和安装角度的定日镜在风作用下的风压分布和受力情况。对于多定日镜之间的风干扰效应，也有不少研究通过数值模拟揭示了定日镜间距、布局方式对风载荷的影响规律。研究发现，当定日镜间距较小时，风干扰效应明显，会导致定日镜风载荷增大；合理的布局方式可以有效减小风干扰，降低定日镜的风载荷。计算风工程仍然面临一些困难和挑战。在模型简化方面，为了提高计算效率，往往需要对实际问题进行一定的简化，但简化模型可能会忽略一些重要因素，影响计算结果的准确性。在模拟复杂地形对风场的影响时，如何准确地考虑地形的粗糙度、地形起伏等因素，仍然是一个有待解决的问题。计算精度也是一个关键问题，虽然数值模拟方法不断发展，但在某些情况下，计算结果与实际情况仍存在一定偏差。湍流模型的选择和参数设置对计算结果有较大影响，不同的湍流模型在不同的流动条件下表现出不同的精度，如何选择合适的湍流模型是提高计算精度的关键之一。计算资源的限制也制约了数值模拟的发展，对于大规模的多定日镜风场模拟，需要大量的计算时间和内存，如何提高计算效率，降低计算成本，是当前研究的热点之一。1.5本文的主要工作本文旨在深入研究多定日镜风场的特性，通过数值模拟方法分析风载荷在定日镜上的分布规律，为定日镜的结构设计和优化提供理论依据。具体工作如下：理论分析：对结构风工程的相关理论进行深入研究，包括风场的基本特性、湍流模型、数值模拟方法等，为后续的数值模拟提供理论基础。模型建立：基于计算流体力学软件，建立多定日镜风场的数值模型，考虑定日镜的几何形状、布局方式、地形等因素，确保模型能够准确反映实际情况。数值模拟：运用建立的数值模型，对不同工况下的多定日镜风场进行模拟计算，分析风场特性，如风速分布、风压分布等，研究定日镜之间的风干扰效应，探讨风载荷与定日镜间距、布局方式等因素的关系。结果验证：将数值模拟结果与现场实测数据或风洞试验结果进行对比分析，验证数值模型的准确性和可靠性，对模拟结果进行修正和完善。结果分析与应用：根据数值模拟结果，总结多定日镜风场的规律，提出定日镜结构设计和布局优化的建议，为塔式太阳能热发电系统的工程实践提供参考。二、数值模拟基本理论和方法2.1流体动力学控制方程2.1.1质量守恒方程质量守恒定律是自然界的基本定律之一，在流体流动中，其体现为质量不会凭空产生或消失，只会在流场中转移。对于一个固定的控制体，单位时间内流入控制体的质量与流出控制体的质量之差，应等于控制体内质量的变化率。其数学表达式在笛卡尔坐标系下为：\frac{\partial\rho}{\partialt}+\frac{\partial(\rhou)}{\partialx}+\frac{\partial(\rhov)}{\partialy}+\frac{\partial(\rhow)}{\partialz}=0其中，\rho为流体密度，t为时间，u、v、w分别为流体在x、y、z方向上的速度分量。该方程的物理含义深刻，\frac{\partial\rho}{\partialt}表示单位时间内控制体内流体密度的变化率，反映了控制体内质量随时间的增减情况；\frac{\partial(\rhou)}{\partialx}+\frac{\partial(\rhov)}{\partialy}+\frac{\partial(\rhow)}{\partialz}则代表了单位时间内通过控制体表面的质量通量的散度，体现了流体流入和流出控制体对质量分布的影响。当流场处于定常状态时，即流场内各点的物理量不随时间变化，\frac{\partial\rho}{\partialt}=0，此时质量守恒方程简化为\nabla\cdot(\rho\vec{V})=0，其中\vec{V}=u\vec{i}+v\vec{j}+w\vec{k}为速度矢量。在不可压缩流体的情况下，流体密度\rho为常数，质量守恒方程进一步简化为\frac{\partialu}{\partialx}+\frac{\partialv}{\partialy}+\frac{\partialw}{\partialz}=0，意味着速度场的散度为零，即流入控制体某一微小体积元的流体体积等于流出该体积元的流体体积，反映了不可压缩流体在流动过程中体积的守恒性。2.1.2动量守恒方程动量守恒方程在描述流体受力和运动变化中起着核心作用，它基于牛顿第二定律，即物体的动量变化率等于作用在物体上的合外力。在流体力学中，动量守恒方程考虑了流体内部的应力以及外部施加的各种力对流体动量的影响。在笛卡尔坐标系下，忽略质量力时，粘性流体的动量守恒方程（Navier-Stokes方程）的分量形式为：\begin{cases}\rho(\frac{\partialu}{\partialt}+u\frac{\partialu}{\partialx}+v\frac{\partialu}{\partialy}+w\frac{\partialu}{\partialz})=-\frac{\partialp}{\partialx}+\mu(\frac{\partial^{2}u}{\partialx^{2}}+\frac{\partial^{2}u}{\partialy^{2}}+\frac{\partial^{2}u}{\partialz^{2}})\\\rho(\frac{\partialv}{\partialt}+u\frac{\partialv}{\partialx}+v\frac{\partialv}{\partialy}+w\frac{\partialv}{\partialz})=-\frac{\partialp}{\partialy}+\mu(\frac{\partial^{2}v}{\partialx^{2}}+\frac{\partial^{2}v}{\partialy^{2}}+\frac{\partial^{2}v}{\partialz^{2}})\\\rho(\frac{\partialw}{\partialt}+u\frac{\partialw}{\partialx}+v\frac{\partialw}{\partialy}+w\frac{\partialw}{\partialz})=-\frac{\partialp}{\partialz}+\mu(\frac{\partial^{2}w}{\partialx^{2}}+\frac{\partial^{2}w}{\partialy^{2}}+\frac{\partial^{2}w}{\partialz^{2}})\end{cases}其中，p为流体压力，\mu为动力粘度。方程左边\rho(\frac{\partialu}{\partialt}+u\frac{\partialu}{\partialx}+v\frac{\partialu}{\partialy}+w\frac{\partialu}{\partialz})表示单位体积流体的动量变化率，包含了非定常项\rho\frac{\partialu}{\partialt}和对流项\rho(u\frac{\partialu}{\partialx}+v\frac{\partialu}{\partialy}+w\frac{\partialu}{\partialz})。非定常项反映了流场内某点速度随时间的变化对动量的影响；对流项则体现了由于流体的宏观运动，使得动量在空间中发生输运和变化。方程右边-\frac{\partialp}{\partialx}+\mu(\frac{\partial^{2}u}{\partialx^{2}}+\frac{\partial^{2}u}{\partialy^{2}}+\frac{\partial^{2}u}{\partialz^{2}})表示作用在单位体积流体上的力，其中-\frac{\partialp}{\partialx}为压力梯度力，它促使流体从高压区域流向低压区域；\mu(\frac{\partial^{2}u}{\partialx^{2}}+\frac{\partial^{2}u}{\partialy^{2}}+\frac{\partial^{2}u}{\partialz^{2}})为粘性力，反映了流体内部不同流层之间由于粘性而产生的相互作用力，其作用是阻碍流体的相对运动。对于理想流体，即不考虑粘性的流体，动力粘度\mu=0，动量守恒方程简化为欧拉方程：\begin{cases}\rho(\frac{\partialu}{\partialt}+u\frac{\partialu}{\partialx}+v\frac{\partialu}{\partialy}+w\frac{\partialu}{\partialz})=-\frac{\partialp}{\partialx}\\\rho(\frac{\partialv}{\partialt}+u\frac{\partialv}{\partialx}+v\frac{\partialv}{\partialy}+w\frac{\partialv}{\partialz})=-\frac{\partialp}{\partialy}\\\rho(\frac{\partialw}{\partialt}+u\frac{\partialw}{\partialx}+v\frac{\partialw}{\partialy}+w\frac{\partialw}{\partialz})=-\frac{\partialp}{\partialz}\end{cases}2.1.3能量守恒方程能量守恒方程在考虑流体能量转化和传递时具有关键应用，它遵循热力学第一定律，即能量既不能被创造也不能被消灭，只能从一种形式转化为另一种形式。在流体流动中，能量的形式主要包括内能、动能和压力能，此外还可能涉及热传导、热辐射以及由于粘性作用产生的耗散能等。一般情况下，可压缩粘性流体的能量守恒方程在笛卡尔坐标系下的表达式为：\rhoc_{p}(\frac{\partialT}{\partialt}+u\frac{\partialT}{\partialx}+v\frac{\partialT}{\partialy}+w\frac{\partialT}{\partialz})=\lambda(\frac{\partial^{2}T}{\partialx^{2}}+\frac{\partial^{2}T}{\partialy^{2}}+\frac{\partial^{2}T}{\partialz^{2}})+\Phi+S_{h}其中，c_{p}为定压比热容，T为温度，\lambda为热导率，\Phi为粘性耗散函数，S_{h}为热源项。方程左边\rhoc_{p}(\frac{\partialT}{\partialt}+u\frac{\partialT}{\partialx}+v\frac{\partialT}{\partialy}+w\frac{\partialT}{\partialz})表示单位体积流体的内能变化率，其中\rhoc_{p}\frac{\partialT}{\partialt}为非定常项，反映了单位时间内由于温度随时间变化引起的内能改变；\rhoc_{p}(u\frac{\partialT}{\partialx}+v\frac{\partialT}{\partialy}+w\frac{\partialT}{\partialz})为对流项，体现了由于流体的宏观流动导致内能在空间中的输运。方程右边\lambda(\frac{\partial^{2}T}{\partialx^{2}}+\frac{\partial^{2}T}{\partialy^{2}}+\frac{\partial^{2}T}{\partialz^{2}})表示通过热传导传递的热量，它描述了由于温度梯度的存在，热量从高温区域向低温区域传递的现象；\Phi为粘性耗散函数，代表了由于流体粘性作用，机械能转化为热能的部分，体现了粘性力对流体能量的耗散效应；S_{h}为热源项，表示单位体积内的外部热源对流体能量的输入，例如化学反应放热、电加热等情况。在一些简化情况下，若忽略粘性耗散和热源项，且假设流体为不可压缩的理想流体，能量守恒方程可简化为：\rhoc_{p}(\frac{\partialT}{\partialt}+u\frac{\partialT}{\partialx}+v\frac{\partialT}{\partialy}+w\frac{\partialT}{\partialz})=\lambda(\frac{\partial^{2}T}{\partialx^{2}}+\frac{\partial^{2}T}{\partialy^{2}}+\frac{\partial^{2}T}{\partialz^{2}})2.2控制方程的离散有限体积法是离散控制方程的常用方法之一，其原理基于积分守恒原理，将计算区域划分为一系列不重叠的控制体积，使每个控制体积都包围一个节点。在每个控制体积上对控制方程进行积分，将偏微分方程转化为关于节点物理量的代数方程，从而实现对控制方程的离散求解。具体步骤如下：首先，对计算区域进行网格划分，生成一系列的控制体积，确定每个控制体积的边界和节点位置。然后，在每个控制体积上对质量守恒方程、动量守恒方程和能量守恒方程进行积分。以质量守恒方程为例，对控制体积V积分：\int_{V}\frac{\partial\rho}{\partialt}dV+\oint_{S}\rho\vec{V}\cdotd\vec{S}=0其中，S为控制体积的表面，\vec{V}为速度矢量，d\vec{S}为表面微元矢量。通过高斯公式将面积分转化为体积分，再利用合适的插值方法，将控制体积界面上的物理量（如密度、速度等）用节点物理量表示，从而得到离散形式的质量守恒方程。常用的离散格式有中心差分格式、一阶迎风格式、混合格式、指数格式和乘方格式等。中心差分格式通过线性插值计算控制体积界面上的物理量，具有二阶精度，在扩散占优的情况下计算精度较高，但当流动为强对流时，离散方程系数可能出现负值，导致计算的收敛性和精度变差。当Peclet数（表征对流与扩散强度之比）小于2时，中心差分格式的计算结果与精确解基本吻合；当Peclet数大于2时，数值解可能失去物理意义。一阶迎风格式在确定界面物理量时考虑了流动方向，将界面上的未知量取为上游节点的值，具有一阶截差，计算稳定，不会引起解的振荡，离散方程系数总是大于零。该格式夸大了扩散项的影响，存在数值扩散问题，使流场中的梯度被抹平。混合格式根据流体的Peclet数在中心差分格式和迎风格式之间进行切换，综合了两种格式的优点，在CFD软件中广泛应用，但只有一阶精度。指数格式是根据一维、稳态、无源的对流扩散方程的精确解建立的离散格式，能对任何Peclet数和网格点都得到精确解，但指数运算耗时，且对于二维或三维问题以及源项不为零的情况，该方案不准确。乘方格式是参考精确解所做的近似，与指数格式精度相近但更省时，与混合格式性质类似，可用作混合格式的替代格式。2.3流场数值计算的主要方法常见的流场数值计算方法包括有限差分法、有限元法和有限体积法等。有限差分法是最早应用的数值分析方法之一，它将求解域划分为差分网格，用有限个网格节点代替连续的求解域。通过泰勒级数展开，把控制方程中的导数用网格节点上的函数值的差商代替进行离散，从而建立以网格节点上的值为未知数的代数方程组。有限差分格式多样，从格式精度可分为一阶格式、二阶格式和高阶格式；从差分的空间形式可分为中心格式和逆风格式；考虑时间因子影响，还可分为显格式、隐格式、显隐交替格式等。该方法数学概念直观，表达简单，适用于有结构网格，在求解双曲型和抛物型偏微分方程时具有优势。在求解波动方程时，有限差分法能够较为准确地模拟波动的传播特性。对于复杂边界条件和椭圆型问题，有限差分法的处理相对困难。有限元法以变分原理和加权余量法为基础，将计算域划分为有限个互不重叠的单元，在每个单元内选择合适的节点作为求解函数的插值点，将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式，借助变分原理或加权余量法将微分方程离散求解。根据权函数和插值函数的不同，有限元法可分为多种计算格式。从权函数选择上，有配置法、矩量法、最小二乘法和伽辽金法等；从计算单元网格形状划分，有三角形网格、四边形网格和多边形网格等；从插值函数精度划分，又分为线性插值函数和高次插值函数等。有限元法能够适应复杂的几何形状和边界条件，通过调整权函数和插值函数的形式，可以生成不同的有限元模型，在结构力学、流体力学等领域得到广泛应用。在模拟复杂形状的飞行器绕流问题时，有限元法能够较好地处理飞行器表面的复杂边界。有限元法的计算过程相对复杂，计算量较大，对计算机内存和计算速度要求较高。有限体积法将区域划分为体积元，通过积分和平衡原理处理方程，主要关注流体流动和传热等问题。它在处理不可压缩流体和包含固态边界的问题时具有优越性，在计算流体力学中应用广泛。有限体积法离散的核心是使用有限个离散点代替原来整个连续的空间，把计算区域分成不重叠的计算网格，确定每个节点位置和节点控制体体积。在控制体上积分控制方程，以便在控制体积的界面上产生离散方程。有限体积法在离散过程中保证了物理量在每个控制体积上的守恒性，这是其重要特点之一。在模拟管道内的流体流动时，有限体积法能够准确地计算管道壁面附近的流动特性。2.4湍流及湍流的数值模拟方法2.4.1湍流及其数学描述湍流是一种高度复杂的流动状态，其特性表现为随机性和不规则性。在湍流中，流体的速度、压力等物理量在空间和时间上都呈现出强烈的脉动，流体质点的运动轨迹极为复杂，形成各种尺度的涡旋结构，这些涡旋相互作用、合并和分裂，使得湍流流动的能量在不同尺度间不断传递和耗散。用于描述湍流的数学模型主要基于雷诺平均Navier-Stokes（RANS）方程。通过对瞬时的Navier-Stokes方程进行时间平均，得到平均流场的运动方程，但此时方程中会出现雷诺应力项，它反映了脉动速度对平均流场的影响。雷诺应力项使得方程不封闭，需要引入湍流模型来建立雷诺应力与平均速度之间的关系，从而实现对方程的求解。常见的湍流模型有零方程模型（如Prandtl混合长模型）、一方程模型（如Spalart-Allmaras模型）和两方程模型（如k-ε模型、k-ω模型）等。Prandtl混合长模型假设湍流切应力与平均速度梯度成正比，通过引入混合长的概念来描述湍流的特性，该模型形式简单，但适用范围有限，仅适用于一些简单的湍流流动。Spalart-Allmaras模型是一种一方程模型，它求解一个关于湍流粘性的输运方程，在航空航天领域中对于模拟边界层流动等具有较好的效果。k-ε模型和k-ω模型是应用广泛的两方程模型，k-ε模型通过求解湍动能k和湍动能耗散率ε的输运方程来封闭RANS方程，适用于各种工程湍流问题的计算；k-ω模型则求解湍动能k和比耗散率ω的输运方程，在近壁区域的计算精度较高。2.4.2湍流的数值模拟方法直接数值模拟（DNS）不引入任何湍流模型，直接对Navier-Stokes方程进行数值求解，能够精确模拟湍流的所有尺度，包括从大尺度的涡旋到小尺度的耗散涡。DNS方法可以获得湍流流场中详细的物理信息，对于研究湍流的基本特性和物理机制具有重要意义。由于湍流包含了极宽范围的尺度，从宏观的大尺度涡旋到分子尺度的耗散，DNS需要极高的计算分辨率来捕捉所有尺度的流动细节，这导致计算量巨大，对计算机的内存和计算速度要求极高。目前，DNS仅适用于简单流动和低雷诺数情况，在实际工程应用中受到很大限制。大涡模拟（LES）直接模拟大尺度湍流结构，对于小尺度湍流结构则采用亚格子模型进行模拟。LES方法通过滤波函数将湍流运动分解为大尺度运动和小尺度运动，大尺度运动直接求解，小尺度运动对大尺度运动的影响通过亚三、单个定日镜表面风压及周围风场研究3.1引言多定日镜风场特性的研究是确保塔式太阳能热发电系统稳定运行和结构安全的关键。在多定日镜风场中，单个定日镜的风场特性是基础，其表面风压分布和周围风场情况直接影响着定日镜的受力状态和稳定性。深入研究单个定日镜风场特性，能够为多定日镜风场模拟提供准确的参数和边界条件，有助于更全面、精确地理解多定日镜系统中风的作用机制，进而为整个定日镜场的优化设计和安全评估提供有力支持。因此，对单个定日镜表面风压及周围风场的研究具有重要的理论和实际意义，是开展多定日镜风场数值模拟研究的重要基础。3.2风对结构的作用3.2.1大气边界层风特性大气边界层是地球表面与自由大气之间直接相互作用的气层，其特性对定日镜的风载荷有着重要影响。在大气边界层中，风速随高度呈现出特定的变化规律。一般来说，在近地面层，由于受到地面摩擦力的影响，风速较小，且风速随高度的增加而逐渐增大。根据相关研究和实际观测，风速剖面通常可以用幂律分布来描述：U(z)=U_{ref}(\frac{z}{z_{ref}})^{\alpha}其中，U(z)为高度z处的风速，U_{ref}为参考高度z_{ref}处的参考风速，\alpha为风速幂指数，其值与地面粗糙度和大气稳定度等因素有关。在平坦地形且大气处于中性稳定状态时，\alpha一般取值在0.16-0.25之间。在城市地区，由于建筑物等障碍物较多，地面粗糙度大，\alpha值可能会偏大；而在开阔的沙漠或海洋地区，地面粗糙度小，\alpha值相对较小。大气边界层中的湍流度也是一个重要特性。湍流度反映了风速的脉动程度，其定义为脉动风速的均方根与平均风速的比值。在近地面层，湍流度较大，随着高度的增加，湍流度逐渐减小。湍流度的大小会影响风对定日镜的作用力，较大的湍流度会使定日镜受到的风载荷产生脉动，增加结构的动力响应。例如，在强风天气且湍流度较高时，定日镜可能会出现明显的振动，这不仅会影响其光学性能，还可能对结构的疲劳寿命产生不利影响。大气边界层中风向也会发生变化。在近地面层，由于受到地形和建筑物等的影响，风向可能会出现不规则的波动；而在较高的高度，风向相对较为稳定。风向的变化会导致定日镜所受风压的方向和大小发生改变，在研究定日镜的风载荷时，需要考虑风向的不确定性。3.2.2风对结构的作用风对定日镜结构产生的作用力主要包括阻力、升力和扭矩等。阻力是风在定日镜表面产生的平行于来流方向的力，其大小与风速的平方、定日镜的迎风面积以及阻力系数有关。阻力系数与定日镜的形状、表面粗糙度和来流的雷诺数等因素相关。对于平板状的定日镜，在一定的雷诺数范围内，其阻力系数大致为1.2-1.4。阻力会使定日镜产生沿风向的位移和变形，当阻力过大时，可能会导致定日镜的支撑结构损坏。升力是风在定日镜表面产生的垂直于来流方向的力，其产生与定日镜表面的压力分布不均匀有关。在一些情况下，如定日镜的安装角度不合理或受到复杂风场的作用时，升力可能会变得较大。升力会使定日镜产生向上或向下的位移，影响其稳定性。如果升力过大，可能会导致定日镜脱离支撑结构，造成严重的安全事故。扭矩是由于风在定日镜表面的作用力分布不均匀而产生的，它会使定日镜绕其支撑轴发生转动。扭矩的大小与定日镜的形状、尺寸、表面风压分布以及支撑方式等因素有关。扭矩会对定日镜的跟踪精度产生影响，进而影响太阳能的收集效率。如果扭矩过大，还可能会损坏定日镜的跟踪机构。这些力对定日镜稳定性的影响是相互关联的。当定日镜受到较大的阻力和升力时，会增加其结构的应力和变形，从而降低结构的稳定性。扭矩的作用也会使定日镜的姿态发生变化，进一步影响其受力状态和稳定性。在设计定日镜时，需要综合考虑这些力的作用，通过合理的结构设计和参数优化，提高定日镜的抗风能力和稳定性。3.3风洞试验简介定日镜风洞试验旨在模拟真实风场环境，研究定日镜在风作用下的力学性能。试验模型通常采用有机玻璃等材料制作，以保证与实物在外形上的几何相似性，缩尺比一般根据试验条件和研究需求确定。在某定日镜风洞试验中，定日镜模型高1.2m，宽0.9m，缩尺比为1/10。模型镜面由8行3列共24块板组成，镜面正反面各布置了144个测压点，以便精确测量定日镜表面的风压分布。测试风剖面和湍流度剖面是风洞试验的重要内容。通过在风洞中设置格栅、尖梯、挡板、粗糙元等装置，可以模拟出符合荷载规范中规定的B类地貌（\alpha=0.16）大气边界层风场。利用美国TSI公司的IFA300热膜风速仪等设备，对模拟风场的风速和湍流度进行测量，确保模拟的风场平均风速剖面符合指数律，紊流度、紊流积分尺度、脉动风速功率谱均满足试验要求。风速谱与Karman谱的对比可以验证模拟风场的准确性。试验工况的设置需要考虑多种因素。风洞试验时，坐标系规定符合右手定则，风向角\beta绕y轴正向在0^{\circ}-180^{\circ}内，每15^{\circ}转动一次；竖向转角x变化绕x轴正向在0^{\circ}-90^{\circ}度，每10^{\circ}进行一次测量。这样可以全面研究不同风向和竖向转角条件下定日镜的风效应，测量工况共有130个。3.4建模和求解过程3.4.1计算模型及网格划分建立定日镜数值模型时，需准确描述定日镜的几何形状和结构特征。利用三维建模软件，根据定日镜的实际尺寸和形状，构建精确的几何模型。在网格划分方面，选择合适的网格类型和尺寸至关重要。常用的网格类型有结构化网格和非结构化网格。结构化网格具有规则的拓扑结构，计算效率高，但对于复杂几何形状的适应性较差；非结构化网格则能更好地适应复杂几何形状，但计算量相对较大。在定日镜数值模拟中，可根据定日镜的形状特点，采用混合网格划分方法，即在定日镜的主体部分使用结构化网格，以提高计算效率；在边界层和复杂几何部位，如镜架、支撑结构等，采用非结构化网格，以提高网格质量和计算精度。网格尺寸对计算精度和效率有显著影响。较小的网格尺寸可以更精确地捕捉流场细节，但会增加计算量和计算时间；较大的网格尺寸虽然计算效率高，但可能会导致计算精度下降。通过网格无关性验证，确定合适的网格尺寸。在某定日镜数值模拟中，通过对比不同网格尺寸下的计算结果，发现当网格尺寸减小到一定程度后，计算结果的变化不再明显，此时的网格尺寸即为合适的网格尺寸。3.4.2边界条件的设定入口边界条件通常设定为速度入口，根据试验或实际工况确定入口风速的大小和方向。在模拟大气边界层风场时，入口风速剖面可采用幂律分布进行设定，以模拟不同高度处的风速变化。出口边界条件一般设定为压力出口，假定出口处的压力为大气压力。壁面边界条件根据定日镜的实际情况进行设定。对于定日镜表面和支撑结构表面，采用无滑移壁面条件，即壁面上的流体速度为零。在近壁区域，流体的流动特性与远离壁面处有很大不同，需要对壁面进行特殊处理，以准确模拟壁面附近的流场。3.4.3湍流模型及壁面处理选择合适的湍流模型是准确模拟定日镜风场的关键。常见的湍流模型有k-ε模型、k-ω模型和SSTk-ω模型等。k-ε模型是应用广泛的两方程模型，通过求解湍动能k和湍动能耗散率ε的输运方程来封闭RANS方程，适用于各种工程湍流问题的计算。k-ω模型则求解湍动能k和比耗散率ω的输运方程，在近壁区域的计算精度较高。SSTk-ω模型结合了k-ε模型和k-ω模型的优点，在近壁区域和远场都能给出较好的计算结果。在定日镜风场模拟中，根据模拟的精度要求和计算资源，选择合适的湍流模型。对于一般的工程计算，k-ε模型通常能满足要求；当需要更精确地模拟近壁区域的流场时，可选择SSTk-ω模型。壁面处理方法常用的有壁面函数法和直接求解壁面附近的流动。壁面函数法通过引入经验公式，将壁面附近的流动与主流区的流动联系起来，从而简化近壁区域的计算。在壁面函数法中，常用的有标准壁面函数和增强壁面函数。标准壁面函数基于对数律假设，适用于高雷诺数流动；增强壁面函数则考虑了更多的物理因素，如压力梯度、曲率等，适用于更复杂的流动情况。直接求解壁面附近的流动需要使用非常细密的网格，计算量巨大，但能得到更准确的壁面附近流场信息。在实际应用中，根据计算精度和计算资源的限制，选择合适的壁面处理方法。3.4.4计算和求解数值计算流程一般包括初始化计算参数、求解控制方程、迭代计算和收敛判断等步骤。在初始化阶段，设定计算的时间步长、收敛精度等参数。选择合适的求解器进行计算，常见的求解器有基于有限体积法的求解器，如FLUENT等。在计算过程中，通过迭代求解控制方程，逐步逼近流场的真实解。当计算结果满足设定的收敛精度时，认为计算收敛，得到最终的流场结果。在计算过程中，需要合理设置参数。时间步长的选择要兼顾计算精度和计算效率，较小的时间步长可以提高计算精度，但会增加计算时间；较大的时间步长则可能导致计算不稳定。收敛精度的设定要根据实际需求确定，过高的收敛精度会增加计算量，过低的收敛精度则会影响计算结果的准确性。3.5计算结果及分析通过数值模拟，得到定日镜表面的原始风压等值线云图。从云图中可以直观地看出定日镜表面风压的分布情况。在迎风面，风压呈现出不均匀分布，靠近来流边缘处的风压较大，而在背风面，风压相对较小，且存在明显的负压区域。对定日镜周围的风环境进行分析，发现风速分布也呈现出复杂的特性。在定日镜的迎风面，风速逐渐减小，在背风面，由于气流的分离和再附，形成了复杂的涡旋结构。这些涡旋会对定日镜产生额外的作用力，影响其稳定性。通过计算风压系数，进一步分析其在定日镜表面的分布规律。风压系数的计算公式为：C_{p}=\frac{p-p_{0}}{\frac{1}{2}\rhoU^{2}}其中，C_{p}为风压系数，p为测点的压力，p_{0}为参考压力，\rho为空气密度，U为参考风速。在定日镜的迎风面，风压系数为正值，且在靠近来流边缘处风压系数较大；在背风面，风压系数为负值，且在某些区域负压较大。风压系数的分布与定日镜的形状、风向和来流条件等因素密切相关。3.6本章小结本章通过理论分析、风洞试验和数值模拟等方法，对单个定日镜表面风压及周围风场进行了深入研究。得到了定日镜表面风压的分布规律，包括迎风面和背风面的风压大小和分布特点。分析了定日镜周围的风环境特性，如风速分布和涡旋情况。通过计算风压系数，进一步明确了风压在定日镜表面的分布规律。这些研究成果对多定日镜风场模拟具有重要的启示。在多定日镜风场模拟中，可以将单个定日镜的风场特性作为基础，考虑定日镜之间的相互干扰效应，进一步研究多定日镜系统的风场特性。单个定日镜的研究结果也为定日镜的结构设计和优化提供了理论依据，通过合理调整定日镜的形状、尺寸和安装角度等参数，可以降低风载荷对定日镜的影响，提高定日镜的抗风能力和稳定性。四、定日镜间相互干扰问题的数值分析4.1引言在塔式太阳能热发电系统中，定日镜数量众多，它们紧密排列组成定日镜场。定日镜间的相互干扰会显著影响风场特性和定日镜所受的风载荷，进而对定日镜的结构安全和整个发电系统的稳定性产生影响。当定日镜间距过小时，前一排定日镜的尾流会影响后一排定日镜的风场，导致后一排定日镜所受的风压分布发生改变，增加结构的受力风险。定日镜的排列方式也会影响风的流动路径和干扰程度，不合理的排列可能会加剧风干扰，降低定日镜的性能。深入研究定日镜间的相互干扰问题，对于优化镜场布置、提高发电效率、保障系统安全稳定运行具有重要意义。它可以为定日镜场的设计提供科学依据，通过合理调整定日镜的间距和排列方式，减小风干扰的影响，降低定日镜的结构成本，提高整个发电系统的经济效益。4.2几种对比工况为了深入研究定日镜间的相互干扰问题，设定了多种对比工况。在定日镜间距方面，设置了小间距工况、中等间距工况和大间距工况。小间距工况下，定日镜间距较小，旨在研究定日镜在紧密排列时风干扰的严重程度以及对风载荷的影响；中等间距工况代表了实际工程中较为常见的间距设置，分析该工况下的风场特性和干扰情况，可为实际工程提供参考；大间距工况下，定日镜间距较大，用于对比不同间距对风干扰的缓解效果，探讨增大间距在减小风干扰方面的作用。在排列方式上，考虑了正方形排列工况、三角形排列工况和交错排列工况。正方形排列工况具有规则的排列形式，便于分析和计算，研究其风干扰特性可以为其他排列方式提供对比基础；三角形排列工况在一些工程中也有应用，分析其风场情况可以了解该排列方式在减少风干扰方面的优势和不足；交错排列工况通过改变定日镜的排列位置，试图打破风的规律性流动，研究其对风干扰的影响，探索更优的排列方式。不同工况的特点和研究目的各不相同。小间距工况的特点是风干扰强烈，能够突出风干扰对定日镜的不利影响，研究目的是确定在极端紧密排列情况下定日镜的受力极限和可能出现的问题；中等间距工况的特点是接近实际工程应用，研究目的是为实际工程中的定日镜间距设计提供直接的参考依据；大间距工况的特点是风干扰相对较小，研究目的是评估增大间距在改善风场和降低风载荷方面的效果，为优化镜场布局提供思路。正方形排列工况的特点是排列规则、易于分析，研究目的是建立基本的排列模型，分析规则排列下的风干扰规律；三角形排列工况的特点是在空间利用和抗风干扰方面可能具有一定优势，研究目的是挖掘该排列方式的潜力，为镜场布局提供更多选择；交错排列工况的特点是排列方式灵活，研究目的是探索通过改变排列方式来减小风干扰的可行性，寻找更有效的抗干扰排列策略。4.3计算工况及网格划分针对不同的对比工况，建立相应的计算模型。在建模过程中，准确描述定日镜的几何形状、尺寸以及它们在镜场中的位置关系。利用专业的三维建模软件，按照实际定日镜的参数构建模型，确保模型的准确性和真实性。网格划分是数值模拟中的关键步骤，它直接影响计算结果的精度和计算效率。采用混合网格划分方法，在定日镜表面和周围流场变化剧烈的区域，使用非结构化网格，以提高对复杂流场的捕捉能力；在远离定日镜的区域，使用结构化网格，以减少计算量。对于定日镜表面，采用细密的非结构化网格，确保能够准确捕捉壁面附近的边界层流动和压力变化；在定日镜之间的间隙以及尾流区域，根据流场的复杂程度，适当调整网格密度，保证对风干扰现象的有效模拟。通过网格无关性验证来确定合理的网格尺寸。选取不同网格尺寸进行模拟计算，对比计算结果。当网格尺寸减小到一定程度后，计算结果的变化不再明显，此时的网格尺寸即为合适的网格尺寸。在某定日镜间相互干扰的数值模拟中，通过逐步减小网格尺寸，发现当网格尺寸减小到一定值后，定日镜表面的风压系数和周围流场的速度分布等计算结果基本保持不变，从而确定了合适的网格尺寸。合理的网格划分既能够保证计算精度，又能控制计算成本，提高计算效率，为准确分析定日镜间的相互干扰问题提供了有力支持。4.4数值模拟与风洞试验的结果对比将数值模拟结果与风洞试验数据进行对比，是验证数值模拟方法准确性的重要手段。通过对比，可以评估数值模拟方法在预测定日镜间风干扰效应方面的可靠性，为进一步的数值研究提供依据。在对比过程中，选取关键参数进行分析，如定日镜表面的风压系数、周围流场的速度分布等。对于风压系数，对比不同工况下数值模拟和试验得到的风压系数分布规律和数值大小。在某工况下，数值模拟得到的定日镜迎风面风压系数与风洞试验结果在趋势上基本一致，都呈现出在边缘区域风压系数较大，中心区域相对较小的特点，但在具体数值上可能存在一定差异。对于流场速度分布，对比不同位置处的速度大小和方向。在定日镜的尾流区域，数值模拟得到的速度衰减趋势与风洞试验结果相符，但在速度的具体数值上可能存在偏差。分析两者差异的原因，主要包括以下几个方面。数值模拟中采用的湍流模型虽然能够模拟湍流的一般特性，但对于复杂的定日镜间流场，可能无法完全准确地描述湍流的细节，从而导致计算结果与试验结果存在差异。网格划分的精度也会影响计算结果。尽管通过网格无关性验证确定了合理的网格尺寸，但在实际计算中，网格仍然是对连续流场的离散近似，可能会引入一定的误差。风洞试验本身也存在一定的误差，如测量仪器的精度、模型制作的误差等，这些因素都会对试验结果产生影响，进而导致与数值模拟结果的差异。边界条件的设定在数值模拟中也具有重要影响。实际风场的边界条件较为复杂，而数值模拟中往往需要进行一定的简化，这种简化可能会导致计算结果与实际情况存在偏差。4.5风压系数干扰效应4.5.1定日镜镜面定日镜镜面风压系数受周围定日镜干扰的变化情况较为复杂。在小间距工况下，由于周围定日镜的尾流影响，镜面风压系数分布明显改变。在前排定日镜的遮挡下，后排定日镜迎风面的风压系数减小，且分布不均匀性增加。在某小间距工况中，后排定日镜迎风面靠近边缘处的风压系数相比无干扰情况下降低了约20%，且在不同位置处的风压系数差异增大，这是因为前排定日镜的尾流导致气流紊乱，使得后排定日镜表面的压力分布更加复杂。干扰对镜面受力产生显著影响。风压系数的变化直接导致镜面所受压力的改变，不均匀的风压分布会使镜面产生不均匀的应力，增加了镜面变形和损坏的风险。当风压系数在镜面上分布不均匀时，镜面不同部位所受的压力不同，可能会导致镜面产生弯曲变形，长期作用下甚至可能使镜面出现裂纹或破裂。干扰还可能使镜面所受的合力方向发生改变，对镜面的支撑结构提出更高的要求，若支撑结构无法承受这种变化的力，可能会导致镜面倾斜或倒塌。4.5.2定日镜背面定日镜背面风压系数的变化与周围定日镜的干扰密切相关。在大间距工况下，定日镜背面的风压系数相对较小，且分布较为均匀。随着间距减小，干扰增强，背面风压系数增大，且分布变得不均匀。在某中等间距工况下，定日镜背面的风压系数相比大间距工况增大了约15%，且在靠近相邻定日镜的一侧，风压系数明显高于其他部位，这是由于相邻定日镜的存在改变了气流的流动路径，使得定日镜背面的气流分离和再附现象更加复杂。分析背风面的压力分布和气流特征，发现背风面存在明显的负压区域。在干扰作用下，负压区域的范围和强度发生变化。当定日镜间距较小时，负压区域范围扩大，强度增强，这会增加定日镜所受的吸力，对定日镜的稳定性产生不利影响。在小间距工况下，定日镜背风面的负压区域范围比大间距工况扩大了约30%，负压强度也增加了约20%，这使得定日镜更容易受到向上的吸力作用，可能导致定日镜脱离支撑结构。气流在背风面形成复杂的涡旋结构，这些涡旋的存在会进一步影响风压系数的分布和定日镜的受力。涡旋的旋转和移动会导致背风面压力的脉动，增加定日镜的振动和疲劳损伤风险。4.6流场干扰效应4.6.1水平剖面在水平剖面上，定日镜间流场的干扰情况十分明显。当定日镜采用正方形排列时，在定日镜之间的间隙处，气流会发生汇合和分离现象。在某正方形排列工况中，相邻定日镜间隙处的气流速度明显增大，这是因为气流在通过狭窄的间隙时，受到收缩作用，流速加快。在定日镜的尾流区域，气流速度逐渐减小，形成低速区。由于相邻定日镜的尾流相互影响，低速区的范围和形状发生改变。在小间距工况下，相邻定日镜的尾流相互叠加，低速区范围扩大，气流的恢复速度减慢，这会影响后续定日镜的风场，增加风干扰的程度。4.6.2垂直剖面研究垂直剖面上流场的变化，发现不同高度处气流的相互作用显著。在靠近地面的区域，由于地面的摩擦作用，气流速度较低，且受到定日镜的影响较大。随着高度的增加，气流速度逐渐增大，定日镜对气流的影响逐渐减小。在某垂直剖面分析中，在