本文由Carl E. Mungan（美国海军学院物理系）和Garth A. Sheldon-Coulson（美国波特兰市）两位作者合作完成，发表于欧洲物理学会期刊《european journal of physics》2021年第42卷，编号025008。这是一篇研究流体力学中非稳态流动问题的学术论文，具体探讨了变截面U型管中液柱振荡的行为。

一、 研究背景与目标

研究的核心科学领域是流体动力学，特别是非稳态伯努利方程的应用。在基础物理教学中，U型管中液体的振荡是一个经典问题，但通常假设管子截面均匀。对于截面变化的U型管，其动力学行为更为复杂，超出了标准伯努利方程的适用范围。本研究旨在填补这一空白，深入分析液体在截面平滑变化的垂直U型管中振荡时的运动方程和压力分布。

研究的学术动机在于展示非稳态伯努利方程在处理此类非定常流动（流场中任一点速度随时间变化）问题上的有效性和正确性。作者指出，此前虽有文献通过拉格朗日方法分析了左右臂截面不同但各自均匀的U型管，但本研究采用非稳态伯努利方程处理截面连续平滑变化的更一般情形，并验证其结果与机械能守恒的一致性。该研究的目的是为高年级本科生流体力学课程提供一个深刻且实用的教学案例，帮助他们超越稳态流假设，理解流体加速流动时的物理本质。

二、 详细研究方法与流程

本研究是一个理论分析与数值计算相结合的工作，不涉及实验室实验，其主要“研究对象”是所建立的数学模型。研究流程可概括为以下几个关键步骤：

步骤一：建立物理模型与坐标系 研究者考虑一个两端开口、垂直放置的U型管，其横截面积A[s]是沿管轴中心线坐标s的平滑函数（s沿逆时针方向绕管轴测量）。管中装有不可压缩、无粘性、无旋的理想液体。在静态平衡时，左右两个自由液面（A和B）位于同一高度（y=0），液柱总长度为L。当液柱发生振荡时，左液面A的垂直位移记为y_A（向下为正），右液面B的位移记为y_B（向上为正）。假设振荡幅度足够小，自由液面始终位于U型管的垂直段内。任一轴线上点P的压力和速度记为p和v。

步骤二：应用非稳态伯努利方程推导运动方程 这是本研究的关键步骤。作者使用了附录中推导的非稳态伯努利方程（公式A9）。该方程在稳态伯努利方程（压力能+动能+势能守恒）的基础上，增加了一项与流动加速度相关的项，其形式为惰性乘以体积流量对时间的导数。 具体应用时，在左右自由液面A和B之间应用此方程（公式2）。由于液体不可压缩，体积流量Q = v A 在整个液柱的任一截面上瞬时相等。在推导对时间导数d(vA)/dt时，需要特别注意，因为即使管的截面积A不随时间显式变化，但跟踪自由液面等随流体运动的质点时，A是时间的隐函数。作者给出了正确的展开式（公式3），并指出忽略其中一项会导致错误。 结合连续方程（v_B = v_A * A[y_A]/A[y_B]）和几何约束（左右臂液体体积变化相等，公式5），最终推导出关于左液面位移y_A[t]的一个积分-微分方程（公式4）。这个方程描述了y_A与其时间导数（速度、加速度）、右液面位移y_B（通过体积约束与y_A关联）以及管子截面积函数A[s]之间的关系。

步骤三：通过能量守恒验证运动方程 为了验证从非稳态伯努利方程推导出的运动方程（公式4）的正确性，研究者采用了另一种独立方法：机械能守恒。他们计算了振荡液柱相对于平衡位置的势能变化（公式6）和整个液柱的动能（公式7，利用连续性将速度v[s]用v_A和面积比表示）。总机械能（U+K）对时间的导数为零。运用莱布尼茨积分法则对含时上下限的积分求导（公式9），分别计算了dU/dt（公式10）和dK/dt（公式11）。将两者相加并设为零，经过化简，最终得到了与公式4完全相同的方程。这一交叉验证确认了理论推导的严谨性。

步骤四：均匀截面U型管的解析解 作为特例和验证，研究者将模型应用于横截面积A为常数的U型管。此时，运动方程大幅简化（公式12），可直接求解。他们得到了左液面（或右液面）做简谐运动的解析解（公式13），其角频率ω = √(2g/L)（公式14）。此外，他们还利用非稳态伯努利方程计算了振荡过程中U型管垂直段内任意深度h处的压力，得到了公式16：p - p_atm = ρ(g + a)h。其中a是该点液体的垂直加速度。这个结果与在垂直加速的电梯中水杯内的压力公式一致，是牛顿第二定律在流体中的体现，而非通常的静水压力公式。这一点纠正了使用稳态伯努利方程会得出的错误结论。

