% 一维对流方程迎风格式、Lax格式、FTCS格式差分法计算% 潭花林清华大学航天航空学院% FTCS格式对于一维对流方程不稳定，最好不用clcclear all% 1.参数定义dx=1;x1=-18;x2=18;x=x1:dx:x2;L1=length(x);% dt=0.5*dx; % 收敛dt=2*dx; % 不收敛t1=0;t2=t1+80*dt;t=t1:dt:t2;L2=length(t);alpha=1;lambda=alpha*dt/dx;geshi=1; % 迎风格式% geshi=2; % Lax格式% geshi=3; % FTCS格式% 2.显式求解zeta=zeros(L1,L2);for kk=1:3geshi=kk;for ii=1:L1if x(ii)>0zeta(ii,1)=1;elseif x(ii)==0zeta(ii,1)=1/2;elseif x(ii)<0zeta(ii,1)=0;endendendendif geshi==1for ii=2:L1for jj=1:(L2-1)zeta(ii,jj+1)=zeta(ii,jj)-lambda*(zeta(ii,jj)-z eta(ii-1,jj));endzeta(1,jj+1)=zeta(2,jj+1);endzeta1=zeta;else if geshi==2for ii=2:(L1-1)for jj=1:(L2-1)zeta(ii,jj+1)=(zeta(ii+1,jj)+zeta(ii-1,jj))/2-. ..lambda/2*(zeta(ii+1,jj)-zeta(ii-1,jj));endzeta(1,jj+1)=zeta(2,jj+1);zeta(L1,jj+1)=zeta(L1,jj)-lambda*(zeta(L1,jj)-z eta(L1-1,jj));endzeta2=zeta;else if geshi==3for ii=2:(L1-1)for jj=1:(L2-1)zeta(ii,jj+1)=zeta(ii,jj)-lambda/2*(zeta(ii+1,j j)-zeta(ii-1,jj));endzeta(1,jj+1)=zeta(2,jj+1);zeta(L1,jj+1)=zeta(L1,jj)-lambda*(zeta(L1,jj)-z eta(L1-1,jj));endzeta3=zeta;endendendend% 3.绘图% 3.1 t=0figure(1)n=1;plot(x,zeta1(1:L1,n),'-k',x,zeta2(1:L1,n),'-.k' ,x,zeta3(1:L1,n),'--k')% 作图% axis equal %%% 是否要求x、y坐标间距相等% grid on %%% 是否要求画网格xlabel('x/m'),ylabel('t/s') % %% x,y轴表示的变量含义%text(1,2,'f(x)') %%% 图中文字标识legend('迎风格式','Lax格式','FTCS格式') %%% 不同曲线的线型区分title('t=0时刻的计算结果') %%% 标题axis([-18,18,-0.2,1.2])% 3.2 t=10figure(2)n=(10-t(1))/dt;plot(x,zeta1(1:L1,n),'-k',x,zeta2(1:L1,n),'-.k' ,x,zeta3(1:L1,n),'--k')% 作图% axis equal %%% 是否要求x、y坐标间距相等% grid on %%% 是否要求画网格xlabel('x/m'),ylabel('t/s') % %% x,y轴表示的变量含义%text(1,2,'f(x)') %%% 图中文字标识legend('迎风格式','Lax格式','FTCS格式') %%% 不同曲线的线型区分title('t=10s时刻的计算结果') %%% 标题% 3.3 t=20figure(3)n=(20-t(1))/dt;plot(x,zeta1(1:L1,n),'-k',x,zeta2(1:L1,n),'-.k' ,x,zeta3(1:L1,n),'--k')% 作图% axis equal %%% 是否要求x、y坐标间距相等% grid on %%% 是否要求画网格xlabel('x/m'),ylabel('t/s') %%% x,y轴表示的变量含义%text(1,2,'f(x)') %%% 图中文字标识legend('迎风格式','Lax格式','FTCS格式') %%% 不同曲线的线型区分title('t=20s时刻的计算结果') %%% 标题% 3.4 t=40figure(4)n=(40-t(1))/dt;plot(x,zeta1(1:L1,n),'-k',x,zeta2(1:L1,n),'-.k' ,x,zeta3(1:L1,n),'--k')% 作图% axis equal %%% 是否要求x、y坐标间距相等% grid on %%% 是否要求画网格xlabel('x/m'),ylabel('t/s') % %% x,y轴表示的变量含义%text(1,2,'f(x)') %%% 图中文字标识legend('迎风格式','Lax格式','FTCS格式') %%% 不同曲线的线型区分title('t=40s时刻的计算结果') %%% 标题聚乙烯（PE）简介1.1聚乙烯化学名称：聚乙烯英文名称：polyethylene，简称PE结构式：聚乙烯是乙烯经聚合制得的一种热塑性树脂，也包括乙烯与少量α-烯烃的共聚物。
聚乙烯是五大合成树脂之一，是我国合成树脂中产能最大、进口量最多的品种。
1.1.1聚乙烯的性能1.一般性能聚乙烯为白色蜡状半透明材料，柔而韧，比水轻，无嗅、无味、无毒，常温下不溶于一般溶剂，吸水性小，但由于其为线性分子可缓慢溶于某些有机溶剂，且不发生溶胀。
工业上为使用和贮存的方便通常在聚合后加入适量的塑料助剂进行造粒，制成半透明的颗粒状物料。
PE易燃，燃烧时有蜡味，并伴有熔融滴落现象。
聚乙烯的性质因品种而异，主要取决于分子结构和密度，也与聚合工艺及后期造粒过程中加入的塑料助剂有关。
2.力学性能PE是典型的软而韧的聚合物。
除冲击强度较高外，其他力学性能绝对值在塑料材料中都是较低的。
PE密度增大，除韧性以外的力学性能都有所提高。
LDPE 由于支化度大，结晶度低，密度小，各项力学性能较低，但韧性良好，耐冲击。
HDPE支化度小，结晶度高，密度大，拉伸强度、刚度和硬度较高，韧性较差些。
相对分子质量增大，分子链间作用力相应增大，所有力学性能，包括韧性也都提高。
几种PE的力学性能见表1-1。
表1-1 几种PE力学性能数据3.热性能PE受热后，随温度的升高，结晶部分逐渐熔化，无定形部分逐渐增多。
其熔点与结晶度和结晶形态有关。
HDPE的熔点约为125～137℃，MDPE的熔点约为126～134℃，LDPE的熔点约为105～115℃。
相对分子质量对PE的熔融温度基本上无影响。
PE的玻璃化温度（T g）随相对分子质量、结晶度和支化程度的不同而异，而且因测试方法不同有较大差别，一般在-50℃以下。
PE在一般环境下韧性良好，耐低温性(耐寒性)优良，PE的脆化温度(T b)约为-80～-50℃，随相对分子质量增大脆化温度降低，如超高相对分子质量聚乙烯的脆化温度低于-140℃。
PE的热变形温度(T HD)较低，不同PE的热变形温度也有差别，LDPE约为38～50℃(0.45MPa，下同)，MDPE约为50～75℃，HDPE约为60～80℃。
PE的最高连续使用温度不算太低，LDPE约为82～100℃，MDPE约为105～121℃，HDPE为121℃，均高于PS和PVC。
PE的热稳定性较好，在惰性气氛中，其热分解温度超过300℃。
PE的比热容和热导率较大，不宜作为绝热材料选用。
PE的线胀系数约在(15～30)×10-5K-1之间，其制品尺寸随温度改变变化较大。
几种PE的热性能见表1-2。
表1-2几种PE热性能4.电性能PE分子结构中没有极性基团，因此具有优异的电性能，几种PE的电性能见表1-3。
PE的体积电阻率较高，介电常数和介电损耗因数较小，几乎不受频率的影响，因而适宜于制备高频绝缘材料。
它的吸湿性很小，小于0.01％（质量分数），电性能不受环境湿度的影响。
尽管PE具有优良的介电性能和绝缘性，但由于耐热性不够高，作为绝缘材料使用，只能达到Y级（工作温度≤90℃）。
表1-3聚乙烯的电性能5.化学稳定性PE是非极性结晶聚合物，具有优良的化学稳定性。
室温下它能耐酸、碱和盐类的水溶液，如盐酸、氢氟酸、磷酸、甲酸、醋酸、氨、氢氧化钠、氢氧化钾以及各类盐溶液(包括具有氧化性的高锰酸钾溶液和重铬酸盐溶液等)，即使在较高的浓度下对PE也无显著作用。
但浓硫酸和浓硝酸及其他氧化剂对聚乙烯有缓慢侵蚀作用。
PE在室温下不溶于任何溶剂，但溶度参数相近的溶剂可使其溶胀。
随着温度的升高，PE结晶逐渐被破坏，大分子与溶剂的作用增强，当达到一定温度后PE可溶于脂肪烃、芳香烃、卤代烃等。
如LDPE能溶于60℃的苯中，HDPE能溶于80～90℃的苯中，超过100℃后二者均可溶于甲苯、三氯乙烯、四氢萘、十氢萘、石油醚、矿物油和石蜡中。
但即使在较高温度下PE仍不溶于水、脂肪族醇、丙酮、乙醚、甘油和植物油中。
PE在大气、阳光和氧的作用下易发生老化，具体表现为伸长率和耐寒性降低，力学性能和电性能下降，并逐渐变脆、产生裂纹，最终丧失使用性能。
为了防止PE的氧化降解，便于贮存、加工和应用，一般使用的PE原料在合成过程中已加入了稳定剂，可满足一般的加工和使用要求。
如需进一步提高耐老化性能，可在PE中添加抗氧剂和光稳定剂等。
6.卫生性PE分子链主要由碳、氢构成，本身毒性极低，但为了改善PE性能，在聚合、成型加工和使用中往往需添加抗氧剂和光稳定剂等塑料助剂，可能影响到它的卫生性。