原理绝密？你不会造氢弹本质上是数学不好！聊聊老专家们牛在哪里

军武数据库

昨天，在一篇讨论氢弹爆炸机制的文章下，看到有人依然在疑惑：氢弹的原理，难道不是最高级别的绝密吗？

说实话，这种问题我已经听得耳朵起茧了。原理级的东西真没什么好保密的，甚至普通人翻翻公开的物理教材、查查学术论文，就能拼凑出一大堆所谓的“绝密”信息。

真正的问题不在于原理是不是绝密，而是在于大多数人的意识中根本分不清“原理”和“工程”的差别。这是一种典型的认知偏差——或者更直白一点，叫认知割裂。

对于核武器的原理级别的东西其实特别简单，归纳一下就这么几条：

E=MC²，这是核心，意思是“质能等价”。所有的核武器无论是裂变武器还是聚变武器，亦或是今后出来的反物质武器都是遵循着这一基本核心原理工作的。在使用过程中，利用物理机制让物质的质量发生变化从而释放出巨大能量。这一点和炸药爆炸利用化学键之间的键能差所释放出的能量不属于一个数量级。因此，核武器威力巨大！

中子平均寿命Λ，，这个部分就有很多人会不理解，实际上kef是指中子的生成率/中子的损失率。也就是，它决定了链式反应一旦发生这个装置能不能自持。也是推算临界质量和临界半径的基本算法。

裂变和聚变的反应截面 σ，这是一个决定了核反应发生的概率。

这就属于统计物理范畴了，主要是衡量两粒子碰撞时发生特定过程的概率的指标。

我们都知道中子击中重核元素（铀235、钚239 等）会发生裂变反应，而氘原子与氘原子碰撞或者氘原子和氚原子碰撞可以产生聚变反应。但是“碰撞”这个事件是由一系列的概率事件合并在一起造成的。在物理学上这个概率叫做“截面”。截面是有单位的，这个单位叫做“靶恩（barn）”，同样这是一个面积单位相当于10⁻²⁸平方米。

经过计算，不同的粒子之间在不同反应条件下有不同的反应截面，这就为核武器材料的选择提供了基本依据。

例如我们常说的氢弹，是以氘-氚聚变反应支持的，但为什么不是氘-氘反应呢？这就是因为，氘-氚聚变反应的截面高，在即便是在较低的温度下，聚变反应就可以引发。

所以说，“σ”决定的是工程上要达到什么条件（温度、密度）才能实现有效反应。同样，“σ”也是各种核弹材料的选择依据。

劳森判据（Lawson Criterion），这是核聚变研究的一条基准，定义聚变反应自持的必要条件 其本质是等离子体内聚变反应产生能量的速率必须超过等离子体能量损失的速率。 这一标准对于实现聚变点火至关重要，因此被称为“判据”。

其实，这个判据就是温度、密度、时间所组成的条件积，更高的温度才能克服带正电的原子核之间的静电排斥力，使它们发生聚变；更高的密度意味着在给定的体积内可以容纳更多的燃料颗粒，从而增加聚变反应的频率；更长的约束时间才能允许发生更多的聚变反应，并可能引发聚变点火。这个判据无论是对我们的可控核聚变还是氢弹的设计都极其重要。简单的说就是可控核聚变目前所做的点火试验是依靠降低密度获得足够的温度和时间（主要是时间）下的聚变反应。而氢弹的爆炸则是极端的提高温度让聚变材料在一瞬间完成聚变和能量输出。

再有就是辐射驱动内爆平衡方程了。

E 是辐射能量密度，F 是辐射能流，S 是源项（裂变初级的能量输出），κ 是吸收系数。这个方程本质上是一个很普遍的能量守恒原理方程。聚变材料的加热算法就需要依靠这个方程来进行推导。可以保证用最小当量的原子弹来引爆最多的聚变材料。

然后就是氢弹聚变材料的状态方程

P 是压力，γ 是比热比，ρ 是密度，e 是内能密度。方程描述了在极端高温高压下，氢弹材料的热力学性质，如压力、密度、内能的关系。这是聚变材料稳定性、分布结构、外壳容器结构和材料的基本计算依据。

所以对于一枚氢弹来说，基本上核心内容就是这些方程了。这些方程其实和F=ma一样，都是最基本的学术性原理级的存在。如果觉得告诉大家是泄密了？大家就会造核弹了？那就太天真了。

F=ma在初中二年级就告诉大伙了，为了让大家记得且深入了解，老师不仅教了这个公式的原理是什么，还要不断在物理试卷里面考试确认大家记得且会用，但讽刺的是——99.9999%的人没有把它真正的用在自己的生活里吧？没有任何人依据这个牛顿运动公式造出一辆车来吧？

怎么氘了氢弹原子弹的时候，大伙就觉得知道了个公式就能造了呢？

难点在哪里呢？——工程计算！

就像绝对的美女只能存在于想象中一样，上面的所有公式也都只是在想象中——理想环境。在工程上，如果只依靠理论计算公式在工程中则是难以为继的。

