CN106570284B - 一种聚焦式太阳能热发电系统的多维分析方法 - Google Patents

一种聚焦式太阳能热发电系统的多维分析方法 Download PDF

Info

Publication number
CN106570284B
CN106570284B CN201610986195.0A CN201610986195A CN106570284B CN 106570284 B CN106570284 B CN 106570284B CN 201610986195 A CN201610986195 A CN 201610986195A CN 106570284 B CN106570284 B CN 106570284B
Authority
CN
China
Prior art keywords
thermal
collecting tube
model
dimensional
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
Application number
CN201610986195.0A
Other languages
English (en)
Other versions
CN106570284A (zh
Inventor
李印实
李陆
孙杰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian Jiaotong University
Original Assignee
Xian Jiaotong University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Xian Jiaotong University filed Critical Xian Jiaotong University
Priority to CN201610986195.0A priority Critical patent/CN106570284B/zh
Publication of CN106570284A publication Critical patent/CN106570284A/zh
Application granted granted Critical
Publication of CN106570284B publication Critical patent/CN106570284B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation

Abstract

本发明公开了一种聚焦式太阳能热发电系统的多维分析方法,并提出一种多维耦合模型。该多维耦合模型针对不同的问题采用不同的子模型,如聚光集热计算借助于蒙特卡洛光线追迹法,热工水力参数采用一维离散模型,集热管温度分布采用二维有限模型,热应力问题采用三维有限元模型,热力循环计算则采用零维模型。这种耦合方式可以实现精度与效率双优,同时也能获得局部详细信息与系统整体性能,该模型的利用具有极大的灵活性。可以根据获得的热工水力参数调节运行系统,可以用来预测各种不确定因素对电厂性能的影响,可以实时定位监测集热管的热应力应变情况。同时该模型方法也很容易应用于其他聚光形式的太阳能热发电系统。

Description

一种聚焦式太阳能热发电系统的多维分析方法
【技术领域】
本发明属于太阳能热利用技术领域,具体涉及一种聚焦式太阳能热发电系统的多维分析方法。
【背景技术】
随着新能源技术的不断提高,人类社会对传统能源的依赖将日趋减小,这既有助于缓解能源紧俏问题,也有助于改善环境污染问题,其中太阳能热发电技术越来越受到人们的青睐。太阳能热发电系统的稳定性,可靠性以及安全性严重制约其商业化发展。聚焦式太阳能热发电系统作为一个复杂的光-热-电转换系统,除了伴随着突出的热结构问题外,一系列不确定因素也对系统平稳性造成冲击,诸如雾霾,云层遮挡等。能够准确有效预测这些因素的影响,对系统的操作运行有积极作用。尽管各国学者在这一领域做了丰富的研究,但成果都过于单一化、分散化,割裂了局部与系统整体之间的联系,同时也缺乏有效的即时预测模型。构建一套精准而高效的光-热-电-结构耦合模型,既能及时监测热工水力参数,又能定位监测热应力-应变突出的管段,还能预测在不同干扰作用下系统性能,将会极大地坚定太阳能热电产业的投资信心,促进太阳能热电产业的发展。
【发明内容】
本发明的目的在于克服上述现有技术的缺点,提供一种聚焦式太阳能热发电系统的多维分析方法,能高效地实现太阳能到电能的转化,同时提供详细的热工水力参数,定向监测集热器热应力应变问题,预测系统性能。
为达到上述目的,本发明采用以下技术方案予以实现:
一种聚焦式太阳能热发电系统的多维分析方法,包括以下步骤:
1)针对太阳能热发电系统聚光器结构建立三维光学计算模型,采用蒙特卡洛光线追迹法计算集热管接受的辐射热流密度分布;
2)沿程对传热工质采用一维离散模型,用于求解热工水力参数;
3)采用有限容积法,对集热器建立二维离散模型,用于求解集热管温度分布;
4)采用有限元法,建立集热管三维离散模型,加载二维温度边界条件,求解集热管的应力应变问题;
5)热力循环采用零维模型,输入主蒸汽参数温度t,压力p和流量得到系统的输出电功W:
本发明进一步的改进在于:
所述步骤1)中,采用蒙特卡洛光线追迹法计算集热管接受的辐射热流密度分布的具体方法如下:
建立聚光器的三维几何模型;设定在抛物槽镜面、真空玻璃管和金属集热管表面可能发生的基本光学事件、初始化光子初始位置及入射方向;启动程序直至计算结束后统计光子在集热管壁面的分布并转化为相应的热流密度分布;基本光学事件包括反射、吸收、透射或折射。
所述步骤2)和3)中求解热工水力参数和集热管温度分布的具体方法如下:
获得集热管热流密度边界后,根据一维离散模型1D和二维离散模型2D建立集热管的热平衡关系,求解热工水力参数和集热管温度分布;其中一维离散模型1D相邻单元温度关系为:
集热管的二维传热控制方程为:
所述步骤4)中,获得集热管温度分布后,根据工质的温度与焓值以及集热管轴向温差的大小判断是否进行三维应力结构分析,即预判是否为集热管承受的应力或变形最大位置;
如果判定是最易出现安全问题的集热管部分,则调用三维离散模型;热结构问题基本控制方程包括三维平衡方程、本构方程和应变位移关系,柱坐标系下分别为:
求解热应力应变问题时,集热管二维温度分布需作为温度边界条件加载给三维模型,用于解决纯三维计算耗时耗内存的缺点。
与现有技术相比,本发明具有以下有益效果:
本发明能够定向监测热应力问题突出的地方,可以模拟各种干扰对系统性能的影响;同时各子模型可以以不同的耦合方式有针对性地解决不同层面的问题。采用该多维耦合模型,使得计算模型详略得当、灵活自如,既有效保留了必要的局部信息和整体性能,又能研究两者之间的相互作用,还实现了计算精度和计算效率的双优,为太阳能热发电系统集成化设计运行提供了高效可靠的分析结果。
【附图说明】
图1为本发明主要子模块模型;
图2为本多维耦合模型程序结构;
图3由MCRT求得的集热管热流密度分布;
图4集热管内外壁面平均温度分布;
图5为集热管应力及变形情况。
【具体实施方式】
下面结合附图对本发明做进一步详细描述:
参见图1,本发明针对太阳能热发电系统聚光器结构建立三维光学计算模型,采用蒙特卡洛光线追迹法(MCRT)计算集热管接受的辐射热流密度分布(见图1);
沿程对传热工质采用一维离散模型(1D),用于求解热工水力参数;采用有限容积法,对集热器建立二维(周向和轴向)离散模型(2D)用于求解集热管温度分布(见图1);
采用有限元法建立集热管三维离散模型(3D),加载[0006]中二维温度边界条件,求解集热管的应力应变问题(见图1);
热力循环采用零维模型(0D),输入[0006]主蒸汽参数温度t,压力p和流量求解系统的输出电功W;即
四个子模型协同运行才能实现目标[0003]。0-1D作为外循环负责热工水力及热电转化计算;1-2D模型共同求解热工水力参数和集热管的温度场;2D温度分布可以预判是否启用3D模型并作为温度边界条件求解热应力应变问题;
一旦获得集热管的热变形数据即可反馈至[0005]更新计算集热管热流密度分布,进一步根据0-1D模型研究其变形对系统性能的影响。
参见图2,该多维耦合模型算法结构为:
采用MCRT求解集热管辐射热流密度分布。该过程主要包括:建立聚光器的三维几何模型;设定在抛物槽镜面,真空玻璃管和金属管表面可能发生的基本光学事件(如反射,吸收,透射,折射);初始化光子初始位置及入射方向;启动程序直至计算结束后统计光子在集热管壁面的分布并转化为相应的热流密度分布(图3)。
获得集热管热流密度边界后,根据一维离散模型1D、二维离散模型2D建立集热管的热平衡关系求解热工水力参数和集热管温度分布。其中一维离散模型1D相邻单元温度关系为:
集热管的二维传热控制方程为:
求解的二维温度为集热管内外壁面温度分布,内外壁面平均温度分布如图4所示。
求解完毕后,一方面热工水力参数可以作为参考来调整电厂的运行,另一方面出口蒸汽参数可以直接传递给零维模型计算净发电量等(见公式(1))
另一方面,获得集热管温度分布后,根据工质的温度与焓值以及集热管轴向温差的大小判断是否进行三维应力结构分析,即预判断是否是集热管承受的应力或变形最大位置。
如果判定是最易出现安全问题的集热管部分,则调用三维离散模型3D。热结构问题基本控制方程包括三维平衡方程,本构方程和应变位移关系,柱坐标系下分别为:
求解热应力应变问题时,集热管二维温度分布需作为温度边界条件加载给三维离散模型3D,这样可以解决纯三维计算耗时耗内存的缺点。
图5为集热管的三维应力及变形情况。该分析结果用于判断集热管是否能够安全运行以及可能的工作寿命,为潜在的危险作预警。
同时集热管弯曲变形数据可以直接反馈至MCRT部分,用于重新计算热流密度分布情况,然后继续调用一维离散模型1D、二维离散模型2D和零维模型
进一步研究弯曲后集热性能的变化。
这里MCRT可以获得镜场的聚光集热性能以及详细的辐射热流密度分布,一维离散模型1D、二维离散模型2D可以获得详细的热工水力参数和集热管的温度分布,零维模型可以获得系统层面的性能参数,三维离散模型3D可以根据需要,定点监测集热管热结构情况,而各子模型的耦合,即本发明所提的多维耦合模型,完全能够灵活高效地解决聚光式太阳能热发电系统多层面的问题。
以上内容仅为说明本发明的技术思想,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本发明权利要求书的保护范围之内。

