Server heat dissipation regulation and control method and device and medium

Abstract

The invention relates to the technical field of heat dissipation regulation and control, in particular to a server heat dissipation regulation and control method and device and a medium. The method comprises the following steps: acquiring a core temperature, an accumulative error and a task vector containing the characteristics of instruction composition, floating point density and the like; predicting an expected heat flux based on the model, calculating a thermal penalty factor positively correlated with an accumulative error, and generating a scheduling scheme by taking the heat flux accounting for penalty in a balanced manner and total execution time as targets; an integral step length inversely proportional to the temperature error change rate is set to calculate an accumulative error; predicting a temperature track through a first-order heat conduction model by using heat flux, and taking an initial slope as a differential item; and calculating a proportional gain associated with the expected heat flux and the accumulative error, and summing each control quantity to generate a heat dissipation control instruction. According to the scheme, local hot spots can be avoided, the control predictability is enhanced, and stable and coordinated temperature regulation is realized.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of heat dissipation regulation. More specifically, the present application relates to a server heat dissipation regulation method, device and medium. BACKGROUND

[0002] With the rapid development of server hardware technology, the integration of chips is constantly improving, which makes the power density of CPU, GPU and other core computing units rise sharply. This high power density operating state causes the server to generate a lot of heat when processing high-load tasks, which will seriously affect the stability and service life of the equipment if not cooled in time and effectively. In order to maintain the normal operation of the system, the mainstream server cooling scheme usually adopts a feedback control mechanism based on temperature sensors. Specifically, this mechanism monitors the temperature of the core area in real time, and when the monitored temperature value exceeds the preset safety threshold, the system will control the cooling fan to increase the speed, aiming to enhance air convection to achieve cooling effect, thereby preventing hardware overheating.

[0003] In order to further optimize the above-mentioned cooling control effect, the proportional-integral-derivative (PID) control algorithm is often introduced into the cooling system in the prior art. The PID control algorithm not only depends on the current temperature error, but also combines the historical cumulative error and the prediction of the future error change trend for comprehensive operation. Compared with simple threshold switch control, the PID algorithm can output more accurate control signals, thereby providing smoother and timely response. This improved control method effectively suppresses the sharp fluctuations of temperature to a certain extent, improves the dynamic adjustment ability of the cooling system to temperature changes, and makes the server maintain a relatively stable operating temperature within a certain range.

[0004] However, the above prior art still has significant defects in practical application. First, the control logic of traditional threshold control and PID control is based on the temperature measurement value after the event, which has obvious hysteresis and cannot make a prediction and intervention before the temperature rises. Second, at the task scheduling level, the existing scheduler usually only takes the computational performance indicators such as throughput and delay as the optimization target, ignoring the heat dissipation factor. This single-dimensional scheduling strategy may inadvertently allocate multiple high-heat-power tasks to a few cores, causing these cores to form local hot spots that are difficult to eliminate through regular fan adjustment in a short period of time. Third, the derivative term of the PID algorithm is calculated based on historical data, which is highly sensitive to delay and easily disturbed by signal noise; and the integral term is prone to integral saturation phenomenon in the case of long-time error accumulation, which will cause the system to overshoot and oscillate, thereby destroying the stability of the system. In summary, the existing technology lacks a perspective of deeply integrating task scheduling and heat dissipation control, and lacks an advanced strategy that can adapt to dynamic changes in load and accurately predict future heat trends, making it difficult to achieve low-energy and high-stability collaborative regulation of server temperature. SUMMARY

[0005] The purpose of the present application is to provide a server heat dissipation regulation method, device and medium to solve the problem of low-energy and high-stability collaborative regulation of server temperature in the prior art. To this end, the present application provides solutions in the following four aspects.

