[未答复] 速度势函数求解,求解出来的结果不正确,请问如何处理?

[复制链接]
zhangtangjia 发表于 2022-9-20 21:15:33
本帖最后由 zhangtangjia 于 2022-9-20 21:24 编辑

我复现论文结果,发现结果与论文结果不一致,怀疑是我求解有点问题,但是实在不知道怎么处理?
代码

clc;clear all;close all
%% 翼型绘制
syms z ksi c0 R theta W_uni ksi0 U alpha

z=ksi+c0^2./ksi;
T_disk=-4*pi*R*U*sin(alpha);
W_uni=U*(ksi-ksi0)*exp(-i*alpha)+(R^2/(ksi-ksi0)*U*exp(i*alpha)-i*T_disk/(2*pi)*log(ksi-ksi0));%速度势
f_conjugat_ksi=diff(W_uni,ksi,1);
f_conjugat_z=diff(W_uni,ksi,1)*diff(z,ksi)^(-1);



% ksi=1/2*(z+-sqrt(z^2-4*c0^2));
ksi0=-0.35;U=0.2;alpha=20;
R=4.35:0.01:10;

ksi=R*exp(i*theta)+ksi0;
c0=R-abs(ksi0);
theta=0:0.1:2*pi;
c0=[];
ksi=[];
z=[];
f_conjugat_ksi_res=[];
f_conjugat_z_res=[];
W_uni_res=[];
theta=0:0.01:2*pi;
% 带入数值求解
for n=1:length(R)
   c0(n)=R(n)-abs(ksi0);
    ksi(n,:)=R(n).*exp(i*theta)+ksi0;
    z(n,:)=ksi(n,:)+c0(n).^2./ksi(n,:);
%     sub [a b],[A B]
    W_uni_res(n,:)=(U*exp(alpha*1i)*R(n).^2)./(ksi(n,:) - ksi0) + U*log(ksi(n,:) - ksi0)*sin(alpha)*R(n).*pi^2*2i + U*exp(-alpha*1i)*(ksi(n,:) - ksi0);
    f_conjugat_ksi_res(n,:)=U*exp(-alpha*1i)-(R(n).^2*U*exp(alpha*1i))./(ksi(n,:) - ksi0).^2 + (R(n).*U*pi^2*sin(alpha)*2i)./(ksi(n,:)-ksi0);
    f_conjugat_z_res(n,:)=-(U*exp(-alpha*1i)-(R(n).^2*U*exp(alpha*1i))./(ksi(n,:) - ksi0).^2 + (R(n).*U*pi^2*sin(alpha)*2i)./(ksi(n,:)-ksi0))./(c0(n)^2./ksi(n,:).^2 - 1);
end


f_conjugat_r=real(z);
f_conjugat_i=imag(z);
abs_f_conjugat=abs(f_conjugat_ksi_res);
figure(2)
pcolor(f_conjugat_r,f_conjugat_i,real(f_conjugat_ksi_res))
shading interp
colorbar
hold on
box on

速度

速度

速度势函数

速度势函数

自己计算结果

自己计算结果

论文结果

论文结果
您需要登录后才可以回帖 登录 | 注册

本版积分规则

相关帖子
相关文章
热门教程
站长推荐
快速回复 返回顶部 返回列表