Method, system and storage medium for measuring filter rod temperature during suction process based on simulation model

Abstract

The present invention provides a method, system and storage medium for measuring the temperature of a filter rod during a smoking process based on a simulation model, and relates to the technical field of tobacco parameter measurement. The present invention obtains design parameters of a filter rod, and performs collective modeling on the filter rod according to the design parameters; performs CFD simulation analysis on the filter rod to obtain analysis results, performs comprehensive analysis on the analysis results to generate a third temperature evaluation coefficient distribution cloud map, generates a filter rod temperature distribution cloud map according to the third temperature evaluation coefficient distribution cloud map and filter rod material parameters, and obtains the temperature of a position of the filter rod to be measured according to the filter rod temperature distribution cloud map.

Description

Technical Field

[0001] The present invention relates to the technical field of tobacco parameter measurement, and in particular to a method, system and storage medium for measuring filter rod temperature during a smoking process based on a simulation model. Background Art

[0002] Smoking is often considered to be associated with lung cancer, cardiovascular disease and other serious respiratory diseases. When tobacco burns, it produces deadly tar toxins, carbon dioxide and addictive nicotine and other harmful substances. According to studies, the emission of toxic substances in the smoke produced by burning cigarettes is 2-3 times higher than that absorbed by the human body. Even if non-smokers are exposed to this smoke environment for a short time, it will pose a potential threat to their health.

[0003] In traditional filter rod designs, smoke cannot reach a high temperature drop, which may cause the temperature of the inlet smoke to rise, thus adversely affecting the user experience. Accurate measurement of filter rod temperature is crucial to optimizing product performance and improving user experience. However, due to the complex heat transfer and material exchange involved in the smoking process, existing measurement methods are difficult to fully capture these dynamic changes. Traditional methods usually rely on direct measurement or simple model simulation, which shows obvious limitations when dealing with high-speed flow of smoke and mixing of multiple substances. In addition, the physical properties of the filter rod material, such as thermal conductivity and specific heat capacity, have not been fully considered, resulting in inaccurate and unreliable temperature measurement results.

[0004] In the prior art, the publication number CN110793669A discloses a method for quickly determining the temperature of the filter rod during the puffing process, by collecting the basic parameters of the puffing process; establishing a finite element analysis model according to the unit division data and the filter rod parameters, performing unit node numbering, determining the heat exchange boundary and calculating the node coordinates; selecting a turbulence model; treating the flow of smoke and nicotine as a steady flow, wherein the nicotine is regarded as a continuous phase, and simulating under the Euler model; determining the mathematical modulus of the calculation; determining the boundary conditions of the calculation, and calculating the change in the temperature of the filter rod. This prior art changes the traditional method of measuring the actual temperature of the filter rod, and can measure the temperature distribution of the cross section at different positions inside the filter rods of different lengths, and has the characteristics of high calculation accuracy; it can play the role of calculating the actual temperature of the filter rod during actual puffing through the basic parameters of the filter rod. However, the prior art still has defects. When measuring the temperature, the prior art only considers the flow and temperature properties of the smoke, and ignores the thermal properties of the filter rod material and the influence of the condensation phenomenon during the smoke puffing process on the measured temperature. Ignoring the thermal properties and condensation phenomenon of the filter rod material will lead to inaccurate temperature measurement, thereby affecting product performance evaluation and optimization. This error may lead to improper material selection and unreasonable structural design, ultimately affecting user experience.

[0005] The above information disclosed in this Background section is only for enhancement of understanding of the background of the present disclosure and therefore it may contain information that does not constitute the prior art that is already known to one of ordinary skill in the art. Summary of the invention

[0006] The object of the present invention is to provide a method, system and storage medium for measuring the temperature of a filter rod during a suction process based on a simulation model, so as to solve the problems raised in the above-mentioned background technology.

[0007] To achieve the above object, the present invention provides the following technical solutions:

[0008] A method for measuring filter rod temperature during a suction process based on a simulation model, the specific steps comprising:

[0009] Step 1: Obtain the design parameters of the filter rod, and perform collective modeling of the filter rod according to the design parameters;

[0010] Step 2: using ICEM to divide the model into hexahedral grids, measuring the material parameters of the filter rod by experimental means, setting boundary conditions and material parameters, performing CFD simulation analysis on the filter rod, obtaining analysis results, and preprocessing the analysis results, wherein the analysis results include a smoke velocity distribution cloud map of the filter rod, a filter rod pressure distribution cloud map, a smoke temperature distribution cloud map, a tow temperature distribution cloud map, a condensed nicotine distribution cloud map, a nicotine mole fraction distribution cloud map, a condensed glycerol distribution cloud map, a glycerol mole fraction distribution cloud map, a condensed water distribution cloud map, and a water mole fraction distribution cloud map;

[0011] Step 3: Fusing the pre-processed flue gas temperature distribution cloud map and the tow temperature distribution cloud map to obtain a first temperature evaluation coefficient distribution cloud map;

[0012] Step 4: Fusing the pre-processed filter rod pressure distribution cloud map, the smoke velocity distribution cloud map and the first temperature evaluation coefficient distribution cloud map to obtain a second temperature evaluation coefficient distribution cloud map;

[0013] Step 5: the condensed nicotine distribution cloud map, the nicotine mole fraction distribution cloud map, the condensed glycerol distribution cloud map, the glycerol mole fraction distribution cloud map, the condensed water distribution cloud map, and the water mole fraction distribution cloud map are merged to obtain a condensed substance evaluation coefficient cloud map, and the second temperature evaluation coefficient distribution cloud map and the condensed substance evaluation coefficient distribution cloud map are merged to obtain a third temperature evaluation coefficient distribution cloud map;

[0014] Step 6: Generate a filter rod temperature distribution cloud map according to the material parameters and the third temperature evaluation coefficient distribution cloud map, and obtain the temperature of the position to be measured on the filter rod according to the filter rod temperature distribution cloud map.

[0015] Furthermore, the preprocessing operation is: to mark the coordinates of the pixel points of the various distribution cloud maps in the analysis results, with the bottom pixel row being the first row and the leftmost pixel column being the first column, and the pixel point positions of the first recognition image and the second recognition image being mapped one by one through the coordinates, so that each pixel point has a unique coordinate value, and each pixel point position of the various distribution cloud maps in the analysis results is mapped one by one.

[0016] Furthermore, the specific logic for obtaining the first temperature evaluation coefficient distribution cloud map is as follows: the first temperature evaluation coefficient of each coordinate pixel point is calculated according to the flue gas temperature distribution cloud map and the tow temperature distribution cloud map, and the first temperature evaluation coefficient distribution cloud map is formed according to the first temperature evaluation coefficient of each coordinate. The specific formula for generating the first temperature evaluation coefficient is as follows:

[0017]

[0018] Among them, ET1 is the first temperature evaluation coefficient, RT is the bundle temperature, and ST is the flue gas temperature.

[0019] Furthermore, the specific logic for obtaining the second temperature evaluation coefficient distribution cloud map is as follows: the second temperature evaluation coefficient of each coordinate pixel point is calculated according to the smoke velocity distribution cloud map, the pressure distribution cloud map and the first temperature evaluation coefficient distribution cloud map, and the second temperature evaluation coefficient distribution cloud map is formed according to the second temperature evaluation coefficient of each coordinate. The specific formula for obtaining the second temperature evaluation coefficient is as follows:

[0020]

[0021] Among them, ET2 is the second temperature evaluation coefficient, ET1 is the first temperature evaluation coefficient, Pa is the pressure inside the filter rod, and V is the flue gas velocity.

[0022] Further, the specific logic for obtaining the third temperature evaluation coefficient distribution cloud map is as follows: according to the condensed nicotine distribution cloud map, the nicotine mole fraction distribution cloud map, the condensed glycerol distribution cloud map, the glycerol mole fraction distribution cloud map, the condensed water distribution cloud map and the water mole fraction distribution cloud map, the condensed substance evaluation coefficient of each coordinate pixel is calculated, and the condensed substance evaluation coefficient distribution cloud map is formed according to the condensed substance evaluation coefficient of each coordinate; then, the third temperature evaluation coefficient of each coordinate pixel is calculated according to the condensed substance evaluation coefficient distribution cloud map and the second temperature evaluation coefficient distribution cloud map, and the third temperature evaluation coefficient distribution cloud map is formed according to the third temperature evaluation coefficient of each coordinate; the specific logic for generating the condensed substance evaluation coefficient is as follows:

[0023] Ln=Yj*M Yj *C Yj +CHO*M CHO *C CHO+HO*M HO *C HO

[0024] Where Ln is the condensate evaluation coefficient, Yj is the nicotine condensation volume, M Yj is the molar fraction of nicotine, C Yj is the specific heat capacity of nicotine, CHO is the condensation volume of glycerol, M CHO is the mole fraction of glycerol, C CHO is the specific heat capacity of glycerol, HO is the condensed volume of water, M HO is the water mole fraction, C HO is the specific heat capacity of water;

[0025] The specific formula for generating the third temperature evaluation coefficient is:

[0026]

[0027] Among them, ET3 is the third temperature evaluation coefficient, and ET2 is the second temperature evaluation coefficient.

[0028] Furthermore, the specific logic for generating the temperature distribution cloud map is: the filter rod temperature at each coordinate pixel point is calculated according to the material parameters and the third temperature evaluation coefficient distribution cloud map, and the filter rod temperature distribution cloud map is formed according to the filter rod temperature at each coordinate. The specific formula for calculating the filter rod temperature is:

[0029]

[0030] Among them, TY is the filter rod temperature, ET3 is the third temperature evaluation coefficient, is the thermal conductivity of the filter rod material, and is the heat exchange coefficient between the filter rod material and the flue gas.

[0031] The present invention further provides a system for measuring the temperature of a filter rod during a suction process based on a simulation model, wherein the system is used to implement any step of the method for measuring the temperature of a filter rod during a suction process based on a simulation model, specifically comprising:

[0032] The geometric modeling module is used to obtain the design parameters of the filter rod and to perform collective modeling of the filter rod according to the design parameters;

