CN113782106B - 一种用于cfd的高焓离解气体表面催化反应速率的求解方法 - Google Patents

一种用于cfd的高焓离解气体表面催化反应速率的求解方法 Download PDF

Info

Publication number
CN113782106B
CN113782106B CN202110862009.3A CN202110862009A CN113782106B CN 113782106 B CN113782106 B CN 113782106B CN 202110862009 A CN202110862009 A CN 202110862009A CN 113782106 B CN113782106 B CN 113782106B
Authority
CN
China
Prior art keywords
adsorption
reaction
atoms
rate
recombination
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.)
Active
Application number
CN202110862009.3A
Other languages
English (en)
Other versions
CN113782106A (zh
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.)
AERODYNAMICS NATIONAL KEY LABORATORY
Original Assignee
AERODYNAMICS NATIONAL KEY LABORATORY
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 AERODYNAMICS NATIONAL KEY LABORATORY filed Critical AERODYNAMICS NATIONAL KEY LABORATORY
Priority to CN202110862009.3A priority Critical patent/CN113782106B/zh
Publication of CN113782106A publication Critical patent/CN113782106A/zh
Application granted granted Critical
Publication of CN113782106B publication Critical patent/CN113782106B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16CCOMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C20/00Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
    • G16C20/10Analysis or design of chemical reactions, syntheses or processes
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/28Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16CCOMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C10/00Computational theoretical chemistry, i.e. ICT specially adapted for theoretical aspects of quantum chemistry, molecular mechanics, molecular dynamics or the like
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/08Fluids

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Computing Systems (AREA)
  • Physics & Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • General Physics & Mathematics (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Physics (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Computer Hardware Design (AREA)
  • Mathematical Optimization (AREA)
  • Fluid Mechanics (AREA)
  • Algebra (AREA)
  • Analytical Chemistry (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Physical Or Chemical Processes And Apparatus (AREA)

Abstract

本发明公开了一种用于CFD的高焓离解气体表面催化反应速率的求解方法,包括基于固体表面存在有限个数的吸附位点的理论,求解气相反应物在固体表面吸附位点上发生的四步反应的反应速率;其中,气相反应物在固体表面吸附位点上发生的四步反应分别为:气相原子的化学吸附、一个气相原子和一个吸附原子的ER复合、两个吸附原子的LH复合和吸附原子的热解附。通过本发明方法获得的反应速率可用以解决化学非平衡流场求解过程遇到的催化壁面表征困难、边界质量能量输运不准确的问题。

Description

一种用于CFD的高焓离解气体表面催化反应速率的求解方法
技术领域
本发明属于数值计算处理领域,尤其涉及可用于计算流体力学(CFD)的、基于多步反应机理和表面微细观参数来求解高焓离解气体各组分在固体材料表面的多相催化反应速率的方法。
背景技术
在飞行器化学非平衡绕流流场和表面气动热的数值计算中,需要考虑表面催化反应对流动和换热的影响,目前各组分的净复合系数是表征各个组分净反应速率的重要参数。
目前多采用给定复合系数为0-1之间的常数或Arrhenius形式的关联式,前者不考虑表面条件的影响,而后者也仅将表面温度的影响考虑在内,因此现有的计算方法无法体现真实的表面物理化学过程,带来表面催化反应速率的不确定性,进而导致流场和气动热计算的不确定性。
因此需要发展可用于计算流体力学(CFD)的、能够表征真实的表面反应过程的方法,包含反应动力学机理、各个反应步骤的速率确定方法等问题。
发明内容
本发明的目的在于,为克服现有技术缺陷,提供了一种用于CFD的高焓离解气体表面催化反应速率的求解方法,获得的反应速率将用以解决化学非平衡流场求解过程遇到的催化壁面表征困难、边界质量能量输运不准确的问题。
本发明目的通过下述技术方案来实现:
一种用于CFD的高焓离解气体表面催化反应速率的求解方法,包括基于固体表面存在有限个数的吸附位点的理论,求解气相反应物在固体表面吸附位点上发生的四步反应的反应速率;其中,气相反应物在固体表面吸附位点上发生的四步反应分别为:气相原子的化学吸附、一个气相原子和一个吸附原子的ER复合、两个吸附原子的LH复合和吸附原子的热解附。所述求解方法包括如下步骤:a:分析当前的流场组分环境,得到气相反应物的浓度;b:根据固体表面参数得到各步反应的速率常数;c:基于速率常数和流场环境获得各步的反应速率;d:以稳定状态下吸附相浓度不随时间改变为约束方程,求解真实的吸附相浓度;e:根据反应物浓度和速率常数获得各步真实反应速率,并得到各反应步共同导致的各组分净反应速率。
根据一个优选的实施方式,,所述步骤a包括:根据非平衡流场求解器获得当前流场,并获得流场网格中壁面法向第一层网格上的气体A的组分浓度cA
根据一个优选的实施方式,所述步骤b包括:令固体表面的吸附位个数为S0,根据固体表面温度Tw、各步反应的概率系数P,活化能E、表面原子的微观运动速率v、吸附位间的平均距离dS获得各步反应的反应速率常数k:
其中,R为通用气体常数,h为普朗克常数,下标Ad、ER、LH、Td分别代表吸附、ER复合、LH复合和热解附四种反应类型;且气相原子的三维平均热运动速度v3D和固体表面原子二维平均热速度v2D分别为:
其中,MA为原子质量。
根据一个优选的实施方式,所述步骤c包括:基于步骤a和步骤b获得的反应物浓度cA和速率常数,针对原子组分A发生的四步反应过程,其中吸附反应生成吸附相原子A(*)、ER/LH复合、热解附消耗吸附相原子A(*),其中Av为阿伏伽德罗常数,其对应的各步反应速率为,RAd=cA×(S0/Av-cA(*))×kAd
RER=cA×cA(*)×kER
RLH=cA(*)×cA(*)×kLH
RTd=cA(*)×kTd
根据一个优选的实施方式,所述步骤d包括:稳定状态下,表面吸附相的浓度不随时间而改变,吸附反应生成吸附相原子A(*)的速率与ER、LH复合和热解附消耗A(*)的速率达到平衡,即是:
从而,基于该约束方程,完成吸附相反应物A(*)的浓度cA(*)的求解。
根据一个优选的实施方式,所述步骤e包括:基于各步反应的速率,获得当前该原子及其相应的双原子分子的净反应速率;
R(A)=-RER-2×RLH
前述本发明主方案及其各进一步选择方案可以自由组合以形成多个方案,均为本发明可采用并要求保护的方案;且本发明,各非冲突选择之间也可以自由组合。本领域技术人员在了解本发明方案后根据现有技术和公知常识可明了有多种组合,均为本发明所要保护的技术方案,在此不做穷举。
本发明的有益效果:本发明方法以真实的多步催化反应过程为理论基础,并根据真实的宏微观表面参数分别确定每一步反应的理论速率,最后考虑各步反应之间的约束和影响,得到真实的各步反应速率。该方法依赖于真实的表面物理化学过程和表面参数,因此更具确定性,同时具有材料依赖性。利用该催化反应速率求解方法作为边界条件计算得到的不同材料表面的气动热Q与实验的对比,证明了该方法的准确性和有效性。
附图说明
图1是本发明方法进行多相催化反应中组分净反应速率的求解流程;
图2是本发明方法涉及的多相催化反应的过程示意图;
图3是本发明方法的气动热求解精度验证。
具体实施方式
以下通过特定的具体实例说明本发明的实施方式,本领域技术人员可由本说明书所揭露的内容轻易地了解本发明的其他优点与功效。本发明还可以通过另外不同的具体实施方式加以实施或应用,本说明书中的各项细节也可以基于不同观点与应用,在没有背离本发明的精神下进行各种修饰或改变。需说明的是,在不冲突的情况下,以下实施例及实施例中的特征可以相互组合。
需要说明的是,为使本发明实施例的目的、技术方案和优点更加清楚,下面对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。
实施例1:
参考图1至图3所示,本发明公开了一种用于CFD的高焓离解气体表面催化反应速率的求解方法。
本发明的求解方法是根据化学反应动力学理论,基于固体表面存在有限个数的吸附位点的理论假设,求解在此吸附位点上发生的四步反应:气相原子的化学吸附、一个气相原子和一个吸附原子的Eley-Rideal(ER)复合、两个吸附原子的Langmuir-Hinshelwood(LH)复合和吸附原子的热解附。
每一步反应消耗或生成原子的速率由各步反应的速率常数和反应物浓度共同决定,各步反应速率常数考虑了表面温度、反应活化能等物理参数的影响,气相反应物的浓度由非平衡流场求解器获得,表面吸附相反应物的浓度由稳态下表面吸附位的覆盖率不随时间改变求解得到。
本发明获得的反应速率将用以解决化学非平衡流场求解过程遇到的催化壁面表征困难、边界质量能量输运不准确的问题。
参考图1所示,本发明求解界面上由催化反应导致的各组分净反应速率的方法流程主要包括以下步骤:
a.根据非平衡流场CFD求解器获得当前流场,提取流场网格中壁面法向第一层网格上的气体A的组分浓度cA(单位mol/m3)。将该气体组分作为多相催化求解方法的输入条件之一。
b.假定固体表面的吸附位个数为S0(单位1/m2),根据固体表面温度TW、活化能E、各步反应的概率系数P,气相原子的三维平均热速度v3D、表面原子的二维平均热速度v2D、吸附位间的平均距离dS等获得各步反应的反应速率常数kAd、kER、kLH、kTd。其中下标Ad、ER、LH、Td分别代表吸附、ER复合、LH复合和热解附四种反应类型:
h为普朗克常数,其中原子的三维和二维平均热速度v3D和v2D通过固体表面温度、通用气体常数R和原子的质量MA计算:
c.针对原子组分A发生的四类反应过程(如图2),基于a和b获得的气相反应物浓度cA、计算过程中得到吸附相原子浓度cA(*)(单位mol/m3)和速率常数,可以得到四类反应的速率,见表1,其中Av为阿伏伽德罗常数。四种反应类型中,吸附反应生成吸附相原子A(*),ER/LH复合、热解附消耗A(*),吸附相反应物的浓度cA(*)作为待解变量;
表1四种反应类型的理论反应速率
d.吸附相反应物浓度cA(*)的求解:稳定状态下,表面吸附相的浓度应不随时间而改变,因此,吸附反应生成吸附相原子A(*)的速率与ER、LH复合和热解附消耗A(*)的速率应该达到平衡(式(7)),该平衡方程为关于cA(*)的一元二次方程,求解即可获得当前吸附相反应物的浓度cA(*)
e.至此,该迭代步下各步反应的速率常数和气相、吸附相反应物浓度均已获得,求得各步反应的速率后,四步反应分别导致的气相原子A变化速率叠加,即获得当前该原子及其相应的双原子分子的净反应速率R(A)和R(A2);
R(A)=-RAd-RER-2×RLH+RTd (8)
f.根据个流动边界上的质量和能量守恒,多相催化反应对流场控制方程组所起的传热和传质作用均可用e步中获得的组分A的净反应速率R(A)表示,因此该方法所求得的净反应速率可以直接用于CFD求解;
本发明方法以真实的多步催化反应过程为理论基础,并根据真实的宏微观表面参数分别确定每一步反应的理论速率,最后考虑各步反应之间的约束和影响,得到真实的各步反应速率。该方法依赖于真实的表面物理化学过程和表面参数,因此更具确定性,同时具有材料依赖性。利用该催化反应速率求解方法作为边界条件计算得到的不同材料表面的气动热Q与实验的对比如图3所示,证明了该方法的准确性和有效性。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。

Claims (1)

1.一种用于CFD的高焓离解气体表面催化反应速率的求解方法,其特征在于,所述求解方法包括基于固体表面存在有限个数的吸附位点的理论,求解气相反应物在固体表面吸附位点上发生的四步反应的反应速率;
其中,气相反应物在固体表面吸附位点上发生的四步反应分别为:气相原子的化学吸附、一个气相原子和一个吸附原子的ER复合、两个吸附原子的LH复合和吸附原子的热解附;
所述求解方法包括如下步骤:
a:分析当前的流场组分环境,得到气相反应物的浓度;
b:根据固体表面参数得到各步反应的速率常数;
c:基于速率常数和流场环境获得各步的反应速率;
d:以稳定状态下吸附相浓度不随时间改变为约束方程,求解真实的吸附相浓度;
e:根据反应物浓度和速率常数获得各步真实反应速率,并得到各反应步共同导致的各组分净反应速率;
所述步骤a包括:根据非平衡流场求解器获得当前流场,并获得流场网格中壁面法向第一层网格上的气体A的组分浓度cA
所述步骤b包括:
令固体表面的吸附位个数为S0,根据固体表面温度Tw、各步反应的概率系数P、活化能E、表面原子的微观运动速率v、吸附位间的平均距离dS获得各步反应的反应速率常数k:
其中,R为通用气体常数,h为普朗克常数,下标Ad、ER、LH、Td分别代表吸附、ER复合、LH复合和热解附四种反应类型;且气相原子的三维平均热运动速度v3D和固体表面原子二维平均热速度v2D分别为:
其中,MA为原子质量;
所述步骤c包括:
基于步骤a和步骤b获得的反应物浓度cA和速率常数,针对原子组分A发生的四步反应过程,其中吸附反应生成吸附相原子A(*)、ER/LH复合、热解附消耗吸附相原子A(*),则吸附、ER复合、LH复合和热解附四步的反应速率分别为:
RAd=cA×(S0/Av-cA(*))×kAd
RER=cA×cA(*)×kER
RLH=cA(*)×cA(*)×kLH
RTd=cA(*)×kTd
其中Av为阿伏伽德罗常数;
所述步骤d包括:
稳定状态下,表面吸附相的浓度不随时间而改变,吸附反应生成吸附相原子A(*)的速率与ER、LH复合和热解附消耗A(*)的速率达到平衡,即是:
从而,基于该约束方程,完成吸附相反应物A(*)的浓度cA(*)的求解;
所述步骤e包括:
基于各步反应的速率,获得当前该原子及其相应的双原子分子的净反应速率;
R(A)=-RER-2×RLH
CN202110862009.3A 2021-07-29 2021-07-29 一种用于cfd的高焓离解气体表面催化反应速率的求解方法 Active CN113782106B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110862009.3A CN113782106B (zh) 2021-07-29 2021-07-29 一种用于cfd的高焓离解气体表面催化反应速率的求解方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110862009.3A CN113782106B (zh) 2021-07-29 2021-07-29 一种用于cfd的高焓离解气体表面催化反应速率的求解方法

Publications (2)

Publication Number Publication Date
CN113782106A CN113782106A (zh) 2021-12-10
CN113782106B true CN113782106B (zh) 2024-04-16

Family

ID=78836578

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110862009.3A Active CN113782106B (zh) 2021-07-29 2021-07-29 一种用于cfd的高焓离解气体表面催化反应速率的求解方法

Country Status (1)

Country Link
CN (1) CN113782106B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2865820A1 (fr) * 2004-01-30 2005-08-05 Novalyst Discovery Procede pour constituer une base de donnees permettant de selectionner au moins un catalyseur adapte a une reaction
CN108573124A (zh) * 2018-05-15 2018-09-25 中国石油大学(华东) 一种定量分析金属嵌入碳纳米管催化剂析氢活性的方法
KR20200001725A (ko) * 2018-06-28 2020-01-07 울산과학기술원 고효율 물산화 촉매, 이의 제조방법 및 광전극
CN111859532A (zh) * 2020-06-16 2020-10-30 空气动力学国家重点实验室 考虑高超声速化学非平衡效应的改进热壁修正方法
CN112800643A (zh) * 2020-12-30 2021-05-14 新源动力股份有限公司 一种波纹流道燃料电池多物理场耦合计算简化方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050251346A1 (en) * 2004-03-29 2005-11-10 Ilie Fishtik Method and apparatus for reaction route graphs for reaction mechanism and kinetics modeling

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2865820A1 (fr) * 2004-01-30 2005-08-05 Novalyst Discovery Procede pour constituer une base de donnees permettant de selectionner au moins un catalyseur adapte a une reaction
CN108573124A (zh) * 2018-05-15 2018-09-25 中国石油大学(华东) 一种定量分析金属嵌入碳纳米管催化剂析氢活性的方法
KR20200001725A (ko) * 2018-06-28 2020-01-07 울산과학기술원 고효율 물산화 촉매, 이의 제조방법 및 광전극
CN111859532A (zh) * 2020-06-16 2020-10-30 空气动力学国家重点实验室 考虑高超声速化学非平衡效应的改进热壁修正方法
CN112800643A (zh) * 2020-12-30 2021-05-14 新源动力股份有限公司 一种波纹流道燃料电池多物理场耦合计算简化方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Alexandra N. Molchanova et al..《PHYSICS OF FLUIDS》 Surface recombination in the direct simulation Monte Carlo method.2018,第30卷第1-18页. *
M Cacciatore et al..《Plasma Sources Sci. Technol.》Dynamics of plasma–surface processes: E–R and L–H atom recombination reactions.2009,第18卷第1-13页. *
杨肖峰 等.《国防科技大学学报 》高焓CO2气流壁面两步催化机制对非平衡气动加热影响的数值模拟.2020,第42卷(第1期),第108-116页. *
梁辉 等.《东南大学学报》纳米γFe2O3催化剂SCR脱除NOx化学反应动力学.2015,第45卷(第5期),第918-922页. *

Also Published As

Publication number Publication date
CN113782106A (zh) 2021-12-10

Similar Documents

Publication Publication Date Title
Nielsen et al. Parallel kinetic Monte Carlo simulation framework incorporating accurate models of adsorbate lateral interactions
Lele et al. Thermal decomposition kinetic of salt hydrates for heat storage systems
Wheeler et al. Modelling of solar thermochemical reaction systems
Zhang et al. Multiscale computational fluid dynamics modeling of thermal atomic layer deposition with application to chamber design
Gakis et al. Detailed investigation of the surface mechanisms and their interplay with transport phenomena in alumina atomic layer deposition from TMA and water
CN111859532B (zh) 考虑高超声速化学非平衡效应的改进热壁修正方法
Partopour et al. Resolved‐particle fixed bed CFD with microkinetics for ethylene oxidation
Ostermeier et al. Numerical calculation of wall-to-bed heat transfer coefficients in Geldart B bubbling fluidized beds with immersed horizontal tubes
Pereira et al. Multi-scale modeling of diffusion and reaction–diffusion phenomena in catalytic porous layers: comparison with the 1D approach
CN111597735A (zh) 机器学习与cvd建模相结合的组分预测方法
Wuerth et al. Development of a continuous fluidized bed reactor for thermochemical energy storage application
Wu et al. Development of methanol steam reforming microreactor based on stacked wave sheets and copper foam for hydrogen production
Coetzee et al. The mechanistic effect over the substrate in a square type atomic layer deposition reactor
CN113782106B (zh) 一种用于cfd的高焓离解气体表面催化反应速率的求解方法
Okoli et al. Application of an equation‐oriented framework to formulate and estimate parameters of chemical looping reaction models
Coetzee et al. The fluid flow effect on the inlet injection of the thin film deposition in a square type atomic layer deposition reactor
Vanka et al. Numerical study of mixed convection flow in an impinging jet CVD reactor for atmospheric pressure deposition of thin films
Nwanna et al. A numerical approach on the selection of the purge flow rate in an atomic layer deposition (ALD) process
De Falco Membrane reactors modeling
Singhal et al. Multiscale modelling of packed bed chemical looping reforming
Othman et al. 3D CFD simulations of MOCVD synthesis system of titanium dioxide nanoparticles
Chutichairattanaphum et al. Effects of Raschig ring packing patterns on pressure drop, heat transfer, methane conversion, and coke deposition on a semi-pilot-scale packed bed reformer
Teixeira et al. A computational fluid dynamics-based sensitivity analysis of the chemical vapor analysis process to synthesize carbon nanotubes
George et al. Optimization of gallium nitride metalorganic chemical vapor deposition process
Li et al. The influences of oxygen concentration and external heating on carbon nanotube growth in diffusion flame

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
GR01 Patent grant
GR01 Patent grant