CN113177374A - 一种燃烧特性获取方法、装置、电子设备及存储介质 - Google Patents
一种燃烧特性获取方法、装置、电子设备及存储介质 Download PDFInfo
- Publication number
- CN113177374A CN113177374A CN202110468022.0A CN202110468022A CN113177374A CN 113177374 A CN113177374 A CN 113177374A CN 202110468022 A CN202110468022 A CN 202110468022A CN 113177374 A CN113177374 A CN 113177374A
- Authority
- CN
- China
- Prior art keywords
- combustion
- flow field
- model
- simulation calculation
- target solid
- 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
Links
Images
Classifications
-
- 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
- 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/08—Thermal analysis or thermal optimisation
-
- 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)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Fluid Mechanics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computing Systems (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及燃烧特性获取方法、装置、电子设备及存储介质,燃烧特性获取方法包括:对与目标固体模型对应的燃烧模型进行网格划分,对网格划分后的所述燃烧模型的燃烧流场进行燃烧模拟计算,得到模拟计算结果,所述燃烧流场为燃料燃烧产生的流场,将所述模拟计算结果与所述目标固体模型关联,对关联后的目标固体模型进行网格划分,对网格划分后的所述目标固体模型在所述燃烧流场作用下的燃烧特性进行模拟计算,得到每个网格内固体结构的位移变化量和等效应力变化量。本发明实施例可以预测目标固体模型对应的零件内部在燃料燃烧时,零件的燃烧特性,进而对产品设计提供了一定的理论依据,更加有助于后期产品设计优化,提高产品性能。
Description
技术领域
本申请涉及计算机技术领域,尤其涉及一种燃烧特性获取方法、装置、电子设备及存储介质。
背景技术
在产品(如:发动机气缸等)设计中,为了准确预测产品结构的合理性,一般需要经过大量的实验,以了解产品内部流场中的热力对于固体结构产生的应力以及应变,进而得到产品在结构上需要改进的方向等。
目前的实验方法是首先生产出实体的产品,对产品施加实际的环境因素,如:向产品内部腔体注入混合气体,并点燃注入腔体的混合气体等,然而,这种产品设计方式,效率非常低,且造成资源的浪费。
发明内容
为了解决上述产品设计需要生产实体产品,再进行试验,效率非常低,且造成资源的浪费的技术问题,本申请提供了一种燃烧特性获取方法、装置、电子设备及存储介质。
第一方面,本申请提供了一种燃烧特性获取方法,包括:
对与目标固体模型对应的燃烧模型进行网格划分;
对网格划分后的所述燃烧模型的燃烧流场进行燃烧模拟计算,得到模拟计算结果,所述燃烧流场为燃料燃烧产生的流场;
将所述模拟计算结果与所述目标固体模型关联,对关联后的目标固体模型进行网格划分;
对网格划分后的所述目标固体模型在所述燃烧流场作用下的燃烧特性进行模拟计算,得到每个网格内固体结构的位移变化量和等效应力变化量。
可选地,对网格划分后的所述燃烧模型的燃烧流场进行燃烧模拟计算,得到模拟计算结果,包括:
设置所述燃烧模型的边界条件、松弛因子和收敛精度;
对网格划分后的所述燃烧模型的冷态流场进行模拟计算,所述冷态流场为燃料与空气混合且燃料未燃烧放热的流场;
基于模拟计算后的所述冷态流场,对网格划分后的所述燃烧模型的燃烧流场进行模拟计算,得到模拟计算结果。
可选地,所述基于模拟计算后的所述冷态流场,对网格划分后的所述燃烧模型的燃烧流场进行模拟计算,得到模拟计算结果,包括:
设置基于模拟计算后的所述冷态流场的初始温度及内部气体浓度;
将已设置初始温度及内部气体浓度的所述冷态流场作为所述燃烧流场的初始流场,模拟计算所述燃烧模型的燃烧流场。
可选地,将所述模拟计算结果与所述目标固体模型关联,对关联后的目标固体模型进行网格划分,包括:
获取所述目标固体模型;
将所述模拟计算结果导入所述目标固体模型;
将已导入所述模拟计算结果的所述目标固体模型进行网格划分。
可选地,对网格划分后的所述目标固体模型在所述燃烧流场作用下的燃烧特性进行模拟计算,得到每个网格内固体结构的位移变化量和等效应力变化量,包括:
设置所述目标固体模型的边界条件;
设置位移求解器及等效应力求解器;
对所述目标固体模型的燃烧特性进行模拟计算,得到每个网格内固体结构的位移变化量和等效应力变化量。
可选地,对与目标固体模型对应的燃烧模型进行网格划分,包括:
使用Gambit网格划分软件采用四面体网格对所述燃烧模型进行网格划分,得到网格文件;
将网格文件导入FLUENT软件。
可选地,在对与目标固体模型对应的燃烧模型进行网格划分之后,所述方法还包括:
检查所述网格文件中网格的质量,将所述网格文件中网格的长度单位修改为预设单位;
设置所述燃烧模型的反应类型;
将所述燃烧模型设置为涡耗散模型;
设置所述燃烧模型考虑湍流流动中的化学反应和混合比率。
第二方面,本申请提供了一种燃烧特性获取装置,包括:
网格划分模块,用于对与目标固体模型对应的燃烧模型进行网格划分;
第一模拟计算模块,用于对网格划分后的所述燃烧模型的燃烧流场进行燃烧模拟计算,得到模拟计算结果,所述燃烧流场为燃料燃烧产生的流场;
关联划分模块,用于将所述模拟计算结果与所述目标固体模型关联,对关联后的目标固体模型进行网格划分;
第二模拟计算模块,用于对网格划分后的所述目标固体模型在所述燃烧流场作用下的燃烧特性进行模拟计算,得到每个网格内固体结构的位移变化量和等效应力变化量。
第三方面,本申请提供了一种电子设备,包括处理器、通信接口、存储器和通信总线,其中,处理器,通信接口,存储器通过通信总线完成相互间的通信;
存储器,用于存放计算机程序;
处理器,用于执行存储器上所存放的程序时,实现第一方面任一所述的燃烧特性获取方法。
第四方面,本申请提供了一种计算机可读存储介质,所述计算机可读存储介质上存储有燃烧特性获取方法的程序,所述燃烧特性获取方法的程序被处理器执行时实现第一方面任一所述的燃烧特性获取方法的步骤。
本申请实施例提供的上述技术方案与现有技术相比具有如下优点:
本申请实施例提供的该方法,通过首先对与目标固体模型对应的燃烧模型进行网格划分,然后对网格划分后的所述燃烧模型的燃烧流场进行燃烧模拟计算,得到模拟计算结果,燃烧流场为燃料燃烧产生的流场,再将所述模拟计算结果与所述目标固体模型关联,对关联后的目标固体模型进行网格划分,最后可以对网格划分后的所述目标固体模型在所述燃烧流场作用下的燃烧特性进行模拟计算,得到每个网格内固体结构的位移变化量和等效应力变化量。
本发明实施例通过模拟燃烧流场,可以快速有效的得到燃烧流场运动规律,利用流场分析结果得到的模拟计算结果对固体结构进行耦合分析,直观表达流场对固体结构产生的冲击影响,可以预测目标固体模型对应的零件内部在燃料燃烧时,零件的燃烧特性,进而对产品设计提供了一定的理论依据,更加有助于后期产品设计优化,提高产品性能。
附图说明
此处的附图被并入说明书中并构成本说明书的一部分,示出了符合本发明的实施例,并与说明书一起用于解释本发明的原理。
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,对于本领域普通技术人员而言,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本申请实施例提供的一种燃烧特性获取方法的流程图;
图2为本申请实施例提供的一种燃烧模型流场网格的示意图;
图3为本申请实施例提供的一种冷态流场甲烷质量分数云图;
图4为本申请实施例提供的一种H2O摩尔分数云图;
图5为本申请实施例提供的一种燃烧流场温度云图;
图6为本申请实施例提供的一种燃烧流场速度云图;
图7为本申请实施例提供的一种燃烧流场密度云图;
图8为本申请实施例提供的一种燃烧中心线温度分布曲线;
图9为本申请实施例提供的一种流场结果加载至热分析;
图10为本申请实施例提供的一种固体结构网格划分图;
图11为本申请实施例提供的一种热分析温度载荷云图;
图12为本申请实施例提供的一种稳态结构分析温度载荷云图;
图13为本申请实施例提供的一种稳态结构分析压力载荷云图;
图14为本申请实施例提供的一种燃烧室位移云图;
图15为本申请实施例提供的一种燃烧室等效应力云图;
图16为本申请实施例提供的一种燃烧特性获取装置的结构图;
图17为本申请实施例提供的一种电子设备的结构图。
具体实施方式
为使本申请实施例的目的、技术方案和优点更加清楚,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本申请的一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都属于本申请保护的范围。
由于目前对于设计完毕的产品的实验方法是首先生产出实体的产品,对产品施加实际的环境因素,如:向产品内部腔体注入混合气体,并点燃注入腔体的混合气体等,然而,这种产品设计方式,效率非常低,且造成资源的浪费。为此,本发明实施例提供了一种燃烧特性获取方法、装置、电子设备及存储介质,本发明实施例通过主要使用流-热-固耦合分析的方式,获取目标固体模型的燃烧特性,流-热-固耦合分析主要针对流体与固体之间的传热以及固体受热之后产生的热应力和热应变的分析。流-热-固耦合分为瞬态和稳态,通过流-热-固耦合分析,可以实时了解流场中的热力对于固体结构产生的应力以及应变,可以较为准确的得到产品在结构上需要改进的方向,可以快速有效的提高产品性能,提升产品市场竞争力。
图1为本申请实施例提供的一种燃烧特性获取方法的流程图,如图1所示,所述燃烧特性获取方法可以包括以下步骤:
步骤S101,对与目标固体模型对应的燃烧模型进行网格划分;
本发明实施例中,目标模型为待获取燃烧特性的零件构建的,示例性的,零件可以指发动机等,燃烧模型为基于目标固体模型内部流体区域构建的模型,用于分析工程中因燃烧流场产生的热力导致机械结构产生形变、断裂破坏等问题。
在该步骤中,可以将燃烧模型的进行网格划分,得到的燃烧模型流场网格如图2所示,以便对每个网格进行模拟计算,具体而言,可以使用Gambit网格划分软件采用四面体网格对所述燃烧模型进行网格划分,得到网格文件,将网格文件导入FLUENT软件。
网格划分的作用是为了使模型成为有限元,划分网格后,单元节点的位移增量是迭代方程计算过程中的基本未知量,使用数值积分方法对网格单元进行模拟求解。
在对与目标固体模型对应的燃烧模型进行网格划分之后,所述方法还包括:
检查所述网格文件中网格的质量,将所述网格文件中网格的长度单位修改为预设单位;设置所述燃烧模型的反应类型;将所述燃烧模型设置为涡耗散模型;设置所述燃烧模型考虑湍流流动中的化学反应和混合比率。
检查所述网格文件中网格的质量的作用是提高计算精度,减少收敛时间,Flunet软件操作流程为:mesh→check。检查结果包括element quality和aspect ratio。elementquality单元数值越大越好,一般不能低于0.7,aspect ratio单元数值越接近1越好。
由于流道模型建模之初是以mm为单位进行建模,导入Fluent软件后默认单位为m,模型单位应与之前保持一致,否则会影响计算结果的准确性,所以需要将所述网格文件中网格的长度单位修改为预设单位。
有限速率/涡耗散模型对湍流燃烧模型既考虑阿列纽斯定律,又考虑混合比率,并且在计算中取两者最小值,此模型考虑混合气体的混合比率,较符合当前工况。
步骤S102,对网格划分后的所述燃烧模型的燃烧流场进行燃烧模拟计算,得到模拟计算结果。
本发明实施例中的燃烧流场为燃料燃烧产生的流场,通过模拟燃烧流场,可以得到燃料在零件内部燃烧时,零件内部的流体区域的流场。
步骤S103,将所述模拟计算结果与所述目标固体模型关联,对关联后的目标固体模型进行网格划分;
该步骤中,可以将模拟计算结果导入目标固体模型,并对目标固体模型本身进行网格划分,以便于分析燃料在目标固体模型对应的零件内燃烧时对零件的影响。
步骤S104,对网格划分后的所述目标固体模型在所述燃烧流场作用下的燃烧特性进行模拟计算,得到每个网格内固体结构的位移变化量和等效应力变化量。
通过对目标固体模型在所述燃烧流场作用下的燃烧特性进行模拟计算,模拟燃料在燃烧时对零件的影响,得到每个网格内固体结构的位移变化量和等效应力变化量。
本发明实施例通过首先对与目标固体模型对应的燃烧模型进行网格划分,然后对网格划分后的所述燃烧模型的燃烧流场进行燃烧模拟计算,得到模拟计算结果,燃烧流场为燃料燃烧产生的流场,再将所述模拟计算结果与所述目标固体模型关联,对关联后的目标固体模型进行网格划分,最后可以对网格划分后的所述目标固体模型在所述燃烧流场作用下的燃烧特性进行模拟计算,得到每个网格内固体结构的位移变化量和等效应力变化量。
本发明实施例通过模拟燃烧流场,可以快速有效的得到燃烧流场运动规律,利用流场分析结果得到的模拟计算结果对固体结构进行耦合分析,直观表达流场对固体结构产生的冲击影响,可以预测目标固体模型对应的零件内部在燃料燃烧时,零件的燃烧特性,进而对产品设计提供了一定的理论依据,更加有助于后期产品设计优化,提高产品性能。
在本发明的又一实施例中,步骤S102对网格划分后的所述燃烧模型的燃烧流场进行燃烧模拟计算,得到模拟计算结果,包括:
步骤201,设置所述燃烧模型的边界条件、松弛因子和收敛精度;
边界条件可以指燃料进口速度、空气进口速度、压力出口以及湍流粘度等。示例性的,燃料进口以及空气进口设置为速度进口,燃料进口速度为5m/s,空气进口为1m/s。
本发明实施例中,还进行燃料进口组分以及空气进口组分为流体成分设置,目的在于模拟真实工况中的气体组成,燃料进口组分为100%甲烷,空气进口为21%氧气,Specification Method均为Intensity and Viscosity Ratio,Turbulent Intensity均为5%,Turbulent Viscosity均为5,出口设置为压力出口,壁面设为绝热壁面,其它保持默认设置。
设置松弛因子,将Momentum设为0.5,Turbulent设为0.5,CO2和O2设为0.8,Energy设为0.85。
收敛精度为e-03,e-03代表10的-3次方,设置收敛精度的目的在于提高计算精度,保证结果的准确性。
步骤202,对网格划分后的所述燃烧模型的冷态流场进行模拟计算,所述冷态流场为燃料与空气混合且燃料未燃烧放热的流场;
本发明实施例中,可以使用默认的收敛精度,从燃料进口进行流场初始化,开始计算冷态流场,冷态流场是没有燃烧放热的流场,只有燃料与空气组分传递和混合扩散,冷态流场常作为燃烧计算的初始流场,对于内部涡流预测以及组分预测具有一定的实际意义,图3为冷态流场甲烷质量分数云图。
冷态流场气体质量分数分布可以对流场内部涡流分布以及组分混合情况进行预测。
步骤203,基于模拟计算后的所述冷态流场,对网格划分后的所述燃烧模型的燃烧流场进行模拟计算,得到模拟计算结果。
在该步骤中,设置基于模拟计算后的所述冷态流场的初始温度及内部气体浓度,将已设置初始温度及内部气体浓度的所述冷态流场作为所述燃烧流场的初始流场,模拟计算所述燃烧模型的燃烧流场。
冷态流场计算结束后,对燃烧流场进行计算,首先要赋予流场较高的温度,初始温度设为2000K,二氧化碳以及氧气浓度设为20%,开始进行计算求解。图4为分析结果得到的H2O摩尔分数云图,可以看出最高为18.3%。图5为燃烧流场的温度云图,最高温度为2640K,图6为燃烧流场内甲烷与空气的流动速度云图,速度最高为6.13m/s,图7为燃烧流场内甲烷与空气的混合气体的密度云图,密度最高为1.17kg/m3,图8为燃烧中心线温度分布曲线,表示各个位置流场稳态温度变化。
图4-图8揭示了流体在运动的过程中,流体温度以及速度等分布情况,可以表现哪些位置受到流体的冲击较严重,哪些位置受到的热应力较大,对于零件的破坏较严重。
本发明实施例能够模拟燃烧模型的燃烧流场,便于后续分析燃烧流场对目标固体模型的影响。
在本发明的又一实施例中,将所述模拟计算结果与所述目标固体模型关联,对关联后的目标固体模型进行网格划分,包括:
步骤301,获取所述目标固体模型;
步骤302,将所述模拟计算结果导入所述目标固体模型;
该步骤中,可以流场分析时得到的包含压力载荷以及温度载荷数据的模拟计算结果导入固体模型。
示例性的,打开结构分析模块,选择import load,加载pressure load和temperature load,然后可以通过图9的方式关联分析数据,图9中results选项,鼠标右键可以显示将数据关联至结构分析模块。
步骤303,将已导入所述模拟计算结果的所述目标固体模型进行网格划分。
如图9所示,可以根据实际情况中零件材料确定采用结构钢材料,然后进行模型网格划分,设置体网格大小10mm,自动生成网格,如图10所示。
在本发明的又一实施例中,对网格划分后的所述目标固体模型在所述燃烧流场作用下的燃烧特性进行模拟计算,得到每个网格内固体结构的位移变化量和等效应力变化量,包括:
步骤401,设置所述目标固体模型的边界条件;
本发明实施例中,边界条件为空气端面设置固定支撑以及对固体结构施加压力载荷以及温度载荷。
步骤402,设置位移求解器及等效应力求解器;
示例性的,分析模块里面有一个solve部分,右键插入deform和Von misesstress,即位移以及等效应力求解器。
步骤403,对所述目标固体模型的燃烧特性进行模拟计算,得到每个网格内固体结构的位移变化量和等效应力变化量。
将流体分析结果导入固体结构,空气口端面设置为Fix Support,设置固定支撑的作用在于根据实际工况限制模型的自由度,防止计算过程中模型受力飞走。导入压力载荷以及温度载荷的目的在于在燃烧过程中模型受到这两种力,加载之后进行分析更加真实可靠。
图11、图12、图13分别为稳态热分析温度载荷云图、稳态结构分析温度载荷云图、稳态结构分析压力载荷云图,温度载荷最高为2355.9℃。最后设置位移、等效应力等求解器,开始进行程序计算。其中,图11-13为导入pressure load和temperature load后自动生成的。
求解完成之后燃烧室位移云图以及等效应力云图分析结果如图14、15所示,最大等效应力为4723.7MPa,最大位移为1.356mm。
得到的应力以及位移计算结果可以看出最大应力及最大变形位置,模型材料有相应的屈服强度,超过屈服强度会使零件被破坏,根据分析结果对零件比较脆弱的部位进行加强,对产品零件进行结构上的优化,提高零件的使用性能。
在本发明的又一实施例中,还提供一种燃烧特性获取装置,如图16所示,所述装置包括:
网格划分模块11,用于对与目标固体模型对应的燃烧模型进行网格划分;
第一模拟计算模块12,用于对网格划分后的所述燃烧模型的燃烧流场进行燃烧模拟计算,得到模拟计算结果,所述燃烧流场为燃料燃烧产生的流场;
关联划分模块13,用于将所述模拟计算结果与所述目标固体模型关联,对关联后的目标固体模型进行网格划分;
第二模拟计算模块14,用于对网格划分后的所述目标固体模型在所述燃烧流场作用下的燃烧特性进行模拟计算,得到每个网格内固体结构的位移变化量和等效应力变化量。
在本发明的又一实施例中,还提供一种电子设备,包括处理器、通信接口、存储器和通信总线,其中,处理器,通信接口,存储器通过通信总线完成相互间的通信;
存储器,用于存放计算机程序;
处理器,用于执行存储器上所存放的程序时,实现前述任一方法实施例所述的燃烧特性获取方法。
本发明实施例提供的电子设备,处理器通过执行存储器上所存放的程序实现了通过首先对与目标固体模型对应的燃烧模型进行网格划分,然后对网格划分后的所述燃烧模型的燃烧流场进行燃烧模拟计算,得到模拟计算结果,燃烧流场为燃料燃烧产生的流场,再将所述模拟计算结果与所述目标固体模型关联,对关联后的目标固体模型进行网格划分,最后可以对网格划分后的所述目标固体模型在所述燃烧流场作用下的燃烧特性进行模拟计算,得到每个网格内固体结构的位移变化量和等效应力变化量。
本发明实施例通过模拟燃烧流场,可以快速有效的得到燃烧流场运动规律,利用流场分析结果得到的模拟计算结果对固体结构进行耦合分析,直观表达流场对固体结构产生的冲击影响,可以预测目标固体模型对应的零件内部在燃料燃烧时,零件的燃烧特性,进而对产品设计提供了一定的理论依据,更加有助于后期产品设计优化,提高产品性能。
上述电子设备提到的通信总线1140可以是外设部件互连标准(PeripheralComponentInterconnect,简称PCI)总线或扩展工业标准结构(ExtendedIndustryStandardArchitecture,简称EISA)总线等。该通信总线1140可以分为地址总线、数据总线、控制总线等。为便于表示,图17中仅用一条粗线表示,但并不表示仅有一根总线或一种类型的总线。
通信接口1120用于上述电子设备与其他设备之间的通信。
存储器1130可以包括随机存取存储器(RandomAccessMemory,简称RAM),也可以包括非易失性存储器(non-volatilememory),例如至少一个磁盘存储器。可选的,存储器还可以是至少一个位于远离前述处理器的存储装置。
上述的处理器1110可以是通用处理器,包括中央处理器(CentralProcessingUnit,简称CPU)、网络处理器(NetworkProcessor,简称NP)等;还可以是数字信号处理器(DigitalSignalProcessing,简称DSP)、专用集成电路(ApplicationSpecificIntegratedCircuit,简称ASIC)、现场可编程门阵列(Field-ProgrammableGateArray,简称FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。
在本发明的又一实施例中,还提供一种计算机可读存储介质,所述计算机可读存储介质上存储有燃烧特性获取方法的程序,所述燃烧特性获取方法的程序被处理器执行时实现前述任一方法实施例所述的燃烧特性获取方法的步骤。
需要说明的是,在本文中,诸如“第一”和“第二”等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
以上所述仅是本发明的具体实施方式,使本领域技术人员能够理解或实现本发明。对这些实施例的多种修改对本领域的技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其它实施例中实现。因此,本发明将不会被限制于本文所示的这些实施例,而是要符合与本文所申请的原理和新颖特点相一致的最宽的范围。
Claims (10)
1.一种燃烧特性获取方法,其特征在于,包括:
对与目标固体模型对应的燃烧模型进行网格划分;
对网格划分后的所述燃烧模型的燃烧流场进行燃烧模拟计算,得到模拟计算结果,所述燃烧流场为燃料燃烧产生的流场;
将所述模拟计算结果与所述目标固体模型关联,对关联后的目标固体模型进行网格划分;
对网格划分后的所述目标固体模型在所述燃烧流场作用下的燃烧特性进行模拟计算,得到每个网格内固体结构的位移变化量和等效应力变化量。
2.根据权利要求1所述的燃烧特性获取方法,其特征在于,对网格划分后的所述燃烧模型的燃烧流场进行燃烧模拟计算,得到模拟计算结果,包括:
设置所述燃烧模型的边界条件、松弛因子和收敛精度;
对网格划分后的所述燃烧模型的冷态流场进行模拟计算,所述冷态流场为燃料与空气混合且燃料未燃烧放热的流场;
基于模拟计算后的所述冷态流场,对网格划分后的所述燃烧模型的燃烧流场进行模拟计算,得到模拟计算结果。
3.根据权利要求2所述的燃烧特性获取方法,其特征在于,所述基于模拟计算后的所述冷态流场,对网格划分后的所述燃烧模型的燃烧流场进行模拟计算,得到模拟计算结果,包括:
设置基于模拟计算后的所述冷态流场的初始温度及内部气体浓度;
将已设置初始温度及内部气体浓度的所述冷态流场作为所述燃烧流场的初始流场,模拟计算所述燃烧模型的燃烧流场。
4.根据权利要求1所述的燃烧特性获取方法,其特征在于,将所述模拟计算结果与所述目标固体模型关联,对关联后的目标固体模型进行网格划分,包括:
获取所述目标固体模型;
将所述模拟计算结果导入所述目标固体模型;
将已导入所述模拟计算结果的所述目标固体模型进行网格划分。
5.根据权利要求1所述的燃烧特性获取方法,其特征在于,对网格划分后的所述目标固体模型在所述燃烧流场作用下的燃烧特性进行模拟计算,得到每个网格内固体结构的位移变化量和等效应力变化量,包括:
设置所述目标固体模型的边界条件;
设置位移求解器及等效应力求解器;
对所述目标固体模型的燃烧特性进行模拟计算,得到每个网格内固体结构的位移变化量和等效应力变化量。
6.根据权利要求1所述的燃烧特性获取方法,其特征在于,对与目标固体模型对应的燃烧模型进行网格划分,包括:
使用Gambit网格划分软件采用四面体网格对所述燃烧模型进行网格划分,得到网格文件;
将网格文件导入FLUENT软件。
7.根据权利要求6所述的燃烧特性获取方法,其特征在于,在对与目标固体模型对应的燃烧模型进行网格划分之后,所述方法还包括:
检查所述网格文件中网格的质量,将所述网格文件中网格的长度单位修改为预设单位;
设置所述燃烧模型的反应类型;
将所述燃烧模型设置为涡耗散模型;
设置所述燃烧模型考虑湍流流动中的化学反应和混合比率。
8.一种燃烧特性获取装置,其特征在于,包括:
网格划分模块,用于对与目标固体模型对应的燃烧模型进行网格划分;
第一模拟计算模块,用于对网格划分后的所述燃烧模型的燃烧流场进行燃烧模拟计算,得到模拟计算结果,所述燃烧流场为燃料燃烧产生的流场;
关联划分模块,用于将所述模拟计算结果与所述目标固体模型关联,对关联后的目标固体模型进行网格划分;
第二模拟计算模块,用于对网格划分后的所述目标固体模型在所述燃烧流场作用下的燃烧特性进行模拟计算,得到每个网格内固体结构的位移变化量和等效应力变化量。
9.一种电子设备,其特征在于,包括处理器、通信接口、存储器和通信总线,其中,处理器,通信接口,存储器通过通信总线完成相互间的通信;
存储器,用于存放计算机程序;
处理器,用于执行存储器上所存放的程序时,实现权利要求1~7任一所述的燃烧特性获取方法。
10.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质上存储有燃烧特性获取方法的程序,所述燃烧特性获取方法的程序被处理器执行时实现权利要求1-7任一所述的燃烧特性获取方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110468022.0A CN113177374A (zh) | 2021-04-28 | 2021-04-28 | 一种燃烧特性获取方法、装置、电子设备及存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110468022.0A CN113177374A (zh) | 2021-04-28 | 2021-04-28 | 一种燃烧特性获取方法、装置、电子设备及存储介质 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113177374A true CN113177374A (zh) | 2021-07-27 |
Family
ID=76926955
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110468022.0A Pending CN113177374A (zh) | 2021-04-28 | 2021-04-28 | 一种燃烧特性获取方法、装置、电子设备及存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113177374A (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070139444A1 (en) * | 2005-12-16 | 2007-06-21 | Dupont Keith A | Methods and apparatus predicting variations in material properties |
CN103235842A (zh) * | 2013-03-29 | 2013-08-07 | 广东电网公司电力科学研究院 | 四角切圆燃煤锅炉燃烧特性的获取方法和系统 |
CN105022900A (zh) * | 2015-08-19 | 2015-11-04 | 电子科技大学 | 基于热固耦合分析的重型数控立车静压转台结构优化方法 |
CN108363857A (zh) * | 2018-02-05 | 2018-08-03 | 燕山大学 | 蓄热式加热炉流场及工件内部温度、热应力分析方法 |
CN112036015A (zh) * | 2020-08-19 | 2020-12-04 | 太原理工大学 | 基于单项热流固耦合模拟仿真的锅体形态优化方法 |
CN112115651A (zh) * | 2020-08-26 | 2020-12-22 | 华南理工大学 | 一种车用气瓶在火烧状态下外部流场温度获取的优化方法 |
CN112329291A (zh) * | 2020-10-28 | 2021-02-05 | 中国电子科技集团公司第五十四研究所 | 一种反射面天线的热固耦合场计算方法 |
-
2021
- 2021-04-28 CN CN202110468022.0A patent/CN113177374A/zh active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070139444A1 (en) * | 2005-12-16 | 2007-06-21 | Dupont Keith A | Methods and apparatus predicting variations in material properties |
CN103235842A (zh) * | 2013-03-29 | 2013-08-07 | 广东电网公司电力科学研究院 | 四角切圆燃煤锅炉燃烧特性的获取方法和系统 |
CN105022900A (zh) * | 2015-08-19 | 2015-11-04 | 电子科技大学 | 基于热固耦合分析的重型数控立车静压转台结构优化方法 |
CN108363857A (zh) * | 2018-02-05 | 2018-08-03 | 燕山大学 | 蓄热式加热炉流场及工件内部温度、热应力分析方法 |
CN112036015A (zh) * | 2020-08-19 | 2020-12-04 | 太原理工大学 | 基于单项热流固耦合模拟仿真的锅体形态优化方法 |
CN112115651A (zh) * | 2020-08-26 | 2020-12-22 | 华南理工大学 | 一种车用气瓶在火烧状态下外部流场温度获取的优化方法 |
CN112329291A (zh) * | 2020-10-28 | 2021-02-05 | 中国电子科技集团公司第五十四研究所 | 一种反射面天线的热固耦合场计算方法 |
Non-Patent Citations (1)
Title |
---|
褚睿峰: "加热炉内燃烧传热及炉管热流固耦合数值模拟研究", 中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑, no. 06, pages 28 - 66 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Li et al. | Numerical prediction of combustion instability limit cycle oscillations for a combustor with a long flame | |
Yılmaz | Effect of swirl number on combustion characteristics in a natural gas diffusion flame | |
Zarvandi et al. | Numerical study of the effects of heat transfer methods on CH4/(CH4+ H2)-AIR pre-mixed flames in a micro-stepped tube | |
Visser | Generic Analysis Methods for Gas Turbine Engine Performance: The development of the gas turbine simulation program GSP | |
CN103364212A (zh) | 煤粉锅炉改烧烟煤变负荷燃烧特性测量方法与系统 | |
Krebs et al. | Comparison of nonlinear to linear thermoacoustic stability analysis of a gas turbine combustion system | |
Tham et al. | Development of a highly reduced mechanism for iso-octane HCCI combustion with targeted search algorithm | |
Zhu et al. | Three-dimensional numerical simulations of spherical flame evolutions in shock and reshock accelerated flows | |
Kang et al. | Transient Thermomechanical Stress Analysis of Hot Surface Ignition Device Using Sequentially Coupled Computational Fluid Dynamics–Finite Element Analysis Approach | |
CN112115630B (zh) | 一种预测炸药热损伤的热力耦合仿真方法 | |
CN113177374A (zh) | 一种燃烧特性获取方法、装置、电子设备及存储介质 | |
Jella et al. | Large eddy simulation of a pressurized, partially premixed swirling flame with finite-rate chemistry | |
Jatale et al. | A validation of flare combustion efficiency predictions from large eddy simulations | |
Rostampour et al. | Numerical investigation of the effect of knock on heat transfer in a turbocharged spark ignition engine | |
El Abbassi et al. | Modelling turbulent combustion coupled with conjugate heat transfer in openfoam | |
Gunadi et al. | TURBULENCE MODELS APPLICATION IN AIR FLOW OF CROSSFLOW TURBINE. | |
Wijeyakulasuriya et al. | Analyzing three-dimensional multiple shock-flame interactions in a constant-volume combustion channel | |
Detomaso et al. | Large Eddy Simulation of a pistonless constant volume combustor: a new concept of pressure gain combustion | |
Shih et al. | Modeling of internal reacting flows and external static stall flows using RANS and PRNS | |
Aboaba et al. | Smart proxy modeling application of artificial intelligence & machine learning in computational fluid dynamics | |
Razdolsky | Mathematical modeling of fire dynamics | |
Weathers et al. | An exercise in model validation: Comparing univariate statistics and Monte Carlo-based multivariate statistics | |
Yu et al. | Large eddy simulation of turbulent premixed piloted flame using artificial thickened flame model coupled with tabulated chemistry | |
Katragadda et al. | Modelling of the curvature term in the flame surface density transport equation: a direct numerical simulations based analysis | |
Wilson | Assessment and analysis of a novel intake for a ramjet engine |
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 |