14 Apr

流体静力平衡的应用

很早以前我就对这个问题感兴趣了,但是一直搁置着,没有怎么研究。最近在阅读《引力与时空》的“潮汐力”那一节时重新回到了这个问题上,决定写点什么东西。在这里不深究流体静力平衡的定义,顾名思义地理解,它就是流体在某个特定的力场下所达到的平衡状态。流体静力学告诉我们:

达到流体静力平衡时,流体的面必定是一个等势面。

这是为什么呢?我们从数学的角度来简单分析一下:只考虑二维情况,假如等势面方程是$U(x,y)=C$,那么两边微分就有
$$0=dU=\frac{\partial U}{\partial x}dx+\frac{\partial U}{\partial y}dy=(\frac{\partial U}{\partial x},\frac{\partial U}{\partial y})\cdot (dx,dy)$$

这意味着向量$(\frac{\partial U}{\partial x},\frac{\partial U}{\partial y})$和向量$(dx,dy)$是垂直的,前者便是力的函数,后者就是一个切向量(三维就是一个切平面)。也就是说合外力必然和流体面垂直,这样才能提供一个相等的方向相反的内力让整个结构体系处于平衡状态!

点击阅读全文...

14 Apr

费曼积分法(8):求高斯积分

自从了解了费曼积分法之后,我就一直想着用费曼积分法来求高斯积分$\int_0^{\infty} e^{-x^2}dx=\frac{\sqrt{\pi}}{2}$这个神奇的积分,但一直无果。在《数学桥》里边,作者是通过将其转变为二重积分来解决的,简洁而巧妙。但是为了显示费曼积分法的威力,我一直想找到高斯积分的其他求法。上星期在《数学物理方法》中看到作者用拉普拉斯变换求出了该积分,眼睛不禁为之一亮,不过这属于积分变换内容,属于“积分符号内取积分”的技巧,在此不作讨论。今天在网上查找资料时,在“赵洁”的一篇论文《含参变量积分》中,看到了一种属于费曼积分法范畴内的方法,特与大家分享。

从“事后分析”来看,高斯积分的结果涉及到了$\sqrt{\pi}$这个量,一般来说我们常见的公式出现$\pi$的不少,可是几乎没有出现$\sqrt{\pi}$的,所以一般来说我们都将它平方。我们引入
$$f(x)=(\int_0^x e^{-t^2}dt)^2$$

点击阅读全文...

18 Apr

纠缠的时空(三):长度收缩和时间延缓

我们之前通过矩阵变换方式推导出了洛伦兹变换以及速度合成公式等结论,不得不说,矩阵推导方式有种引人入胜的魅力。今天,在讲述相对论(包括电动力学、广义相对论)的书籍里边,在数学形式上取而代之了张量这一工具,这实际上是对矩阵的一个推广(之前已经提到过,二阶张量相当于矩阵)。采用这样的形式在于它充分体现了相对论的对称和变换关系。本文将来谈及狭义相对论的一些基本结论,包括同时性、长度收缩、时间延缓等。

本文的光速$c=1$。

同时的相对性

在同一时空中,采取两个时空坐标进行洛伦兹变换,再作差,我们得到:
\begin{equation}\left[\begin{array}{c} \Delta x\\ \Delta t \end{array}\right]=\frac{1}{\sqrt{1-v^2}}\left[\begin{array}{c c}1 & v\\ v & 1 \end{array}\right]\left[\begin{array}{c}\Delta x'\\ \Delta t' \end{array}\right]\end{equation}

点击阅读全文...

25 Apr

学习场论(电磁场、重力场)

本博客的文章其实一定程度上反映了我在该时期的学习研究,所以我觉得写blog是一件很惬意的事情,它记录着我的成长历程。读者可能留意到,我上学期说对量子力学很感兴趣,也算是入了一点点门。这学期开学初表示对摄动理论方面的知识很感兴趣,也研究了一两个星期。再后来就将学习重点放在了相对论上面了。现在呢?我在学习朗道的《场论》,主要先学习电磁场(电动力学)。

有的读者可能比较无语:你怎么变来变去,学习不是贵在精而不在多吗?

点击阅读全文...

22 May

当Matlab遇上牛顿法

牛顿法是求方程近似根的一个相当有用而且快捷的方法,我们最近科学计算软件课程(Matlab)的一个作业就是编写求方程近似解的程序,其中涉及到牛顿法。我们要实现的目标是,用户输入一道方程,脚本就自动求出根来。这看起来是一个挺简单的循环迭代程序,但是由于Matlab本身的特殊性,却产生了不少困难。

