用于尾气后处理的通孔式载体阻力系数计算方法技术

本发明专利技术提供一种用于尾气后处理的通孔式载体阻力系数计算方法，包括以下步骤：获取各温度、流量条件下对应的背压试验数据记为Pti,qj,其中i≥2，j≥2；根据通孔式载体的规格，获取通孔式载体内气体的实际流通面积s；根据各温度、流量试验条件下对应的背压试验数据，由气体状态方程计算得出各条件下流经通孔式载体的气体密度，记为ρti,qj；计算各温度条件下的动力粘度，记为μti；计算出各温度、流量试验条件下通孔式载体内部气流通道的平均气体流速Vti,qj；得到各温度、流量试验条件下通孔式载体单位长度压力损失值dpti,qj/dx；基于Forchheimer模型计算相同温度条件下每两组数据对应的阻力系数α和β值；求解最终的通孔式载体阻系数。该方法的误差小。

Calculation method of drag coefficient of through-hole carrier for tail gas aftertreatment

The invention provides a method for calculating the resistance coefficient of a through hole carrier for tail gas treatment, including the following steps: obtaining the corresponding back pressure test data of each temperature and flow condition is Pti, QJ, in which I is more than 2, and j is more than 2; according to the specification of the through hole carrier, the actual circulation area of the gas in the through hole carrier s is obtained; The corresponding back pressure test data under the temperature and flow test conditions are calculated by the gas state equation. The gas density of the flow through the through hole carrier is calculated and recorded as P Ti, QJ. The dynamic viscosity of each temperature condition is recorded as Mu Ti, and the air passage of the through hole carrier is calculated under the conditions of each temperature and flow test. The flow velocity of Vti and QJ is Vti, QJ, and the pressure loss value of the per unit length of through hole carrier is dpti and qj/dx under the conditions of each temperature and flow test. Based on the Forchheimer model, the resistance coefficient alpha and beta of every two groups of data under the same temperature are calculated, and the final resistance coefficient of the through hole carrier is solved. The error of this method is small.

