CN104679953B - 户内变电站主变室流场温度场数值模拟快速计算方法、系统 - Google Patents
户内变电站主变室流场温度场数值模拟快速计算方法、系统 Download PDFInfo
- Publication number
- CN104679953B CN104679953B CN201510076326.7A CN201510076326A CN104679953B CN 104679953 B CN104679953 B CN 104679953B CN 201510076326 A CN201510076326 A CN 201510076326A CN 104679953 B CN104679953 B CN 104679953B
- Authority
- CN
- China
- Prior art keywords
- mrow
- module
- msup
- msub
- temperature
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 28
- 238000004088 simulation Methods 0.000 title claims abstract description 18
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 72
- 238000000034 method Methods 0.000 claims abstract description 25
- 230000004888 barrier function Effects 0.000 claims description 12
- 125000006850 spacer group Chemical group 0.000 claims description 9
- 238000000638 solvent extraction Methods 0.000 claims description 8
- 230000004907 flux Effects 0.000 claims description 4
- 230000005540 biological transmission Effects 0.000 claims description 3
- 238000009792 diffusion process Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 claims description 3
- 238000012546 transfer Methods 0.000 claims description 3
- 241000208340 Araliaceae Species 0.000 claims description 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 claims description 2
- 235000008434 ginseng Nutrition 0.000 claims description 2
- 230000009466 transformation Effects 0.000 claims 1
- 230000008569 process Effects 0.000 abstract description 9
- 238000005516 engineering process Methods 0.000 description 6
- 239000012530 fluid Substances 0.000 description 6
- 230000008901 benefit Effects 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000012800 visualization Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000012067 mathematical method Methods 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000009423 ventilation Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Abstract
本发明公开了一种户内变电站主变室流场温度场数值模拟快速计算方法，首先是建立输入参数表和输出参数表；再进行几何建模和网格划分；之后将网格划分结果输入到速度场快速求解算法以及温度场快速求解算法中进行计算，结果存储到输出参数表中。本发明还公开了一种户内变电站主变室流场温度场数值模拟快速计算系统，包括了参数化外观建模模块、智能网格生成模块、湍流方程快速求解模块、温度场快速求解模块、速度场快速求解模块、输入输出模块和用户输入模块。本发明实现了户内变电站主变室流场温度场数值的快速计算，使得建模和网格划分过程变得快速、准确，且占用的内存大大的减小，以方便实际工程应用。
Description
技术领域
本发明涉及一种能源系统动态特性及控制领域，特别是一种户内变电站主变室流场温度场数值模拟快速计算方法、系统。
背景技术
计算流体力学技术：计算流体力学(CFD)是近代流体力学，数值数学和计算机科学结合的产物，它以电子计算机为工具，应用各种离散化的数学方法，对流体力学的各类问题进行数值实验、计算机模拟和分析研究，以解决各种实际问题。
由于计算流体力学技术具有成本低和能模拟较复杂或较理想的过程等优点，因此采用计算流体力学的方法对户内变电站主变室的流场和温度场进行分析具有重要的工程实用价值。但是常规的CFD计算中，前期的外观建模、绘制网格以及后期计算中所涉及到的求解方法需要耗费大量的计算时间和计算内存，不仅对使用人员的建模能力、网格绘制经验的要求较高，而且对于计算机的内存空间和计算性能也有很高的要求，计算耗时长、内存占用量大、计算结果受人为因素影响严重等问题已称为阻碍CFD计算推广应用的主要瓶颈。
发明内容
本发明的其中一个目的就是提供一种户内变电站主变室流场温度场数值模拟快速计算方法，可以实现智能的网格划分过程，大幅缩短计算求解时长以及减小计算机所占内存，为户内变电站主变室流场与温度场的计算提供便捷、快速且满足工程精度需要的数值计算手段。
本发明的目的之一是通过这样的技术方案实现的，它包括有如下的具体步骤：
1)根据参数化外观建模算法以及智能网格划分算法、速度场快速求解算法以及温度场快速求解算法所需的输入参数和输出参数，建立输入、输出参数表；
2)根据参数化外观建模算法、智能网格划分算法对户内变电站主变室的计算空间进行外观建模和网格划分；
3)将上述网格划分的结果输入到速度场快速求解算法以及温度场快速求解算法中进行计算，并将流场的计算结果存储到输出参数表；
4)将输出参数表中的数据输入到流场可视化程序中；
步骤1)中所述的输入参数表用于存储流场温度场数值模拟计算的参数，包括户内变电站主变室参数化建模和网格划分所需的整体尺寸参数、户内变电站主变室用于内部的障碍物及边界尺寸参数、速度场快速求解算法以及温度场快速求解算法所需的参数；所述的输出参数表用于存储计算结果；
步骤2)中所述的参数化外观建模算法，具体步骤为：
2-1-1)读取用户输入的变电站主变室整体的尺寸信息；
2-1-2)经系统自动检查无误后，依次读入障碍物的位置及尺寸信息，进出口位置、尺寸及方向信息，风扇的位置、尺寸及方向信息；
2-1-3)自动完成变电站主变室的外观建模，并且将外观建模的信息传递到智能网格生成算法中；
步骤2)中所述的智能网格划分算法的具体步骤为：
2-2-1)读入由参数化外观建模算法生成的信息，包括障碍物位置及尺寸信息、进出口位置以及尺寸和方向信息、风扇的位置以及尺寸和方向信息、还有从输入参数表中读取的X、Y、Z三个坐标轴的网格整体数量；
2-2-2)采用障碍物和边界的坐标值将X、Y、Z三个坐标轴分成若干段；
2-2-3)将分成的每一段的长度逐一求出；
步骤3)中所述的速度场快速求解算法包括密度简化算法、微分方程分布投影算法以及湍流方程求解快速算法；
所述的密度简化法采用的是Boussinesq近似，即将密度以及流场的物性参数均视为常数，其公式为：其中R为气体常数，取值为R＝287.1J/(Kg·K)，公式中T0为参考温度，取值为室外温度值，P0为参考压力值，取值为一个大气压；
进一步所述的微分方程分布投影算法，将时间步从n推进到n+1的方法是：首先从时间步n推进到对流-扩散中间步*，应用的公式为：
其中u*为中间步*的速度项，un为当前时间步n的速度项，un-1为上一个时间步的速度项，V代表网格单元面中心的速度矢量，用来计算面元的通量，中间步的推进完成之后，V*通过格点中心的u*的插值计算，Vn-1为上一个时间步的V，而Vn+1为下一个时间步长的V，vt为湍动强度,△t为时间步长，Tn为当前时间步的温度项，S为控制面积,为拉普拉斯算子；
所述的湍流方程快速求解算法采用的公式为：vT＝0.03874VL，其中u1为流场中任意一点x沿X轴的速度，u2为点x沿Y轴的速度，u3为点x沿Z轴的速度，L为长度尺度，取值为点x与壁面之间的最短距离；
所述的温度场快速求解算法的主要公式为：
其中k为传热系数，T*为中间时间步的温度项，Tn为当前时间步的温度项,Tn-1为上一个时间步的温度项，ST为源项，Vn+1为下一个时间步长的V。
进一步，所述的输入参数表可以自动完成对用户输入界面输入的数据存储。
本发明的另一个目的在于提供了一种户内变电站主变室流场温度场数值模拟快速计算系统，它可以自动完成几何建模和网格划分，并且温度场快速求解模块、速度场快速求解模块以及湍流方程快速求解模块的使用，大幅缩减了计算量。
本发明是通过这样的技术方案实现的，它包括有参数化外观建模模块、智能网格划分模块、湍流方程快速求解模块、温度场快速求解模块以及速度场快速求解模块、流场可视化模块、输入参数模块、输出参数模块、用户输入模块；
输入参数模块接收输入到用户输入模块中的户内变电站主变室几何参数，然后经参数化外观建模模块进行几何建模，参数化建模数据输入智能网格划分模块自动完成网格划分，将智能网格划分模块中输出的边界信息以及障碍物信息输入到湍流方程快速求解模块进行计算，将其余划分结果输入到速度场快速求解模块以及温度场快速求解模块中进行计算，再将计算结果输入到输出参数模块中，最后将输出参数表中接收到的计算结果输入到流场可视化模块。
进一步，所述的输入参数模块与用户输入模块以及参数化外观建模模块相连，用于存储用户输入模块中输入的户内变电站主变室流场温度场计算所需的尺寸、障碍物及边界信息；所述的输出参数模块用于记录智能网格划分模块的网格划分结果、户内变电站主变室的尺寸、温度场快速求解模块以及速度场快速求解模块求解出的流场温度场信息结果。
进一步，所述的流场温度场信息结果包括：变电站主变室流场中所有点的温度值、压力值以及沿X、Y、Z轴向方向的速度。
由于采用了上述技术方案，本发明具有如下的优点：
1)现有的诸多CFD技术和系统需要使用者自行建模以及网格划分，而这一部分的工作量往往很大，据统计，外观建模以及网格划分需要占用使用者接近80％的工作时间，而真正从系统开始计算一直到算出最终结果只占到整体开发过程的10％，另外在后处理中使用了10％的时间。变电站主变室本身具有空间大、内部结构复杂，扰流元件多等特点，决定了其前处理过程,即外观建模以及网格生成更加的耗时耗力，而本算法和系统的开发是的前处理过程，即建模和网格划分过程可以通过用户参数化输入之后由参数化建模模块以及智能网格划分模块自动完成，十分方便快捷，用户只需按照附件中表1.1的参数输入表输入变电站主变室整体的尺寸、障碍物的位置尺寸以及边界的位置和信息就可以进行自动生成高质量的外观模型以及网格。
2)在完成了外观建模和网格划分之后的计算过程，现有的CFD技术也因为需要求解大量的偏微分方程而使得整体的计算时间长和计算占用内存大，而本发明中提到的温度场快速求解模块和速度场快速求解模块以及湍流方程快速求解模块由于使用了更为简单的求解过程，大大缩减了偏微分方程的求解量，使得求解过程十分快速。
本发明的其他优点、目标和特征在某种程度上将在随后的说明书中进行阐述，并且在某种程度上，基于对下文的考察研究对本领域技术人员而言将是显而易见的，或者可以从本发明的实践中得到教导。本发明的目标和其他优点可以通过下面的说明书和权利要求书来实现和获得。
附图说明
本发明的附图说明如下。
图1为本发明的整体算法流程图。
图2为参数化外观建模算法流程图。
图3为智能网格生成算法流程图。
图4为湍流方程快速计算流程图。
图5为本发明的结构示意图。
具体实施方式
下面结合附图和实施例对本发明作进一步说明。
本发明提供的是一种户内变电站主变室流场温度场数值模拟快速计算方法、系统，目的在于针对所有的计算流体力学技术的前期外观建模复杂、网格划分过程不智能、计算求解过程耗时长以及占用内存大的问题，为户内变电站主变室流场温度场的计算提供了便捷快速且满足工程精度的数值计算方法。
一种户内变电站主变室流场温度场数值模拟快速计算方法，先根据参数化外观建模算法以及智能网格划分算法、速度场快速求解算法、温度场快速求解算法所需参数以及输出参数，建立输入参数表和输出参数表，所述输入参数表存储该算法下用于流场温度场模拟计算的参数，包括变电站主变室参数化建模和网格生成所需的整体尺寸参数和户内变电站主变室内部的障碍物以及边界尺寸参数以及速度场快速求解算法、温度场快速求解算法所需要的参数，所述输出表用于存储该算法的计算结果，并且将计算所得数据用于流场的可视化分析之中；再根据参数化外观建模算法以及智能网格划分算法对户内变电站主变室的计算空间进行外观建模和网格划分，将网格划分的结果输入到速度场快速求解算法以及温度场快速求解算法中进行计算，将流场的计算结果存储到对应的输出参数表，然后将输出参数表中的数据输入到可视化程序中，以此完成流场的可视化。其中所述的参数输入表可见于附件中表1.1，参数输出表可见于附件表1.2。
上述的参数化外观建模算法中，首先读取用户输入的变电站主变室整体的尺寸信息，经检查无误，依次读入障碍物信息、进出口信息、风扇信息，然后自动将变电站主变室的外观建模完成。并且将外观建模的信息传递到智能网格生成算法中。
上述的网格自动生成算法分为三个步骤，首先读入由参数化外观建模算法生成的信息，这些信息包括障碍物信息、进出口信息、风扇信息以及从输入参数表中读取X、Y、Z三个坐标轴的网格整体数量；然后使用障碍物和边界的坐标值将X、Y、Z轴分成若干段；最后将每一段的长度逐一求出。
上述的速度场快速求解算法包含密度简化算法、微分方程分布投影算法以及湍流方程求解快速算法。
所述的密度简化法采用的是Boussinesq近似，即将密度和流场的物性参数均看作常数，其公式为：其中R为气体常数，取值为R＝287.1J/(Kg·K)，公式中T0为参考温度，取值为室外温度值，P0为参考压力值，取值为一个大气压。
所述的微分方程分布投影算法，将时间步从n推进到n+1的方法是：首先从时间步n推进到对流-扩散中间步*，应用的公式为：
其中u*为中间步*的速度项，un为当前时间步n的速度项，un-1为上一个时间步的速度项，V代表网格单元面中心的速度矢量，用来计算面元的通量，中间步的推进完成之后，V*通过格点中心的u*的插值计算，Vn-1为上一个时间步的V，而Vn+1为下一个时间步长的V，vt为湍动强度,△t为时间步长，Tn为当前时间步的温度项，S为控制面积,为拉普拉斯算子。
进一步，所述的湍流方程快速求解算法采用的公式为：vT＝0.03874VL，其中u1为流场中任意一点x沿X轴的速度，u2为点x沿Y轴的速度，u3为点x沿Z轴的速度，L为长度尺度，取值为点x与壁面之间的最短距离。
进一步，所述的温度场快速求解算法的主要公式为：
其中k为传热系数，T*为中间时间步的温度项，Tn为当前时间步的温度项,Tn-1为上一个时间步的温度项，ST为源项，Vn+1为下一个时间步长的V。
根据上文所述的户内变电站主变室流场温度场数值模拟的快速计算方法，其特征在于，利用建立的输入参数表可以自动完成对用户输入界面输入数据的存储，其中所述的速度场快速求解算法、温度场快速求解算法接收输入参数表中的内容并且进行计算，能够快速地得到整体流场的速度、温度以及压力场。
根据前文所述的一种户内变电站主变室通风换热性能计算分析方法的整体流程图可见于附件中图1.
一种户内变电站主变室流场温度场数值模拟快速计算系统，其特征在于，包括依次连接的参数化外观建模模块、智能网格划分模块、湍流方程快速求解模块、温度场快速求解模块以及速度场快速求解模块，还包括与输入输出模块相连的流场可视化模块以及用户输入模块，模块之间的连接方式如附图5中所示。所述输入输出模块包括输入参数表和输出参数表，所述输入参数表和用户输入模块以及参数化外观建模模块相连，其用于存储用户输入模块中输入的关于户内变电站主变室流场温度场计算所需的尺寸、障碍物以及边界的信息。所述的输入表如附件表1.1所示。所述的输出参数表记录了智能网格划分模块的网格划分结果、关于户内变电站主变室的尺寸，温度场快速求解模块以及速度场快速求解模块求解的流场温度场信息结果，所述的流场温度场信息结果包括了变电站主变室流场中所有点的温度值、压力值、沿X、Y、Z轴向方向的速度。所述输出参数表见附件表1.2。
所述的参数化外观建模模块流程图如附件中图2所示，按照流程图中的顺序依次从输入参数表中将所需的参数输入模块并且按流程图开始计算。
所述的智能网格划分模块计算流程图如附件中图3所示，按照流程图中的顺序依次从输入参数表以及参数化建模算法生成的数据中进行数据输入并且按照流程开始计算。
所述的湍流方程快速求解模块计算流程图如附件中图4所示，按照流程图中的顺序依次从输入智能网格划分模块中输出的边界以及障碍物点的信息读入开始，按照流程进行计算。
表1.1输入参数表
表1.2输出参数表
最后说明的是，以上实施例仅用以说明本发明的技术方案而非限制，尽管参照较佳实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，可以对本发明的技术方案进行修改或者等同替换，而不脱离本技术方案的宗旨和范围，其均应涵盖在本发明的权利要求范围当中。
Claims (5)
1.一种户内变电站主变室流场温度场数值模拟快速计算方法，其特征在于，具体步骤如下：
1)根据参数化外观建模算法、智能网格划分算法、速度场快速求解算法、温度场快速求解算法建立输入参数表和输出参数表；
2)根据参数化外观建模算法以及智能网格划分算法对户内变电站主变室的计算空间进行外观建模和网格划分；
3)将上述网格划分的结果输入到速度场快速求解算法以及温度场快速求解算法中进行计算，将流场的计算结果存储到输出参数表中；
4)将输出参数表中的数据输入到流场可视化程序中；
步骤1)中所述的输入参数表存储用于流场温度场数值模拟计算的参数，包括户内变电站主变室参数化建模和网格划分所需的整体尺寸参数、户内变电站主变室内部的障碍物及边界尺寸参数，速度场快速求解算法以及温度场快速求解算法所需的参数；所述的输出参数表用于存储计算结果；
步骤2)中所述的参数化外观建模算法，具体步骤为：
2-1-1)读取用户输入的变电站主变室整体的尺寸信息；
2-1-2)经系统自动检查无误后，依次读入障碍物的位置及尺寸信息，进出口位置、尺寸及方向信息，风扇的位置、尺寸及方向信息；
2-1-3)自动完成变电站主变室的外观建模，并且将外观建模的信息传递到智能网格生成算法中；
步骤2)中所述的智能网格划分算法的具体步骤为：
2-2-1)读入由参数化外观建模算法生成的信息，包括障碍物位置及尺寸信息，进出口位置、尺寸及方向信息，风扇的位置、尺寸及方向信息，还有从输入参数表中读取的X、Y、Z三个坐标轴的网格整体数量；
2-2-2)采用障碍物和边界的坐标值将X、Y、Z三个坐标轴分成若干段；
2-2-3)将分成的每一段的长度逐一求出；
步骤3)中所述的速度场快速求解算法包括密度简化算法、微分方程分布投影算法以及湍流方程快速求解算法；
所述的密度简化法采用的是Boussinesq近似，即将密度看作常数，流场的物性参数也看作常数，其公式为：其中R为气体常数，取值为R＝287.1J/(Kg·K)，公式中T0为参考温度，取值为室外温度值，P0为参考压力值，取值为一个大气压；
所述的微分方程分布投影算法，将时间步从n推进到n+1的方法是：首先从时间步n推进到对流-扩散中间步*，应用的公式为：
<mrow>
<msub>
<mo>&Integral;</mo>
<mrow>
<mi>c</mi>
<mi>v</mi>
</mrow>
</msub>
<mfrac>
<mrow>
<msup>
<mi>u</mi>
<mo>*</mo>
</msup>
<mo>-</mo>
<msup>
<mi>u</mi>
<mi>n</mi>
</msup>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mi>d</mi>
<mi>V</mi>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mo>(</mo>
<mrow>
<msub>
<mo>&Integral;</mo>
<mrow>
<mi>c</mi>
<mi>s</mi>
</mrow>
</msub>
<mn>3</mn>
<msup>
<mi>u</mi>
<mi>n</mi>
</msup>
<msup>
<mi>V</mi>
<mi>n</mi>
</msup>
<mi>d</mi>
<mi>S</mi>
<mo>-</mo>
<msub>
<mo>&Integral;</mo>
<mrow>
<mi>c</mi>
<mi>s</mi>
</mrow>
</msub>
<msup>
<mi>u</mi>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msup>
<mi>V</mi>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mi>d</mi>
<mi>S</mi>
</mrow>
<mo>)</mo>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mo>(</mo>
<mrow>
<msub>
<mo>&Integral;</mo>
<mrow>
<mi>c</mi>
<mi>s</mi>
</mrow>
</msub>
<msub>
<mi>v</mi>
<mi>t</mi>
</msub>
<mo>&dtri;</mo>
<msup>
<mi>u</mi>
<mo>*</mo>
</msup>
<mi>d</mi>
<mi>S</mi>
<mo>+</mo>
<msub>
<mo>&Integral;</mo>
<mrow>
<mi>c</mi>
<mi>s</mi>
</mrow>
</msub>
<msub>
<mi>v</mi>
<mi>t</mi>
</msub>
<mo>&dtri;</mo>
<msup>
<mi>u</mi>
<mi>n</mi>
</msup>
<mi>d</mi>
<mi>S</mi>
</mrow>
<mo>)</mo>
<mo>+</mo>
<msub>
<mo>&Integral;</mo>
<mrow>
<mi>c</mi>
<mi>s</mi>
</mrow>
</msub>
<mo>&lsqb;</mo>
<mo>-</mo>
<mi>g</mi>
<mi>&beta;</mi>
<mrow>
<mo>(</mo>
<msup>
<mi>T</mi>
<mi>n</mi>
</msup>
<mo>-</mo>
<msub>
<mi>T</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mi>d</mi>
<mi>V</mi>
</mrow>
其中u*为中间步*的速度项，un为当前时间步n的速度项，un-1为上一个时间步的速度项，V代表网格单元面中心的速度矢量，用来计算面元的通量，中间步的推进完成之后，V*通过格点中心的u*的插值计算，Vn-1为上一个时间步的V，而Vn+1为下一个时间步长的V，vt为湍动强度,△t为时间步长，Tn为当前时间步的温度项，S为控制面积,为拉普拉斯算子；
所述的湍流方程快速求解算法采用的公式为：vT＝0.03874VL，其中u1为流场中任意一点x沿X轴的速度，u2为点x沿Y轴的速度，u3为点x沿Z轴的速度，L为长度尺度，取值为点x与壁面之间的最短距离；
所述的温度场快速求解算法的主要公式为：
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<msup>
<mi>T</mi>
<mo>*</mo>
</msup>
<mo>-</mo>
<msup>
<mi>T</mi>
<mi>n</mi>
</msup>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>u</mi>
<mi>j</mi>
</msub>
<mfrac>
<mrow>
<mo>&part;</mo>
<msup>
<mi>T</mi>
<mi>n</mi>
</msup>
</mrow>
<mrow>
<mo>&part;</mo>
<msub>
<mi>x</mi>
<mi>j</mi>
</msub>
</mrow>
</mfrac>
</mrow>
</mtd>
<mtd>
<mrow>
<mfrac>
<mrow>
<msup>
<mi>T</mi>
<mrow>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>-</mo>
<msup>
<mi>T</mi>
<mo>*</mo>
</msup>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<mfrac>
<msub>
<mi>v</mi>
<mi>T</mi>
</msub>
<msub>
<mi>Pr</mi>
<mi>T</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<msup>
<mi>T</mi>
<mrow>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msup>
</mrow>
<mrow>
<mo>&part;</mo>
<msubsup>
<mi>x</mi>
<mi>j</mi>
<mn>2</mn>
</msubsup>
</mrow>
</mfrac>
<mo>+</mo>
<msub>
<mi>S</mi>
<mi>T</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>;</mo>
</mrow>
其中k为传热系数，T*为中间时间步的温度项，Tn为当前时间步的温度项,Tn-1为上一个时间步的温度项，ST为源项，Vn+1为下一个时间步长的V。
2.如权利要求1所述的户内变电站主变室流场温度场数值模拟快速计算方法，其特征在于，所述的输入参数表可以自动完成对用户输入界面输入数据的存储。
3.实现权利要求1和2所述方法的一种户内变电站主变室流场温度场数值模拟快速计算系统，其特征在于，包括：参数化外观建模模块、智能网格划分模块、湍流方程快速求解模块、温度场快速求解模块以及速度场快速求解模块、流场可视化模块、输入参数模块、输出参数模块、用户输入模块；
输入参数模块接收输入到用户输入模块中的户内变电站主变室几何参数，然后经参数化外观建模模块进行几何建模，参数化建模数据输入智能网格划分模块自动完成网格划分，将智能网格划分模块中输出的边界信息以及障碍物信息输入到湍流方程快速求解模块进行计算，将其余划分结果输入到速度场快速求解模块以及温度场快速求解模块中进行计算，再将计算结果输入到输出参数模块中，最后将输出参数表中接收到的计算结果输入到流场可视化模块。
4.如权利要求3所述的户内变电站主变室流场温度场数值模拟快速计算系统，其特征在于，所述的输入参数模块与用户输入模块以及参数化外观建模模块相连，用于存储用户输入模块中输入的户内变电站主变室流场温度场计算所需的尺寸、障碍物及边界信息；所述的输出参数模块用于记录智能网格划分模块的网格划分结果、户内变电站主变室的尺寸、温度场快速求解模块以及速度场快速求解模块求解出的流场温度场信息结果。
5.如权利要求4所述的户内变电站主变室流场温度场数值模拟快速计算系统，其特征在于，所述的流场温度场信息结果包括：变电站主变室流场中所有点的温度值、压力值和沿X、Y、Z三个轴向方向的速度。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510076326.7A CN104679953B (zh) | 2015-02-12 | 2015-02-12 | 户内变电站主变室流场温度场数值模拟快速计算方法、系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510076326.7A CN104679953B (zh) | 2015-02-12 | 2015-02-12 | 户内变电站主变室流场温度场数值模拟快速计算方法、系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104679953A CN104679953A (zh) | 2015-06-03 |
| CN104679953B true CN104679953B (zh) | 2017-12-15 |
Family
ID=53314992
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510076326.7A Active CN104679953B (zh) | 2015-02-12 | 2015-02-12 | 户内变电站主变室流场温度场数值模拟快速计算方法、系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104679953B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105354438B (zh) * | 2015-12-09 | 2018-09-04 | 国网重庆市电力公司电力科学研究院 | 一种考虑地形影响的uhvdc线路地面三维合成电场计算方法 |
| CN110046440B (zh) * | 2019-04-22 | 2021-03-12 | 国电联合动力技术有限公司 | 基于cfd数据库的流场模拟计算方法及装置 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103166133A (zh) * | 2013-02-27 | 2013-06-19 | 重庆市电力公司电力科学研究院 | 户内变电站主变压器室通风换热优化控制方法 |
| CN104036125A (zh) * | 2014-05-30 | 2014-09-10 | 河北省电力建设调整试验所 | 一种油浸式变压器内部温度场的精确计算方法 |
-
2015
- 2015-02-12 CN CN201510076326.7A patent/CN104679953B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103166133A (zh) * | 2013-02-27 | 2013-06-19 | 重庆市电力公司电力科学研究院 | 户内变电站主变压器室通风换热优化控制方法 |
| CN104036125A (zh) * | 2014-05-30 | 2014-09-10 | 河北省电力建设调整试验所 | 一种油浸式变压器内部温度场的精确计算方法 |
Non-Patent Citations (4)
| Title |
|---|
| The CFD Numerical Simulation on Ventilation of Indoor Transformer Substation;XU Luwen 等;《Hydromechatronics Engineering》;20130930;第41卷(第18期);第141-144页 * |
| 户内散热器室自然通风数值分析;邵志伟 等;《河南科技大学学报：自然科学版》;20141231;第35卷(第6期);第66-70页 * |
| 辐射对 220kV 主变压器室温度场分布的模拟研究;王香丽 等;《发电与空调》;20120115;第33卷(第1期);第65-68页 * |
| 进风口对变压器室通风效果的影响;高学平 等;《土木建筑与环境工程》;20121231;第34卷;第117-122页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104679953A (zh) | 2015-06-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Daniels et al. | A suite of computationally expensive shape optimisation problems using computational fluid dynamics | |
| Shen et al. | Investigation of response surface methodology for modelling ventilation rate of a naturally ventilated building | |
| CN108647370A (zh) | 基于双环迭代的无人直升机气动外形优化设计方法 | |
| CN104933225B (zh) | 实现计算流体力学大规模实时模拟的方法 | |
| CN109614638B (zh) | 一种非直接建模的城市风环境cfd模拟方法 | |
| CN110348059A (zh) | 一种基于结构化网格的通道内流场重构方法 | |
| CN109948214B (zh) | 城市多尺度风环境数值模拟方法 | |
| CN113177370B (zh) | 考虑环境湿度因素风流-粉尘气固两相流动数值模拟方法 | |
| Ivanell et al. | Micro-scale model comparison (benchmark) at the moderately complex forested site Ryningsnäs | |
| Amoiralis et al. | Freeform deformation versus B-spline representation in inverse airfoil design | |
| CN104679953B (zh) | 户内变电站主变室流场温度场数值模拟快速计算方法、系统 | |
| CN101436303B (zh) | 从物体三维图像中获取四面体网格的方法 | |
| CN103778098A (zh) | 基于格子Boltzmann理论实现协同计算大涡模拟系统及方法 | |
| US20220414282A1 (en) | Boundary layer mesh generation method based on anisotropic volume harmonic field | |
| Zhang et al. | Three-dimensional wind field reconstruction using tucker decomposition with optimal sensor placement | |
| Zhu et al. | Hydrodynamic design of a circulating water channel based on a fractional-step multi-objective optimization | |
| CN112528538A (zh) | 一种室内新风仿真模拟方法及室内设计系统 | |
| Andreev et al. | A computational technology for constructing the optimal shape of a power plant blade assembly taking into account structural constraints | |
| CN100447816C (zh) | Ct投影数据三维解析模拟方法 | |
| CN113673177B (zh) | 网格空隙空间识别和自动种子设定检测 | |
| Ma | Multimedia simulation-based architecture cad system model | |
| CN103778289B (zh) | 一种风速流线可视化的建筑设计方法 | |
| CN113505544A (zh) | 基于有限体积法自行车运动虚拟数值风洞系统 | |
| CN113887147A (zh) | 一种基于mri数据的脑血流自动化分析系统 | |
| Su et al. | AWSD: An Aircraft Wing Dataset Created by an Automatic Workflow for Data Mining in Geometric Processing. |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |