CN107121363A - 一种结晶工艺过程颗粒尺寸分布预测的方法 - Google Patents

一种结晶工艺过程颗粒尺寸分布预测的方法 Download PDF

Info

Publication number
CN107121363A
CN107121363A CN201710173148.9A CN201710173148A CN107121363A CN 107121363 A CN107121363 A CN 107121363A CN 201710173148 A CN201710173148 A CN 201710173148A CN 107121363 A CN107121363 A CN 107121363A
Authority
CN
China
Prior art keywords
mrow
msub
particle size
mfrac
size distribution
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201710173148.9A
Other languages
English (en)
Inventor
徐三
李秀喜
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
South China University of Technology SCUT
Original Assignee
South China University of Technology SCUT
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by South China University of Technology SCUT filed Critical South China University of Technology SCUT
Priority to CN201710173148.9A priority Critical patent/CN107121363A/zh
Publication of CN107121363A publication Critical patent/CN107121363A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/02Investigating particle size or size distribution
    • G01N15/0205Investigating particle size or size distribution by optical means

Landscapes

  • Chemical & Material Sciences (AREA)
  • Dispersion Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Organic Low-Molecular-Weight Compounds And Preparation Thereof (AREA)

Abstract

本发明属于结晶过程控制技术领域,公开了一种结晶工艺过程颗粒尺寸分布预测的方法。所述方法为:通过超声颗粒粒度分析仪获取实时颗粒尺寸分布数据,衰减全反射傅里叶红外光谱获取实时浓度数据,采用矩量法计算成核和生长速率,然后带入到如下常微分方程组中,运用MATLAB中的ODE45求解器计算,最终获得结晶工艺过程的颗粒尺寸分布。本发明的预测方法集成了矩量法和离散法的优点,提高计算求解的效率,对计算机的硬件要求降低,具有良好的成本优势;并解决了颗粒尺寸分布预测的离散问题,使结晶工艺过程得以从实验室向实际生产过程进行放大,应用范围广泛。

Description