Claims (2)

1.一种聚焦式太阳能热发电系统的多维分析方法,其特征在于,包括以下步骤:
1)针对太阳能热发电系统聚光器结构建立三维光学计算模型,采用蒙特卡洛光线追迹法计算集热管接受的辐射热流密度分布;
2)沿程对传热工质采用一维离散模型,用于求解热工水力参数;
3)采用有限容积法,对集热器建立二维离散模型,用于求解集热管温度分布;求解热工水力参数和集热管温度分布的具体方法如下:
获得集热管热流密度边界后,根据一维离散模型1D和二维离散模型2D建立集热管的热平衡关系,求解热工水力参数和集热管温度分布;其中一维离散模型1D相邻单元温度关系为:
集热管的二维传热控制方程为:
4)采用有限元法,建立集热管三维离散模型,加载二维温度边界条件,求解集热管的应力应变问题;获得集热管温度分布后,根据工质的温度与焓值以及集热管轴向温差的大小判断是否进行三维应力结构分析,即预判是否为集热管承受的应力或变形最大位置;
如果判定是最易出现安全问题的集热管部分,则调用三维离散模型;热结构问题基本控制方程包括三维平衡方程、本构方程和应变位移关系,柱坐标系下分别为:
求解热应力应变问题时,集热管二维温度分布需作为温度边界条件加载给三维模型,用于解决纯三维计算耗时耗内存的缺点;
5)热力循环采用零维模型,输入主蒸汽参数温度t,压力p和流量得到系统的输出电功W:
2.根据权利要求1所述的聚焦式太阳能热发电系统的多维分析方法,其特征在于,所述步骤1)中,采用蒙特卡洛光线追迹法计算集热管接受的辐射热流密度分布的具体方法如下:
建立聚光器的三维几何模型;设定在抛物槽镜面、真空玻璃管和金属集热管表面可能发生的基本光学事件、初始化光子初始位置及入射方向;启动程序直至计算结束后统计光子在集热管壁面的分布并转化为相应的热流密度分布;基本光学事件包括反射、吸收、透射或折射。
CN201610986195.0A 2016-11-09 2016-11-09 一种聚焦式太阳能热发电系统的多维分析方法 Active CN106570284B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610986195.0A CN106570284B (zh) 2016-11-09 2016-11-09 一种聚焦式太阳能热发电系统的多维分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610986195.0A CN106570284B (zh) 2016-11-09 2016-11-09 一种聚焦式太阳能热发电系统的多维分析方法

