纤维元素的质量守恒可以表示为
dzd(ρfufAfYs)=−πdfms˙(22.1)
式中
ρf=fiber density
uf=fiber velocity vector
Af= 与流动方向平行的纤维表面积矢量
Ys= 纤维中溶剂 s 的质量分数
m˙s′′= 溶剂 s 的蒸发质量流率
df= 纤维直径
z= 沿纤维的坐标
m˙s′′ 是使用薄膜理论计算的。
m˙s′′=Mscβln(1−ψs,I1−ψs,g)(22.2)
其中,
β= 从适当的相关性估计的质量转移系数(参见动量、热量和质量转移的相关性(第952页))
Ms= 溶剂的分子量
c= 周围气体的摩尔浓度
ψs,g= 周围气体中溶剂蒸气的摩尔分数
在纤维表面,气体中溶剂的摩尔分数 ψs,l 与纤维中溶剂的质量分数 Ys 通过Flory [185](第1067页)给出的气液平衡方程相关联,
ψs,I=pps, vap Yse1−Ys+χ(1−Ys)2(22.3)
何处
χ=Flory-Huggins parameter
p= 周围流动中的绝对压力
ps,vap= 溶剂的饱和蒸汽压
这些方程仅在选择了干纺纤维时使用。
纤维的形成基于在卷取点施加于纤维的张力,这些张力导致纤维的拉伸和延伸。
对纤维微元进行力的平衡分析,得出纤维动量变化的方程。
dzd(ρfufufAf)=dzdF+Ffriction +Fgravitation (22.4)
由于纤维的加速度、与周围气体的摩擦力以及重力作用,纤维中的张力会发生改变。
与周围气体的摩擦力是通过计算得出的。
Ffriction =21ρcf,axπdfuf−upar (uf−upar )(22.5)
其中
ρ= 气体密度
cf,ax= 沿纤维轴向的摩擦系数
upar = 沿纤维方向的气体速度
重力由以下公式计算得出
Fgravitation =ρf4πdf2g⋅nf(22.6)
其中,nf 表示纤维元素的方向向量。
拉伸力 F 与应力张量的分量相关联。
F=Af(τzz−τrr)(22.7)
忽略粘弹性效应并假设牛顿流动,便可以得到
τzz=2ηfdzduf(22.8)
τrr=−ηfdzduf(22.9)
导致
F=3Afηfdzduf(22.10)
通过将零剪切粘度 ηf 乘以三,来估算拉伸粘度。
为了计算沿纺丝线方向的纤维温度,对进入和传递至微分纤维单元的热焓进行平衡处理。
dzd(ρfuf⋅Afhf)=dzd(λfAfdzdTf)+πdf(α(T−Tf)−m˙s′′hs,ν)(22.11)
+Q˙viscous heating +Q˙radiation, abs −Q˙radiation, emission
式中
hf=fiber enthalpy
λf=fiber thermal conductivity
Tf=fiber temperature
hs,v=enthalpy of the solvent vapor
α=heat transfer coefficient
在熔体纺丝过程中,m˙s′′ 为零,因为没有质量转移。由于粘性加热产生的热量项是从圆柱形流动的流体力学中推导出来的。
Q˙viscous heating =4πdf2(4(dzduf)2−32m˙s′′′2)(22.12)
辐射热交换由最后两项考虑。
Q˙radiation, abs =dfεfG(22.13)
Q˙radiation, emission =πdfεfσTf4(22.14)
式中
G=thermal irradiation
εf=fiber’s emissivity
σ=Boltzman constant
纤维焓值 hf 与纤维温度 Tf 之间的关系如下:
hf=∫TrefTf((1−Ys)Cpp+YsCps)dT(22.15)
它利用了聚合物的比热容 Cpp 和纤维中溶剂的比热容 Cps。
在给定温度 Tv 下,溶剂蒸气的焓值取决于在蒸发温度 Tvap 下给出的蒸发热 Δhs,并通过以下公式计算得出:
hs,v=∫TrefTvapCps,ldT+Δhs∣Tvap+∫TvapTvCps,vdT(22.16)
以下是翻译后的中文文本:
Cps,l= 溶剂液体的比热容
Cps,v= 溶剂蒸气的比热容