Self-correction optimization method and system for heat transfer boundary conditions of water jacket heater

Abstract

The invention discloses a self-correction optimization method and system for heat transfer boundary conditions of a water jacket heater, and particularly relates to the technical field of condition optimization control. The method comprises the following steps: constructing a thermal fatigue deformation trend estimation model by collecting a furnace wall surface temperature sequence, furnace wall surface heat flux distribution information and a furnace wall thermal stress response curve of a furnace body structure component, and generating a geometric parameter drift vector; performing structure equivalent mapping analysis on the geometric parameter drift vector, and extracting contact thermal resistance change information; performing furnace wall heat flow field reconstruction on the geometric parameter drift vector to obtain a dynamic correction vector; and the contact thermal resistance change information and the heat flow correction vector are comprehensively analyzed, a time-varying boundary condition inversion model is constructed, and a multi-parameter coupling boundary correction set is output and used for updating the boundary conditions in the water jacket heating furnace heat transfer simulation model, so that identification and self-adaptive correction of the boundary conditions are realized, and the heat transfer control precision and the heat efficiency are improved.

Description

Technical Field

[0001] The present invention relates to the technical field of condition optimization control, and more specifically, the present invention relates to a method and system for self-correcting and optimizing the heat transfer boundary conditions of a water jacket heating furnace. Background Art

[0002] During the long-term operation of a water jacket heating furnace, due to the continuous thermal cycling load on structural components such as furnace walls or pipes, irreversible thermal expansion accumulation and geometric deformation occur, which in turn causes dynamic drift of heat transfer boundary conditions such as contact thermal resistance and heat flux distribution.

[0003] Existing self-correcting and optimizing methods generally calculate based on the initial design geometric parameters, lacking real-time perception and dynamic response to structural deformation and boundary condition changes during operation. Summary of the Invention

[0004] In order to overcome the above-mentioned defects of the prior art, embodiments of the present invention provide a method and system for self-correcting and optimizing the heat transfer boundary conditions of a water jacket heating furnace to solve the problems raised in the above background art.

[0005] To achieve the above object, the present invention provides the following technical solutions: A method for self-correcting and optimizing the heat transfer boundary conditions of a water jacket heating furnace, comprising the following steps: S1. Collect the furnace wall surface temperature sequence, furnace wall surface heat flux distribution information, and furnace wall thermal stress response curve of the furnace body structural components under thermal cycling conditions, and obtain the comprehensive state parameters of the furnace body structural components under thermal cycling conditions; S2. Based on the comprehensive state parameters of the furnace body structural components under thermal cycling conditions, construct a thermal fatigue deformation trend estimation model, and generate a geometric parameter drift vector for evaluating the deformation of the furnace body structural components; S3. Perform structural equivalent mapping analysis on the geometric parameter drift vector, and extract the contact thermal resistance change information in the heat transfer path; S4. Based on the geometric parameter drift vector, reconstruct the furnace wall heat flow field, and obtain the dynamic correction vector of the actual heat flux distribution; S5. Comprehensively analyze the contact thermal resistance change information in the heat transfer path and the dynamic correction vector of the actual heat flux distribution, construct a time-varying boundary condition inversion model, and output a multi-parameter coupled boundary correction set; S6. Update the boundary conditions in the heat transfer simulation model of the water jacket heating furnace according to the multi-parameter coupled boundary correction set.

[0006] In a preferred embodiment, S1 is specifically: Collect the furnace wall surface temperature sequence, furnace wall surface heat flux distribution information, and furnace wall thermal stress response curve of the furnace body structural components under long-term thermal cycling conditions during the operation of the water jacket heating furnace; Perform a time-series variation analysis on the surface temperature sequence of the furnace wall to obtain the temperature change characteristics of the furnace wall surface; Perform a spatial distribution analysis on the heat flux distribution information of the furnace wall surface to obtain the spatial distribution characteristics of the heat flux on the furnace wall surface; Perform feature recognition on the thermal stress response curve of the furnace wall to obtain the response characteristics of the thermal stress change of the furnace wall; Combine the temperature change characteristics of the furnace wall surface, the spatial distribution characteristics of the heat flux on the furnace wall surface, and the response characteristics of the thermal stress change of the furnace wall to obtain the comprehensive state parameters of the furnace body structure components under thermal cycling conditions.

[0007] In a preferred embodiment, S2 is specifically as follows: Based on the comprehensive state parameters of the furnace body structure components under thermal cycling conditions, establish the correlation between thermal expansion displacement and structural deformation; According to the correlation between thermal expansion displacement and structural deformation, determine the thermal fatigue deformation trend of the furnace body structure components; According to the thermal fatigue deformation trend of the furnace body structure components, calculate the geometric shape drift amount generated by the furnace body structure components under thermal cycling conditions; Based on the geometric shape drift amount generated by the furnace body structure components under thermal cycling conditions, generate a geometric parameter drift vector for evaluating the deformation of the furnace body structure components.

[0008] In a preferred embodiment, S3 is specifically as follows: Adopt the structural equivalent mapping method to convert the geometric shape drift vector generated by the furnace body structure components under thermal cycling conditions into the structural change characteristics of the heat transfer path; According to the structural change characteristics of the heat transfer path, analyze the change law of the heat conduction interface contact area of the furnace body structure components, and determine the change amount of the heat conduction interface contact area; According to the change amount of the heat conduction interface contact area, analyze the change of the heat transfer capacity of the thermal contact interface of the furnace body structure components, and obtain the change information of the contact thermal resistance of the thermal contact interface of the furnace body structure components.

[0009] In a preferred embodiment, S4 is specifically as follows: Based on the geometric shape drift vector generated by the furnace body structure components under thermal cycling conditions, establish a numerical reconstruction model of the furnace wall heat flow field; According to the numerical reconstruction model of the furnace wall heat flow field, perform a reconstruction analysis on the spatial distribution characteristics of the heat flux on the furnace wall surface, and obtain the change law of the heat flux distribution on the furnace wall surface under the actual operating conditions of the furnace wall structure components; According to the change law of the heat flux distribution on the furnace wall surface, determine the position and heat flux change amplitude of the local offset area of the heat flux of the furnace wall structure components, and generate a dynamic correction vector of the actual heat flow distribution of the furnace wall structure components.

[0010] In a preferred embodiment, S5 specifically includes: Based on the contact thermal resistance change information of the thermal contact interface of the furnace body structure components, establish a time-varying model of the change of the contact thermal resistance parameters; Based on the dynamic correction vector of the actual heat flux distribution of the furnace wall structure components, establish a time-varying model of the dynamic correction parameters of the heat flow; According to the time-varying model of the change of the contact thermal resistance parameters and the time-varying model of the dynamic correction parameters of the heat flow, conduct a multi-parameter comprehensive analysis of the heat transfer boundary conditions of the furnace body structure components, and construct an inverse model of the boundary conditions for the dynamic change of the heat transfer characteristics of the furnace body structure components; Through the inverse model of the boundary conditions, output a multi-parameter coupling boundary correction set of the heat transfer boundary conditions under the thermal cycle conditions of the furnace body structure components.

