Minnaert 假设气泡关闭时发生周期性膨胀和收缩,周围水也跟着振动,就嗷地一声叫了出来!设有个半径为 \(r\) 的泡形成后开始简谐振动,半径有

\[r=r_0+a\sin\frac{2\pi t}T
\]

下面考虑势能和动能的守恒,先算势能。当气泡半径减小 \(x\) 时,体积减小的比例是 \(\left(\dfrac{r-x}r\right)^3\)。假设瞬间的体积压缩是绝热过程(合理,毕竟嗷的一瞬间传热哪有那么快的),则压力满足

\[\frac p{p_0}=\frac{v_0^k}{v^k}=\left(\frac r{r-x}\right)^{3\kappa},\frac{c_p}{c_v}=\kappa
\]

假设气泡的脉动不太大(合理,气泡也不可能一瞬间变多小),则可以直接一阶近似

\[p-p_0=\frac{3\kappa px}r
\]

在压缩到最小体积下的势能是:

\[-\int_{V_0}^V\left(p-p_0\right)\mathrm dV=\int_0^a\frac{3\kappa px}r4\pi r^2\mathrm dx=6\pi\kappa pra^2
\]

再算动能。谐振子的动能由水的运动决定,假设水的运动是径向的。气泡壁的速度为:

\[\frac{\mathrm dr}{\mathrm dt}=\frac{2\pi a}T\cos\frac{2\pi T}T
\]

距离 \(R\) 处的一个小水团子的速度为

\[\frac{r^2}{R^2}\frac{2\pi a}T\cos\frac{2\pi T}T
\]

气泡脉动时所有水团子的总动能我们直接积分只因到无穷远。

\[\frac 12\int\left(\frac{\mathrm dr}{\mathrm dt}\right)_{\max}^2\mathrm dm=\frac\rho2\int_r^\infty\left(\frac{r^2}{R^2}\frac{2\pi a}T\right)^24\pi R^2\mathrm dR=\frac{8\pi^3\rho r^3a^2}{T^2}
\]

简谐运动让动能等于势能就有

\[6\pi\kappa pra^2=\frac{8\pi^3\rho r^3a^2}{T^2}
\]

算出气泡的脉动周期和频率,即我们要求的 Minnaert 共振频率!

\[T^2=\frac{4\pi^2\rho r^2}{3\kappa p},f=\frac 1T=\frac1{2\pi r}\sqrt\frac{3\kappa p}\rho
\]

代入室温下的气体数据

\[f=\frac1{2\pi r}\sqrt\frac{3\times1.400\times100000}{1000}\approx\frac{3.26\text{(m/s)}}r
\]

Minnaert的更多相关文章

  1. Opengles 管线编程介绍

      OpenGL ES 2.0可编程管道 上图橙色部分(Vertex Shader和Fragment Shader)为此管道的可编程部分.整个管道包含以下两个规范: 1)         OpenGL ...

随机推荐

  1. 实验11.ACL实验

    # 实验11.ACL实验 本实验用于测试ACL,类似于防火墙. 拓扑 要求阻塞PC1到PC2和server的全部协议,阻塞client1到server1的icmp协议 具体配置 首先利用ospf协议实 ...

  2. Python str 转 b’二进制串

    用raw_unique_escape来编码无义意的二进制串 bytes(strtext, encoding='raw_unique_escape')

  3. 在WPF中使用WriteableBitmap对接工业相机及常用操作

    写作背景 写这篇文章主要是因为工业相机(海康.大恒等)提供的.NET开发文档和示例程序都是用WinForm项目来说明举例的,而在WPF项目中对图像的使用和处理与在WinForm项目中有很大不同.在Wi ...

  4. 《刚刚问世》系列初窥篇-Java+Playwright自动化测试-1-环境准备与搭建

    1.简介 Python+Playwright系列的文章还没有结束,就有好的小伙伴或者童鞋们私信公众号留言,问宏哥什么时候出Java语言的Playwright的自动化测试文章.本来想趁热打铁将Pytho ...

  5. power bi 如何删除敏感度标签

    经验证,此方法不够彻底,我的office excel打开后还是要添加敏感度标签,即使我把敏感度标签删掉也不行. 当我把创建敏感度标签的管理员账户删掉之后,虽然打开excel还是会显示敏感度标签,但是已 ...

  6. 【Python】基于动态规划和K聚类的彩色图片压缩算法

    引言 当想要压缩一张彩色图像时,彩色图像通常由数百万个颜色值组成,每个颜色值都由红.绿.蓝三个分量组成.因此,如果我们直接对图像的每个像素进行编码,会导致非常大的数据量.为了减少数据量,我们可以尝试减 ...

  7. ubuntu16.04 python2&3 pip升级后报错:sys.stderr.write(f"ERROR: {exc}")

    ubuntu16.04 python2&3 pip升级后报错: sys.stderr.write(f"ERROR: {exc}") 描述 最近使用ubuntu16.04上的 ...

  8. VUE手稿1

  9. SpringCloud 微服务简介

    一.认识微服务随着互联网行业的发展,对服务的要求也越来越高,服务架构也从单体架构逐渐演变为现在流行的微服务架构.这些架构之间有怎样的差别呢? 1.单体架构:将业务的所有功能集中在一个项目中开发,打成一 ...

  10. 解决方案 | 在 Tkinter 中导入 pywinauto/pyautogui 时窗口大小发生变化

    上面问题也可以换一个说法,pywinauto/pyautogui 时改变了tkinter的原有的窗口大小.这个问题困扰了我好几天而且网上有这样的问题但是并没有答案,今天摸索出答案给大家分享下.解决方法 ...