CN106156460B - 获得核燃料棒内部燃料芯块温度分布的方法 - Google Patents

获得核燃料棒内部燃料芯块温度分布的方法 Download PDF

Info

Publication number
CN106156460B
CN106156460B CN201510142261.1A CN201510142261A CN106156460B CN 106156460 B CN106156460 B CN 106156460B CN 201510142261 A CN201510142261 A CN 201510142261A CN 106156460 B CN106156460 B CN 106156460B
Authority
CN
China
Prior art keywords
fuel
temperature
fuel rod
pellets
theoretical
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
CN201510142261.1A
Other languages
English (en)
Other versions
CN106156460A (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.)
State Power Investment Group Science and Technology Research Institute Co Ltd
Original Assignee
National Nuclear Beijing Science And Technology Research Institute Co Ltd
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 National Nuclear Beijing Science And Technology Research Institute Co Ltd filed Critical National Nuclear Beijing Science And Technology Research Institute Co Ltd
Priority to CN201510142261.1A priority Critical patent/CN106156460B/zh
Publication of CN106156460A publication Critical patent/CN106156460A/zh
Application granted granted Critical
Publication of CN106156460B publication Critical patent/CN106156460B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Monitoring And Testing Of Nuclear Reactors (AREA)

Abstract

本发明提供了一种获得核燃料棒内部燃料芯块温度分布的方法。该方法包括：获得燃料棒有效温度TFdff和燃料芯块理论温度分布
Figure DDA0000689569770000011
的函数表达式；获得燃料棒燃耗深度BU；计算燃料棒包壳平均温度
Figure DDA0000689569770000012
由燃料棒燃耗深度BU和燃料棒功率p，计算燃料棒有效温度TFdff；由燃料棒包壳平均温度
Figure DDA0000689569770000013
和燃料棒有效温度TFdff，计算燃料芯块的平均温度
Figure DDA0000689569770000014
由燃料棒燃耗深度BU和燃料棒功率p，计算燃料芯块理论平均温度
Figure DDA0000689569770000015
利用燃料芯块平均温度
Figure DDA0000689569770000016
和燃料芯块理论平均温度
Figure DDA0000689569770000017
对燃料芯块理论温度分布进行校正，得到校正后的燃料芯块温度分布Tpellet(r)。本发明简化了核燃料棒内部燃料芯块温度分布的计算流程，计算速度更快，使得工程人员对燃料芯块温度的实时监控成为可能。

Description