[0006] In the first aspect, the present application provides a server heat dissipation regulation method, comprising the following steps: obtaining the current temperature of each computing core of the server, the cumulative temperature error integral value of the last control period, and the task feature vector of the task set to be scheduled, wherein the task feature vector includes the instruction type composition ratio, the floating point operation density and the memory access frequency of the task; based on the task feature vector and the preset instruction set heat power consumption model, predicting the expected heat flux time series of the task set to be scheduled on each core; calculating a thermal penalty factor positively correlated with the cumulative temperature error integral value of each computing core in the last control period; generating and executing a task scheduling scheme to balance the expected heat flux of each core with the thermal penalty factor and the total execution time as the target; calculating the temperature error of the current temperature and the target temperature; setting an integral time step inversely proportional to the absolute value of the temperature error rate, and calculating the cumulative temperature error integral value of the current period by using the trapezoidal integral method; predicting the temperature change trajectory in the future time window by using the expected heat flux time series and the current temperature through the first-order heat conduction model of the server, and taking the initial slope of the trajectory as the derivative term; Calculate a proportional gain coefficient that is positively correlated with the absolute value of the expected heat flux and the cumulative temperature error integral value of the current cycle; based on the proportional gain coefficient, the cumulative temperature error integral value and the differential term, generate proportional, integral and differential control quantities and sum them to generate a heat dissipation control command for adjusting the speed of the server cooling fan.

[0007] Preferably, the step of predicting the expected heat flux time series of the scheduled task set executed on each core includes: decomposing each task in the task set into a sequence of micro-operations within a unit time slice; querying a preset instruction set heat power model to obtain the basic heat power value corresponding to each micro-operation; weighting and summing the basic heat power values ​​of each micro-operation according to the instruction type composition ratio, floating-point operation density, and memory access frequency in the task feature vector to obtain the instantaneous heat power within the time slice; and concatenating the instantaneous heat power values ​​in consecutive time slices to obtain the expected heat flux time series executed on the corresponding core.

[0008] Preferably, the formula for calculating the thermal penalty factor is: ,in, For the first One core thermal penalty factor, For the first The cumulative temperature error integral value of the previous cycle of each core. This is the preset penalty gain coefficient.

[0009] Preferably, the step of generating and executing a task scheduling scheme with the goal of balancing the expected heat flux of each core component and the total execution time, which are included in the thermal penalty factor, includes: Construct the optimization objective function ,in, The total completion time for the task set. For the first One core thermal penalty factor, To be allocated to the first The expected total heat production of all tasks in each core, and For preset weighting coefficients, and These are the normalized time and normalized heat used to make the components of the objective function dimensionless; the particle swarm optimization algorithm is used to iteratively optimize the allocation scheme of tasks on each core until the objective function is achieved. The value of converges to the minimum, thus obtaining the task scheduling scheme.

[0010] Preferably, the step of calculating the cumulative temperature error integral value for the current period includes: Calculate the current temperature error Temperature error compared to the previous cycle the difference between the temperature error and the temperature error of the previous period ; According to the formula Calculate the integral time step , wherein is the basic time step, is the preset integral step adjustment coefficient; The cumulative temperature error integral value of the current period is calculated by the formula , wherein is the cumulative temperature error integral value of the previous period.

[0011] Preferably, the temperature change trajectory in the future time window is predicted by the server first-order heat conduction model using the expected heat flux time series and the current temperature, and the initial slope of the trajectory is taken as the differential term, comprising: Establish the first-order heat conduction differential equation of the server core: , wherein is the core temperature, is the expected heat flux, is the ambient temperature, and is the pre-calibrated model coefficient; Substitute the temperature at the current time and the initial value of the expected heat flux time series into the equation to calculate the initial slope ; The calculated initial slope is taken as the differential term.

[0012] Preferably, the calculation formula of the proportional gain coefficient is: , wherein is the proportional gain coefficient, is the basic proportional gain coefficient, is the average value of the expected heat flux in the future time window, is the absolute value of the cumulative temperature error integral value of the current period, and is the preset gain adjustment coefficient.