Publications (2)

Publication Number Publication Date
CN106570284A CN106570284A (zh) 2017-04-19
CN106570284B true CN106570284B (zh) 2019-11-08

Family

ID=58540663

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610986195.0A Active CN106570284B (zh) 2016-11-09 2016-11-09 一种聚焦式太阳能热发电系统的多维分析方法

Country Status (1)

Country Link
CN (1) CN106570284B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108871731A (zh) * 2017-05-11 2018-11-23 天津滨海光热跟踪技术有限公司 集热器漏光光线捕捉测量方法
CN109376450B (zh) * 2018-11-09 2023-07-07 上海电气集团股份有限公司 一种太阳能吸热管温度场分析建模方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101278292A (zh) * 2004-12-14 2008-10-01 施卢默格控股有限公司 用于油藏模拟器中的耦合应力/液体流动的有限容积方法
CN102270251A (zh) * 2011-04-18 2011-12-07 西安交通大学 一种多阶多表面复杂太阳能聚焦集热系统设计方法
CN102609567A (zh) * 2012-01-13 2012-07-25 西安交通大学 一种太阳能聚焦集热系统光-热-流体耦合优化设计方法
CN103727687A (zh) * 2013-12-30 2014-04-16 北京工业大学 一种双熔盐太阳能热发电传热蓄热系统

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5005360A (en) * 1990-02-22 1991-04-09 Mcmurtry J A Solar energy system for generating electricity

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101278292A (zh) * 2004-12-14 2008-10-01 施卢默格控股有限公司 用于油藏模拟器中的耦合应力/液体流动的有限容积方法
CN102270251A (zh) * 2011-04-18 2011-12-07 西安交通大学 一种多阶多表面复杂太阳能聚焦集热系统设计方法
CN102609567A (zh) * 2012-01-13 2012-07-25 西安交通大学 一种太阳能聚焦集热系统光-热-流体耦合优化设计方法
CN103727687A (zh) * 2013-12-30 2014-04-16 北京工业大学 一种双熔盐太阳能热发电传热蓄热系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
聚焦式太阳能热发电中的蓄热技术及系统;尹辉斌 等;《热能动力工程》;20130131;第1-6,15页 *

