对我来说,我认为这将是Rayleigh-Plesset的任何变化方程,它描述了体积下的泡沫液体振荡压力场的影响。
下面,看到的一个例子Keller-Miksis方程(类似于Rayleigh-Plesset,但占温和的液体压缩性)模拟的径向振动惯性collpasing泡沫。
清楚,clc
T0 = 298;%环境温度,K
频率= 26500;%驱动频率,赫兹
r0 = 4.5;%初始气泡半径,微米
p0 = 1;%环境压力,自动取款机
pA = 1.2;%声压,atm
rho0 = 997;%的液体密度、公斤m3
sigmaL = 0.072;%表面张力,N m - 1
muL = 1 e - 3;%的液体粘度、公斤m - 1 s - 1
c = 1485;%声速,m s - 1
γ= 1.4;%具体加热空气比率,-
tspan =[0 1 /频率);
IC = [r0 * 1 e-6 0 p0 * 101325 + 2 * sigmaL / (r0 * 1 e-6) T0);
odefun = @ (t, w) keller_miksis (t, w, r0 * 1 e-6 T0, rho0,频率,p0 * 101325 * 101325, sigmaL, muL, c,γ);
(时间,索尔)= ode15s (odefun tspan, IC);
t = *频率;
R =索尔(:1)* 1 e + 6;
S =索尔(:,2);
p =索尔(:,3)/ 101325;
T =索尔(:,4);
集(0,“defaulttextinterpreter”,“乳胶”)
集(0,“defaultAxesTickLabelInterpreter”,“乳胶”)
集(0,“defaultLegendInterpreter”,“乳胶”)
图(1)
次要情节(2、2、1)
情节(t, R)
标题(气泡半径的)
包含(“\ωt美元”)
ylabel (' R美元,美元\μm美元的)
次要情节(2 2 2)
semilogy (t, p)
标题(“气压”)
包含(“\ωt美元”)
ylabel (“$ p $, atm机”)
次要情节(2,2,3)
情节(t, t (: 1))
标题(气体温度的)
包含(“\ωt美元”)
ylabel (‘T美元,K’)
次要情节(2,2,4)
情节(t、abs (S / c))
标题(“界面马赫数”)
包含(“\ωt美元”)
ylabel (“马美元,”)
函数f = keller_miksis (t w r0、T0 rho0,频率,p0, pA, sigmaL, muL, c,γ)
R = w (1);
博士= w (2);
p = w (3);
T = w (4);
dpdt = 3 *γ* p * / R博士;
dTdt = 3 * (1-gamma) * T * / R博士;
pinf = p0-pA * sin(2 *π*频率* t);
pL = (p0 + 2 * sigmaL / r0) * (r0 / R) ^(3 *γ)2 * sigmaL / R-4 * muL * / R博士;
ddR =((1 +博士/ c) * (pL-pinf) / (rho0) -1.5 * (1-dR / (3 * c)) *(^ 2)博士+…
(R / c (rho0 *)) * (dpdt + (2 * sigmaL * + 4 * muL * ^ 2)博士/ (R ^ 2))) / ((1-dR / c) * R + 4 * muL / (rho0 * c));
f =[博士;…
ddR;…
dpdt;…
dTdt];
结束