[0033] A suction simulation module is used to divide the model into hexahedral grids using ICEM, measure the material parameters of the filter rod by experimental means, set boundary conditions and material parameters, perform CFD simulation analysis on the filter rod, obtain analysis results, and preprocess the analysis results, wherein the analysis results include a smoke velocity distribution cloud map of the filter rod, a filter rod pressure distribution cloud map, a smoke temperature distribution cloud map, a tow temperature distribution cloud map, a condensed nicotine distribution cloud map, a nicotine mole fraction distribution cloud map, a condensed glycerol distribution cloud map, a glycerol mole fraction distribution cloud map, a condensed water distribution cloud map, and a water mole fraction distribution cloud map;

[0034] A temperature analysis module, used for fusing the pre-processed flue gas temperature distribution cloud map and the tow temperature distribution cloud map to obtain a first temperature evaluation coefficient distribution cloud map;

[0035] A flue gas analysis module is used to fuse the filter rod pressure distribution cloud map, the flue gas velocity distribution cloud map and the first temperature evaluation coefficient distribution cloud map after pretreatment to obtain a second temperature evaluation coefficient distribution cloud map;

[0036] A condensation analysis module, for fusing the condensed nicotine distribution cloud map, the nicotine mole fraction distribution cloud map, the condensed glycerol distribution cloud map, the glycerol mole fraction distribution cloud map, the condensed water distribution cloud map, and the water mole fraction distribution cloud map to obtain a condensed substance evaluation coefficient cloud map, and fusing the second temperature evaluation coefficient distribution cloud map and the condensed substance evaluation coefficient distribution cloud map to obtain a third temperature evaluation coefficient distribution cloud map;

[0037] The temperature measurement module is used to generate a filter rod temperature distribution cloud map according to the material parameters and the third temperature evaluation coefficient distribution cloud map, and obtain the temperature of the filter rod at the position to be measured according to the filter rod temperature distribution cloud map.

[0038] The present invention further provides a computer-readable storage medium, characterized in that the storage medium internally stores a computer program that can be executed by a processor, and when the computer program is executed by the processor, it can implement the method for measuring the filter rod temperature during the suction process based on the simulation model.

[0039] Compared with the prior art, the present invention has the following beneficial effects:

[0040] The present invention simulates the condensation phenomenon of the main condensate in the suction process, generates a condensation material evaluation coefficient and a condensation material evaluation coefficient cloud map according to the simulation results, and can directly see the influence of the condensation phenomenon of each part of the filter rod on the temperature measurement, so that the temperature measurement accuracy is higher;

[0041] The present invention also takes into account the material properties of the filter rod and obtains a temperature distribution cloud map based on the analysis results of the material properties of the filter rod and the simulation results, which not only improves the measurement accuracy, but also provides a scientific basis for the material selection and structural optimization of the filter rod, thereby improving product performance and enhancing the user's smoking experience. BRIEF DESCRIPTION OF THE DRAWINGS

[0042] Figure 1 It is a schematic diagram of the overall method flow of the present invention.

[0043] Figure 2 This is the flue gas temperature distribution cloud diagram within 2 seconds.

[0044] Figure 3 This is the condensation distribution cloud diagram of propylene glycol within 2 seconds.

[0045] Figure 4 It is a schematic diagram of the overall system structure of the present invention. DETAILED DESCRIPTION

[0046] In order to make the objectives, technical solutions and advantages of the present invention more clearly understood, the present invention is further described in detail below in conjunction with specific embodiments.

[0047] It should be noted that, unless otherwise defined, the technical terms or scientific terms used in the present invention should be understood by people with ordinary skills in the field to which the present invention belongs. The words "first", "second" and similar words used in the present invention do not indicate any order, quantity or importance, but are only used to distinguish different components. "Include" or "comprise" and similar words mean that the elements or objects appearing before the word include the elements or objects listed after the word and their equivalents, without excluding other elements or objects. "Connect" or "connected" and similar words are not limited to physical or mechanical connections, but may include electrical connections, whether direct or indirect. "Up", "down", "left", "right" and the like are only used to indicate relative positional relationships. When the absolute position of the described object changes, the relative positional relationship may also change accordingly.

[0048] Example:

[0049] See also Figures 1 to 3 , the present invention provides a technical solution:

[0050] A method for measuring filter rod temperature during a suction process based on a simulation model, the specific steps comprising:

[0051] Step 1: Obtain the design parameters of the filter rod, and perform collective modeling of the filter rod according to the design parameters;

[0052] Step 2: using ICEM to divide the model into hexahedral grids, measuring the material parameters of the filter rod by experimental means, setting boundary conditions and material parameters, performing CFD simulation analysis on the filter rod, obtaining analysis results, and preprocessing the analysis results, wherein the analysis results include a smoke velocity distribution cloud map of the filter rod, a filter rod pressure distribution cloud map, a smoke temperature distribution cloud map, a tow temperature distribution cloud map, a condensed nicotine distribution cloud map, a nicotine mole fraction distribution cloud map, a condensed glycerol distribution cloud map, a glycerol mole fraction distribution cloud map, a condensed water distribution cloud map, and a water mole fraction distribution cloud map;

