核工程专业公式大全（终极版）

核工程专业公式大全目录TOC\o"1-1"\h\z\u前言：说明与阅读指南 1第一类：核物理基础公式 2第二类：核衰变与放射性公式 3第三类：中子截面与反应率公式 5第四类：中子慢化与能谱公式 7第五类：中子扩散方程 9第六类：反应堆临界理论 11第七类：反应堆动力学与反应性控制 13第八类：热工水力学基础公式 15第九类：核燃料与燃料循环公式 18第十类：辐射防护与屏蔽公式 20第十一类：核安全与事故分析公式 22第十二类：概率安全评价（PSA）公式 23第十三类：核聚变物理公式 25附录：补充说明 27核工程专业公式大全前言：说明与阅读指南本汇编涵盖核工程领域的全部核心公式体系，按照核物理基础、中子物理、反应堆理论、热工水力、燃料循环、辐射防护、安全分析、聚变物理等学科模块展开。每个公式包含：标准形式、符号含义、应用场景、限制条件、计算举例等五大要素。符号统一说明如下：符号含义单位σ微观截面b（1b=10⁻²⁴cm²），cm²Σ宏观截面cm⁻¹φ中子通量密度n·cm⁻²·s⁻¹v中子速度cm/sn中子数密度n·cm⁻³N核密度cm⁻³λ衰变常数/中子代时间s⁻¹/sΛ中子代时间sβ缓发中子份额无量纲ρ反应性pcm（1pcm=10⁻⁵Δk/k）k_eff有效增殖因子无量纲Q热功率/反应热W，MeVT温度K，℃k热导率W·m⁻¹·K⁻¹h换热系数W·m⁻²·K⁻¹D扩散系数cmL扩散长度cmτ慢化年龄/约束时间cm²/sE能量eV，MeVm质量kg，uρ_m密度g·cm⁻³A活度Bq，Cid厚度/距离cm，mμ线衰减系数cm⁻¹ω角频率rad/sα多普勒系数pcm/Kγ空泡系数pcm/%void

第一类：核物理基础公式1.1质量亏损与结合能爱因斯坦质能方程：Δ符号：ΔE为结合能（MeV），Δm为质量亏损（u），c为光速（2.9979×10⁸m/s）。应用场景：计算原子核稳定性、评估核反应释放能量。限制条件：适用于核反应前后质量差异计算。举例：U-235的Δm≈0.186u，ΔE≈173.7MeV。比结合能：E符号：E_B/A为平均每核子的结合能（MeV/u），A为质量数。应用场景：原子核稳定性排序、裂变与聚变能量来源解释。1.2核反应Q值计算Q符号：m_i为入射粒子质量，m_t为靶核质量，m_p为出射粒子质量，m_f为剩余核质量。应用场景：核反应能量释放评估、靶核计算、聚变与裂变反应设计。限制条件：需精确质能当量换算（1u=931.494MeV/c²）。举例：D+T→⁴He+n反应中，m(D)=2.014102u，m(T)=3.016049u，m(⁴He)=4.002603u，m(n)=1.008665u，Q=(5.030151-5.011268)×931.5=17.6MeV。1.3核半径公式R符号：R为核半径（cm），r₀为核半径常数（1.2×10⁻¹³cm或1.25fm），A为质量数。应用场景：估算核尺寸、反应截面计算。限制条件：近似公式，实际核具有弥散边界。举例：U-235核半径R=1.2×10⁻¹³×235^{1/3}≈7.4×10⁻¹³cm。

第二类：核衰变与放射性公式2.1放射性衰变定律衰变速率：A放射性衰变定律：N活度：A符号：N为放射性核数，N₀为初始核数，λ为衰变常数（s⁻¹），A为活度（Bq）。应用场景：放射性核素衰变预测、放射性废物处置计算、医学同位素剂量估算、辐照后冷却时间计算。限制条件：适用于无限小时间间隔和独立衰变事件，统计规律要求大样本。举例：Cs-137（λ=7.29×10⁻¹⁰s⁻¹，半衰期30.17年），初始活度10⁷Bq，10年后A=10⁷×e⁻⁷・²⁹×¹⁰⁻¹⁰׳・¹⁵⁶×¹⁰⁸≈7.3×10⁶Bq。2.2半衰期与平均寿命半衰期：T平均寿命：τ符号：T_{1/2}为半衰期（s），τ为平均寿命（s）。应用场景：放射性废物分类（短半衰期/长半衰期）、燃料设计。举例：I-131半衰期8.02天，平均寿命约11.6天。2.3放射性平衡与衰变链母-子体衰变：N长期平衡（λ₁≪λ₂）：N2.4衰变热计算（Wigner-Way公式）P符号：P为衰变功率（W），t为冷却时间（s），t₀为辐照时间。应用场景：乏燃料池设计、反应堆停堆后冷却安全分析。限制条件：经验公式，用于近似估算，不同堆型存在差异。2.5核素活度与质量换算A符号：m为核素质量（g），M为摩尔质量（g/mol），N_A为阿伏伽德罗常数（6.022×10²³mol⁻¹）。应用场景：放射性源项制备、剂量估算、燃料裂变产物分析。

