← 返回 冯康 简介档案 深度档案
计算数学数值分析

冯康(Feng Kang)

中国计算数学的开创者——他在 1960 年代与西方几乎同时、却完全独立地发明了有限元方法;在 1980 年代又单枪匹马提出哈密顿系统的「辛几何算法」,开创了「保结构算法」这一新领域,让长期数值模拟不再因能量漂移而失真。从三峡大坝的应力计算到分子动力学、天体力学的长时间轨道追踪,其方法已渗入现代科学计算的每一根血管。

一、生平与时代背景

冯康 1920 年 9 月(一说 9 月 9 日)生于江苏南京,童年在苏州度过,就读于苏州中学。1939 年春,他以同等学历考入福建协和学院数理系;同年秋又重新考入中央大学(今南京大学)电机工程系,两年后转入物理系。他在大学期间兼修了电机、物理、数学三系的主要课程,这一宽厚的跨学科基础对他日后的研究起了独特作用。1944 年毕业后,他患了脊椎结核,在家休养期间坚持自学数学,方向逐步转向基础数学。

1945 年,冯康任复旦大学数学物理系助教;1946 年到清华大学任物理系助教,1951 年转任数学系助教。1951 年他调入中国科学院数学研究所任助理研究员;1951—1953 年赴苏联莫斯科斯捷克洛夫数学研究所进修,师从拓扑学大师列夫·庞特里亚金(Lev Pontryagin),研究拓扑群与李群。1957 年他调入新成立的中国科学院计算技术研究所,从此转向应用数学与计算数学,成为中国计算数学与科学计算的奠基人和领导者。

1978 年,中科院计算中心成立,冯康出任首任主任(至 1987 年任名誉主任)。1980 年当选中国科学院学部委员(院士)。1990 年他倡导并亲自筹备成立了中国科学院科学与工程计算国家重点实验室。1993 年 8 月 17 日,冯康病逝,享年 72 岁。他去世后,中科院计算数学与科学工程计算研究所于 1994 年设立「冯康科学计算奖」(Feng Kang Prize),1995 年首次颁发,用以表彰 45 岁以下在科学计算领域做出突出贡献的海内外华人学者。

冯康的研究横跨拓扑群、广义函数、计算数学与科学工程计算。他主持的「有限元方法」获 1982 年国家自然科学奖二等奖,「地震勘探数值方法」获 1987 年国家科技进步奖二等奖,「哈密尔顿系统的辛几何算法」获 1997 年国家自然科学奖一等奖(追授)。

二、学术生涯年表

年份 事件
1920 9 月(一说 9 月 9 日)生于江苏南京,童年在苏州
1939 春 考入福建协和学院数理系(存疑:维基百科仅记中央大学,此处据中科院 AMSS 官方介绍)
1939 秋 重新考入中央大学电机工程系
1941 转入中央大学物理系
1944 物理系毕业;后患脊椎结核,在家自学数学
1945 任复旦大学数学物理系助教
1946 任清华大学物理系助教
1951 转清华数学系助教;同年调中科院数学所助理研究员
1951—1953 赴苏联莫斯科斯捷克洛夫数学研究所进修,师从庞特里亚金
1957 调中科院计算技术研究所,转向计算数学
1950 年代末—1960 年代 提出「基于变分原理的差分格式」,独立创始有限元方法(1965 年发表)
1978 任中科院计算中心首任主任
1980 当选中国科学院学部委员(院士)
1982 「有限元方法」获国家自然科学奖二等奖
1984 起 转向动力系统,提出哈密顿系统的辛几何算法
1987 「地震勘探数值方法」获国家科技进步奖二等奖;任计算中心名誉主任
1990 倡导成立中科院科学与工程计算国家重点实验室
1993 8 月 17 日病逝,享年 72 岁
1994 中科院设立「冯康科学计算奖」
1997 「哈密顿系统的辛几何算法」获国家自然科学奖一等奖(追授)

三、核心成果(逐条技术解析)

1. 基于变分原理的差分格式——独立创立有限元方法(1950 年代末—1965)

问题背景:1950 年代中国大规模开展水利工程(如大坝建设),需要计算弹性结构的应力分布,归结为求解椭圆型偏微分方程的边值问题。传统差分方法在复杂几何与不规则网格上难以处理。

关键思想:冯康没有走直接差分离散的老路,而是从变分原理出发:把椭圆边值问题写成弱形式(虚功原理 / 伽辽金形式)——求 $u\in H^1_0(\Omega)$ 使得

$$a(u,v):=\int_\Omega \nabla u\cdot\nabla v\,dx=\int_\Omega fv\,dx=:(f,v),\quad \forall v\in H^1_0(\Omega),$$