[0053] The specific logic of using ICEM to divide the model into hexahedral meshes is as follows: ICEM is used to perform two O-shaped cuts on the three-dimensional model, and the corresponding line association is established to divide the block into two blocks, obtaining a mesh model with a basic mesh number of 463644, and the mesh division quality is above 0.7.

[0054] The simulation mainly analyzes the flow field distribution of cigarette smoke after 2s of smoking in the cooling filter. The cycle is 0.2s and the total cycle is 2s. When the boundary conditions are added, the smoke inlet temperature is 260℃ and the temperature of the reflux generated at the smoke outlet is 42.85℃. The boundary conditions are set as velocity inlet and pressure outlet. The pressure transition of the smoke outlet pressure is loaded by UDF. The wall surface is formed paper with a width of 25mm. The filter rod and the external environment realize thermal convection. The wall boundary condition is used. The natural convection heat transfer coefficient is 10W / (m 2 ·K).

[0055] The structure adopts a multiphase flow mixture model. The phase material of the first phase is mixture-template, the second phase is waterliquid, the third phase is C10H14N2, and the fourth phase is C3H8O2. The material parameters are the density of acetate rod, specific heat capacity of acetate fiber, thermal conductivity of acetate fiber, flue gas density, specific heat capacity of flue gas, thermal conductivity of flue gas and heat transfer coefficient between flue gas and acetate fiber measured by experimental means; in addition, initial standardization is used during initialization, the initialization parameter assignment UDS is 27, and the remaining parameters are all 0.

[0056] During simulation, the flue gas temperature distribution cloud map and the fiber bundle temperature distribution cloud map can be directly obtained.

[0057] During simulation, the smoke temperature distribution cloud map, condensed nicotine distribution cloud map, nicotine mole fraction distribution cloud map, condensed glycerol distribution cloud map, glycerol mole fraction distribution cloud map, condensed water distribution cloud map and water mole fraction distribution cloud map can be obtained by performing condensation analysis on the inside of the filter rod.

[0058] The theory of the condensation process mainly involves two important concepts, namely saturated vapor pressure and saturated temperature. The so-called saturated vapor pressure refers to the gas pressure when the gas is in dynamic equilibrium with the corresponding saturated liquid at a specific temperature, while the saturated temperature refers to the lowest temperature at which the gas is converted into a liquid. When the temperature of the gas drops below the saturation temperature, its saturated vapor pressure will also decrease accordingly. When the saturated vapor pressure of the gas is less than the external pressure, the gas molecules will begin to gather to form droplets, and then condensation will occur. Condensation is an important part of the design of the cooling structure of the filter rod, which involves the process of converting gas or steam into liquid. This process is based on the phase change law of matter. When the gas molecules gather to a certain extent, they will form liquid particles, thereby undergoing a phase transition. In this way, condensation can absorb the temperature of the high-temperature flue gas very well, reduce the sensory temperature of the flue gas, and increase the smoking taste.

[0059] Specifically include: calculating the saturation temperature of the gas, calculating the mass of the condensed liquid through the saturation temperature, and the specific logic for calculating the saturation temperature is:

[0060] T SAT =73.213+7980*ρ gas *Y i,gas

[0061] Among them, T SAT is the saturation temperature, ρ gas is the smoke density, Y i,gas is the molar volume of condensable substances in smoke. Condensable substances in smoke include water, nicotine and glycerol;

[0062] The specific formula used to calculate the mass of condensed liquid is:

[0063]

[0064] Among them, m l is the mass of condensed liquid, VOF gas is the volume fraction of smoke, T gas is the flue gas temperature; this process is compiled by CFD. During the simulation, the flue gas temperature distribution cloud map, condensed nicotine distribution cloud map, nicotine mole fraction distribution cloud map, condensed glycerol distribution cloud map, glycerol mole fraction distribution cloud map, condensed water distribution cloud map and water mole fraction distribution cloud map can be obtained by performing condensation analysis on the inside of the filter rod;

[0065] The flue gas velocity distribution cloud map and filter rod pressure distribution cloud map can be obtained by analyzing the unbalanced heat transfer process;

[0066] Non-equilibrium heat transfer describes the energy exchange process of a system when the heat is not balanced. In non-equilibrium heat conduction, the transfer of heat can be achieved mainly through three core pathways: transfer, convection, and radiation. Unlike the mechanism of heat transfer in equilibrium, the speed of heat transfer in non-equilibrium is not simply set according to the principle of heat conduction, but is actually affected by a combination of factors such as temperature gradient, fluid velocity, and changes in physical properties. In the non-equilibrium state of the heat transfer process, the speed of heat transfer is not only affected by the temperature difference, but also by various other thermodynamic parameters. For example, when flowing liquids flow in pipes due to factors such as friction and turbulence, they generate heat energy. At this moment, heat propagation is affected not only by the temperature difference, but also by multiple variables such as the speed of the fluid movement, the appearance of the fluid pipe, and the material used. In the process of non-equilibrium heat transfer, phase changes and various chemical changes of matter may occur, and the changes will further change the speed and mechanism of heat transport.