步骤五：变截面U型管的数值解——以线性渐变管为例 对于更一般的变截面情况，运动方程（公式4）通常没有解析解，需要进行数值求解。研究者选择了一个具体例子进行演示：一个截面为圆形、半径r随s线性增加的U型管。他们假设半径r = (s + H)t，其中t是渐变率，H是左臂允许液面上升的最大高度参数。 在此假设下，面积函数A[s] = πr^2已知。体积约束方程（公式5）可以化为一个关于y_B和y_A的三次方程，并求出了显式解（公式19）。同时，惰性积分（公式4右侧积分）也可以求出解析表达式（公式20）。将这些表达式代入运动方程（公式4），就得到了一个关于y_A[t]的常微分方程。 研究者使用Mathematica软件对该方程进行数值求解。他们设定了一套参数（g=9.8 m/s², L=0.2 m, H=0.1 m, t=0.1），并考虑了两种初始条件：1）小振幅（初始将左液面下推1毫米后释放）；2）大振幅（初始将左液面上拉8厘米后释放）。通过NDSolve命令求解方程，并绘制了左右液面位移随时间变化的曲线（图2和图3）。

三、 主要研究结果

1. 理论方程的建立与验证：成功推导出描述变截面U型管内液面振荡的积分-微分方程（公式4），并利用能量守恒原理完成了严谨的交叉验证，确保了理论框架的自洽性。
2. 均匀管的结果：确认了在均匀截面U型管中，液面做严格的简谐振动，角频率由√(2g/L)给出，与经典结论一致。更重要的是，明确了振荡过程中液体内部压力的正确分布由p - p_atm = ρ(g + a)h给出，其中a是局部加速度，这澄清了常见的概念误区。
3. 变截面管的数值结果：
   • 小振幅振荡：对于线性渐变的U型管，在微小扰动下（初始位移1 mm），液面的振荡（图2蓝色曲线）非常接近正弦曲线（图2绿色虚线）。其振荡频率比均匀管的固有频率（约14 rad/s）高出约29%，说明截面变化影响了系统的等效“劲度”或惯性。
   • 大振幅振荡：当初始位移很大时（初始上拉8 cm，接近左臂极限H），振荡表现出明显的非线性（图3）。位移-时间曲线不再是完美的正弦波，波形发生畸变。这表明在变截面系统中，大幅振荡的运动方程本质上是非线性的，无法用简单的简谐运动来描述。右液面的位移（图2、图3中红色曲线，纵坐标放大了5倍）形状与左液面不同，直观展示了由于左右臂截面积不同导致的运动不对称性。

四、 研究结论与价值

本研究的主要结论是：非稳态伯努利方程是分析U型管等系统中液体非定常流动的有效理论工具。它能够准确预测变截面U型管中液面的振荡运动，并给出振荡过程中液体内部正确的压力分布，而传统的稳态伯努利方程在此类动态问题中不适用。 其科学价值在于： 1. 理论层面：清晰展示了非稳态伯努利方程的推导与应用，强调了在计算流量时间导数时需考虑面积随流体质点运动而变化的隐式时间依赖性，这一点容易被忽略。 2. 教学层面：为本科生流体力学课程提供了一个内容深刻、逻辑完整、兼具解析与数值方法的优秀教学案例。它引导学生从简单的均匀管解析解过渡到复杂的变截面管数值解，理解线性与非线性振荡的区别，并深化对伯努利方程适用条件的认识。 3. 方法论层面：演示了如何将理论建模（非稳态伯努利方程、能量守恒）、解析求解（特例）和数值计算（一般情况）相结合来解决物理问题。

五、 研究亮点

1. 研究方法的创新性与严谨性：采用非稳态伯努利方程处理变截面U型管振荡这一非标准问题，并罕见地通过机械能守恒对推导出的运动方程进行了严格的交叉验证，体现了理论工作的严谨性。
2. 对经典概念的澄清与深化：明确指出并纠正了在动态振荡问题中滥用静水压力公式或稳态伯努利方程的常见错误，得出了符合牛顿第二定律的压力分布公式，加深了对流体动态行为的理解。
3. 从线性到非线性的完整展示：研究不仅给出了均匀管线性简谐运动的解析解作为基准，更重要的是通过数值方法探讨了变截面管非线性振荡的特性（波形畸变），展示了物理系统的复杂性。
4. 详尽的附录推导：论文附录提供了非稳态伯努利方程的完整推导过程，从牛顿第二定律应用于流体微元出发，清晰引出了“惰性”这一概念，使其不仅是一篇研究报告，也是一个自包含的教学资料。

六、 其他有价值内容

论文的附录部分本身具有很高的教学价值。它从流体微元受力分析出发，系统推导了非稳态伯努利方程，并解释了“惰性”概念的来源——类似于电路中的电感。这有助于学生理解该方程与稳态方程的本质区别，以及其物理图像的来源。 此外，研究者在讨论中提及了相关但不同的问题（如潜艇注水下潜），并将自己的结果与已知的物理现象（加速电梯中的流体压力）进行类比，体现了知识的联系与迁移。文中还指出了其他文献中曾出现的错误或疏忽，例如在计算d(vA)/dt时漏项的问题，这对科研工作者和学生都具有警示意义。