WO2022099714A1 - Dynamic mesh method-based method for tow heating performance three-dimensional simulation in high temperature carbonization furnace - Google Patents

Dynamic mesh method-based method for tow heating performance three-dimensional simulation in high temperature carbonization furnace Download PDF

Info

Publication number
WO2022099714A1
WO2022099714A1 PCT/CN2020/129173 CN2020129173W WO2022099714A1 WO 2022099714 A1 WO2022099714 A1 WO 2022099714A1 CN 2020129173 W CN2020129173 W CN 2020129173W WO 2022099714 A1 WO2022099714 A1 WO 2022099714A1
Authority
WO
WIPO (PCT)
Prior art keywords
tow
carbonization furnace
option
temperature carbonization
mesh
Prior art date
Application number
PCT/CN2020/129173
Other languages
French (fr)
Chinese (zh)
Inventor
张毅鹏
孙中心
Original Assignee
天华化工机械及自动化研究设计院有限公司
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 天华化工机械及自动化研究设计院有限公司 filed Critical 天华化工机械及自动化研究设计院有限公司
Priority to PCT/CN2020/129173 priority Critical patent/WO2022099714A1/en
Publication of WO2022099714A1 publication Critical patent/WO2022099714A1/en

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
    • 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/08—Thermal analysis or thermal optimisation

Definitions

  • the invention relates to the technical field of high-temperature carbonization furnace design simulation analysis method used in carbon fiber production.
  • Carbon fiber production is a high-energy-consuming industry, among which the high-temperature carbonization furnace is one of the largest energy consumers in carbon fiber production equipment. At the same time, the high-temperature carbonization furnace is also the key equipment for carbon fiber production. Converted to carbon fiber with a carbon content greater than 90%.
  • High temperature carbonization furnace is the integration of high temperature technology and high temperature equipment, and the operating temperature is generally 1000°C-1600°C. At present, in the production of carbon fiber, the quality of the product is mainly tested after the production process is completed, resulting in that the quality of the product cannot be dynamically adjusted during production, and the product qualification rate is difficult to improve.
  • the uniformity of the temperature field in the high-temperature carbonization furnace is very important to the performance of the product, and the temperature in the furnace is above 1000 °C, and the conventional testing methods cannot meet the requirements. Therefore, a reasonable design method needs to be selected to ensure that the furnace wall can meet the specifications. Surface temperature, so that the tow is heated evenly in the furnace cavity, and the unit energy consumption is reduced at the same time.
  • the present invention is aimed at the lack of a three-dimensional heating performance that can ensure that the furnace wall can reach a surface temperature that meets the specification in the design process of the existing high-temperature carbonization furnace, so that the tow is heated evenly in the furnace cavity, and at the same time, the unit energy consumption is reduced.
  • the technology of the simulation method is insufficient, and a three-dimensional simulation method of the heating performance of the tow in the high temperature carbonization furnace based on the dynamic grid method is proposed.
  • a three-dimensional simulation method for tow heating performance in a high-temperature carbonization furnace based on a dynamic grid method characterized in that the method comprises the following steps:
  • step (2) transfer the fluid computational domain of the muffle cavity of the high-temperature carbonization furnace and the three-dimensional simulation model of the heated tow established in step (1) to the Global Mesh Setup module of the ICEM software, respectively.
  • the Global Mesh Setup module Volume Meshing Parameters Set the meshing type to Tetra/Mixed, and mesh the fluid computational domain of the muffle cavity of the high-temperature carbonization furnace and the 3D simulation model of the heated tow.
  • the entire computational domain uses unstructured meshes. The mesh is refined at the wall surface.
  • the mesh quality of the overall structure is guaranteed to be greater than 0.5, and the inlet and outlet and wall boundary names of all 3D simulation models are defined, including the tow wall surface and the inlet and outlet of the furnace cavity. and the name of the wall boundary of the furnace cavity;
  • v, ⁇ , ⁇ are the flow velocity, density and viscosity coefficient of the fluid, respectively, and d is the characteristic length;
  • the turbulence model is selected as the laminar model;
  • ⁇ -fluid density; t-time; V-velocity vector where u, v, w are the components of V in the three directions of x, y and z.
  • ⁇ is the dynamic viscosity
  • F b is the volume force on the micro-element
  • C p specific heat capacity
  • T temperature
  • k fluid heat transfer coefficient
  • S T viscous dissipation term
  • the parameters set for the fluid computational domain of the muffle cavity of the high-temperature carbonization furnace and the three-dimensional simulation model of the heated tow include: the geometry and size of the muffle cavity, the geometry and size of the inlet and outlet seals, Inlet and outlet sealing nitrogen pipe inlet size, nitrogen pipe outlet size, geometry and geometry of heated tow.
  • step (4) the center point of the three-dimensional simulation model of the fluid computational domain of the muffle cavity of the high-temperature carbonization furnace is selected as the detection surface, the detection surface is the X-direction plane passing through the center point, the surface of the tow is selected as the detection surface, and the detection surface is The X-direction plane passing through the center point.
  • the present invention can intuitively visualize the temperature change process of the tow in the carbon fiber production process by simulating the temperature field characteristics of the tow surface during the heating process in the high-temperature carbonization furnace, and then determine the furnace cavity. Whether the heating capacity of the tow meets the process requirements, thereby providing support for subsequent structural optimization.
  • the invention can reduce experimental cost, optimize product design, shorten product development cycle, provide theoretical support for reducing carbon fiber production energy consumption, and also provide basis for relevant numerical simulation research.
  • FIG. 1 is a schematic diagram of a three-dimensional model established in the simulation method of the present invention.
  • FIG. 2 is a flow chart of a three-dimensional simulation method for tow heating performance in a high-temperature carbonization furnace based on a dynamic grid method of the present invention.
  • FIG. 3 is a schematic diagram of the temperature distribution of the tow in the furnace cavity at any time of the present invention during the heating process.
  • Figure 4 is a schematic diagram of the temperature distribution on the surface of the tow of the present invention.
  • FIG. 5 is a schematic diagram of the temperature distribution in the furnace chamber of the present invention.
  • a three-dimensional simulation method for tow heating performance in a high-temperature carbonization furnace based on a dynamic grid method disclosed in the present invention specifically includes the following steps, as shown in FIG. 2 :
  • step (2) transfer the fluid computational domain of the muffle cavity of the high-temperature carbonization furnace and the three-dimensional simulation model of the heated tow established in step (1) to the Global Mesh Setup module of the ICEM software, respectively.
  • the Global Mesh Setup module Volume Meshing Parameters Set the meshing type to Tetra/Mixed, and mesh the fluid computational domain of the muffle cavity of the high-temperature carbonization furnace and the 3D simulation model of the heated tow.
  • the entire computational domain uses an unstructured mesh, that is, a tetrahedral mesh. grid. In order to ensure the accuracy of the flow field calculation results, the meshes near the surface of the tow are refined.
  • the mesh quality of the overall structure is guaranteed to be greater than 0.5.
  • v, ⁇ , ⁇ are the flow velocity, density and viscosity coefficient of the fluid, respectively, and d is the characteristic length;
  • the turbulence model is selected as the laminar model;
  • ⁇ -fluid density; t-time; V-velocity vector where u, v, w are the components of V in the three directions of x, y and z.
  • ⁇ is the dynamic viscosity
  • F b is the volume force on the micro-element
  • C p specific heat capacity
  • T temperature
  • k fluid heat transfer coefficient
  • S T viscous dissipation term
  • the detection surface is the X-direction plane passing through the center point, and the surface of the tow is selected as the detection surface, and the detection surface is the detection surface.
  • the X-direction plane of the center point it is preferable to select the center point of the three-dimensional simulation model of the fluid computational domain of the muffle cavity of the high-temperature carbonization furnace as the detection surface, the detection surface is the X-direction plane passing through the center point, and the surface of the tow is selected as the detection surface, and the detection surface is the detection surface.
  • the X-direction plane of the center point is preferable to select the center point of the three-dimensional simulation model of the fluid computational domain of the muffle cavity of the high-temperature carbonization furnace as the detection surface, the detection surface is the X-direction plane passing through the center point, and the surface of the tow is selected as the detection surface, and the detection surface is the detection surface.
  • the X-direction plane of the center point it is preferable to select the center point of the three-dimensional simulation model
  • the post-processing software CFD-POST is used to realize the visibility of the temperature field in the furnace cavity and the surface temperature field of the tow during the tow heating process, as shown in Figure 4 and Figure 5. It can be seen from the figure that not only the accuracy of the three-dimensional calculation results can fully meet the requirements, but also the types of output data are diverse, and the output results are more intuitive.
  • the three-dimensional simulation method for tow heating performance in a high-temperature carbonization furnace based on the dynamic grid method provided by the present invention has the following advantages:
  • the traditional method is usually to test the performance of the tow after heating.
  • the invention creatively uses the moving grid method to simulate the motion state of the tow, and judges the heating capacity of the furnace through the cloud image of the temperature change on the surface of the tow, so that the structure design can be better measured.
  • the invention defines the motion law of the tow through UDF, and by modifying the airflow velocity in the inlet boundary velocity-inlet, the changing law of the temperature distribution on the surface of the tow and the heat storage capacity in the muffle furnace under different tow motion states can be obtained.
  • the temperature distribution on the surface of the tow is analyzed. The temperature of the tow increases gradually with the length of the tow entering the furnace cavity, and the temperature of the furnace cavity remains constant at all times.
  • the surface temperature of the tow gradually increases during the process of entering the furnace cavity, the tow is heated, the furnace wall has a good heating effect, and its heat storage performance is superior; as shown in Figure 3, from any time the tow is in the furnace
  • the schematic diagram of the temperature distribution in the heating process in the cavity and the schematic diagram of the temperature distribution on the surface of the tow as shown in Figure 4 show that the temperature of the tow surface is heated in a gradient increase, and the temperature distribution on the surface of the tow is uniform; as shown in Figure 5, the present invention
  • the schematic diagram of the temperature distribution in the furnace cavity shows that the temperature distribution in the muffle cavity is uniform, indicating that the air flow in the furnace cavity is reasonably distributed.
  • multiple simulations should be performed, and the experimental results should be compared and analyzed to obtain the best solution for the speed of the tow movement and the heat storage performance of the muffle furnace cavity.

