4-偏微分方程及其求解实例材料.ppt

当B2-4AC0时，方程为椭圆型（elliptic）PDE 当B2-4AC=0时，方程为抛物型（parabolic）PDE 当B2-4AC0时，方程为双曲型（hyperbolic）PDE;(1) 导热方程:;偏微分方程的边界条件;偏微分方程数值求解方法;偏微分方程数值求解方法;偏微分方程数值求解方法;偏微分方程数值求解方法;有限差分法求解偏微分方程;有限差分法求解偏微分方程;有限差分法求解偏微分方程; 1、显式差分法 2、隐式差分法 ;j;function PDE1_FDM_im % 显式差分法求解一维热传导方程 clear all;clc n=6;% 空间节点数 m=8;% 时间节点数 T(1,1:m)=100;T(n,1:m)=0; % 边界条件 T(1:n,1)=0; % 初始条件 alpha=2; % 热扩散系数,m/s a=10; % for x坐标长度,cm b=8; % for 时间长度,s h=a/(n-1); % 空间步长 k=b/(m-1); % 时间步长 r=alpha.*k./h.^2;;for j=2:m % for time for i=2:n-1 % for x T(i,j+1)=(1-2.*r).*T(i,j)+r.*(T(i+1,j)+T(i-1,j)); end end T8=T(1:n,m);j;i=1;偏微分方程的求解实例1: 恒靠近速度时两等直径液滴形成的液膜内流  体排液速率的模拟—问题描述(h(r,t)~t ?);偏微分方程的求解实例1: 恒靠近速度时两等直径液滴形成的液膜内流  体排液速率的模拟—模型建立;偏微分方程的求解实例1: 恒靠近速度时两等直径液滴形成的液膜内流  体排液速率的模拟—方程离散;r;function CVFD clear all;clc; format long e h0c=2.8e-6; % 初始液膜中心厚度, m mu=0.3; % 液相流体的粘度, Pa.s rb=0.0015; % 液滴半径,m theta=0.03; % 液膜的表面张力, N/m ra=0.5*rb; % 膜出口压力为零区域 V=2e-6; % 膜出口ra处靠近速度m/s m=100; r=[0,0,linspace(0,ra,m)]; n=m+2; dr=ra./(m-1); % 节点的步长,Δr,m h0=h0c+r.^2./rb; % 初始膜厚度m h0(1)=h0(5);h0(2)=h0(4);;hold on k=10; tf=12; for i=1:k [t,h]=ode15s(@Non,[0,tf],h0,[],mu,theta,rb,V,dr,n) k=length(t); % 不同时刻的液膜厚度 plot(r(3:n)./ra,h(k,3:n)./h0c) h0=h(k,:); end hold off function dhdt=Non(t,h,mu,theta,rb,V,dr,n) h(n+1)=(4.*dr.^2./rb+2.*h(n)-h(n-1).*(1-1./(2.*(n-3))))./(1./(2.*(n-3))+1); dhdt(3)=-theta./9./mu.*(h(3).^3.*(2.*h(5)-8.*h(4)+6.*h(3))./dr.^4); h(1)=h(5);h(2)=h(4);;for i=4:n-1 Dhr(i)=(h(i+1)-h(i-1))/(2*dr); D2hr(i)=(h(i+1)-2*h(i)+h(i-1))/(dr^2); D3hr(i)=(h(i+2)-2.*h(i+1)+2.*h(i-1)-… h(i-2))./(2.*dr.^3); D4hr(i)=(h(i+2)-4.*h(i+1)+6.*h(i)-… 4.*h(i-1)+h(i-2))./(dr.^4); H1(i)=D4hr(i)+2./((i-1).*dr).*D3hr(i)-… 1./((i-3).^2.*dr.^2).*D2hr(i)+… 1./((i-3).^3.*dr.^3).*Dhr(i); H2(i)=Dhr(i).*(D3hr(i)+1./((i-3).*dr).*… D2hr(i)-1./((i-3).^2.*dr.^2).*Dhr(i)); dhdt(i)=-theta./8./mu.*(h(i).^3.*H1(i)./3+…