Matlab是为了数值计算(尤其是矩阵运算)而生的,因此它并不擅长处理符号计算。这就给我们编程带来了困难。在网上随便一搜,就可以发现,网上的Matlab牛顿法程序都是要求用户同时输入方程及其导函数,这显然是不方便的,因为Matlab本身就具备了求导功能。下面我们来分析一下困难在哪里。

我们要实现的最基本功能是定义一个函数,然后可以根据该函数求具体的函数值,并且自动求该函数的导数,接着求导数值。这些看起来很基本的功能在Matlab中却很难调和,因为Matlab的“函数”定义很广,一个具有特定功能的M文件叫“函数”,一个运算式$f(x)$也可能是一个函数,显然后者是可以求导的,前者却不行,所以Matlab一刀砍——不能对函数求导!!

点击阅读全文...

11 May

电的相对论效应——磁“子虚乌有”?

也许大家会觉得,相对论中有一个因子
$$\gamma =\frac{1}{\sqrt{1-\frac{v^2}{c^2}}}$$
因此,相对论的效应只有在高速情况下,即v比较接近于c的情况下才会凸显出来。这在一般情况下是正确的,但是却不全对。因为存在相当明显的、速度低于1mm/s的相对论效应——那就是几乎人尽皆知的“磁”。

之前已经提及过,磁场可以解释为电场的相对论效应,因此所有电磁现象都可以归因为电场和相对论。事实上,这是正确的,只是教科书上并没有明确说出这一点而已。于是我们就不难理解“为什么电磁学的麦克斯韦方程组会与相对论协调”、“为什么电场与磁场的表现如此相似”等等问题了,因为它们的探究本身就在相对论的框架下,磁场和电场都是一个东西的结果。

点击阅读全文...

7 Jun

《虚拟的实在(2)》——为什么引力如此复杂?

上一篇文章里我已经从我自己的理解角度简单说了一下场论的必要性,这次让我们再次谈到这个话题,企图在文字层面上得到更深入的认识。

上一两周的时间,我一直在找资料,主要是线性引力的资料,并且发现了很多有趣的东西,在此一并与大家分享一下。首先,当我在Google中输入“线性引力”时,我发现了一本“奇书”,一本名副其实的“巨著”——《引力论》!洋洋1300多页的大作,三位“超级巨星”——C.W.麦思纳(Charles W.Misner)、K.S.索恩(Kip S.Thorne)、J.A.惠勒(John Archibald Wheeler)——联合编写,恐怕再也找不到哪本书可以PK它的“全明星阵容”了。该书英文名为Gravitation,中文是由台湾翻译的,繁体中文版。全书讲述了引力的研究历史和发展情况,更重要的是几乎每一处历史都给出了数学论证!最最重要的,作者惠勒还是跟爱因斯坦同一个研究时代的人,我们可以最真实的感受到那年代的研究。看到这里,我就迫不及待地想买了,由于各种原因,我们很难买到,到图书馆找,发现有英文版的,就马上借过来了,另外因为买不到中文版,我只好到网上买了电子版,然后打印出来了。不过不是很清晰,而且自我感觉中文翻译不是很好(当然,已经够我们阅读了)。

点击阅读全文...

1 Jul

与首都机场的“难分难舍”

无尽的等待

无尽的等待

上个月的最后三天(06.28-06.30),我去国家天文台参加了第三届宇宙驿站的站长联谊会及科普研讨会。会议在河北兴隆天文台举行,我们按照计划是先到北京总部,然后去兴隆,然后回到北京总部解散。具体的故事我会另写文章与大家分享,本文主要想说一下我与北京首都国际机场的“难分难舍”的返程之旅。

按照计划,我是昨晚9点的飞机,今天凌晨应该可以到广州。我七点多到机场,八点左右就办完了登记手续,然而,我们等了两三个钟,最终得到的结果是:由于雷暴雨的影响(北京并没有下雨,估计是途中某个地方的上空天气太糟糕),该航班取消,补到第二天七点......这对我来说可真是个大考验。虽说航空公司会为我们联系宾馆,但是效率之低让不少人在机场抗议,于是乎冰冷的机场一下子就热闹起来的(取消的不知我们一趟航班,还有很多其他航班)。而我虽然来过好几次北京,毕竟还属于“异客”,自然经验不足,但我做出了一个很大胆的决定:在机场过夜!

点击阅读全文...