CN114444413A - 一种板状燃料堆芯亚通道级三维热工水力分析方法 - Google Patents
一种板状燃料堆芯亚通道级三维热工水力分析方法 Download PDFInfo
- Publication number
- CN114444413A CN114444413A CN202210073056.4A CN202210073056A CN114444413A CN 114444413 A CN114444413 A CN 114444413A CN 202210073056 A CN202210073056 A CN 202210073056A CN 114444413 A CN114444413 A CN 114444413A
- Authority
- CN
- China
- Prior art keywords
- fuel
- channel
- grid
- coolant
- equation
- 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.)
- Granted
Links
- 239000000446 fuel Substances 0.000 title claims abstract description 172
- 238000004458 analytical method Methods 0.000 title claims abstract description 26
- 239000002826 coolant Substances 0.000 claims abstract description 88
- 239000007787 solid Substances 0.000 claims abstract description 73
- 238000000034 method Methods 0.000 claims abstract description 40
- 238000012546 transfer Methods 0.000 claims abstract description 27
- 238000004364 calculation method Methods 0.000 claims abstract description 23
- 230000008878 coupling Effects 0.000 claims abstract description 12
- 238000010168 coupling process Methods 0.000 claims abstract description 12
- 238000005859 coupling reaction Methods 0.000 claims abstract description 12
- 239000012530 fluid Substances 0.000 claims description 77
- 238000005253 cladding Methods 0.000 claims description 53
- 238000012937 correction Methods 0.000 claims description 13
- 239000008188 pellet Substances 0.000 claims description 13
- 238000004134 energy conservation Methods 0.000 claims description 10
- 230000014509 gene expression Effects 0.000 claims description 10
- 230000008569 process Effects 0.000 claims description 10
- 238000004422 calculation algorithm Methods 0.000 claims description 8
- 238000001727 in vivo Methods 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000005192 partition Methods 0.000 claims description 6
- 230000000704 physical effect Effects 0.000 claims description 5
- 229910052770 Uranium Inorganic materials 0.000 claims description 4
- JFALSRSLKYAFGM-UHFFFAOYSA-N uranium(0) Chemical compound [U] JFALSRSLKYAFGM-UHFFFAOYSA-N 0.000 claims description 4
- 238000012935 Averaging Methods 0.000 claims description 3
- 230000009471 action Effects 0.000 claims description 3
- 238000013459 approach Methods 0.000 claims description 3
- 238000009835 boiling Methods 0.000 claims description 3
- 230000000295 complement effect Effects 0.000 claims description 3
- 238000006073 displacement reaction Methods 0.000 claims description 3
- 230000008030 elimination Effects 0.000 claims description 3
- 238000003379 elimination reaction Methods 0.000 claims description 3
- 230000004907 flux Effects 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 238000012986 modification Methods 0.000 claims description 3
- 230000004048 modification Effects 0.000 claims description 3
- 239000003758 nuclear fuel Substances 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 238000012216 screening Methods 0.000 claims description 3
- 238000004781 supercooling Methods 0.000 claims description 3
- 230000001131 transforming effect Effects 0.000 claims description 3
- 238000004088 simulation Methods 0.000 abstract description 8
- 238000011160 research Methods 0.000 abstract description 4
- 238000005516 engineering process Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 7
- 230000001052 transient effect Effects 0.000 description 3
- 230000000712 assembly Effects 0.000 description 2
- 238000000429 assembly Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- ZOXJGFHDIHLPTG-UHFFFAOYSA-N Boron Chemical compound [B] ZOXJGFHDIHLPTG-UHFFFAOYSA-N 0.000 description 1
- 241000270295 Serpentes Species 0.000 description 1
- 239000006096 absorbing agent Substances 0.000 description 1
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 1
- 229910052782 aluminium Inorganic materials 0.000 description 1
- 229910052796 boron Inorganic materials 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000003381 stabilizer Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- 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
-
- 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/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- 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
- Y02E30/00—Energy generation of nuclear origin
- Y02E30/30—Nuclear fission reactors
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Fluid Mechanics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computing Systems (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Monitoring And Testing Of Nuclear Reactors (AREA)
Abstract
本发明公开了一种板状燃料堆芯亚通道级三维热工水力分析方法，该方法主要内容如下：1、建立由通道级网格与亚网格组成的亚通道级网格；2、对生成的两类网格进行标记；3、在控制方程中附加源项、修正几何与引入湍流搅混以构建冷却剂流动换热分析模型；4、针对原固体域设置离散节点，构造燃料板热工特性数值求解模型；5、建立流固耦合传热模型。本方法通过对板状燃料全堆芯进行亚通道级网格划分，在精确计算冷却剂温度场与流场、板状燃料温度分布的基础上，简化网格规模、显著降低计算资源，同时能够捕捉冷却剂流动与换热细节，为在工程中使用数值模拟技术对板状燃料全堆芯进行热工水力特性分析提供了一种兼具效率与精度的新型研究方法。
Description
技术领域
本发明属于板状燃料反应堆数值模拟技术领域，具体涉及到一种板状燃料堆芯亚通道级三维热工水力分析方法。
背景技术
板状燃料由于其功率密度高、燃耗较深、排布紧凑、换料周期长的特点，被广泛应用于小型堆、测试堆与一体化实验堆中。其核心结构单元为由燃料芯块、包壳与冷却剂矩形窄缝流道组成的闭式通道。在进行板状燃料反应堆正常或事故工况下的瞬态安全分析时，燃料板最高温度、包壳表面温度、冷却剂流场等热工水力参数极为重要。但由于燃料组件排布的复杂性与矩形窄缝通道数量的繁多，实现针对板状燃料反应堆全堆芯的数值模拟目前仍具有较大挑战性。
现阶段国内外针对板状燃料反应堆的热工水力分析方法主要包括以下三类：
1、一维系统分析程序——如RELAP5、MITARS、TRACE，此类程序能够构建核电厂宏观回路与各类热力构件模型，同时满足在较短时间内计算核电厂系统各处稳态参数分布与瞬态响应的要求。但该类程序使用离散的一维节点描述堆内各类系统级构件，如蒸汽发生器、稳压器、阀门、主泵，在程序中均采用单个节点表示，无法刻画构件内物理过程与参数分布。对于复杂的堆芯结构，仅使用单节点的集总值显然无法满足精细尺度下堆芯安全分析的需求。
2、子通道分析程序——如COBRA、COSINE，此类程序将堆芯划分为相邻的并联流道，并采用各类描述热工水力参数间关系的经验表达式计算子通道内冷却剂流动特性值，相比于一维系统程序能够较精确地反映堆芯各流道内部的流动换热细节，但此类程序对于瞬变工况模拟能力仍具有较大不足，且目前现有的子通道程序大多基于压水堆棒束通道进行开发，在板状燃料反应堆中应用较少。
3、计算流体力学(CFD)仿真软件——如FLUENT、CFX、STAR-CCM+，此类软件要求针对堆芯几何模型进行微尺度的网格划分，通过求解控制体质量、动量、能量守恒方程获得流场参数与温度分布。CFD仿真精度较高，结果可视化能力高，在采用恰当计算模型与方法的条件下能够捕捉冷却剂精细流动与传热现象。但同 CFD方法要求在几何建模软件中构建精细的堆芯结构，所需求的建模工作量较大；且计算速度缓慢，计算资源耗费量大。
发明内容
本发明的目的在于提供一种板状燃料堆芯亚通道级三维热工水力研究方法，该方法通过全堆芯亚通道级网格划分，固体域节点温度场差分计算与流固耦合模型的建立，能够在保证计算精度的条件下，显著降低计算资源与计算时间。
为了达到上述目的，本发明采用如下技术方案：
一种板状燃料堆芯亚通道级三维热工水力研究方法，包括如下步骤：
步骤1：建立由全堆芯通道级网格与流体域亚网格两部分组成的亚通道级网格
(1)：通过几何建模软件SOLIDWORKS建立板状燃料堆芯三维几何模型，此堆芯三维几何模型由众多相同的燃料芯块板、燃料包壳与冷却剂通道三部分组成；
(2)：在网格划分软件中基于步骤1-(1)中建立的板状燃料堆芯三维几何模型，在水平剖面上沿各燃料板中心对称轴进行纵向网格划分，生成既包含燃料板固体域同时也包含冷却剂流体域的通道级网格。需要特别说明的是，考虑到计算流体力学(CFD)软件的求解性质，对于网格内部介质的有限体积法数值求解必须针对流体；因此在后续处理中通道级网格内部原燃料板固体域被简化为与冷却剂性质相同的流体域，而由该简化方式引入的人工修正在步骤3中做了详细介绍；
(3)：为了精确捕捉矩形窄缝通道内冷却剂流动与换热细节，针对通道级网格中原流体域沿纵向均匀剖分若干节点，划分的网格中左右两列网格与毗邻固体域一同被合并为新的通道级网格，而流体域其余网格即为生成的亚网格，最后沿堆芯高度方向均匀划分预设数量节点；至此，完成全堆芯亚通道级网格的建立，通过该方式生成的网格数量远小于常规CFD模拟所需要的数目，为建模与网格划分过程节省了可观的工作量；
步骤2：对建立的全堆芯通道级网格与流体域亚网格进行网格标记
将生成的全堆芯亚通道级网格导入计算流体力学软件中，通过用户自定义函数分别对通道级网格与流体域亚网格进行网格标记，具体步骤如下：
设堆芯高度方向为z轴，垂直于燃料板长度方向为x轴，沿燃料板长度方向为y轴，假设每个标准燃料组件中燃料板的个数为NOP，则根据计算需要确定的通道级网格节点个数也为NOP，使用N表示某处通道级网格节点沿x方向所处的相对位置，1≤N≤NOP，N＝1时x＝Thc0+Thw+0.5ThU，其中x为节点与燃料组件边界横向距离，Thc0、Thw、ThU分别为组件边界包壳厚度、矩形窄缝厚度与铀燃料板厚度；当1＜N≤NOP时，DeltaN＝Δx(N,N-1)＝Thw+2Thc+ThU， Thc为燃料板包壳厚度，即相邻燃料板中心距为一定值；
对于流体域亚网格，采用补集筛选法对其进行网格标记，设x方向每两块燃料板之间对冷却剂通道均匀划分亚网格节点数为n，则第N个通道网格节点处网格坐标为：
另设一变量E，令满足上述坐标的节点处E＝1，则其余节点均为亚通道级网格节点，赋值E＝0；至此，对亚通道级网格中两套不同网格的标记完成，在后续分类过程中只需引入简单的判断语句遴选变量E的取值即可；
步骤3：通过针对标记的两类网格在控制守恒方程中附加源项、修正几何参数与引入通道间湍流搅混以构建冷却剂流动换热分析模型
由于通道级网格包络了原燃料板固体域，因此在对控制体进行有限体积法积分求解时作以下几何参数修正：
控制体体积修正：
控制体侧面积修正：
其中：
V'、V0分别为修正后与未修正的控制体体积；
Sy'、Sy0分别为修正后与未修正的控制体沿y方向表面积；
L、D分别为控制体沿垂直板排列方向的高度与宽度；
经过修正后，通道级网格内进行守恒方程的数值求解时所使用的几何参数均映射为真实的流体域参数值；
为代替全堆芯亚通道级网格划分方式对板状燃料堆芯几何固体域结构的简化带来的冷却剂流动特性的影响，需在守恒方程中引入额外的动量、能量源项，并满足以下假设：
1)冷却剂流动阻力使用由摩擦压降转换得到的动量源项
表示：
式中：
ρ——冷却剂密度/kg·m-3；
De——冷却剂通道水力直径/m，对于矩形通道通过
计算得出，Δx为控制体径向厚度；
f——摩擦阻力因子；
在计算因流动摩擦引入的动量损失时，特征长度取通道级网格的水力直径；针对层流、湍流分别采用不同的阻力摩擦因子计算关系式，引入无量纲数Re，Re
为雷诺数，通过下式计算得出：
其中μ为冷却剂的动力粘性系数/kg·m-1·s-1；
当Re<2300时，冷却剂流动视为层流，摩擦阻力因子由修正后Hagen-Poiseuill 公式计算得出；而当Re≥2300时，流动转为湍流，摩擦阻力因子通过Kakac经验关系式计算，具体表达式如下：
2)计算域中燃料板与包壳不发生结构形变与相对位移；
3)计算过程中冷却剂不发生过冷沸腾，始终保持单相流动；
基于上述假设，对流体控制方程中的动量方程和能量方程中添加源项：在动量方程中添加包含粘性损失、惯性损失的压降作为动量源项
在能量方程中添加燃料板表面的对流换热量Qhf作为能量源项SE；
分析通道级网格与流体域亚网格的划分特性后发现，其毗邻界面为纯流体域，因此在交界面处考虑相邻通道P、N间湍流搅混带来的动量与能量交换，具体形式为：
其中：
GP、GN分别为P、N通道内流体质量流密度/kg·m-2·s-1；
DP、DN分别为P、N通道宽度/m；
hP、hN为相邻P、N控制体内流体比焓/kJ·kg-1；
ΔxP,N为P、N两点间距离/m；
β为相邻通道间湍流搅混系数，为雷诺数与通道几何参数的函数，β通过以下关系式计算得出：
最终，考虑由通道级网格内部固体域的简化带来的动量损失、燃料板释热引入的能量源项与通道间湍流搅混作用后，给出亚通道网格与通道网格内流体域的守恒方程如下形式：
质量守恒方程：
动量守恒方程
至此，建立了完整的冷却剂流动换热特性分析模型；
步骤4：针对原固体域设置离散节点与构建燃料板热工特性数值求解模型
分析步骤1中通道级网格的划分特性发现，其内部同时包含由燃料板组成的原固体域与冷却剂流体域，虽然在步骤3中已经讨论了固体域的简化引入的模型修正，但简化后的原固体域仍不适用于流体守恒方程；因此在通道级网格内部引入独立的板状燃料元件热工特性分析模型对其温度场进行求解，具体如下：
首先提出以下假设：
1)反应堆燃料板结构均匀分布在堆芯中；
2)单个控制体内所有燃料板具有相同的几何特性；
3)单个控制体内所有燃料板具有相同的热工特性；
根据燃料板与包壳组成的固体域几何特性沿径向划分N个节点，其中最外层两个节点分布在包壳内、外侧边界区域，其余节点从燃料板对称中心向外侧依次排布，忽略燃料板的轴向导热，即对于单个控制体，其中添加的固体域温度求解算法是一维的；对于第i个计算节点，其能量守恒方程应满足以下形式：
式中：
ρi为计算节点i处控制体密度/kg·m-3；
Cpi为定压比热容kJ·kg-1·K-1；
Δxi为控制体径向宽度/m；
Ti为节点i处的温度/K；
qi-1,i与qi+1,i分别为i-1节点与i+1节点传递到i处的面热流密度/kW·m-2； Qi为节点i处控制体内体积释热率/kW·m-3。
其中，两侧相邻计算节点传递给中心节点的面热流密度通过傅立叶导热定律求得：
qi-1,i＝Ki-1,i(Ti-1-Ti)
qi+1,i＝Ki+1,i(Ti+1-Ti)
式中：
Ti-1、Ti+1分别为计算节点i-1、i+1处的温度/K；
Ki-1,i为计算节点i-1和i之间的传热系数/kW·m-2·K-1，等同于Ki,i-1；Ki+1,i为计算节点i+1和i之间的传热系数/kW·m-2·K-1，等同于Ki,i+1，均表示为控制体内燃料板几何参数与导热系数λ二者的函数：
式中xi为节点i距芯块中心的径向距离，某两节点间的平均导热系数由积分平均值得出；
根据由燃料板及包壳组成的固体域以及与冷却剂交界面的传热特点，建立如下边界条件：
1)燃料板中心为对称边界，因此对于计算节点1的控制方程为绝热形式: K0,1＝Ki-1,i＝0；
2)在燃料芯块与包壳之间存在接触导热，即在i＝N-1处qi,i-1＝hs(Ti-Ti-1)，其中hs为燃料芯块与包壳间的接触导热系数/kW·m-2·K-1，由下式计算得出：
式中：
ψ为带有量纲的常数，取10/m1/2；
λU、λc分别为燃料芯块、包壳的导热系数/W·m-1·K-1；
Pi为表面接触压力/Pa；
Θ为燃料芯块的Meyer硬度数/N·m-2；
δU、δc分别为在交界面处燃料芯块、包壳的表面粗糙度/m；
3)在燃料板包壳外表面，即i＝N-1处，Ki+1,i＝H，Ti+1＝Tf，则有 qi+1,i＝H(Tf-Ti)；
其中：Ki+1,i为计算节点i+1和i之间的传热系数kW·m-2·K-1；
H为燃料板包壳表面的对流换热系数kW·m-2·K-1；
qi+1,i为计算节点i+1传导到计算节点i的热量kW·m-2；
Tf为矩形窄缝内冷却剂主流温度/K；
Ti+1为计算节点i+1处的温度/K；
Ti为计算节点i处的温度/K；
为了实现温度场的程序求解，将能量守恒方程转化为差分方程形式：
式中：Ti n为计算节点i处在第n个时间步下的温度/K，Ti n-1为计算节点i处在第n-1个时间步下的温度/K，Ti n+1为计算节点i处在第n+1个时间步下的温度/K；
将差分方程变形为固体域温度场差分方程：
式中：系数Ai＝-Ki-1,i，系数
系数Ci＝-Ki+1,i，系数
对于每个通道级网格中固体域内N个节点都存在上述关于其自身节点温度与相邻两节点温度的三元差分方程，整理为一N+1元一次方程组，并在初始、最末方程中添加了附加的本构边界条件作为方程闭合的必要条件，对矩阵方程的化简采用高斯消去法，求解算法采用高斯迭代法，最终得到n时刻各个节点的温度值，至此，完成了针对原固体域离散节点的设置并构建了燃料板热工特性数值求解模型。
步骤5：在通道级网格内固体域与流体域交界面处通过换热关系建立流固耦合换热模型
为了进一步求解冷却剂热工水力参数，将通道网格中固体域边界，即燃料包壳最外层控制体与亚通道级网格中流体域节点联立并构建流固耦合换热模型：
对于流固边界，热流密度满足Dittus-Boleter公式：
式中：λ为流体导热系数/W·m-1·K-1；Pr为普朗特数；上式中包含受温度场影响而变化的流体物性参数；因此对于冷却剂流场的计算需要通过 Dittus-Boleter关系式与固体域温度场进行耦合求解，具体实现步骤如下：
1)：假定亚通道级网格内初始冷却剂流场与温度场并赋值，通过检索得到对应温度下的物性参数值，通过考虑附加的动量、能量源项和湍流搅混项的冷却剂守恒方程求得该时刻冷却剂流场；
2)：将通道网格内冷却剂流场参数代入Dittus-Boleter公式求解燃料包壳外表面对流换热系数H；
3)：以燃料包壳外表面对流换热系数H为边界条件求解步骤4中固体域温度场差分方程，获得燃料板及包壳温度分布；
4)：根据燃料包壳、矩形窄缝通道内冷却剂温度、换热面积与表面对流换热系数求得表面对流换热量；
5)：将得到的表面对流换热量作为通道网格内冷却剂控制方程中附加能量源项，根据冷却剂流动换热特性分析模型求解其下一时间步亚通道级网格内流场与温度场；
6)：获得该时刻下流体物性参数，返回步骤①进行通道级网格流场参数的迭代求解，最终得到所需时间内板状燃料反应堆全堆芯热工水力参数分布。
和现有技术相比较，本发明具有以下特点：
1、堆芯控制体划分方式简单，计算网格数量少；
2、在通道级网格的基础上，采用的进一步亚通道级网格划分能够精确捕捉矩形窄缝通道内流动换热细节；
3、克服了一维系统分析程序无法获得板状燃料堆芯固体域准确的温度分布、对功率分布不均匀瞬态计算过程响应模拟较差的问题；
4、针对燃料板温度场求解设计的算法模型普适性强，可在简单参数调节后应用于其他堆型的固体域传热计算；
5、耗费计算资源少，计算速度快且计算精度高；
6、构建的流固耦合换热模型可在获得冷却剂热工水力参数的同时计算得到燃料板、包壳的温度场，对板状燃料反应堆安全分析具有重要意义。
本发明方法已通过程序完成实践验证，在该方法下针对板状燃料全堆芯的数值仿真能够准确地获得燃料板、包壳温度场以及冷却剂流场、温度场等热工水力特性参数，具有精确捕捉冷却剂通道内流体域流动换热细节的能力，且同时具有计算资源耗费量少，并行能力高，计算速度快的特点。
附图说明
图1为相邻燃料板间通道级/亚通道级网格划分示意图。
图2为燃料板与包壳组成固体域内温度场差分算法节点选取示意图。
图3为流固耦合模型流程图。
图4为本发明板状燃料堆芯亚通道级三维热工水力分析方法系统框图。
具体实施方式
为了更深层次地介绍本发明技术路线与实施意义，以下根据图4所示系统框图为主要流程，结合附图对本发明作进一步的详细描述。
本发明提出了一种板状燃料堆芯亚通道级三维热工水力研究方法，包括以下步骤：
步骤1：建立由全堆芯通道级网格与流体域亚网格两部分组成的全堆芯亚通道级网格
(1)：通过几何建模软件建立板状燃料堆芯三维几何模型，此堆芯三维几何模型由众多相同的燃料板、燃料包壳与冷却剂通道三部分组成；常见的板状燃料反应堆燃料组件分别为板状燃料反应堆标准组件盒与控制组件盒结构示意图，主要区别在于控制组件盒两侧布置了硼材料制成的中子吸收体，而标准组件盒仅由铝制包壳包裹的铀燃料芯块板组成；
(2)：在网格划分软件中基于步骤1-(1)中建立的板状燃料堆芯三维几何模型，在水平剖面上沿各燃料板中心对称轴进行纵向网格划分，生成既包含燃料板固体域同时也包含冷却剂流体域的通道级网格；考虑到计算流体力学软件的求解性质，对于网格内部介质的有限体积法数值求解必须针对流体；因此在后续处理中通道级网格内部原燃料板固体域被简化为与冷却剂性质相同的流体域，而由该简化方式引入的人工修正在步骤3中做了详细介绍；
(3)：为了精确捕捉矩形窄缝通道内冷却剂流动与换热细节，针对通道级网格中原流体域沿纵向均匀剖分若干节点，划分的网格中左右两列网格与毗邻固体域一同被合并为新的通道级网格，而流体域其余网格即为生成亚网格，最后沿堆芯高度方向均匀划分预设数量节点；至此，完成全堆芯亚通道级网格的建立，建立的网格示意如图1所示，通过该方式生成的网格数量远小于常规CFD模拟所需要的数目，为建模与网格划分过程节省了可观的工作量；
步骤2：对建立的全堆芯通道级网格与流体域亚网格进行网格标记
将生成的全堆芯亚通道级网格导入计算流体力学软件中，通过用户自定义函数分别对通道级网格与流体域亚网格进行网格标记，具体步骤如下：
设堆芯高度方向为z轴，垂直于燃料板长度方向为x轴，沿燃料板长度方向为y轴，假设每个标准燃料组件中燃料板的个数为NOP，则根据计算需要确定的通道级网格节点个数也为NOP，使用N表示某处通道级网格节点沿x方向所处的相对位置，1≤N≤NOP，N＝1时x＝Thc0+Thw+0.5ThU，其中x为节点与燃料组件边界横向距离，Thc0、Thw、ThU分别为组件边界包壳厚度、矩形窄缝厚度与铀燃料板厚度；当1＜N≤NOP时，DeltaN＝Δx(N,N-1)＝Thw+2Thc+ThU， Thc为燃料板包壳厚度，即相邻燃料板中心距为一定值；
对于流体域亚网格，采用补集筛选法对其进行网格标记，设x方向每两块燃料板之间对冷却剂通道均匀划分亚网格节点数为n，则第N个通道网格节点处网格坐标为：
另设一变量E，令满足上述坐标的节点处E＝1，则其余节点均为亚通道级网格节点，赋值E＝0；至此，对亚通道级网格中两套不同网格的标记完成，在后续分类过程中只需引入简单的判断语句遴选变量E的取值即可；
步骤3：通过针对标记的两类网格在控制守恒方程中附加源项、修正几何参数与引入通道间湍流搅混以构建冷却剂流动换热分析模型
由于通道级网格包络了原燃料板固体域，因此在对控制体进行有限体积法积分求解时作以下几何参数修正：
控制体体积修正：
控制体侧面积修正：
其中：
V'、V0分别为修正后与未修正的控制体体积；
Sy'、Sy0分别为修正后与未修正的控制体沿y方向表面积；
L、D分别为控制体沿垂直板排列方向的高度与宽度；
经过修正后，通道级网格内进行守恒方程的数值求解时所使用的几何参数均映射为真实的流体域参数值；
为代替全堆芯亚通道级网格划分方式对板状燃料堆芯几何固体域结构的简化带来的冷却剂流动特性的影响，需在守恒方程中引入额外的动量、能量源项，并满足以下假设：
1)冷却剂流动阻力使用由摩擦压降转换得到的动量源项
表示：
式中：
ρ——冷却剂密度/kg·m-3；
De——冷却剂通道水力直径/m，对于矩形通道通过
计算得出，Δx为控制体径向厚度；
f——摩擦阻力因子；
在计算因流动摩擦引入的动量损失时，特征长度取通道级网格的水力直径；针对层流、湍流分别采用不同的阻力摩擦因子计算关系式，引入无量纲数Re，Re
为雷诺数，通过下式计算得出：
其中μ为冷却剂的动力粘性系数/kg·m-1·s-1；
当Re<2300时，冷却剂流动视为层流，摩擦阻力因子由修正后Hagen-Poiseuill 公式计算得出；而当Re≥2300时，流动转为湍流，摩擦阻力因子通过Kakac经验关系式计算，具体表达式如下：
2)计算域中燃料板与包壳不发生结构形变与相对位移；
3)计算过程中冷却剂不发生过冷沸腾，始终保持单相流动；
基于上述假设，对流体控制方程中的动量方程和能量方程中添加源项：在动量方程中添加包含粘性损失、惯性损失的压降作为动量源项
在能量方程中添加燃料板表面的对流换热量Qhf作为能量源项SE；
分析通道级网格与流体域亚网格的划分特性后发现，其毗邻界面为纯流体域，因此在交界面处考虑相邻通道P、N间湍流搅混带来的动量与能量交换，具体形式为：
其中：
GP、GN分别为P、N通道内流体质量流密度/kg·m-2·s-1；
DP、DN分别为P、N通道宽度/m；
hP、hN为相邻P、N控制体内流体比焓/kJ·kg-1；
ΔxP,N为P、N两点间距离/m；
β为相邻通道间湍流搅混系数，为雷诺数与通道几何参数的函数，β通过以下关系式计算得出：
最终，考虑由通道级网格内部固体域的简化带来的动量损失、燃料板释热引入的能量源项与通道间湍流搅混作用后，给出亚通道网格与通道网格内流体域的守恒方程如下形式：
质量守恒方程：
动量守恒方程
至此，建立了完整的冷却剂流动换热特性分析模型；
步骤4：针对原固体域设置离散节点与构建燃料板热工特性数值求解模型
分析步骤1中通道级网格的划分特性发现，其内部同时包含由燃料板组成的原固体域与冷却剂流体域，虽然在步骤3中已经讨论了固体域的简化引入的模型修正，但简化后的原固体域仍不适用于流体守恒方程；因此在通道级网格内部引入独立的板状燃料元件热工特性分析模型对其温度场进行求解，具体如下：
首先提出以下假设：
1)反应堆燃料板结构均匀分布在堆芯中；
2)单个控制体内所有燃料板具有相同的几何特性；
3)单个控制体内所有燃料板具有相同的热工特性；
根据燃料板与包壳组成的固体域几何特性沿径向划分N个节点，其中最外层两个节点分布在包壳内、外侧边界区域，其余节点从燃料板对称中心向外侧依次排布,最终的节点划分示意图如图2；忽略燃料板的轴向导热，即对于单个控制体，其中添加的固体域温度求解算法是一维的；对于第i个计算节点，其能量守恒方程应满足以下形式：
式中：
ρi为计算节点i处控制体密度/kg·m-3；
Cpi为定压比热容kJ·kg-1·K-1；
Δxi为控制体径向宽度/m；
Ti为节点i处的温度/K；
qi-1,i与qi+1,i分别为i-1节点与i+1节点传递到i处的面热流密度/kW·m-2； Qi为节点i处控制体内体积释热率/kW·m-3。
其中，两侧相邻计算节点传递给中心节点的面热流密度通过傅立叶导热定律求得：
qi-1,i＝Ki-1,i(Ti-1-Ti)
qi+1,i＝Ki+1,i(Ti+1-Ti)
式中：
Ti-1、Ti+1分别为计算节点i-1、i+1处的温度/K；
Ki-1,i为计算节点i-1和i之间的传热系数/kW·m-2·K-1，等同于Ki,i-1；Ki+1,i为计算节点i+1和i之间的传热系数/kW·m-2·K-1，等同于Ki,i+1，均表示为控制体内燃料板几何参数与导热系数λ二者的函数：
式中xi为节点i距芯块中心的径向距离，某两节点间的平均导热系数由积分平均值得出；
根据由燃料板及包壳组成的固体域以及与冷却剂交界面的传热特点，建立如下边界条件：
1)燃料板中心为对称边界，因此对于计算节点1的控制方程为绝热形式: K0,1＝Ki-1,i＝0；
2)在燃料芯块与包壳之间存在接触导热，即在i＝N-1处qi,i-1＝hs(Ti-Ti-1)，其中hs为燃料芯块与包壳间的接触导热系数/kW·m-2·K-1，由下式计算得出：
式中：
ψ为带有量纲的常数，取10/m1/2；
λU、λc分别为燃料芯块、包壳的导热系数/W·m-1·K-1；
Pi为表面接触压力/Pa；
Θ为燃料芯块的Meyer硬度数/N·m-2；
δU、δc分别为在交界面处燃料芯块、包壳的表面粗糙度/m；
3)在燃料板包壳外表面，即i＝N-1处，Ki+1,i＝H，Ti+1＝Tf，则有qi+1,i＝H(Tf-Ti)；
其中：Ki+1,i为计算节点i+1和i之间的传热系数kW·m-2·K-1；
H为燃料板包壳表面的对流换热系数kW·m-2·K-1；
qi+1,i为计算节点i+1传导到计算节点i的热量kW·m-2；
Tf为矩形窄缝内冷却剂主流温度/K；
Ti+1为计算节点i+1处的温度/K；
Ti为计算节点i处的温度/K；
为了实现温度场的程序求解，将能量守恒方程转化为差分方程形式：
式中：Ti n为计算节点i处在第n个时间步下的温度/K，Ti n-1为计算节点i处在第n-1个时间步下的温度/K，Ti n+1为计算节点i处在第n+1个时间步下的温度/K；
将差分方程变形为固体域温度场差分方程：
式中：系数Ai＝-Ki-1,i，系数
系数Ci＝-Ki+1,i，系数
对于每个通道级网格中固体域内N个节点都存在上述关于其自身节点温度与相邻两节点温度的三元差分方程，整理为一N+1元一次方程组，并在初始、最末方程中添加了附加的本构边界条件作为方程闭合的必要条件。对于初始节点与最末节点，分别不考虑前、后计算节点，即：
i＝1时，A1＝0
i＝n时，Cn＝0
将边界条件代入差分方程中整理可得：
左侧N×N系数为中心三对角矩阵，对矩阵方程的化简采用高斯消去法得到对角矩阵：
化简后对角矩阵求解算法采用高斯迭代法，最终得到n时刻各个节点的温度值；至此，完成了针对原固体域离散节点的设置并构建了燃料板热工特性数值求解模型；
步骤5：在通道级网格内固体域与流体域交界面处通过换热关系建立流固耦合换热模型
为了进一步求解冷却剂热工水力参数，将通道网格中固体域边界，即燃料包壳最外层控制体与亚通道级网格中流体域节点联立并构建流固耦合换热模型：
对于流固边界，热流密度满足Dittus-Boleter公式：
式中：λ为流体导热系数/W·m-1·K-1；Pr为普朗特数；上式中包含受温度场影响而变化的流体物性参数；因此对于冷却剂流场的计算需要通过 Dittus-Boleter关系式与固体域温度场进行耦合求解，求解流程如附图3所示，具体实现步骤如下：
1)：假定亚通道级网格内初始冷却剂流场与温度场并赋值，通过检索得到对应温度下的物性参数值，通过考虑附加的动量、能量源项和湍流搅混项的冷却剂守恒方程求得该时刻冷却剂流场；
2)：将通道网格内冷却剂流场参数代入Dittus-Boleter公式求解燃料包壳外表面对流换热系数H；
3)：以燃料包壳外表面对流换热系数H为边界条件求解步骤4中固体域温度场差分方程，获得燃料板及包壳温度分布；
4)：根据燃料包壳、矩形窄缝通道内冷却剂温度、换热面积与表面对流换热系数求得表面对流换热量；
5)：将得到的表面对流换热量作为通道网格内冷却剂控制方程中附加能量源项，根据冷却剂流动换热特性分析模型求解其下一时间步亚通道级网格内流场与温度场；
6)：获得该时刻下流体物性参数，返回步骤①进行通道级网格流场参数的迭代求解，最终得到所需时间内板状燃料反应堆全堆芯热工水力参数分布。
Claims (1)
1.一种板状燃料堆芯亚通道级三维热工水力分析方法，其特征在于：包括如下步骤：
步骤1：建立由全堆芯通道级网格与流体域亚网格两部分组成的全堆芯亚通道级网格
(1)：通过几何建模软件建立板状燃料堆芯三维几何模型，此堆芯三维几何模型由众多相同的燃料板、燃料包壳与冷却剂通道三部分组成；
(2)：在网格划分软件中基于步骤1-(1)中建立的板状燃料堆芯三维几何模型，在水平剖面上沿各燃料板中心对称轴进行纵向网格划分，生成既包含燃料板固体域同时也包含冷却剂流体域的通道级网格；考虑到计算流体力学软件的求解性质，对于网格内部介质的有限体积法数值求解必须针对流体；因此在后续处理中通道级网格内部原燃料板固体域被简化为与冷却剂性质相同的流体域，而由该简化方式引入的人工修正在步骤3中做了详细介绍；
(3)：为了精确捕捉矩形窄缝通道内冷却剂流动与换热细节，针对通道级网格中原流体域沿纵向均匀剖分若干节点，划分的网格中左右两列网格与毗邻固体域一同被合并为新的通道级网格，而流体域其余网格即为生成亚网格，最后沿堆芯高度方向均匀划分预设数量节点；至此，完成全堆芯亚通道级网格的建立；
步骤2：对建立的全堆芯通道级网格与流体域亚网格进行网格标记
将生成的全堆芯亚通道级网格导入计算流体力学软件中，通过用户自定义函数分别对通道级网格与流体域亚网格进行网格标记，具体步骤如下：
设堆芯高度方向为z轴，垂直于燃料板长度方向为x轴，沿燃料板长度方向为y轴，假设每个标准燃料组件中燃料板的个数为NOP，则根据计算需要确定的通道级网格节点个数也为NOP，使用N表示某处通道级网格节点沿x方向所处的相对位置，1≤N≤NOP，N＝1时x＝Thc0+Thw+0.5ThU，其中x为节点与燃料组件边界横向距离，Thc0、Thw、ThU分别为组件边界包壳厚度、矩形窄缝厚度与铀燃料板厚度；当1＜N≤NOP时，
Thc为燃料板包壳厚度，即相邻燃料板中心距为一定值；
对于流体域亚网格，采用补集筛选法对其进行网格标记，设x方向每两块燃料板之间对冷却剂通道均匀划分亚网格节点数为n，则第N个通道网格节点处网格坐标为：
另设一变量E，令满足上述坐标的节点处E＝1，则其余节点均为亚通道级网格节点，赋值E＝0；至此，对亚通道级网格中两套不同网格的标记完成，在后续分类过程中只需引入简单的判断语句遴选变量E的取值即可；
步骤3：通过针对标记的两类网格在控制守恒方程中附加源项、修正几何参数与引入通道间湍流搅混以构建冷却剂流动换热分析模型
由于通道级网格包络了原燃料板固体域，因此在对控制体进行有限体积法积分求解时作以下几何参数修正：
控制体体积修正：
控制体侧面积修正：
其中：
V'、V0分别为修正后与未修正的控制体体积；
Sy'、Sy0分别为修正后与未修正的控制体沿y方向表面积；
L、D分别为控制体沿垂直板排列方向的高度与宽度；
经过修正后，通道级网格内进行守恒方程的数值求解时所使用的几何参数均映射为真实的流体域参数值；
为代替全堆芯亚通道级网格划分方式对板状燃料堆芯几何固体域结构的简化带来的冷却剂流动特性的影响，需在守恒方程中引入额外的动量、能量源项，并满足以下假设：
1)冷却剂流动阻力使用由摩擦压降转换得到的动量源项
表示：
式中：
ρ——冷却剂密度/kg·m-3；
De——冷却剂通道水力直径/m，对于矩形通道通过
计算得出，Δx为控制体径向厚度；
f——摩擦阻力因子；
在计算因流动摩擦引入的动量损失时，特征长度取通道级网格的水力直径；针对层流、湍流分别采用不同的阻力摩擦因子计算关系式，引入无量纲数Re，Re为雷诺数，通过下式计算得出：
其中μ为冷却剂的动力粘性系数/kg·m-1·s-1；
当Re<2300时，冷却剂流动视为层流，摩擦阻力因子由修正后Hagen-Poiseuill公式计算得出；而当Re≥2300时，流动转为湍流，摩擦阻力因子通过Kakac经验关系式计算，具体表达式如下：
2)计算域中燃料板与包壳不发生结构形变与相对位移；
3)计算过程中冷却剂不发生过冷沸腾，始终保持单相流动；
基于上述假设，对流体控制方程中的动量方程和能量方程中添加源项：在动量方程中添加包含粘性损失、惯性损失的压降作为动量源项
在能量方程中添加燃料板表面的对流换热量Qhf作为能量源项SE；
分析通道级网格与流体域亚网格的划分特性后发现，其毗邻界面为纯流体域，因此在交界面处考虑相邻通道P、N间湍流搅混带来的动量与能量交换，具体形式为：
其中：
GP、GN分别为P、N通道内流体质量流密度/kg·m-2·s-1；
DP、DN分别为P、N通道宽度/m；
hP、hN为相邻P、N控制体内流体比焓/kJ·kg-1；
ΔxP,N为P、N两点间距离/m；
β为相邻通道间湍流搅混系数，为雷诺数与通道几何参数的函数，β通过以下关系式计算得出：
最终，考虑由通道级网格内部固体域的简化带来的动量损失、燃料板释热引入的能量源项与通道间湍流搅混作用后，给出亚通道网格与通道网格内流体域的守恒方程如下形式：
质量守恒方程：
动量守恒方程
至此，建立了完整的冷却剂流动换热特性分析模型；
步骤4：针对原固体域设置离散节点与构建燃料板热工特性数值求解模型
分析步骤1中通道级网格的划分特性发现，其内部同时包含由燃料板组成的原固体域与冷却剂流体域，虽然在步骤3中已经讨论了固体域的简化引入的模型修正，但简化后的原固体域仍不适用于流体守恒方程；因此在通道级网格内部引入独立的板状燃料元件热工特性分析模型对其温度场进行求解，具体如下：
首先提出以下假设：
1)反应堆燃料板结构均匀分布在堆芯中；
2)单个控制体内所有燃料板具有相同的几何特性；
3)单个控制体内所有燃料板具有相同的热工特性；
根据燃料板与包壳组成的固体域几何特性沿径向划分N个节点，其中最外层两个节点分布在包壳内、外侧边界区域，其余节点从燃料板对称中心向外侧依次排布，忽略燃料板的轴向导热，即对于单个控制体，其中添加的固体域温度求解算法是一维的；对于第i个计算节点，其能量守恒方程应满足以下形式：
式中：
ρi为计算节点i处控制体密度/kg·m-3；
Δxi为控制体径向宽度/m；
Ti为节点i处的温度/K；
qi-1,i与qi+1,i分别为i-1节点与i+1节点传递到i处的面热流密度/kW·m-2；
Qi为节点i处控制体内体积释热率/kW·m-3。
其中，两侧相邻计算节点传递给中心节点的面热流密度通过傅立叶导热定律求得：
qi-1,i＝Ki-1,i(Ti-1-Ti)
qi+1,i＝Ki+1,i(Ti+1-Ti)
式中：
Ti-1、Ti+1分别为计算节点i-1、i+1处的温度/K；
Ki-1,i为计算节点i-1和i之间的传热系数/kW·m-2·K-1，等同于Ki,i-1；Ki+1,i为计算节点i+1和i之间的传热系数/kW·m-2·K-1，等同于Ki,i+1，均表示为控制体内燃料板几何参数与导热系数λ二者的函数：
式中xi为节点i距芯块中心的径向距离，某两节点间的平均导热系数由积分平均值得出；
根据由燃料板及包壳组成的固体域以及与冷却剂交界面的传热特点，建立如下边界条件：
1)燃料板中心为对称边界，因此对于计算节点1的控制方程为绝热形式:
K0,1＝Ki-1,i＝0；
2)在燃料芯块与包壳之间存在接触导热，即在i＝N-1处qi,i-1＝hs(Ti-Ti-1)，其中hs为燃料芯块与包壳间的接触导热系数/kW·m-2·K-1，由下式计算得出：
式中：
ψ为带有量纲的常数，取10/m1/2；
λU、λc分别为燃料芯块、包壳的导热系数/W·m-1·K-1；
Pi为表面接触压力/Pa；
Θ为燃料芯块的Meyer硬度数/N·m-2；
δU、δc分别为在交界面处燃料芯块、包壳的表面粗糙度/m；
3)在燃料板包壳外表面，即i＝N-1处，Ki+1,i＝H，Ti+1＝Tf，则有qi+1,i＝H(Tf-Ti)；
其中：Ki+1,i为计算节点i+1和i之间的传热系数kW·m-2·K-1；
H为燃料板包壳表面的对流换热系数kW·m-2·K-1；
qi+1,i为计算节点i+1传导到计算节点i的热量kW·m-2；
Tf为矩形窄缝内冷却剂主流温度/K；
Ti+1为计算节点i+1处的温度/K；
Ti为计算节点i处的温度/K；
为了实现温度场的程序求解，将能量守恒方程转化为差分方程形式：
式中：Ti n为计算节点i处在第n个时间步下的温度/K，Ti n-1为计算节点i处在第n-1个时间步下的温度/K，Ti n+1为计算节点i处在第n+1个时间步下的温度/K；
将差分方程变形为固体域温度场差分方程：
式中：系数Ai＝-Ki-1,i，系数
系数Ci＝-Ki+1,i，系数
对于每个通道级网格中固体域内N个节点都存在上述关于其自身节点温度与相邻两节点温度的三元差分方程，整理为一N+1元一次方程组，并在初始、最末方程中添加了附加的本构边界条件作为方程闭合的必要条件，对矩阵方程的化简采用高斯消去法，求解算法采用高斯迭代法，最终得到n时刻各个节点的温度值，至此，完成了针对原固体域离散节点的设置并构建了燃料板热工特性数值求解模型；
步骤5：在通道级网格内固体域与流体域交界面处通过换热关系建立流固耦合换热模型
为了进一步求解冷却剂热工水力参数，将通道网格中固体域边界，即燃料包壳最外层控制体与亚通道级网格中流体域节点联立并构建流固耦合换热模型：
对于流固边界，热流密度满足Dittus-Boleter公式：
式中：λ为流体导热系数/W·m-1·K-1；Pr为普朗特数；上式中包含受温度场影响而变化的流体物性参数；因此对于冷却剂流场的计算需要通过Dittus-Boleter关系式与固体域温度场进行耦合求解，具体实现步骤如下：
1)：假定亚通道级网格内初始冷却剂流场与温度场并赋值，通过检索得到对应温度下的物性参数值，通过考虑附加的动量、能量源项和湍流搅混项的冷却剂守恒方程求得该时刻冷却剂流场；
2)：将通道网格内冷却剂流场参数代入Dittus-Boleter公式求解燃料包壳外表面对流换热系数H；
3)：以燃料包壳外表面对流换热系数H为边界条件求解步骤4中固体域温度场差分方程，获得燃料板及包壳温度分布；
4)：根据燃料包壳、矩形窄缝通道内冷却剂温度、换热面积与表面对流换热系数求得表面对流换热量；
5)：将得到的表面对流换热量作为通道网格内冷却剂控制方程中附加能量源项，根据冷却剂流动换热特性分析模型求解其下一时间步亚通道级网格内流场与温度场；
6)：获得该时刻下流体物性参数，返回步骤①进行通道级网格流场参数的迭代求解，最终得到所需时间内板状燃料反应堆全堆芯热工水力参数分布。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210073056.4A CN114444413B (zh) | 2022-01-21 | 2022-01-21 | 一种板状燃料堆芯亚通道级三维热工水力分析方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210073056.4A CN114444413B (zh) | 2022-01-21 | 2022-01-21 | 一种板状燃料堆芯亚通道级三维热工水力分析方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114444413A true CN114444413A (zh) | 2022-05-06 |
| CN114444413B CN114444413B (zh) | 2023-10-24 |
Family
ID=81367404
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210073056.4A Active CN114444413B (zh) | 2022-01-21 | 2022-01-21 | 一种板状燃料堆芯亚通道级三维热工水力分析方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114444413B (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115169265A (zh) * | 2022-07-28 | 2022-10-11 | 中国核动力研究设计院 | 基于数值分析的搅混系数分析方法、系统、设备和介质 |
| CN116415449A (zh) * | 2023-06-12 | 2023-07-11 | 天津渤化工程有限公司 | 顺酐反应器的设计方法、设计系统及信息数据处理终端 |
| CN116933553A (zh) * | 2023-08-02 | 2023-10-24 | 上海交通大学 | 数值反应堆中子学的非结构网格体积修正方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112699620A (zh) * | 2021-01-08 | 2021-04-23 | 西安交通大学 | 基于计算流体力学的反应堆堆芯热工水力特性分析方法 |
| CN112989651A (zh) * | 2021-02-06 | 2021-06-18 | 西安交通大学 | 反应堆堆芯多物理场耦合方法 |
| EP3896703A1 (en) * | 2020-04-14 | 2021-10-20 | Westinghouse Electric Sweden AB | Computer-based simulation methods for boiling water reactors (bwr) |
-
2022
- 2022-01-21 CN CN202210073056.4A patent/CN114444413B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP3896703A1 (en) * | 2020-04-14 | 2021-10-20 | Westinghouse Electric Sweden AB | Computer-based simulation methods for boiling water reactors (bwr) |
| CN112699620A (zh) * | 2021-01-08 | 2021-04-23 | 西安交通大学 | 基于计算流体力学的反应堆堆芯热工水力特性分析方法 |
| CN112989651A (zh) * | 2021-02-06 | 2021-06-18 | 西安交通大学 | 反应堆堆芯多物理场耦合方法 |
Non-Patent Citations (2)
| Title |
|---|
| 曾和义;郭;: "湍流模型在堆芯热工水力性能分析中的应用", 原子能科学技术, no. 01 * |
| 罗磊;陈文振;陈志云;朱倩;: "单个燃料元件热工水力三维数值模拟", 海军工程大学学报, no. 01 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115169265A (zh) * | 2022-07-28 | 2022-10-11 | 中国核动力研究设计院 | 基于数值分析的搅混系数分析方法、系统、设备和介质 |
| CN115169265B (zh) * | 2022-07-28 | 2023-09-05 | 中国核动力研究设计院 | 基于数值分析的搅混系数分析方法、系统、设备和介质 |
| CN116415449A (zh) * | 2023-06-12 | 2023-07-11 | 天津渤化工程有限公司 | 顺酐反应器的设计方法、设计系统及信息数据处理终端 |
| CN116415449B (zh) * | 2023-06-12 | 2023-08-22 | 天津渤化工程有限公司 | 顺酐反应器的设计方法、设计系统及信息数据处理终端 |
| CN116933553A (zh) * | 2023-08-02 | 2023-10-24 | 上海交通大学 | 数值反应堆中子学的非结构网格体积修正方法 |
| CN116933553B (zh) * | 2023-08-02 | 2024-02-13 | 上海交通大学 | 数值反应堆中子学的非结构网格体积修正方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114444413B (zh) | 2023-10-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114444413A (zh) | 一种板状燃料堆芯亚通道级三维热工水力分析方法 | |
| CN111261232B (zh) | 反应堆一回路冷却剂流场、温度场和应力场耦合计算方法 | |
| CN112699620B (zh) | 基于计算流体力学的反应堆堆芯热工水力特性分析方法 | |
| CN114282460B (zh) | 一种铅铋快堆堆芯热工水力特性分析方法 | |
| CN115859851B (zh) | 一种液态金属耦合超临界二氧化碳共轭传热的计算方法 | |
| CN111027205B (zh) | 一种液态金属冷却快堆堆芯盒间流数值模拟方法 | |
| CN115659908B (zh) | 一种印刷电路板换热器的三单元非平衡多孔介质方法 | |
| CN115862779B (zh) | 一种基于两方程热湍流模型的液态金属共轭传热计算方法 | |
| CN114154262B (zh) | 一种交叉流印刷电路板式换热器正向设计方法 | |
| Hung et al. | Optimal design of a concentric heat exchanger for high-temperature systems using CFD simulations | |
| Wang et al. | Numerical study on the purge gas flow and heat transfer characteristics in helium cooled solid breeder blanket of CFETR | |
| CN116562153B (zh) | 一种液态金属热分层特性的计算方法 | |
| Zhu et al. | Validation of the Porous Medium Approximation for Hydrodynamics Analysis in Compact Heat Exchangers | |
| Razzaghpanah et al. | Natural convection heat transfer from a bundle of heated circular cylinders with staggered arrangement immersed in molten solar salt | |
| Ye et al. | Numerical analysis of the flow behavior in a helically coiled once through steam generator | |
| Lee et al. | Dual-fluid topology optimization of printed-circuit heat exchanger with low-pumping-power design | |
| Zahid et al. | Experimental and CFD simulation study of shell and tube heat exchangers with different baffle segment configurations | |
| Feng et al. | Validation of the TWOPORFLOW code for the core analysis of liquid metal-cooled reactor with selected experiments | |
| Li et al. | Wenxi Tian, and Suizheng Qiu Xi'an Jiaotong University, Xi'an, Shaanxi, China dlzhang@ mail. xjtu. edu. cn | |
| Peng et al. | Study on anti-disturbance capacity of LBE cooled double loop natural circulation system under asymmetric conditions | |
| Wang et al. | CFD-Based Optimization of a Shell-and-Tube Heat Exchanger. | |
| Yildiz et al. | Large eddy simulation of 5-tube bundle helical coil steam generator test section | |
| Wu et al. | Numerical study of the flow and heat transfer characteristics within the annular fuel element with double-sided cooling | |
| Aldlemy et al. | Numerical simulation on the effect of pipe roughness in turbulent flow | |
| Yerramalle et al. | Pore scale numerical investigation of mixed convection from an isolated heat source in a channel with a porous insert |
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 |