24小时热门版块排行榜

返回列表

晓嘉嘉

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 54.5
散金: 100
红花: 1
帖子: 63
在线: 15.4小时
虫号: 6058623
注册: 2017-03-19
性别: GG
专业: 机械工程

[求助] 气泡动力学

根据别人论文的方程和数据提供对单气泡的动力学做了相应的仿真，但是老是出现警告: 在 t=6.355348e-06 处失败。在时间 t 处，若不将步长降至允许的最小值(1.355253e-20)以下，积分公差要求无法满足。希望老师们、大神们给予菜鸟的我一些帮助和指导。编码如下：
function jjdcxformat('long');
global f R0 p0 p_a Ki pi rho sigma mu w p_v m p_g I c n
f=0.2e6;%输入频率
R0=6e-6;%初始半径
time=1/f;%时间
m=20;%气体摩尔数
n=8.314;%气体常数
p0=1.013e5;%作用在气泡壁上的流体静力，取大气压
p_g=1.01e5;%水温下的蒸汽压
p_a=1.5*p0;%气泡内的原始压力
Ki=1;%气体的多变指数
rho=1000;%水的密度
p_v=0.0233e5;%气泡内的蒸汽压
sigma=7.25e-2;%液体的表面张力
mu=1.0e-3;%液体的运动粘度
c=1500;%声速
pi=3.1415926;
w=2*pi*f;
ttime=5*time;
options = odeset('RelTol',1e-10,'AbsTol',1e-8);
[t,y]=ode45('qipao3',[0 ttime],[R0,0],options);
plot(t/time,y(:,1)/R0,'b');
hold on;function dy=qipao3(t,y)
global f R0 Ki p0 p_a pi rho sigma mu h p_v
h=R0/8.54;
f11=p_a*sin(2*pi*f*t);
dy=[y(2);((p0+2*sigma/R0-p_v)*(R0/y(1))^(3*Ki)-2*sigma/y(1)-4*mu*y(2)/y(1)+p_v-p0-f11)/(rho*y(1))-(3*y(2).^2)/(2*y(1))];
end

回复此楼