[0013] In the second aspect, a server heat dissipation control device comprises the following modules: ​The acquisition module is configured to acquire the current temperature of each computing core of the server, the accumulated temperature error integral value of the last control period, and a task feature vector of a task set to be scheduled, the task feature vector including a task instruction type composition ratio, a floating point operation density, and a memory access frequency; the generation module is configured to predict an expected heat flux time sequence of the task set to be scheduled executing on each core based on the task feature vector and a preset instruction set thermal power consumption model; calculate a thermal penalty factor positively correlated with the accumulated temperature error integral value of each computing core in the last control period; generate and execute a task scheduling scheme with the goal of balancing the expected heat flux of each core with the total execution time of the thermal penalty factor; the prediction module is configured to calculate a temperature error of the current temperature and a target temperature; set an integral time step inversely proportional to the absolute value of the temperature error rate, and calculate the accumulated temperature error integral value of the current period using the trapezoidal integration method; predict a temperature change trajectory in a future time window using the first-order heat conduction model of the server using the expected heat flux time sequence and the current temperature, and use the initial slope of the trajectory as a differential term; the calculation module is configured to calculate a proportional gain coefficient positively correlated with the absolute value of the expected heat flux and the accumulated temperature error integral value of the current period; generate a proportional, integral, and differential control amount based on the proportional gain coefficient, the accumulated temperature error integral value, and the differential term, and sum them up to generate a heat dissipation control instruction for regulating the rotation speed of the server heat dissipation fan.

[0014] In a third aspect, a computer-readable storage medium comprises: The computer program stored on the readable storage medium is executed by a processor to implement the server heat dissipation regulation method described above.

[0015] The beneficial effects of the present application are: the present application inputs the prediction of the thermal characteristics of the task in the task scheduling stage, and combines the historical thermal state of the core to generate a scheduling scheme that can balance the expected thermal load of each computing core, thereby avoiding the local hot spot problem caused by uneven task allocation from the root, and improving the uniformity of the overall temperature field. The differential term is constructed by predicting the future temperature change trend based on the heat conduction model, replacing the traditional differential link, enhancing the predictability and anti-interference ability of the control; by inputting the proportional gain associated with the expected heat flux and using the trapezoidal integration method to calculate the accumulated error, the accuracy of the control response and the stability of the system are improved, and the temperature overshoot and oscillation are suppressed. The closed-loop feedback control of the task-level thermal behavior prediction is combined to achieve more stable cooperative regulation of the server temperature. BRIEF DESCRIPTION OF DRAWINGS

[0016] Figure 1 A step flowchart of the server heat dissipation regulation method in the embodiment is schematically shown. DETAILED DESCRIPTION

[0017] The technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention.

[0018] like Figure 1 As shown, a server heat dissipation control method in this embodiment includes the following steps: Step S1: Obtain the current temperature of each computing core of the server, the cumulative temperature error integral value of the previous control cycle, and the task feature vector of the task set to be scheduled. The task feature vector includes the instruction type composition ratio, floating-point operation density, and memory access frequency of the task.

[0019] Real-time temperature data is acquired by reading the digital thermal sensors (DTS) built into each computing core; the cumulative temperature error integral value stored at the end of the previous control cycle is read from memory; for each task in the set of tasks to be scheduled, offline or online analysis is performed through static code analysis or hardware performance counter (PMU) to count the number of different types of instructions such as integer operations, floating-point operations, memory accesses, and branches during execution, calculate the proportion of each type of instruction to the total number of instructions, and count the number of floating-point operations and memory accesses per unit time to obtain the task feature vector.

[0020] In an optional embodiment, the step of predicting the expected heat flux time series of the set of tasks to be scheduled executed on each core includes: Each task in the task set is decomposed into a sequence of micro-operations within a unit time slice; a preset instruction set thermal power consumption model is queried to obtain the basic thermal power consumption value corresponding to each micro-operation; based on the proportion of instruction types, floating-point operation density, and memory access frequency in the task feature vector, the basic thermal power consumption value of each micro-operation is weighted and summed to obtain the instantaneous thermal power consumption within the time slice; the instantaneous thermal power consumption values ​​in consecutive time slices are concatenated to obtain the expected heat flux time sequence executed on the corresponding core.

[0021] Specifically, this prediction process divides a scheduled computational task into multiple unit time slices. For example, a video transcoding task is divided into multiple 10ms unit time slices based on its execution time. Within the first 10ms time slice, the execution of the task is decomposed into a series of micro-operations, such as integer addition, floating-point multiplication, and memory reads. A pre-established power consumption model database is queried, which records the basic thermal power consumption of different micro-operations, for example, integer addition is 0.1W, floating-point multiplication is 0.3W, and memory reads are 0.2W.