Abstract

A dynamic mesh method-based method for tow heating performance three-dimensional simulation in a high temperature carbonization furnace, which comprises: (1), constructing a three-dimensional simulation model; (2), transmitting to a Global Mesh Setup module of ICEM software, and performing meshing; (3), transmitting to a FLUENT module of ANSYS software, and configuring a boundary condition; (4), performing simulation computation and obtaining a result, and using same to serve as an indicator for determining tow heating performance in a high temperature carbonization furnace; and (5), repeating steps (1)-(4), achieving a result having visual, directly observable output of a temperature field within a furnace cavity and a tow surface temperature field during a tow heating process via a temperature change cloud map of a surface of measurement. Costs are reduced, a product development cycle is shortened, and a basis is provided for numerical simulation research.

Description

一种基于动网格方法的高温碳化炉中丝束加热性能三维模拟方法A three-dimensional simulation method of tow heating performance in high temperature carbonization furnace based on dynamic grid method 技术领域technical field
本发明涉及到碳纤维生产中所用的高温碳化炉设计模拟分析方法技术领域。The invention relates to the technical field of high-temperature carbonization furnace design simulation analysis method used in carbon fiber production.
背景技术Background technique
碳纤维生产属于高耗能产业，其中的高温碳化炉是碳纤维生产设备中的耗能大户之一，同时，高温碳化炉也是碳纤维生产的关键设备，主要用于对预氧丝进行高温碳化，使其转化为碳元素含量大于90％的碳纤维。高温碳化炉是高温技术和高温设备的集成，使用温度一般在1000℃-1600℃。目前在碳纤维生产中，产品的质量主要在生产流程完成后进行检测，导致产品的质量不能在生产中进行动态调整，产品合格率难以提高。高温碳化炉内温度场的均匀性对产品的性能至关重要，并且炉内温度在1000℃以上，常规的测试手段不能满足要求，因此需要选择合理的设计方法，保证炉壁能达到符合规范的表面温度，使丝束在炉腔内受热均匀，同时降低单位能耗。Carbon fiber production is a high-energy-consuming industry, among which the high-temperature carbonization furnace is one of the largest energy consumers in carbon fiber production equipment. At the same time, the high-temperature carbonization furnace is also the key equipment for carbon fiber production. Converted to carbon fiber with a carbon content greater than 90%. High temperature carbonization furnace is the integration of high temperature technology and high temperature equipment, and the operating temperature is generally 1000℃-1600℃. At present, in the production of carbon fiber, the quality of the product is mainly tested after the production process is completed, resulting in that the quality of the product cannot be dynamically adjusted during production, and the product qualification rate is difficult to improve. The uniformity of the temperature field in the high-temperature carbonization furnace is very important to the performance of the product, and the temperature in the furnace is above 1000 °C, and the conventional testing methods cannot meet the requirements. Therefore, a reasonable design method needs to be selected to ensure that the furnace wall can meet the specifications. Surface temperature, so that the tow is heated evenly in the furnace cavity, and the unit energy consumption is reduced at the same time.
发明内容SUMMARY OF THE INVENTION
综上所述，本发明是针对现有高温碳化炉设计过程中，缺少能保证炉壁能达到符合规范的表面温度，使丝束在炉腔内受热均匀，同时降低单位能耗的加热性能三维模拟方法的技术不足，而提出一种基于动网格方法的高温碳化炉中丝束加热性能三维模拟方法。To sum up, the present invention is aimed at the lack of a three-dimensional heating performance that can ensure that the furnace wall can reach a surface temperature that meets the specification in the design process of the existing high-temperature carbonization furnace, so that the tow is heated evenly in the furnace cavity, and at the same time, the unit energy consumption is reduced. The technology of the simulation method is insufficient, and a three-dimensional simulation method of the heating performance of the tow in the high temperature carbonization furnace based on the dynamic grid method is proposed.
为解决本发明所提出的技术问题，采用的技术方案为：In order to solve the technical problem proposed by the present invention, the technical scheme adopted is:
一种基于动网格方法的高温碳化炉中丝束加热性能三维模拟方法，其特征在于所述方法包括有如下步骤：A three-dimensional simulation method for tow heating performance in a high-temperature carbonization furnace based on a dynamic grid method, characterized in that the method comprises the following steps:
(1)、采用三维CAD软件SOLIDWORKS建立高温碳化炉马弗腔体流体计算域和被加热丝束的三维仿真模型；(1) Using the 3D CAD software SOLIDWORKS to establish the 3D simulation model of the fluid computational domain of the muffle cavity of the high temperature carbonization furnace and the heated tow;
(2)、将步骤(1)建立的高温碳化炉马弗腔体流体计算域和被加热丝束的三维仿真模型分别传递到ICEM软件的Global Mesh Setup模块中，在Global Mesh Setup模块中Volume Meshing Parameters设置网格划分类型为Tetra/Mixed，对高温碳化炉马弗腔体流体计算域和被加热丝束的三维仿真模型进行网格划分，整个计算域采用非结构网格，对靠近丝束表面壁面处进行网格加密，根据ICEM软件中网格质量的判断标准，保证整体结构的网格质量大于0.5，定义所有三维仿真 模型的进出口与壁面边界名称，包括丝束壁面、炉腔进出口和炉腔壁面边界名称；(2), transfer the fluid computational domain of the muffle cavity of the high-temperature carbonization furnace and the three-dimensional simulation model of the heated tow established in step (1) to the Global Mesh Setup module of the ICEM software, respectively. In the Global Mesh Setup module, Volume Meshing Parameters Set the meshing type to Tetra/Mixed, and mesh the fluid computational domain of the muffle cavity of the high-temperature carbonization furnace and the 3D simulation model of the heated tow. The entire computational domain uses unstructured meshes. The mesh is refined at the wall surface. According to the grid quality judgment standard in the ICEM software, the mesh quality of the overall structure is guaranteed to be greater than 0.5, and the inlet and outlet and wall boundary names of all 3D simulation models are defined, including the tow wall surface and the inlet and outlet of the furnace cavity. and the name of the wall boundary of the furnace cavity;
(3)、将划分好网格的流体计算域和被加热丝束的三维仿真模型传递到ANSYS软件的FLUENT模块，并进行设置边界条件；(3) Transfer the meshed fluid computational domain and the three-dimensional simulation model of the heated tow to the FLUENT module of the ANSYS software, and set the boundary conditions;
(4)、在ANSYS软件中的FLUENT模块里设置温度检测面，并进行仿真运算得到结果，以此作为判定高温碳化炉中丝束加热性能的指标；ANSYS软件中的FLUENT模块进行设置的过程如下：(4) Set the temperature detection surface in the FLUENT module in the ANSYS software, and perform simulation operations to obtain the results, which are used as indicators to determine the heating performance of the tow in the high-temperature carbonization furnace; the process of setting the FLUENT module in the ANSYS software is as follows :
(4.1)、在User Defined选项导入根据工艺参数编制的自定义丝束运动参数，通过UDF控制被加热丝束的运动状态，实现丝束在高温碳化炉中运动过程；(4.1) Import the custom tow motion parameters compiled according to the process parameters in the User Defined option, and control the motion state of the heated tow through UDF to realize the motion process of the tow in the high temperature carbonization furnace;
(4.2)、在General选项中，将y方向Gravitational Acceleration根据实际设定为预设值，time选项设置为Transient瞬态传热；(4.2) In the General option, set the Gravitational Acceleration in the y direction to the preset value according to the actual situation, and set the time option to Transient transient heat transfer;
(4.3)、将Models选项中的Energy勾选Energy Equation，Viscous Models选项中选取laminar模型，为了判断炉腔内气流运动状态，引入雷诺数进行描述，雷诺数的计算公式为：(4.3) Check the Energy Equation in the Models option, and select the laminar model in the Viscous Models option. In order to judge the airflow movement state in the furnace cavity, the Reynolds number is introduced for description. The calculation formula of the Reynolds number is:
Figure PCTCN2020129173-appb-000001
Figure PCTCN2020129173-appb-000001
其中，其中v、ρ、μ分别为流体的流速、密度与黏性系数，d为特征长度；通过雷诺数的计算，进而选择湍流模型为laminar模型；Among them, v, ρ, μ are the flow velocity, density and viscosity coefficient of the fluid, respectively, and d is the characteristic length; through the calculation of the Reynolds number, the turbulence model is selected as the laminar model;
(4.4)在Materials Fluid选项部分选择空气和氮气；在Materials Solid选项建立丝束的物理参数，包括密度、比热容和热传导率参数；(4.4) Select air and nitrogen in the Materials Fluid option section; establish the physical parameters of the tow in the Materials Solid option, including density, specific heat capacity and thermal conductivity parameters;
(4.5)、在Cell Zone Conditions选项中将Fluid1部分设为氮气，Fluid2部分设为氧气；将Solid1部分设为丝束；(4.5), in the Cell Zone Conditions option, set the Fluid1 part to nitrogen, and the Fluid2 part to oxygen; set the Solid1 part to tow;
(4.6)、在Boundary Conditions选项中设置入口边界条件为Pressure-inlet，并将Velocity Magnitude根据实际要求设置为预设值，Thermal选项设置为udf tm-inlet，设置出口边界条件为Pressure-outlet，将一侧壁面设置为对流换热面，UDF定义每小时内炉壁空气综合温度值，对流换热系数根据实际要求设置为预设值，其他壁面设置为绝热壁面，将丝束与气体接触面设置为Coupled；(4.6) In the Boundary Conditions option, set the inlet boundary condition to Pressure-inlet, set Velocity Magnitude to the preset value according to the actual requirements, set the Thermal option to udf tm-inlet, set the outlet boundary condition to Pressure-outlet, and set the One side wall is set as the convective heat transfer surface, UDF defines the comprehensive temperature value of the furnace wall air per hour, the convective heat transfer coefficient is set to the preset value according to the actual requirements, the other walls are set as the adiabatic wall surface, and the contact surface between the tow and the gas is set is Coupled;
(4.7)、在Dynamic Mesh选项，并勾选激活动网格选项，勾选Mesh Methods选项中的Smoothing，Remeshing选项；在Dynamic Mesh Zones中定义运动区域，在Zone Names下选择rigid wall，在Type中选择Rigid Body，在Preview Mesh  Motion中设置预览网格运动；(4.7), in the Dynamic Mesh option, and check the active mesh option, check the Smoothing and Remeshing options in the Mesh Methods option; define the motion area in Dynamic Mesh Zones, select rigid wall under Zone Names, and in Type Select Rigid Body and set preview mesh motion in Preview Mesh Motion;
(4.8)、选择Check Case后进行计算，计算时间根据实际工程中参数设定；其中选择Check Case后进行计算所使用的高温碳化炉全流场计算的三维数学模型，包括有连续性方程、动量方程和能量方程，分别如下所示：(4.8) Calculate after selecting Check Case, and the calculation time is set according to the parameters in the actual project; among them, the three-dimensional mathematical model of the full flow field calculation of the high-temperature carbonization furnace used for calculation after selecting Check Case, including continuity equation, momentum The equation and the energy equation, respectively, are as follows:
Figure PCTCN2020129173-appb-000002
Figure PCTCN2020129173-appb-000002
式中，ρ-流体密度；t-时间；V-速度矢量，其中u、v、w为V在x、y及z三个方向上的分量。In the formula, ρ-fluid density; t-time; V-velocity vector, where u, v, w are the components of V in the three directions of x, y and z.
动量方程的Navier-Stokes方程：Navier-Stokes equation for momentum equation:
Figure PCTCN2020129173-appb-000003
Figure PCTCN2020129173-appb-000003
其中，μ是动力黏度，F b是微元上的体积力； Among them, μ is the dynamic viscosity, F b is the volume force on the micro-element;
能量守恒方程：Energy conservation equation:
Figure PCTCN2020129173-appb-000004
Figure PCTCN2020129173-appb-000004
其中，C p—比热容，T—温度，k—流体传热系数，S T—粘性耗散项； Among them, C p —specific heat capacity, T—temperature, k—fluid heat transfer coefficient, S T —viscous dissipation term;
(5)、在相同设置条件下，通过将高温碳化炉马弗腔体流体计算域和被加热丝束的三维仿真模型设置不同参数，并重复步骤(1)-(4)，以进行多次模拟计算，由检测面的温度变化云图，利用后处理软件CFD-POST实现丝束加热过程中炉腔内温度场、丝束表面温度场的可视性，直观输出的结果。(5) Under the same setting conditions, set different parameters by setting the fluid computational domain of the muffle cavity of the high temperature carbonization furnace and the three-dimensional simulation model of the heated tow, and repeat steps (1)-(4) for several times. For simulation calculation, from the temperature change cloud map of the detection surface, the post-processing software CFD-POST is used to realize the visibility of the temperature field in the furnace cavity and the surface temperature field of the tow during the heating process of the tow, and the results are intuitively output.
作为对本发明作进一步限定的技术方案包括有：As the technical scheme further limiting the present invention includes:
步骤(1)中对高温碳化炉马弗腔体流体计算域和被加热丝束的三维仿真模型设定的参数包括：马弗腔体几何形状和几何尺寸、进出口密封几何形状和几何尺寸、进出口密封氮气管进口尺寸、氮气管出口尺寸、被加热丝束的几何形状和几何尺寸。In step (1), the parameters set for the fluid computational domain of the muffle cavity of the high-temperature carbonization furnace and the three-dimensional simulation model of the heated tow include: the geometry and size of the muffle cavity, the geometry and size of the inlet and outlet seals, Inlet and outlet sealing nitrogen pipe inlet size, nitrogen pipe outlet size, geometry and geometry of heated tow.
步骤(4)中，选择高温碳化炉马弗腔体流体计算域的三维仿真模型的中心点作为检测面，检测面为过中心点的X方向平面，选择丝束表面作为检测面，检测面为过中心点的X方向平面。In step (4), the center point of the three-dimensional simulation model of the fluid computational domain of the muffle cavity of the high-temperature carbonization furnace is selected as the detection surface, the detection surface is the X-direction plane passing through the center point, the surface of the tow is selected as the detection surface, and the detection surface is The X-direction plane passing through the center point.
本发明的有益效果为：本发明通过对高温碳化炉中在加热过程中的丝束表面 温度场特性进行模拟，直观的可以将碳纤维生产过程中丝束的温度变化过程实现可视化，然后判断炉腔对丝束的加热能力是否满足工艺要求，从而为后续的结构优化提供支持。本发明可以降低实验成本，优化产品设计，缩短产品开发周期，为降低碳纤维生产能耗提供理论支持，也为相关的数值模拟研究提供依据。The beneficial effects of the present invention are as follows: the present invention can intuitively visualize the temperature change process of the tow in the carbon fiber production process by simulating the temperature field characteristics of the tow surface during the heating process in the high-temperature carbonization furnace, and then determine the furnace cavity. Whether the heating capacity of the tow meets the process requirements, thereby providing support for subsequent structural optimization. The invention can reduce experimental cost, optimize product design, shorten product development cycle, provide theoretical support for reducing carbon fiber production energy consumption, and also provide basis for relevant numerical simulation research.
附图说明Description of drawings
图1是本发明模拟方法中建立的三维模型示意图。FIG. 1 is a schematic diagram of a three-dimensional model established in the simulation method of the present invention.
图2是本发明一种基于动网格方法的高温碳化炉中丝束加热性能三维模拟方法的流程图。FIG. 2 is a flow chart of a three-dimensional simulation method for tow heating performance in a high-temperature carbonization furnace based on a dynamic grid method of the present invention.
图3是本发明任意时刻丝束在炉腔内加热过程的温度分布示意图。FIG. 3 is a schematic diagram of the temperature distribution of the tow in the furnace cavity at any time of the present invention during the heating process.
图4是本发明丝束表面的温度分布示意图。Figure 4 is a schematic diagram of the temperature distribution on the surface of the tow of the present invention.
图5是本发明炉腔内温度分布示意图。FIG. 5 is a schematic diagram of the temperature distribution in the furnace chamber of the present invention.
具体实施方式Detailed ways
以下结合附图和本发明优选的具体实施例对本发明的方法作进一步地说明。The method of the present invention will be further described below with reference to the accompanying drawings and preferred specific embodiments of the present invention.
本发明公开的一种基于动网格方法的高温碳化炉中丝束加热性能三维模拟方法，具体包括有如下步骤，如图2中所示：A three-dimensional simulation method for tow heating performance in a high-temperature carbonization furnace based on a dynamic grid method disclosed in the present invention specifically includes the following steps, as shown in FIG. 2 :
(1)、采用三维CAD(Computer Aided Design计算机辅助设计)软件SOLIDWORKS建立高温碳化炉马弗腔体流体计算域和被加热丝束的三维仿真模型，如图1中所示；并设定的相关参数包括：马弗腔体几何形状和几何尺寸、进出口密封几何形状和几何尺寸、进出口密封氮气管进口尺寸、氮气管出口尺寸、被加热丝束的几何形状和几何尺寸。(1) Use 3D CAD (Computer Aided Design) software SOLIDWORKS to establish a 3D simulation model of the fluid computational domain of the muffle cavity of the high temperature carbonization furnace and the heated tow, as shown in Figure 1; and set the relevant Parameters include: muffle cavity geometry and geometry, inlet and outlet seal geometry and geometry, inlet and outlet seal nitrogen pipe inlet size, nitrogen pipe outlet size, and heated tow geometry and geometry.
(2)、将步骤(1)建立的高温碳化炉马弗腔体流体计算域和被加热丝束的三维仿真模型分别传递到ICEM软件的Global Mesh Setup模块中，在Global Mesh Setup模块中Volume Meshing Parameters设置网格划分类型为Tetra/Mixed，对高温碳化炉马弗腔体流体计算域和被加热丝束的三维仿真模型进行网格划分，整个计算域采用非结构网格，也就是四面体网格。为了保证流场计算结果的准确性，对靠近丝束表面壁面处进行网格加密，根据ICEM软件中网格质量的判断标准，保证整体结构的网格质量大于0.5，时为了便于后期设置计算条件，定义所有三维仿真模型的进出口与壁面边界名称，包括丝束壁面、炉腔进出口和炉腔壁面边界名称。(2), transfer the fluid computational domain of the muffle cavity of the high-temperature carbonization furnace and the three-dimensional simulation model of the heated tow established in step (1) to the Global Mesh Setup module of the ICEM software, respectively. In the Global Mesh Setup module, Volume Meshing Parameters Set the meshing type to Tetra/Mixed, and mesh the fluid computational domain of the muffle cavity of the high-temperature carbonization furnace and the 3D simulation model of the heated tow. The entire computational domain uses an unstructured mesh, that is, a tetrahedral mesh. grid. In order to ensure the accuracy of the flow field calculation results, the meshes near the surface of the tow are refined. According to the grid quality judgment standard in ICEM software, the mesh quality of the overall structure is guaranteed to be greater than 0.5. In order to facilitate the later setting of calculation conditions , to define the boundary names of inlet, outlet and wall of all 3D simulation models, including the boundary names of the tow wall, inlet and outlet of the furnace cavity, and the wall surface of the furnace cavity.
(3)、将划分好网格的流体计算域和被加热丝束的三维仿真模型传递到ANSYS软件的FLUENT模块，并进行设置边界条件。(3) Transfer the meshed fluid computational domain and the three-dimensional simulation model of the heated tow to the FLUENT module of the ANSYS software, and set the boundary conditions.
(4)、在ANSYS软件中的FLUENT模块里设置温度检测面，并进行仿真运算得到结果，以此作为判定高温碳化炉中丝束加热性能的指标；ANSYS软件中的FLUENT模块进行设置的过程如下：(4) Set the temperature detection surface in the FLUENT module in the ANSYS software, and perform simulation operations to obtain the results, which are used as indicators to determine the heating performance of the tow in the high-temperature carbonization furnace; the process of setting the FLUENT module in the ANSYS software is as follows :
(4.1)、在User Defined选项导入根据工艺参数编制的自定义丝束运动参数，通过UDF控制被加热丝束的运动状态，实现丝束在高温碳化炉中运动过程；(4.1) Import the custom tow motion parameters compiled according to the process parameters in the User Defined option, and control the motion state of the heated tow through UDF to realize the motion process of the tow in the high temperature carbonization furnace;
(4.2)、在General选项中，将y方向Gravitational Acceleration根据实际设定为预设值，time选项设置为Transient瞬态传热；(4.2) In the General option, set the Gravitational Acceleration in the y direction to the preset value according to the actual situation, and set the time option to Transient transient heat transfer;
(4.3)、将Models选项中的Energy勾选Energy Equation，Viscous Models选项中选取laminar模型，为了判断炉腔内气流运动状态，引入雷诺数进行描述，雷诺数的计算公式为：(4.3) Check the Energy Equation in the Models option, and select the laminar model in the Viscous Models option. In order to judge the airflow movement state in the furnace cavity, the Reynolds number is introduced for description. The calculation formula of the Reynolds number is:
Figure PCTCN2020129173-appb-000005
Figure PCTCN2020129173-appb-000005
其中，其中v、ρ、μ分别为流体的流速、密度与黏性系数，d为特征长度；通过雷诺数的计算，进而选择湍流模型为laminar模型；Among them, v, ρ, μ are the flow velocity, density and viscosity coefficient of the fluid, respectively, and d is the characteristic length; through the calculation of the Reynolds number, the turbulence model is selected as the laminar model;
(4.4)在Materials Fluid选项部分选择空气和氮气；在Materials Solid选项建立丝束的物理参数，包括密度、比热容和热传导率参数；(4.4) Select air and nitrogen in the Materials Fluid option section; establish the physical parameters of the tow in the Materials Solid option, including density, specific heat capacity and thermal conductivity parameters;
(4.5)、在Cell Zone Conditions选项中将Fluid1部分设为氮气，Fluid2部分设为氧气；将Solid1部分设为丝束；(4.5), in the Cell Zone Conditions option, set the Fluid1 part to nitrogen, and the Fluid2 part to oxygen; set the Solid1 part to tow;
(4.6)、在Boundary Conditions选项中设置入口边界条件为Pressure-inlet，并将Velocity Magnitude根据实际要求设置为预设值，Thermal选项设置为udf tm-inlet，设置出口边界条件为Pressure-outlet，将一侧壁面设置为对流换热面，UDF定义每小时内炉壁空气综合温度值，对流换热系数根据实际要求设置为预设值，其他壁面设置为绝热壁面，将丝束与气体接触面设置为Coupled；(4.6) In the Boundary Conditions option, set the inlet boundary condition to Pressure-inlet, set Velocity Magnitude to the preset value according to the actual requirements, set the Thermal option to udf tm-inlet, set the outlet boundary condition to Pressure-outlet, and set the One side wall is set as the convective heat transfer surface, UDF defines the comprehensive temperature value of the furnace wall air per hour, the convective heat transfer coefficient is set to the preset value according to the actual requirements, the other walls are set as the adiabatic wall surface, and the contact surface between the tow and the gas is set is Coupled;
(4.7)、在Dynamic Mesh选项，并勾选激活动网格选项，勾选Mesh Methods选项中的Smoothing，Remeshing选项；在Dynamic Mesh Zones中定义运动区域，在Zone Names下选择rigid wall，在Type中选择Rigid Body，在Preview Mesh Motion中设置预览网格运动；(4.7), in the Dynamic Mesh option, and check the active mesh option, check the Smoothing and Remeshing options in the Mesh Methods option; define the motion area in Dynamic Mesh Zones, select rigid wall under Zone Names, and in Type Select Rigid Body and set preview mesh motion in Preview Mesh Motion;
(4.8)、选择Check Case后进行计算，计算时间根据实际工程中参数设定；其中选择Check Case后进行计算所使用的高温碳化炉全流场计算的三维数学模型，包括有连续性方程、动量方程和能量方程，分别如下所示：(4.8) Calculate after selecting Check Case, and the calculation time is set according to the parameters in the actual project; among them, the three-dimensional mathematical model of the full flow field calculation of the high-temperature carbonization furnace used for calculation after selecting Check Case, including continuity equation, momentum The equation and the energy equation, respectively, are as follows:
Figure PCTCN2020129173-appb-000006
Figure PCTCN2020129173-appb-000006
式中，ρ-流体密度；t-时间；V-速度矢量，其中u、v、w为V在x、y及z三个方向上的分量。In the formula, ρ-fluid density; t-time; V-velocity vector, where u, v, w are the components of V in the three directions of x, y and z.
动量方程的Navier-Stokes方程：Navier-Stokes equation for momentum equation:
Figure PCTCN2020129173-appb-000007
Figure PCTCN2020129173-appb-000007
其中，μ是动力黏度，F b是微元上的体积力； Among them, μ is the dynamic viscosity, F b is the volume force on the micro-element;
能量守恒方程：Energy conservation equation:
Figure PCTCN2020129173-appb-000008
Figure PCTCN2020129173-appb-000008
其中，C p—比热容，T—温度，k—流体传热系数，S T—粘性耗散项； Among them, C p —specific heat capacity, T—temperature, k—fluid heat transfer coefficient, S T —viscous dissipation term;
具体实施过程中，优选选择高温碳化炉马弗腔体流体计算域三维仿真模型的中心点作为检测面，检测面为过中心点的X方向平面，选择丝束表面作为检测面，检测面为过中心点的X方向平面。In the specific implementation process, it is preferable to select the center point of the three-dimensional simulation model of the fluid computational domain of the muffle cavity of the high-temperature carbonization furnace as the detection surface, the detection surface is the X-direction plane passing through the center point, and the surface of the tow is selected as the detection surface, and the detection surface is the detection surface. The X-direction plane of the center point.
(5)、在相同设置条件下，通过将高温碳化炉马弗腔体流体计算域和被加热丝束的三维仿真模型设置不同参数，并重复步骤(1)-(4)，以进行多次模拟计算，由检测面的温度变化云图，利用后处理软件CFD-POST实现丝束加热过程中炉腔内温度场、丝束表面温度场的可视性，直观输出的结果。(5) Under the same setting conditions, set different parameters by setting the fluid computational domain of the muffle cavity of the high temperature carbonization furnace and the three-dimensional simulation model of the heated tow, and repeat steps (1)-(4) for several times. For simulation calculation, from the temperature change cloud map of the detection surface, the post-processing software CFD-POST is used to realize the visibility of the temperature field in the furnace cavity and the surface temperature field of the tow during the heating process of the tow, and the results are intuitively output.
利用后处理软件CFD-POST实现丝束加热过程中炉腔内温度场、丝束表面温度场的可视性如图4、图5。由图可见三维计算结果不仅精度完全可以满足要求，而且输出数据的类型多样，输出的结果也更加直观。与现有技术相比，本发明提供的一种基于动网格方法的高温碳化炉中丝束加热性能三维模拟方法，具有如下优势：The post-processing software CFD-POST is used to realize the visibility of the temperature field in the furnace cavity and the surface temperature field of the tow during the tow heating process, as shown in Figure 4 and Figure 5. It can be seen from the figure that not only the accuracy of the three-dimensional calculation results can fully meet the requirements, but also the types of output data are diverse, and the output results are more intuitive. Compared with the prior art, the three-dimensional simulation method for tow heating performance in a high-temperature carbonization furnace based on the dynamic grid method provided by the present invention has the following advantages:
(1)通过在CFD计算软件中引入丝束在高温碳化炉中的运动状态，更接近实际。(1) By introducing the motion state of the tow in the high-temperature carbonization furnace into the CFD calculation software, it is closer to reality.
(2)关于高温碳化炉丝束加热性能的问题，传统的方法通常都是对加热后的丝束性能进行检测。本发明创造性地利用动网格方法模拟丝束的运动状态，通过丝束表面温度变化云图判定炉膛的加热能力，从而可以更好的衡量结构设计。(2) Regarding the heating performance of the tow in the high temperature carbonization furnace, the traditional method is usually to test the performance of the tow after heating. The invention creatively uses the moving grid method to simulate the motion state of the tow, and judges the heating capacity of the furnace through the cloud image of the temperature change on the surface of the tow, so that the structure design can be better measured.
(3)可以直观动态的计算丝束表面在高温碳化炉内运动过程中任意时刻的温度分布特性。(3) The temperature distribution characteristics of the tow surface at any time during the movement of the tow surface in the high temperature carbonization furnace can be intuitively and dynamically calculated.
(4)可以进一步用来研究高温碳化炉对丝束的加热性能，从而为高温碳化炉中在设计时提供参考。(4) It can be further used to study the heating performance of the high temperature carbonization furnace to the tow, so as to provide a reference for the design of the high temperature carbonization furnace.
本发明通过UDF定义丝束运动规律，入口边界velocity-inlet中修改气流速度大小可以得出不同丝束运动状态下丝束表面温度分布与马弗炉内蓄热能力的变化规律。分析丝束表面温度分布，丝束温度随着进入炉腔内长度逐渐增加，炉腔温度时刻保持恒定。综上，丝束在进入炉腔的过程中表面温度逐渐升高，丝束被加热，炉膛壁面具有良好的加热效果，其蓄热性能优越；如图3所示，从任意时刻丝束在炉腔内加热过程的温度分布示意图和如图4所示丝束表面的温度分布示意图显示出丝束表面被加热的过程温度呈现梯度增加，丝束表面温度分布均匀；如图5所示，本发明炉腔内温度分布示意图显示马弗腔内温度分布均匀，说明该炉腔内的气流组织合理分布。为验证仿真结果，应多次模拟，比较分析实验结果，得出丝束运动速度与马弗炉腔蓄热性能的最佳方案。The invention defines the motion law of the tow through UDF, and by modifying the airflow velocity in the inlet boundary velocity-inlet, the changing law of the temperature distribution on the surface of the tow and the heat storage capacity in the muffle furnace under different tow motion states can be obtained. The temperature distribution on the surface of the tow is analyzed. The temperature of the tow increases gradually with the length of the tow entering the furnace cavity, and the temperature of the furnace cavity remains constant at all times. To sum up, the surface temperature of the tow gradually increases during the process of entering the furnace cavity, the tow is heated, the furnace wall has a good heating effect, and its heat storage performance is superior; as shown in Figure 3, from any time the tow is in the furnace The schematic diagram of the temperature distribution in the heating process in the cavity and the schematic diagram of the temperature distribution on the surface of the tow as shown in Figure 4 show that the temperature of the tow surface is heated in a gradient increase, and the temperature distribution on the surface of the tow is uniform; as shown in Figure 5, the present invention The schematic diagram of the temperature distribution in the furnace cavity shows that the temperature distribution in the muffle cavity is uniform, indicating that the air flow in the furnace cavity is reasonably distributed. In order to verify the simulation results, multiple simulations should be performed, and the experimental results should be compared and analyzed to obtain the best solution for the speed of the tow movement and the heat storage performance of the muffle furnace cavity.
本发明所述的实施例仅仅是对本发明的优选实施方式进行的描述，并非对本发明构思和范围进行限定，在不脱离本发明设计思想的前提下，本领域中工程技术人员对本发明的技术方案作出的各种变型和改进，均应落入本发明的保护范围，本发明请求保护的技术内容，已经全部记载在权利要求书中。The embodiments of the present invention are only descriptions of the preferred embodiments of the present invention, and do not limit the concept and scope of the present invention. Various modifications and improvements made should fall within the protection scope of the present invention, and the technical content claimed in the present invention has been fully recorded in the claims.