[0067] The unbalanced heat transfer process is obtained by defining the velocity and pressure boundaries. Specifically, the flow time is obtained. When the flow time is greater than 2s, the velocity is set to 0. When the flow time is less than 2s, the velocity is:

[0068] V=0.5936×sin(1.57×t)

[0069] Among them, V is the flue gas velocity and t is the flow time;

[0070] When the flow time is greater than 2s, the pressure is set to 0, and when the flow time is less than 2s, the pressure is:

[0071] Pa=14.42-864.69×t+432.35×t 2

[0072] Wherein, Pa is the filter rod pressure.

[0073] The preprocessing operation is: to mark the coordinates of the pixel points of various distribution cloud maps in the analysis results, the bottom pixel row is the first row, the leftmost pixel column is the first column, the pixel point positions of the first recognition image and the second recognition image are mapped one by one through the coordinates, so that each pixel point has a unique coordinate value, and each pixel point position of the various distribution cloud maps in the analysis results is mapped one by one.

[0074] Step 3: Fusing the pre-processed flue gas temperature distribution cloud map and the tow temperature distribution cloud map to obtain a first temperature evaluation coefficient distribution cloud map;

[0075] See also Figure 2 The smoke temperature distribution cloud diagram within 2s is as follows: Figure 2 As shown;

[0076] The specific logic for obtaining the first temperature evaluation coefficient distribution cloud map is as follows: the first temperature evaluation coefficient of each coordinate pixel point is calculated according to the flue gas temperature distribution cloud map and the tow temperature distribution cloud map, and the first temperature evaluation coefficient distribution cloud map is formed according to the first temperature evaluation coefficient of each coordinate. The specific formula for generating the first temperature evaluation coefficient is as follows:

[0077]

[0078] Among them, ET1 is the first temperature evaluation coefficient, RT is the bundle temperature, and ST is the flue gas temperature. The first temperature evaluation coefficient ET1 reflects the influence of the internal heat source on the filter rod temperature. The larger the value, the greater the influence. The generation of this coefficient can provide an important basis for the determination of the filter rod temperature.

[0079] Step 4: Fusing the pre-processed filter rod pressure distribution cloud map, the smoke velocity distribution cloud map and the first temperature evaluation coefficient distribution cloud map to obtain a second temperature evaluation coefficient distribution cloud map;

[0080] The specific logic for obtaining the second temperature evaluation coefficient distribution cloud map is as follows: the second temperature evaluation coefficient of each coordinate pixel point is calculated according to the smoke velocity distribution cloud map, the pressure distribution cloud map and the first temperature evaluation coefficient distribution cloud map, and the second temperature evaluation coefficient distribution cloud map is formed according to the second temperature evaluation coefficient of each coordinate. The specific formula for obtaining the second temperature evaluation coefficient is as follows:

[0081]

[0082] Among them, ET2 is the second temperature evaluation coefficient, ET1 is the first temperature evaluation coefficient, Pa is the pressure inside the filter rod, and V is the smoke velocity. The second temperature evaluation coefficient ET2 reflects the influence of the internal heat source and smoke flow on the filter rod temperature. The larger its value, the greater the influence. The generation of this coefficient can provide an important basis for the determination of the filter rod temperature. The greater the pressure Pa inside the filter rod, the greater the effect of the gas flow inside the filter rod on the increase of the filter rod temperature. The smoke velocity V reflects the flow velocity of the smoke inside the filter rod. The larger its value, the greater the effect of the gas flow inside the filter rod on the reduction of the filter rod temperature.

[0083] Step 5: the condensed nicotine distribution cloud map, the nicotine mole fraction distribution cloud map, the condensed glycerol distribution cloud map, the glycerol mole fraction distribution cloud map, the condensed water distribution cloud map, and the water mole fraction distribution cloud map are merged to obtain a condensed substance evaluation coefficient cloud map, and the second temperature evaluation coefficient distribution cloud map and the condensed substance evaluation coefficient distribution cloud map are merged to obtain a third temperature evaluation coefficient distribution cloud map;

[0084] See also Figure 3 , glycerol condenses within 2 seconds. The glycerol distribution cloud diagram is as follows Figure 3 shown.

[0085] The specific logic for obtaining the third temperature evaluation coefficient distribution cloud map is as follows: the condensed substance evaluation coefficient of each coordinate pixel point is calculated according to the condensed nicotine distribution cloud map, the nicotine mole fraction distribution cloud map, the condensed glycerol distribution cloud map, the glycerol mole fraction distribution cloud map, the condensed water distribution cloud map and the water mole fraction distribution cloud map, and the condensed substance evaluation coefficient distribution cloud map is formed according to the condensed substance evaluation coefficient of each coordinate; the third temperature evaluation coefficient of each coordinate pixel point is calculated according to the condensed substance evaluation coefficient distribution cloud map and the second temperature evaluation coefficient distribution cloud map, and the third temperature evaluation coefficient distribution cloud map is formed according to the third temperature evaluation coefficient of each coordinate; the specific logic for generating the condensed substance evaluation coefficient is as follows:

[0086] Ln=Yj*M Yj *C Yj +CHO*M CHI *C CHO +HO*M HO *C HO