一种结晶工艺过程颗粒尺寸分布预测的方法
技术领域
本发明属于结晶过程控制技术领域,具体涉及一种结晶工艺过程颗粒尺寸分布预测的方法。
背景技术
结晶作为化学工业中重要的单元操作,它能够以低成本获取高附加值的产品。而结晶产品质量的关键指标之一是颗粒尺寸分布,对其准确预测能极大改善产品的质量。目前在结晶领域,用于描述颗粒物制造过程并预测颗粒尺寸分布的粒数衡算模型是一个在时间域和尺寸域上多尺度的偏微分方程。但是它的解析解只有在非常理想状态下才可获得,因此如何寻求高效的数值解成为了当今的一个技术难点。文献中报道的关于此类方程的求解模式主要有矩量法、离散法等。单纯的采用矩量法或者离散法所产生的非线性模型问题非常巨大,无论求解效率还是在求解可靠性上都受到巨大挑战,特别是工艺放大更是难以办到。
粒数恒算模型代表着一些我们感兴趣的沿着一定的性质坐标分布的物理量,比如种群密度,它是关于个体在一定空间范围内的数量分布(亦或者质量、体积分布)(Wokd,S.Nonlinear partial least squares modellingⅡ.Spline inner relation[J].Chemometrics and Intelligent Laboratory Systems,1992,14(1):71-84)。该模型类似于描述某种机械系统中相空间概率守恒的经典刘维尔方程。在上个世纪,粒数恒算模型就已经被建立来描述各种颗粒物体系:生物细胞、液滴和气泡、造粒物、晶体、聚合物、凝溶胶甚至宇宙中天体颗粒物。到了20世纪60年代,粒数恒算模型开始用于求解化学工程问题,其通用形式:
粒数恒算模型描述的是一个粒子相空间而不是所谓简单的欧式空间,粒子相坐标由内部坐标和外部坐标构成,前者规定粒子的质量指标比如尺寸、晶型等,后者指定粒子在过程单元中运动的物理位置。式1-1中n代表粒子数密度函数,它是单位体积晶浆中单位粒度的晶体的个数,是颗粒尺寸的函数(式1-2给出了其具体定义式,以人口作比喻,晶体粒数密度相当于人口密度,即等价于某一地区在某个年龄范围(比如20-30岁)内的人数,其单位为#·m-1·m-3另外在上式1-3中引入的晶体粒数N,其单位是#·m-3。另外,式1-1中的ξ代表相空间变量(包括物理空间坐标变量和区分粒子性质的内部坐标变量,如粒径),t代表时间,h是体系中新粒子引入时的净速率,在一个完整的系统中它一般由三大机理构成:成核、聚结和破碎。当加入这三种机理时就会形成一个复杂的积分-微分方程体系。即便提供合适的初边值条件,在数学分析层面上,这个完整的表达式也是难以处理的。尹秋响等(尹秋响,仲维正,王胜春,王静康.小粒度维生素C冷却结晶动力学[J].天津大学学报,2004,33(1):23-27)采用矩量法对粒数衡算模型求解,获得工业结晶中常用的成核和生长经验式:
式中mj表示粒数密度的j阶矩。
Ali Abbas等(Ali Abbas,Jose A.Romagnoli.Multiscale modeling,simulationand validation of batch cooling crystallization)采用离散法对粒数衡算模型进行离散降阶,获得一系列容易求解的常微分方程组:
式中φ表示相邻尺寸区间长度,φ=Li-Li-1
上述方法中矩量法一方面无法求解与尺寸相关的晶体生长过程,另一方面从单一的矩量方程组无法预测不同的颗粒尺寸分布。离散法由于生长和成核速率模型存在,极大增加模型预测的计算负荷和数值稳定性。
发明内容
为了解决以上现有技术的缺点和不足之处,本发明的目的在于提供一种结晶工艺过程颗粒尺寸分布预测的方法。
本发明目的通过以下技术方案实现:
一种结晶工艺过程颗粒尺寸分布预测的方法,包括如下步骤:
通过超声颗粒粒度分析仪获取实时颗粒尺寸分布数据,衰减全反射傅里叶红外光谱获取实时浓度数据,采用矩量法计算成核和生长速率,然后带入到如下常微分方程组(1-7)~(1-9)中,运用MATLAB中的ODE45求解器计算,最终获得结晶工艺过程的颗粒尺寸分布;
利用上述方法对阿司匹林结晶工艺过程中产品颗粒尺寸分布进行预测的方法,具体步骤如下:
将250g阿司匹林溶解在1000mL的乙醇溶液中,以0.5℃/min的降温速率进行降温,同时打开超声颗粒粒度分析仪和衰减全反射傅里叶红外光谱,在线采集颗粒尺寸分布和浓度数据,采用矩量法计算阿司匹林在乙醇中冷却结晶的成核和生长速率,然后带入常微分方程组(1-7)~(1-9)中,运用MATLAB中的ODE45求解器计算,最终获得阿司匹林在乙醇中结晶的颗粒尺寸分布。
本发明的预测方法具有如下优点及有益效果:
(1)本发明的预测方法集成了矩量法和离散法的优点,提高计算求解的效率,对计算机的硬件要求降低,具有良好的成本优势;
(2)本发明解决了颗粒尺寸分布预测的离散问题,使结晶工艺过程得以从实验室向实际生产过程进行放大;
(3)本发明采用描述颗粒物制造过程(沉淀过程、乳液聚合过程等)的粒数衡算模型进行预测,应用范围更广。
附图说明
图1为本发明实施例对法对阿司匹林结晶工艺过程中产品颗粒尺寸分布进行预测的装置结构示意图;图中编号说明如下:①-衰减全反射傅里叶红外光谱,②-超声颗粒粒度分析仪,③-浊度仪,④-循环水浴控制器,⑤-结晶器,⑥-搅拌桨。
图2为本发明实施例中结晶过程的降温曲线图。
图3为本发明实施例中超声颗粒粒度分析仪采集的颗粒尺寸分布图。
图4为本发明实施例中衰减全反射傅里叶红外光谱采集计算的浓度变化曲线图。
图5为本发明实施例中阿司匹林在乙醇中冷却结晶的成核速率随过饱和度的变化曲线图。
图6为本发明实施例中阿司匹林在乙醇中冷却结晶的生长速率随过饱和度的变化曲线图。
图7为本发明实施例最终获得的阿司匹林在乙醇中结晶的颗粒尺寸分布图。
具体实施方式
下面结合实施例及附图对本发明作进一步详细的描述,但本发明的实施方式不限于此。
实施例
采用图1所示的装置对法对阿司匹林结晶工艺过程中产品颗粒尺寸分布进行预测:
图1装置由衰减全反射傅里叶红外光谱①,超声颗粒粒度分析仪②,浊度仪③,循环水浴控制器④,结晶器⑤和搅拌桨⑥组成。搅拌桨搅拌条件(150rpm/min)下在结晶器中将250g阿司匹林溶解在1L的乙醇溶液中(溶解温度为35℃),循环水浴控制器,在体系降温之前设置好降温程序(第一阶段,40℃恒温搅拌30分钟;第二阶段40℃降温至35℃,用时10分钟;第三阶段35℃恒温搅拌30分钟;第四阶段35℃降温至25℃,用时20分钟;第五阶段,25℃恒温搅拌60分钟,其结晶过程的降温曲线如图2所示)。当循环水浴控制器第一阶段结束时打开超声颗粒粒度分析仪和衰减全反射傅里叶红外光谱在线采集颗粒尺寸分布(如图3所示)和浓度数据(图4)。同时打开浊度仪,浊度仪在线监控成核时刻(表现为浊度值的突变)。采用矩量法计算阿司匹林在乙醇中冷却结晶的成核和生长速率(图5和图6),然后带入常微分方程组(1-7)~(1-9)中,运用MATLAB中的ODE45求解器计算。最终获得10分钟、90分钟及150分钟阿司匹林在乙醇中结晶的颗粒尺寸分布模拟值(如图7所示)。
上述实施例为本发明较佳的实施方式,但本发明的实施方式并不受上述实施例的限制,其它的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化,均应为等效的置换方式,都包含在本发明的保护范围之内。

Claims (2)

1.一种结晶工艺过程颗粒尺寸分布预测的方法,其特征在于包括如下步骤:
通过超声颗粒粒度分析仪获取实时颗粒尺寸分布数据,衰减全反射傅里叶红外光谱获取实时浓度数据,采用矩量法计算成核和生长速率,然后带入到如下常微分方程组(1-7)~(1-9)中,运用MATLAB中的ODE45求解器计算,最终获得结晶工艺过程的颗粒尺寸分布;
<mrow> <mfrac> <mrow> <msub> <mi>dN</mi> <mn>1</mn> </msub> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <mo>=</mo> <mo>-</mo> <mfrac> <mi>G</mi> <mrow> <mn>2</mn> <msub> <mi>&amp;phi;</mi> <mn>1</mn> </msub> </mrow> </mfrac> <msub> <mi>N</mi> <mn>1</mn> </msub> <mo>+</mo> <mi>B</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <mfrac> <mrow> <msub> <mi>dN</mi> <mi>i</mi> </msub> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <mo>=</mo> <mo>-</mo> <mfrac> <mi>G</mi> <mrow> <mn>2</mn> <msub> <mi>&amp;phi;</mi> <mi>i</mi> </msub> </mrow> </mfrac> <msub> <mi>N</mi> <mi>i</mi> </msub> <mo>+</mo> <mfrac> <mi>G</mi> <mrow> <mn>2</mn> <msub> <mi>&amp;phi;</mi> <mrow> <mi>i</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mrow> </mfrac> <msub> <mi>N</mi> <mrow> <mi>i</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>,</mo> <mi>i</mi> <mo>=</mo> <mn>2</mn> <mo>,</mo> <mn>3</mn> <mo>...</mo> <mo>,</mo> <mi>n</mi> <mo>-</mo> <mn>1</mn> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <mfrac> <mrow> <msub> <mi>dN</mi> <mi>n</mi> </msub> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <mo>=</mo> <mo>-</mo> <mfrac> <mi>G</mi> <mrow> <mn>2</mn> <msub> <mi>&amp;phi;</mi> <mi>n</mi> </msub> </mrow> </mfrac> <msub> <mi>N</mi> <mi>n</mi> </msub> <mo>+</mo> <mfrac> <mi>G</mi> <mrow> <mn>2</mn> <msub> <mi>&amp;phi;</mi> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mrow> </mfrac> <msub> <mi>N</mi> <mrow> <mi>n</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mn>9</mn> <mo>)</mo> </mrow> <mo>.</mo> </mrow>
2.采用权利要求1所述的方法对阿司匹林结晶工艺过程中产品颗粒尺寸分布进行预测,具体步骤如下:
将250g阿司匹林溶解在1000mL的乙醇溶液中,以0.5℃/min的降温速率进行降温,同时打开超声颗粒粒度分析仪和衰减全反射傅里叶红外光谱,在线采集颗粒尺寸分布和浓度数据,采用矩量法计算阿司匹林在乙醇中冷却结晶的成核和生长速率,然后带入常微分方程组(1-7)~(1-9)中,运用MATLAB中的ODE45求解器计算,最终获得阿司匹林在乙醇中结晶的颗粒尺寸分布。
CN201710173148.9A 2017-03-22 2017-03-22 一种结晶工艺过程颗粒尺寸分布预测的方法 Pending CN107121363A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710173148.9A CN107121363A (zh) 2017-03-22 2017-03-22 一种结晶工艺过程颗粒尺寸分布预测的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710173148.9A CN107121363A (zh) 2017-03-22 2017-03-22 一种结晶工艺过程颗粒尺寸分布预测的方法