[0011] In a preferred embodiment, S6 specifically includes: Based on the multi-parameter coupling boundary correction set, determine the contact thermal resistance correction parameters and the furnace wall heat flux distribution correction parameters; Update the contact thermal resistance correction parameters to the contact thermal resistance boundary conditions of the thermal contact interface in the heat transfer simulation model of the water jacket heating furnace; Update the furnace wall heat flux distribution correction parameters to the spatial distribution boundary conditions of the furnace wall surface heat flux in the heat transfer simulation model of the water jacket heating furnace; The updated contact thermal resistance boundary conditions of the thermal contact interface and the spatial distribution boundary conditions of the furnace wall surface heat flux form the current heat transfer boundary conditions of the heat transfer simulation model of the water jacket heating furnace.

[0012] On the other hand, the present invention provides a self-correcting and optimizing system for the heat transfer boundary conditions of a water jacket heating furnace, including: A state acquisition module, which acquires the furnace wall surface temperature sequence, the furnace wall surface heat flux distribution information, and the furnace wall thermal stress response curve of the furnace body structure components under thermal cycle conditions, and obtains the comprehensive state parameters of the furnace body structure components under thermal cycle conditions; A trend estimation module, which constructs a thermal fatigue deformation trend estimation model based on the comprehensive state parameters of the furnace body structure components under thermal cycle conditions, and generates a geometric parameter drift vector used to evaluate the deformation of the furnace body structure components; A thermal resistance analysis module, which conducts a structural equivalent mapping analysis on the geometric parameter drift vector and extracts the contact thermal resistance change information in the heat transfer path; A heat flux reconstruction module, which reconstructs the furnace wall heat flow field based on the geometric parameter drift vector and obtains the dynamic correction vector of the actual heat flux distribution; A boundary inversion module, which comprehensively analyzes the contact thermal resistance change information in the heat transfer path and the dynamic correction vector of the actual heat flux distribution, constructs a time-varying boundary condition inversion model, and outputs a multi-parameter coupling boundary correction set; The model update module updates the boundary conditions in the heat transfer simulation model of the water jacket heating furnace according to the multi-parameter coupling boundary correction set.

[0013] The technical effects and advantages of a method and system for self-correcting and optimizing the heat transfer boundary conditions of a water jacket heating furnace according to the present invention are as follows: By collecting the furnace wall surface temperature sequence, heat flux distribution information, and thermal stress response curve, the operating state of the furnace body structure under thermal cycling conditions can be comprehensively perceived, improving the accuracy of condition monitoring; by constructing a thermal fatigue deformation trend estimation model, the cumulative drift characteristics of the geometric structure can be identified in a timely manner, effectively revealing the correlation between structural deformation and boundary changes; through structural equivalent mapping analysis and furnace wall heat flow field reconstruction, the change characteristics of contact thermal resistance and heat flow distribution can be accurately extracted, enhancing the quantification ability of boundary condition changes; by constructing a time-varying boundary condition inversion model, the accurate deduction of heat transfer boundary conditions can be achieved; by using a multi-parameter coupling boundary correction set to update the boundary conditions in the heat transfer simulation model of the water jacket heating furnace, problems such as local hot spots and reduced efficiency caused by model mismatch can be effectively avoided, improving the simulation accuracy, optimization stability, and operation safety of the heating process. Brief Description of the Drawings

[0014] Figure 1 It is a schematic diagram of a method for self-correcting and optimizing the heat transfer boundary conditions of a water jacket heating furnace according to the present invention; Figure 2 It is a schematic structural diagram of a system for self-correcting and optimizing the heat transfer boundary conditions of a water jacket heating furnace according to the present invention. Detailed Embodiments

[0015] Next, the technical solutions in the embodiments of the present invention will be clearly and completely described in conjunction with the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are only a part of the embodiments of the present invention, rather than all of the embodiments. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of the present invention without creative efforts shall fall within the protection scope of the present invention.

[0016] Embodiment 1

[0017] Figure 1 A method for self-correcting and optimizing the heat transfer boundary conditions of a water jacket heating furnace according to the present invention is given, which includes the following steps: S1. Collect the furnace wall surface temperature sequence, furnace wall surface heat flux distribution information, and furnace wall thermal stress response curve of the furnace body structure components under thermal cycling conditions to obtain the comprehensive state parameters of the furnace body structure components under thermal cycling conditions; S2. Based on the comprehensive state parameters of the furnace body structure components under thermal cycling conditions, construct a thermal fatigue deformation trend estimation model and generate a geometric parameter drift vector for evaluating the deformation of the furnace body structure components; S3. Conduct a structural equivalent mapping analysis on the geometric parameter drift vector to extract the information on the change of contact thermal resistance in the heat transfer path; S4. Based on the geometric parameter drift vector, reconstruct the heat flux field of the furnace wall to obtain the dynamic correction vector of the actual heat flux distribution; S5. Conduct a comprehensive analysis on the information on the change of contact thermal resistance in the heat transfer path and the dynamic correction vector of the actual heat flux distribution, construct a time-varying boundary condition inversion model, and output a multi-parameter coupled boundary correction set; S6. Update the boundary conditions in the heat transfer simulation model of the water jacket furnace according to the multi-parameter coupled boundary correction set.

[0018] S1. Collect the furnace wall surface temperature sequence, the furnace wall surface heat flux distribution information, and the furnace wall thermal stress response curve of the furnace body structure components under thermal cycling conditions, and obtain the comprehensive state parameters of the furnace body structure components under thermal cycling conditions, including: Collect the furnace wall surface temperature sequence, the furnace wall surface heat flux distribution information, and the furnace wall thermal stress response curve of the furnace body structure components of the water jacket furnace under long-term thermal cycling conditions during the operation of the water jacket furnace; Use the temperature measurement sensors installed on the furnace wall surface of the furnace body structure components of the water jacket furnace to continuously collect in real time the temperature change of the furnace wall surface of the furnace body structure components under long-term thermal cycling conditions, and arrange the collected continuous temperature data in sequence at fixed time intervals to form the furnace wall surface temperature sequence; adopt the heat flux measurement devices distributed at different spatial positions on the furnace wall surface to collect in real time the magnitude of the heat flux on the furnace wall surface, and correspond it with the position coordinate information to form the furnace wall surface heat flux distribution information; at the same time, adopt the thermal stress detection device to monitor and record in real time the change of the thermal stress suffered by the furnace body structure components during the thermal cycling operation, and plot the recorded data as the response curve of the thermal stress changing with time.

[0019] Conduct a time series change analysis on the furnace wall surface temperature sequence to obtain the temperature change characteristics of the furnace wall surface; In order to obtain the temperature change characteristics in the furnace wall surface temperature sequence, adopt the time series analysis algorithm, use the furnace wall surface temperature sequence as the input, calculate the difference between the temperature value at each moment in the furnace wall surface temperature sequence and the temperature value at the adjacent time before it, and calculate the change trend of the difference with time, and then determine the change law of the furnace wall surface temperature by statistically analyzing the frequency, amplitude, and duration of the appearance of the difference change trend, and finally output the time series change characteristics of the furnace wall surface temperature.