[0087] Where Ln is the condensate evaluation coefficient, Yj is the nicotine condensation volume, M Yj is the molar fraction of nicotine, C Yj is the specific heat capacity of nicotine, CHO is the condensation volume of glycerol, M CHO is the mole fraction of glycerol, C CHO is the specific heat capacity of glycerol, HO is the condensed volume of water, M HO is the water mole fraction, C HO is the specific heat capacity of water; the condensed material evaluation coefficient Ln reflects the thermodynamic properties of the condensate in the filter rod. The larger its value, the higher the degree of the filter rod temperature drop caused by the condensed material in the filter rod; Yj, CHO, and HO reflect the condensed volume of each substance. The larger the volume, the more the component is in the mixture, and the easier it is to cause the filter rod temperature to drop; M Yj 、M CHO and M HO Reflects the proportion of each component. The larger the molar fraction, the higher the proportion of the component in the mixture; C Yj , C CHO and C HO Reflects the ability of each component to absorb or release heat. The larger the specific heat capacity, the more heat the component can absorb or release when the temperature changes, and the lower the filter rod temperature.

[0088] The specific formula for generating the third temperature evaluation coefficient is:

[0089]

[0090] Among them, ET3 is the third temperature evaluation coefficient, and ET2 is the second temperature evaluation coefficient. The third temperature evaluation coefficient ET3 reflects the influence of internal condensed matter, heat source and flue gas flow on the filter rod temperature. The larger the value, the greater the influence. The generation of this coefficient can provide an important basis for the determination of the filter rod temperature.

[0091] Step 6: Generate a filter rod temperature distribution cloud map according to the material parameters and the third temperature evaluation coefficient distribution cloud map, and obtain the temperature of the position to be measured on the filter rod according to the filter rod temperature distribution cloud map.

[0092] The specific logic for generating the temperature distribution cloud map is as follows: the filter rod temperature at each coordinate pixel is calculated based on the material parameters and the third temperature evaluation coefficient distribution cloud map, and the filter rod temperature distribution cloud map is constructed based on the filter rod temperature at each coordinate. The specific formula for calculating the filter rod temperature is as follows:

[0093]

[0094] Among them, TY is the filter rod temperature, ET3 is the third temperature evaluation coefficient, DR is the thermal conductivity of the filter rod material, EC is the heat transfer coefficient between the filter rod material and the flue gas, CA The filter rod temperature TY is the filter rod temperature that needs to be measured. The thermal conductivity coefficient DR of the filter rod material reflects the thermal conductivity of the filter rod. The larger the value, the stronger the thermal conductivity of the filter rod material, and the more effective the heat transfer through the material. The heat transfer coefficient EC between the filter rod material and the flue gas reflects the heat exchange capacity between the filter rod material and the flue gas. The larger the value, the faster the heat exchange between the filter rod material and the flue gas, and the easier it is for the filter rod to reach the equilibrium temperature. The specific heat capacity C of the filter rod material A It reflects the ability to absorb heat when the temperature rises. The larger the value, the more heat the filter rod material needs to heat up, and it is relatively difficult to reach a high temperature state.

[0095] See also Figure 4 The present invention further provides a system for measuring the temperature of a filter rod during a suction process based on a simulation model, wherein the system is used to implement any step of the method for measuring the temperature of a filter rod during a suction process based on a simulation model, specifically comprising:

[0096] The geometric modeling module is used to obtain the design parameters of the filter rod and to perform collective modeling of the filter rod according to the design parameters;

[0097] A suction simulation module is used to divide the model into hexahedral grids using ICEM, measure the material parameters of the filter rod by experimental means, set boundary conditions and material parameters, perform CFD simulation analysis on the filter rod, obtain analysis results, and preprocess the analysis results, wherein the analysis results include a smoke velocity distribution cloud map of the filter rod, a filter rod pressure distribution cloud map, a smoke temperature distribution cloud map, a tow temperature distribution cloud map, a condensed nicotine distribution cloud map, a nicotine mole fraction distribution cloud map, a condensed glycerol distribution cloud map, a glycerol mole fraction distribution cloud map, a condensed water distribution cloud map, and a water mole fraction distribution cloud map;

[0098] A temperature analysis module, used for fusing the pre-processed flue gas temperature distribution cloud map and the tow temperature distribution cloud map to obtain a first temperature evaluation coefficient distribution cloud map;

[0099] A flue gas analysis module is used to fuse the filter rod pressure distribution cloud map, the flue gas velocity distribution cloud map and the first temperature evaluation coefficient distribution cloud map after pretreatment to obtain a second temperature evaluation coefficient distribution cloud map;

[0100] A condensation analysis module, for fusing the condensed nicotine distribution cloud map, the nicotine mole fraction distribution cloud map, the condensed glycerol distribution cloud map, the glycerol mole fraction distribution cloud map, the condensed water distribution cloud map, and the water mole fraction distribution cloud map to obtain a condensed substance evaluation coefficient cloud map, and fusing the second temperature evaluation coefficient distribution cloud map and the condensed substance evaluation coefficient distribution cloud map to obtain a third temperature evaluation coefficient distribution cloud map;