Publications (1)

Publication Number Publication Date
CN107121363A true CN107121363A (zh) 2017-09-01

Family

ID=59718129

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710173148.9A Pending CN107121363A (zh) 2017-03-22 2017-03-22 一种结晶工艺过程颗粒尺寸分布预测的方法

Country Status (1)

Country Link
CN (1) CN107121363A (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109283104A (zh) * 2018-11-15 2019-01-29 北京理工大学 一种rdx制备过程中结晶液中产品粒径分布在线监测方法
CN110243825A (zh) * 2019-06-27 2019-09-17 青岛科技大学 一种基于在线图像法测量结晶过程溶解度、超溶解度及溶液浓度的方法
CN115684484A (zh) * 2022-10-20 2023-02-03 天津大学 一种预测晶体颗粒临界结块周期的方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000225302A (ja) * 1999-02-04 2000-08-15 Mitsubishi Chemicals Corp 粒度分布の特徴量の推定方法、粒子の製造方法及びこれらに用いる装置
JP2001056282A (ja) * 1999-06-08 2001-02-27 Solt Industry Center Of Japan 粒径分布測定装置
CN1696654A (zh) * 2005-06-06 2005-11-16 清华大学 一种颗粒粒径分布的测量方法
CN101769847A (zh) * 2010-01-20 2010-07-07 中昊晨光化工研究院 一种测量聚四氟乙烯粉体平均粒径及其粒径分布的方法
CN102460113A (zh) * 2009-05-11 2012-05-16 西安大略大学 监测介质的颗粒尺寸分布的超声法方法
CN105247342A (zh) * 2013-02-28 2016-01-13 尼尔·M·戴 用于确定颗粒尺寸的方法及设备

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000225302A (ja) * 1999-02-04 2000-08-15 Mitsubishi Chemicals Corp 粒度分布の特徴量の推定方法、粒子の製造方法及びこれらに用いる装置
JP2001056282A (ja) * 1999-06-08 2001-02-27 Solt Industry Center Of Japan 粒径分布測定装置
CN1696654A (zh) * 2005-06-06 2005-11-16 清华大学 一种颗粒粒径分布的测量方法
CN102460113A (zh) * 2009-05-11 2012-05-16 西安大略大学 监测介质的颗粒尺寸分布的超声法方法
CN101769847A (zh) * 2010-01-20 2010-07-07 中昊晨光化工研究院 一种测量聚四氟乙烯粉体平均粒径及其粒径分布的方法
CN105247342A (zh) * 2013-02-28 2016-01-13 尼尔·M·戴 用于确定颗粒尺寸的方法及设备

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
尹秋响,仲维正,王胜春,王静康: "小粒度维生素C冷却结晶动力学", 《天津大学学报》 *
江波: "阿司匹林结晶过程多尺度建模与模拟", 《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109283104A (zh) * 2018-11-15 2019-01-29 北京理工大学 一种rdx制备过程中结晶液中产品粒径分布在线监测方法
CN110243825A (zh) * 2019-06-27 2019-09-17 青岛科技大学 一种基于在线图像法测量结晶过程溶解度、超溶解度及溶液浓度的方法
CN110243825B (zh) * 2019-06-27 2021-09-10 青岛科技大学 一种基于在线图像法测量结晶过程溶解度、超溶解度及溶液浓度的方法
CN115684484A (zh) * 2022-10-20 2023-02-03 天津大学 一种预测晶体颗粒临界结块周期的方法

Similar Documents

Publication Publication Date Title
CN107121363A (zh) 一种结晶工艺过程颗粒尺寸分布预测的方法
CN111028893B (zh) 一种晶体生长预测方法及装置
Chung et al. Optimal model-based experimental design in batch crystallization
Mesbah et al. A control oriented study on the numerical solution of the population balance equation for crystallization processes
WO2021068901A1 (zh) 基于指数时间差分格式求解的材料微观结构演化模拟
Hemalatha et al. Multiobjective optimization and experimental validation for batch cooling crystallization of citric acid anhydrate
Sun et al. Temperature dependence of the solid-liquid interface free energy of Ni and Al from molecular dynamics simulation of nucleation
Liu et al. Dynamic relationships, regional differences, and driving mechanisms between economic development and carbon emissions from the farming industry: empirical evidence from rural China
McKenzie et al. Hybrid Monte Carlo technique for modeling of crystal nucleation and application to lithium disilicate glass-ceramics
Su et al. Modeling and Bayesian parameter estimation for semibatch pH‐shift reactive crystallization of l‐glutamic acid
Bonnin-Paris et al. Determination of the metastable zone width of glycine aqueous solutions for batch crystallizations
Makoveeva et al. Mathematical modeling of crystallization process from a supercooled binary melt
Xu et al. Multiscale modeling and simulation of directional solidification process of Ni-based superalloy turbine blade casting
Jiang et al. Fast prediction of the quasi phase equilibrium in phase field model for multicomponent alloys based on machine learning method
CN114819407A (zh) 湖泊蓝藻水华动态预测方法及装置
Dli et al. Intelligent control system architecture for phosphorus production from apatite-nepheline ore waste
Lou et al. Recessive transition mechanism of arable land use based on the perspective of coupling coordination of input–output: A case study of 31 provinces in China
Unno et al. Partial seeding policy for controlling the crystal quality in batch cooling crystallization
Chaudhury et al. Integrated population balance model development and validation of a granulation process
CN105260503B (zh) 结晶度与晶化温度关系模型建立方法
Matveev et al. Oscillating stationary distributions of nanoclusters in an open system
Simão et al. A case study for decentralized heat storage solutions in the agroindustry sector using phase change materials
CN103472721A (zh) 自适应机器学习的农药废液焚烧炉炉温最佳化系统及方法
Faria et al. Modelling agglomeration degree in sucrose crystallisation
Bajcinca et al. Method of moments over orthogonal polynomial bases

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20170901