[0020] Conduct a spatial distribution analysis on the furnace wall surface heat flux distribution information to obtain the spatial distribution characteristics of the heat flux on the furnace wall surface; To obtain the spatial distribution characteristics of the heat flux on the furnace wall surface, a spatial distribution analysis method is adopted. The heat flux distribution information on the furnace wall surface is spatially divided according to the geometric coordinates of the furnace wall. By statistically analyzing the heat flux data in each divided area, the average value, peak value of the heat flux in each area, and the degree of non-uniformity of the heat flux distribution (the difference between the average value and the peak value of the heat flux) are calculated, and the spatial distribution law of the heat flux on the surface of the furnace wall structural components is determined. Finally, the spatial distribution characteristics of the heat flux on the furnace wall surface are output.

[0021] Perform feature recognition on the furnace wall thermal stress response curve to obtain the response characteristics of the furnace wall thermal stress change; To obtain the response characteristics of the furnace wall thermal stress change, a feature recognition method is adopted to process the furnace wall thermal stress response curve. By identifying the peak value, valley value, and thermal stress change slope in the thermal stress response curve, the typical states during the thermal stress change process are determined. And by calculating the amplitude difference and occurrence frequency between the peak and valley values of the thermal stress, the characteristic states of the furnace wall thermal stress changing with the thermal cycle conditions are determined. Finally, the response characteristics of the furnace wall thermal stress change are output.

[0022] Combine the temperature change characteristics on the furnace wall surface, the spatial distribution characteristics of the heat flux on the furnace wall surface, and the response characteristics of the furnace wall thermal stress change to obtain the comprehensive state parameters of the furnace body structural components under thermal cycle conditions; Combine the time-series change characteristics of the furnace wall surface temperature, the spatial distribution characteristics of the heat flux on the furnace wall surface, and the response characteristics of the furnace wall thermal stress change to establish the comprehensive state parameters of the furnace body structural components under thermal cycle conditions, including the temperature change characteristics, heat flux spatial distribution characteristics, and thermal stress response characteristics.

[0023] S2. Based on the comprehensive state parameters of the furnace body structural components under thermal cycle conditions, construct a thermal fatigue deformation trend estimation model and generate a geometric parameter drift vector used to evaluate the deformation of the furnace body structural components, including: Based on the comprehensive state parameters of the furnace body structural components under thermal cycle conditions, establish the correlation relationship between thermal expansion displacement and structural deformation; Select the temperature change characteristics on the furnace wall surface as the independent variable. The temperature change characteristics are composed of the time-series difference generated by the change of the furnace wall surface temperature with the thermal cycle conditions, expressed as the amplitude, frequency, and duration of the temperature change. Select the response characteristics of the furnace wall thermal stress change as another independent variable. The thermal stress response characteristics are composed of the peak value, valley value, and thermal stress change slope in the thermal stress response curve. Finally, select the spatial distribution characteristics of the heat flux on the furnace wall surface as the third independent variable. The spatial distribution characteristics are composed of the average value, peak value of the heat flux in each area of the furnace wall, and the degree of non-uniformity of the heat flux distribution.

[0024] Using a finite element analysis tool, a structural model of the furnace body structure components is constructed. The structural model is established based on the geometric parameters, material parameters, and initial state conditions in the design stage of the water jacket heating furnace. The structural model specifies the thermophysical properties, elastoplastic properties, and geometric dimensions of the furnace wall material to ensure that the structural model can truly reflect the actual structural characteristics and material properties of the furnace body structure components.

[0025] Taking the surface temperature change characteristics, thermal stress response characteristics, and heat flux spatial distribution characteristics of the furnace wall as the input to the structural model, a thermal-structural coupling analysis is carried out. The analysis process is as follows: Establish a mapping relationship between the furnace wall temperature change characteristics and the temperature field boundary conditions in the structural model, and convert the temperature amplitude change and duration represented by the temperature change characteristics into temperature load parameters in the structural model; Establish a mapping relationship between the furnace wall thermal stress response characteristics and the boundary constraint conditions applied in the structural model, that is, convert the difference and frequency of the thermal stress peak and valley into the corresponding structural load boundary conditions in the structural model; Establish a mapping relationship between the furnace wall heat flux spatial distribution characteristics and the heat flow boundary of the structural model, that is, convert the average value, peak value, and uneven degree of the heat flux distribution into the heat flow load parameters of the structural model.

[0026] After the mapping is completed and the boundary conditions of the structural model are applied, a thermal-structural coupling calculation of the structural model is carried out. A nonlinear finite element analysis algorithm is adopted, specifically selecting the elastoplastic deformation analysis and thermal-mechanical coupling solution method to accurately simulate and analyze the actual deformation state of the furnace body structure components.

[0027] After completing the thermal-structural coupling calculation of the furnace body structure components, a thermal expansion displacement distribution nephogram and a deformation state nephogram of the furnace body structure components are obtained. Correlation analysis is carried out on the calculated thermal expansion displacement data and deformation data. Specifically: Calculate the spatial correspondence relationship between the thermal expansion displacement field and the structural deformation field, that is, calculate the spatial comparison relationship between the thermal expansion displacement value and the structural deformation value at each specific position point on the furnace wall structural model, and obtain the spatial correlation relationship between the thermal expansion displacement and the structural deformation in the furnace wall structure components. Then, through numerical statistics, obtain the sensitivity degree and response law of the overall structural deformation of the furnace wall structure components to the thermal expansion displacement, so as to establish the correlation relationship between the thermal expansion displacement and the structural deformation of the furnace body structure components.

[0028] According to the correlation relationship between the thermal expansion displacement and the structural deformation, determine the thermal fatigue deformation trend of the furnace body structure components; Taking the multi-cycle thermal cycle conditions during the operation of the furnace wall structure components as the background, record the cumulative change of the structural deformation amount in each thermal cycle, and fit the cumulative change trend of the structural deformation amount into a thermal fatigue deformation trend curve. The thermal fatigue deformation trend curve reflects the structural deformation law of the furnace body structure components gradually accumulated with the thermal cycle conditions, and accurately describes the long-term change trend and change speed of the structural deformation.

[0029] Calculate the geometric shape drift of the furnace body structure components under thermal cycling conditions according to the thermal fatigue deformation trend of the furnace body structure components; Calculate the geometric shape drift of the furnace body structure components under thermal cycling conditions according to the thermal fatigue deformation trend curve. The calculation method is as follows: calculate the geometric difference between the initial design state and the deformed state of the furnace wall structure components after long-term thermal cycling conditions. The calculation process includes the spatial difference between the structural coordinate positions of each key part and each specific position of the furnace wall structure components after thermal cycling and the initial design structural coordinate positions. Calculate the cumulative change law of the spatial difference with the increase of the thermal cycling period to obtain the geometric shape drift of the furnace body structure components.