[0101] The temperature measurement module is used to generate a filter rod temperature distribution cloud map according to the material parameters and the third temperature evaluation coefficient distribution cloud map, and obtain the temperature of the filter rod at the position to be measured according to the filter rod temperature distribution cloud map.

[0102] The present invention further provides a computer-readable storage medium, characterized in that the storage medium internally stores a computer program that can be executed by a processor, and when the computer program is executed by the processor, it can implement the method for measuring the filter rod temperature during the suction process based on the simulation model.

[0103] The above formulas are all dimensionless and numerical calculations. The formula is a formula for the most recent real situation obtained by collecting a large amount of data and performing software simulation. The preset parameters in the formula are set by technicians in this field according to actual conditions.

[0104] The above embodiments may be implemented in whole or in part by software, hardware, firmware or any other combination thereof. When implemented by software, the above embodiments may be implemented in whole or in part in the form of a computer program product. Those skilled in the art may appreciate that the units and algorithm steps of each example described in conjunction with the embodiments disclosed herein may be implemented by electronic hardware, or a combination of computer software and electronic hardware. Whether these functions are performed by hardware or software methods depends on the specific application and design constraints of the technical solution.

[0105] The units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, and may be located in one place or distributed on multiple network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.

[0106] The above description is only a specific implementation manner of the present application, but the protection scope of the present application is not limited thereto. Any technician familiar with the technical field can easily think of changes or substitutions within the technical scope disclosed in the present application, which should be included in the protection scope of the present application.

Claims

1. A method for measuring filter rod temperature during suction based on a simulation model, characterized in that: The specific steps include: Step 1: Obtain the design parameters of the filter rod, and perform collective modeling of the filter rod according to the design parameters; Step 2: using ICEM to divide the model into hexahedral grids, measuring the material parameters of the filter rod by experimental means, setting boundary conditions and material parameters, performing CFD simulation analysis on the filter rod, obtaining analysis results, and preprocessing the analysis results, wherein the analysis results include a smoke velocity distribution cloud map of the filter rod, a filter rod pressure distribution cloud map, a smoke temperature distribution cloud map, a tow temperature distribution cloud map, a condensed nicotine distribution cloud map, a nicotine mole fraction distribution cloud map, a condensed glycerol distribution cloud map, a glycerol mole fraction distribution cloud map, a condensed water distribution cloud map, and a water mole fraction distribution cloud map; Step 3: Fusing the pre-processed flue gas temperature distribution cloud map and the tow temperature distribution cloud map to obtain a first temperature evaluation coefficient distribution cloud map; Step 4: Fusing the pre-processed filter rod pressure distribution cloud map, the smoke velocity distribution cloud map and the first temperature evaluation coefficient distribution cloud map to obtain a second temperature evaluation coefficient distribution cloud map; Step 5: the condensed nicotine distribution cloud map, the nicotine mole fraction distribution cloud map, the condensed glycerol distribution cloud map, the glycerol mole fraction distribution cloud map, the condensed water distribution cloud map, and the water mole fraction distribution cloud map are merged to obtain a condensed substance evaluation coefficient cloud map, and the second temperature evaluation coefficient distribution cloud map and the condensed substance evaluation coefficient distribution cloud map are merged to obtain a third temperature evaluation coefficient distribution cloud map; Step 6: Generate a filter rod temperature distribution cloud map according to the material parameters and the third temperature evaluation coefficient distribution cloud map, and obtain the temperature of the position to be measured on the filter rod according to the filter rod temperature distribution cloud map.

2. The method for measuring the temperature of a filter rod during the suction process based on a simulation model according to claim 1, characterized in that: The preprocessing operation is: to mark the coordinates of the pixel points of various distribution cloud maps in the analysis results, the bottom pixel row is the first row, the leftmost pixel column is the first column, the pixel point positions of the first recognition image and the second recognition image are mapped one by one through the coordinates, so that each pixel point has a unique coordinate value, and each pixel point position of the various distribution cloud maps in the analysis results is mapped one by one.

3. The method for measuring the temperature of a filter rod during the suction process based on a simulation model according to claim 1, characterized in that: The specific logic for obtaining the first temperature evaluation coefficient distribution cloud map is as follows: the first temperature evaluation coefficient of each coordinate pixel point is calculated according to the flue gas temperature distribution cloud map and the tow temperature distribution cloud map, and the first temperature evaluation coefficient distribution cloud map is formed according to the first temperature evaluation coefficient of each coordinate. The specific formula for generating the first temperature evaluation coefficient is as follows: Among them, ET1 is the first temperature evaluation coefficient, RT is the bundle temperature, and ST is the flue gas temperature.

4. The method for measuring the temperature of a filter rod during the suction process based on a simulation model according to claim 1, characterized in that: The specific logic for obtaining the second temperature evaluation coefficient distribution cloud map is as follows: the second temperature evaluation coefficient of each coordinate pixel point is calculated according to the smoke velocity distribution cloud map, the pressure distribution cloud map and the first temperature evaluation coefficient distribution cloud map, and the second temperature evaluation coefficient distribution cloud map is formed according to the second temperature evaluation coefficient of each coordinate. The specific formula for obtaining the second temperature evaluation coefficient is as follows: Among them, ET2 is the second temperature evaluation coefficient, ET1 is the first temperature evaluation coefficient, Pa is the pressure inside the filter rod, and Y is the flue gas velocity.