Also Published As

Publication number Publication date
CN106570284A (zh) 2017-04-19

Similar Documents

Publication Publication Date Title
Qiu et al. Thermal performance analysis of a parabolic trough solar collector using supercritical CO2 as heat transfer fluid under non-uniform solar flux
Wu et al. Effect of cooling channel position on heat transfer characteristics and thermoelectric performance of air-cooled PV/T system
Das et al. 3D numerical study on estimating flow and performance parameters of solar updraft tower (SUT) plant: Impact of divergent angle of chimney, ambient temperature, solar flux and turbine efficiency
Kumar et al. Numerical investigation of natural convection heat loss in modified cavity receiver for fuzzy focal solar dish concentrator
CN106936627B (zh) 一种基于大数据分析挖掘的火电设备性能监测方法
CN106570284B (zh) 一种聚焦式太阳能热发电系统的多维分析方法
CN105973505B (zh) 一种太阳能腔式吸热器开口处热流密度测定方法
Qaisrani et al. Heat losses and thermal stresses of an external cylindrical water/steam solar tower receiver
CN105720917B (zh) 基于红外图像处理的光伏阵列热斑快速检测系统及方法
CN110442959B (zh) 太阳能碟式聚光系统聚光性能的预测与保持方法
CN106369838B (zh) 一种槽式聚光太阳能集热系统设计方法
CN108131848A (zh) 一种dsg槽式光热电站集热场工质参数预测方法
Abdunnabi et al. Experimental validation of forced circulation of solar water heating systems in TRNSYS
CN107947182A (zh) 一种动态潮流分析方法及动态潮流分析系统
Hu et al. Numerical simulation on convective thermal loss of a cavity receiver in a solar tower power plant
CN108446465A (zh) 通过工质分解在线测算火电机组厂用蒸汽量的方法
Kribus Thermal integral micro-generation systems for solar and conventional use
Wang et al. Numerical simulation of the heat flux distribution in a solar cavity receiver
Gobereit et al. Cost analysis of different operation strategies for falling particle receivers
Guo et al. The dynamic behavior of once-through direct steam generation parabolic trough solar collector row under moving shadow conditions
Zhou et al. Experimental and numerical evaluation of a lab-scale external solar receiver
Ozalp et al. Effect of cameralike aperture in quest for maintaining quasi-constant radiation inside a solar reactor
CN108416372B (zh) 基于多种传感器的风电场风速获取方法
Kerzmann et al. Flow rate optimization of a linear concentrating photovoltaic system
Henrique Rossa et al. Simulation of Energy Production in Grid-Connected Photovoltaic Systems From Measured and Calculated Data From Clear-Sky Radiation Model

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant