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. 10-Python进程与线程

    Python进程 创建新进程 from multiprocessing import Process import time def run_proc(name): #子进程要执行的代码 for i ...

  2. P8594 「KDOI-02」一个仇的复

    我会组合数! 首先发现同一列只有被不同的横块填或被一个相同的竖块填,且用竖块填完1列之后,会分成两个封闭的长方形,而长方形内部则用横块来填充. 先考虑一个子问题,某个 \(2 \times n\) 长 ...

  3. 解决:编译安卓源码时 JDK 报错 error='Not enough space' (errno=12)

    背景 在编译 Android 10 代码的时候,OpenJDK发现报错: OpenJDK 64-Bit Server VM warning: INFO: os::commit_memory(.., . ...

  4. qt中的 connect 函数

    1.connect()函数实现的是信号与槽的关联. 注意:只有QO bject类及其派生的类才能使用信号和槽的机制 2.函数原型 static QMetaObject::Connection conn ...

  5. NXP i.MX 8M Mini开发板规格书(四核ARM Cortex-A53 + 单核ARM Cortex-M4,主频1.6GHz)

    1 评估板简介 创龙科技TLIMX8-EVM是一款基于NXP i.MX 8M Mini的四核ARM Cortex-A53 + 单核ARM Cortex-M4异构多核处理器设计的高性能评估板,由核心板和 ...

  6. VS License Header Manager 插件

    /************************************************************************************* * * 文 件 名: %F ...

  7. Asp .Net Core 系列:基于 Castle DynamicProxy + Autofac 实践 AOP 以及实现事务、用户填充功能

    目录 什么是 AOP ? .Net Core 中 有哪些 AOP 框架? 基于 Castle DynamicProxy 实现 AOP IOC中使用 Castle DynamicProxy 实现事务管理 ...

  8. Springboot+Shiro+Mybatis+mysql实现权限安全认证

    Shiro是Apache 的一个强大且易用的Java安全框架,执行身份验证.授权.密码学和会话管理.Shiro 主要分为两个部分就是认证和授权两部分 一.介绍 Subject代表了当前用户的安全操作 ...

  9. Oracle plsql中文字段乱码,where条件中文字段搜不到结果集

    设置系统环境变量 变量名:NLS_LANG 变量值:AMERICAN_AMERICA.ZHS16GBK

  10. CF372C

    思路 根据题意可以得到dp转移方程是 \(f_{i,j}=\max\{f_{i-1,k}+b_i-|a_i-j|\}\) 而且 \(j-(t_{i}-t_{i-1})\times d\le k\le ...