CN112528572B - Low-temperature carbonization furnace tow heating process three-dimensional simulation method based on OVERSET model - Google Patents
Low-temperature carbonization furnace tow heating process three-dimensional simulation method based on OVERSET model Download PDFInfo
- Publication number
- CN112528572B CN112528572B CN202011275956.4A CN202011275956A CN112528572B CN 112528572 B CN112528572 B CN 112528572B CN 202011275956 A CN202011275956 A CN 202011275956A CN 112528572 B CN112528572 B CN 112528572B
- Authority
- CN
- China
- Prior art keywords
- setting
- model
- low
- carbonization furnace
- temperature carbonization
- 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
Links
- 238000004088 simulation Methods 0.000 title claims abstract description 50
- 238000003763 carbonization Methods 0.000 title claims abstract description 47
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000010438 heat treatment Methods 0.000 title claims abstract description 31
- 238000004364 calculation method Methods 0.000 claims abstract description 47
- 239000012530 fluid Substances 0.000 claims abstract description 31
- 230000033001 locomotion Effects 0.000 claims abstract description 24
- 238000009826 distribution Methods 0.000 claims abstract description 17
- 238000013178 mathematical model Methods 0.000 claims abstract description 10
- 238000012800 visualization Methods 0.000 claims abstract description 7
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 claims description 24
- 230000008569 process Effects 0.000 claims description 15
- 238000001514 detection method Methods 0.000 claims description 13
- 229910052757 nitrogen Inorganic materials 0.000 claims description 9
- 230000000903 blocking effect Effects 0.000 claims description 6
- 229910001873 dinitrogen Inorganic materials 0.000 claims description 6
- 238000004134 energy conservation Methods 0.000 claims description 6
- 239000000463 material Substances 0.000 claims description 6
- 238000007781 pre-processing Methods 0.000 claims description 6
- 238000012546 transfer Methods 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 5
- 239000007789 gas Substances 0.000 claims description 4
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 3
- 238000009413 insulation Methods 0.000 claims description 3
- 230000003993 interaction Effects 0.000 claims description 3
- 239000001301 oxygen Substances 0.000 claims description 3
- 229910052760 oxygen Inorganic materials 0.000 claims description 3
- 239000007787 solid Substances 0.000 claims description 3
- 230000001052 transient effect Effects 0.000 claims description 3
- 238000009825 accumulation Methods 0.000 claims description 2
- 229920000049 Carbon (fiber) Polymers 0.000 abstract description 18
- 239000004917 carbon fiber Substances 0.000 abstract description 18
- 238000013461 design Methods 0.000 abstract description 9
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 abstract description 9
- 238000007380 fibre production Methods 0.000 abstract description 5
- 238000005265 energy consumption Methods 0.000 abstract description 3
- 238000004458 analytical method Methods 0.000 abstract description 2
- 238000004519 manufacturing process Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 238000005338 heat storage Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 239000000498 cooling water Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 239000012774 insulation material Substances 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000013041 optical simulation Methods 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 229910001220 stainless steel Inorganic materials 0.000 description 1
- 239000010935 stainless steel Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
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
-
- 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)
- Furnace Housings, Linings, Walls, And Ceilings (AREA)
- Waste-Gas Treatment And Other Accessory Devices For Furnaces (AREA)
Abstract
A three-dimensional simulation method for a heating process of a low-temperature carbonization furnace tow based on an OVERSET model relates to the technical field of design analysis methods of low-temperature carbonization furnaces. Constructing a three-dimensional mathematical model of a low-temperature carbonization furnace full flow field; establishing a three-dimensional simulation model of a fluid calculation domain and a heated filament bundle in a muffle cavity of the low-temperature carbonization furnace by using SOLIDWORKS, and setting geometric parameters; carrying out mesh division on the three-dimensional simulation model; transmitting the three-dimensional simulation model to an ANSYS FLUENT calculation module and setting boundary conditions; transmitting the result obtained by the simulation operation to POST-processing software CFD-POST to realize the visualization of physical parameters; different working temperatures and tow movement speeds are set, and simulation calculation is repeatedly carried out for a plurality of times. The distribution state of the surface temperature of the tows at different furnace chamber temperatures and tow movement speeds is judged, so that the heating performance can be better predicted, the experimental cost of furnace body structure design is reduced, and theoretical support is provided for reducing the energy consumption of carbon fiber production.
Description
Technical Field
The invention relates to the technical field of design analysis methods of low-temperature carbonization furnaces used in the production process of carbon fibers.
Background
In the production process of carbon fiber production, a low-temperature carbonization furnace working in an environment of 300-1000 ℃ for a long time is one of key equipment, and mainly comprises a frame body (furnace shell), a heat insulation material, a stainless steel muffle, a muffle counterweight device, a heating element, a sensing element, an inlet nitrogen seal, an outlet cooling water tank, an outlet nitrogen seal, a muffle internal gas detection device, a high-purity nitrogen pipeline system, an electric appliance temperature control system and the like. Through development and development in the last forty years, although the carbon fibers in China have prepared products with the level similar to that of foreign T300-grade carbon fibers, the annual yield and the product performance of the carbon fibers still can not meet the requirements of domestic markets for the carbon fibers. Compared with the international advanced carbon fiber product technical level, the main problems of the domestic carbon fiber are obviously reflected in the poor uniformity and stability of the carbon fiber, and the main reasons are that the carbon fiber process is limited abroad, and the other important influencing factor is that the main production equipment in the carbon fiber production process lags behind abroad. If the current passive situation needs to be changed, the industrial safety threat is broken, the autonomous development of the carbon fiber production line equipment for realizing the localization and industrialization needs to be solved urgently, and professional researchers related to carbon fibers are also needed to continuously improve and research.
The temperature uniformity in the muffle cavity influences the production quality of carbon fibers, the quality of finished carbon fibers is generally detected, factors in the production process are difficult to control, the process cannot be continuously tracked through test equipment, and the quality of the carbon fibers cannot be guaranteed. Therefore, a reasonable design method needs to be selected, the furnace wall can reach the surface temperature meeting the specification, and the unit energy consumption is reduced. The aim of the research is to simulate the heating performance of the tows in the furnace chamber in the design process so as to judge whether the temperature of the surfaces of the tows meets the requirements of the production process, thereby optimally designing the structure of the furnace body and reducing the manufacturing cost of the low-temperature carbonization furnace on the premise of not reducing the existing heating effect.
Disclosure of Invention
In conclusion, the invention aims to reduce the experimental cost of furnace body structure design, provide theoretical support for reducing the energy consumption of carbon fiber production, and provide a basis for related numerical simulation research.
The technical scheme adopted for realizing the purpose of the invention is as follows:
a three-dimensional simulation method for a low-temperature carbonization furnace tow heating process based on an OVERSET model is characterized by comprising the following steps:
(1) constructing a three-dimensional mathematical model of the total flow field of the low-temperature carbonization furnace;
the three-dimensional continuity equation, the momentum equation and the energy conservation equation contained in the three-dimensional mathematical model for calculating the total flow field of the low-temperature carbonization furnace are respectively shown in formulas (1), (2) and (3):
three-dimensional continuity equation:
where ρ -fluid density; t-time; v-velocity vector, where u, V, w are the components of V in the three x, y and z directions;
the momentum equation:
wherein μ is dynamic viscosity, F b Is the volume force on the infinitesimal;
energy conservation equation:
wherein, C p Specific heat capacity, T-temperature, k-coefficient of heat transfer of the fluid, S T -a viscous dissipation term;
(2) according to the geometric parameters of the low-temperature carbonization furnace in the actual engineering, establishing a three-dimensional simulation model of a fluid calculation domain and a heated filament bundle in a muffle cavity of the low-temperature carbonization furnace by using SOLIDWORKS software, and setting related geometric parameters;
(3) respectively transmitting the fluid calculation domain established in the last step and the three-dimensional simulation model of the heated filament bundle to a Blocking module of ICEM software, and performing grid division on the fluid calculation domain and the three-dimensional simulation model of the heated filament bundle in the Blocking module in an O-Block mode; carrying out grid encryption on the position close to the surface wall surface of the filament bundle, ensuring that the number of grid layers is not less than five according to the requirement of a standard wall surface function, adopting a BiGeometric mode for a grid division strategy, setting a control ratio factor to be a default value of 1.2, ensuring that the grid quality of the whole structure is more than 0.85 according to the standard of the grid quality in ICEM software, and simultaneously defining the names of an inlet, an outlet and a wall surface boundary of a fluid calculation domain and a three-dimensional simulation model of the heated filament bundle;
(4) transmitting the fluid calculation domain divided with the grid and the three-dimensional simulation model of the heated filament bundle to a FLUENT calculation module in ANSYS software, and setting boundary conditions; the procedure for setting the FLUENT module in ANSYS software is as follows:
(4.1) introducing custom tow motion parameters compiled according to process parameters into a User Defined option, and controlling the motion state of the heated tows through UDF to realize the motion process of the tows in a low-temperature carbonization furnace;
(4.2) importing a background grid and a foreground grid which are required by calculation and are divided by ICEM software in the File option;
(4.3) in the General option, setting the y-direction probabilistic accumulation option as a preset value according to the actual situation, and setting the time option as Transient heat transfer;
(4.4) selecting an Energy Equation from the Models options, selecting a laminar model from the Viscous Models options, and introducing a Reynolds number for judging the motion state of the airflow in the furnace cavity for description, wherein the Reynolds number has a calculation formula as follows:
wherein v, rho and mu are respectively the flow velocity, density and viscosity coefficient of the fluid, and d is the characteristic length; selecting a turbulence model as a laminar model through the calculation of Reynolds number;
(4.5) selecting air and nitrogen in Materials Fluid option section; establishing physical parameters of the tows in a Materials Solid option, wherein the physical parameters comprise density, specific heat capacity and heat conductivity physical parameters;
(4.6) in the Cell Zone Conditions option, part of Fluid1 is set to nitrogen and part of Fluid2 is set to oxygen; setting Solid1 part as a tow;
(4.7) setting an inlet Boundary condition as Pressure-unlet in the Boundary Conditions option, setting the Velocity magnet as a preset value according to actual requirements, setting the Thermal option as UDF tm-unlet, setting an outlet Boundary condition as Pressure-outlet, setting one side wall surface as a convection heat exchange surface, defining a comprehensive temperature value of furnace wall air in each hour by UDF, setting a convection heat exchange coefficient as a preset value according to actual requirements, setting other wall surfaces as heat insulation wall surfaces, and setting a contact surface between the filament bundle and gas as Coupled;
(4.8) creating an interaction Interface in an Overset Interface option, selecting a Background grid in Background Zones option, selecting a foreground grid in Component Zones, and preprocessing an Overset model in calculation, wherein the preprocessing comprises point searching, hole digging and interpolation relation establishment;
(4.9) after selecting the Check Case, calculating based on the three-dimensional mathematical model of the low-temperature carbonization furnace full flow field in the step (1), and setting the calculation time according to parameters in actual engineering;
(5) transmitting the result obtained by the simulation operation to POST-processing software CFD-POST to realize the visualization of physical parameters in the heating process of the tows in the carbonization furnace;
(6) and under the same setting condition, setting different working temperatures and different tow movement speeds and repeating the steps (2) to (5) to perform multiple times of simulation calculation, judging the distribution state of the surface temperature of the tows at different furnace chamber temperatures and tow movement speeds, and predicting the heating performance of the furnace chamber on the tows, so that the distribution state is used as a basis for designing the muffle cavity structure and the operation process parameters of the low-temperature carbonization furnace.
The related parameters set according to the low-temperature carbonization furnace in the actual engineering in the step (2) comprise: the geometrical shape and the geometrical size of the muffle cavity, the geometrical shape and the geometrical size of the inlet and outlet seal, the inlet size of the inlet and outlet seal nitrogen gas pipe, the size of the outlet of the nitrogen gas pipe, and the geometrical shape and the geometrical size of the tows.
And (3) defining the boundary names of the inlet, the outlet and the wall surfaces of the three-dimensional simulation model in the step (3) to comprise the boundary names of the tow wall surface, the inlet and the outlet of the furnace chamber and the wall surface of the furnace chamber.
In the step (4), the central point of the three-dimensional simulation model of the fluid calculation domain is selected as a detection surface, the detection surface is an X-direction plane passing through the central point, the tow surface is selected as a detection surface, and the detection surface is an X-direction plane passing through the central point.
The simulation result in the step (5) comprises: and (3) detecting a temperature change cloud picture of the surface, and realizing visualization of a temperature field in a furnace cavity and a temperature field on the surface of the tows in the heating process of the tows by using POST-processing software CFD-POST.
The invention has the beneficial effects that: the invention is based on the flow field numerical calculation method, the heating capacity research is carried out on the surface temperature field characteristics of the filament bundle in the hearth under different production process conditions, the defects of the prior art at the design stage of the low-temperature carbonization furnace are overcome, and compared with the prior art, the invention has the following advantages:
(1) the calculation result is close to the actual result by introducing the moving state of the tows in the low-temperature carbonization furnace in CFD calculation software.
(2) Regarding the problem of heating performance of the tow in the low-temperature carbonization furnace, the conventional method is to measure the performance of the heated tow. The invention creatively utilizes the OVERSET model to simulate the movement state of the tows, and judges the heating capacity of the hearth through the cloud chart of the surface temperature change of the tows, thereby being capable of better measuring the structural design.
(3) The temperature distribution characteristics of the surface of the filament bundle at any moment in the process of moving in the low-temperature carbonization furnace can be intuitively and dynamically calculated.
(4) The method can be further used for researching the heating performance of the low-temperature carbonization furnace on the tows, thereby providing a reference for the design of the low-temperature carbonization furnace.
Drawings
FIG. 1 is a schematic diagram of a three-dimensional simulation model established in the simulation method of the present invention.
Fig. 2 is a flow chart of the implementation of the OVERSET model in numerical calculations.
FIG. 3 is a schematic representation of the temperature profile of the tow during heating in the furnace chamber at any point in the present invention.
FIG. 4 is a schematic representation of the temperature profile of the surface of the tow of the present invention.
Fig. 5 is a schematic view showing the temperature distribution in the cavity of the furnace according to the present invention.
Detailed Description
The structure of the present invention will be further described with reference to the accompanying drawings and preferred embodiments of the present invention.
The invention discloses a three-dimensional simulation method for a low-temperature carbonization furnace filament bundle heating process based on an OVERSET model, which comprises the following steps:
(1) constructing a three-dimensional mathematical model of the total flow field of the low-temperature carbonization furnace, and providing a theoretical basis for subsequent simulation calculation;
the three-dimensional continuity equation, the momentum equation and the energy conservation equation contained in the three-dimensional mathematical model for calculating the total flow field of the low-temperature carbonization furnace are respectively shown in formulas (1), (2) and (3):
three-dimensional continuity equation:
where ρ -fluid density; t-time; v-velocity vector, where u, V, w are the components of V in the three x, y and z directions;
the momentum equation:
wherein μ is dynamic viscosity, F b Is the volume force on the infinitesimal;
energy conservation equation:
wherein, C p Specific heat capacity, T-temperature, k-coefficient of heat transfer of the fluid, S T -a viscous dissipation term.
(2) According to the geometric parameters of the low-temperature carbonization furnace in the actual engineering, a three-dimensional simulation model of a fluid calculation domain and a heated filament bundle in a muffle cavity of the low-temperature carbonization furnace is established by adopting three-dimensional aided design software SOLIDWORKS (solid-state optical simulation system), as shown in figure 1, and relevant geometric parameters are set; at least comprises the following steps: the geometrical shape and the geometrical size of the muffle cavity, the geometrical shape and the geometrical size of the inlet and outlet seal, the inlet size of the inlet and outlet seal nitrogen gas pipe, the size of the outlet of the nitrogen gas pipe, and the geometrical shape and the geometrical size of the tows.
(3) Respectively transmitting the fluid calculation domain established in the last step and the three-dimensional simulation model of the heated filament bundle to a Blocking module of ICEM software, and performing grid division on the three-dimensional simulation model in the Blocking module in an O-Block mode; carrying out grid encryption on the position close to the surface wall surface of the filament bundle, ensuring that the number of grid layers is not less than five according to the requirement of a standard wall surface function, adopting a BiGeometric mode for a grid division strategy, setting a control ratio factor to be a default value of 1.2, ensuring that the grid quality of the integral structure is more than 0.85 according to the standard of the grid quality in ICEM software, and simultaneously defining the names of a fluid calculation domain and the inlet and outlet and wall surface boundary of a heated filament bundle three-dimensional simulation model in order to facilitate the later-stage setting of calculation conditions; the method mainly comprises the names of the tow wall surface, the inlet and the outlet of the furnace chamber and the boundary of the furnace chamber wall surface.
(4) Respectively transmitting the fluid calculation domain with the divided grids and the three-dimensional simulation model of the heated filament bundle to an ANSYS FLUENT calculation module and setting boundary conditions as shown in FIG. 1; the method comprises the steps of selecting a central point of a three-dimensional simulation model of a fluid calculation domain as a detection surface, wherein the detection surface is an X-direction plane passing through the central point, selecting a tow surface as the detection surface, and the detection surface is the X-direction plane passing through the central point.
The process of setting the ANSYS FLUENT module is as follows:
(4.1) introducing custom tow motion parameters compiled according to process parameters into a User Defined option, and controlling the motion state of the heated tows through UDF to realize the motion process of the tows in a low-temperature carbonization furnace;
(4.2) importing a background grid and a foreground grid which are required by calculation and are divided by ICEM software in the File option;
(4.3) in the General option, setting the y-direction visual access option as a preset value according to the actual condition, and setting the time option as Transient heat transfer;
(4.4) selecting an Energy Equation from the Models options, selecting a laminar model from the Viscous Models options, and introducing a Reynolds number for judging the motion state of the airflow in the furnace cavity for description, wherein the Reynolds number has a calculation formula as follows:
wherein v, rho and mu are respectively the flow velocity, density and viscosity coefficient of the fluid, and d is the characteristic length; selecting a turbulence model as a laminar model through the calculation of Reynolds number;
(4.5) selecting air and nitrogen in Materials Fluid option section; establishing physical parameters of the tows in a Materials Solid option, wherein the physical parameters comprise density, specific heat capacity and thermal conductivity physical parameters;
(4.6) in the Cell Zone Conditions option, part of Fluid1 is set to nitrogen and part of Fluid2 is set to oxygen; setting Solid1 part as a tow;
(4.7) setting an inlet Boundary condition as Pressure-inlet in a Boundary Conditions option, setting a Velocity map as a preset value according to an actual requirement, setting a Thermal option as UDF tm-inlet, setting an outlet Boundary condition as Pressure-outlet, setting a side wall surface as a convection heat exchange surface, defining a comprehensive temperature value of furnace wall air in each hour by UDF, setting a convection heat exchange coefficient as a preset value according to an actual requirement, setting other wall surfaces as heat insulation wall surfaces, and setting a contact surface of a filament bundle and gas as Coupled;
(4.8) creating an interaction Interface in an Overset Interface option, selecting a Background grid in Background Zones option, selecting a foreground grid in Component Zones, and preprocessing an Overset model in calculation, wherein the preprocessing comprises point searching, hole digging and interpolation relationship establishment, and is shown in FIG. 2;
and (4.9) after selecting the Check Case, calculating based on the three-dimensional mathematical model in the step (1), and setting the calculation time according to parameters in actual engineering.
(5) Transmitting the result obtained by the simulation operation to POST-processing software CFD-POST to realize the visualization of physical parameters in the heating process of the tows in the carbonization furnace; the simulation result comprises: and (3) detecting a temperature change cloud chart of the surface, and realizing visualization of a temperature field in a furnace cavity and a temperature field on the surface of the tows in the heating process of the tows by using POST-processing software CFD-POST (computational fluid dynamics-POST) as shown in figures 4 and 5. The three-dimensional calculation result can completely meet the requirement on precision, the types of output data are various, and the output result is more visual.
(6) Under the same setting condition, different working temperatures and different tow movement speeds are set, and the steps (2) to (5) are repeated to carry out simulation calculation for multiple times, the distribution state of the surface temperature of the tows is judged when the temperatures of the furnace chambers and the tow movement speeds are different, the heating performance of the furnace chambers to the tows can be well predicted, and the distribution state is used as a basis for designing the muffle cavity structure and the operation process parameters of the low-temperature carbonization furnace.
According to the invention, the tow motion rule is defined by udf, and the change rule of tow surface temperature distribution and the heat storage capacity in the muffle furnace in different tow motion states can be obtained by modifying the air flow speed in the inlet boundary velocity-inlet. And (3) analyzing the distribution of the surface temperature of the tows, wherein the temperature of the tows gradually increases along with the length of the tows entering the furnace chamber, and the temperature of the furnace chamber is kept constant all the time.
In conclusion, the surface temperature of the tows gradually rises in the process of entering the furnace chamber, the tows are heated, the wall surface of the furnace chamber has a good heating effect, and the heat storage performance is excellent; the schematic diagram of the temperature distribution of the filament bundle in the heating process in the furnace chamber at any moment is shown in fig. 3, and the schematic diagram of the temperature distribution of the surface of the filament bundle is shown in fig. 4, which shows that the temperature of the heating process of the surface of the filament bundle is increased in a gradient manner, and the temperature distribution of the surface of the filament bundle is uniform; the schematic diagram of the temperature distribution in the furnace chamber is shown in fig. 5, which shows that the temperature distribution in the muffle chamber is uniform, and the reasonable distribution of the airflow organization in the furnace chamber is illustrated. In order to verify the simulation result, simulation is carried out for multiple times, and the experimental result is compared and analyzed to obtain the optimal scheme of the movement speed of the tows and the heat storage performance of the muffle furnace chamber.
The embodiments of the present invention are described only for the preferred embodiments of the present invention, and not for the limitation of the concept and scope of the present invention, and various modifications and improvements made to the technical solution of the present invention by those skilled in the art without departing from the design concept of the present invention shall fall into the protection scope of the present invention, and the technical content of the present invention which is claimed is fully set forth in the claims.
Claims (5)
1. A three-dimensional simulation method for a low-temperature carbonization furnace tow heating process based on an OVERSET model is characterized by comprising the following steps:
(1) constructing a three-dimensional mathematical model of the whole flow field of the low-temperature carbonization furnace;
the three-dimensional continuity equation, the momentum equation and the energy conservation equation contained in the three-dimensional mathematical model for calculating the total flow field of the low-temperature carbonization furnace are respectively shown in formulas (1), (2) and (3):
three-dimensional continuity equation:
where ρ -fluid density; t-time; v-velocity vector, where u, V, w are the components of V in the three x, y and z directions;
the momentum equation:
wherein μ is dynamic viscosity, F b Is the volume force on the infinitesimal;
energy conservation equation:
wherein, C p Specific heat capacity, T-temperature, k-coefficient of heat transfer of the fluid, S T -a viscous dissipation term;
(2) according to the geometric parameters of the low-temperature carbonization furnace in the actual engineering, establishing a three-dimensional simulation model of a fluid calculation domain and a heated filament bundle in a muffle cavity of the low-temperature carbonization furnace by using SOLIDWORKS software, and setting related geometric parameters;
(3) respectively transmitting the fluid calculation domain established in the last step and the three-dimensional simulation model of the heated filament bundle to a Blocking module of ICEM software, and performing grid division on the fluid calculation domain and the three-dimensional simulation model of the heated filament bundle in the Blocking module in an O-Block mode; carrying out grid encryption on the position close to the surface wall surface of the filament bundle, ensuring that the number of grid layers is not less than five according to the requirement of a standard wall surface function, adopting a BiGeometric mode for a grid division strategy, setting a control ratio factor to be a default value of 1.2, ensuring that the grid quality of the whole structure is more than 0.85 according to the standard of the grid quality in ICEM software, and simultaneously defining the names of an inlet, an outlet and a wall surface boundary of a fluid calculation domain and a three-dimensional simulation model of the heated filament bundle;
(4) transmitting the fluid calculation domain with the divided grids and the three-dimensional simulation model of the heated filament bundle to a FLUENT calculation module in ANSYS software, and setting boundary conditions; the procedure for setting the FLUENT module in the ANSYS software is as follows:
(4.1) introducing custom tow motion parameters compiled according to the process parameters into a User Defined option, and controlling the motion state of the heated tow through UDF to realize the motion process of the tow in a low-temperature carbonization furnace;
(4.2) importing a background grid and a foreground grid which are required by calculation and are divided by ICEM software in the File option;
(4.3) in the General option, setting the y-direction probabilistic accumulation option as a preset value according to the actual situation, and setting the time option as Transient heat transfer;
(4.4) selecting an Energy Equation from the Energy options in the Models options, selecting a Laminar model from the Viscous Models options, and introducing a Reynolds number for describing in order to judge the airflow motion state in the furnace cavity, wherein the Reynolds number has a calculation formula as follows:
wherein v, rho and mu are respectively the flow velocity, density and viscosity coefficient of the fluid, and d is the characteristic length; selecting a turbulence model as a laminar model through the calculation of Reynolds number;
(4.5) selecting air and nitrogen in Materials Fluid option section; establishing physical parameters of the tows in a Materials Solid option, wherein the physical parameters comprise density, specific heat capacity and thermal conductivity physical parameters;
(4.6) in the Cell Zone Conditions option, part of Fluid1 is set to nitrogen and part of Fluid2 is set to oxygen; setting Solid1 part as a tow;
(4.7) setting an inlet Boundary condition as Pressure-unlet in the Boundary Conditions option, setting the Velocity magnet as a preset value according to actual requirements, setting the Thermal option as UDF tm-unlet, setting an outlet Boundary condition as Pressure-outlet, setting one side wall surface as a convection heat exchange surface, defining a comprehensive temperature value of furnace wall air in each hour by UDF, setting a convection heat exchange coefficient as a preset value according to actual requirements, setting other wall surfaces as heat insulation wall surfaces, and setting a contact surface between the filament bundle and gas as Coupled;
(4.8) creating an interaction Interface in an Overset Interface option, selecting a Background grid in Background Zones option, selecting a foreground grid in Component Zones, and preprocessing an Overset model in calculation, wherein the preprocessing comprises point searching, hole digging and interpolation relation establishment;
(4.9) after selecting the Check Case, calculating based on the three-dimensional mathematical model of the low-temperature carbonization furnace full flow field in the step (1), and setting the calculation time according to parameters in actual engineering;
(5) transmitting the result obtained by the simulation operation to POST-processing software CFD-POST to realize the visualization of physical parameters in the heating process of the tows in the carbonization furnace;
(6) and under the same setting condition, setting different working temperatures and different tow movement speeds and repeating the steps (2) to (5) to perform multiple times of simulation calculation, judging the distribution state of the surface temperature of the tows at different furnace chamber temperatures and tow movement speeds, and predicting the heating performance of the furnace chamber on the tows, so that the distribution state is used as a basis for designing the muffle cavity structure and the operation process parameters of the low-temperature carbonization furnace.
2. The OVERSET model-based three-dimensional simulation method for the tow heating process of the low-temperature carbonization furnace, according to claim 1, is characterized in that: the relevant parameters set according to the low-temperature carbonization furnace in the practical engineering in the step (2) comprise: the geometrical shape and the geometrical size of the muffle cavity, the geometrical shape and the geometrical size of the inlet and outlet seal, the inlet size of the inlet and outlet seal nitrogen gas pipe, the size of the outlet of the nitrogen gas pipe, and the geometrical shape and the geometrical size of the tows.
3. The OVERSET model-based three-dimensional simulation method for the tow heating process of the low-temperature carbonization furnace, according to claim 1, is characterized in that: and (4) defining the boundary names of the inlet, the outlet and the wall surfaces of the three-dimensional simulation model in the step (3) to comprise the boundary names of the tow wall surface, the inlet, the outlet and the wall surface of the furnace chamber.
4. The OVERSET model-based three-dimensional simulation method for the tow heating process of the low-temperature carbonization furnace, according to claim 1, is characterized in that: in the step (4), the central point of the three-dimensional simulation model of the fluid calculation domain is selected as a detection surface, the detection surface is an X-direction plane passing through the central point, the tow surface is selected as a detection surface, and the detection surface is an X-direction plane passing through the central point.
5. The OVERSET model-based three-dimensional simulation method for the tow heating process of the low-temperature carbonization furnace, according to claim 1, is characterized in that: the simulation result in the step (5) comprises the following steps: and (3) detecting a temperature change cloud picture of the surface, and realizing visualization of a temperature field in a furnace cavity and a temperature field on the surface of the tows in the heating process of the tows by using POST-processing software CFD-POST.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011275956.4A CN112528572B (en) | 2020-11-16 | 2020-11-16 | Low-temperature carbonization furnace tow heating process three-dimensional simulation method based on OVERSET model |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011275956.4A CN112528572B (en) | 2020-11-16 | 2020-11-16 | Low-temperature carbonization furnace tow heating process three-dimensional simulation method based on OVERSET model |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112528572A CN112528572A (en) | 2021-03-19 |
| CN112528572B true CN112528572B (en) | 2022-09-16 |
Family
ID=74982531
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011275956.4A Active CN112528572B (en) | 2020-11-16 | 2020-11-16 | Low-temperature carbonization furnace tow heating process three-dimensional simulation method based on OVERSET model |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112528572B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112800700B (en) * | 2021-04-13 | 2021-06-25 | 中国空气动力研究与发展中心计算空气动力研究所 | Low-temperature surface dry mode frosting simulation method and device, electronic equipment and medium |
| CN113268824B (en) * | 2021-04-23 | 2023-07-11 | 上海卫星工程研究所 | Fluid loop numerical simulation method, system and medium based on flow field in temperature control valve |
| CN113128097B (en) * | 2021-04-29 | 2023-11-17 | 浙江理工大学 | Method for simulating and predicting heat transfer performance of porous nanofiber medium |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105930585A (en) * | 2016-04-21 | 2016-09-07 | 厦门大学 | CFD-based simulation method for flow field and temperature field of Shell gasifier |
| CN108363857A (en) * | 2018-02-05 | 2018-08-03 | 燕山大学 | Recuperative heater flow field and inside workpiece temperature, thermal-stress analysis method |
| CN111651908A (en) * | 2020-07-17 | 2020-09-11 | 天华化工机械及自动化研究设计院有限公司 | ANSYS-based multi-field coupling stress distribution simulation method for high-temperature carbonization furnace |
| CN111931363A (en) * | 2020-07-28 | 2020-11-13 | 南京钢铁股份有限公司 | Motion simulation method for plate blank of roller-hearth heating furnace |
-
2020
- 2020-11-16 CN CN202011275956.4A patent/CN112528572B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105930585A (en) * | 2016-04-21 | 2016-09-07 | 厦门大学 | CFD-based simulation method for flow field and temperature field of Shell gasifier |
| CN108363857A (en) * | 2018-02-05 | 2018-08-03 | 燕山大学 | Recuperative heater flow field and inside workpiece temperature, thermal-stress analysis method |
| CN111651908A (en) * | 2020-07-17 | 2020-09-11 | 天华化工机械及自动化研究设计院有限公司 | ANSYS-based multi-field coupling stress distribution simulation method for high-temperature carbonization furnace |
| CN111931363A (en) * | 2020-07-28 | 2020-11-13 | 南京钢铁股份有限公司 | Motion simulation method for plate blank of roller-hearth heating furnace |
Non-Patent Citations (4)
| Title |
|---|
| Fluent对碳纤维预氧化炉优化设计的可行性分析;闫亮等;《建材世界》;20090215(第01期);97-100页 * |
| 基于ANSYS CFX的反应炉燃烧传热数值模拟研究;闫媛媛等;《陕西科技大学学报(自然科学版)》;20131025(第05期);134-138页 * |
| 基于FLUENT软件的热压氧化高压釜流场数值模拟;吴国玉等;《黄金》;20201015(第10期);64-67页 * |
| 碳纤维材料高温碳化过程研究进展概述;欧阳艳艳等;《计量科学与技术》;20201031;68-72页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112528572A (en) | 2021-03-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112528572B (en) | Low-temperature carbonization furnace tow heating process three-dimensional simulation method based on OVERSET model | |
| CN111400958B (en) | ANSYS-based simulation method for oxygen content distribution of high-temperature carbonization furnace | |
| WO2022099713A1 (en) | Three-dimensional simulation method for tow heating process in low temperature carbonization furnace based on overset model | |
| CN111400934B (en) | Method for simulating oxygen content distribution of low-temperature carbonization furnace based on WORKBENCH | |
| CN111651908B (en) | ANSYS-based multi-field coupling stress distribution simulation method for high-temperature carbonization furnace | |
| CN112446177B (en) | Simulation method for heat insulation performance of external heat insulation material of high-temperature carbonization furnace | |
| CN111680432B (en) | Low-temperature carbonization furnace multi-coupling field stress distribution simulation method based on WORKBENCH | |
| Yang et al. | Experimental and numerical investigations on the thermal performance of a horizontal spiral-coil ground heat exchanger | |
| WO2022011723A1 (en) | Ansys-based multi-field coupling stress distribution simulation method for high temperature carbonization furnace | |
| CN110489835B (en) | Natural ventilation and wall heat storage coupling simulation method based on Ansys software | |
| WO2022099714A1 (en) | Dynamic mesh method-based method for tow heating performance three-dimensional simulation in high temperature carbonization furnace | |
| CN112270109B (en) | Method for simulating heating performance of graphite rod in high-temperature carbonization furnace | |
| WO2021207953A1 (en) | Workbench-based low-temperature carbonization furnace oxygen content distribution simulation method | |
| CN106897520A (en) | A kind of heat transfer system analysis method for reliability containing fuzzy parameter | |
| CN105087882A (en) | Partitioning method for heat treatment stages of vertical quenching furnace | |
| Shen et al. | Soft sensor modeling of blast furnace wall temperature based on temporal–spatial dimensional finite-element extrapolation | |
| CN110737935B (en) | Indoor thermal environment modeling method based on digital twinning | |
| CN108090293A (en) | A kind of blast furnace hearth and bottom corrodes enveloping surface and determines method | |
| WO2021207952A1 (en) | High-temperature carbonization furnace oxygen content distribution simulation method based on ansys | |
| CN110489912B (en) | Method for simulating numerical value of solar cross-season soil heat storage layered slice | |
| CN110728030B (en) | Method for simulating coupling of two sides of steam and water based on numerical value of surface condenser | |
| CN115659875B (en) | Rod bundle fuel assembly test device and test parameter acquisition method thereof | |
| WO2022099712A1 (en) | Simulation method for heating performance of graphite rod in high-temperature carbonization furnace | |
| CN116227060A (en) | Tube-shell heat exchanger structure optimization analysis method based on orthogonal test | |
| CN112435716B (en) | Visual simulation method for dynamic distribution characteristics of oxygen concentration in high-temperature carbonization furnace |
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 | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20240205 Address after: 100728 No. 22 North Main Street, Chaoyang District, Beijing, Chaoyangmen Patentee after: CHINA PETROLEUM & CHEMICAL Corp. Country or region after: China Patentee after: Sinopec Shangai Petrochemical Co.,Ltd. Address before: 730060 Gansu city of Lanzhou province Xigu District Heshui Road No. 3 Patentee before: TIANHUA INSTITUTE OF CHEMICAL MACHINERY AND AUTOMATION Co.,Ltd. Country or region before: China |
|
| TR01 | Transfer of patent right |