[0030] Generate a geometric parameter drift vector for evaluating the geometric parameter drift caused by the deformation of the furnace body structure components based on the geometric shape drift of the furnace body structure components under thermal cycling conditions; Vectorize the geometric shape drift of the furnace body structure components according to the position coordinates. Specifically, take the initial position of the structure as the coordinate origin, decompose the drift amount in the coordinate axis directions, that is, determine the projection sizes of the drift amount in the horizontal, vertical, and radial three orthogonal directions respectively; use the projection sizes to form a multi-dimensional geometric parameter drift vector. Each element in the drift vector represents the drift amount of a certain position of the furnace body structure components in a specific direction, indicating the geometric parameter drift characteristics of the furnace body structure components under multiple thermal cycling conditions.

[0031] S3. Conduct a structural equivalent mapping analysis on the geometric parameter drift vector to extract the contact thermal resistance change information in the heat transfer path, including: Adopt the structural equivalent mapping method to convert the geometric shape drift vector of the furnace body structure components under thermal cycling conditions into the structural change characteristics of the heat transfer path; Use the geometric parameter drift vector as the input parameter of the structural equivalent mapping method. The geometric parameter drift vector of the furnace body structure components contains the geometric drift direction and drift amplitude information of each specific position of the furnace body structure components under thermal cycling conditions, that is, including the displacement amounts in the horizontal, vertical, and radial directions. Map the drift amount at each specific position in the geometric parameter drift vector of the furnace body structure components to the corresponding heat transfer path position respectively to obtain the structural change characteristics of the heat conduction path of the furnace body structure components.

[0032] For example, the drift amount of the furnace body structure component in the vertical direction is mapped to the change in the interface position of the heat conduction path in the vertical direction, that is, the actual contact interface position in the vertical direction of the heat conduction path has a corresponding vertical offset; the drift amount of the furnace body structure component in the horizontal direction is mapped to the change in the contact interface position of the heat conduction path in the horizontal direction; the drift amount of the furnace body structure component in the radial direction is mapped to the change in the interface position of the heat conduction path in the radial direction. Through the above mapping process, the structural change characteristics of the heat conduction path of the furnace body structure component are formed, reflecting the structural change of the heat transfer path caused by the structural drift in each direction and position of the furnace body structure component.

[0033] According to the structural change characteristics of the heat transfer path, analyze the change law of the heat conduction interface contact area of the furnace body structure component, and determine the change amount of the heat conduction interface contact area; Adopt the interface reconstruction analysis method, using the structural change characteristics as input parameters, and calculate the influence of the relative offset of each interface position in the structural change characteristics on the heat conduction interface contact area.

[0034] Specifically, the interface reconstruction analysis method includes the following: determine the position and area size of the initial interface contact area; according to the offset amount of each position in the structural change characteristics, determine the specific offset value at each specific interface contact position; by analyzing the coincidence degree between the offset value of the interface contact position and the initial interface position, calculate the area of the interface coincidence area. The calculation method is: calculate the intersection area size of the initial interface contact area and the actual contact area after the interface offset, that is, the area of the interface coincidence area; calculate the change in the area of the interface coincidence area for each heat cycle period, and obtain the change law of the interface contact area within the heat cycle period.

[0035] After calculation by the above interface reconstruction analysis method, a change curve of the interface contact area with the change of heat cycle conditions is formed. The change curve expresses the trend and specific change speed of the change in the heat conduction interface contact area of the furnace body structure component gradually decreasing or increasing with the number of heat cycles, and finally determines the change amount of the heat conduction interface contact area.

[0036] According to the change amount of the heat conduction interface contact area, analyze the change of the heat transfer capacity of the heat contact interface of the furnace body structure component, and obtain the change information of the contact thermal resistance of the heat contact interface of the furnace body structure component; Specifically, the formula for calculating the heat conduction capacity of the thermal contact interface is given by the interface heat conduction theory: the heat conduction capacity of the thermal contact interface is equal to the product of the material thermal conductivity and the contact area of the thermal conduction interface divided by the thickness of the thermal resistance in contact with the thermal conduction interface. The formula for calculating the heat conduction capacity shows that the change in the contact area of the thermal conduction interface will directly affect the magnitude of the heat conduction capacity; the change amount of the contact area of the thermal conduction interface is input into the formula for calculating the heat conduction capacity to obtain the heat conduction capacity value under each state of the change in the contact area of the thermal conduction interface; finally, the change difference of the heat transfer capacity of the thermal contact interface relative to the initial state under each heat cycle is calculated, which is the change information of the thermal contact conduction capacity of the thermal contact interface.

[0037] According to the change information of the heat transfer capacity of the thermal contact interface of the furnace body structure components, the change information of the contact thermal resistance of the thermal contact interface of the furnace body structure components is obtained. Using the thermal resistance calculation method, the formula for calculating the contact thermal resistance is: the contact thermal resistance of the thermal contact interface is equal to the temperature difference across the thermal contact interface divided by the heat transfer capacity of the thermal contact interface. According to the change information of the thermal contact conduction capacity of the thermal contact interface and the actually measured temperature difference data of the thermal contact interface, the change amount of the contact thermal resistance of the thermal contact interface with the heat cycle is calculated. The contact thermal resistance increases as the heat conduction capacity of the thermal contact interface decreases.

[0038] S4. Based on the geometric parameter drift vector, reconstruct the heat flow field of the furnace wall to obtain the dynamic correction vector of the actual heat flow distribution, including: Based on the geometric shape drift vector generated by the furnace body structure components under thermal cycling conditions, establish a numerical reconstruction model of the heat flow field of the furnace wall; The geometric shape drift vector represents the drift amplitude and position of the furnace wall structure components in three directions: the horizontal direction, the vertical direction, and the radial direction. The establishment of the numerical reconstruction model of the heat flow field of the furnace wall adopts the numerical simulation method of heat flow field reconstruction. The finite difference method is selected for the numerical simulation method. The finite difference method divides the furnace wall surface into multiple continuous grid cells, forms a finite difference calculation grid by discretizing the furnace wall surface, and specifies the size and distribution density of the grid cells to ensure that the numerical simulation accuracy meets the requirements of actual engineering applications.

[0039] Based on the initial geometric structure of the furnace wall, an initial numerical model of the furnace wall heat flux field is established. The numerical model of the furnace wall heat flux field defines the initial heat flux boundary conditions on the furnace wall surface, that is, the heat flux density distribution of each grid cell on the furnace wall surface in the undeformed state. The boundary position of the initial numerical model of the furnace wall heat flux field is corrected by the geometric shape drift vector generated by the furnace structure components under thermal cycling conditions, that is, the positions of each grid cell in the horizontal, vertical, and radial directions are adjusted to form the corrected grid cell positions of the heat flux field. The correction method is as follows: the initial position of each grid cell on the furnace wall is respectively added with the amplitude of the drift vector at the corresponding position to realize the mapping transformation from the initial grid cell to the grid cell position in the actual drift state, and a numerical reconstruction model of the furnace wall heat flux field is established.

