分享到:

从裂解汽油中回收C_9芳烃的模拟探讨

前言裂介汽油是煤柴油高温裂介后的认~q。部份的液相产品,约有普来自汽油气提塔(包括压缩段间冷凝液,,主要为c。以下较重组分。约有告来自脱丁烷塔塔底,主要为c。以_匕较轻组分。 ’从某厂装置试运转数据看来,最明显的是c。芳烃的变化,仅为原设计值的一半左右,且其中三甲苯含量显著减少,茹和茹满,尤其是丙基苯增加较多,因此给‘c。芳烃回收增加了一定困难。这种严重差距,就对芳烃回收提出了怎样才能回收以作为歧化原料的研究课题。 要解决这个问题,_必须搞清下列几个疑点: (1)是否由于原料不同而引起的? (2)是否由于设计或操作不合理而引起的? (3)C。芳烃究竟在哪里损失的或者如何才能回收到C。芳烃? 对上述这些问题的回答,我们和厂方商定:用实验和计算机模拟手段予以逐步解决。本文是用计算机对生产装置模拟计算的结果,仅为本课题的开始阶段。数学模型和框图 对于多组分连续蒸馏塔的计算,是通过物料和能量衡算以及每块理论板上的相平衡方程来建立数学模型...  (本文共12页) 阅读全文>>

《数学学习与研究》2014年03期
数学学习与研究

加边三对角矩阵特征值及其应用

引言加边三对角矩阵在图像处理、随机过程和数值分析等学科中有广泛的应用[1]-[3],该类矩阵的有关研究一直受到人们的关注[4],本文考虑加边三对角矩阵的特征值问题,并给出了一个例子.定义1 n-阶对称矩阵形式如下:Βn=α1b1b2……bn-2bn-1b1α2β2……0 0b2β2α3…0 0bn-20 0…βn-2αn-1βn-1bn-10 0…0βn-1α??????????????????n∈Rn×n.MERGEFORMAT(1)若βi≥0(i=2,3,…,n-1),则称Βn为加边非负雅克比矩阵;若βi≠0(i=2,3,…,n-1),则称Βn为加边不可约三对角矩阵.本文首先给出了有关加边对角矩阵的若干性质,由此引出加边三对角矩阵的特征值问题.引理1[5]假设λ20,满足Ax=ρ(A)x.定理2若加边三对角矩阵Βn的元素bi0(i=1,2,…,n-1),μn为其最大特征值,那么存在实向量x,使得...  (本文共2页) 阅读全文>>

《计算机工程与应用》2012年17期
计算机工程与应用

求分块周期三对角矩阵逆矩阵的新算法

1引言分块周期三对角矩阵是一类重要的特殊矩阵,它在科学与工程计算方面应用广泛,研究求其逆矩阵的相关算法一直是人们研究的热点问题,近年来得到不少求该类矩阵逆的新算法[1-4]。在文献[1]中,IKEBEY证明了三对角矩阵的逆矩阵可以由四组向量确定。已有学者通过三对角矩阵的LU和UL分解,给出了求此四组向量的方法,进而给出了求三对角矩阵逆矩阵的显式算法。本文在前人工作的基础上,利用Sheman-Morrison-Woodbuy公式,通过选用两个参数矩阵得到求分块周期三对角矩阵的逆矩阵的新算法,并据此得到三对角矩阵和对称三对角矩阵求逆矩阵的新算法,新算法计算量要比传统算法小,计算效率有显著提高。设分块周期三对角矩阵有如下形式:A=è÷÷÷÷÷÷÷a1b100tc2a2b2000c3a300000an-1bn-1s00cnan(1)其中,A的元素aibicist都是m阶方阵。若s=t=0,则矩阵A为分块三对角...  (本文共3页) 阅读全文>>

《炼油技术与工程》2006年04期
炼油技术与工程

改进三对角矩阵法在粗甲醇精馏塔计算中的应用

甲醇精馏是用精馏的方法除去粗甲醇中的各 类杂质,以达到一定的质量标准,获得工业级精甲 醇的过程。粗甲醇中含有多种组分〔‘],因为它们 的沸点各异,工业上多采用二塔或三塔流程。二 塔常压流程为最早,也是目前使用最多的技术,基 于节能考虑,加压的三塔流程得到了迅速发展。 甲醇精馏塔向多塔式的发展,使其计算复杂 化,对精馏塔的计算机模拟有必要从简化计算发 展到精确模拟。同时,对精馏塔的精确模拟计算, 也是寻找更为合理的工艺操作条件的手段之一。 1数学模型 本文采用5组分体系描述粗甲醇精馏系统, 即甲醇一水一正辛烷一二甲醚一l一丁醇。其中1一丁醇 是多元醇的代表组分,正辛烷是高级烃的代表组 分。由于甲醇精馏体系是多元组分的汽液平衡, 至今尚未见实验数据的报导。实践证明:UMFAC 法可适用于粗甲醇体系的汽液平衡计算川,且计 算结果与生产实际相符。本文采用UNIFAC法计 算粗甲醇体系的活度系数,采用维里方程计算气 体的逸度系数,而精馏塔...  (本文共5页) 阅读全文>>

《大学数学》2003年03期
大学数学

一类对称三对角矩阵的合同对角化算法的实现

1 提出问题我们先考虑对称三对角矩阵A=a b 0… 0b a b… 00 b a 00 0 b0 0… b a n× n,  a≥ 2 | b| 0 ,对应的二次型X′AX=a( x21+ x22 +… + x2n) + 2 b( x1x2 + x1x3+… + xn- 1xn) .要对其标准化的方法很多 ,我们介绍一种比较方便易懂的方法 ,该方法便于在计算机上编程实现 .2 解决问题在矩阵 A中 ,用第一行乘上 - ba加到第二行 ,得a b 0… 00 a- b2a b… 00 b a 00 0 b0 0… b a n× n.令λ1=a,λ2 =a- b2a,则用第二行乘上 - bλ2加到第三行得主对角线上的元素为λ3=a- b2λ2.如此反复 ,可得λn=a- b2λn- 1,其中λi( i=1 ,2 ,… ,n)分别为第 i行上主对角线上的元素 ,则原矩阵 A化为一个上三角矩阵 .同理 ,该化简过程作用于上...  (本文共3页) 阅读全文>>

《广东财经职业学院学报》2002年02期
广东财经职业学院学报

实对称三对角矩阵的特征值的估算

形如δ1 ε1 0……… 0ε1 δ2 ε2 0ε2 δ3   0 εn- 10……… 0εn- 1 δn的实方阵称为n阶实对称三对角矩阵 ,本文以字母T特记之 对T的特征值进行估算是具有实际应用意义的 ,本文的主要目的是希望找到一种对T的特征值进行有效估算的方法 为了叙述的方便 ,我们约定以下几个符号 :1 .m表示集合 {1 ,2 ,3,… ,m};2 .T 1 2…k1 2…k 表示方阵T的k阶顺序主子阵 ,以 φk(λ)表示其特征多项式 ;3.Vc 表示实系数多项式序列fn(x) ,fn - 1 (x) ,… ,f1 (x) ,f0 (x)在x =c时的变号数 ;4.ρ(A)表示n阶方阵A =(aij)的谱半径 ,即若λ1 ,λ2 ,λ3,… ,λn 是A在复数域的特征值 ,则 :  ρ(A) =maxi {|λi| },i∈n .一、关于T的特征值的个数问题T作为特...  (本文共4页) 阅读全文>>