goodisok 的博客

皮埃尔-西蒙·拉普拉斯 (Pierre-Simon Laplace, 1749-1827)一、出身——从诺曼底农村到巴黎科学院1749年3月23日，拉普拉斯出生在法国诺曼底的博蒙昂诺日（Beaumont-en-Auge）的一个农民家庭。 与牛顿的剑桥精英教育不同，拉普拉斯的教育条件极为有限——他的父亲是个苹果酒农。但他从小展现出惊人的数学天赋。16岁进入卡昂大学，18岁带着一封推荐信前往巴黎，找到了达朗贝尔（d’Alembert，当时法国最权威的数学家）。 面试未果，拉普拉斯写一封信给达朗贝尔，证明了自己的水平。 达朗贝尔读后回答： “我现在不用推荐信就能了解你的能力了。你不需要别的推荐信，你自己就是最好的推荐信。” 拉普拉斯随即被任命为巴黎军事学院的数学教授——当时他只有19岁。 二、拉普拉斯的核心贡献2.1 拉普拉斯方程 ∇2ϕ=0\nabla^2\phi=0∇2ϕ=0——势流理论的数学心脏这是拉普拉斯在天体力学中推导出来的方程（1785年），但它的应用远超出了天文范围。 方程的形式： ∂2ϕ∂x2+∂2ϕ∂y2+∂2ϕ∂z2=0 \frac{\partial^2\p...

莱昂哈德·欧拉 (Leonhard Euler, 1707-1783)一、出生与成长——瑞士教堂里的数学家1707年4月15日，瑞士巴塞尔一个牧师家庭迎来了长子——莱昂哈德·欧拉。 他的父亲原本期望他继承教职。少年欧拉确实在巴塞尔大学学习神学——但他同时也跟着约翰·伯努利（Johann Bernoulli）学数学。伯努利很快就发现这个学生的天分非同一般，劝说欧拉的父亲：”让他学数学。” 欧拉听从了这个建议。 1727年，20岁的欧拉参加了巴黎科学院的悬赏论文竞赛——题目是帆船桅杆的最优布置。他拿了第二名。第一名是他的老师伯努利。伯努利的评论是：”27岁的我教20岁的欧拉——但最多三年他就会超过我。” 结果不到三年就应验了。 1727年的另一件事1727年，牛顿去世。欧拉同年移居圣彼得堡科学院。这个时间衔接颇有象征意味：牛顿开头的时代，由欧拉来接续。 二、惊人的产出：886篇论文欧拉的产出量是人类历史上无与伦比的。 数据 数值 已发表的学术论文 886篇 出版的著作 45本 去世后留下的手稿 够科学院发47年 平均每年发表论文数量 约15-20篇 在...

“我不知道在世人眼中我是什么样的人。但在我自己看来，我不过是一个在海边玩耍的孩童，不时为发现一块更光滑的卵石或更美丽的贝壳而沾沾自喜，而真理的浩瀚海洋尚未被我发现。”—— 艾萨克·牛顿，晚年自述 艾萨克·牛顿 (Isaac Newton, 1643-1727)一、牛顿时空定位——一个不该出生的人1.1 诞生的时间点1643年1月4日，英格兰林肯郡伍尔索普庄园。一个早产儿，小到可以放进一夸脱的壶里。接生婆说这孩子活不了几个小时。 这个婴儿活到了84岁。 他出生的那一年，英国正在进行内战。伽利略在11个月前去世。欧洲还陷入在三十年战争的泥潭中。科学正在经历从”亚里士多德的权威”到”亲自实验”的艰难转型。 1.2 历史的转折在牛顿之前，人们对运动的理解基于亚里士多德（公元前4世纪）的框架： 概念 亚里士多德（~350 BC） 牛顿之前（17世纪中期） 牛顿之后（1687年之后） 运动 需要力才能维持 伽利略提出惯性，未完成 力改变运动，非维持运动 天体运动 圆周运动是”完美”的 开普勒提出椭圆轨道，无物理解释 万有引力解释了开普勒定律 数学描述 几何 解析几...

本文是 John D. Anderson《Fundamentals of Aerodynamics》（第六版）读书笔记系列的第十二篇（完结篇），覆盖教材第13-14章「湍流模型与现代数值方法」的核心内容。从工程CFD的奠基到高保真仿真前沿，为这个系列画上句号。 第十二篇：湍流模型与现代数值方法——从工程 CFD 到高保真仿真一、CFD 的百年演变1.1 概念引入如果你在 1950 年问一位空气动力学家如何计算一个翼型的升力和阻力，答案会是：去风洞。如果你在 2026 年问同样的问题，答案会是：先用 CFD 算一遍，再用风洞验证几个关键工况。 计算流体力学（Computational Fluid Dynamics, CFD） 已经从学术工具变成了工程标准。 1.2 历史节点 年份 事件 影响 1910 Richardson 提出”数值天气预报”概念 CFD 最早的思想 1930s 手算松弛法处理势流 电脑出现前的数值方法 1946 第一台电子计算机 ENIAC 诞生 开始有了”计算”的硬件 1960s NASA Ames 的团队开发隐式格式 可压缩流 N...

本文是 John D. Anderson《Fundamentals of Aerodynamics》（第六版）读书笔记系列的第十一篇，覆盖教材第12章《黏流与边界层》的核心内容。这是理解真实流动物理的关键——解释了为什么势流理论预测的”零阻力”在现实中并不存在。 第十一篇：黏流与边界层——普朗特的开创一、边界层——空气动力学史上最重要的概念1.1 概念引入1904 年，在海德堡举行的第三届国际数学大会上，一位 29 岁的德国学者提出了一个看似简单的想法，却永远改变了空气动力学的面貌。 Ludwig Prandtl 指出： 对于黏性很小的流体（如空气），黏性效应局限在物体表面附近一个很薄的区域内。在这个区域以外，流动可以看作无黏的势流。 这个”薄区域”就是边界层。 1.2 边界层的定义边界层是紧贴在物体壁面上的流体层，其中： 速度梯度很大——从壁面零速度（无滑移条件）变化到外部势流速度 黏性剪切力不可忽略 厚度 δ(x)\delta(x)δ(x) 沿流向增大 边界层厚度 δ\deltaδ 的定义通常为：从壁面到速度达到 0.99Ve0.99V_e0.99Ve​ 的距离...

本文是 John D. Anderson《Fundamentals of Aerodynamics》（第六版）读书笔记系列的第十篇，覆盖教材第11章「超声速翼型设计」的核心内容。阅读本文前需掌握斜激波与膨胀波（第九讲）。 第十篇：超声速翼型设计——激波-膨胀波理论与超临界翼型一、超声速翼型与亚声速翼型的根本区别 1.1 概念引入在亚声速环境下设计翼型时，我们追求薄翼理论中的 Cl=2παC_l = 2\pi\alphaCl​=2πα 和低 CdC_dCd​。但在超声速环境中，翼型设计的哲学被彻底颠覆了。 根本区别： 特性 亚声速翼型 超声速翼型 前缘形状 圆弧形（半径有限） 尖锐（半径很小） 厚度 中（8-15%） 薄（3-8%） 最大厚度位置 30-50% chord 靠近 50% chord 弯度 有弯度更好 弯度几乎不增加升力 激波系统 无（Ma<0.7Ma<0.7Ma<0.7） 前缘+后缘激波、膨胀波 阻力主要来源 摩擦阻力 波阻（占主导） 1.2 超声速翼型的基本类型 类型 截面形状 特点 菱形（双...

本文是 John D. Anderson《Fundamentals of Aerodynamics》（第六版）读书笔记系列的第九篇，覆盖教材第10章「超声速线化理论」的核心内容。阅读本文前需掌握正激波（第八讲）的基本概念。 第九篇：超声速线化理论——从斜激波到普朗特-迈耶膨胀波一、超声速流动的新物种1.1 概念引入当物体以超声速飞行时，它周围的气体行为与亚声速时完全不同： 信息向下游传播（亚声速时信息可以向上游传播） 激波从正前方变成了倾斜的 出现了膨胀波——一种”负激波” 升力和阻力都需要全新的理论框架 1.2 超声速的”信息锥”在超声速流动中，扰动以声速传播，物体比声速快。因此扰动的传播范围局限在一个圆锥内——马赫锥（Mach Cone）。 马赫锥的半顶角： μ=arcsin⁡(1Ma) \mu = \arcsin\left(\frac{1}{Ma}\right) μ=arcsin(Ma1​) Ma=1.2Ma = 1.2Ma=1.2：μ=56.4∘\mu = 56.4^\circμ=56.4∘——宽锥 Ma=2.0Ma = 2.0Ma=2.0：μ=30∘\mu...

本文是 John D. Anderson《Fundamentals of Aerodynamics》（第六版）读书笔记系列的第八篇，覆盖教材第9章《可压缩流基础》。这是全书从低速空气动力学走向高速空气动力学的转折点。 第八篇：可压缩流基础——从亚声速到超声速的入口一、为什么可压缩性在 Ma > 0.3 时必须考虑？1.1 概念引入密度的相对变化 Δρ/ρ∝Ma2\Delta \rho / \rho \propto Ma^2Δρ/ρ∝Ma2。当 Ma=0.3Ma = 0.3Ma=0.3 时，密度变化约 4.5%，这已经是工程上可接受的近似误差上限。 一旦进入 Ma>0.3Ma > 0.3Ma>0.3 的区域，以下现象开始出现： 伯努利方程失效（因为它假设 ρ=const\rho = \text{const}ρ=const） 温度开始显著变化 激波——流动中的一种不连续面——出现 阻力出现一个新的分量：波阻 1.2 热力学基础可压缩流需要引入热力学语言。Anderson 第9章从四个基本热力学参数开始： 内能 e,焓 h=e+pρ,熵 s,比热比 γ...

本文是 John D. Anderson《Fundamentals of Aerodynamics》（第六版）读书笔记系列的第七篇，覆盖教材第8章部分内容和有限翼理论的工程扩展。在掌握了普朗特升力线理论的基础上，现在我们来看真实机翼的复杂性。 第七篇：有限翼实验与修正——从升力线到真实机翼一、升力线理论的局限性 普朗特升力线理论是空气动力学的杰作，但它有几个关键假设限制了适用范围： 小展弦比失效：AR<4AR < 4AR<4 时，升力线理论的精度显著下降 后掠翼不适用：升力线理论假设无后掠 跨声速/超声速失效：可压缩性效应难以处理 非线性大迎角失效：失速附近的流动 Anderson 在第8章中详细讨论了这些局限性以及对应的修正方法。 二、后掠翼的空气动力学2.1 为什么需要后掠翼在 Ma≈0.7−1.2Ma \approx 0.7-1.2Ma≈0.7−1.2 的跨声速飞行中，如果机翼垂直于来流，局部流速会在机翼上表面达到 Ma>1Ma > 1Ma>1，形成激波，带来剧烈的阻力增长（波阻）。 后掠翼的原理：只有垂直于前缘的...

本文是 John D. Anderson《Fundamentals of Aerodynamics》（第六版）读书笔记系列的第六篇，覆盖教材第7章《有限翼理论》（Prandtl’s Lifting-Line Theory）。这是Anderson全书中最具工程实用价值的一章之一。 第六篇：有限翼理论与升力线——从二维到三维的跨越一、为什么二维理论不够？ 1.1 概念引入薄翼理论给出了一个完美的二维翼型：Cl=2παC_l = 2\pi\alphaCl​=2πα，CdC_dCd​ 很小，唯一的力是升力。但现实中的机翼是有限展长的——有两个端部（翼尖），空气可以从下表面通过翼尖流到上表面。 这个”翼端泄漏”效应会带来三个重要后果： 升力降低：在翼尖处升力降为零 诱导阻力出现：额外的阻力，在二维中不存在 翼尖涡形成：真实机翼后面拖出两条集中的涡 Anderson 将普朗特（Prandtl）在 1918-1919 年左右提出的升力线理论描述为”空气动力学史上最优雅的理论之一”。 1.2 物理机制升力沿着展向分布——翼根处大，翼尖处小。展向升力梯度意味着环量沿着展向变化。 根据亥...