第三类：中子截面与反应率公式3.1核反应率与截面定义反应率密度：R宏观截面：Σ总反应率：R符号：R为反应率（cm⁻³・s⁻¹），N为单位体积原子核数（cm⁻³），σ为微观截面（b，1b=10⁻²⁴cm²），Φ为中子通量密度（n・cm⁻²・s⁻¹），Σ为宏观截面（cm⁻¹）。应用场景：反应堆功率计算、屏蔽设计、燃料燃耗估算。限制条件：截面与中子能量强相关，需能谱平均。举例：Φ=3×10¹³n·cm⁻²·s⁻¹，Σ_f=0.1cm⁻¹，R_f=3×10¹³×0.1=3×10¹²fissions·cm⁻³·s⁻¹。3.2布勒特-魏格纳公式（共振截面）σ符号：σ_0为峰值截面，E₀为共振能量，Γ为总宽度，Γ_γ为辐射俘获宽度。应用场景：共振吸收截面计算、反应堆能谱分析、多普勒效应分析。限制条件：适用于单能级窄共振，多能级需叠加。举例：U-238在6.67eV处有强共振，σ_0约20000b。3.3多普勒展宽公式Γσ符号：k_B为玻尔兹曼常数（8.617×10⁻⁵eV/K），M为靶核质量，T为温度（K）。应用场景：燃料多普勒反馈分析、事故分析。限制条件：适用于低能中子共振区。3.4有效增殖因子（六因子公式）kk符号：ε为快中子裂变因子，p为逃脱共振俘获概率，f为热中子利用系数，η为每次裂变平均热中子产额，P\{FNL}为快中子不泄露概率，P\{TNL}为热中子不泄露概率。各因子详细定义：εpfη应用场景：反应堆临界计算、装料方案设计、控制棒效应分析。限制条件：适用于热中子反应堆，快堆需修正。举例：PWR典型值：ε≈1.03，p≈0.85，f≈0.7，η≈2.0，P\{FNL}≈0.97，P\{TNL}≈0.97，k_{eff}≈1.03×0.85×0.7×2.0×0.97×0.97≈1.15。3.5四因子公式k

第四类：中子慢化与能谱公式4.1弹性散射能量变化最大能量损失（质心系）：αE最小散射后能量：E符号：A为靶核质量数，θ_c为质心系散射角，E、E'为散射前后中子能量。应用场景：慢化剂选择、能谱设计、反应性控制。举例：H-1（A=1）：α=0；D-2（A=2）：α=1/9；C-12（A=12）：α≈0.716。4.2对数能降μ平均对数能降：ξ轻核近似：ξ4.3慢化能力与慢化比慢化能力：ξ慢化比：ξ应用场景：不同慢化剂性能比较。举例：H₂O：ξΣ_s≈1.35cm⁻¹，慢化比≈70；D₂O：ξΣ_s≈0.18cm⁻¹，慢化比≈5000。4.4中子慢化能谱（费米谱）Φ热中子能谱（麦克斯韦分布）：N裂变中子能谱（麦克斯韦形式）：N瓦特谱：N应用场景：堆芯能谱分析、多群常数制备、探测器响应计算。限制条件：热堆中1/E谱适用于慢化区，有强吸收时偏离。