然后在计算域上作三角剖分,取分片多项式(通常分片线性)构成的有限维子空间 $V_h\subset H^1_0$,求解离散问题 $a(u_h,v_h)=(f,v_h),\ \forall v_h\in V_h$。他证明了收敛性与误差估计,并指出由此导出的计算格式等价于一种差分格式——故他称之为「基于变分原理的差分格式」。

结果:这正是西方所称的有限元方法(finite element method, FEM)。冯康在 1950 年代末到 1960 年代初、几乎与西方(Turner–Clough–Martin–Topp 1956 年的波音刚度法、Argyris、Zienkiewicz 等人的工作)同时但完全独立地发展了这一方法,1965 年在《应用数学与计算数学》发表论文《基于变分原理的差分格式》。

与前人关系:变分思想可追溯到 Courant 1943 年;工程界的刚度法由 Turner 等(1956)提出,Zienkiewicz 在 1960 年前后推广了「有限元」这一名称。冯康的独立工作证明了该方法的数学收敛性,并给出了系统的理论框架。Peter Lax 在冯康去世后的 SIAM News 讣告中确认了其独立创始地位。

意义:有限元方法被公认为 20 世纪计算数学的里程碑之一,是固体力学、流体力学、电磁场等工程仿真的核心工具。冯康因此被视为有限元法的奠基人之一。

局限性:当时中国与西方学术隔绝,冯康的中文论文未及时为国际所知,导致其优先权在西方文献中长期未被充分承认;直到 1980 年代后才随国际交流恢复而获得公认。

2. 自然边界归化与自然边界元方法(1970 年代)

问题背景:对无界区域(如外部问题)的椭圆方程,有限元需要截断区域并施加人工边界条件,误差较大。

关键思想:冯康系统研究把椭圆型方程边值问题归化为边界积分方程(自然边界归化),发展了自然边界元方法。其核心是通过格林函数把区域上的微分方程转化为边界上的积分方程,从而降低维数、精确处理无界区域。

结果:他提出的「自然边界归化」与经典的直接边界归化、间接边界归化并列为边界元方法的三大主要途径;他还发展了有限元与边界元的自然耦合方法,兼具两者优点。

3. 间断有限元空间中的嵌入理论(1970 年代)

冯康在 1970 年代研究了间断有限元空间中的嵌入(embedding)理论,把经典的椭圆型偏微分方程理论推广到不同维数的组合情形,为弹性组合结构(如板—梁、壳—加强筋耦合)提供了严格的数学基础。

4. 哈密顿系统的辛几何算法——开创保结构算法(1984 起)

问题背景:哈密顿系统 $\dot z=J\nabla H(z)$($J$ 为标准辛矩阵)描述了天体力学、分子动力学、经典力学等大量守恒系统,其精确流保持相空间的辛结构(辛 2-形式 $\omega=dp\wedge dq$,即刘维尔定理)。但传统数值积分(如普通 Runge–Kutta 法)不保持辛结构,长时间模拟会出现能量单调漂移、轨道扭曲,使长期预测失效。

关键思想:冯康 1984 年提出,数值积分器本身应当是一个辛变换(保持辛 2-形式的正则变换),这样才能在长时间计算中保持系统的几何结构。他系统构造了辛算法(symplectic algorithm),如辛 Euler 法、Störmer–Verlet 法(蛙跳法)等,并证明它们在长时间跟踪中远优于传统算法。

精确表述:一个一步法 $z_{n+1}=\Phi_h(z_n)$ 称为辛的,若其雅可比满足 $\Phi_h'^T J\Phi_h'=J$(即保持辛形式)。辛算法保持相空间体积与能量的有界振荡(而非漂移),能正确再现哈密顿系统的不变环面(KAM 环面)。

推广:冯康带领研究组把思想从有限维哈密顿系统推广到无穷维哈密顿 PDE(如波动方程、非线性薛定谔方程),并进一步推广到具有李代数结构的动力系统——切触系统(contact systems)、无源系统(divergence-free systems)等,形成「结构保持算法」(structure-preserving algorithms / geometric numerical integration)的统一框架。

意义:辛几何算法开创了几何数值积分这一活跃领域,广泛应用于天体力学(太阳系长期演化模拟)、分子动力学(蛋白质折叠模拟)、等离子体物理等。「哈密顿系统的辛几何算法」获 1997 年国家自然科学奖一等奖(追授)。

与前人关系:辛积分器的个别例子(如 Verlet 法)在物理中早已使用,但冯康是第一个系统地从辛几何观点建立理论、提出构造方法并证明其长期优势的人,被国际公认为该领域的开创者。

5. 地震勘探数值方法

冯康将有限元与波动方程数值方法应用于地震勘探,发展了地震波传播的数值模拟技术,为中国油气勘探提供了关键计算工具,获 1987 年国家科技进步奖二等奖。

四、代表著作与文献

五、师承、学派与同时代交往

六、争议、公案与历史评价

七、身后影响与学术谱系

八、延伸阅读与权威来源