A self-correction optimization method and system for heat transfer boundary conditions of a water-jacketed heating furnace

Abstract

The present invention discloses a self-correction optimization method and system for the heat transfer boundary conditions of a water-jacketed heating furnace, which specifically relates to the technical field of condition optimization control. The method is to construct a thermal fatigue deformation trend estimation model by collecting the furnace wall surface temperature sequence of the furnace body structural components, the furnace wall surface heat flux distribution information and the furnace wall thermal stress response curve, and generate a geometric parameter drift vector; perform structural equivalent mapping analysis on the geometric parameter drift vector to extract contact thermal resistance change information; perform furnace wall thermal flow field reconstruction on the geometric parameter drift vector to obtain a dynamic correction vector; perform comprehensive analysis on the contact thermal resistance change information and the heat flow correction vector, construct a time-varying boundary condition inversion model, output a multi-parameter coupled boundary correction set, and use it to update the boundary conditions in the water-jacketed heating furnace heat transfer simulation model, realize the identification and adaptive correction of the boundary conditions, and improve the heat transfer control accuracy and thermal efficiency.

Description

Technical Field

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

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

[0003] Existing self-correcting optimization methods are generally based on the calculation of initial design geometric parameters, and lack real-time perception and dynamic response to structural deformation and changes in boundary conditions during operation. Summary of the Invention

[0004] In order to overcome the above-mentioned defects of the prior art, an embodiment of the present invention provides a method and system for self-correction optimization of heat transfer boundary conditions of a water-jacketed heating furnace to solve the problems raised in the above-mentioned background technology.

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

[0006] A self-correction optimization method for heat transfer boundary conditions of a water-jacket heating furnace comprises the following steps:

[0007] S1. Collecting 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 cycle conditions to obtain the comprehensive state parameters of the furnace body structure components under thermal cycle conditions;

[0008] S2. Based on the comprehensive state parameters of the furnace structure components under thermal cycle conditions, a thermal fatigue deformation trend estimation model is constructed, and a geometric parameter drift vector is generated for evaluating the deformation of the furnace structure components;

[0009] S3. Perform structural equivalent mapping analysis on the geometric parameter drift vector to extract the contact thermal resistance change information in the heat transfer path;

[0010] S4. Reconstruct the furnace wall heat flow field based on the geometric parameter drift vector to obtain the dynamic correction vector of the actual heat flow distribution;

[0011] 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, build a time-varying boundary condition inversion model, and output a multi-parameter coupled boundary correction set;

[0012] S6. Update the boundary conditions in the water jacket heating furnace heat transfer simulation model according to the multi-parameter coupled boundary correction set.

[0013] In a preferred embodiment, S1 is specifically:

[0014] 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 cycle conditions during the operation of the water jacket heating furnace;

[0015] The temperature variation characteristics of the furnace wall surface are obtained by analyzing the time series variation of the furnace wall surface temperature.

[0016] Perform spatial distribution analysis on the heat flux distribution information on the furnace wall surface to obtain the heat flux spatial distribution characteristics on the furnace wall surface;

[0017] Perform feature recognition on the furnace wall thermal stress response curve to obtain the response characteristics of the furnace wall thermal stress change;

[0018] 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 changes of the furnace wall are combined to obtain the comprehensive state parameters of the furnace structure components under thermal cycle conditions.

[0019] In a preferred embodiment, S2 is specifically:

[0020] Based on the comprehensive state parameters of the furnace structure components under thermal cycle conditions, the correlation between thermal expansion displacement and structural deformation is established;

[0021] Determine the thermal fatigue deformation trend of the furnace structure components based on the correlation between thermal expansion displacement and structural deformation;

[0022] According to the thermal fatigue deformation trend of the furnace structure components, calculate the geometric shape drift of the furnace structure components under thermal cycle conditions;

[0023] Based on the geometric shape drift of the furnace structure components under thermal cycle conditions, a geometric parameter drift vector is generated for evaluating the deformation of the furnace structure components.

[0024] In a preferred embodiment, S3 is specifically:

[0025] The structural equivalent mapping method is used to convert the geometric shape drift vectors of the furnace structure components under thermal cycling conditions into the structural change characteristics of the heat transfer path;

[0026] According to the structural change characteristics of the heat transfer path, the change law of the heat conduction interface contact area of ​​the furnace structure components is analyzed to determine the change amount of the heat conduction interface contact area;

[0027] According to the change in the contact area of ​​the heat conduction interface, the change in the heat transfer capacity of the thermal contact interface of the furnace structure component is analyzed, and the contact thermal resistance change information of the thermal contact interface of the furnace structure component is obtained.

[0028] In a preferred embodiment, S4 is specifically:

[0029] Based on the geometric shape drift vectors generated by the furnace structure components under thermal cycling conditions, a numerical reconstruction model of the furnace wall thermal flow field is established;

[0030] Based on the numerical reconstruction model of the furnace wall thermal flow field, the spatial distribution characteristics of the heat flux on the furnace wall surface are reconstructed and analyzed to obtain the distribution variation law of the heat flux on the furnace wall surface under the actual operating conditions of the furnace wall structural components;

[0031] According to the variation law of heat flux distribution on the furnace wall surface, the position of the local heat flux offset area and the heat flux variation amplitude of the furnace wall structural component are determined, and the dynamic correction vector of the actual heat flux distribution of the furnace wall structural component is generated.

[0032] In a preferred embodiment, S5 is specifically:

[0033] Based on the contact thermal resistance variation information of the thermal contact interface of the furnace structure components, a time-varying model of the contact thermal resistance parameter variation is established;

[0034] Based on the dynamic correction vector of the actual heat flux distribution of the furnace wall structural components, a time-varying model of the heat flux dynamic correction parameters is established;

[0035] Based on the time-varying model of the contact thermal resistance parameters and the time-varying model of the heat flow dynamic correction parameters, a multi-parameter comprehensive analysis of the heat transfer boundary conditions of the furnace structure components is conducted, and an inversion model of the boundary conditions of the dynamic changes in the heat transfer characteristics of the furnace structure components is constructed;

