冯康(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 年国家科技进步奖二等奖。
四、代表著作与文献
- 《基于变分原理的差分格式》(《应用数学与计算数学》,1965):独立创始有限元方法的标志性论文。
- 《数值计算方法》(合著,1970 年代):中国计算数学的早期教材。
- 《哈密顿系统的辛算法》(论文集):辛几何算法的系列工作。
- 论文《差分格式与辛几何》(1980 年代):辛算法领域的奠基性论文之一。
- 《冯康先生纪念文集》(袁亚湘编,科学出版社,2020):生平、贡献与学术传承。
- 传记《冯康传》(宁肯、汤涛著,浙江教育出版社,2019,ISBN 9787553698496)。
五、师承、学派与同时代交往
- 师承:中央大学物理系;1946 年起在清华大学受教;博士阶段(存疑:冯康未取得正式博士学位,苏联斯捷克洛夫研究所时期受庞特里亚金指导)。
- 庞特里亚金(Lev Pontryagin):在莫斯科的导师,把他带入拓扑群与李群的世界,这一纯数学训练为其后来的辛几何工作埋下伏笔。
- 学派归属:中国计算数学学派的奠基人与领导者;国际上是有限元方法的独立奠基人之一、几何数值积分(保结构算法)的开创者。
- 学生与合作者:带领秦孟兆等学者发展了哈密顿系统辛算法;培养了中国几代计算数学人才。
- 国际交往:1978 年后在法国、意大利、日本、美国等十余所大学讲学,使中国计算数学走向世界;与 Peter Lax 等国际计算数学大师有交往,Lax 为其撰写了权威讣告。
六、争议、公案与历史评价
- 有限元方法的优先权:冯康在 1950 年代末—1960 年代初独立发展有限元方法,但因当时中国与西方学术隔绝、论文以中文发表,其贡献在西方文献中长期未获承认。1980 年代后,经 Peter Lax 等国际权威学者的确认,冯康的独立创始地位才得到国际公认。这是 20 世纪「多重独立发现」的典型案例,并非学术争议,但涉及历史叙述的公平性。
- 辛算法的归属:辛积分器的个别格式(如 Verlet 法)在冯康之前已在物理中使用,因此辛算法的「优先权」表述需谨慎:准确的说法是,冯康是第一个系统建立辛几何算法理论、提出通用构造方法并证明其长期优势的学者,开创了「保结构算法」领域,这一点在国际上无争议。
- 历史评价:冯康是中国计算数学无可争议的奠基人。他在两个相隔二十年的方向(有限元、辛算法)上都做出了奠基性贡献,且每个方向都成为科学计算的主流。他以工程需求为起点、以严格数学理论为支撑的研究风格,影响了整个中国计算数学界。
七、身后影响与学术谱系
- 有限元方法:已成为工程仿真(ANSYS、ABAQUS 等商业软件)的核心,覆盖结构力学、流体、电磁、热传导等几乎所有工程领域;中国计算力学的建立直接源于冯康的工作。
- 保结构算法 / 几何数值积分:辛算法是分子动力学(GROMACS、NAMD)、天体力学长期轨道计算的标准方法;保结构思想已推广到李群积分器、体积保持、接触系统等,形成「几何数值积分」这一独立学科方向。
- 中国计算数学学科:创办中科院计算中心(后发展为计算数学与科学工程计算研究所)与科学与工程计算国家重点实验室,培养了几代人才,使中国计算数学跻身世界前列。
- 命名遗存:冯康科学计算奖(CAS/LSEC 设立);冯康楼(中科院)。
- 对国际计算数学的影响:冯康是国际计算数学界公认的大师,其工作被写入 Ciarlet《The Finite Element Method》、Leimkuhler–Reich《Simulating Hamiltonian Dynamics》、Hairer–Lubich–Wanner《Geometric Numerical Integration》等经典著作。
八、延伸阅读与权威来源
- 中国科学院数学与系统科学研究院:冯康人物介绍:中科院官方传记,本档案主要史实出处(生平、研究领域、奖项)。
- LSEC: 冯康纪念页 Home page of later Professor Feng Kang:计算中心撰写的纪念文章,含各阶段贡献细节。
- Wikipedia: Feng Kang:生平、教育与贡献的交叉核对。
- Peter Lax, "Feng Kang", SIAM News, 1993:国际计算数学大师 Peter Lax 撰写的权威讣告,确认有限元独立创始地位。
- 冯康科学计算奖官网:以冯康命名的奖项信息。
- Wikipedia: Finite element method:有限元法的西方发展脉络,用于对比冯康的独立工作。
- Wikipedia: Symplectic integrator:辛积分器与保结构算法的技术背景。
- Wikipedia: Boundary element method:边界元法,含自然边界元的地位。
- 传记《冯康传》(宁肯、汤涛,浙江教育出版社,2019):权威中文传记。