[0022] The feature vector of the video transcoding task is extracted, showing a high density of floating-point operations and a moderate frequency of memory accesses. Based on these features, the base thermal power consumption is adjusted; for example, the power consumption of floating-point multiplication is multiplied by a weighting factor of 1.2, resulting in 0.36W. The weighted thermal power consumption values ​​of all micro-operations within this time slice are summed to obtain a total instantaneous thermal power consumption value, such as 25W. This process is repeated for each subsequent 10ms time slice, potentially yielding values ​​such as 28W, 22W, etc. Continuing these consecutive instantaneous thermal power consumption values ​​constitutes the expected heat flux time series for the task executed on this core, such as sequence data points 25, 28, 22, etc.

[0023] In this step, not only is the current physical state (temperature) obtained, but the historical heat accumulation is also quantified by the integral value of the previous cycle. More importantly, by extracting task features (such as floating-point density), the system is given the ability to predict the future load thermal characteristics, providing a data foundation for subsequent feedforward control.

[0024] Step S2: Based on the task feature vector and the preset instruction set thermal power consumption model, predict the expected heat flux time series of the task set to be scheduled on each core; calculate the thermal penalty factor that is positively correlated with the cumulative temperature error integral value of each computing core in the previous control cycle; generate and execute the task scheduling scheme with the goal of balancing the expected heat flux of each core and the total execution time of the thermal penalty factor.

[0025] Beforehand, a series of benchmark programs are run on the target server, and power consumption and instruction execution are measured simultaneously. Using methods such as multiple linear regression, an instruction set thermal power consumption model is established, which maps task feature vectors to core power consumption. For example, thermal power consumption equals weight coefficient A multiplied by floating-point operation density, plus weight coefficient B multiplied by memory access frequency, plus a basic power consumption constant. For any task in the task set to be scheduled, the task feature vector is substituted into the model to calculate the expected thermal power consumption value when executed on a single core. This value is the expected heat flux. If the task execution time is known, a step signal with the amplitude of this heat flux value can be generated during the task duration as the expected heat flux time series.

[0026] For each computing core in the server, the cumulative temperature error integral value stored in the previous cycle is obtained, and the integral value is multiplied by a preset proportional constant to obtain the thermal penalty factor of the core. For example, if the cumulative temperature error integral value of a core in the previous cycle is large, it indicates that the core is in a high-temperature state, and the calculated thermal penalty factor is also correspondingly large.

[0027] Construct a multi-objective optimization function with two objectives: first, to minimize the expected heat flux of the largest task in the core after considering the thermal penalty factor; and second, to minimize the total completion time of the task set, i.e., the maximum completion time. The expected heat flux after considering the thermal penalty factor is equal to the sum of the expected heat fluxes of all tasks allocated to that core multiplied by one plus the thermal penalty factor of that core. Use heuristic algorithms such as genetic algorithms or particle swarm optimization to solve the optimization problem and obtain a mapping scheme for allocating tasks in the task set to specific computing cores. Call the task scheduling interface at the operating system level to bind tasks to the specified cores for execution according to this scheme.

[0028] To add a scheduling penalty to the continuously overheating core, in one optional embodiment, the thermal penalty factor is calculated as follows: ,in, For the first One core thermal penalty factor, For the first The cumulative temperature error integral value of the previous cycle of each core. This is the preset penalty gain coefficient.

[0029] For example, suppose there are two cores, core 1 and core 2, and a preset penalty gain coefficient. The value is 0.05. After a control cycle ends, check the cumulative temperature error integral value of each core from the previous cycle. If core 1 has been operating under high load for a long time, the cumulative temperature error integral value from the previous cycle will be 0.05. A value of -50 degrees Celsius per second indicates that the temperature remains above the target value.

[0030] Calculate the thermal penalty factor of core 1 according to the formula. The thermal penalty factor of core 1 is obtained. The value is 3.5. For core 2, if the integral value of the previous cycle is +10 degrees per second, it indicates that the temperature is lower than the target value, and the penalty factor for core 2 is calculated accordingly. A penalty factor of 1.5 will reduce the priority of new tasks being assigned to core 1 in subsequent task scheduling decisions.

