CN105468832B - 一种超临界水堆堆芯核热耦合迭代方法 - Google Patents
一种超临界水堆堆芯核热耦合迭代方法 Download PDFInfo
- Publication number
- CN105468832B CN105468832B CN201510801836.6A CN201510801836A CN105468832B CN 105468832 B CN105468832 B CN 105468832B CN 201510801836 A CN201510801836 A CN 201510801836A CN 105468832 B CN105468832 B CN 105468832B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- power density
- density distribution
- neutronics calculation
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Monitoring And Testing Of Nuclear Reactors (AREA)
Abstract
超临界水堆堆芯核热耦合迭代方法,执行完m次中子学计算后执行一次热工水力计算,m为设定的中子学计算迭代次数,m大于5;在执行中子学计算中,引入自适应松弛因子ω并利用功率密度分布迭代函数对第2次至第m次执行中子学计算得到的功率密度分布进行调整,调整后功率密度分布用于下一次执行中子学计算中或热工水力计算中。本发明解决了SCWR堆芯核热耦合计算中迭代计算难以收敛的问题,提高核热耦合计算的计算效率和数值稳定性。
Description
技术领域
本发明涉及核反应堆设计技术领域,具体地,涉及一种超临界水堆堆芯核热耦合迭代方法。
背景技术
超临界水堆(SCWR)的运行工况在水的热力学临界点(374℃,22.1MPa)以上,具备较高的***热效率且反应堆***简化、紧凑。但是,超临界水的密度、比热以及热导率在超临界区域变化剧烈,在通常定义的运行压力25MPa、堆芯冷却剂进口/出口温度280℃/500℃的超临界水堆中,冷却剂的密度由堆芯入口的约777kg/m3变化到堆芯出口的约90kg/m3。剧烈变化的水密度显著影响堆芯内中子慢化能力,慢化能力变化引起的功率密度分布变化反过来又将显著影响堆芯内水密度分布。这种物理-热工耦合特性显著强于传统的压水堆,要求在堆芯稳态性能分析及方案设计时必须考虑核热耦合。
SCWR堆芯中,冷却剂密度、慢化剂密度分布不均匀且变化剧烈,功率密度分布的微小变化可能引起冷却剂密度、慢化剂密度的剧烈变化,进而显著改变中子截面参数,反过来使功率密度分布也产生剧烈变化。上述强烈核热耦合特性,使得SCWR堆芯的核热耦合计算面临迭代收敛性、数值稳定性以及计算效率的挑战。
现有核热耦合迭代技术中,耦合迭代方法通常采用解耦迭代方法和逐次迭代方法。解耦迭代方法要求每个迭代步中,所有计算完全收敛后才开始下一步的迭代计算。对于一个给定迭代步,基于当前热工-水力参数和中子截面参数的预估值,计算至中子裂变源完全收敛。解耦迭代方法的计算效率较低,且对于迭代计算参量变化剧烈的分析对象,如冷却剂密度变化剧烈的SCWR堆芯,甚至会导致迭代计算不收敛。逐次迭代方法假设一个功率密度分布,通过迭代公式不断进行迭代,直至相邻两次的功率密度分布最大相对偏差小于收敛精度要求。对于迭代计算参量变化剧烈的SCWR堆芯,相邻两次的功率密度分布容易出现功率振荡,导致无法收敛。特别是对于多流程SCWR堆芯,在堆芯顶部和底部可能分别形成功率峰,这将进一步导致堆芯内功率密度分布难以收敛。
发明内容
本发明的目的就在于克服上述现有技术的缺点和不足,提供一种超临界水堆堆芯核热耦合迭代方法,该方法解决了超临界水堆堆芯核热耦合计算过程中迭代计算难以收敛的问题,并提高核热耦合计算的计算效率和数值稳定性。
本发明解决上述问题所采用的技术方案是:
一种超临界水堆堆芯核热耦合迭代方法,执行完m次中子学计算后执行一次热工水力计算,m为设定的中子学计算迭代次数,m大于5。本技术方案中,在执行每步热工水力计算前并不要求外迭代完全收敛才执行热工水力计算,也不要求每执行一次外迭代就进行一次热工水力计算;核热耦合计算的效率得到显著提高。
作为本发明的进一步改进,在执行中子学计算中,引入自适应松弛因子ω并利用功率密度分布迭代函数对第2次至第m次执行中子学计算得到的功率密度分布进行调整,得到的调整后功率密度分布用于下一次执行中子学计算中或执行热工水力计算中,具体包括以下步骤:
S1、设定中子学计算迭代次数m,构建功率密度分布迭代函数,初始化执行中子学计算次数n;
S2、执行第一次中子学计算,并统计中子学计算次数后跳转到步骤S3;
S3、利用执行第一次中子学计算得到的功率密度分布执行第二次中子学计算,并统计中子学计算次数后跳转到步骤S4;
S4、计算最近一次执行中子学计算得到的功率密度分布的相对偏差,并根据该相对偏差选取自适应松弛因子ω的取值;
S5、将步骤S4中选取的自适应松弛因子ω的取值和最近一次执行中子学计算得到的功率密度分布代入功率密度分布迭代函数计算出调整后功率密度分布,判断中子学计算次数是否等于中子学计算迭代次数,是则将中子学计算次数设为初始值后跳转到步骤S7,否则跳转到步骤S6;
S6、利用步骤S5中的调整后功率密度分布再一次执行中子学计算,统计中子学计算次数,跳转到步骤S4;
S7、利用步骤S5中的调整后功率密度分布执行热工水力计算,判断执行热工水力计算后得到的功率密度分布是否收敛,是则结束计算,否则跳转到步骤S2。本技术方案中,在功率密度分布迭代计算过程中引入松弛因子对功率密度分布进行调整,可以避免由于相邻两次迭代计算结果差值过大而引起迭代过程发散,自适应松弛因子随着迭代的收敛情况而自行调整取值,在保证迭代收敛性和数值稳定性的前提下,尽可能提高计算效率。
进一步,所述功率密度分布迭代函数为:
式(1)中,n为中子学计算次数;ωn为对第n次执行中子学计算得到的功率密度分布进行调整时使用的自适应松弛因子ω的取值;Pn-1为第n-1次执行中子学计算得到的功率密度分布,Pn为第n次执行中子学计算得到的功率密度分布,P′n-1为第n-1次执行中子学计算得到的功率密度分布的调整后功率密度分布,P′n为采用ωn对Pn进行调整得到的调整后功率密度分布。
进一步,步骤S4中所述选取自适应松弛因子ω的取值方法为:
式中,εp n为第n次执行中子学计算得到的功率密度分布的相对偏差,其中:
式(3)中,Pn-1为第n-1次执行中子学计算得到的功率密度分布,P′n-1为第n-1次执行中子学计算得到的功率密度分布的调整后功率密度分布,Pn为第n次执行中子学计算得到的功率密度分布。
进一步,所述中子学计算迭代次数m为20-40。
优选的,所述中子学计算迭代次数m为30。
进一步,所述步骤S6中判断执行热工水力计算后得到的功率密度分布是否收敛的具体方法为:判断执行热工水力计算后的功率密度分布和执行热工水力计算前的功率密度分布是否满足公式(2),满足则收敛,不满足不收敛:
式(4)中,f(q′m)为执行热工水力计算后的功率密度分布,q'm为第m次执行中子学计算得到的功率密度分布的调整后功率密度分布,δ为设定的收敛值。
综上,本发明的有益效果是:
1、本发明通过采用改进的外迭代耦合方法和匹配的自适应松弛因子,有效解决了SCWR堆芯核热耦合计算时的迭代收敛性问题,并确保了耦合计算的数值稳定性;
2、本发明将SCWR堆芯核热耦合计算的计算效率提高了三倍以上。
附图说明
图1是本发明的超临界水堆堆芯核热耦合迭代方法;
图2是超临界水堆堆芯核热耦合计算流程。
具体实施方式
本发明通过采用改进的外迭代耦合方法和匹配的自适应松弛因子,提高核热耦合计算的迭代收敛性、计算效率和数值稳定性。其主要思想是执行完设定次数(m次)的中子学计算后执行一次热工水力计算,m为设定的中子学计算迭代次数,m大于5。在执行中子学计算中,引入自适应松弛因子ω并利用功率密度分布迭代函数对每轮第2次至第m次执行中子学计算得到的功率密度分布进行调整得到调整后功率密度分布,该调整后功率密度分布用于下一次执行中子学计算中或执行热工水力计算中,保证迭代收敛性和效率。下面结合实施例及附图,对本发明作进一步地的详细说明,但本发明的实施方式不限于此。
实施例1:
如图1所示,一种超临界水堆堆芯核热耦合迭代方法,包括以下的步骤S1-S7:
S1、设定中子学计算迭代次数m,构建功率密度分布迭代函数;初始化执行中子学计算次数n,本实施例中初始化n时将n的值设定为0;
S2、执行第一次中子学计算,并统计中子学计算次数后跳转到步骤S3;本实施例中,统计中子学计算次数的方法为:使中子学计算次数自加1,即把n+1的值赋给n,本步骤中,中子学计算次数自加1后n的值为1;
S3、利用执行第一次中子学计算得到的功率密度分布P1执行第二次中子学计算,并统计中子学计算次数后(此时n的值为2)跳转到步骤S4;本步骤中,统计中子学计算次数的方法为:使中子学计算次数自加1,即将n+1的值赋给n;
S4、计算最近一次执行中子学计算得到的功率密度分布的相对偏差,并根据该相对偏差选取自适应松弛因子ω的取值;其中,各次执行中子学计算得到的功率密度分布的相对偏差εp n的计算公式如下:
式(3)中,Pn-1为第n-1次执行中子学计算得到的功率密度分布,Pn为第n次执行中子学计算得到的功率密度分布,P′n-1为第n-1次执行中子学计算得到的功率密度分布的调整后功率密度分布,实际上,P′n-1也是采用自适应松弛因子ω对Pn-1进行调整得到的调整后功率密度分布。在执行第n次执行中子学计算后、执行第n+1次执行中子学计算前的最近一次执行中子学计算即第n次执行中子学计算,最近一次执行中子学计算得到的功率密度分布即Pn。本发明的调整后功率密度分布是指采用相应的自适应松弛因子的取值对执行中子学计算得到的功率密度分布进行调整后的功率密度分布。
本实施例中将每执行m次中子学计算后执行一次热工水力计算称为一轮热耦合计算;则在每轮热耦合计算中,第一次执行中子学计算后的功率密度分布不进行调整,因此不用计算相对偏差;第二次执行中子学计算后的功率密度分布的相对偏差εp 2为第二次执行中子学计算得到的功率密度P2与第一次执行中子学计算得到的功率密度P1之差的绝对值除以P1;第三次及之后的执行中子学计算后的功率密度分布的相对偏差εp n(n=3,4,...,m)为当次执行中子学计算得到的功率密度Pn与上次执行中子学计算得到的功率密度分布的调整后功率密度分布P′n-1之差的绝对值除以上次执行中子学计算得到的功率密度分布的调整后功率密度分布P′n-1。
选取自适应松弛因子ω的取值公式如下:
式中,ωn即为对第n次执行中子学计算得到的功率密度分布进行调整时使用的自适应松弛因子ω的取值,εp n为第n次执行中子学计算得到的功率密度分布的相对偏差。
S5、将步骤S4中选取的自适应松弛因子ω的取值ωn和最近一次执行中子学计算得到的功率密度分布代入功率密度分布迭代函数计算出调整后功率密度分布,判断中子学计算次数是否等于中子学计算迭代次数,是则将中子学计算次数设为初始值后跳转到步骤S7,否则跳转到步骤S6;所述功率密度分布迭代函数主要用于对第2至第m次执行中子学计算得到的功率密度分布进行调整,在第2至第m次执行中子学计算中,每次执行完中子学计算后对得到的功率密度分布进行一次调整,第2至第m-1次的调整后的功率密度分布用于下一次的执行中子学计算中和下一次的功率密度分布的偏差计算中,第m次的调整后的功率密度分布用于执行热工水力计算中,功率密度分布迭代函数公式如下:
式(1)中,n为中子学计算次数;Pn-1为第n-1次执行中子学计算得到的功率密度分布,Pn为第n次执行中子学计算得到的功率密度分布,P′n-1为第n-1次执行中子学计算得到的功率密度分布的调整后功率密度分布,P′n为采用ωn对Pn进行调整得到的调整后功率密度分布,将ωn、Pn-1、Pn、P′n-1代入式(1)就能计算出Pn的调整后功率密度分布P′n。
本实施例中,将中子学计算次数设为初始值即进行清零操作也即将n设置为0,统计中子学计算次数即把n+1的值赋给n,新的中子学计算次数等于上一次的中子学计算次数加1;本实施例中,执行一次中子学计算即进行一次外迭代。在其他实施方式中,初始化是也可以将n设定为其他值,例如但不限于等于m,这种情况下本步骤中,判断执行中子学计算次数是否达到中子学计算迭代次数时判断n是否等于0,将执行中子学计算次数设定为初始值时将n设置为m;本步骤、步骤S2和S3中的统计中子学计算次数的方式是:使中子学计算次数自减1。
S6、利用步骤S5中的调整后功率密度分布再一次执行中子学计算,统计中子学计算次数,跳转到步骤S4。
S7、利用步骤S5中的调整后功率密度分布(实际上此时步骤S5中的调整后功率密度分布就是对第m次执行中子学计算得到的功率密度分布Pm进行调整后得到的调整后功率密度分布P′m)执行热工水力计算,判断执行热工水力计算后得到的功率密度分布是否收敛,是则结束计算,否则跳转到步骤S2。
本实施例中,所述中子学计算迭代次数m为30次,实际应用中,中子学计算迭代次数m可以在20-40次范围内调整。但中子学计算迭代次数m为30次时,相比传统的解耦迭代方法耦合计算耗时减少了2/3,显著提高了计算效率。
步骤S4中,式(2)表明当相邻两次外迭代(中子学计算)的功率密度分布相对偏差大于0.1时,自适应松弛因子取值为0.2;当相邻两次外迭代的功率密度分布相对偏差小于或等于0.1时,自适应松弛因子取值为0.9和“0.2-0.3×lg(εp n)二者间的较小值。实验表明,采用匹配的自适应松弛因子,一方面兼顾了计算效率,另一方面有效解决了堆芯核热耦合计算时的迭代收敛性问题,并确保了耦合计算的数值稳定性。
步骤S7中,判断执行热工水力计算后得到的功率密度分布是否收敛的具体方法为:判断执行热工水力计算后的功率密度分布和执行热工水力计算前的功率密度分布是否满足公式(4),满足则收敛,不满足不收敛:
式中,f(q'm)为执行热工水力计算后的功率密度,q'm为第m次执行中子学计算得到的功率密度分布的调整后功率密度分布,δ为预先设定的收敛值。本实施例中,第m次执行中子学计算得到的功率密度分布的调整后功率密度分布P′m也即对执行第30次中子学计算后的功率密度q30进行调整得到的调整后功率密度分布P′30。
前述执行堆芯中子学计算作为整个超临界水堆堆芯核热耦合迭代方法中的外迭代,在超临界水堆堆芯核热耦合计算过程中,若采用传统的解耦迭代方法,需要等到外迭代完全收敛后,才执行热工水力计算,则每次热工水力计算前需至少执行100次堆芯中子学计算,计算效率较低且当冷却剂密度变化剧烈时,迭代计算不收敛。若直接采用外迭代耦合方法,需要每执行一次外迭代,就进行一次热工水力计算,由于执行一次热工水力计算耗时较长,整个迭代过程中热工水力计算部分耗时过多,导致整体计算效率较低。
本实施例中,一种超临界水堆堆芯核热耦合迭代方法采用改进的外迭代耦合方法,在执行每步热工水力计算前并不要求外迭代完全收敛才执行热工水力计算,也不要求每执行一次外迭代就进行一次热工水力计算。本实施例中的改进的外迭代耦合方法是每执行设定次数(m次)的外迭代计算后执行一次热工水力计算。采用这种改进的外迭代耦合方法,核热耦合计算的收敛性得到显著提高,并未出现功率或密度不收敛情形,同时核热耦合计算的效率得到显著提高。本实施例中,每执行30次外迭代计算后执行一次热工水力计算,相比传统的解耦迭代方法耦合计算耗时减少了2/3,显著提高了计算效率。
在核热耦合计算过程中,若采用传统的逐次迭代方法,由于SCWR堆芯冷却剂密度、慢化剂密度分布不均匀且变化剧烈,功率密度分布的微小变化可能引起冷却剂密度、慢化剂密度的剧烈变化,进而显著改变中子截面参数,反过来使功率密度分布也产生剧烈变化,将极易出现功率振荡,导致无法收敛。本实施例中,还在功率密度分布迭代计算过程中引入松弛因子,可以避免由于相邻两次迭代计算结果差值过大而引起迭代过程发散。若采用固定的松弛因子,当松弛因子较大时,虽然计算耗时相对较少,但不易实现迭代计算收敛,当松弛因子较小时,虽然容易实现迭代计算收敛的同时,但将花费更多的计算时间。一般而言,针对特定问题,如SCWR堆芯核热耦合计算,存在一个最合适的松弛因子。为此,本实施例中提出了匹配的自适应松弛因子。匹配的自适应松弛因子随着迭代的收敛情况而自行调整,在保证迭代收敛性和数值稳定性的前提下,尽可能提高计算效率。
综上,上述超临界水堆堆芯核热耦合迭代方法,通过采用改进的外迭代耦合方法和匹配的自适应松弛因子,有效解决了SCWR堆芯核热耦合计算时的迭代收敛性问题,并确保了耦合计算的数值稳定性,同时将SCWR堆芯核热耦合计算的计算效率提高了三倍以上。
本实施例中的执行中子学计算和执行热工水力计算均采用现有技术中的常用方法进行计算,一般由相应的计算机程序执行计算,本实施例中不再赘述其具体方法和流程。
本实施例中初始功率密度分布、执行中子学计算后的功率密度分布、执行热工水力学计算后的功率密度分布的计算方法为现有技术,本实施例中不再赘述,本发明的目的是通过这些功率密度分布来选取自适应松弛因子ω取值。
实施例2:
本实施例中提供一种基于实施例1或实施例2中的迭代方法的超临界水堆堆芯核热耦合计算方法,如图2所示,包括以下步骤:
(一)、计算组件不同的燃耗工况;
(二)、执行截面拟合处理;
(三)、初始控制棒棒位和水密度分布;
(四)、调整控制棒位置;
(五)、堆芯核热耦合迭代,直到执行热工水力计算后得到的功率密度分布收敛;
(六)、判断是否到临界状态,是则跳转到步骤(七),否则跳转到步骤(四);
(七)、根据最大功率密度分布计算最大包壳壁面温度;
(八)、执行燃耗计算;
(九)、判断燃料是否到达寿期末,是则结束计算,否则跳转到步骤(十);
(十)、燃耗步加1,跳转到步骤(四)。
该计算方法中,步骤五即实施例1中的超临界水堆堆芯核热耦合迭代方法,其余步骤全部采用现有技术中的处理方法,本实施例中不再赘述各步骤的处理方法。
以上仅是本发明的优选实施方式,本发明的保护范围并不仅局限于上述实施例,凡属于本发明思路下的技术方案均属于本发明的保护范围。应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理前提下的若干改进和润饰,应视为本发明的保护范围。
Claims (5)
1.一种超临界水堆堆芯核热耦合迭代方法,其特征在于,执行完m次中子学计算后执行一次热工水力计算,m为设定的中子学计算迭代次数,m大于5;
在执行中子学计算中,引入自适应松弛因子ω并利用功率密度分布迭代函数对第2次至第m次执行中子学计算得到的功率密度分布进行调整,得到的调整后功率密度分布用于下一次执行中子学计算或执行热工水力计算中,所述超临界水堆堆芯核热耦合迭代方法具体包括以下步骤:
S1、设定中子学计算迭代次数m,构建功率密度分布迭代函数,初始化执行中子学计算次数n;
S2、执行第一次中子学计算,并统计中子学计算次数后跳转到步骤S3;
S3、利用执行第一次中子学计算得到的功率密度分布执行第二次中子学计算,并统计中子学计算次数后跳转到步骤S4;
S4、计算最近一次执行中子学计算得到的功率密度分布的相对偏差,并根据该相对偏差选取自适应松弛因子ω的取值;
S5、将步骤S4中选取的自适应松弛因子ω的取值和最近一次执行中子学计算得到的功率密度分布代入功率密度分布迭代函数计算出调整后功率密度分布,判断中子学计算次数是否等于中子学计算迭代次数,是则将中子学计算次数设为初始值后跳转到步骤S7,否则跳转到步骤S6;
S6、利用步骤S5中的调整后功率密度分布再一次执行中子学计算,统计中子学计算次数,跳转到步骤S4;
S7、利用步骤S5中的调整后功率密度分布执行热工水力计算,判断执行热工水力计算后得到的功率密度分布是否收敛,是则结束计算,否则跳转到步骤S2;
步骤S4中所述选取自适应松弛因子ω的取值方法为:
<mrow>
<msub>
<mi>&omega;</mi>
<mi>n</mi>
</msub>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>min</mi>
<mrow>
<mo>(</mo>
<mn>0.9</mn>
<mo>,</mo>
<mn>0.2</mn>
<mo>-</mo>
<mn>0.3</mn>
<mo>&times;</mo>
<mi>lg</mi>
<mo>(</mo>
<mrow>
<msup>
<msub>
<mi>&epsiv;</mi>
<mi>p</mi>
</msub>
<mi>n</mi>
</msup>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msup>
<msub>
<mi>&epsiv;</mi>
<mi>p</mi>
</msub>
<mi>n</mi>
</msup>
<mo>&le;</mo>
<mn>0.1</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0.2</mn>
<mo>,</mo>
<msup>
<msub>
<mi>&epsiv;</mi>
<mi>p</mi>
</msub>
<mi>n</mi>
</msup>
<mo>></mo>
<mn>0.1</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mi>n</mi>
<mo>=</mo>
<mn>2</mn>
<mo>,</mo>
<mn>4</mn>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<mi>m</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
式中,εp n为第n次执行中子学计算得到的功率密度分布的相对偏差,其中:
<mrow>
<msup>
<msub>
<mi>&xi;</mi>
<mi>p</mi>
</msub>
<mi>n</mi>
</msup>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mo>|</mo>
<mrow>
<msub>
<mi>P</mi>
<mi>n</mi>
</msub>
<mo>-</mo>
<msub>
<msup>
<mi>P</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
<mo>|</mo>
</mrow>
<mrow>
<msub>
<msup>
<mi>P</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>,</mo>
<mi>n</mi>
<mo>=</mo>
<mn>3</mn>
<mo>,</mo>
<mn>4</mn>
<mo>,</mo>
<mn>5</mn>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<mi>m</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mo>|</mo>
<mrow>
<msub>
<mi>P</mi>
<mi>n</mi>
</msub>
<mo>-</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
<mo>|</mo>
</mrow>
<msub>
<mi>P</mi>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mfrac>
<mo>,</mo>
<mi>n</mi>
<mo>=</mo>
<mn>2</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
式(3)中,Pn-1为第n-1次执行中子学计算得到的功率密度分布,P′n-1为第n-1次执行中子学计算得到的功率密度分布的调整后功率密度分布,Pn为第n次执行中子学计算得到的功率密度分布。
2.根据权利要求1所述的一种超临界水堆堆芯核热耦合迭代方法,其特征在于,所述功率密度分布迭代函数为:
<mrow>
<msubsup>
<mi>P</mi>
<mi>n</mi>
<mo>&prime;</mo>
</msubsup>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&omega;</mi>
<mi>n</mi>
</msub>
<mo>&times;</mo>
<msub>
<mi>P</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>&omega;</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>&prime;</mo>
</msubsup>
<mo>,</mo>
<mi>n</mi>
<mo>=</mo>
<mn>3</mn>
<mo>,</mo>
<mn>4</mn>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<mi>m</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&omega;</mi>
<mi>n</mi>
</msub>
<mo>&times;</mo>
<msub>
<mi>P</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>&omega;</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>,</mo>
<mi>n</mi>
<mo>=</mo>
<mn>2</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
式(1)中,n为中子学计算次数;ωn为对第n次执行中子学计算得到的功率密度分布进行调整时使用的自适应松弛因子ω的取值;Pn-1为第n-1次执行中子学计算得到的功率密度分布,Pn为第n次执行中子学计算得到的功率密度分布,P′n-1为第n-1次执行中子学计算得到的功率密度分布的调整后功率密度分布,P′n为采用ωn对Pn进行调整得到的调整后功率密度分布。
3.根据权利要求1或2所述的一种超临界水堆堆芯核热耦合迭代方法,其特征在于,所述中子学计算迭代次数m为20-40。
4.根据权利要求1或2所述的一种超临界水堆堆芯核热耦合迭代方法,其特征在于,所述中子学计算迭代次数m为30。
5.根据权利要求1或2所述的一种超临界水堆堆芯核热耦合迭代方法,其特征在于,所述步骤S6中判断执行热工水力计算后得到的功率密度分布是否收敛的具体方法为:判断执行热工水力计算后的功率密度分布和执行热工水力计算前的功率密度分布是否满足公式(4),满足则收敛,不满足不收敛:
<mrow>
<mfrac>
<mrow>
<mo>|</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<msub>
<msup>
<mi>q</mi>
<mo>&prime;</mo>
</msup>
<mi>m</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<msup>
<mi>q</mi>
<mo>&prime;</mo>
</msup>
<mi>m</mi>
</msub>
<mo>|</mo>
</mrow>
<mrow>
<msub>
<msup>
<mi>q</mi>
<mo>&prime;</mo>
</msup>
<mi>m</mi>
</msub>
</mrow>
</mfrac>
<mo>&le;</mo>
<mi>&delta;</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
式中,f(q'm)为执行热工水力计算后的功率密度分布,q'm为第m次执行中子学计算得到的功率密度分布的调整后功率密度分布,δ为设定的收敛值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510801836.6A CN105468832B (zh) | 2015-11-19 | 2015-11-19 | 一种超临界水堆堆芯核热耦合迭代方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510801836.6A CN105468832B (zh) | 2015-11-19 | 2015-11-19 | 一种超临界水堆堆芯核热耦合迭代方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105468832A CN105468832A (zh) | 2016-04-06 |
CN105468832B true CN105468832B (zh) | 2018-03-30 |
Family
ID=55606528
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510801836.6A Active CN105468832B (zh) | 2015-11-19 | 2015-11-19 | 一种超临界水堆堆芯核热耦合迭代方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105468832B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108875213B (zh) * | 2018-06-19 | 2022-04-12 | 哈尔滨工程大学 | 一种反应堆堆芯热工水力多尺度分析方法 |
CN112052579A (zh) * | 2020-08-28 | 2020-12-08 | 上海核工程研究设计院有限公司 | 一种基于浮动网格的核-热-力多物理耦合计算方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4711753A (en) * | 1986-03-19 | 1987-12-08 | Westinghouse Electric Corp. | Calibration of a nuclear reactor core parameter predictor |
CN103902784A (zh) * | 2014-04-11 | 2014-07-02 | 华北电力大学 | 用于超临界水堆瞬态核热耦合的安全分析计算装置 |
-
2015
- 2015-11-19 CN CN201510801836.6A patent/CN105468832B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4711753A (en) * | 1986-03-19 | 1987-12-08 | Westinghouse Electric Corp. | Calibration of a nuclear reactor core parameter predictor |
CN103902784A (zh) * | 2014-04-11 | 2014-07-02 | 华北电力大学 | 用于超临界水堆瞬态核热耦合的安全分析计算装置 |
Non-Patent Citations (3)
Title |
---|
SCWR 堆芯三维瞬态物理-热工水力耦合程序开发;王连杰等;《核动力工程》;20141231;第35卷(第S2期);186-189 * |
基于不动点理论的核热耦合计算收敛特性研究;刘仕倡等;《核动力工程》;20140630;第35卷(第3期);152-157 * |
超临界水堆反应堆物理-热工水力耦合程序超临界水堆反应堆物理-热工水力耦合程序***MCATHAS的开发;安萍 等;《核动力工程》;20101231;第31卷(第6期);52-74 * |
Also Published As
Publication number | Publication date |
---|---|
CN105468832A (zh) | 2016-04-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107122546A (zh) | 一种压水堆稳态计算的多物理耦合方法 | |
CN105653869B (zh) | 一种超临界水堆堆芯三维瞬态性能分析方法 | |
CN103177154B (zh) | 一种获得核燃料组件共振参数的方法 | |
Ma et al. | RMC/CTF multiphysics solutions to VERA core physics benchmark problem 9 | |
Mao et al. | Topological susceptibility to the one-loop order in chiral perturbation theory | |
CN103985058B (zh) | 一种基于改进多中心校正内点法的可用输电能力计算方法 | |
CN103020468B (zh) | 获得核反应堆热耦合的方法 | |
CN110427588B (zh) | 用于堆芯功率在线重构的本征正交分解样本选取方法 | |
CN105468832B (zh) | 一种超临界水堆堆芯核热耦合迭代方法 | |
CN104573364B (zh) | 一种乏燃料水池及贮存格架热工水力耦合计算方法 | |
Zhu et al. | The Implementation and Analysis of the MOC and CMFD Adjoint Capabilities in the 2D-1D Code MPACT | |
Lee et al. | Development status of Monte Carlo code at UNIST | |
Shen et al. | A regularized Newton method for degenerate unconstrained optimization problems | |
CN107103189B (zh) | 一种进行反应堆临界曲率搜索的方法 | |
Yang et al. | Core design study on CANDU-SCWR with 3D neutronics/thermal-hydraulics coupling | |
Wade et al. | Optimal control of nuclear reactor depletion | |
Gill | Inline thermal and xenon feedback iterations in Monte Carlo reactor calculations | |
García et al. | A subchannel coarsening method for Serpent2-SUBCHANFLOW applied to a full-core VVER problem | |
CN105389424A (zh) | 一种超临界水堆组件截面拟合方法 | |
JP2020046343A (ja) | 原子炉評価装置、原子炉評価装置方法及び原子炉評価プログラム | |
Kim et al. | Depletion analysis of fuel assembly with the IDTMC method | |
Pierre et al. | Analysis of the Superphenix start-up tests with APOLLO-3: From zero power isothermal conditions to dynamic power transient analysis | |
Matsuoka et al. | Radial electric field formation including electron radial drift for a core electron-root confinement (CERC) plasma in LHD | |
CN113792500B (zh) | 一种铅铋堆***热耦合计算方法 | |
Choi et al. | Improvement of resonance interference treatment in STREAM |
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 |