| 查看: 3143 | 回复: 3 | ||
[求助]
气泡动力学
|
|
根据别人论文的方程和数据提供对单气泡的动力学做了相应的仿真，但是老是出现警告: 在 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 |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有10人回复
-
什么时候能放榜呀？
已经有3人回复
-
一个有机合成实验室都需要哪些设备？
已经有8人回复
-
UV压敏胶开发
已经有5人回复
-
面上提前没消息，有中的吗
已经有16人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有25人回复
-
2026年国自然面上资助率
已经有22人回复
-
基础研究怎么拉横向，学校到款任务越来越多，难以完成 拉横向，都有哪些途径啊
已经有12人回复
-
微信指数没变化，科研之友没阅读
已经有19人回复
2楼2018-01-14 12:32:26
3楼2018-01-17 12:30:27
4楼2019-10-03 21:47:16