[0031] In an optional embodiment, a task scheduling scheme is generated and executed with the goal of balancing the expected heat flux of each core factor that includes the thermal penalty factor with the total execution time, including: Construct the optimization objective function ,in, The total completion time for the task set. For the first One core thermal penalty factor, To be allocated to the first The expected total heat production of all tasks in each core, and For preset weighting coefficients, and These are the normalized time and normalized heat, respectively, used to make the components of the objective function dimensionless; The particle swarm optimization algorithm is used to iteratively optimize the allocation scheme of tasks on each core until the objective function is achieved. The value of converges to the minimum, thus obtaining the task scheduling scheme.

[0032] Specifically, assuming three tasks need to be scheduled to two cores, the weighting coefficients are... and The values ​​are set to 0.6 and 0.4 respectively. Based on the thermal penalty factors calculated above (thermal penalty factor of core 1),... It is 3.5, with 2 cores. (1) The particle swarm optimization algorithm randomly generates multiple scheduling schemes as initial particles. For example, scheme A assigns tasks 1 and 2 to core 1, and task 3 to core 2. The total completion time is calculated based on this scheme. The total heat generated by core 1 is allocated within 100ms. For 500 Total heat production of core 2 For 200 .

[0033] Calculate the objective function of scheme A Value. Assuming normalized time. 120ms, normalized heat For 1000 The objective function can be calculated. The value is 1.28. The algorithm will continue to evaluate other options, such as option B, which assigns task 1 to core 1 and tasks 2 and 3 to core 2, and calculate... Value. Through multiple iterations, the particle swarm optimization algorithm continuously adjusts task allocation to find the optimal value. The option with the smallest value is selected and adopted as the implementation option because it represents the best balance between speed and thermal management.

[0034] In this step, by predicting heat flux and introducing a thermal penalty factor, the scheduling algorithm no longer focuses solely on performance but actively avoids "hot" cores, distributing the thermal load evenly across all computing units. This effectively prevents localized hotspots caused by task backlog, reduces the pressure on subsequent cooling systems, and improves chip reliability.

[0035] Step S3: Calculate the temperature error between the current temperature and the target temperature; set an integral time step that is inversely proportional to the absolute value of the temperature error change rate, and calculate the cumulative temperature error integral value of the current period using the trapezoidal integral method; use the expected heat flux time series and the current temperature to predict the temperature change trajectory within the future time window using the server's first-order heat conduction model, and use the initial slope of the trajectory as the differential term.

[0036] The steps for calculating the cumulative temperature error integral value for the current period include: Calculate the current temperature error Temperature error compared to the previous cycle The difference is used to obtain the change in temperature error. ; According to the formula Calculate the integration time step ,in, Based on the time step, This is the preset integral step size adjustment coefficient; Through formula Calculate the cumulative temperature error integral value for the current period. ,in, This is the cumulative temperature error integral value from the previous cycle.

[0037] For example, assume a base time step The integral step size adjustment coefficient is 10ms. It is 0.5. Period, measured temperature error It is 2℃. At present... During the cycle, due to the start of a high-load task, the temperature rose rapidly, resulting in a measured temperature error. The temperature is 8℃. At this point, the change in temperature error... The temperature is 6℃. Due to the drastic temperature changes, a smaller integration step size is needed to improve accuracy.

[0038] The integration time step can be calculated using the formula. The integration time step is 2.5ms, which is relatively short. The cumulative temperature error integral value is updated using the trapezoidal integration method. If the cumulative value from the previous cycle... If the time interval is 20 milliseconds, then the cumulative value for the current period is... It is 32.5 degrees per millisecond. If the temperature change is gradual, for example... The calculated value is 0.2℃. This will be close to 10ms, thus reducing unnecessary calculations while maintaining accuracy.

[0039] Using the expected heat flux time series and the current temperature, a server-side first-order heat conduction model is used to predict the temperature change trajectory within a future time window, with the initial slope of the trajectory as the differential term, including: Establish the first-order heat conduction differential equation for the server core: ,in, For the core temperature, For the expected heat flux, For ambient temperature, and These are the pre-calibrated model coefficients; The current temperature Initial values ​​of the expected heat flux time series Substituting into the equation, the initial slope is calculated. ; The calculated initial slope As the differential term.