获得核燃料棒内部燃料芯块温度分布的方法
技术领域
本发明涉及核技术领域，尤其涉及一种获得核燃料棒内部燃料芯块温度分布的方法。
背景技术
反应堆以装载的核燃料产生热量，将热量带出后传递给汽轮机，推动汽轮机发电。反应堆发热的主要原件是核燃料，核燃料一般被制作成很小的柱状体(俗称芯块)，填充在燃料棒包壳内，并将多组燃料棒捆绑成束状，组成核燃料组件。工程实践中，人们非常关注核燃料棒的温度参数，因为反应堆运行过程中，燃料棒的温度很高，很可能导致燃料棒融化，从而产生事故，威胁核电厂的安全。目前传统的做法是计算出燃料棒的平均功率，通过一个简单的折算，得到燃料棒的平均温度，这种处理方式简单，但缺乏精度，并且往往对燃料棒的温度估计过低，从而造成安全隐患。
在现代反应堆物理计算中，先进节块法可以高速准确地求解三维中子扩散方程，得到了工业界广泛的应用。然而有的节块法采用了横向积分技术对空间进行近似，这类节块法(例如节块展开法NEM、解析节块法ANM等)只能直接得到组件的平均功率，而在核设计和安全分析中往往需要提供组件内燃料棒的三维功率分布。因此需要根据节块法提供的节块平均量等信息对节块内的燃料棒进行精细功率重构。
精细功率重构方法可分为两类，一类是嵌入法(Embedded Method)，即利用均匀化介质组件内的功率分布给出非均匀介质组件计算所需的边界条件，对感兴趣的组件进行非均匀介质组件输运计算，得出组件内功率分布。这一类方法具有较高的精度，但是计算复杂，计算时间很长，而且通常只限于两群模型。另一类是调制法(Modulation Method)，即事先算出单个组件内部的非均匀功率分布(使用简单边界条件，如反射边界条件)作为形状函数，由均匀化介质组件内功率分布乘以该形状函数就得出非均匀介质组件内的功率分布。这一类方法的精度与嵌入法相近，但是计算效率更高，因此几乎程序都采用了调制法进行精细功率重构。
调制法也根据通量展开式的形式不同而不同，但是原理基本一致。在均匀节块内，将第g群二维中子通量密度展开成多项式Pi(x，y)或/和双曲(或三角)函数Fj(x，y)的形式：
Figure BDA0000689569750000021
其中：ai和bi为展开系数。
然后根据展开系数的个数，确定N个约束条件。根据这些约束条件求出的展开式系数ai和bi，进而可以求得均匀节块内的精细通量分布
Figure BDA0000689569750000022
在均匀的精细中子通量密度基础上，乘以事先由组件计算程序计算得到的组件形状因子fg(x，y)即可获得非均匀节块内的精细通量分布
Figure BDA0000689569750000023
Figure BDA0000689569750000024
最后根据式(3)获得非均匀节块内的精细功率分布phet(x，y)。
Figure BDA0000689569750000025
根据传热模型，可以由外至内计算出燃料棒的温度。图1为燃料棒横截面上温度场分区的示意图。如图1所示，燃料棒film层外为冷却剂区，认为film层外层温度等于冷却剂温度Tmod。
(1)Film层温差ΔTfilm
Figure BDA0000689569750000026
其中：Qco为包壳外表面的热流密度，该项由其他程序计算给出(如反应堆堆芯分析程序)；hfilm为film层导热系数。
(2)化学沉积层温差ΔTcrud
Figure BDA0000689569750000027
其中：kcrud为化学沉积层厚度，会随着燃耗的变化而不同；tcrud为化学沉积层导热系数，约等于0.5。
(3)包壳外表面温度T(Rco)
T(Rco)＝Tmod+ΔTcrud+ΔTfilm (6)
若T(Rco)大于饱和冷却剂饱和温度，则认为泡核沸腾将包壳外表面的热量通过化学沉积层带走，此时可以忽略ΔTcrud，包壳外表面的温度可以用Jens-Lottes沸腾传热方程计算。
(4)包壳层温差ΔTclad
Figure BDA0000689569750000031
其中：kclad为包壳导热系数；Rci为包壳内径；Rco为包壳外径。
(5)气隙层温差ΔTgap
Figure BDA0000689569750000032
其中：hgap为气隙导热系数；Rhd为热肿胀后芯块半径。
(6)芯块温度计算
将芯块分为J个子环区域，第j个子环区域区的温差ΔTj通过公式(9)计算得到：
Figure BDA0000689569750000033
其中：Qj为第j个子环区域的线功率密度；qj为第j个子环区域的体积功率密度；kj为第j个子环区域的导热系数；rj为第j个子环区域的半径。
则第j个子环区域的平均温度可由外至内累加得到，如公式(10)所示：
Figure BDA0000689569750000034
上述过程虽然是燃料棒温度计算的标准流程，但在实际过程可操作性低，例如燃料的导热系数很难精确给出，并且计算过程复杂，计算效率低。
发明内容
(一)要解决的技术问题
鉴于上述技术问题，本发明提供了一种获得核燃料棒内部燃料芯块温度分布的方法，以简化燃料棒内部燃料芯块温度分布的计算流程，满足实际工程需求。
(二)技术方案
根据本发明的一个方面，提供了一种获得核燃料棒内部燃料芯块温度分布的方法。该方法包括：步骤A：获得燃料棒有效温度TFeff和燃料芯块理论温度分布
Figure BDA0000689569750000041
的函数表达式；步骤B：获得燃料棒燃耗深度BU；步骤C：计算燃料棒包壳平均温度
Figure BDA0000689569750000042
步骤D：由燃料棒燃耗深度BU和燃料棒功率p，依据燃料棒有效温度TFeff的函数表达式，计算燃料棒有效温度TFeff；步骤E：由燃料棒包壳平均温度
Figure BDA0000689569750000043
和燃料棒有效温度TFeff，计算燃料芯块的平均温度
Figure BDA0000689569750000044
步骤F：由燃料棒燃耗深度BU和燃料棒功率p，依据燃料芯块理论温度分布
Figure BDA0000689569750000045
的函数表达式，计算燃料芯块理论平均温度
Figure BDA0000689569750000046
以及步骤G：利用燃料芯块平均温度
Figure BDA0000689569750000047
和燃料芯块理论平均温度
Figure BDA0000689569750000048
对燃料芯块理论温度分布进行校正，得到校正后的燃料芯块温度分布Tpellet(r)。
(三)有益效果
从上述技术方案可以看出，本发明获得核燃料棒内部燃料芯块温度分布的方法具有以下有益效果：
(1)通过利用函数拟合的方式计算燃料棒内芯块温度分布，简化了核燃料棒内部燃料芯块温度分布的计算流程，计算速度更快，使得工程人员对燃料芯块温度的实时监控成为可能；
(2)传统的燃料棒内部芯块平均温度认为与燃料棒功率成正比，这种处理方式基于经验，精确度低，本发明可计算出芯块内部的温度分布信息，不仅计算速度更快，并且计算精度更高。
附图说明
图1为燃料棒横截面上温度场分区的示意图；
图2为根据本发明实施例获得核燃料棒内部燃料芯块温度分布的方法的流程图；
图3为图2所示方法中燃料芯块分区示意图计算燃料棒燃料芯块平均温度的步骤中对燃料芯块进行分区的示意图。
具体实施方式
本发明提供一种快速准确的获得核燃料棒内部燃料芯块温度分布的方法(Intra-Pin Temperature Distribution Calculation Method，简称IPT技术)，可方便地应用于多个相关领域，为安全分析和燃料性能分析提供数据支撑。
为使本发明的目的、技术方案和优点更加清楚明白，以下结合具体实施例，并参照附图，对本发明进一步详细说明。
在本发明的一个示例性实施例中，提供了一种获得核燃料棒内部燃料芯块温度分布的方法。利用本实施例的方法，可以计算每个时间点的核燃料棒内部的温度分布。
图2为根据本发明实施例获得核燃料棒内部燃料芯块温度分布的方法的流程图。如图2所示，本实施例获得核燃料棒内部燃料芯块温度分布的方法包括：
步骤A：采用拟合方法获得燃料棒有效温度TFeff和燃料芯块理论温度分布
Figure BDA0000689569750000051
的函数表达式；
首先定义燃料棒有效温度TFeff，如式(11)所示：
Figure BDA0000689569750000052
其中，
Figure BDA0000689569750000053
为燃料棒包壳平均温度，
Figure BDA0000689569750000054
为燃料棒燃料芯块平均温度，θf为加权系数，为经验系数，一般取0＜θf≤02。
本发明认为，TFeff是燃料棒功率p和燃耗深度BU的函数，以多项式的形式表示，本发明以四阶多项式为例，如下所示：
TFeff(p，BU)＝Tmod+c1(BU)p+c2(BU)p2+c3(BU)p3+c4(BU)p4 (12)
其中，Tmod为慢化剂温度，为已知量。
由燃料性能计算程序计算出不同的燃料棒功率pn(n＝1，2，3，4)、不同的燃料棒燃耗深度BUm(m＝1，2，…，M)的TFeff(p，BU)，组成方程组，如公式(13)所示：
Figure BDA0000689569750000061
对于每一燃料棒燃耗深度BUm，计算四个不同的燃料棒功率pn下的燃料棒有效温度TFeff(p，BU)，联立式(13)，可得系数c1(BUm)、c2(BUm)、c3(BUm)、c4(BUm)。
本发明中，认为燃料棒内部燃料芯块理论温度分布
Figure BDA0000689569750000066
与燃料棒功率p、燃料棒燃耗深度BU、及燃料芯块中心点的距离r有关，燃料芯块温度分布可写成三者的函数形式，如公式(14)所示：
Figure BDA0000689569750000062
函数类型可以用多项式、指数函数、高斯多项式、幂函数、三角函数等，视实际工程需要而定。本实施例中以多项式为例进行说明，即将燃料芯块理论温度分布展开为二阶多项式的形式，如公式(15)所示：
Figure BDA0000689569750000063
用燃料性能分析程序计算多组标准工况下的燃料芯块温度分布，组成方程组，如公式(16)所示。其中，标准工况是指不同的燃料棒燃耗深度BUm(m＝1，2，…，M)，不同的燃料棒功率pi(i＝1，2，…，I)。
Figure BDA0000689569750000064
根据公式(16)可以求得燃料芯块理论温度分布
Figure BDA0000689569750000067
的函数表达式的展开系数a0(BUm)、a1(BUm)、a2(BUm)、a3(BUm)、a4(BUm)、a5(BUm)。
可见，经过多项式拟合，可以获得燃料棒有效温度TFeff的函数表达式和燃料芯块理论温度分布
Figure BDA0000689569750000065
的函数表达式，这些函数表达式将用于后续步骤。
步骤B：计算当前时刻，即第i个时间段终止点，的燃料棒燃耗深度BUi，EOS；
当最初计算时，i＝1，而后，随着i的逐步递增，就可以计算后续各个时刻的燃料棒燃耗深度BUi，EOS。
由于燃料的燃耗过程是一个时间概念，因此在本发明中，人为将时间划分成若干个时间段，每一时间段的起始点用BOS表示，每一时间段的终止点用EOS表示。
在该步骤中，利用第i个时间段长度ti和燃料棒功率pi计算出该时间段的燃耗深度变化量ΔBUi，与初始燃耗深度BUi，BOS加和，得到第i时间段终止点的燃料棒燃耗深度BUi，EOS，如式17所示：
BUi，EOS＝BUi，BOS+C·pi，BOS·ti·θ+C·pi，EOS·ti·(1-θ) (17)
其中，BUi，EOS为第i个时间段终止点的燃耗深度(同时也是第i+1个时间段起始点的燃耗深度，即BUi+，BOS＝BUi，EOS)，pi，BOS为第i个时间段起始点的燃料棒功率，pi，EOS第i个时间段终止点的燃料棒功率，θ为预估校正因子，一般取0.5，C为燃耗深度转换因子，为常量。对于第i个时间段，BUi，BOS、pi，BOS、pi，EOS、ti都是已知的，因此根据式(17)可以很方便的求得BUi，EOS。
后续步骤(步骤C-步骤F)均是对每一时间步结束点的处理，为描述方便，在相应的参数中会省略时间段的标记i。其中，上述的第i个时间段终止点，的燃料棒燃耗深度BUi，EOS也简写为BU。
步骤C：计算第i个时间段终止点的燃料棒包壳平均温度
Figure BDA0000689569750000071
本实施例中，采用如下公式计算燃料棒包壳平均温度：
Figure BDA0000689569750000072
其中：Tmod为第i个时间段终止点的包壳外侧的慢化剂的温度，ΔTclad为燃料棒包壳温升，计算公式为
Figure BDA0000689569750000073
其中：Qco为包壳外侧的慢化剂的热流密度；Rci和Rco分别为包壳内径和外径；kclad为包壳导热系数。
对于当前时刻，Qco，i是已知的，而kclad、Rci、Rco为常量，因此可求得第i时刻的包壳温升ΔTclad，i。而第i个时间段终止点的包壳外侧的慢化剂温度Tmod也是已知的，因此由式(18)可求得当前时刻的燃料包壳平均温度
Figure BDA0000689569750000081
步骤D：依据步骤A得到的燃料棒有效温度TFeff的函数表达式，由第i个时间段终止点的燃料棒燃耗深度BU和燃料棒功率p，计算第i个时间段终止点的燃料棒有效温度TFeff；
第i个时间段终止点的燃料棒有效温度TFeff为：
TFeff＝Tmod+c1p+c2p2+c3p3+c4p4 (19)
其中，对于第i个时间段终止点而言，燃料棒燃耗深度BU和燃料棒功率p是已知的。
其中，cj(BU)是由步骤A得到的cj(BUm)获得的多个燃料棒燃耗深度BUM下函数表达式系数的二维表经由线性插值得到，即：
Figure BDA0000689569750000082
其中，BU的取值介于BUz和BUz+1之间，k＝1、2、3、4。
步骤E：由步骤C得到的第i个时间段终止点的燃料棒包壳平均温度
Figure BDA0000689569750000083
和步骤D得到的第i个时间段终止点的燃料棒有效温度TFeff计算第i个时间段终止点的燃料芯块的平均温度
Figure BDA0000689569750000084
燃料棒由燃料包壳与燃料芯块组成，工程应用中，燃料棒有效温度可由燃料棒包壳平均温度与燃料棒燃料芯块平均温度加权组成，如公式(21)所示：
Figure BDA0000689569750000085
其中，θf为加权系数，定义如前所述。
由公式(21)可得第i个时间段终止点的燃料芯块的平均温度
Figure BDA0000689569750000086
计算式：
Figure BDA0000689569750000091
其中，
Figure BDA0000689569750000092
由步骤C获得，TFeff由步骤D获得。
步骤F：由燃料棒燃耗深度BU和燃料棒功率p，依据步骤A获得的燃料芯块理论温度分布Tpellet(r)的函数表达式，计算第i个时间段终止点的燃料芯块理论平均温度
Figure BDA0000689569750000093
该步骤F具体又可以包括：
子步骤F1：将燃料棒燃料芯块根据距离中心点的近分为N个燃料环，如图3所示；
将燃料棒燃料芯块根据距离中心点的远近分为N个燃料环，对于第n个燃料环，其内、外径分别用rn和rn+1表示，其中：n＝1、2、……、N。
子步骤F2：对于第n个燃料环，依据步骤A获得的燃料芯块理论温度分布Tpellet(r)的函数表达式计算该燃料环的理论平均温度
Figure BDA0000689569750000094
其中：n＝1、2、……、N：
燃料芯块理论温度分布的函数表达式由步骤A中式(15)表示，即：
Figure BDA0000689569750000095
其中，aj(BU)是由步骤A得到的aj(BUm)获得的二维表经由线性插值得到，即：
Figure BDA0000689569750000096
其中，BU的取值介于BUz和BUz+1之间，j＝0、1、2、3、4、5。
根据式(23)，分别求出每一燃料环上内径的理论温度
Figure BDA0000689569750000097
和外径的理论温度
Figure BDA0000689569750000098
则第n个燃料环的理论平均温度
Figure BDA0000689569750000099
用公式(25)计算：
Figure BDA00006895697500000910
子步骤F3：利用燃料芯块中每一燃料环的平均温度
Figure BDA0000689569750000101
合并出燃料棒燃料芯块的平均温度
Figure BDA0000689569750000102
如公式26所示：
Figure BDA0000689569750000103
其中，rpellet为整个燃料芯块的外径。
步骤G：利用步骤E求得的第i个时间段终止点的燃料芯块平均温度
Figure BDA0000689569750000104
和步骤F求得的燃料芯块理论平均温度
Figure BDA0000689569750000105
对燃料芯块理论温度分布Tpellet(r)进行校正，得到校正后的燃料芯块温度分布，如式(27)所示。
Figure BDA0000689569750000106
与燃料芯块理论温度分布Tpellet(r)相比，校正后的燃料芯块温度分布Tpellet(r)更为准确，更能反映燃料芯块温度分布的实际情况。
步骤H：令i＝i+1，执行步骤B，计算第i+1时间段的燃料芯块温度分布，直至用户设定的时间终止点。
通过本实施例的方法，工程师可以快速方便的获得反应堆内每一时刻每一燃料棒的芯块温度分布，获知最高芯块温度的数值，可帮助工程师快速做出反应堆是否安全的判断，有利于核反应堆的安全。
需要说明的是，本实施例是一个完整的获取每一时刻燃料棒芯块温度分布的方法，而在实际工程中，如果步骤B中的燃料棒燃耗深度BUi，EOS已知的情况下，本发明可以直接计算某一个时间点的燃料棒芯块温度分布，而不关心其前一时间段燃料棒芯块温度分布的状态。本领域技术人员很清楚其实现方式，此处不再赘述。
至此，已经结合附图对本实施例进行了详细描述。依据以上描述，本领域技术人员应当对本发明获得核燃料棒内部燃料芯块温度分布的方法有了清楚的认识。
此外，在附图或说明书描述中，相似或相同的部分都使用相同的图号，附图中未绘示或描述的实现方式，为所属技术领域中普通技术人员所知的形式。并且，上述对各元件和方法的定义并不仅限于实施例中提到的各种具体结构、形状或方式，本领域普通技术人员可对其进行简单地更改或替换，例如：
(1)步骤A中多项式表达形式还可以用幂函数、三角函数、双曲函数等形式表达；
(2)低阶展开式可以用高阶展开式来代替；
(3)实施例可提供包含特定值的参数的示范，但参数无需确切等于相应的值，而是可在可接受的误差容限或设计约束内近似于相应的值。
综上所述，本发明可快速精确计算出与反应堆安全关系最大的燃料芯块温度分布，更好地保障了反应堆安全，提升了核电站的安全性，可广泛应用于反应堆设计分析、安全分析、模拟仿真、反应堆监测等诸多领域。
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限制本发明，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (10)