[0036] Through the boundary condition inversion model, a multi-parameter coupled boundary correction set of heat transfer boundary conditions under thermal cycle conditions of furnace structure components is output.

[0037] In a preferred embodiment, S6 is specifically:

[0038] Based on the multi-parameter coupled boundary correction set, the correction parameters of contact thermal resistance and furnace wall heat flux distribution are determined;

[0039] Update the contact thermal resistance correction parameters to the thermal contact interface contact thermal resistance boundary conditions in the water jacket heating furnace heat transfer simulation model;

[0040] The furnace wall heat flux distribution correction parameters are updated to the furnace wall surface heat flux spatial distribution boundary conditions in the water jacket heating furnace heat transfer simulation model;

[0041] The updated thermal contact interface thermal resistance boundary conditions and the furnace wall surface heat flux spatial distribution boundary conditions form the current heat transfer boundary conditions of the water jacket heating furnace heat transfer simulation model.

[0042] In another aspect, the present invention provides a water-jacket heating furnace heat transfer boundary condition self-correction optimization system, comprising:

[0043] The state acquisition module collects 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 cycle conditions to obtain the comprehensive state parameters of the furnace body structure components under thermal cycle conditions;

[0044] The trend estimation module builds a thermal fatigue deformation trend estimation model based on the comprehensive state parameters of the furnace structure components under thermal cycle conditions, and generates a geometric parameter drift vector for evaluating the deformation of the furnace structure components;

[0045] Thermal resistance analysis module, performs structural equivalent mapping analysis on geometric parameter drift vectors and extracts contact thermal resistance change information in the heat transfer path;

[0046] The heat flow reconstruction module reconstructs the furnace wall heat flow field based on the geometric parameter drift vector and obtains the dynamic correction vector of the actual heat flow distribution;

[0047] The boundary inversion module 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 coupled boundary correction set;

[0048] The model update module updates the boundary conditions in the water jacket heating furnace heat transfer simulation model according to the multi-parameter coupling boundary correction set.

[0049] The technical effects and advantages of the self-correction optimization method and system for heat transfer boundary conditions of a water jacket heating furnace of the present invention are as follows:

[0050] By collecting the furnace wall surface temperature series, heat flux distribution information and thermal stress response curve, the operating status of the furnace structure under thermal cycle conditions can be fully perceived, thereby improving the accuracy of status monitoring. By constructing a thermal fatigue deformation trend estimation model, the cumulative drift characteristics of the geometric structure can be timely identified, effectively revealing the relationship between structural deformation and boundary changes. Through structural equivalent mapping analysis and furnace wall thermal flow field reconstruction, the changing characteristics of contact thermal resistance and heat flux distribution can be accurately extracted, enhancing the ability to quantify boundary condition changes. By constructing a time-varying boundary condition inversion model, the heat transfer boundary conditions can be accurately deduced. Through a multi-parameter coupled boundary correction set, the boundary conditions in the water-jacket heating furnace heat transfer simulation model can be updated, effectively avoiding local hot spots and efficiency reduction problems caused by model mismatch, and improving the simulation accuracy of the heating process, optimizing stability and operational safety. BRIEF DESCRIPTION OF THE DRAWINGS

[0051] Figure 1 A schematic diagram of a self-correction optimization method for heat transfer boundary conditions of a water-jacket heating furnace according to the present invention;

[0052] Figure 2 The present invention is a structural schematic diagram of a water-jacket heating furnace heat transfer boundary condition self-correction optimization system. DETAILED DESCRIPTION

[0053] The following will provide a clear and complete description of the technical solutions in the embodiments of the present invention in conjunction with the accompanying drawings. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. All other embodiments obtained by ordinary technicians in this field based on the embodiments of the present invention without making any creative efforts shall fall within the scope of protection of the present invention.

[0054] Example 1

[0055] Figure 1 A method for self-correcting and optimizing heat transfer boundary conditions of a water-jacketed heating furnace is provided, which comprises the following steps:

[0056] S1. Collecting 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 cycle conditions to obtain the comprehensive state parameters of the furnace body structure components under thermal cycle conditions;

[0057] S2. Based on the comprehensive state parameters of the furnace structure components under thermal cycle conditions, a thermal fatigue deformation trend estimation model is constructed, and a geometric parameter drift vector is generated for evaluating the deformation of the furnace structure components;

[0058] S3. Perform structural equivalent mapping analysis on the geometric parameter drift vector to extract the contact thermal resistance change information in the heat transfer path;

[0059] S4. Reconstruct the furnace wall heat flow field based on the geometric parameter drift vector to obtain the dynamic correction vector of the actual heat flow distribution;

[0060] 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, build a time-varying boundary condition inversion model, and output a multi-parameter coupled boundary correction set;

[0061] S6. Update the boundary conditions in the water jacket heating furnace heat transfer simulation model according to the multi-parameter coupled boundary correction set.

[0062] 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 structure components under thermal cycle conditions to obtain the comprehensive state parameters of the furnace structure components under thermal cycle conditions, including:

[0063] 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 cycle conditions during the operation of the water jacket heating furnace;

[0064] A temperature measuring sensor installed on the furnace wall surface of the water-jacketed heating furnace body structure assembly is used to continuously collect the temperature changes of the furnace wall surface of the furnace body structure assembly under long-term thermal cycle conditions in real time, and the collected continuous temperature data are arranged in sequence at fixed time intervals to form a furnace wall surface temperature sequence; a heat flux measuring device distributed at different spatial positions on the furnace wall surface is used to collect the heat flux size of the furnace wall surface in real time, and the position coordinate information is corresponded to form the heat flux distribution information of the furnace wall surface; at the same time, a thermal stress detection device is used to monitor and record the thermal stress changes of the furnace body structure assembly during the thermal cycle operation in real time, and the recorded data is plotted as a response curve of thermal stress changes over time.

[0065] The temperature variation characteristics of the furnace wall surface are obtained by analyzing the time series variation of the furnace wall surface temperature.

