| 查看: 2313 | 回复: 1 | |||
sas4444木虫 (小有名气)
|
[交流]
rayleigh-plesset仿真 已有1人参与
|
|
function Rx=RP(t,R) P0=100000; %初始压强P0=0.1MPa T0=298.16; %初始温度T0=273.16+25 Hv=2257.2; %焓Hv=2257.2KJ/Kg R0=0.01; %初始半径R0=0.02m rho=1000; %液体密度ρ=1000Kg/cm2 nu=1*10^-6; %运动粘度ν=1.0*10^-6m2/s gama=0.0725; %表面张力γ=0.0725N/m m=0.018; %摩尔质量m=0.018Kg/mol Rg=8.3144; %通用气体常数Rg=8.3144KJ/(Kg·mol·K) alpha=R0/8.86; %常数α=R0/8.86 r=0.008; %距离r=8*10^-3m k=1.66; %绝热系数κ=1.66 Rx=zeros(2,1); Rx(1)=R(2); Rx(2)=(P0*exp((Hv*m/Rg)*(1/T0-(R(1)^3-alpha^3)^(k-1)/(T0*R0^(3*(k-1)))) ... +2*gama/R0*(R0^(3*(k-1))/(R(1)^3-alpha^3)^(k-1))*(R0/R(1))^(3*k)-(P0-rho/2*(R(1)^2/r^2*R(2))^2))/rho ... -(1.5*R(2)^2+4*nu*R(2)/R(1)+2*gama/(rho*R(1))))/R(1); end 然后 [t,R]=ode45('RP',[0,100],[0.000007,2000]); 如何避免的出来的是复数或者是无穷大 |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有6人回复
-
明天应该可查了！？
已经有5人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有7人回复
-
没有任何消息-是不是就凉了
已经有8人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
人气不行了
已经有10人回复
-
我面上完蛋了
已经有10人回复
-
能否退出参与的面上项目解除限项
已经有24人回复
-
filecode，4个jtjc了
已经有17人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 多时滞线性系统仿真
已经有6人回复
- adams2012 最新出版的书
已经有180人回复
- 流体力学公式中存在的问题
已经有5人回复
- 发些CFD基础理论的书
已经有83人回复
- ANSYS二维电场模拟电阻变化，求指教
已经有4人回复
- Rayleigh信道与AWGN信道的BER
已经有6人回复
- matlab 32位和64位 区别
已经有7人回复
- 怎么把图中的数据转成表格导出来
已经有8人回复
- 【原创】两种热虹吸塔在HYSYS中的动态实现及分析
已经有7人回复
2楼2018-06-28 09:58:20