【技术实现步骤摘要】
用于尾气后处理的通孔式载体阻力系数计算方法
本专利技术涉及一种数据计算方法，尤其是用于尾气后处理产品背压计算的通孔式载体阻力系数计算方法。
技术介绍
随着国家内燃机尾气污染物排放标准越来越严格，仅凭借内燃机的机内净化技术已经无法满足日益严苛的排放标准。在国五、国六阶段，根据不同的净化技术路线，需要安装一系列的后处理装置，这会在一定程度上增大背压，而背压对发动机扭矩、功率、油耗及污染物成分等均有较大影响。目前后处理装置系统背压评价基于仿真分析计算和台架试验，相较于台架试验，仿真分析计算具有节约成本、开发周期短等优点。其中CFD(计算流体动力学)仿真分析计算能考虑产品结构等因素对净化器的影响，该评价方式与台架试验的结果往往偏差较小。后处理装置系统背压的仿真分析评价的关键是算准净化器载体的背压，目前除DPF及GPF载体以外，其余净化器载体均为通孔式载体。市面上的CFD分析软件有很多，如ANSYSFluent、AVLFIRE与Star-CCM+等，但由于大部分是通用CFD软件，不能准确的考虑通孔式载体内真实的流场分布，且现有的通孔式载体阻力计算多依赖于经验参数，难以满足对尾气后处理装置系统背压计算的分析需求。
技术实现思路
本专利技术的目的克服现有技术中存在的不足，提供一种用于尾气后处理的通孔式载体阻力系数计算方法，该计算方法能准确的考虑通孔式载体内流场分布，计算得到的通孔式载体阻力系数运用于AVLFIRE计算流体分析软件，能准确的计算出含通孔式载体的后处理装置的系统背压及流场分布，满足对机动车后处理分析的需求。本专利技术采用的技术方案是：一种用于尾气后处理的通孔式载体阻力系数计算方法，包括以下步骤：步骤S1，获取i个不同温度点，各温度点下j个不同质量流量下的通孔式载体背压，各温度、流量条件下的对应的背压试验数据记为Pti,qj,其中i≥2，j≥2；步骤S2，根据通孔式载体的规格，获取通孔式载体内气体的实际流通面积s；步骤S3，根据各温度、流量试验条件下对应的背压试验数据，由气体状态方程计算得出各条件下流经通孔式载体的气体密度，记为ρti,qj；步骤S4，计算各温度条件下的动力粘度，记为μti；步骤S5，根据速度公式V＝q/(s·ρ)计算出各温度、流量试验条件下通孔式载体内部气流通道的平均气体流速Vti,qj；q为质量流量，s为通孔式载体内部气体的实际流通面积，ρ为流经通孔式载体的气体密度；步骤S6，各温度、流量试验条件下对应的背压试验数据除以通孔式载体长度，得到各温度、流量试验条件下通孔式载体单位长度压力损失值dpti,qj/dx,其中x表示通孔式载体的轴向长度；步骤S7，通孔式载体阻力计算基于Forchheimer模型其中α为粘性阻力系数，β为惯性阻力系数；在该Forchheimer模型方程中，在对各温度、流量试验条件下，dpti,qj/dx，μti，Vti,qj及ρti,qj均已在步骤S1～S6中获取，故仅α和β为待求量；将多组试验数据按照试验温度点的不同将数据分为i类，每类内含j组数据，将每类中的其中两组数据两两进行求逆矩阵，求出每两组数据对应的阻力系数α和β值；根据各α和β值的差异，筛选有效数据；步骤S8，对于步骤S7中的有效数据，将各温度、流量试验条件下获取的dpti,qj/dx，μti，Vti,qj及ρti,qj信息分别代入Forchheimer模型，则会产生多个二元一次方程的超静定方程组，方程组的未知量为x′＝(α，β)，求解该方程组的最小二乘解，该最小二乘解即为该超静定方程组最优解，该最优解即为计算得到的最终阻力系数α和β值。进一步地，步骤S1中，通孔式载体根据外观尺寸、目数、壁厚、催化剂涂敷状态选定规格；步骤S2中，将通孔式载体的外观尺寸、目数、壁厚、催化剂涂敷状态信息输入到一维热力学循环软件AVLBOOST中，以获取通孔式载体内气体的实际流通面积s；进一步地，步骤S4中，根据以下两方程μ·107＝0.3875·t+180.5-ABS((t-200)·8.5/200)，Pa·s，范围t＝0～400℃μ·107＝0.2725·t+220.5-ABS((t-600)·2.5/200)，Pa·s，范围t＝400～800℃计算出各温度条件下的动力粘度，记为μti；t为温度，μ为动力粘度。进一步地，步骤S7中，某类数据中所求得到的阻力系数α和β值存在偏差大于相应预设阈值的情况，则剔除该温度条件下的相关数据。本专利技术的优点在于：根据通孔式载体背压试验数据，计算得到阻力系数α和β值，将所获取的阻力系数代入AVLFIRE计算流体分析软件进行仿真计算分析，所得到的通孔式载体仿真计算背压结果与试验数据的误差能控制在5％以内，能满足工程运用的需求，该方法获取的通孔式载体阻力系数使得背压的计算值与试验数据有很好的一致性。目前载体技术也在日益发展，该方法能运用于各种形式的通孔式载体，规避了使用经验参数带来的不可控误差这一问题。附图说明图1为本专利技术的流程图。具体实施方式下面结合具体附图和实施例对本专利技术作进一步说明。按照本专利技术提供的技术方案，用于尾气后处理的通孔式载体阻力系数计算方法，所述计算方法包括如下步骤：步骤S1，进行所需某一外观尺寸、目数、壁厚、催化剂涂敷状态规格(例如φ101.6×152.4/300-6.5，催化剂涂敷量260g/L的规格)的通孔式载体背压测试试验，该背压测试试验在热流试验台上进行，不同温度、流量下的通孔式载体背压数据应尽可能多。假设试验在5个温度点t1、t2、t3、t4、t5(温度间隔100℃)下进行，在各温度点测量4个质量流量q1、q2、q3、q4(流量间隔200kg/h)下的通孔式载体背压，共计20组背压试验数据，各温度、流量下的对应的背压试验数据记为Pti,qj,其中i＝1,2,3,4,5；j＝1,2,3,4。步骤S2，将通孔式载体的外观尺寸、目数、壁厚、催化剂涂敷状态规格信息输入到一维热力学循环软件AVLBOOST中，在summary的GEOMETRY信息中获取通孔式载体内气体的实际流通面积s；步骤S3，根据各温度、流量试验条件下对应的背压试验数据，由气体状态方程计算得出各条件下流经通孔式载体的气体密度，记为ρti,qj；步骤S4，根据以下两方程μ·107＝0.3875·t+180.5-ABS((t-200)·8.5/200)，Pa·s，范围t＝0～400℃μ·107＝0.2725·t+220.5-ABS((t-600)·2.5/200)，Pa·s，范围t＝400～800℃计算出各温度条件下的动力粘度，记为μti；Pa·s是单位，帕·秒；t为温度，μ为动力粘度；步骤S5，根据速度公式V＝q/(s·ρ)计算出各温度、流量试验条件下通孔式载体内部气流通道的平均气体流速Vti,qj；q为质量流量，s为通孔式载体内部气体的实际流通面积，ρ为流经通孔式载体的气体密度；步骤S6，各温度、流量试验条件下对应的背压试验数据除以通孔式载体长度，可得到各温度、流量试验条件下通孔式载体单位长度压力损失值dpti,qj/dx,其中x表示通孔式载体的轴向长度；步骤S7，通孔式载体阻力计算基于Forchheimer模型其中α为粘性阻力系数，β为惯性阻力系数；在该Forchheimer模型方程中，在对各温度、流量试验本文档来自技高网...