[0040] According to the numerical reconstruction model of the furnace wall heat flux field, the spatial distribution characteristics of the heat flux on the furnace wall surface are reconstructed and analyzed to obtain the variation law of the heat flux distribution on the furnace wall surface under the actual operating conditions of the furnace wall structure components; Taking the spatial distribution characteristics of the heat flux on the furnace wall surface as input data, which includes the average value, peak value, and non-uniformity degree of the heat flux at different spatial positions on the furnace wall surface; mapping the above input data to the corrected grid cell positions of the numerical reconstruction model of the furnace wall heat flux field to realize the mapping distribution of the spatial distribution characteristics of the heat flux on the furnace wall surface under actual drift conditions. The mapping method is as follows: the heat flux characteristic values on the furnace wall surface are assigned to each grid cell according to the actual position of the furnace wall grid cell; then, numerical solutions are carried out according to the heat conduction equation to calculate the heat conduction relationship between each grid cell, so as to obtain the actual distribution of the heat flux on the furnace wall surface at each corrected grid cell position.

[0041] In order to obtain the variation law of the heat flux distribution on the furnace wall surface under the actual operating conditions of the furnace wall structure components, a heat flux variation trend analysis method is adopted. The implementation process is as follows: calculate the difference value between the heat flux of each grid cell in the actual drift state and the heat flux of the corresponding grid cell in the initial state; statistically analyze the spatial distribution of the difference values on the entire furnace wall surface; then, based on the spatial distribution of the difference values, analyze the overall heat flux distribution variation trend of the furnace wall structure components, which is manifested as the regional distribution of the increase or decrease of the heat flux of each grid cell, the variation trend of the heat flux distribution non-uniformity, and the positions with significant heat flux changes in local areas, and finally form the variation law of the heat flux distribution on the furnace wall surface under the actual operating conditions of the furnace wall structure components.

[0042] According to the variation law of the heat flux distribution on the furnace wall surface, determine the position and heat flux variation amplitude of the local offset area of the heat flux of the furnace wall structure components, and generate a dynamic correction vector of the actual heat flow distribution of the furnace wall structure components; In the variation law of the heat flux distribution on the furnace wall surface, identify the local areas where there are significant changes in heat flux. The areas with significant heat flux changes are defined as the positions where the difference in heat flux from the initial state exceeds a preset threshold; the preset threshold is set as a multiple of the average change value of the overall heat flux of the furnace wall; locate the positions of the grid cells that exceed the preset threshold, so as to identify the positions of the local offset areas of the heat flux of the furnace wall structural components; at the same time, calculate the magnitude of the heat flux change in the local offset area, and the calculation method is the difference between the actual heat flux value in the local offset area of the heat flux and the heat flux value at the corresponding position in the initial state, and the positive or negative of the difference indicates the trend of heat flux increase or decrease.

[0043] Vectorially represent the positions of the local offset areas of the heat flux of the furnace wall structural components and the corresponding magnitudes of the heat flux changes, that is, taking the position coordinates of each local offset area on the furnace wall surface as the base points and the magnitude of the heat flux change as the vector length to form a dynamic heat flow correction vector; then arrange all the dynamic heat flow correction vectors according to the positions of the furnace wall geometric structure to form a dynamic correction vector of the actual heat flow distribution of the furnace wall structural components. The dynamic correction vector expresses the change characteristics of the actual heat flow distribution of the furnace wall structural components with the structural drift under the action of the thermal cycle working condition.

[0044] S5. Comprehensively analyze the information on the change of the contact thermal resistance in the heat transfer path and the dynamic correction vector of the actual heat flow distribution, construct a time-varying boundary condition inversion model, and output a multi-parameter coupled boundary correction set, including: Based on the information on the change of the contact thermal resistance of the thermal contact interface of the furnace body structural components, establish a time-varying model of the change of the contact thermal resistance parameters; The information on the change of the contact thermal resistance is the thermal resistance change characteristics caused by the geometric drift of the structure and the change of the interface contact area at the thermal contact interface of the furnace body structural components under long-term thermal cycle conditions, which is manifested as a trend curve of the cumulative change of the thermal resistance with the thermal cycle period. In order to reflect the dynamic process of the change of the contact thermal resistance with time, a time series analysis method is selected to construct a time-varying model of the change of the contact thermal resistance parameters.

[0045] The autoregressive moving average model is selected for the time series analysis method: taking the numerical values of the contact thermal resistance in each thermal cycle period in the information on the change of the contact thermal resistance as the basic data sequence, preprocess the basic data sequence, including detrending processing and standardization processing; select the order of the autoregressive moving average model, and the order selection is determined by the Bayesian information criterion; after determining the order, estimate the model parameters of the basic data sequence, and the parameter estimation method is the maximum likelihood estimation method to obtain a time-varying model of the change of the contact thermal resistance parameters with the thermal cycle period. The time-varying model is expressed as the numerical value of the contact thermal resistance in the current period is the sum of the products of the numerical values of the contact thermal resistance in several historical periods and the model estimation parameters plus the model error term. Among them, the model error term is manifested as a random disturbance term, indicating the influence of unpredictable small disturbances on the change of the thermal resistance.

[0046] Based on the dynamic correction vector of the actual heat flux distribution of the furnace wall structure components, a time-varying model of the dynamic correction parameters of heat flow is established; The dynamic correction vector of the actual heat flux distribution represents the heat flux change caused by geometric drift of the furnace wall structure components under the action of the thermal cycle condition, including the positions of the local offset regions on the furnace wall surface and the corresponding heat flux change amplitudes. In order to reflect the dynamic change process of the actual heat flux change with the thermal cycle period, a state space model is selected to establish a time-varying model of the dynamic correction parameters of heat flow.

[0047] The process of establishing the state space model is as follows: construct a state vector, which includes the heat flux correction amplitudes of each local region of the furnace wall; define a state transition equation, which is expressed as the current heat flux correction state vector being equal to the heat flux correction state vector of the previous cycle multiplied by the state transition matrix plus the process disturbance vector, and the state transition matrix is defined as a matrix reflecting the evolution of the heat flux correction state vector with the thermal cycle period, and the process disturbance vector is a random vector subject to Gaussian distribution; construct an observation equation, which is expressed as the actual observed heat flux distribution correction data vector being equal to the state vector multiplied by the observation matrix plus the observation noise vector; the observation matrix is defined as a matrix reflecting the mapping relationship between the actual measured heat flux data and the state vector, and the observation noise vector is the random error existing in the actual measurement process; finally, the Kalman filter algorithm is used for state space model parameter estimation to obtain a time-varying model of the dynamic correction parameters of heat flow of the furnace wall structure components.