[0040] Specifically, a predicted value about future temperature trends is provided to the PID controller, assuming that the thermal model coefficients of a server core have been pre-calibrated, and the heat generation coefficient... The heat dissipation coefficient is 0.15. The value is 0.03. Current ambient temperature. The core's current temperature is 25°C, as measured by a sensor. The temperature is 75℃. Based on the upcoming scheduled task, the expected heat flux time series for a period of time is predicted, with the first value being... It is 40W.

[0041] The initial slope is calculated by substituting the real-time and predicted data into the first-order heat conduction differential equation. .get The value is 4.5. The unit of the obtained value 4.5 is °C / s, which represents the rate of change of the core temperature under the current state and expected load. The differential term... This input is fed into the derivative of the PID controller, enabling the controller to respond to the expected temperature rise trend in advance, rather than waiting for the temperature error to actually occur before making adjustments.

[0042] In this step, by adopting an integration step size that is inversely proportional to the rate of change of error, details can be captured when the temperature fluctuates drastically, while integration saturation can be avoided in steady state. The differential term is calculated using a first-order heat conduction model, and the trend is predicted by using physical laws rather than simply data differences, which greatly improves the ability to predict future temperature changes and enhances the ability to resist noise interference.

[0043] Step S4: Calculate the proportional gain coefficient that is positively correlated with the absolute value of the expected heat flux and the cumulative temperature error integral value of the current cycle; generate proportional, integral and differential control quantities based on the proportional gain coefficient, the cumulative temperature error integral value and the differential term, and sum them to generate a heat dissipation control command for adjusting the speed of the server cooling fan.

[0044] Specifically, a base proportional gain coefficient is set, plus a positive coefficient multiplied by the expected heat flux of the task to be executed on the current core, plus another positive coefficient multiplied by the absolute value of the integral value of the cumulative temperature error in the current cycle, and the sum of the three is used as the proportional gain coefficient; when a high-power task is about to be executed or the historical cumulative temperature difference is large, the proportional gain coefficient will increase accordingly, making the control effect stronger.

[0045] Multiply the calculated proportional gain coefficient by the current temperature error to obtain the proportional term; multiply a fixed preset integral gain coefficient by the cumulative temperature error integral value of the current period to obtain the integral term; multiply a fixed preset differential gain coefficient by the differential term to obtain the differential term. Summate the calculated proportional, integral, and differential control values ​​algebraically to obtain a comprehensive control output value; convert the comprehensive control output value into a pulse width modulation (PWM) duty cycle signal for the corresponding fan using a preset mapping function or lookup table method, for example, an integer from 0 to 255; send this PWM signal to the fan controller via the server's baseboard management controller (BMC) to adjust the actual fan speed.

[0046] In an optional embodiment, the formula for calculating the proportional gain coefficient is: ,in, This is the proportional gain coefficient. The basic proportional gain coefficient, This represents the average expected heat flux within the future time window. This represents the absolute value of the integral of the accumulated temperature error for the current period. and This is the preset gain adjustment coefficient.

[0047] Specifically, adjust the proportional gain in the PID controller. It can adapt to different workloads and system thermal states. Set the base proportional gain coefficient. The gain adjustment coefficient is 0.8. It is 0.01. The value is 0.002. Within one control period, the average heat flux for a future time window is predicted. A value of 50W indicates an upcoming high-load scenario. Calculate the cumulative temperature error integral value for the current cycle. -60 degrees per second, absolute value The value of 60 indicates that the core has been deviating from the target temperature for some time.

[0048] Calculate the proportional gain coefficient according to the formula. Larger This value allows the PID controller to respond more strongly and quickly to the current temperature error. An aggressive control strategy is necessary under high loads and existing temperature deviations to prevent temperature runaway. Conversely, under low loads and stable temperatures, the calculated... The value will be smaller, making the control effect more gentle.

[0049] This step achieves dynamic matching of control intensity. The proportional gain coefficient is positively correlated with the predicted heat flux and cumulative error, meaning that when facing an upcoming high load or significant historical heat buildup, the system will output a more aggressive control input to quickly suppress temperature rise; while maintaining gentle control under low loads. This non-linear adjustment mechanism optimizes fan efficiency and reduces unnecessary noise and vibration while ensuring that the temperature does not exceed the limit.