【技术保护点】
1.一种用于尾气后处理的通孔式载体阻力系数计算方法，其特征在于，包括以下步骤：步骤S1，获取i个不同温度点，各温度点下j个不同质量流量下的通孔式载体背压，各温度、流量条件下的对应的背压试验数据记为Pti,qj,其中i≥2，j≥2；步骤S2，根据通孔式载体的规格，获取通孔式载体内气体的实际流通面积s；步骤S3，根据各温度、流量试验条件下对应的背压试验数据，由气体状态方程计算得出各条件下流经通孔式载体的气体密度，记为ρti,qj；步骤S4，计算各温度条件下的动力粘度，记为μti；步骤S5，根据速度公式V＝q/(s·ρ)计算出各温度、流量试验条件下通孔式载体内部气流通道的平均气体流速Vti,qj；q为质量流量，s为通孔式载体内部气体的实际流通面积，ρ为流经通孔式载体的气体密度；步骤S6，各温度、流量试验条件下对应的背压试验数据除以通孔式载体长度，得到各温度、流量试验条件下通孔式载体单位长度压力损失值dpti,qj/dx,其中x表示通孔式载体的轴向长度；步骤S7，通孔式载体阻力计算基于Forchheimer模型

【技术特征摘要】
1.一种用于尾气后处理的通孔式载体阻力系数计算方法，其特征在于，包括以下步骤：步骤S1，获取i个不同温度点，各温度点下j个不同质量流量下的通孔式载体背压，各温度、流量条件下的对应的背压试验数据记为Pti,qj,其中i≥2，j≥2；步骤S2，根据通孔式载体的规格，获取通孔式载体内气体的实际流通面积s；步骤S3，根据各温度、流量试验条件下对应的背压试验数据，由气体状态方程计算得出各条件下流经通孔式载体的气体密度，记为ρti,qj；步骤S4，计算各温度条件下的动力粘度，记为μti；步骤S5，根据速度公式V＝q/(s·ρ)计算出各温度、流量试验条件下通孔式载体内部气流通道的平均气体流速Vti,qj；q为质量流量，s为通孔式载体内部气体的实际流通面积，ρ为流经通孔式载体的气体密度；步骤S6，各温度、流量试验条件下对应的背压试验数据除以通孔式载体长度，得到各温度、流量试验条件下通孔式载体单位长度压力损失值dpti,qj/dx,其中x表示通孔式载体的轴向长度；步骤S7，通孔式载体阻力计算基于Forchheimer模型其中α为粘性阻力系数，β为惯性阻力系数；在该Forchheimer模型方程中，在对各温度、流量试验条件下，dpti,qj/dx，μti，Vti,qj及ρti,qj均已在步骤S1～S6中获取，故仅α和β为待求量；将多组试验数据按照试验温度点的不同将数据分为i类，每类内含j组数据，将每类中的其中两组数据两两进行求逆矩阵，求出每两组数据对应的阻力系数...

【专利技术属性】
技术研发人员：付细平，陈增响，马相雪，何伟娇，盛杰，何龙，张俊超，郭鑫，宫婷婷，
申请(专利权)人：无锡威孚力达催化净化器有限责任公司，
类型：发明
国别省市：江苏,32