Claims (3)

  1. 一种基于动网格方法的高温碳化炉中丝束加热性能三维模拟方法，其特征在于所述方法包括有如下步骤：A three-dimensional simulation method for tow heating performance in a high-temperature carbonization furnace based on a dynamic grid method, characterized in that the method comprises the following steps:
    (1)、采用三维CAD软件SOLIDWORKS建立高温碳化炉马弗腔体流体计算域和被加热丝束的三维仿真模型；(1) Using the 3D CAD software SOLIDWORKS to establish the 3D simulation model of the fluid computational domain of the muffle cavity of the high temperature carbonization furnace and the heated tow;
    (2)、将步骤(1)建立的高温碳化炉马弗腔体流体计算域和被加热丝束的三维仿真模型分别传递到ICEM软件的Global Mesh Setup模块中，在Global Mesh Setup模块中Volume Meshing Parameters设置网格划分类型为Tetra/Mixed，对高温碳化炉马弗腔体流体计算域和被加热丝束的三维仿真模型进行网格划分，整个计算域采用非结构网格，对靠近丝束表面壁面处进行网格加密，根据ICEM软件中网格质量的判断标准，保证整体结构的网格质量大于0.5，定义所有三维仿真模型的进出口与壁面边界名称，包括丝束壁面、炉腔进出口和炉腔壁面边界名称；(2), transfer the fluid computational domain of the muffle cavity of the high-temperature carbonization furnace and the three-dimensional simulation model of the heated tow established in step (1) to the Global Mesh Setup module of the ICEM software, respectively. In the Global Mesh Setup module, Volume Meshing Parameters Set the meshing type to Tetra/Mixed, and mesh the fluid computational domain of the muffle cavity of the high-temperature carbonization furnace and the 3D simulation model of the heated tow. The entire computational domain uses unstructured meshes. The mesh is refined at the wall surface. According to the grid quality judgment standard in the ICEM software, the mesh quality of the overall structure is guaranteed to be greater than 0.5, and the inlet and outlet and wall boundary names of all 3D simulation models are defined, including the tow wall surface and the inlet and outlet of the furnace cavity. and the name of the wall boundary of the furnace cavity;
    (3)、将划分好网格的流体计算域和被加热丝束的三维仿真模型传递到ANSYS软件的FLUENT模块，并进行设置边界条件；(3) Transfer the meshed fluid computational domain and the three-dimensional simulation model of the heated tow to the FLUENT module of the ANSYS software, and set the boundary conditions;
    (4)、在ANSYS软件中的FLUENT模块里设置温度检测面，并进行仿真运算得到结果，以此作为判定高温碳化炉中丝束加热性能的指标；ANSYS软件中的FLUENT模块进行设置的过程如下：(4) Set the temperature detection surface in the FLUENT module in the ANSYS software, and perform simulation operations to obtain the results, which are used as indicators to determine the heating performance of the tow in the high-temperature carbonization furnace; the process of setting the FLUENT module in the ANSYS software is as follows :
    (4.1)、在User Defined选项导入根据工艺参数编制的自定义丝束运动参数，通过UDF控制被加热丝束的运动状态，实现丝束在高温碳化炉中运动过程；(4.1) Import the custom tow motion parameters compiled according to the process parameters in the User Defined option, and control the motion state of the heated tow through UDF to realize the motion process of the tow in the high temperature carbonization furnace;
    (4.2)、在General选项中，将y方向Gravitational Acceleration根据实际设定为预设值，time选项设置为Transient瞬态传热；(4.2) In the General option, set the Gravitational Acceleration in the y direction to the preset value according to the actual situation, and set the time option to Transient transient heat transfer;
    (4.3)、将Models选项中的Energy勾选Energy Equation，Viscous Models选项中选取laminar模型，为了判断炉腔内气流运动状态，引入雷诺数进行描述，雷诺数的计算公式为：(4.3) Check the Energy Equation in the Models option, and select the laminar model in the Viscous Models option. In order to judge the airflow movement state in the furnace cavity, the Reynolds number is introduced for description. The calculation formula of the Reynolds number is:
    Figure PCTCN2020129173-appb-100001
    Figure PCTCN2020129173-appb-100001
    其中，其中v、ρ、μ分别为流体的流速、密度与黏性系数，d为特征长度；通过雷诺数的计算，进而选择湍流模型为laminar模型；Among them, v, ρ, μ are the flow velocity, density and viscosity coefficient of the fluid, respectively, and d is the characteristic length; through the calculation of the Reynolds number, the turbulence model is selected as the laminar model;
    (4.4)在Materials Fluid选项部分选择空气和氮气；在Materials Solid选项建立丝束的物理参数，包括密度、比热容和热传导率参数；(4.4) Select air and nitrogen in the Materials Fluid option section; establish the physical parameters of the tow in the Materials Solid option, including density, specific heat capacity and thermal conductivity parameters;
    (4.5)、在Cell Zone Conditions选项中将Fluid1部分设为氮气，Fluid2部分设为氧气；将Solid1部分设为丝束；(4.5), in the Cell Zone Conditions option, set the Fluid1 part to nitrogen, and the Fluid2 part to oxygen; set the Solid1 part to tow;
    (4.6)、在Boundary Conditions选项中设置入口边界条件为Pressure-inlet，并将Velocity Magnitude根据实际要求设置为预设值，Thermal选项设置为udftm-inlet，设置出口边界条件为Pressure-outlet，将一侧壁面设置为对流换热面，UDF定义每小时内炉壁空气综合温度值，对流换热系数根据实际要求设置为预设值，其他壁面设置为绝热壁面，将丝束与气体接触面设置为Coupled；(4.6) In the Boundary Conditions option, set the inlet boundary condition to Pressure-inlet, set the Velocity Magnitude to the preset value according to the actual requirements, set the Thermal option to udftm-inlet, set the outlet boundary condition to Pressure-outlet, and set a The side wall surface is set as the convective heat transfer surface, UDF defines the comprehensive temperature value of the furnace wall air per hour, the convective heat transfer coefficient is set to a preset value according to the actual requirements, other walls are set as adiabatic wall surfaces, and the contact surface between the tow and the gas is set as Coupled;
    (4.7)、在Dynamic Mesh选项，并勾选激活动网格选项，勾选Mesh Methods选项中的Smoothing，Remeshing选项；在Dynamic Mesh Zones中定义运动区域，在Zone Names下选择rigid wall，在Type中选择Rigid Body，在Preview Mesh Motion中设置预览网格运动；(4.7), in the Dynamic Mesh option, and check the active mesh option, check the Smoothing and Remeshing options in the Mesh Methods option; define the motion area in Dynamic Mesh Zones, select rigid wall under Zone Names, and in Type Select Rigid Body and set preview mesh motion in Preview Mesh Motion;
    (4.8)、选择Check Case后进行计算，计算时间根据实际工程中参数设定；其中选择Check Case后进行计算所使用的高温碳化炉全流场计算的三维数学模型，包括有连续性方程、动量方程和能量方程，分别如下所示：(4.8) Calculate after selecting Check Case, and the calculation time is set according to the parameters in the actual project; among them, the three-dimensional mathematical model of the full flow field calculation of the high-temperature carbonization furnace used for calculation after selecting Check Case, including continuity equation, momentum The equation and the energy equation, respectively, are as follows:
    Figure PCTCN2020129173-appb-100002
    Figure PCTCN2020129173-appb-100002
    式中，ρ-流体密度；t-时间；V-速度矢量，其中u、v、w为V在x、y及z三个方向上的分量。In the formula, ρ-fluid density; t-time; V-velocity vector, where u, v, w are the components of V in the three directions of x, y and z.
    动量方程的Navier-Stokes方程：Navier-Stokes equation for momentum equation:
    Figure PCTCN2020129173-appb-100003
    Figure PCTCN2020129173-appb-100003
    其中，μ是动力黏度，F b是微元上的体积力； Among them, μ is the dynamic viscosity, F b is the volume force on the micro-element;
    能量守恒方程：Energy conservation equation:
    Figure PCTCN2020129173-appb-100004
    Figure PCTCN2020129173-appb-100004
    其中，C p—比热容，T—温度，k—流体传热系数，S T—粘性耗散项； Among them, C p —specific heat capacity, T—temperature, k—fluid heat transfer coefficient, S T —viscous dissipation term;
    (5)、在相同设置条件下，通过将高温碳化炉马弗腔体流体计算域和被加热丝束的三维仿真模型设置不同参数，并重复步骤(1)-(4)，以进行多次模拟计算，由检测面的温度变化云图，利用后处理软件CFD-POST实现丝束加热过程 中炉腔内温度场、丝束表面温度场的可视性，直观输出的结果。(5) Under the same setting conditions, set different parameters by setting the fluid computational domain of the muffle cavity of the high temperature carbonization furnace and the three-dimensional simulation model of the heated tow, and repeat steps (1)-(4) for several times. For simulation calculation, from the temperature change cloud map of the detection surface, the post-processing software CFD-POST is used to realize the visibility of the temperature field in the furnace cavity and the surface temperature field of the tow during the heating process of the tow, and the results are intuitively output.
  2. 根据权利要求1所述的一种基于动网格方法的高温碳化炉中丝束加热性能三维模拟方法，其特征在于：步骤(1)中对高温碳化炉马弗腔体流体计算域和被加热丝束的三维仿真模型设定的参数包括：马弗腔体几何形状和几何尺寸、进出口密封几何形状和几何尺寸、进出口密封氮气管进口尺寸、氮气管出口尺寸、被加热丝束的几何形状和几何尺寸。A three-dimensional simulation method of tow heating performance in a high temperature carbonization furnace based on a dynamic grid method according to claim 1, characterized in that: in step (1), the fluid calculation domain and the heated temperature of the muffle cavity of the high temperature carbonization furnace are calculated The parameters set by the three-dimensional simulation model of the tow include: muffle cavity geometry and geometry, inlet and outlet seal geometry and geometry, inlet and outlet seal nitrogen pipe inlet size, nitrogen pipe outlet size, and the geometry of the heated tow shape and geometry.
  3. 根据权利要求1所述的一种基于动网格方法的高温碳化炉中丝束加热性能三维模拟方法，其特征在于：步骤(4)中，选择高温碳化炉马弗腔体流体计算域三维仿真模型的中心点作为检测面，检测面为过中心点的X方向平面，选择丝束表面作为检测面，检测面为过中心点的X方向平面。A three-dimensional simulation method for tow heating performance in a high-temperature carbonization furnace based on a dynamic grid method according to claim 1, characterized in that: in step (4), a three-dimensional simulation of the fluid computational domain of the muffle cavity of the high-temperature carbonization furnace is selected The center point of the model is used as the detection surface, and the detection surface is the X direction plane passing through the center point. The surface of the tow is selected as the detection surface, and the detection surface is the X direction plane passing through the center point.