[0050] In a second embodiment, the present invention also provides a server heat dissipation control device, comprising the following modules: The acquisition module is used to acquire the current temperature of each computing core of the server, the cumulative temperature error integral value of the previous control cycle, and the task feature vector of the task set to be scheduled. The task feature vector includes the instruction type composition ratio, floating-point operation density and memory access frequency of the task. The generation module is used to predict the expected heat flux time series of the task set to be scheduled on each core based on the task feature vector and the preset instruction set thermal power consumption model; calculate the thermal penalty factor that is positively correlated with the cumulative temperature error integral value of each computing core in the previous control cycle; and generate and execute the task scheduling scheme with the goal of balancing the expected heat flux of each core and the total execution time, which includes the thermal penalty factor. The prediction module is used to calculate the temperature error between the current temperature and the target temperature; it sets an integral time step that is inversely proportional to the absolute value of the temperature error change rate, and uses the trapezoidal integral method to calculate the cumulative temperature error integral value of the current period; using the expected heat flux time series and the current temperature, it predicts the temperature change trajectory within the future time window through the server's first-order heat conduction model, and uses the initial slope of the trajectory as the differential term. The calculation module is used to calculate a proportional gain coefficient that is positively correlated with the absolute value of the expected heat flux and the cumulative temperature error integral value of the current period; based on the proportional gain coefficient, the cumulative temperature error integral value and the differential term, it generates proportional, integral and differential control quantities and sums them to generate a heat dissipation control command for adjusting the speed of the server cooling fan.

[0051] The present invention also provides a computer-readable storage medium storing a computer program, which, when executed by a processor, implements the server heat dissipation control method described above.

[0052] The system also includes other components well known to those skilled in the art, such as communication buses and communication interfaces, the settings and functions of which are known in the art and therefore will not be described in detail here.

[0053] In the description of this specification, "multiple" means at least two, such as two, three or more, etc., unless otherwise expressly and specifically defined.

[0054] While various embodiments of the invention have been shown and described in this specification, it will be apparent to those skilled in the art that such embodiments are provided by way of example only. Many modifications, alterations, and alternatives will occur to those skilled in the art without departing from the spirit and essence of the invention.

Claims

1. A method for regulating server heat dissipation, characterized in that, Includes the following steps: The current temperature of each computing core of the server, the cumulative temperature error integral value of the previous control cycle, and the task feature vector of the task set to be scheduled are obtained. The task feature vector includes the instruction type composition ratio, floating-point operation density and memory access frequency of the task. Based on the task feature vector and the preset instruction set thermal power consumption model, predict the expected heat flux time series of the task set to be scheduled on each core. Calculate the thermal penalty factor that is positively correlated with the cumulative temperature error integral value of each computing core in the previous control cycle; With the goal of balancing the expected heat flux of each core component and the total execution time by taking into account the thermal penalty factor, a task scheduling scheme is generated and executed. Calculate the temperature error between the current temperature and the target temperature; set an integral time step that is inversely proportional to the absolute value of the rate of change of temperature error, and calculate the cumulative temperature error integral value of the current period using the trapezoidal integral method; use the expected heat flux time series and the current temperature to predict the temperature change trajectory within the future time window through the server's first-order heat conduction model, and use the initial slope of the trajectory as the differential term. Calculate a proportional gain coefficient that is positively correlated with the absolute value of the expected heat flux and the cumulative temperature error integral value of the current cycle; based on the proportional gain coefficient, the cumulative temperature error integral value and the differential term, generate proportional, integral and differential control quantities and sum them to generate a heat dissipation control command for adjusting the speed of the server cooling fan.

2. The server heat dissipation control method according to claim 1, characterized in that, The step of predicting the expected heat flux time series of the scheduled task set executed on each core includes: Each task in the task set is broken down into a sequence of micro-operations within a unit time slice; Query the preset instruction set thermal power consumption model to obtain the basic thermal power consumption value corresponding to each micro-operation; Based on the proportion of instruction types, floating-point operation density, and memory access frequency in the task feature vector, the basic thermal power consumption values ​​of each micro-operation are weighted and summed to obtain the instantaneous thermal power consumption within the time slice. By concatenating the instantaneous thermal power values ​​within consecutive time slices, the expected thermal flux time series executed on the corresponding core is obtained.