[0048] According to the time-varying model of the contact thermal resistance parameter change and the time-varying model of the dynamic correction parameters of heat flow, a multi-parameter comprehensive analysis of the heat transfer boundary conditions of the furnace body structure components is carried out, and an inversion model of the boundary conditions for the dynamic change of the heat transfer characteristics of the furnace body structure components is constructed; A multi-parameter collaborative analysis method is adopted, specifically: synchronize the output values of the time-varying model of the contact thermal resistance parameter and the time-varying model of the dynamic correction parameters of heat flow; the synchronization method is to determine a unified time scale of the thermal cycle period, and align the output values of the time-varying model of the contact thermal resistance parameter and the time-varying model of the dynamic correction parameters of heat flow according to the unified period to form two parameter sequences; then, taking the parameter sequences as inputs, a parameter collaborative analysis matrix is established, and each element of the collaborative analysis matrix represents the mutual influence intensity between the contact thermal resistance parameter and the heat flux correction parameter within the corresponding thermal cycle period; the calculation method is to normalize the product of the changes of the two parameters within the same period; finally, through the eigenvalue analysis of the collaborative analysis matrix, a set of dynamic change characteristic parameters of the heat transfer boundary conditions of the furnace body structure components is obtained.

[0049] Based on the dynamically varying characteristic parameter set, an inverse model for the boundary conditions of the dynamic variation of the heat transfer characteristics of the furnace structure components is constructed. The inverse model for the boundary conditions adopts the inverse problem solving method: set the forward heat transfer model of the furnace wall structure components, and the forward heat transfer model expresses the calculation formula for the temperature field distribution of the furnace wall of the furnace structure components under specific heat transfer boundary condition parameters; with the dynamically varying characteristic parameter set as the input, construct an objective function, and the objective function is expressed as the sum of squares of the differences between the actually measured temperature field and the temperature field calculated by the forward model; use the genetic algorithm optimization method to solve the minimum value of the objective function to obtain the optimal estimation result of the heat transfer boundary condition parameters; define the initial population, crossover rate, mutation rate, and selection rule in the genetic algorithm, and use the fitness function to evaluate the optimization degree of each generation of the population to obtain the estimated value of the optimal heat transfer boundary condition parameters under actual conditions.

[0050] Through the inverse model for the boundary conditions, output the multi-parameter coupled boundary correction set for the heat transfer boundary conditions under the thermal cycle conditions of the furnace structure components; The multi-parameter coupled boundary correction set includes the contact thermal resistance correction parameter and the furnace wall heat flux distribution correction parameter, and is presented in the form of a multi-dimensional vector.

[0051] S6. Update the boundary conditions in the heat transfer simulation model of the water jacket heating furnace according to the multi-parameter coupled boundary correction set, including: Based on the multi-parameter coupled boundary correction set, determine the contact thermal resistance correction parameter and the furnace wall heat flux distribution correction parameter; Extract the contact thermal resistance correction parameter and the furnace wall heat flux distribution correction parameter from the multi-parameter coupled boundary correction set.

[0052] The contact thermal resistance correction parameter is expressed as the value after correction of the contact thermal resistance parameter at each position in the thermal contact interface, and the heat flux distribution correction parameter is expressed as the value after correction of the heat flux at each local position of the furnace wall. In order to ensure that the contact thermal resistance correction parameter and the furnace wall heat flux distribution correction parameter are applicable to the heat transfer simulation model of the water jacket heating furnace, a parameter mapping process is implemented. The parameter mapping method is: according to the initial design structure positions of each contact interface position and the furnace wall area in the heat transfer simulation model of the water jacket heating furnace, establish the mapping relationship between the contact thermal resistance correction parameter and the furnace wall heat flux distribution correction parameter and the positions in the heat transfer simulation model of the water jacket heating furnace, that is, map each correction parameter to the corresponding geometric position in the heat transfer simulation model of the water jacket heating furnace.

[0053] Update the contact thermal resistance boundary condition of the thermal contact interface in the heat transfer simulation model of the water jacket heating furnace with the contact thermal resistance correction parameter; A heat transfer simulation model of the water jacket heater is established. The heat transfer simulation model of the water jacket heater is constructed based on the initial geometric model, material property parameters, and initial heat transfer boundary conditions at the design stage of the water jacket heater. The numerical value of the initial contact thermal resistance of the thermal contact interface is defined. Specifically, the initial contact thermal resistance of the thermal contact interface is the thermal resistance value set according to the ideal contact conditions at the design stage.

[0054] The update method of the contact thermal resistance correction parameter is as follows: Replace and update each corresponding correction value in the contact thermal resistance correction parameter with the initial contact thermal resistance parameter of the thermal contact interface of the heat transfer simulation model of the water jacket heater one by one: Replace the contact thermal resistance value at each position in the boundary condition of the thermal contact interface of the water jacket heater simulation model with the corresponding contact thermal resistance correction parameter to ensure that the boundary condition of the thermal contact interface in the water jacket heater simulation model accurately represents the dynamic characteristics of the contact thermal resistance under the actual operating conditions.

[0055] Update the furnace wall heat flux distribution correction parameter to the spatial distribution boundary condition of the furnace wall surface heat flux in the heat transfer simulation model of the water jacket heater; The initial boundary condition of the spatial distribution of the furnace wall surface heat flux in the heat transfer simulation model of the water jacket heater is defined as the initial heat flux density distribution state set according to the ideal working conditions at the design stage, which is expressed as the initial heat flux density value of each local area on the furnace wall surface.

[0056] The update method is as follows: Perform spatial grid division on the boundary condition of the furnace wall surface heat flux of the simulation model, and divide the furnace wall surface into several local area grids; Each local area grid corresponds to an initial heat flux density value in the initial boundary condition. Map the furnace wall heat flux distribution correction parameter to the grid area on the furnace wall surface. The position mapping method is to determine the grid position corresponding to each heat flux correction value in the heat flux distribution correction parameter according to the initial position coordinates of each local area grid on the furnace wall surface; Replace the initial heat flux density value of each grid area on the furnace wall surface with the correction value at the corresponding position in the heat flux distribution correction parameter under the actual operating conditions. The update process ensures that the spatial distribution boundary condition of the furnace wall surface heat flux reflects the actual dynamic change state of the heat flow.

[0057] The updated contact thermal resistance boundary condition of the thermal contact interface and the spatial distribution boundary condition of the furnace wall surface heat flux form the current heat transfer boundary condition of the heat transfer simulation model of the water jacket heater; The current heat transfer boundary condition includes the updated boundary condition of the thermal contact interface and the boundary condition of the furnace wall surface heat flux. In order to verify the effectiveness and accuracy of the updated heat transfer simulation model of the water jacket heater, verification is required: The surface temperature distribution of the furnace wall and the thermal stress response data measured under the actual operating conditions are used as the comparison benchmark for the output results of the heat transfer simulation model of the water jacket heating furnace; the simulation process is to use the updated heat transfer boundary conditions to perform numerical simulation calculations of the furnace wall temperature field and thermal stress field; the finite element numerical solution method is adopted to perform the simulation solution of the furnace wall temperature field and thermal stress field; the simulation model solution adopts the thermal-structural coupling analysis method to ensure that the simulation results accurately reflect the actual thermal cycle conditions; the simulation results of the furnace wall surface temperature field and thermal stress field are obtained through calculation.

