计算程序计算流体力学对流方程有限差分法Lax格式迎风格式FTCS格式.docx
- 文档编号:10195466
- 上传时间:2023-02-09
- 格式:DOCX
- 页数:14
- 大小:34.66KB
计算程序计算流体力学对流方程有限差分法Lax格式迎风格式FTCS格式.docx
《计算程序计算流体力学对流方程有限差分法Lax格式迎风格式FTCS格式.docx》由会员分享,可在线阅读,更多相关《计算程序计算流体力学对流方程有限差分法Lax格式迎风格式FTCS格式.docx(14页珍藏版)》请在冰豆网上搜索。
计算程序计算流体力学对流方程有限差分法Lax格式迎风格式FTCS格式
%一维对流方程迎风格式、Lax格式、FTCS格式差分法计算
%潭花林清华大学航天航空学院
%FTCS格式对于一维对流方程不稳定,最好不用
clc
clearall
%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);
forkk=1:
3
geshi=kk;
forii=1:
L1
ifx(ii)>0
zeta(ii,1)=1;
else
ifx(ii)==0
zeta(ii,1)=1/2;
else
ifx(ii)<0
zeta(ii,1)=0;
end
end
end
end
ifgeshi==1
forii=2:
L1
forjj=1:
(L2-1)
zeta(ii,jj+1)=zeta(ii,jj)-lambda*(zeta(ii,jj)-zeta(ii-1,jj));
end
zeta(1,jj+1)=zeta(2,jj+1);
end
zeta1=zeta;
elseifgeshi==2
forii=2:
(L1-1)
forjj=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));
end
zeta(1,jj+1)=zeta(2,jj+1);
zeta(L1,jj+1)=zeta(L1,jj)-lambda*(zeta(L1,jj)-zeta(L1-1,jj));
end
zeta2=zeta;
elseifgeshi==3
forii=2:
(L1-1)
forjj=1:
(L2-1)
zeta(ii,jj+1)=zeta(ii,jj)-lambda/2*(zeta(ii+1,jj)-zeta(ii-1,jj));
end
zeta(1,jj+1)=zeta(2,jj+1);
zeta(L1,jj+1)=zeta(L1,jj)-lambda*(zeta(L1,jj)-zeta(L1-1,jj));
end
zeta3=zeta;
end
end
end
end
%3.绘图
%3.1t=0
figure
(1)
n=1;
plot(x,zeta1(1:
L1,n),'-k',x,zeta2(1:
L1,n),'-.k',x,zeta3(1:
L1,n),'--k')%作图
%axisequal%%%是否要求x、y坐标间距相等
%gridon%%%是否要求画网格
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.2t=10
figure
(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')%作图
%axisequal%%%是否要求x、y坐标间距相等
%gridon%%%是否要求画网格
xlabel('x/m'),ylabel('t/s')%%%x,y轴表示的变量含义
%text(1,2,'f(x)')%%%图中文字标识
legend('迎风格式','Lax格式','FTCS格式')%%%不同曲线的线型区分
title('t=10s时刻的计算结果')%%%标题
%3.3t=20
figure(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')%作图
%axisequal%%%是否要求x、y坐标间距相等
%gridon%%%是否要求画网格
xlabel('x/m'),ylabel('t/s')%%%x,y轴表示的变量含义
%text(1,2,'f(x)')%%%图中文字标识
legend('迎风格式','Lax格式','FTCS格式')%%%不同曲线的线型区分
title('t=20s时刻的计算结果')%%%标题
%3.4t=40
figure(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')%作图
%axisequal%%%是否要求x、y坐标间距相等
%gridon%%%是否要求画网格
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力学性能数据
性能
LDPE
LLDPE
HDPE
超高相对分子质量聚乙烯
邵氏硬度(D)
拉伸强度/MPa
拉伸弹性模量/MPa
压缩强度/MPa
缺口冲击强度/kJ·m-2
弯曲强度/MPa
41~46
7~20
100~300
12.5
80~90
12~17
40~50
15~25
250~550
—
>70
15~25
60~70
21~37
400~1300
22.5
40~70
25~40
64~67
30~50
150~800
—
>100
—
3.热性能
PE受热后,随温度的升高,结晶部分逐渐熔化,无定形部分逐渐增多。
其熔点与结晶度和结晶形态有关。
HDPE的熔点约为125~137℃,MDPE的熔点约为126~134℃,LDPE的熔点约为105~115℃。
相对分子质量对PE的熔融温度基本上无影响。
PE的玻璃化温度(Tg)随相对分子质量、结晶度和支化程度的不同而异,而且因测试方法不同有较大差别,一般在-50℃以下。
PE在一般环境下韧性良好,耐低温性(耐寒性)优良,PE的脆化温度(Tb)约为-80~-50℃,随相对分子质量增大脆化温度降低,如超高相对分子质量聚乙烯的脆化温度低于-140℃。
PE的热变形温度(THD)较低,不同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热性能
性能
LDPE
LLDPE
HDPE
超高相对分子质量聚乙烯
熔点/℃
热降解温度(氮气)/℃
热变形温度(0.45MPa)/℃
脆化温度/℃
线性膨胀系数/(×10-5K-1)
比热容/J·(kg·K)-1
热导率/W·(m·K)-1
105~115
>300
38~50
-80~-50
16~24
2218~2301
0.35
120~125
>300
50~75
-100~-75
—
—
—
125~137
>300
60~80
-100~-70
11~16
1925~2301
0.42
190~210
>300
75~85
-140~-70
—
—
—
4.电性能
PE分子结构中没有极性基团,因此具有优异的电性能,几种PE的电性能见表1-3。
PE的体积电阻率较高,介电常数和介电损耗因数较小,几乎不受频率的影响,因而适宜于制备高频绝缘材料。
它的吸湿性很小,小于0.01%(质量分数),电性能不受环境湿度的影响。
尽管PE具有优良的介电性能和绝缘性,但由于耐热性不够高,作为绝缘材料使用,只能达到Y级(工作温度≤90℃)。
表1-3聚乙烯的电性能
性能
LDPE
LLDPE
HDPE
超高相对分子质量聚乙烯
体积电阻率/Ω·cm
介电常数/F·m-1(106Hz)
介电损耗因数(106Hz)
介电强度/kV·mm-1
≥1016
2.25~2.35
<0.0005
>20
≥1016
2.20~2.30
<0.0005
45~70
≥1016
2.30~2.35
<0.0005
18~28
≥1017
≤2.35
<0.0005
>35
5.化学稳定性
PE是非极性结晶聚合物,具有优良的化学稳定性。
室温下它能耐酸、碱和盐类的水溶液,如盐酸、氢氟酸、磷酸、甲酸、醋酸、氨、氢氧化钠、氢氧化钾以及各类盐溶液(包括具有氧化性的高锰酸钾溶液和重铬酸盐溶液等),即使在较高的浓度下对PE也无显著作用。
但浓硫酸和浓硝酸及其他氧化剂对聚乙烯有缓慢侵蚀作用。
PE在室温下不溶于任何溶剂,但溶度参数相近的溶剂可使其溶胀。
随着温度的升高,PE结晶逐渐被破坏,大分子与溶剂的作用增强,当达到一定温度后PE可溶于脂肪烃、芳香烃、卤代烃等。
如LDPE能溶于60℃的苯中,HDPE能溶于80~90℃的苯中,超过100℃后二者均可溶于甲苯、三氯乙烯、四氢萘、十氢萘、石油醚、矿物油和石蜡中。
但即使在较高温度下PE仍不溶于水、脂肪族醇、丙酮、乙醚、甘油和植物油中。
PE在大气、阳光和氧的作用下易发生老化,具体表现为伸长率和耐寒性降低,力学性能和电性能下降,并逐渐变脆、产生裂纹,最终丧失使用性能。
为了防止PE的氧化降解,便于贮存、加工和应用,一般使用的PE原料在合成过程中已加入了稳定剂,可满足一般的加工和使用要求。
如需进一步提高耐老化性能,可在PE中添加抗氧剂和光稳定剂等。
6.卫生性
PE分子链主要由碳、氢构成,本身毒性极低,但为了改善PE性能,在聚合、成型加工和使用中往往需添加抗氧剂和光稳定剂等塑料助剂,可能影响到它的卫生性。
树脂生产厂家在聚合时总是选用无毒助剂,且用量极少,一般树脂不会受到污染。
PE长期与脂肪烃、芳香烃、卤代烃类物质接触容易引起溶胀,PE中有些低相对分子质量组分可能会溶于其中,因此,长期使用PE容器盛装食用油脂会产生一种蜡味,影响食用效果。
1.1.2聚乙烯的分类
聚乙烯的生产方法不同,其密度及熔体流动速率也不同。
按密度大小主要分为低密度聚乙烯(LDPE)、线型低密度聚乙烯(LLDPE)、中密度聚乙烯(MDPE)、高密度聚乙烯(HDPE)。
其中线性低密度聚乙烯属于低密度聚乙烯中的一种,是工业上常用的聚乙烯,其他分类法有时把MDPE归类于HDPE或LLDPE。
按相对分子质量可分为低相对分子质量聚乙烯、普通相对分子质量聚乙烯、超高相对分子质量聚乙烯。
按生产方法可分为低压法聚乙烯、中压法聚乙烯和高压法聚乙烯。
1.低密度聚乙烯
英文名称:
Lowdensitypolyethylene,简称LDPE
低密度聚乙烯,又称高压聚乙烯。
无味、无臭、无毒、表面无光泽、乳白色蜡状颗粒,密度0.910~0.925g/cm3,质轻,柔性,具有良好的延伸性、电绝缘性、化学稳定性、加工性能和耐低温性(可耐-70℃),但力学强度、隔湿性、隔气性和耐溶剂性较差。
分子结构不够规整,结晶度较低(55%~65%),熔点105~115℃。
LDPE可采用热塑性成型加工的各种成型工艺,如注射、挤出、吹塑、旋转成型、涂覆、发泡工艺、热成型、热风焊、热焊接等,成型加工性好。
主要用作农膜、工业用包装膜、药品与食品包装薄膜、机械零件、日用品、建筑材料、电线、电缆绝缘、吹塑中空成型制品、涂层和人造革等。
2.高密度聚乙烯
英文名称:
HighDensityPolyethylene,简称HDPE
高密度聚乙烯,又称低压聚乙烯。
无毒、无味、无臭,白色颗粒,分子为线型结构,很少有支化现象,是典型的结晶高聚物。
力学性能均优于低密度聚乙烯,熔点比低密度聚乙烯高,约125~137℃,其脆化温度比低密度聚乙烯低,约-100~-70℃,密度为0.941~0.960g/cm3。
常温下不溶于一般溶剂,但在脂肪烃、芳香烃和卤代烃中长时间接触时能溶胀,在70℃以上时稍溶于甲苯、醋酸中。
在空气中加热和受日光影响发生氧化作用。
能耐大多数酸碱的侵蚀。
吸水性小,具有良好的耐热性和耐寒性,化学稳定性好,还具有较高的刚性和韧性,介电性能、耐环境应力开裂性亦较好。
HDPE可采用注射、挤出、吹塑、滚塑等成型方法,生产薄膜制品、日用品及工业用的各种大小中空容器、管材、包装用的压延带和结扎带,绳缆、鱼网和编织用纤维、电线电缆等。
3.线性低密度聚乙烯
英文名称:
LinearLowDensityPolyethylene,简称LLDPE
线形低密度聚乙烯被认为是“第三代聚乙烯”的新品种,是乙烯与少量高级α-烯烃(如丁烯-1、己烯-1、辛烯-1、四甲基戊烯-1等)在催化剂作用下,经高压或低压聚合而成的一种共聚物,为无毒、无味、无臭的乳白色颗粒,密度0.918~0.935g/cm3。
与LDPE相比,具有强度大、韧性好、刚性大、耐热、耐寒性好等优点,且软化温度和熔融温度较高,还具有良好的耐环境应力开裂性,耐冲击强度、耐撕裂强度等性能。
并可耐酸、碱、有机溶剂等。
LLDPE可通过注射、挤出、吹塑等成型方法生产农膜、包装薄膜、复合薄膜、管材、中空容器、电线、电缆绝缘层等。
由于不存在长支链,LLDPE的65%~70%用于制作薄膜。
4.中密度聚乙烯
英文名称:
Mediumdensitypolyethylene,简称MDPE
中密度聚乙烯是在合成过程中用α-烯烃共聚,控制密度而成。
MDPE的密度为0.926~0.953g/cm3,结晶度为70%~80%,平均相对分子质量为20万,拉伸强度为8~24MPa,断裂伸长率为50%~60%,熔融温度126~135℃,熔体流动速率为0.1~35g/10min,热变形温度(0.46MPa)49~74℃。
MDPE最突出的特点是耐环境应力开裂性及强度的长期保持性。
MDPE可用挤出、注射、吹塑、滚塑、旋转、粉末成型加工方法,生产工艺参数与HDPE和LDPF相似,常用于管材、薄膜、中空容器等。
5.超高相对分子质量聚乙烯
英文名称:
ultra-highmolecularweightpolyethylene,简称UHMWPE
超高相对分子质量聚乙烯冲击强度高,耐疲劳,耐磨,是一种线型结构的具有优异综合性能的热塑性工程塑料。
其相对分子质量达到300~600万,密度0.936~0.964g/cm3,热变形温度(0.46MPa)85℃,熔点130~136℃。
UHMWPE因相对分子质量高而具有其他塑料无可比拟的优异性能,如耐冲击、耐磨损、自润滑性、耐化学腐蚀等性能,广泛应用于机械、运输、纺织、造纸、矿业、农业、化工及体育运动器械等领域,其中以大型包装容器和管道的应用最为广泛。
另外,由于超高相对分子质量聚乙烯优异的生理惰性,已作为心脏瓣膜、矫形外科零件、人工关节等在临床医学上使用,而且,超高相对分子质量聚乙烯耐低温性能优异,在-40℃时仍具有较高的冲击强度,甚至可在-269℃下使用。
超高相对分子质量聚乙烯纤维的复合材料在军事上已用作装甲车辆的壳体、雷达的防护罩壳、头盔等;体育用品上已制成弓弦、雪橇和滑水板等。
由于超高相对分子质量聚乙烯熔融状态的粘度高达108Pa·s,流动性极差,其熔体流动速率几乎为零,所以很难用一般的机械加工方法进行加工。
近年来,通过对普通加工设备的改造,已使超高相对分子质量聚乙烯由最初的压制-烧结成型发展为挤出、吹塑和注射成型以及其他特殊方法的成型。
6.茂金属聚乙烯
茂金属聚乙烯(mPE)是近年来迅速发展的一类新型高分子树脂,其相对分子质量分布窄,分子链结构和组成分布均一,具有优异的力学性能和光学性能,已被广泛应用于包装、电气绝缘制品等。
1.1.3聚乙烯的成型加工
PE的熔体粘度比PVC低,流动性能好,不需加入增塑剂已具有很好的成型加工性能。
前文已介绍了各类聚乙烯可采用的成型加工方法,下面主要介绍在成型过程中应注意的几个问题。
①聚乙烯属于结晶性塑料,吸湿小,成型前不需充分干燥,熔体流动性极好,流动性对压力敏感,成型时宜用高压注射,料温均匀,填充速度快,保压充分。
不宜用直接浇口,以防收缩不均,内应力增大。
注意选择浇口位置,防止产生缩孔和变形。
②PE的热容量较大,但成型加工温度却较低,成型加工温度的确定主要取决于相对分子质量、密度和结晶度。
LDPE在180℃左右,HDPE在220℃左右,最高成型加工温度一般不超过280℃。
③熔融状态下,PE具有氧化倾向,因而,成型加工中应尽量减少熔体与空气的接触及在高温下的停留时间。
④PE的熔体粘度对剪切速率敏感,随剪切速率的增大下降得较多。
当剪切速率超过临界值后,易出现熔体破裂等流动缺陷。
⑤制品的结晶度取决于成型加工中对冷却速率的控制。
不论采取快速冷却还是缓慢冷却,应尽量使制品各部分冷却速率均匀一致,以免产生内应力,降低制品的力学性能。
⑥收缩范围和收缩值大(一般成型收缩率为1.5%~5.0%),方向性明显,易变形翘曲,冷却速度宜慢,模具设冷料穴,并有冷却系统。
⑦软质塑件有较浅的侧凹槽时,可强行脱模。
1.1.4聚乙烯的改性
聚乙烯属非极性聚合物,与无机物、极性高分子相容性弱,因此其功能性较差,采用改性可提高PE的耐热老化性、高速加工性、冲击强度、粘接性、生物相容性等性质。
常用的改性方法包括物理改性和化学改性。
1.物理改性
物理改性是在PE基体中加入另一组分(无机组分、有机组分或聚合物等)的一种改性方法。
常用的方法有增强改性、共混改性、填充改性。
(1)增强改性增强改性是指填充后对聚合物有增强效果的改性。
加入的增强剂有玻璃纤维、碳纤维、石棉纤维、合成纤维、棉麻纤维、晶须等。
自增强改性也属于增强改性的一种。
①自增强改性。
所谓自增强就是使用特殊的加工成型方法,使得材料内部组织形成伸直链晶体,材料内部大分子晶体沿应力方向有序排列,材料的宏观强度得到大幅度提高,同时分子链有序排列将使结晶度提高,从而使材料的强度进一步提高,由于所形成的增强相与基体相的分子结构相同,因而不存在外增强材料中普遍存在的界面问题。
如采用超高相对分子质量聚乙烯(UHMPE)纤维增强LDPE,在加热加压成型的条件下,可以形成良好的界面,最大限度发挥基体和纤维的强度。
②纤维增强改性。
纤维增强聚合物基复合材料由于具有比强度高、比刚度高等优点而得到广泛应用。
如采用经KH-550偶联剂处理的长玻璃纤维(LGF)与PE复合制备的PE/LGF复合材料,当LGF加入量为3O%(质量分数)、长度约为35mm时,复合材料的拉伸强度和冲击强度分别为52.5MPa和52kJ/m。
③晶须改性。
晶须的加入能够大幅度提高HDPE材料的力学性能,包括短期力学性能及耐长期蠕变性能。
晶须对HDPE材料的增强作用主要归因于它们之间的良好界面粘接,同时刚性的晶须则能够承担较大的外界应力使复合材料的模量得到提高。
④纳米粒子增强改性。
少量无机刚性粒子填充PE可同时起到增韧与增强的作用。
如将表面处理过的纳米SiO2粒子填充mLLDPE-LDPE,SiO2纳米粒子均匀分散于基材中,与基材形成牢固的界面结合,当填充质量分数为2%时,拉伸强度、断裂伸长率分别提高了13.7MPa和174.9%。
(2)共混改性共混改性主要目的是改善PE的韧性、冲击强度、粘接性、高速加工性等各种缺陷,使其具有较好的综合性能。
共混改性主要是向PE基体中加入另一种聚合物,如塑料类、弹性体类等聚合物,以及不同种类的PE之间进行共混。
①PE系列的共混改性。
单一组分的PE往往很难满足加工要求,而通过不同种类PE之间的共混改性可以获得性能优良的PE材料。
如通过LDPE与LLDPE共混,解决了LDPE因大量添加阻燃剂和抗静电剂等助剂造成力学性能急剧降低的问题;LLDPE与HDPE共混后可以提高产品的综合性能。
②PE与弹性体的共混改性。
弹性体具有低的表面张力、较强的极性、突出的增韧作用,因此与PE共混后,既能保持PE的原有性能,同时也可以制备出具有综合优良性能的PE。
如LDPE-聚烯烃弹性体(POE)共混物,当POE的质量分数为3O%时,共混体系的拉伸强度达到最大值,为21.5MPa。
③PE与塑料的共混改性。
聚乙烯具有良好的韧性,但制品的强度和模量较低,与工程塑料等共混可提高复合体系的综合力学性能。
但PE和这类高聚物的界面问题也是影响其共混物性能的主要原因,因此通常需要加入界面相容剂以提高共混物的力学性能。
(3)填充改性填充改性是在PE基质中加入无机填料或有机填料,一方面可以降低成本达到增重的目的,另一方面可提高PE的功能性,如电性能、阻燃性能等,但同时对复合材料的力学性能和加工性能带来一定程度的影响。
无论是无机填料还是有机填料,填料与PE基体的相容性和界面粘接强度是PE填充改性必须面临的问题,而PE是非极性化合物,与填料相容性差,因此,必须对填料进行表面处理。
填料的表面处理
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 计算 程序 流体力学 对流 方程 有限 差分法 Lax 格式 迎风 FTCS
![提示](https://static.bdocx.com/images/bang_tan.gif)