3. The server heat dissipation control method according to claim 1, characterized in that, The formula for calculating the thermal penalty factor is as follows: ,in, For the first One core thermal penalty factor, For the first The cumulative temperature error integral value of the previous cycle of each core. This is the preset penalty gain coefficient.

4. The server heat dissipation control method according to claim 1, characterized in that, The process of generating and executing a task scheduling scheme with the goal of balancing the expected heat flux of each core component and the total execution time, which are included in the thermal penalty factor, includes: Construct the optimization objective function ,in, The total completion time for the task set. For the first One core thermal penalty factor, To be allocated to the first The expected total heat production of all tasks in each core, and For preset weighting coefficients, and These are the normalized time and normalized heat, respectively, used to make the components of the objective function dimensionless; The particle swarm optimization algorithm is used to iteratively optimize the allocation scheme of tasks on each core until the objective function is achieved. The value of converges to the minimum, thus obtaining the task scheduling scheme.

5. The server heat dissipation control method according to claim 1, characterized in that, The step of calculating the cumulative temperature error integral value for the current period includes: Calculate the current temperature error Temperature error compared to the previous cycle The difference is used to obtain the change in temperature error. ; According to the formula Calculate the integration time step ,in, Based on the time step, This is the preset integral step size adjustment coefficient; Through formula Calculate the cumulative temperature error integral value for the current period. ,in, This is the cumulative temperature error integral value from the previous cycle.

6. The server heat dissipation control method according to claim 1, characterized in that, The step of using the expected heat flux time series and the current temperature to predict the temperature change trajectory within a future time window through a server-side first-order heat conduction model, and using the initial slope of the trajectory as the differential term, includes: Establish the first-order heat conduction differential equation for the server core: ,in, For the core temperature, For the expected heat flux, For ambient temperature, and These are the pre-calibrated model coefficients; The current temperature Initial values ​​of the expected heat flux time series Substituting into the equation, the initial slope is calculated. ; The calculated initial slope As the differential term.

7. The server heat dissipation control method according to claim 1, characterized in that, The formula for calculating the proportional gain coefficient is as follows: ,in, This is the proportional gain coefficient. The basic proportional gain coefficient, This represents the average expected heat flux within the future time window. This represents the absolute value of the integral of the accumulated temperature error for the current period. and This is the preset gain adjustment coefficient.

8. A server heat dissipation control device, characterized in that, Includes the following modules: The acquisition module is used to acquire the current temperature of each computing core of the server, the cumulative temperature error integral value of the previous control cycle, and the task feature vector of the task set to be scheduled. The task feature vector includes the instruction type composition ratio, floating-point operation density and memory access frequency of the task. The generation module is used to predict the expected heat flux time series of the task set to be scheduled on each core based on the task feature vector and the preset instruction set thermal power consumption model. Calculate the thermal penalty factor that is positively correlated with the cumulative temperature error integral value of each computing core in the previous control cycle; With the goal of balancing the expected heat flux of each core component and the total execution time by taking into account the thermal penalty factor, a task scheduling scheme is generated and executed. The prediction module is used to calculate the temperature error between the current temperature and the target temperature; it sets an integral time step that is inversely proportional to the absolute value of the temperature error change rate, and uses the trapezoidal integral method to calculate the cumulative temperature error integral value of the current period; using the expected heat flux time series and the current temperature, it predicts the temperature change trajectory within the future time window through the server's first-order heat conduction model, and uses the initial slope of the trajectory as the differential term. The calculation module is used to calculate a proportional gain coefficient that is positively correlated with the absolute value of the expected heat flux and the cumulative temperature error integral value of the current period; based on the proportional gain coefficient, the cumulative temperature error integral value and the differential term, it generates proportional, integral and differential control quantities and sums them to generate a heat dissipation control command for adjusting the speed of the server cooling fan.

9. A computer-readable storage medium storing a computer program thereon, characterized in that, When the computer program is executed by the processor, it implements the server heat dissipation control method as described in any one of claims 1-7.

Images