CN111859531B - 一种求解化学非平衡流动化学反应源项的预估修正方法 - Google Patents

一种求解化学非平衡流动化学反应源项的预估修正方法 Download PDF

Info

Publication number
CN111859531B
CN111859531B CN202010540454.3A CN202010540454A CN111859531B CN 111859531 B CN111859531 B CN 111859531B CN 202010540454 A CN202010540454 A CN 202010540454A CN 111859531 B CN111859531 B CN 111859531B
Authority
CN
China
Prior art keywords
round
chemical reaction
sub
moment
component density
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
CN202010540454.3A
Other languages
English (en)
Other versions
CN111859531A (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 CN202010540454.3A priority Critical patent/CN111859531B/zh
Publication of CN111859531A publication Critical patent/CN111859531A/zh
Application granted granted Critical
Publication of CN111859531B publication Critical patent/CN111859531B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/15—Vehicle, aircraft or watercraft design
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
    • G16C20/10—Analysis or design of chemical reactions, syntheses or processes
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Mathematical Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Optimization (AREA)
  • Chemical & Material Sciences (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Computing Systems (AREA)
  • Analytical Chemistry (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Mathematical Physics (AREA)
  • Fluid Mechanics (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Algebra (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Automation & Control Theory (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Computational Mathematics (AREA)
  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明涉及流体力学数值计算技术领域，公开了一种求解化学非平衡流动化学反应源项的预估修正方法，包括以下步骤：a.根据当前时刻的流场解获得当前时刻的组分密度和化学反应生成率，预设该计算轮次N＝1；b.开展下一轮预估，通过二分法剖分时间离散间隔；依次对二分点的参数值进行预估，并对上轮各点的参数值进行修正；c.计算第N+1轮和第N轮的预估值之差，若相对值较大，则增加一轮预估修正，即N＝N+1，进行步骤b操作，否则跳转至下一步骤；d.将第N+1轮的计算结果作为最终结果。本方法采用核心控制方程隐式处理和化学反应源项半隐式处理，在保持原有LUSGS算法的优势的同时，解决了化学反应源项的刚性问题。

Description

一种求解化学非平衡流动化学反应源项的预估修正方法
技术领域
本发明涉及流体力学数值计算技术领域，尤其涉及一种求解化学非平衡流动化学反应源项的预估修正方法。
背景技术
在利用计算流体力学技术求解高超声速化学非平衡流动研究中，需要解决因化学反应特征时间尺度与流动不匹配而引起的刚性问题，其中化学反应生成率及各组分密度是所关注的重要参数。时间推进的LUSGS隐式算法是可压缩流动计算的常规方法，而将化学反应源项作隐式处理的全隐式算法虽然有较高的稳定性和收敛性，但会造成系数矩阵的非对角化，破坏原有算法的优越性。因此，针对高超声速化学非平衡流动的数值计算，需要发展具有较高收敛性和计算效率的求解化学非平衡流动化学反应源项的方法。
发明内容
本发明所要解决的技术问题是：针对上述存在的问题，提供了一种求解化学非平衡流动化学反应源项的预估修正方法，本方法借鉴流动与化学反应解耦的物理分裂思想，假设化学非平衡过程是由系列局部等压非平衡子过程组成，冻结组分连续性方程的对流和扩散项，分层次对时间步作子过程剖分，并在不同层次对预估值进行修正，进而获得下一时刻的化学反应源项。
本发明采用的技术方案如下：一种求解化学非平衡流动化学反应源项的预估修正方法，包括以下步骤：
步骤1：根据当前时刻的流场，解获得到当前时刻的参数值，设定该计算轮次为初始轮次；
步骤2：根据当前时刻已知的参数值，开展下一轮次的预估，先通过二分法剖分时间离散间隔，然后依次对二分点的参数值进行预估，并对上一轮次各二分点的参数值进行修正；
步骤3：计算第N+1轮和第N轮的预估值之差，若差值大于设定的阈值，则返回步骤2，增加一轮预估修正；否则跳转至下一步骤，其中N为大于等于1的自然数；
步骤4：将第N+1轮计算得到的参数值作为最终结果。
这样，通过综合采用核心控制方程隐式处理和化学反应源项半隐式处理的方式，在保持原有算法优势的情形下，同时解决了化学反应源项的刚性问题，且程序实现简单。
进一步的，所述参数值包括组分密度和化学反应生成率。
进一步的，所述步骤2具体包括以下步骤：
步骤21：根据当前时刻已知的组分密度和化学反应生成率，计算下一时刻的组分密度；
步骤22：进入下一轮次预估修正；
步骤23：对上一轮次2N-1个子过程作二分，形成2N个子过程；
步骤24：设定子过程时刻i＝1，开始按子过程时刻进行推进；
步骤25：求解第N+1轮二分点时刻2i-1的组分密度和化学反应生成率，其中，组分密度根据第N+1轮上一子时刻2i-2和第N轮对应下一子时刻i的相应值共同求得；对应的化学反应生成率由化学反应动力学获得；
步骤26：求解第N+1轮整点2i时刻的组分密度和化学反应生成率，其中，组分密度根据第N+1轮上一子时刻2i-1和第N轮子时刻i的相应值共同求得，对应的化学反应生成率由化学反应动力学获得；
步骤27：若当前子过程为该轮次最后子过程，即i＝2N-1，则进入下一步；否则子过程时刻继续推进，i＝i+1，并返回步骤25；
步骤28：计算完当前轮次中的所有子过程，获得基于本轮次预估修正的下一时刻组分密度和化学反应生成率。
进一步的，所述步骤25中，第N+1轮二分点位置的组分密度计算公式如下：
其中，表示第N+1轮前2i-2时刻的组分密度，/>表示第(N+1)轮前2i-2时刻的化学反应生成率，/>表示第N轮对应下一子时刻i的化学反应生成率。
进一步的，所述步骤26中，第N+1轮整点位置的组分密度计算公式如下：
其中，表示第N+1轮上一子时刻2i-1的组分密度，/>表示第N+1轮上一子时刻2i-1的化学反应生成率，/>表示第N轮对应下一子时刻i的化学反应生成率。
与现有技术相比，采用上述技术方案的有益效果为：本发明的求解化学反应源项的预估修正方法，对核心控制方程隐式处理、化学反应源项半隐式处理，保持原有LUSGS算法的优势，同时高效准确地求解化学反应源项，程序实现简单。
附图说明
图1是本发明的求解化学反应源项的预估修正实施示意图。
图2是本发明的求解化学非平衡流动化学反应源项的预估修正方法的流程图。
具体实施方式
本发明提供一种求解化学非平衡流动化学反应源项的预估修正方法。本发明的求解化学非平衡流动化学反应源项的预估修正方法借鉴流动与化学反应解耦的物理分裂思想，假设化学非平衡过程是由系列局部等压非平衡子过程组成，冻结组分连续性方程的对流和扩散项，分层次对时间步作子过程剖分，并在不同层次对预估值进行修正，进而获得下一时刻的化学反应源项。
下面结合附图进一步对本发明的求解化学非平衡流动化学反应源项的预估修正方法做详细说明。
如图1所示，在已知当前时刻t0的流场，借鉴流动与化学反应解耦的物理分裂思想，假设化学非平衡过程是由系列局部等压非平衡子过程组成，并按照二分法对时间间隔Δt做均匀剖分，第1轮生成1个子过程(21-1)，第2轮生成2个子过程(22-1)，……，第N轮生成2N-1个子过程。按照预估修正规则对各轮次的子过程作计算。
对照如图2所示的流程图给出具体实施步骤。
步骤01：基于当前时刻流场准备求解化学反应源项所需的变量。
步骤02：获取当前时刻t0的组分密度ρ0和化学反应生成率ω0。
步骤03：设定初始轮次N＝1。
步骤04：计算下一时刻t0+Δt的组分密度ρ1 (1)，即
对应的化学反应生成率ω1 (1)由化学反应动力学获得。
步骤05：在已知第N轮各子过程的解情况下，可对第N+1轮各子过程进行预估修正。所以，进入下一轮预估修正，N＝N+1。
步骤06：对上一轮次2N-1个Δt子过程作二分，形成2N个Δt子过程。
步骤07：按子过程时刻顺序开始推进，设定子过程初始时刻i＝1。
步骤08：已获得第(N+1)轮前(2i-2)时刻的组分密度和化学反应生成率求解第N+1轮二分点时刻(2i-1)的组分密度和化学反应生成率，其中组分密度根据第(N+1)轮上一子时刻(2i-2)和第N轮对应下一子时刻(i)的相应值共同求得，如公式(1)
对应的化学反应生成率由化学反应动力学获得。
步骤09：对于第N+1轮整点2i时刻，其组分密度根据第N+1轮上一子时刻(2i-1)和第N轮该子时刻(i)的相应值共同求得，如公式(2)
对应的化学反应生成率由化学反应动力学获得。
步骤10：进行判断
若当前子过程为该轮次最后子过程，即i＝2N-1，则进入步骤11；否则，i＝i+1，并返回步骤08计算下一子过程。
步骤11：计算完成当前轮次中所有子过程后，获得基于本轮次预估修正的下一时刻组分密度和化学反应生成率。
步骤12：进行判断
若第N+1轮和第N轮的化学反应生成率之差足够小，则进入步骤13；否则，N＝N+1，增加一轮预估修正。
步骤13：最终获得下一时刻的组分密度和化学反应生成率，将其作为最终计算结果，结束化学反应源项的预估修正计算。
本发明的求解化学非平衡流动化学反应源项的预估修正方法综合采用核心控制方程隐式处理、化学反应源项半隐式处理，保持原有算法的优势，同时解决化学反应源项的刚性问题，程序实现简单。本发明的求解化学非平衡流动化学反应源项的预估修正方法可快速高效求解化学非平衡流动问题，具有较好的计算精度和效率。
本发明并不局限于前述的具体实施方式。本发明扩展到任何在本说明书中披露的新特征或任何新的组合，以及披露的任一新的方法或过程的步骤或任何新的组合。如果本领域技术人员，在不脱离本发明的精神所做的非实质性改变或改进，都应该属于本发明权利要求保护的范围。

Claims (1)

1.一种求解化学非平衡流动化学反应源项的预估修正方法，其特征在于，包括以下步骤：
步骤1：根据当前时刻的流场，计算得到当前时刻的参数值，设定该计算轮次为初始轮次，所述参数值包括组分密度和化学反应生成率；
步骤21：根据当前时刻已知的组分密度和化学反应生成率，计算下一时刻的组分密度；
步骤22：进入下一轮次预估修正；
步骤23：对上一轮次2N-1个子过程作二分，形成2N个子过程；
步骤24：设定子过程时刻i＝1，开始按子过程时刻进行推进；
步骤25：求解第N+1轮二分点时刻2i-1的组分密度和化学反应生成率，其中，组分密度根据第N+1轮上一子时刻2i-2和第N轮对应下一子时刻i的相应值共同求得；对应的化学反应生成率由化学反应动力学获得；
步骤26：求解第N+1轮整点2i时刻的组分密度和化学反应生成率，其中，组分密度根据第N+1轮二分点时刻2i-1和第N轮对应下一子时刻i的相应值共同求得，对应的化学反应生成率由化学反应动力学获得；
步骤27：若当前子过程为该轮次最后子过程，即i＝2N-1，则进入下一步；否则子过程时刻继续推进，i＝i+1，并返回步骤25；
步骤28：计算完当前轮次中的所有子过程，获得基于本轮次预估修正的下一时刻组分密度和化学反应生成率；
所述步骤25中，第N+1轮二分点位置的组分密度计算公式如下：
其中，表示第N+1轮上一子时刻2i-2的组分密度，/>表示第N+1轮上一子时刻2i-2的化学反应生成率，/>表示第N轮对应下一子时刻i的化学反应生成率；
所述步骤26中，第N+1轮整点位置的组分密度计算公式如下：
其中，表示第N+1轮二分点位置2i-1的组分密度，/>表示第N+1轮上一子时刻2i-1的化学反应生成率，/>表示第N轮对应下一子时刻i的化学反应生成率；
步骤3：计算第N+1轮和第N轮的预估值之差，若差值大于设定的阈值，则返回步骤2，增加一轮预估修正；否则跳转至下一步骤，其中N为大于等于1的自然数；
步骤4：将第N+1轮计算得到的参数值作为最终结果。
CN202010540454.3A 2020-06-15 2020-06-15 一种求解化学非平衡流动化学反应源项的预估修正方法 Active CN111859531B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010540454.3A CN111859531B (zh) 2020-06-15 2020-06-15 一种求解化学非平衡流动化学反应源项的预估修正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010540454.3A CN111859531B (zh) 2020-06-15 2020-06-15 一种求解化学非平衡流动化学反应源项的预估修正方法

Publications (2)

Publication Number Publication Date
CN111859531A CN111859531A (zh) 2020-10-30
CN111859531B true CN111859531B (zh) 2023-10-13

Family

ID=72987405

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010540454.3A Active CN111859531B (zh) 2020-06-15 2020-06-15 一种求解化学非平衡流动化学反应源项的预估修正方法

Country Status (1)

Country Link
CN (1) CN111859531B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116306384B (zh) * 2023-05-23 2023-08-25 中国空气动力研究与发展中心计算空气动力研究所 基于热化学非平衡源项约束的流场获取方法及其相关组件
CN117393062B (zh) * 2023-12-13 2024-02-23 上海交通大学四川研究院 刚性化学反应流回退自适应半隐半显耦合时间的模拟方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107368661A (zh) * 2017-07-27 2017-11-21 中国空气动力研究与发展中心计算空气动力研究所 一种高超声速飞行器热气动弹性特性的耦合计算方法
CN108520144A (zh) * 2018-04-08 2018-09-11 中国人民解放军国防科技大学 基于压力的特征线法中无激波压力最速上升曲线的求解方法
CN111222241A (zh) * 2020-01-06 2020-06-02 中国人民解放军国防科技大学 热化学非平衡条件下流场数据的数值计算方法和装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3182371B1 (en) * 2015-12-17 2018-09-26 Stmicroelectronics Sa Threshold determination in for example a type ransac algorithm

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107368661A (zh) * 2017-07-27 2017-11-21 中国空气动力研究与发展中心计算空气动力研究所 一种高超声速飞行器热气动弹性特性的耦合计算方法
CN108520144A (zh) * 2018-04-08 2018-09-11 中国人民解放军国防科技大学 基于压力的特征线法中无激波压力最速上升曲线的求解方法
CN111222241A (zh) * 2020-01-06 2020-06-02 中国人民解放军国防科技大学 热化学非平衡条件下流场数据的数值计算方法和装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Z. Nussinov 等.Emergent quantization and possible measurements of Planck's constant from the saturation of extrapolated high temperature liquid viscosity and other transport properties.《arXiv》.2015,1-27. *
陈博 等.离子流场模型中高压直流导线表面粗糙系数的表征方法.《高电压技术》.2019,第45卷(第5期),1443-1450. *

Also Published As

Publication number Publication date
CN111859531A (zh) 2020-10-30

Similar Documents

Publication Publication Date Title
CN111859531B (zh) 一种求解化学非平衡流动化学反应源项的预估修正方法
CN107644630B (zh) 基于神经网络的旋律生成方法、装置及存储介质
CN105300384B (zh) 一种用于卫星姿态确定的交互式滤波方法
CN111260201B (zh) 一种基于分层随机森林的变量重要性分析方法
CN103439887B (zh) 低阶系统itae最优的pi控制器参数整定方法及系统
CN109447133B (zh) 一种基于svr算法的位置信息野值的剔除方法
CN111027732A (zh) 一种多风电场出力场景的生成方法及系统
CN111679643B (zh) 基于过程迁移模型的间歇过程二阶修正自适应优化方法
CN112001596B (zh) 一种时间序列数据异常点检测方法及系统
CN108803528A (zh) 基于多变量相关性及时滞性的流程工业系统预测模型
CN111176810B (zh) 基于优先级的气象水文数据处理调度系统
JP2020102082A (ja) 予測モデル生成装置、予測モデル生成方法、および、コンピュータプログラム
CN107871180A (zh) 基于牛顿‑柯特斯公式构造背景值的gm（1，1）模型预测方法
JP2007173334A (ja) 半導体熱処理装置
CN111241749B (zh) 一种基于储备池计算的永磁同步电动机混沌预测方法
CN112985688B (zh) 压力传感器的输出校准方法、装置及可读存储介质
CN112631214B (zh) 基于改进入侵杂草优化算法的柔性作业车间分批调度方法
CN108388218A (zh) 基于潜变量过程迁移模型的修正自适应批次过程优化方法
Villaverde et al. High-confidence predictions in systems biology dynamic models
Fourment et al. Mono and multi‐objective optimization techniques applied to a large range of industrial test cases using Metamodel assisted Evolutionary Algorithms
CN108197809B (zh) 基于维度最优转换的权值共享深度网络的实时调度方法
CN114819398A (zh) 基于灰色马尔可夫链的北斗卫星钟差序列组合预测方法
CN109884893B (zh) 一种多工序变量间动态时滞估计方法
CN111950123A (zh) 一种陀螺仪误差系数曲线拟合预测方法及系统
CN113868887A (zh) 水洗塔模拟计算方法

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