第五类：中子扩散方程5.1中子输运方程（玻尔兹曼方程）1符号：Φ为中子角通量密度，Ω→为方向向量，Σ_t为总宏观截面，S应用场景：中子输运解析、全堆芯物理计算、屏蔽设计、事故分析。限制条件：七维方程（时间×空间3×角度2×能量），求解极其困难。5.2中子扩散方程瞬态扩散方程：1稳态扩散方程：D有源稳态扩散方程：∇符号：D为扩散系数（cm），L为扩散长度（cm），φ为中子通量密度。应用场景：临界计算、通量分布计算。限制条件：扩散近似成立条件——各向同性散射、弱吸收、远离边界与强吸收区。5.3扩散系数DΣ符号：Σ_{tr}为输运截面，μ0应用场景：扩散系数计算、群常数制备。5.4扩散长度与慢化长度扩散长度平方：L慢化长度：τ徙动长度：M符号：L为扩散长度（cm），τ为慢化年龄（cm²），M为徙动长度（cm）。应用场景：泄露效应分析、堆芯尺寸设计。举例：H₂O慢化剂，L²≈2.5cm²，τ≈28cm²，M≈5.5cm。5.5费米年龄方程∇符号：q为慢化密度（慢化到τ以下的中子数/cm³・s），τ为慢化年龄（cm²）。

应用场景：慢化过程解析、中子年龄计算。5.6边界条件（齐次化条件）外推距离：d真空边界：ϕ通量连续条件：ϕ5.7中子守恒与中子流中子流密度：J→