5. The method for measuring the temperature of a filter rod during the suction process based on a simulation model according to claim 1, characterized in that: The specific logic for obtaining the third temperature evaluation coefficient distribution cloud map is as follows: according to the condensed nicotine distribution cloud map, the nicotine mole fraction distribution cloud map, the condensed glycerol distribution cloud map, the glycerol mole fraction distribution cloud map, the condensed water distribution cloud map and the water mole fraction distribution cloud map, the condensed substance evaluation coefficient of each coordinate pixel point is calculated, and the condensed substance evaluation coefficient distribution cloud map is formed according to the condensed substance evaluation coefficient of each coordinate; then, the third temperature evaluation coefficient of each coordinate pixel point is calculated according to the condensed substance evaluation coefficient distribution cloud map and the second temperature evaluation coefficient distribution cloud map, and the third temperature evaluation coefficient distribution cloud map is formed according to the third temperature evaluation coefficient of each coordinate; the specific logic for generating the condensed substance evaluation coefficient is as follows: Ln=Yj*M Yj *C Yj +CHO*M CHO *C CHO +HO*M HO *C HO Where Ln is the condensate evaluation coefficient, Yj is the nicotine condensation volume, M Yj is the molar fraction of nicotine, C Yj is the specific heat capacity of nicotine, CHO is the condensation volume of glycerol, M CHO is the mole fraction of glycerol, C CHO is the specific heat capacity of glycerol, HO is the condensed volume of water, M HO is the water mole fraction, C HO is the specific heat capacity of water; The specific formula for generating the third temperature evaluation coefficient is: Among them, ET3 is the third temperature evaluation coefficient, and ET2 is the second temperature evaluation coefficient.

6. The method for measuring filter rod temperature during suction based on a simulation model according to claim 1, characterized in that: The specific logic for generating the temperature distribution cloud map is as follows: the filter rod temperature at each coordinate pixel is calculated based on the material parameters and the third temperature evaluation coefficient distribution cloud map, and the filter rod temperature distribution cloud map is constructed based on the filter rod temperature at each coordinate. The specific formula for calculating the filter rod temperature is as follows: Among them, TY is the filter rod temperature, ET3 is the third temperature evaluation coefficient, DR is the thermal conductivity of the filter rod material, EC is the heat transfer coefficient between the filter rod material and the flue gas, C A is the specific heat capacity of the filter rod material.

7. A filter rod temperature measurement system during the suction process based on a simulation model, characterized in that: The system is used to implement any step of the method for measuring the temperature of the filter rod during the suction process based on the simulation model according to claims 1-6, specifically comprising: The geometric modeling module is used to obtain the design parameters of the filter rod and to perform collective modeling of the filter rod according to the design parameters; A suction simulation module is used to divide the model into hexahedral grids using ICEM, measure the material parameters of the filter rod by experimental means, set boundary conditions and material parameters, perform CFD simulation analysis on the filter rod, obtain analysis results, and preprocess the analysis results, wherein the analysis results include a smoke velocity distribution cloud map of the filter rod, a filter rod pressure distribution cloud map, a smoke temperature distribution cloud map, a tow temperature distribution cloud map, a condensed nicotine distribution cloud map, a nicotine mole fraction distribution cloud map, a condensed glycerol distribution cloud map, a glycerol mole fraction distribution cloud map, a condensed water distribution cloud map, and a water mole fraction distribution cloud map; A temperature analysis module, used for fusing the pre-processed flue gas temperature distribution cloud map and the tow temperature distribution cloud map to obtain a first temperature evaluation coefficient distribution cloud map; A flue gas analysis module is used to fuse the filter rod pressure distribution cloud map, the flue gas velocity distribution cloud map and the first temperature evaluation coefficient distribution cloud map after pretreatment to obtain a second temperature evaluation coefficient distribution cloud map; A condensation analysis module, for fusing the condensed nicotine distribution cloud map, the nicotine mole fraction distribution cloud map, the condensed glycerol distribution cloud map, the glycerol mole fraction distribution cloud map, the condensed water distribution cloud map, and the water mole fraction distribution cloud map to obtain a condensed substance evaluation coefficient cloud map, and fusing the second temperature evaluation coefficient distribution cloud map and the condensed substance evaluation coefficient distribution cloud map to obtain a third temperature evaluation coefficient distribution cloud map; The temperature measurement module is used to generate a filter rod temperature distribution cloud map according to the material parameters and the third temperature evaluation coefficient distribution cloud map, and obtain the temperature of the filter rod at the position to be measured according to the filter rod temperature distribution cloud map.

8. A computer-readable storage medium, characterized in that: The storage medium internally stores a computer program that can be executed by a processor, and when the computer program is executed by the processor, the method for measuring the temperature of a filter rod during a suction process based on a simulation model as described in any one of claims 1 to 6 can be implemented.

Images