1.一种获得核燃料棒内部燃料芯块温度分布的方法，其特征在于，包括：
步骤A：获得燃料棒有效温度TFeff和燃料芯块理论温度分布
Figure FDA0002226135890000011
的函数表达式；
步骤B：获得燃料棒燃耗深度BU；
步骤C：计算燃料棒包壳平均温度
Figure FDA0002226135890000012
步骤D：由燃料棒燃耗深度BU和燃料棒功率p，依据燃料棒有效温度TFeff的函数表达式，计算燃料棒有效温度TFeff；
步骤E：由燃料棒包壳平均温度
Figure FDA0002226135890000013
和燃料棒有效温度TFeff，计算燃料芯块的平均温度
Figure FDA0002226135890000014
步骤F：由燃料棒燃耗深度BU和燃料棒功率p，依据燃料芯块理论温度分布
Figure FDA0002226135890000015
的函数表达式，计算燃料芯块理论平均温度
Figure FDA0002226135890000016
以及
步骤G：利用燃料芯块平均温度
Figure FDA0002226135890000017
和燃料芯块理论平均温度
Figure FDA0002226135890000018
对燃料芯块理论温度分布进行校正，得到校正后的燃料芯块温度分布Tpellet(r)。
2.根据权利要求1所述的方法，其特征在于：
所述步骤A中，设定多个燃料棒燃耗深度BUM，对于每一燃料棒燃耗深度BUm，通过拟合的方式获得在该燃耗深度BUm下燃料棒有效温度TFeff和燃料芯块理论温度分布
Figure FDA0002226135890000019
的函数表达式；
所述步骤D中，燃料棒有效温度TFeff函数表达式中的各项系数，由步骤A中获得的多个燃料棒燃耗深度下的燃料棒有效温度TFeff函数表达式系数经由当前燃料棒燃耗深度BU差值得到。
3.根据权利要求2所述的方法，其特征在于，所述步骤A中，采用多项式拟合、幂函数拟合、三角函数拟合或双曲函数拟合的方式来获得燃料棒有效温度TFeff和燃料芯块理论温度分布
Figure FDA00022261358900000110
4.根据权利要求1所述的方法，其特征在于，所述步骤B中获得燃料棒燃耗深度BU的方法为：
直接获得；或
利用如下公式计算第i个时间段的终止点燃料棒燃耗深度BU：
BU＝BUi，BOS+C·pi，BOS·ti·θ+C·pi，EOS·ti·(1-θ)
其中：
BUi，BOS为第i个时间段起始点的燃料棒燃耗深度；
pi，BOS为第i个时间段起始点的燃料棒功率；
pi，EOS第i个时间段终止点的燃料棒功率；
C为燃耗深度转换因子，θ为预估校正因子，ti为第i个时间段长度。
5.根据权利要求1所述的方法，其特征在于，所述步骤C中，利用如下公式计算燃料棒包壳平均温度
Figure FDA0002226135890000021
Figure FDA0002226135890000022
其中：
Tmod为包壳外侧的慢化剂的温度；
ΔTclad为燃料棒包壳温升，计算公式为
Figure FDA0002226135890000023
其中：Qco为包壳外侧的慢化剂的热流密；Rco和Rci为包壳外径和内径；kclad为包壳导热系数。
6.根据权利要求1所述的方法，其特征在于，所述步骤E中，采用如下公式计算燃料芯块的平均温度
Figure FDA0002226135890000024
Figure FDA0002226135890000025
其中，θf为加权系数。
7.根据权利要求1所述的方法，其特征在于，所述步骤F包括：
子步骤F1：将燃料棒燃料芯块根据距离中心点远近分为N个燃料环；
子步骤F2：对于第n个燃料环，依据燃料芯块理论温度分布
Figure FDA0002226135890000026
的函数表达式计算其理论平均温度
Figure FDA0002226135890000031
其中：n＝1、2、……、N；以及
子步骤F3：利用燃料芯块中每一燃料环的平均温度
Figure FDA0002226135890000032
合并出燃料棒燃料芯块的平均温度
Figure FDA0002226135890000033
8.根据权利要求7所述的方法，其特征在于，所述子步骤F2中，对于第n个燃料环：
依据燃料芯块理论温度分布
Figure FDA0002226135890000034
的函数表达式分别计算其内径和外径的理论温度，而后经平均获得该第n个燃料环的理论平均温度
Figure FDA0002226135890000035
其中，所述燃料芯块理论温度分布
Figure FDA0002226135890000036
的函数表达式中的各项系数，由步骤A中获得的多个燃料棒燃耗深度下的燃料芯块理论温度分布
Figure FDA0002226135890000037
的函数表达式系数经由当前燃料棒燃耗深度BU差值得到。
9.根据权利要求1所述的方法，其特征在于，所述步骤G中，按照以下公式得到校正后的燃料芯块温度分布Tpellet(r)：
Figure FDA0002226135890000038
10.根据权利要求1至9中任一项所述的方法，其特征在于，设定时间段序号i，其初始值为1；
所述步骤B～步骤G均为针对时间段i的终止点而言，所述步骤G之后还包括：
步骤H：令i＝i+1，重新执行步骤B，计算新的第i个时间段终止点的校正后的燃料芯块温度分布Tpellet(r)。
CN201510142261.1A 2015-03-27 2015-03-27 获得核燃料棒内部燃料芯块温度分布的方法 Active CN106156460B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510142261.1A CN106156460B (zh) 2015-03-27 2015-03-27 获得核燃料棒内部燃料芯块温度分布的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510142261.1A CN106156460B (zh) 2015-03-27 2015-03-27 获得核燃料棒内部燃料芯块温度分布的方法