第六类：反应堆临界理论6.1临界条件k临界反应堆：k\{eff}=1；次临界：k\{eff}<1；超临界：k_{eff}>1。反应性：ρ6.2平板反应堆临界方程B材料曲率：B临界条件：B符号：B_g²为几何曲率（cm⁻²），B_m²为材料曲率。应用场景：临界尺寸计算、燃料富集度确定。举例：圆柱堆B_g²=(π/H+d)²+(2.405/R+d)²。6.3圆柱反应堆临界方程d贝塞尔函数解：ϕ几何曲率：B临界条件：(6.4球面反应堆临界方程ϕ几何曲率：B6.5反应堆功率与通量关系功率密度：q堆芯总功率：P符号：E_f为每次裂变释放能量（~200MeV≈3.2×10⁻¹¹J），q'''为体积释热率（W・cm⁻³）。举例：P=3000MWth，需约9.4×10¹⁹次裂变/秒。6.6中子注量与功率换算Φ应用场景：热工水力-中子学耦合设计。

第七类：反应堆动力学与反应性控制7.1点堆动力学方程缓发中子点堆方程（带i组缓发中子）：dnd符号：n(t)为中子密度，C_i(t)为第i组缓发中子先驱核浓度，β_i为第i组缓发中子份额，β=Σβ_i为总缓发中子份额，Λ为中子代时间（s），λ_i为第i组衰变常数。应用场景：核反应堆瞬态分析、控制棒提棒/落棒分析、事故分析、功率调节系统设计。限制条件：适用于堆芯空间分布均匀的近似（点堆近似），空间效应显著时不适用。举例：阶跃反应性ρ=0.5β（β≈0.0065），功率短期响应显著。7.2中子倍增公式（阶跃反应性输入近似解）：n7.3反应性方程（倒时方程）ρ缓发中子近似：ρ符号：l为中子寿命（s），ω为角频率（rad/s）。应用场景：反应性测定、动态参数提取、功率变化频率分析。7.4Nordheim-Fuchs模型事故分析中功率与反应性的关系：dρ应用场景：反应性事故分析（如控制棒弹出事故、临界事故）。7.5中子代时间Λ符号：l_{mod}为中子代平均寿命（含慢化时间），Λ为中子代时间（主要受缓发中子控制）。应用场景：点堆动力学参数计算、反应堆周期估算。7.6反应堆周期T符号：T_p为反应堆周期（s）。

应用场景：提升功率速率控制、安全分析。7.7反应性系数多普勒系数：α慢化剂温度系数：α空泡系数：α功率系数：α反应性反馈总和：ρ应用场景：反应堆固有安全性评估、功率运行稳定性分析。举例：PWR中α_D为负（-2～-3pcm/K），空泡系数为负（-50～-200pcm/%void），自稳性良好。

第八类：热工水力学基础公式8.1质量守恒方程一维稳态：d一维瞬态：∂ρ两相流质量守恒：∂符号：ρ为密度（kg/m³），u为流速（m/s），A为截面积（m²），α为空泡份额。应用场景：堆芯通道热工水力设计、两相流分析。限制条件：一维模型需通过子通道方法近似三维效应。8.2动量守恒方程（一维）稳态：dP两相动量方程：-符号：P为压力（Pa），f为摩擦系数，D_h为水力直径（m）。应用场景：冷却剂循环分析、泵扬程设计、LOCA事故分析。8.3能量守恒方程（一维）稳态单相：W瞬态：ρ符号：c_p为比热容（J・kg⁻¹・K⁻¹），W为质量流率（kg/s），P_h为加热周长（m），q''为热通量密度（W・m⁻²）。8.4热传导方程圆柱坐标（燃料棒径向）：1稳态：1应用场景：燃料中心温度计算、热应力分析。限制条件：k通常随温度变化，需数值求解。8.5燃料棒温度分布（稳态）定热导率时中心温度：T变热导率积分形式：T符号：T_s为包壳表面温度（K），T_0为中心温度，R为燃料半径（m）。举例：UO₂燃料（T_s=600K，R=0.005m，q'''=5×10⁸W/m³，k≈4.5W/m・K）：ΔT≈695K，T_0≈1295K。8.6间隙热传导包壳表面热通量：q有效传热系数：1应用场景：燃料温度精确计算、芯块-包壳力学相互作用分析。8.7换热系数公式努塞尔数：Nu雷诺数：Re普朗特数：Pr应用场景：对流换热计算、热交换器设计、失水事故后传热分析。8.8Dittus-Boelter公式（湍流强制对流）Nun=0.4（加热）、0.3（冷却）。应用场景：单相流动沸腾前传热计算、冷却剂通道设计。限制条件：适用于光滑圆管，棒束修正：Nu=Nu\{tube}×F\{geometry}。8.9压降计算摩擦压降：Δ局部压降：Δ加速压降：Δ重力压降：Δ总压降：Δ符号：f为范宁摩擦因子，G为质量流速（kg・m⁻²・s⁻¹），K为形阻系数。应用场景：冷却剂循环设计、泵选型、自然循环分析。8.10两相流基本参数滑速比：S空泡份额：α质量含气率：x漂移流模型：u符号：j=αu_g+(1-α)u_f为体积通量密度，C₀为分布参数（\1.1\1.3），u_gj为漂移速度。应用场景：沸腾两相流动分析、反应堆热工水力瞬态。8.11临界热流密度（CHF）经验公式CHF与DNBR定义：DNBR应用场景：燃耗限值、安全I类准则、正常运行工况安全评估。设计准则：DNBR≥1.3（正常运行），DNBR≥1.16（事故工况）。8.12膜态沸腾传热Nu应用场景：再淹没阶段热工分析、ECCS性能验证。

第九类：核燃料与燃料循环公式9.1燃耗深度计算BU燃耗百分数：F符号：BU为燃耗深度（MWd/tU），P为热功率（MW），t为满功率运行时间（day），M_U为初始铀装量（tU）。应用场景：燃料更换周期规划、乏燃料后处理设计。举例：P=3000MW，M_U=80t，运行365天，BU=3000×365/80≈13688MWd/tU。9.2同位素燃耗方程d符号：N\i为第i种核素核密度，F\{yield}为裂变产额项。应用场景：燃料演变计算、反应堆物理-燃耗耦合计算。9.3钚再生/转换比转换比：CR增殖比（BR）：BR9.4裂变产物毒物效应平衡氙-135毒物：Xe平衡氙浓度：X氙振荡条件：ρ氙毒物反应性：ρ符号：γ_X为氙产额（约0.003），λ_X为Xe-135衰变常数（2.09×10⁻⁵s⁻¹）。应用场景：功率变化后碘坑效应分析、反应堆启动后氙毒物计算。举例：高功率停堆后约5-10小时达到最大氙毒物（负反应性约-0.1~-0.3Δk/k）。9.5钐-149毒物ρ符号：γ_Sm为钐产额（约0.013），λ_Sm为Pm-149衰变常数。应用场景：长寿命毒物效应分析（约40小时达到饱和）。

第十类：辐射防护与屏蔽公式10.1辐射强度衰减窄束γ射线衰减：I宽束γ射线衰减（含累积因子）：I符号：I为射线强度，μ为线衰减系数（cm⁻¹），x为屏蔽厚度（cm），B(x)为累积因子。应用场景：屏蔽设计、剂量估算、辐射防护最优化。限制条件：窄束公式仅对理想平行束有效，实际宽束必须使用B(x)。举例：²⁴¹Am（59.5keVγ），混凝土μ≈0.4cm⁻¹，20cm后衰减I/I₀≈0.00034。10.2累积因子经验公式（泰勒公式）B伯杰公式：B等比级数公式（G-P公式）：B符号：A、α₁、α₂、a、b、K均为拟合参数，与辐射能量和屏蔽材料相关。应用场景：γ射线宽束屏蔽计算、医院及核设施屏蔽设计、累积因子快速估算。限制条件：参数由经验拟合获得，不同能区和材料参数不同。10.3剂量与剂量率吸收剂量率（γ）：D剂量当量率：H有效剂量：E符号：μ_en/ρ为质能吸收系数，w_R为辐射权重因子，w_T为组织权重因子。应用场景：剂量计算、事故后辐射后果评估。10.4中子剂量转换中子注量至有效剂量：E符号：C_n(E_n)为中子注量-剂量转换系数（Sv・cm²），ICRP74给出推荐值。

应用场景：中子屏蔽设计、中子辐射防护。10.5点核积分法点源剂量率：D符号：r→p为计算点坐标，r→s为源点坐标，S_i为源强，应用场景：γ射线屏蔽快速计算（3D辐射场表征）、复杂几何体剂量场估算。限制条件：忽略多次散射和光子能量变化的精细处理，适合工程设计初期估算。10.6中子屏蔽（快中子）快中子衰减（距离平方反比修正）：Φ符号：Σ_{rem}为移除截面（cm⁻¹），表征快中子被慢化/吸收的总能力。应用场景：快中子屏蔽层设计（如压力容器外石蜡或含氢材料层）。10.7布拉格-格雷空腔原理D符号：J为通过空腔的粒子注量，(S/ρ)

第十一类：核安全与事故分析公式11.1锆-水反应速率ZrBaker-Just公式：Δ符号：Q为反应热（约5.8MJ/kgZr），E_a为活化能（约125kJ/mol），C为指前因子。应用场景：严重事故中氢产量、包壳氧化评估（如福岛事故分析）。限制条件：仅适用于温度约1000-1500℃范围。11.2安全壳直接加热（DCH）双隔间平衡（TCE）模型dd应用场景：严重事故下安全壳超压载荷计算、安全壳失效概率分析。限制条件：需要大量经验模型参数。11.3包壳峰值温度（PCT）准则10CFR50.46准则：PCT≤1477K（1204℃），锆-水反应率≤17%。11.4源项方程（放射性核素释放）d符号：R_i为迁移/释放率（s⁻¹）。应用场景：事故后辐射后果评价、应急计划区划分、裂变产物释放率计算。

第十二类：概率安全评价（PSA）公式12.1故障树分析顶事件概率（与门）：P顶事件概率（或门）：P含共因失效修正：P应用场景：PSA1级分析（堆芯损伤概率）、系统可靠性分析、核电站安全审评。限制条件：假设基本事件相互独立（实际可能含有依赖关系，需共因失效模型修正）。12.2事件树分析事故序列概率：P总风险：Risk符号：P\{IE}为始发事件频率，P\{branch,i}为第i个分支成功/失效概率，C_{seq}为后果度量。应用场景：PSA2级、3级分析（安全壳响应、场外剂量后果）、核电站风险评价。12.3设备可靠性参数失效概率：P平均无故障时间（MTBF）：MTBF不可用度：U=λλ+μ应用场景：设备可靠性数据库建立、维修周期优化、全厂断电（SBO）概率估算。12.4马尔可夫模型d符号：P(t)应用场景：多状态系统可靠性分析（如氢爆事件树与传统马尔可夫模型对比）。限制条件：状态数增加时计算复杂度迅速增大。12.5始发事件频率确定方法多种方法：工程判断（EngineeringJudgment）、参考已有PSA结果、EPRI始发事件清单、系统化逻辑分类、电厂能量平衡故障树、运行经验分析、失效模式与效应分析（FMEA）等。12.6易损性函数与贝叶斯网络（多灾害PSA）P符号：a_m为地震动峰值加速度，a_c为中值抗震能力，β为对数标准差。应用场景：外部灾害（地震/洪水/极端天气）风险评价。

第十三类：核聚变物理公式13.1劳森判据能量收支平衡条件：P净功率输出条件：nT聚变三乘积：nT符号：n为等离子体数密度（m⁻³），T为离子温度（keV），τ_E为能量约束时间（s）。应用场景：聚变装置设计（托卡马克/仿星器/惯性约束），点火条件判定、聚变增益因子Q计算。限制条件：不同聚变燃料（D-T/D-D/p-B¹¹）判据不同；磁约束与惯性约束热力学条件不同。举例：ITER设计n=1×10²⁰m⁻³，T=15keV，τ_E