[0066] In order to obtain the temperature change characteristics in the furnace wall surface temperature sequence, a time series analysis algorithm is adopted. The furnace wall surface temperature sequence is used as input. The temperature value at each moment in the furnace wall surface temperature sequence is calculated with the temperature value at the previous adjacent time. The trend of the difference over time is also calculated. The frequency, amplitude and duration of the difference change trend are statistically analyzed to determine the change law of the furnace wall surface temperature. Finally, the time series change characteristics of the furnace wall surface temperature are output.

[0067] Perform spatial distribution analysis on the heat flux distribution information on the furnace wall surface to obtain the heat flux spatial distribution characteristics on the furnace wall surface;

[0068] In order to obtain the spatial distribution characteristics of the heat flux on the furnace wall surface, the 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 statistics of the heat flux data in each divided area, the average value, peak value and unevenness of the heat flux distribution in each area (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.

[0069] Perform feature recognition on the furnace wall thermal stress response curve to obtain the response characteristics of the furnace wall thermal stress change;

[0070] In order to obtain the response characteristics of the thermal stress changes in the furnace wall, the feature recognition method is used to process the thermal stress response curve of the furnace wall. By identifying the peak values, valley values ​​and thermal stress change slope in the thermal stress response curve, the typical state of the thermal stress change process is determined. By calculating the amplitude difference and occurrence frequency between the thermal stress peaks and valley values, the characteristic state of the furnace wall thermal stress changing with the thermal cycle conditions is determined, and finally the response characteristics of the furnace wall thermal stress change are output.

[0071] The temperature variation 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 are combined to obtain the comprehensive state parameters of the furnace structure components under thermal cycle conditions;

[0072] The temporal variation characteristics of the furnace wall surface temperature, the spatial distribution characteristics of the furnace wall surface heat flux and the response characteristics of the furnace wall thermal stress change are combined to establish the comprehensive state parameters of the furnace structure components under thermal cycle conditions, which include temperature variation characteristics, heat flux spatial distribution characteristics and thermal stress response characteristics.

[0073] S2. Based on the comprehensive state parameters of the furnace structure components under thermal cycling conditions, a thermal fatigue deformation trend estimation model is constructed, and a geometric parameter drift vector is generated for evaluating the deformation of the furnace structure components, including:

[0074] Based on the comprehensive state parameters of the furnace structure components under thermal cycle conditions, the correlation between thermal expansion displacement and structural deformation is established;

[0075] The temperature variation characteristics of the furnace wall surface were selected as the independent variable. The temperature variation characteristics are composed of the time series differences in the furnace wall surface temperature as the thermal cycle conditions change, expressed as the amplitude, frequency, and duration of the temperature variation. The response characteristics of the furnace wall thermal stress variation were selected as another independent variable. The thermal stress response characteristics are composed of the peak and valley values ​​in the thermal stress response curve and the slope of the thermal stress variation. Finally, the spatial distribution characteristics of the heat flux on the furnace wall surface were selected as the third independent variable. The spatial distribution characteristics are composed of the average value, peak value, and unevenness of the heat flux distribution in each area of ​​the furnace wall.

[0076] Finite element analysis tools were used to construct a structural model of the furnace's structural components. This model was established based on the geometric parameters, material parameters, and initial conditions used during the water-jacketed furnace design phase. The structural model specified the thermophysical properties, elastic-plastic properties, and geometric dimensions of the furnace wall material, ensuring that the model accurately reflects the actual structural characteristics and material properties of the furnace's structural components.

[0077] A thermal-structural coupling analysis was conducted using the furnace wall surface temperature variation characteristics, thermal stress response characteristics, and heat flux spatial distribution characteristics as input to the structural model. The analysis process is as follows: a mapping relationship is established between the furnace wall temperature variation characteristics and the temperature field boundary conditions in the structural model, converting the temperature amplitude change and duration represented by the temperature variation characteristics into temperature load parameters in the structural model; a mapping relationship is established between the furnace wall thermal stress response characteristics and the boundary constraints imposed in the structural model, that is, converting the difference and frequency of the thermal stress peak and valley values ​​into the corresponding structural load boundary conditions in the structural model; and a mapping relationship is established between the furnace wall heat flux spatial distribution characteristics and the heat flow boundary of the structural model, that is, converting the average value, peak value, and uneven distribution of the heat flux into the heat flow load parameters of the structural model.

[0078] After mapping and applying the structural model boundary conditions, thermal-structural coupling calculations are performed on the structural model. Using a nonlinear finite element analysis algorithm, specifically using elastic-plastic deformation analysis and thermal-mechanical coupling solutions, the actual deformation state of the furnace structure components is accurately simulated and analyzed.

[0079] After completing the thermal-structural coupling calculation of the furnace structure components, a cloud diagram of the thermal expansion displacement distribution and deformation state of the furnace structure components is obtained. A correlation analysis is performed on the calculated thermal expansion displacement data and deformation data. Specifically, the spatial correspondence between the thermal expansion displacement field and the structural deformation field is calculated. That is, the spatial comparison relationship between the thermal expansion displacement value and the structural deformation value at each specific location on the furnace wall structure model is calculated to obtain the spatial correlation between the thermal expansion displacement and structural deformation in the furnace wall structure components. Numerical statistics are then used to determine the sensitivity and response pattern of the overall structural deformation of the furnace wall structure components to thermal expansion displacement, thereby establishing a correlation between the thermal expansion displacement and structural deformation of the furnace structure components.

[0080] Determine the thermal fatigue deformation trend of the furnace structure components based on the correlation between thermal expansion displacement and structural deformation;

[0081] Using the multi-cycle thermal cycling conditions of the furnace wall structural components during operation as the background, the cumulative changes in structural deformation within each thermal cycle were recorded, and the cumulative change trend of the structural deformation was fitted into a thermal fatigue deformation trend curve. The thermal fatigue deformation trend curve reflects the gradual accumulation of structural deformation of the furnace structure components under thermal cycling conditions and accurately describes the long-term trend and rate of change of structural deformation.

[0082] According to the thermal fatigue deformation trend of the furnace structure components, calculate the geometric shape drift of the furnace structure components under thermal cycle conditions;

[0083] Based on the thermal fatigue deformation trend curve, the geometric drift of the furnace wall components under thermal cycling conditions is calculated. This calculation method involves calculating the geometric differences between the initial design state of the furnace wall components and their deformed state after long-term thermal cycling. The calculation process includes the spatial differences between the structural coordinate positions of key parts and specific locations of the furnace wall components after thermal cycling and the initial design coordinate positions. The cumulative change in these spatial differences with increasing thermal cycling cycles is calculated to determine the geometric drift of the furnace wall components.

[0084] Based on the geometric shape drift of the furnace structure components under thermal cycle conditions, a geometric parameter drift vector for evaluating the deformation of the furnace structure components is generated;

[0085] The geometric drift of the furnace structure components is vectorized according to positional coordinates. Specifically, the drift is decomposed along the coordinate axes, taking the initial position of the structure as the coordinate origin. This means that the projected magnitudes of the drift in the horizontal, vertical, and radial directions are determined. These projected magnitudes then form a multidimensional geometric parameter drift vector. Each element in the drift vector represents the drift of a specific position in a specific direction, characterizing the geometric parameter drift characteristics of the furnace structure components under multiple thermal cycles.

[0086] S3. Perform structural equivalent mapping analysis on the geometric parameter drift vector to extract the contact thermal resistance change information in the heat transfer path, including:

[0087] The structural equivalent mapping method is used to convert the geometric shape drift vectors of the furnace structure components under thermal cycling conditions into the structural change characteristics of the heat transfer path;

[0088] The geometric parameter drift vector is used as the input parameter of the structural equivalent mapping method. The geometric parameter drift vector of the furnace structure component contains the geometric drift direction and drift amplitude information at each specific location of the furnace structure component under thermal cycling conditions, namely, the displacement in the horizontal, vertical, and radial directions. The drift value at each specific location in the geometric parameter drift vector of the furnace structure component is mapped to the corresponding heat transfer path position, thereby obtaining the structural change characteristics of the heat conduction path of the furnace structure component.

[0089] For example, the vertical drift of the furnace structure is mapped to a 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 experiences a corresponding vertical offset. The horizontal drift of the furnace structure is mapped to a change in the contact interface position of the heat conduction path in the horizontal direction. The radial drift of the furnace structure is mapped to a 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 structure component are formed, reflecting the structural changes of the heat transfer path caused by the structural drift of the furnace structure component in various directions and positions.

[0090] According to the structural change characteristics of the heat transfer path, the change law of the heat conduction interface contact area of ​​the furnace structure components is analyzed to determine the change amount of the heat conduction interface contact area;

[0091] The interface reconstruction analysis method is adopted, and the structural change characteristics are used as input parameters to calculate the influence of the relative offset of each interface position in the structural change characteristics on the contact area of ​​the heat conduction interface.

[0092] Specifically, the interface reconstruction analysis method includes the following: determining the position and area of ​​the initial interface contact area; determining the specific offset value at each specific interface contact position based on the offset value of each position in the structural change characteristics; and calculating the interface overlap area by analyzing the degree of overlap between the offset value of the interface contact position and the initial interface position. The calculation method is to calculate the intersection area of ​​the initial interface contact area and the actual contact area after the interface offset, i.e., the interface overlap area; and calculating the change in the interface overlap area during each thermal cycle to obtain the pattern of interface contact area change during the thermal cycle.

[0093] The interface reconstruction analysis method described above generates a curve for how the interface contact area changes with thermal cycling conditions. This curve depicts the trend and specific rate of change in the contact area of ​​the heat conduction interface of the furnace structure components as it decreases or increases with the number of thermal cycles, ultimately determining the change in the heat conduction interface contact area.

[0094] According to the change in the contact area of ​​the heat conduction interface, the change in the heat transfer capacity of the thermal contact interface of the furnace structure component is analyzed to obtain the contact thermal resistance change information of the thermal contact interface of the furnace structure component;

[0095] Specifically, the theory of interfacial heat conduction is used to develop a formula for calculating the thermal conductivity of a thermal contact interface: the thermal conductivity of a thermal contact interface is equal to the product of the material's thermal conductivity and the interface's contact area, divided by the thermal resistance thickness of the interface. This formula shows that changes in the interface's contact area directly affect the thermal conductivity. The change in the interface's contact area is input into the formula to obtain the thermal conductivity value for each change in contact area. Finally, the difference in thermal conductivity relative to the initial state during each thermal cycle is calculated to provide information on the change in thermal conductivity of the interface.

[0096] Based on the change in heat transfer capacity of the thermal contact interface of the furnace structural components, the change in contact thermal resistance of the thermal contact interface of the furnace structural components is obtained. Using the thermal resistance calculation method, the contact thermal resistance is calculated as follows: 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. Based on the change in thermal contact conductivity of the thermal contact interface and the actual measured temperature difference data of the thermal contact interface, the change in thermal contact resistance of the thermal contact interface over the thermal cycle is calculated. The contact thermal resistance increases as the thermal conductivity of the thermal contact interface decreases.

[0097] S4. Based on the geometric parameter drift vector, the furnace wall heat flow field is reconstructed to obtain the dynamic correction vector of the actual heat flow distribution, including:

[0098] Based on the geometric shape drift vectors generated by the furnace structure components under thermal cycling conditions, a numerical reconstruction model of the furnace wall thermal flow field is established;

[0099] The geometric drift vector represents the drift amplitude and position of the furnace wall structural components in the horizontal, vertical, and radial directions. The numerical reconstruction model of the furnace wall thermal flow field is established using a numerical simulation method for thermal flow field reconstruction. The finite difference method is selected as the numerical simulation method. This method divides the furnace wall surface into multiple continuous grid cells. By discretizing the furnace wall surface, a finite difference calculation grid is formed. The size and distribution density of the grid cells are specified to ensure that the numerical simulation accuracy meets the requirements of practical engineering applications.

[0100] Based on the initial geometric structure of the furnace wall, an initial numerical model of the furnace wall thermal flow field is established. The numerical model of the furnace wall thermal flow field defines the initial thermal flow boundary conditions of the furnace wall surface, that is, the heat flux density distribution of each grid unit on the furnace wall surface in the undeformed state. The boundary position of the initial numerical model of the furnace wall thermal flow field is corrected using the geometric shape drift vector generated by the furnace body structural components under thermal cycle conditions, that is, the position of each grid unit in the horizontal, vertical and radial directions is adjusted to form a corrected grid unit position of the thermal flow field. The correction method is: the initial position of each grid unit of the furnace wall is added with the drift vector amplitude of the corresponding position, to realize the mapping conversion from the initial grid unit to the grid unit position under the actual drift state, and to establish a numerical reconstruction model of the furnace wall thermal flow field.

[0101] Based on the numerical reconstruction model of the furnace wall thermal flow field, the spatial distribution characteristics of the heat flux on the furnace wall surface are reconstructed and analyzed to obtain the distribution variation law of the heat flux on the furnace wall surface under the actual operating conditions of the furnace wall structural components;

[0102] The spatial distribution characteristics of the heat flux on the furnace wall surface are used as input data, including the average value, peak value, and degree of heat flux unevenness at different spatial locations on the furnace wall surface. The above input data are mapped to the positions of each corrected grid unit in the numerical reconstruction model of the furnace wall thermal flow field to achieve 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 value of the furnace wall surface is assigned to each grid unit according to the actual position of the furnace wall grid unit; then, the heat conduction equation is numerically solved to calculate the heat flux conduction relationship between each grid unit, thereby obtaining the actual distribution of the heat flux on the furnace wall surface at each corrected grid unit position.

[0103] In order to obtain the heat flux distribution variation pattern of the furnace wall structure components under the actual operating conditions, a heat flux variation trend analysis method is adopted. The implementation process is as follows: the difference 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 is calculated; the spatial distribution of the difference value on the entire furnace wall surface is statistically analyzed; and based on the spatial distribution of the difference value, the overall heat flux distribution variation trend of the furnace wall structure components is analyzed. This is manifested in the regional distribution of the increase or decrease in the heat flux of each grid cell, the trend of the change in the heat flux distribution unevenness, and the location of the significant heat flux change in the local area. Finally, the heat flux distribution variation pattern of the furnace wall structure components under the actual operating conditions is obtained.

[0104] According to the variation law of heat flux distribution on the furnace wall surface, the position of the local heat flux offset area and the heat flux variation amplitude of the furnace wall structural component are determined, and the dynamic correction vector of the actual heat flux distribution of the furnace wall structural component is generated;

[0105] Within the heat flux distribution variation pattern of the furnace wall surface, local areas with significant heat flux changes are identified. These areas are defined as locations 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 in the overall heat flux of the furnace wall. The grid cells exceeding the preset threshold are located to identify the locations of local heat flux offsets of the furnace wall structural components. The magnitude of the heat flux change in these local offset areas is calculated as the difference between the actual heat flux value of the local offset area and the heat flux value of the corresponding location in the initial state. The positive or negative value of this difference indicates an increasing or decreasing trend in the heat flux.

[0106] The positions of the localized heat flux offsets and the corresponding heat flux variation amplitudes of the furnace wall structural components are represented by vectors. This is done by using the coordinates of each localized offset region on the furnace wall surface as the base point and the heat flux variation amplitude as the vector length to form a dynamic heat flux correction vector. All these dynamic heat flux correction vectors are then arranged according to the position of the furnace wall geometry to form a dynamic correction vector for the actual heat flux distribution of the furnace wall structural components. This dynamic correction vector expresses the characteristics of the actual heat flux distribution of the furnace wall structural components as it changes with structural drift under thermal cycling conditions.

[0107] 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, build a time-varying boundary condition inversion model, and output a multi-parameter coupled boundary correction set, including:

[0108] Based on the contact thermal resistance variation information of the thermal contact interface of the furnace structure components, a time-varying model of the contact thermal resistance parameter variation is established;

[0109] Contact thermal resistance variation characterizes the thermal resistance variation at the thermal contact interface of furnace components under long-term thermal cycling conditions, resulting from structural geometric drift and changes in interface contact area. This variation is expressed as a trend curve showing the cumulative change in thermal resistance over the thermal cycling period. To reflect the dynamic process of contact thermal resistance variation over time, a time series analysis method was used to construct a time-varying model of contact thermal resistance parameter changes.

[0110] The time series analysis method uses an autoregressive moving average model: the contact thermal resistance values ​​within each thermal cycle are used as the basic data sequence. The basic data sequence is preprocessed, including detrending and standardization. The order of the autoregressive moving average model is selected, and the order selection is determined using the Bayesian Information Criterion. After the order is determined, the model parameters of the basic data sequence are estimated using the maximum likelihood estimation method to obtain a time-varying model of the contact thermal resistance parameters as they change over the thermal cycle. The time-varying model is expressed as the value of the contact thermal resistance in the current cycle being the sum of the products of the contact thermal resistance values ​​in several historical cycles and the model estimated parameters, plus the model error term. The model error term is a random disturbance term, representing the impact of unpredictable small disturbances on the thermal resistance change.

[0111] Based on the dynamic correction vector of the actual heat flux distribution of the furnace wall structural components, a time-varying model of the heat flux dynamic correction parameters is established;

[0112] The dynamic correction vector for the actual heat flux distribution represents the heat flux variation caused by geometric drift of the furnace wall components under thermal cycling conditions. It includes the location of each localized offset region on the furnace wall surface and the corresponding heat flux variation. To reflect the dynamic changes in actual heat flux over the thermal cycling period, a state-space model was used to establish a time-varying model for the dynamic heat flux correction parameters.

[0113] The process of establishing the state space model is as follows: constructing a state vector, which contains the heat flux correction amplitude of each local area of ​​the furnace wall; defining the state transfer equation, which is expressed as the current heat flux correction state vector is equal to the heat flux correction state vector of the previous cycle multiplied by the state transfer matrix plus the process disturbance vector. The state transfer 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 obeying the Gaussian distribution; constructing the observation equation, which is expressed as the actual observed heat flux distribution correction data vector is 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 to estimate the state space model parameters to obtain a time-varying model of the dynamic correction parameters of the heat flux of the furnace wall structural components.

[0114] Based on the time-varying model of the contact thermal resistance parameters and the time-varying model of the heat flow dynamic correction parameters, a multi-parameter comprehensive analysis of the heat transfer boundary conditions of the furnace structure components is conducted, and an inversion model of the boundary conditions of the dynamic changes in the heat transfer characteristics of the furnace structure components is constructed;

[0115] A multi-parameter collaborative analysis method is adopted, specifically: the time-varying model output value of the contact thermal resistance parameter is synchronized with the time-varying model output value of the heat flow dynamic correction parameter; the synchronization method is to determine a unified thermal cycle time scale, and the time-varying model output value of the contact thermal resistance parameter is aligned with the time-varying model output value of the heat flow dynamic correction parameter according to a unified period to form two parameter sequences; then, the parameter sequence is used as input to establish a parameter collaborative analysis matrix, and each element of the collaborative analysis matrix represents the mutual influence intensity between the contact thermal resistance parameter and the heat flow correction parameter in the corresponding thermal cycle; the calculation method is to normalize the product of the two parameter changes in the same period; finally, by analyzing the eigenvalues ​​of the collaborative analysis matrix, the dynamically changing characteristic parameter set of the heat transfer boundary conditions of the furnace structure components is obtained.

[0116] Based on the dynamically changing characteristic parameter set, an inverse boundary condition model for the dynamic changes in the heat transfer characteristics of the furnace body structural components is constructed. The boundary condition inverse model adopts an inverse problem solving method: a heat transfer forward model of the furnace wall structural components is set, and the heat transfer forward model expresses the calculation formula for the temperature field distribution of the furnace wall structural components under specific heat transfer boundary condition parameters; using the dynamically changing characteristic parameter set as input, an objective function is constructed, and the objective function is expressed as the sum of the squares of the difference between the actual measured temperature field and the temperature field calculated by the forward model; a genetic algorithm optimization method is used to solve the minimum value of the objective function to obtain the optimal estimation result of the heat transfer boundary condition parameters; the initial population, crossover rate, mutation rate, and selection rules are defined in the genetic algorithm, and the fitness function is used to evaluate the optimization degree of each generation of population to obtain the optimal heat transfer boundary condition parameter estimation value under actual conditions.

[0117] Through the boundary condition inversion model, the multi-parameter coupled boundary correction set of heat transfer boundary conditions under the thermal cycle conditions of the furnace structure components is output;

[0118] The multi-parameter coupled boundary correction set includes contact thermal resistance correction parameters and furnace wall heat flux distribution correction parameters, which are expressed in the form of multi-dimensional vectors.

[0119] S6. Update the boundary conditions in the water-jacketed heating furnace heat transfer simulation model based on the multi-parameter coupled boundary correction set, including:

[0120] Based on the multi-parameter coupled boundary correction set, the correction parameters of contact thermal resistance and furnace wall heat flux distribution are determined;

[0121] The contact thermal resistance correction parameters and furnace wall heat flux distribution correction parameters are extracted from the multi-parameter coupled boundary correction set.

[0122] The contact thermal resistance correction parameter is expressed as the corrected value of the contact thermal resistance parameter at each location on the thermal contact interface, and the heat flux distribution correction parameter is expressed as the corrected value of the heat flux at each local location on the furnace wall. To ensure that the contact thermal resistance correction parameter and the furnace wall heat flux distribution correction parameter are applicable to the water-jacketed furnace heat transfer simulation model, a parameter mapping process is implemented. The parameter mapping method is as follows: Based on the initial design structural position of each contact interface location and the furnace wall area in the water-jacketed furnace heat transfer simulation model, a mapping relationship is established between the contact thermal resistance correction parameter and the furnace wall heat flux distribution correction parameter and the water-jacketed furnace heat transfer simulation model position. In other words, each correction parameter is mapped to the corresponding geometric position in the water-jacketed furnace heat transfer simulation model.

[0123] Update the contact thermal resistance correction parameters to the thermal contact interface contact thermal resistance boundary conditions in the water jacket heating furnace heat transfer simulation model;

[0124] A heat transfer simulation model for a water-jacketed furnace is established. This model is based on the initial geometry, material properties, and heat transfer boundary conditions from the water-jacketed furnace design phase. This model defines the initial contact resistance of the thermal contact interface. Specifically, this initial contact resistance is the thermal resistance value set during the design phase based on ideal contact conditions.

[0125] The method for updating the contact thermal resistance correction parameters is: the correction value corresponding to each position in the contact thermal resistance correction parameters is replaced one by one with the initial contact thermal resistance parameters of the thermal contact interface of the water jacket heating furnace heat transfer simulation model: the contact thermal resistance value of each position in the thermal contact interface boundary conditions of the water jacket heating furnace simulation model is replaced with the contact thermal resistance correction parameters of the corresponding position to ensure that the thermal contact interface boundary conditions in the water jacket heating furnace simulation model accurately express the dynamic characteristics of the contact thermal resistance under actual operating conditions.

[0126] The furnace wall heat flux distribution correction parameters are updated to the furnace wall surface heat flux spatial distribution boundary conditions in the water jacket heating furnace heat transfer simulation model;

[0127] The initial boundary condition of the heat flux spatial distribution on the furnace wall surface of the water-jacketed heating furnace heat transfer simulation model is defined as the initial heat flux density distribution state set according to the ideal operating conditions in the design stage, which is expressed as the initial heat flux density value of each local area on the furnace wall surface.

[0128] The update method involves performing spatial meshing of the heat flux boundary conditions on the furnace wall surface of the simulation model, dividing the furnace wall surface into several localized grid regions. Each localized grid region is assigned an initial heat flux density value in the initial boundary conditions, and the heat flux distribution correction parameters are mapped to the grid regions on the furnace wall surface. The position mapping method determines the grid positions corresponding to each heat flux correction value in the heat flux distribution correction parameter based on the initial position coordinates of each localized grid region on the furnace wall surface. The initial heat flux density value for each grid region on the furnace wall surface is replaced with the correction value for the corresponding position in the heat flux distribution correction parameter under actual operating conditions. The update process ensures that the spatial heat flux distribution boundary conditions on the furnace wall surface reflect the dynamic changes in the actual heat flux.

[0129] The updated thermal contact interface thermal resistance boundary conditions and the furnace wall surface heat flux spatial distribution boundary conditions form the current heat transfer boundary conditions of the water jacket heating furnace heat transfer simulation model;

[0130] The current heat transfer boundary conditions include the updated thermal contact interface boundary conditions and the furnace wall surface heat flux boundary conditions. In order to verify the effectiveness and accuracy of the updated water jacket heating furnace heat transfer simulation model, the following verifications are required:

[0131] The furnace wall surface temperature distribution and thermal stress response data measured under actual operating conditions are used as comparison benchmarks for the output results of the water-jacketed heating furnace heat transfer simulation model. The simulation process uses 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 used to simulate and solve 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 operating conditions. The simulation results of the furnace wall surface temperature field and thermal stress field are obtained by calculation.

[0132] The simulation results were compared and analyzed with the actual measured data, and the error between the two was calculated. At each measurement location on the furnace wall, the difference between the simulated temperature and thermal stress values ​​and the actual measured temperature and thermal stress values ​​was calculated. All differences were statistically analyzed to obtain the distribution of the errors. Based on the results of the error statistical analysis, the accuracy of the updated heat transfer simulation model boundary conditions was evaluated. If the error met the preset error threshold, it indicated that the updated heat transfer boundary conditions accurately described the furnace wall heat transfer state under actual operating conditions. The error threshold was derived from the allowable deviation standard for actual engineering applications of water-jacketed heating furnaces. The error threshold was determined based on the maximum temperature and thermal stress deviations that the furnace wall structural components could tolerate under safe operating conditions, and was determined based on the operating specifications, design specifications, performance index requirements, and long-term operating experience data of the water-jacketed heating furnace.

[0133] Example 2

[0134] The difference between Example 2 of the present invention and Example 1 is that this example introduces a self-correction optimization system for heat transfer boundary conditions of a water-jacket heating furnace.

[0135] Figure 2 A schematic structural diagram of a water-jacket heating furnace heat transfer boundary condition self-correction optimization system is provided. The water-jacket heating furnace heat transfer boundary condition self-correction optimization system comprises:

[0136] The state acquisition module collects 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 cycle conditions to obtain the comprehensive state parameters of the furnace body structure components under thermal cycle conditions;

[0137] The trend estimation module builds a thermal fatigue deformation trend estimation model based on the comprehensive state parameters of the furnace structure components under thermal cycle conditions, and generates a geometric parameter drift vector for evaluating the deformation of the furnace structure components;

[0138] Thermal resistance analysis module, performs structural equivalent mapping analysis on geometric parameter drift vectors and extracts contact thermal resistance change information in the heat transfer path;

[0139] The heat flow reconstruction module reconstructs the furnace wall heat flow field based on the geometric parameter drift vector and obtains the dynamic correction vector of the actual heat flow distribution;

[0140] The boundary inversion module 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 coupled boundary correction set;

[0141] The model update module updates the boundary conditions in the water jacket heating furnace heat transfer simulation model according to the multi-parameter coupling boundary correction set.

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

[0143] The above embodiments can be implemented in whole or in part via 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 comprises one or more computer instructions or computer programs. When loaded or executed on a computer, the processes or functions described in the embodiments of this application are fully or partially performed. The computer can be a general-purpose computer, a special-purpose computer, a computer network, or other programmable device. The computer instructions can be stored in a computer-readable storage medium or transmitted from one computer-readable storage medium to another. For example, the computer instructions can be transmitted from one website, computer, server, or data center to another website, computer, server, or data center via wired means (e.g., infrared, wireless, microwave, etc.). The computer-readable storage medium can be any available medium accessible by a computer or a data storage device such as a server or data center that contains a collection of one or more available media. The available medium can be magnetic media (e.g., floppy disks, hard disks, tapes), optical media (e.g., DVDs), or semiconductor media. The semiconductor media can be a solid-state drive.

[0144] Those skilled in the art will appreciate that the modules and algorithm steps of each example described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, or a combination of computer software and electronic hardware. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Professional and technical personnel can use different methods to implement the described functions for each specific application, but such implementation should not be considered beyond the scope of this application.

[0145] Those skilled in the art will 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 aforementioned method embodiments and will not be repeated here.

[0146] In the several embodiments provided in this 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 schematic. For example, the division of the modules is only a logical function division. In actual implementation, there may be other division methods, such as 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 mutual coupling or direct coupling or communication connection shown or discussed can be through some interfaces, indirect coupling or communication connection of devices or modules, which can be electrical, mechanical or other forms.

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

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

[0149] If the 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 this understanding, the technical solution of the present application, or the part that contributes to the prior art, or the part of the technical solution, can be embodied in the form of a software product. The computer software product is stored in a storage medium and includes several instructions for enabling a computer device (which can be a personal computer, server, or network device, etc.) to execute all or part of the steps of the method described in each embodiment of the present application. The aforementioned storage medium includes various media that can store program codes, such as a USB flash drive, a mobile hard disk, a read-only memory (ROM), a random access memory (RAM), a magnetic disk, or an optical disk.

[0150] The above description is merely a specific embodiment of the present application, but the scope of protection of the present application is not limited thereto. Any changes or substitutions that can be easily conceived by a person skilled in the art within the technical scope disclosed in this application should be included in the scope of protection of this application. Therefore, the scope of protection of this application should be based on the scope of protection of the claims.

[0151] Finally: The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present invention should be included in the scope of protection of the present invention.

Claims

1. A self-correction optimization method for heat transfer boundary conditions of a water-jacket heating furnace, characterized in that: The steps include: S1. Collecting 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 cycle conditions to obtain the comprehensive state parameters of the furnace body structure components under thermal cycle conditions; S2. Based on the comprehensive state parameters of the furnace structure components under thermal cycle conditions, a thermal fatigue deformation trend estimation model is constructed, and a geometric parameter drift vector is generated for evaluating the deformation of the furnace structure components; S3. Perform structural equivalent mapping analysis on the geometric parameter drift vector to extract the contact thermal resistance change information in the heat transfer path; S4. Reconstruct the furnace wall heat flow field based on the geometric parameter drift vector to obtain the dynamic correction vector of the actual heat flow 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, build a time-varying boundary condition inversion model, and output a multi-parameter coupled boundary correction set, specifically: Based on the contact thermal resistance variation information of the thermal contact interface of the furnace structure components, a time-varying model of the contact thermal resistance parameter variation is established; Based on the dynamic correction vector of the actual heat flux distribution of the furnace wall structural components, a time-varying model of the heat flux dynamic correction parameters is established; Based on the time-varying model of the contact thermal resistance parameters and the time-varying model of the heat flow dynamic correction parameters, a multi-parameter comprehensive analysis of the heat transfer boundary conditions of the furnace structure components is conducted, and an inversion model of the boundary conditions of the dynamic changes in the heat transfer characteristics of the furnace structure components is constructed; Through the boundary condition inversion model, the multi-parameter coupled boundary correction set of heat transfer boundary conditions under the thermal cycle conditions of the furnace structure components is output; S6. Update the boundary conditions in the water jacket heating furnace heat transfer simulation model according to the multi-parameter coupled boundary correction set.

2. The method for self-correction optimization of heat transfer boundary conditions of a water-jacket heating furnace according to claim 1, characterized in that: S1, 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 cycle conditions during the operation of the water jacket heating furnace; The temperature variation characteristics of the furnace wall surface are obtained by analyzing the time series variation of the furnace wall surface temperature. Perform spatial distribution analysis on the heat flux distribution information on the furnace wall surface to obtain the heat flux spatial distribution characteristics on the furnace wall surface; Perform feature recognition on the furnace wall thermal stress response curve to obtain the response characteristics of the furnace wall thermal stress change; 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 changes of the furnace wall are combined to obtain the comprehensive state parameters of the furnace structure components under thermal cycle conditions.

3. The method for self-correction optimization of heat transfer boundary conditions of a water jacket heating furnace according to claim 2, characterized in that: S2, specifically: Based on the comprehensive state parameters of the furnace structure components under thermal cycle conditions, the correlation between thermal expansion displacement and structural deformation is established; Determine the thermal fatigue deformation trend of the furnace structure components based on the correlation between thermal expansion displacement and structural deformation; According to the thermal fatigue deformation trend of the furnace structure components, calculate the geometric shape drift of the furnace structure components under thermal cycle conditions; Based on the geometric shape drift of the furnace structure components under thermal cycle conditions, a geometric parameter drift vector is generated for evaluating the deformation of the furnace structure components.

4. The method for self-correction optimization of heat transfer boundary conditions of a water-jacket heating furnace according to claim 3, characterized in that: S3, specifically: The structural equivalent mapping method is used to convert the geometric shape drift vectors of the furnace 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, the change law of the heat conduction interface contact area of ​​the furnace structure components is analyzed to determine the change amount of the heat conduction interface contact area; According to the change in the contact area of ​​the heat conduction interface, the change in the heat transfer capacity of the thermal contact interface of the furnace structure component is analyzed, and the contact thermal resistance change information of the thermal contact interface of the furnace structure component is obtained.

5. The method for self-correction optimization of heat transfer boundary conditions of a water-jacket heating furnace according to claim 4, characterized in that: S4, specifically: Based on the geometric shape drift vectors generated by the furnace structure components under thermal cycling conditions, a numerical reconstruction model of the furnace wall thermal flow field is established; Based on the numerical reconstruction model of the furnace wall thermal flow field, the spatial distribution characteristics of the heat flux on the furnace wall surface are reconstructed and analyzed to obtain the distribution variation law of the heat flux on the furnace wall surface under the actual operating conditions of the furnace wall structural components; According to the variation law of heat flux distribution on the furnace wall surface, the position of the local heat flux offset area and the heat flux variation amplitude of the furnace wall structural component are determined, and the dynamic correction vector of the actual heat flux distribution of the furnace wall structural component is generated.

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

7. A water-jacket heating furnace heat transfer boundary condition self-correction optimization system, used to implement a water-jacket heating furnace heat transfer boundary condition self-correction optimization method according to any one of claims 1 to 6, characterized in that: include: The state acquisition module collects 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 cycle conditions to obtain the comprehensive state parameters of the furnace body structure components under thermal cycle conditions; The trend estimation module builds a thermal fatigue deformation trend estimation model based on the comprehensive state parameters of the furnace structure components under thermal cycle conditions, and generates a geometric parameter drift vector for evaluating the deformation of the furnace structure components; Thermal resistance analysis module, performs structural equivalent mapping analysis on geometric parameter drift vectors and extracts contact thermal resistance change information in the heat transfer path; The heat flow reconstruction module reconstructs the furnace wall heat flow field based on the geometric parameter drift vector and obtains the dynamic correction vector of the actual heat flow distribution; The boundary inversion module 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 coupled boundary correction set, specifically: Based on the contact thermal resistance variation information of the thermal contact interface of the furnace structure components, a time-varying model of the contact thermal resistance parameter variation is established; Based on the dynamic correction vector of the actual heat flux distribution of the furnace wall structural components, a time-varying model of the heat flux dynamic correction parameters is established; Based on the time-varying model of the contact thermal resistance parameters and the time-varying model of the heat flow dynamic correction parameters, a multi-parameter comprehensive analysis of the heat transfer boundary conditions of the furnace structure components is conducted, and an inversion model of the boundary conditions of the dynamic changes in the heat transfer characteristics of the furnace structure components is constructed; Through the boundary condition inversion model, the multi-parameter coupled boundary correction set of heat transfer boundary conditions under the thermal cycle conditions of the furnace structure components is output; The model update module updates the boundary conditions in the water jacket heating furnace heat transfer simulation model according to the multi-parameter coupling boundary correction set.

Images