原因很简单：工程计算必须面对真实的材料。

在理想情况下，大多数材料被当成绝对刚体、完美光滑、性质均匀——这种假设是为了让公式可以优雅地写在黑板上，也方便理论推导。只要按这种假设来算，结果永远和公式吻合得漂亮无比。

但到了真实世界，材料不是这样。它们会有缺陷、晶格畸变、表面粗糙、内部夹杂，温度和应力变化还会让材料性质瞬间改变。这意味着在实际工程计算中，公式必须被修正，甚至修正项有时候比原公式本身还重要。

这就是为什么，理论模型在实验室的白板上看上去像左边的“理想情况”，而一旦进入真实工程环境，就会变成右边的“实际情况”——形变、扰动、不可控的微观差异，全部压在计算模型的背上，让你不得不面对复杂得多的数值求解问题。

同样材料的配比也成了一个脱离公式的存在，公式计算中材料也是100%的纯度无丝毫瑕疵的。以钚弹来说，钚239的临界质量是10公斤，球体临界直径为99毫米。这是一个经过理论计算的数值，就是通过前面中子平均寿命Λ进行推算得到的一个结果。但是天底下就没有100%纯度的钚239！

钚是一种重要的人造元素，是通过核反应堆中一颗中子轰击一颗铀238原子后嬗变生成的。但这个事情任何人都难以保证被中子轰击的铀238原子只接收到了一颗中子。这样就不仅仅有钚239生成，同时还会生成钚240和钚241。

以钚240为例，它会自发衰变释放中子。这样一个钚核心内的中子通量就会因为钚240的存在而导致中子通量和计算不符。理论上的 10 公斤，到了工程里也许需要 11 公斤、12 公斤，或者是9.5公斤、8公斤，甚至在理论的安全装置之外还需要更复杂的临界安全冗余设计。

同样，核心外部的反射层大家知道是用金属铍来反射中子，但是要知道的是金属铍是一种相当坚固难以加工成型的的金属。

为了加工金属铍作为核心外壳，就需要向金属铍中增加铝元素使之具备可加工性。但要知道的是铝和金属铍的中子吸收、透射特性又不一样，因此，这个部分也需要进行修正。

类似的在工程上，各种理想公式到实际运用的工程算式中有大量的修正常数需要计算。这也是为什么明明都有了理想公式，还得再做大量计算的原因。举个例子：

在核弹的工程计算中，往往会有大量的弹性力学的平衡方程（应力应变）、热传导方程（温度场）、流体力学方程（速度、压力场）。这些方程都是这样的形式：

例如，我们要制作一个装载包覆聚变材料的壳体，就需要设定一个结构方程，以便于控制这个壳体在任何方向上都能有效的承载内部的聚变材料。这样：就要讨论壳体的位移场了

在弹性力学中，位移场和应力场的关系由应变–位移关系式和本构方程确定，例如：

结合材料的弹性常数（杨氏模量 E、泊松比ν）就能得到应力–应变关系式：

这些方程联立后，才能进入工程可用的平衡方程：

但是这样的计算实际上仅仅算了材料上的一个点，没错这是一个连续的偏导数到偏微分再到全微分的推导计算过程。

也只有这种计算我们才可以用尽可能少的泡沫塑料包裹聚变材料，使之不至于影响核弹的直径。

当然了，壳体的材料并非理想刚体，形状也不是完美球体，材料各向异性、温度梯度、装配误差都会影响位移场的分布，就需要在核武器的设计结构上进行多点采样，然后一个点一个点的计算。这也是为什么核武器计算的工程量巨大的原因。

所以，数学学不好造什么核弹啊！想都别想！

但为什么说现在计算方便了呢？这是在计算机系统发展之后引入了有限元分析（FEA，Finite Element Analysis）的计算方法之后产生的计算优势。

通过计算机有序的把一个连续的结构体切割成更小的单元，通过高密度的计算更小单元的特性以达到之前手工计算的结果。这个方法从速度上来说要比人工计算快了几万甚至几十万倍。在这种条件下，量变就会产生质变，之前的核武器计算是依靠经验+推导+测试来进行的，现在是可以在计算空间内不断的试错、对比不同方案找到计算结果的可行性。

在这种方法下武器设计师可以在计算机里“虚拟引爆”成百上千个方案，在计算空间里快速试错、对比不同设计，找到更优解，而不必像过去那样依赖消耗昂贵材料的实爆试验。在完成一部分关键基础数据采集后，就可以借助计算机进行大规模的试算与解析。这样反复迭代，最终就能在有限元分析的结果基础上，收敛出一个大差不差、具备工程可行性的核心设计方案。

这也正是为什么说，今天的核武器设计早已不是什么“绝对的机密”。在具备足够的基础数据后，通过高精度的计算分析，就能推导出可行的设计方案。说白了，核武器的理论和设计方法早在全球范围内就已不再是无法触碰的禁区。

