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. Python str 转 b’二进制串

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

  2. 高通mm-camera平台 Camera移植

    高通Cam-X平台 Camera移植 注:此文档以在高通8916平台移植OV5648为例,给大家讲解Android SOC的底层Camera. Reference: https://blog.csdn ...

  3. 嵌入式HLS 案例开发步骤分享——基于Zynq-7010/20工业开发板(1)

    目 录 前 言 3 1 HLS 开发流程说明 5 1.1 HLS 工程导入 5 1.2 编译与仿真 6 1.3 综合 8 1.4 IP 核封装 10 1.5 IP 核测试 14 前 言 本文主要介绍 ...

  4. 转-tomcat在控制台打印时乱码

    描述 使用windows的控制台启动tomcat (8.5),Tomcat终端打印出现乱码,在CSDN上找到一个简单有效的解决办法,在此记录. 解决方案 找到Tomcat的根目录下的 /conf/lo ...

  5. spring-关于组件的注入及获取流程

    一.组件注入的基本流程: 容器初始化: Spring应用启动时,会读取配置(如XML配置.注解配置等),并根据这些配置创建Bean定义(BeanDefinition). 根据Bean定义,Spring ...

  6. yb课堂之用户登陆校验拦截器开发 《十一》

    开发对应的登陆拦截器 开发loginInterceptor 登陆校验成功放行 登陆不成功返回json数据 LoginInterceptor.java package net.ybclass.onlin ...

  7. Mybatis 一级缓存

    Mybatis一级缓存介绍 什么是缓存 程序经常要调用的对象存在内容中,方法其使用时可以快速调用,不必去数据库或者其他持久化设备中查询,主要就是提高性能 Mybatis一级缓存 简介:一级缓存的作用域 ...

  8. oeasy教您玩转vim - 43 - # 替换模式

    ​ 替换模式 回忆上节课内容 上次我们学到的最有用的就是c 他和d的区别就是删除之后进入到插入模式 c可以配合motion 可以用ciw来快速删除当前光标所在的单词 可i和a 配合的文字块 w wor ...

  9. oeasy教您玩转vim - 77 - # 保留环境viminfo

    ​ 保留环境viminfo 回忆组合键映射的细节 上次我们定义了session :mks 还可以加载会话session :source Session.vim vim -S Session.vim 基 ...

  10. TIER 2: Archetype

    TIER 2: Archetype 扫描 nmap 使用 nmap 进行扫描目标 IP,发现目标是 Windows 服务器,开放 SMB 和 SQL Server 服务. SMB SMB 之前已经接触 ...