PCT/CN2020/129173 2020-11-16 2020-11-16 Dynamic mesh method-based method for tow heating performance three-dimensional simulation in high temperature carbonization furnace WO2022099714A1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
PCT/CN2020/129173 WO2022099714A1 (en) 2020-11-16 2020-11-16 Dynamic mesh method-based method for tow heating performance three-dimensional simulation in high temperature carbonization furnace

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/CN2020/129173 WO2022099714A1 (en) 2020-11-16 2020-11-16 Dynamic mesh method-based method for tow heating performance three-dimensional simulation in high temperature carbonization furnace

Publications (1)

Publication Number Publication Date
WO2022099714A1 true WO2022099714A1 (en) 2022-05-19

Family

ID=81602082

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2020/129173 WO2022099714A1 (en) 2020-11-16 2020-11-16 Dynamic mesh method-based method for tow heating performance three-dimensional simulation in high temperature carbonization furnace

Country Status (1)

Country Link
WO (1) WO2022099714A1 (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116502564A (en) * 2023-06-27 2023-07-28 江铃汽车股份有限公司 Parameter acquisition method, system and equipment for face wind sensation evaluation
CN116822021A (en) * 2023-06-28 2023-09-29 东方空间技术(山东)有限公司 Determination method and device for diversion trench parameters of trapezoid launching pad and computing equipment
CN116894372A (en) * 2023-09-11 2023-10-17 潍柴动力股份有限公司 Temperature prediction method and device for frozen sand mold based on variable temperature environment
CN117113777A (en) * 2023-09-08 2023-11-24 武汉理工大学 Drill string vortex-induced vibration calculation method considering internal flow

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100661054B1 (en) * 2005-09-07 2006-12-26 요업기술원 Apparatus for designing a furnace and method therefor
CN201099715Y (en) * 2007-07-05 2008-08-13 聚川企业股份有限公司 Perpendicular wind path type carbon fibre oxidation oven
CN102051711A (en) * 2009-11-17 2011-05-11 安徽鑫丰碳纤维有限公司 Process for producing polyacrylonitrile based carbon fibers
CN111400958A (en) * 2020-04-15 2020-07-10 天华化工机械及自动化研究设计院有限公司 ANSYS-based simulation method for oxygen content distribution of high-temperature carbonization furnace
CN111651908A (en) * 2020-07-17 2020-09-11 天华化工机械及自动化研究设计院有限公司 ANSYS-based multi-field coupling stress distribution simulation method for high-temperature carbonization furnace
CN111680432A (en) * 2020-07-17 2020-09-18 天华化工机械及自动化研究设计院有限公司 Low-temperature carbonization furnace multi-coupling field stress distribution simulation method based on WORKBENCH

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100661054B1 (en) * 2005-09-07 2006-12-26 요업기술원 Apparatus for designing a furnace and method therefor
CN201099715Y (en) * 2007-07-05 2008-08-13 聚川企业股份有限公司 Perpendicular wind path type carbon fibre oxidation oven
CN102051711A (en) * 2009-11-17 2011-05-11 安徽鑫丰碳纤维有限公司 Process for producing polyacrylonitrile based carbon fibers
CN111400958A (en) * 2020-04-15 2020-07-10 天华化工机械及自动化研究设计院有限公司 ANSYS-based simulation method for oxygen content distribution of high-temperature carbonization furnace
CN111651908A (en) * 2020-07-17 2020-09-11 天华化工机械及自动化研究设计院有限公司 ANSYS-based multi-field coupling stress distribution simulation method for high-temperature carbonization furnace
CN111680432A (en) * 2020-07-17 2020-09-18 天华化工机械及自动化研究设计院有限公司 Low-temperature carbonization furnace multi-coupling field stress distribution simulation method based on WORKBENCH

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
OUYANG YANYAN, YAXUAN YAO, ZHENG LIU, HUI CHEN, LINGLING ＲEN: "Overview of Ｒesearch Progresses in High-Temperature Carbonization Process of Carbon Fibers", METROLOGY SCIENCE AND TECHNOLOGY, vol. 10, 31 October 2020 (2020-10-31), pages 68 - 72, XP055930294, ISSN: 2096-9015, DOI: 10.3969/j.issn.2096-9015.2020.10.17 *
ZHANG HUANLEI, CHENG LI-SHENG;ZHANG ZHENG-HE;TAN JING;YANG WEI-MIN;AN YING: "The Simulation Study on the Cooperative Graphitization of Carbon Fibers via Microwave and Laser", INFRARED, vol. 41, no. 2, 29 February 2020 (2020-02-29), pages 37 - 43, XP055930296, ISSN: 1672-8785, DOI: 10.3969/j.issn.1672-8785.2020.02.006 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116502564A (en) * 2023-06-27 2023-07-28 江铃汽车股份有限公司 Parameter acquisition method, system and equipment for face wind sensation evaluation
CN116502564B (en) * 2023-06-27 2023-10-31 江铃汽车股份有限公司 Parameter acquisition method, system and equipment for face wind sensation evaluation
CN116822021A (en) * 2023-06-28 2023-09-29 东方空间技术(山东)有限公司 Determination method and device for diversion trench parameters of trapezoid launching pad and computing equipment
CN116822021B (en) * 2023-06-28 2024-02-02 东方空间技术(山东)有限公司 Determination method and device for diversion trench parameters of trapezoid launching pad and computing equipment
CN117113777A (en) * 2023-09-08 2023-11-24 武汉理工大学 Drill string vortex-induced vibration calculation method considering internal flow
CN117113777B (en) * 2023-09-08 2024-04-09 武汉理工大学 Drill string vortex-induced vibration calculation method considering internal flow
CN116894372A (en) * 2023-09-11 2023-10-17 潍柴动力股份有限公司 Temperature prediction method and device for frozen sand mold based on variable temperature environment
CN116894372B (en) * 2023-09-11 2023-12-15 潍柴动力股份有限公司 Temperature prediction method and device for frozen sand mold based on variable temperature environment

Similar Documents

Publication Publication Date Title
WO2022099714A1 (en) Dynamic mesh method-based method for tow heating performance three-dimensional simulation in high temperature carbonization furnace
WO2022099713A1 (en) Three-dimensional simulation method for tow heating process in low temperature carbonization furnace based on overset model
Afshari et al. On numerical methods; optimization of CFD solution to evaluate fluid flow around a sample object at low Re numbers
CN111400958A (en) ANSYS-based simulation method for oxygen content distribution of high-temperature carbonization furnace
CN112446177B (en) Simulation method for heat insulation performance of external heat insulation material of high-temperature carbonization furnace
CN111651908B (en) ANSYS-based multi-field coupling stress distribution simulation method for high-temperature carbonization furnace
Mohotti et al. Comparison of Reynolds Averaging Navier-Stokes (RANS) turbulent models in predicting wind pressure on tall buildings
Zhang et al. Numerical and experimental study on hydraulic performance of emitters with arc labyrinth channels
Liu et al. An extensive comparison of modified zero-equation, standard k-ε, and LES models in predicting urban airflow
WO2022099715A1 (en) Simulation method for visualizing dynamic distribution characteristics of oxygen concentration in high-temperature carbonization furnace
CN113177370B (en) Wind flow-dust gas-solid two-phase flow numerical simulation method considering environmental humidity factor
WO2022011723A1 (en) Ansys-based multi-field coupling stress distribution simulation method for high temperature carbonization furnace
CN112528572B (en) Low-temperature carbonization furnace tow heating process three-dimensional simulation method based on OVERSET model
WO2022160850A1 (en) Method and apparatus for determining fire resistance data, and electronic device
CN109829210A (en) A kind of S type coefficient of Pitot tube scaling method calculated based on CFD numerical value
Care et al. On the impact of anemometer size on the velocity field in a closed wind tunnel
Xu et al. A review of simplified numerical approaches for fast urban airflow simulation
WO2021207952A1 (en) High-temperature carbonization furnace oxygen content distribution simulation method based on ansys
Hernández-Castillo et al. Heat transfer by natural convection and radiation in three dimensional differentially heated tall cavities
CN105808860A (en) Analyzing method based on computational fluid dynamics
CN109657368B (en) Optimization method of air conditioner radiator
WO2022099712A1 (en) Simulation method for heating performance of graphite rod in high-temperature carbonization furnace
CN113536640B (en) Optimization design method for internal flow channel structure of air distributor based on orthogonal test
CN112270109B (en) Method for simulating heating performance of graphite rod in high-temperature carbonization furnace
WO2022099716A1 (en) Method for simulating thermal insulation property of external thermal insulation material of high-temperature carbonization furnace

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 20961268

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 20961268

Country of ref document: EP

Kind code of ref document: A1