因此，国际核查组织的重点早已不再是拦截所谓的“核武器技术”，而是集中火力监控原料与生产设备——尤其是浓缩铀、武器级钚以及离心机等关键环节。凭心而论，这些机构早在二十多年前就已无力阻止核武器设计技术的扩散，也就只好好去不断收紧对原材料获取和关键设备制造的管控口子。

但矛盾的是——目前核武器真正的技术阶梯其实也比不在材料和生产设备上，而是一个国家的计算能力。例如美国，1992年9月23日，美国进行了最后一次核试验，并于1996年签署了《全面禁止核试验条约》。这代表着美国从1992年9月23日就不进行核武器研究了吗？当然不是！

只不过从那之后，美国的核武器研发重心，已经从实爆试验转移到了大规模计算机模拟。在橡树岭国家实验室等机构，大型计算机集群不断翻新升级，其中有一项最重要的任务，就是进行核武器的全流程模拟计算——从材料的微观物理特性，到内爆的流体动力学过程，再到聚变反应的辐射输运，全部通过超级计算机反复推演。

不仅仅是美国，我们自己也有天河、太湖之光等一系列的超级计算机。

大部分超级计算机的主要用途之一就是用来模拟和验证核武器的设计。

一些小的国家虽然没有超级计算机集群，但同样可以通过数十、数百台服务器组成的计算阵列，完成类似的核武器模拟任务。虽然这样的运算效率比顶级超算差上几千倍，但依然可以在几天甚至一两周内完成一次完整的工程参数验证。

这就是为什么 W 君一直强调：很多国家完全有能力在很短时间内完成核武器的理论设计与工程验证——关键不再是“能不能”，而是“想不想”。

然而，越是在今天核武器设计门槛不断降低、甚至出现“白菜化”趋势的时候，我们反而更应该敬佩当年的先驱者——因为他们不仅没有超级计算机，甚至连成熟的数值模拟软件都没有，却硬是在纸面推导、穿孔卡片机和实验数据之间，把一颗颗原子弹、氢弹做了出来。

邓稼先，在我们的两弹一星项目组，负责的核心工作是工程化可行性计算。他不仅要从理论计算确定内爆方式，还要把这些结果转化为具体的装置结构参数——炸药透镜的几何形状、炸药分块数目、引信位置与延时曲线。这是一个高度复杂的流体力学+爆炸波聚焦计算问题，在当时没有成熟软件，全靠团队编写数值计算程序在机房用穿孔卡片跑结果。邓稼先的贡献，是把这些计算做到可在工厂精确加工制造的程度。

彭桓武，作为理论物理学家，他在核武器工程阶段的突出工作是临界参数计算与结构优化。他推导并计算了不同材料、不同几何形状下的临界质量、临界半径公式，并将其转化为工程上的核装药优化设计方案。这让研制团队在材料利用率和当量提升之间找到最佳平衡，减少了宝贵裂变材料的浪费。

于敏，虽然大家记住的是“构型”，但于敏在工程层面的硬功夫是辐射传输计算与燃料压缩优化。他主导建立了辐射驱动内爆的数学模型，并用有限差分法计算辐射能量在反射器腔体内的分布与时序。这类计算直接决定了二级燃料的压缩速度、均匀性以及压缩后的温度场分布，使得燃料在最短时间内达到点火条件。

郭永怀，他最大的贡献是把空气动力学和冲击波物理的方法引入核武器结构设计。他负责研究高超声速冲击波在炸药、金属壳层、反射器材料中的传播特性，并建立了一整套数值模拟与实验验证体系。这使得理论计算不再是“真空中的理想模型”，而能精准地匹配现实中的材料与结构响应，减少了试验次数加速了我们两弹一星工作的进程。

王淦昌，他带领着带领团队开展了大量中子物理实验与数据测量工作，尤其是精确测定铀-235、钚-239 在不同能量段的中子截面数据。这些一手数据被直接输入到临界质量、反应速率等核动力学计算中，是整个计算体系的物理基础。没有这些实验数据，所有的数值模拟都会停留在纸面假设。

程开甲，这是咱们核试验与数据分析体系的核心人物。他设计并实施了核试验中的数据采集与反演计算流程，从高压传感器布置、X射线闪光摄影，到高速数据采集系统的同步控制，确保能在毫秒甚至微秒级别捕捉爆轰波传播和壳层压缩的真实参数。这些反演数据又反过来修正理论计算，提高了下一轮设计的精度。

对他们的赞美与敬佩，多少字都嫌少。

如果你不研究核技术，可能会以为他们的“牛”，只是因为他们造出了原子弹、氢弹——那只是一个结果。真正让人敬畏的，是他们做到这个结果的过程。

打个比方：当你第一次见到高精度 CNC 机床加工出的严丝合缝的工件，会忍不住惊叹这背后的科技与精度。可他们却能在几十年前，没有 CNC、没有数控、甚至没有现代化检测设备的年代，直接拿出一模一样的东西，然后风轻云淡地说一句：“手搓的。”

没错，这些元勋们就是在没有超级计算机的时候依靠纸、笔、计算尺手工的算出来的！