Publications (2)

Publication Number Publication Date
CN106156460A CN106156460A (zh) 2016-11-23
CN106156460B true CN106156460B (zh) 2020-03-17

Family

ID=57340623

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510142261.1A Active CN106156460B (zh) 2015-03-27 2015-03-27 获得核燃料棒内部燃料芯块温度分布的方法

Country Status (1)

Country Link
CN (1) CN106156460B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108256212B (zh) * 2018-01-16 2021-08-31 中南大学 多棒燃料元件并行性能分析方法、系统及存储介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4054487A (en) * 1974-06-04 1977-10-18 Tokyo Shibaura Denki Kabushiki Kaisha Nuclear fuel rods
SU1198323A1 (ru) * 1984-01-20 1985-12-15 Purmal Modris Ya Способ работы факельной топки
CN101809674A (zh) * 2007-07-30 2010-08-18 原子燃料工业株式会社 多普勒反应性系数测量方法
CN103020468A (zh) * 2012-12-26 2013-04-03 中山大学 核反应堆核热耦合计算方法
CN103902784A (zh) * 2014-04-11 2014-07-02 华北电力大学 用于超临界水堆瞬态核热耦合的安全分析计算装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4054487A (en) * 1974-06-04 1977-10-18 Tokyo Shibaura Denki Kabushiki Kaisha Nuclear fuel rods
SU1198323A1 (ru) * 1984-01-20 1985-12-15 Purmal Modris Ya Способ работы факельной топки
CN101809674A (zh) * 2007-07-30 2010-08-18 原子燃料工业株式会社 多普勒反应性系数测量方法
CN103020468A (zh) * 2012-12-26 2013-04-03 中山大学 核反应堆核热耦合计算方法
CN103902784A (zh) * 2014-04-11 2014-07-02 华北电力大学 用于超临界水堆瞬态核热耦合的安全分析计算装置

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
"Numerical investigation of heat transfer in upward flows of supercritical water in circular tubes and tight fuel rod bundles";Jue Yang，Yoshiaki Oka.et.al.;《Nuclear Engineering and Design》;20071231;第420-430页 *
"Temperature Distribution on Fuel Rods: A Study on the Effect of Eccentricity in the Position of UO2 Pellets";Gaspar Jr., J. C. A.;《Nuclear Energy For New Europe》;20111231;第1-11页 *
"微观燃耗模型在ＣＯＳＩＮＥ软件包堆芯程序中的研究与应用";胡啸宇，王苏等;《原子能科学技术》;20131231;第47卷;第503-506页 *