[0058] The simulation results are compared and analyzed with the actual measurement data, and the error between the simulation results and the actual measurement data is calculated: at each measurement position on the furnace wall surface, the differences between the temperature values and thermal stress values calculated by simulation and the actual measurement temperature values and thermal stress values are calculated respectively; statistical analysis is performed on all the differences to obtain the distribution of the errors. According to the results of the error statistical analysis, the accuracy of the updated heat transfer simulation model boundary conditions is evaluated; if the error meets the preset error threshold, it indicates that the updated heat transfer boundary conditions can accurately describe the furnace wall heat transfer state under the actual operating conditions. The error threshold comes from the allowable deviation standard of the actual engineering application of the water jacket heating furnace. The determination method of the error threshold is: based on the maximum temperature and thermal stress deviations that the furnace wall structural components can accept under safe operating conditions, it is given based on the operating specifications, design specifications, performance index requirements and long-term operating experience data of the water jacket heating furnace.

[0059] Embodiment 2 The difference between Embodiment 2 and Embodiment 1 of the present invention is that this embodiment introduces a self-correcting and optimizing system for the heat transfer boundary conditions of a water jacket heating furnace.

[0060] Figure 2 The structural schematic diagram of a self-correcting and optimizing system for the heat transfer boundary conditions of a water jacket heating furnace according to the present invention is given. A self-correcting and optimizing system for the heat transfer boundary conditions of a water jacket heating furnace includes: A state acquisition module, which acquires the furnace wall surface temperature sequence, the furnace wall surface heat flux distribution information and the furnace wall thermal stress response curve of the furnace body structural components under the thermal cycle conditions, and obtains the comprehensive state parameters of the furnace body structural components under the thermal cycle conditions; A trend estimation module, which constructs a thermal fatigue deformation trend estimation model based on the comprehensive state parameters of the furnace body structural components under the thermal cycle conditions, and generates a geometric parameter drift vector used to evaluate the deformation of the furnace body structural components; A thermal resistance analysis module, which performs a structural equivalent mapping analysis on the geometric parameter drift vector and extracts the contact thermal resistance change information in the heat transfer path; A heat flow reconstruction module, which reconstructs the furnace wall heat flow field based on the geometric parameter drift vector and obtains a dynamic correction vector of the actual heat flow distribution; The boundary inversion module comprehensively analyzes the information on the change of contact thermal resistance in the heat transfer path and the dynamic correction vector of the actual heat flux distribution, constructs a time-varying boundary condition inversion model, and outputs a multi-parameter coupled boundary correction set. The model update module updates the boundary conditions in the heat transfer simulation model of the water jacket furnace according to the multi-parameter coupled boundary correction set.

[0061] The above formulas are all dimensionless and take their numerical values for calculation. The formulas are obtained by collecting a large amount of data for software simulation to obtain a formula that is closest to the actual situation. The preset parameters and threshold selection in the formulas are set by those skilled in the art according to the actual situation.

[0062] The above embodiments can be implemented in whole or in part by software, hardware, firmware, or any other combination. When implemented using software, the above embodiments can be implemented in whole or in part in the form of a computer program product. The computer program product includes one or more computer instructions or computer programs. When the computer instructions or computer programs are loaded or executed on a computer, the processes or functions described in the embodiments of the present application are generated in whole or in part. The computer can be a general-purpose computer, a special-purpose computer, a computer network, or other programmable devices. The computer instructions can be stored in a computer-readable storage medium, or transmitted from one computer-readable storage medium to another computer-readable storage medium. For example, the computer instructions can be transmitted from one website, computer, server, or data center to another website, computer, server, or data center by wire (such as infrared, wireless, microwave, etc.). The computer-readable storage medium can be any available medium that the computer can access, or a data storage device such as a server or data center that includes one or more collections of available media. The available media can be magnetic media (such as floppy disks, hard disks, magnetic tapes), optical media (such as DVDs), or semiconductor media. The semiconductor media can be a solid-state drive.

[0063] Those of ordinary skill in the art can realize that the modules and algorithm steps of the examples described in combination with the embodiments disclosed herein can be implemented by electronic hardware, or a combination of computer software and electronic hardware. Whether these functions are executed in a hardware or software manner depends on the specific application and design constraints of the technical solution. Skilled professionals can use different methods to implement the described functions for each specific application, but such implementation should not be considered to exceed the scope of the present application.

[0064] Those skilled in the art can clearly understand that for the convenience and brevity of description, the specific working processes of the systems, devices, and modules described above can refer to the corresponding processes in the foregoing method embodiments, and will not be described herein again.

[0065] In several embodiments provided in the present application, it should be understood that the disclosed systems, devices, and methods can be implemented in other ways. For example, the device embodiments described above are merely illustrative. For example, the division of the modules is only a logical function division. In actual implementation, there can be other division methods. For example, multiple modules or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the displayed or discussed coupling or direct coupling or communication connection between each other can be through some interfaces. The indirect coupling or communication connection of the devices or modules can be in electrical, mechanical, or other forms.

[0066] The modules described as separate components may or may not be physically separated. The components displayed as modules may or may not be physical modules. They can be located in one place or distributed to multiple network modules. Some or all of the modules can be selected according to actual needs to achieve the purpose of the solution of this embodiment.

[0067] In addition, in each embodiment of the present application, the functional modules can be integrated into one processing module, or each module can exist physically alone, or two or more modules can be integrated into one module.

[0068] If the above functions are implemented in the form of software function modules and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on such an understanding, the technical solution of the present application, in essence, or the part that contributes to the prior art, or a part of this technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to enable a computer device (which can be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the methods described in each embodiment of the present application. The foregoing storage medium includes: various media such as USB flash drives, mobile hard disks, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical discs that can store program codes.

[0069] As described above, the above is only the specific implementation manner of the present application, but the protection scope of the present application is not limited thereto. Any person skilled in the art within the technical scope disclosed in the present application can easily think of changes or substitutions, which should all be covered within the protection scope of the present application. Therefore, the protection scope of the present application should be subject to the protection scope of the claims.

[0070] Finally, the above are only the preferred embodiments of the present invention and are not intended to limit the present invention. Any modifications, equivalent replacements, improvements, etc. made within the spirit and principles of the present invention shall be included within the protection scope of the present invention.

Claims

1. A self-correcting and optimizing method for the heat transfer boundary conditions of a water jacket heating furnace, characterized in that, It includes the following steps: S1. Collect the furnace wall surface temperature sequence, furnace wall surface heat flux distribution information, and furnace wall thermal stress response curve of the furnace body structure components under thermal cycling conditions to obtain the comprehensive state parameters of the furnace body structure components under thermal cycling conditions; S2. Based on the comprehensive state parameters of the furnace body structure components under thermal cycling conditions, construct a thermal fatigue deformation trend estimation model and generate a geometric parameter drift vector for evaluating the deformation of the furnace body structure components; S3. Conduct a structural equivalent mapping analysis on the geometric parameter drift vector to extract the contact thermal resistance change information in the heat transfer path; S4. Based on the geometric parameter drift vector, reconstruct the furnace wall heat flow field to obtain a dynamic correction vector of the actual heat flow distribution; S5. Conduct a comprehensive analysis of the contact thermal resistance change information in the heat transfer path and the dynamic correction vector of the actual heat flow distribution, construct a time-varying boundary condition inversion model, and output a multi-parameter coupled boundary correction set; S6. Update the boundary conditions in the heat transfer simulation model of the water jacket heating furnace according to the multi-parameter coupled boundary correction set.

2. A self-correction optimization method for the heat transfer boundary condition of a water jacket heating furnace according to claim 1, characterized in that S1 is specifically as follows: Collect the furnace wall surface temperature sequence, furnace wall surface heat flux distribution information, and furnace wall thermal stress response curve of the furnace body structure components under long-term thermal cycling conditions during the operation of the water jacket heating furnace; Conduct a time-series change analysis on the furnace wall surface temperature sequence to obtain the temperature change characteristics of the furnace wall surface; Conduct a spatial distribution analysis on the furnace wall surface heat flux distribution information to obtain the spatial distribution characteristics of the furnace wall surface heat flux; Conduct a feature recognition on the furnace wall thermal stress response curve to obtain the response characteristics of the furnace wall thermal stress change; Combine the temperature change characteristics of the furnace wall surface, the spatial distribution characteristics of the furnace wall surface heat flux, and the response characteristics of the furnace wall thermal stress change to obtain the comprehensive state parameters of the furnace body structure components under thermal cycling conditions.

3. A method for self-correcting and optimizing the heat transfer boundary condition of a water jacket heating furnace according to claim 2, characterized in that S2 is specifically as follows: Based on the comprehensive state parameters of the furnace body structure components under thermal cycling conditions, establish the correlation between thermal expansion displacement and structural deformation; According to the correlation between thermal expansion displacement and structural deformation, determine the thermal fatigue deformation trend of the furnace body structure components; According to the thermal fatigue deformation trend of the furnace body structure components, calculate the geometric shape drift amount generated by the furnace body structure components under thermal cycling conditions; Based on the geometric shape drift amount generated by the furnace body structure components under thermal cycling conditions, generate a geometric parameter drift vector for evaluating the deformation of the furnace body structure components.

4. A method for self-correcting and optimizing the heat transfer boundary conditions of a water jacket heating furnace according to claim 3, characterized in that, S3 is specifically as follows: Adopt the structural equivalent mapping method to convert the geometric shape drift vector generated by the furnace body structure components under thermal cycling conditions into the structural change characteristics of the heat transfer path; According to the structural change characteristics of the heat transfer path, analyze the change law of the heat conduction interface contact area of the furnace body structure components, and determine the change amount of the heat conduction interface contact area; According to the change amount of the heat conduction interface contact area, analyze the change of the heat transfer capacity of the heat contact interface of the furnace body structure components to obtain the contact thermal resistance change information of the heat contact interface of the furnace body structure components.

5. A method for self-correcting and optimizing the heat transfer boundary condition of a water jacket heating furnace according to claim 4, characterized in that S4 is specifically as follows: Based on the geometric shape drift vector generated by the furnace body structure components under thermal cycling conditions, establish a numerical reconstruction model of the furnace wall heat flow field; Based on the numerical reconstruction model of the furnace wall heat flux field, the spatial distribution characteristics of the furnace wall surface heat flux are reconstructed and analyzed to obtain the variation law of the furnace wall surface heat flux distribution under the actual operating conditions of the furnace wall structural components; According to the variation law of the furnace wall surface heat flux distribution, determine the position of the local offset area of the heat flux of the furnace wall structural components and the variation range of the heat flux, and generate a dynamic correction vector of the actual heat flux distribution of the furnace wall structural components.

6. A method for self-correcting and optimizing the heat transfer boundary conditions of a water jacket heating furnace according to claim 5, characterized in that, S5, specifically: Based on the information on the change of the contact thermal resistance at the thermal contact interface of the furnace body structural components, establish a time-varying model of the change of the contact thermal resistance parameters; Based on the dynamic correction vector of the actual heat flux distribution of the furnace wall structural components, establish a time-varying model of the dynamic correction parameters of the heat flow; According to the time-varying model of the change of the contact thermal resistance parameters and the time-varying model of the dynamic correction parameters of the heat flow, conduct a multi-parameter comprehensive analysis of the heat transfer boundary conditions of the furnace body structural components, and construct an inversion model of the boundary conditions for the dynamic change of the heat transfer characteristics of the furnace body structural components; Through the boundary condition inversion model, output a multi-parameter coupling boundary correction set for the heat transfer boundary conditions under the thermal cycle conditions of the furnace body structural components.

7. A method for self-correcting and optimizing the heat transfer boundary conditions of a water jacket heating furnace according to claim 6, characterized in that, S6, specifically: Based on the multi-parameter coupling boundary correction set, determine the contact thermal resistance correction parameters and the furnace wall heat flux distribution correction parameters; Update the contact thermal resistance correction parameters to the contact thermal resistance boundary conditions of the thermal contact interface in the heat transfer simulation model of the water jacket heating furnace; Update the furnace wall heat flux distribution correction parameters to the spatial distribution boundary conditions of the furnace wall surface heat flux in the heat transfer simulation model of the water jacket heating furnace; The updated contact thermal resistance boundary conditions of the thermal contact interface and the spatial distribution boundary conditions of the furnace wall surface heat flux form the current heat transfer boundary conditions of the heat transfer simulation model of the water jacket heating furnace.

8. A heat transfer boundary condition self-correction and optimization system for a water jacket heating furnace, which is used to implement the heat transfer boundary condition self-correction and optimization method for a water jacket heating furnace according to any one of claims 1-7, characterized in that, Including: A state acquisition module that acquires the furnace wall surface temperature sequence, the furnace wall surface heat flux distribution information, and the furnace wall thermal stress response curve of the furnace body structural components under thermal cycle conditions, and obtains the comprehensive state parameters of the furnace body structural components under thermal cycle conditions; A trend estimation module that constructs a thermal fatigue deformation trend estimation model based on the comprehensive state parameters of the furnace body structural components under thermal cycle conditions, and generates a geometric parameter drift vector used to evaluate the deformation of the furnace body structural components; A thermal resistance analysis module that conducts a structural equivalent mapping analysis of the geometric parameter drift vector and extracts the information on the change of the contact thermal resistance in the heat transfer path; A heat flux reconstruction module that reconstructs the furnace wall heat flux field based on the geometric parameter drift vector and obtains a dynamic correction vector of the actual heat flux distribution; A boundary inversion module that comprehensively analyzes the information on the change of the contact thermal resistance in the heat transfer path and the dynamic correction vector of the actual heat flux distribution, constructs a time-varying boundary condition inversion model, and outputs a multi-parameter coupling boundary correction set; A model update module that updates the boundary conditions in the heat transfer simulation model of the water jacket heating furnace according to the multi-parameter coupling boundary correction set.

Images