Also Published As

Publication number Publication date
CN106156460A (zh) 2016-11-23

Similar Documents

Publication Publication Date Title
Conner et al. CFD methodology and validation for single-phase flow in PWR fuel assemblies
Ryu et al. Solution of the BEAVRS benchmark using the nTRACER direct whole core calculation code
CN107066745B (zh) 获取快中子堆堆芯瞬态过程三维中子通量密度分布的方法
Heidmann et al. An experimental study of the effect of wake passing on turbine blade film cooling
CN112699620A (zh) 基于计算流体力学的反应堆堆芯热工水力特性分析方法
CN106156460B (zh) 获得核燃料棒内部燃料芯块温度分布的方法
Hawkes et al. Thermal predictions of the AGR-3/4 experiment with post irradiation examination measured time-varying gas gaps
CN112231960B (zh) 二维移动热传导模型和模型建立方法、应用方法
Chen et al. Development of Subchannel Code SUBSC for high-fidelity multi-physics coupling application
Meehan et al. Sensitivity analysis of transient critical heat flux by RIA under high-pressure flow boiling conditions in TRTL
Maekawa et al. Prediction of weld residual stress in a pressurized water reactor pressurizer surge nozzle
Li et al. A Combined Method for Predicting the Boron Deposited Mass and the CIPS Risk
Inaba et al. Evaluation of maximum fuel temperature in HTTR
Kendrick et al. CASL multiphysics modeling of crud deposition in PWRs
Ostanek et al. Establishing a methodology for resolving convective heat transfer from complex geometries
Rose et al. Redwing: a MOOSE Application for Coupling MPACT and BISON
Ren et al. Implementation of the Lax-Wendroff Method in COBRA-TF for Solving Two-Phase Flow Transport Equations
Vasiliev et al. Application of thermal hydraulic and severe accident code SOCRAT/V3 to bottom water reflood experiment QUENCH-LOCA-1
Canepa et al. Preparation of nuclear libraries with deterministic and stochastic methods for LWR reflectors
CN112182849B (zh) 再淹没临界后换热分析方法及装置
CN117473873B (zh) 基于DeepM&Mnet神经网络的核热耦合实现方法
Musa et al. An Experimental Reversed Heat Flux Investigation of the Helium-Cooled Modular Divertor with Multiple Jets
Zhang et al. Point kinetics calculations with fully coupled thermal fluids reactivity feedback
Kyriakopoulos et al. Validation of MOOSE’s subchannel module using the AREVA FCTF heated bundle benchmark
Kiss et al. Further details of a numerical analysis on the thermal hydraulic effect of wrapped wire spacers in fuel bundle

Legal Events

Date Code Title Description
C06 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
CP03 Change of name, title or address

Address after: 102209 Room 401, 4th floor, No. C, Ya'an commercial building, north of beiqijiasi Road, Changping District, Beijing

Patentee after: STATE POWER INVESTMENT CORPORATION Research Institute

Country or region after: China

Address before: 102209 Room 401, 4th floor, No. C, Ya'an commercial building, north of beiqijiasi Road, Changping District, Beijing

Patentee before: STATE NUCLEAR POWER Research Institute

Country or region before: China

CP03 Change of name, title or address