Dynamic Temperature Distribution Research Articles

Transient thermal management of laser systems using Plate-Fin phase change heat Exchangers: Experimental and computational study

To mitigate transient thermal shocks in lasers and reduce thermal stresses caused by temperature fluctuations, the use of phase change materials (PCMs) in thermal management systems is a viable solution. This study proposes an innovative two-dimensional transient heat transfer model specifically designed for plate-fin phase change heat exchangers (PFPCHEs). The model meticulously simulates the complex heat transfer phenomena within the heat exchanger, including fluid convection, solid thermal conduction, and the phase change processes of PCM. The convective heat transfer coefficient between fluid and plate is calculated using the Wieting correlation. An advanced pulse-mode experimental test platform was constructed to validate the model, dynamically monitoring and recording outlet temperatures and heat exchange performance for stringent experimental comparison. The research explores the impact of key operating parameters such as initial temperature, flow rate, and inlet temperature of the cooling cycle on the performance of the heat exchanger, providing valuable design and control strategies for transient thermal management of laser systems. The experimental results confirm that the model accurately predicts the dynamic response characteristics and temperature distribution of PFPCHEs under multi-cycle pulse loads, with 96% of the predictions within a 10% error margin. A significant finding is that the initial temperature has negligible influence on the heat transfer characteristics when it is below the solid phase temperature of the PCM. Moreover, by meticulously adjusting the flow rate and inlet temperature of cooling cycle, it is possible to effectively maintain the stability of the outlet temperature throughout the entire pulse cycle. This is crucial for minimizing temperature fluctuations within the thermal management system and extending the service life of electronic components.

Calculation method for the dynamic temperature distribution of winding area in dry-type on-board traction transformer

Dry-type on-board traction transformers have attracted considerable attention because of their lightweight and flameproof characteristics. The temperature distribution of winding areas is a critical indicator for representing operation performance. This study developed a calculation model for the dynamic temperature distribution of winding areas in dry-type on-board traction transformers under time-varying environments. First, a dynamic temperature calculation model was proposed based on the principle of thermoelectric analogy. Next, for addressing the transient outlet temperature of air ducts in the model, the influences of environment parameters and air duct size on heat dissipation were investigated, and a calculation model for the transient outlet temperature was proposed. Finally, the validity of the proposed dynamic temperature calculation method was verified experimentally, and the method was compared with conventional computational fluid dynamics numerical simulation using line data. The results revealed that for various structural parameters and time-varying operation scenarios, the proposed method could accurately obtain the temperature distribution in the winding area and could save more than 99% of the time cost compared with the conventional method.© 2017 Elsevier Inc. All rights reserved.

Dynamic temperature distribution characteristics of a large glasshouse with cooling system during the start-stop stage

Excessive heat accumulation leads to high temperature in the glasshouse, the high temperature directly affects the normal growth of crops. This research based on the cooling system as the cooling method of glasshouse, experimental and computational fluid dynamics technique were used to investigate the cooling effect in the start-stop stage. CFD was used to simulate and optimize the temperature dynamics during the start-stop phase of the glasshouse cooling system in the most common summer climates (rainy and sunny). A glasshouse test platform was built in Wuhan in order to conduct temperature experiments to verify the effectiveness of the simulation. The results show that starting the cooling system in summer can effectively reduce the temperature inside the glasshouse in a short time (when it is rainy, the temperature is reduced by 8.1 °C; when it is sunny, the temperature is reduced by 5 °C). During the start-stop stage, the temperature trend of the horizontal and vertical planes is the same. Until the system is running stably, the indoor temperature is maintained uniform and stable, with little fluctuation. In order to create the best growing environment for crops, the air velocity of the air supply outlet is increased to 6.0 m/s on sunny days, which is helpful to increase crop yield and reduce greenhouse gas emissions.

Heat Transfer Analysis of MgB2 Coil in Heat Treatment Process for Future Fusion Reactor

State of the art MgB <sub xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">2</sub> is reviewed as a potential material for the poloidal field (PF) coils of the future fusion reactor due to its high critical temperature and low material cost. The heat treatment process is a crucial step in the development of MgB <sub xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">2</sub> magnets. The temperature lag in heat treatment of large magnets can lead to insufficient thermal reaction time. It may be infeasible to control the temperature of a magnet according to the heat treatment scheme recommended for the MgB <sub xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">2</sub> wire. Hence, the heat treatment process of a large magnet needs to be evaluated. Therefore, the dynamic temperature distribution of a MgB <sub xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">2</sub> PF coil is obtained by simulating the heat transfer in heat treatment process. A suitable heat treatment schedule for a large magnet is proposed and the experimental results of a sub-size Cable-In-Conduit Conductor manufactured with MgB <sub xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">2</sub> strand confirmed the feasibility of the newly proposed heat treatment process. The results provide a reference for the heat treatment method of a future larger MgB <sub xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">2</sub> coil.

Numerical analytical model of metal temperature of steam boiler drum

The problem of developing and identifying a temperature distribution model in the wall of the boiler drum was solved. The numerical model is based on the analytical solution of the one-dimensional heat conduction equation with heat exchange at the boundaries at constant heat transfer coefficients. The model is implemented in a dynamic systems simulation package. Using the bounded series of eigenfunctions of the analytic solution, a representation of the plant in the state space was obtained. A universal algorithm for calculating the eigennumbers of the analytic solution was proposed. The eigenvalues dependences of the analytical solution on the values of heat transfer coefficients and thermal conductivity was investigated. Based on the results of the study, data tables were obtained to approximate the dependences of eigennumbers on heat transfer coefficients in a limited range of variation. The methodology for implementing a dynamic temperature distribution model, which takes into account changes in heat transfer coefficients during simulation, was tested. The time of change of coefficients is perceived as the beginning of a new transient from the current state of the model. Presented structural diagrams of the implementation of the model in dynamic systems simulation packages. The proposed structures allow the introduction of parametric disturbances at the time of their occurrence in the simulation. According to the experimental data of steam boiler firing, identification of the dependence of the heat transfer coefficient on the temperature of the metal wall of the boiler drum was carried out. The results showed that the heat transfer coefficients are influenced not only by the metal temperature, but also by additional parameters of the steam boiler operation. As a result of experiments, it was found that the accuracy of the model can be improved by taking into account the dependence of the heat transfer coefficient on the internal surface of the drum on the flow rate of feed water into the boiler drum. The results of the numerical experiment are presented.

Thermal circuit model of prismatic lithium cell considering dynamic non-uniform characteristics during charging-discharging in energy storage

Current analytical and simulation models for lithium battery thermal behaviour encounter efficiency or accuracy challenges in energy storage applications. In this paper, an analytical thermal analysis approach for prismatic lithium cells considering dynamic non-uniform characteristics is proposed to calculate the dynamic temperature distribution. A thermal circuit model and analytical differential equation group were formulated to capture the dynamic non-uniform thermal characteristics during charging/discharging. A solving algorithm was provided to balance accuracy and efficiency by establishing an iterative mechanism over the modelling parameters. The proposed model and algorithm were applied to three commercial lithium cells. The experiments under various charge-discharge conditions verified the millisecond-level computational cost and less than 3 % calculation accuracy. The parameter calibration, iteration mechanism, and reproducible codes enhanced the engineering applicability. The performance of the approach suggested its excellent potential in battery thermal management for energy storage.

Non-uniform heat generation model of pouch lithium-ion battery based on regional heat generation rate

The non-uniform distribution of temperature seriously affects the cycle performance, cycle life and safety of the pouch lithium-ion batteries. In this paper, a dynamic non-uniform heat generation model is developed to the high complexity and low precision of the modeling and characterization of dynamic non-uniform thermal characteristics of the pouch battery. Firstly, the regional heat generation rate of different regions is estimated by the particle swarm optimization based on the cell surface temperature to build a non-uniform heat generation model. Then, the battery's temperature distribution characteristics and heat generation distribution characteristics under the effect of ambient temperature, state of charge and discharge rate are analyzed. Furthermore, combined with the three-dimensional heat transfer equation, a three-dimensional heat model with non-uniform a heat source is modeled by COMSOL, the dynamic non-uniform temperature distribution of the battery is simulated at different ambient temperatures and discharge rates, and then the local temperature at any location of the battery is obtained. Finally, the simulated value of the local temperature is compared with the experimental value, and the average error under constant current discharge condition is less than 1 °C, and the relative error under dynamic current condition is less than 3 %. The result shows that the model can accurately estimate pouch lithium-ion batteries' dynamic non-uniform heat generation distribution under different working conditions.

Applications of distributed fiber Bragg gratings to the measurements of in-tube fluid temperature distribution

The temperature distribution information of the two-phase fluid inside a tube can effectively reflect the heat transfer of the fluid, which is the key information in the study of the heat transfer of flowing fluid in a tube. This article aims to develop a method for measuring the temperature distribution of fluid inside a tube based on distributed fiber Bragg gratings (DFBGs). It is more convenient than traditional measurement methods using thermocouples that are impossible to be installed to the inner tube of a double tube heat exchanger. By using a stainless steel seamless tube, ϕ1.0 × 0.2 mm in size, and a hose to seal and protect the optical fiber with Bragg grating sensors inside, we successfully shielding the stress effects from the thermal effects, and also enhance the durability of the fiber. Since the seal leads to extra thermal resistance to the DFBGs, obtaining the response performance of the sealed DFBGs to a temperature change is essential for reliable applications. The performance characterization shows that the temperature measurement repeatability of the sealed DFBGs is within 0.2 °C for a period up to 35 months. The time constants of the DFBGs are not longer than 0.25 s, capable of measuring the temperature oscillations with frequency less than 2 Hz. Comparison with commercial thermocouples supports the feasibility and the reliability that the sealed DFBGs can be used as temperature sensors. To obtain the equivalent thermal conductivity between the tube that sealed the fiber and the fiber, κ¯air, the response to a temperature jump was also calculated based on a simplified model, which can comprehensively reflect the heat transfer of the air inside the sealed tube and provide the reference for design optimization of a sealed FBG for dynamic measurements. Finally, the capillary tube with sealed DFBGs was placed in the condenser tube of a two-phase circulation loop to exhibit the feasibility of the steady and dynamic temperature distribution measurements of the fluid inside a tube.

Moore–Gibson–Thompson Stability Model in a Two-Temperature Photonic Semiconductor Excited Medium Affected by Rotation and Initial Stress

In this paper, the two-temperature theory is used to examine a novel model that generalizes the Moore–Gibson–Thompson (MGT) effect according to two-dimensional electronic/thermoelastic deformation. The main equations for a semiconductor medium in the context of the impact of rotation are explained in terms of the impact of the initial hydrostatic stress at the free surface. The normal-mode approach is used to derive the precise formulae for the fundamental physical quantities (i.e., normal displacement, normal load stress, electronic diffusion (carrier density), dynamic and conductive temperature distribution) under the influence of the two-temperature coefficient. The comparison with the base state is performed using linear stability analysis. To make some comparisons based on the various values of thermal memories, the influence of a number of novel parameters is applied to each of our primary physical quantities, such as the rotation parameter and the initial stress. An example of the main fields’ perturbation is also obtained and graphically described.

Real-time temperature predictions via state-space model and parameters identification within rack-based cooling data centers

Current rack-based cooling architecture of data centers (DCs) is a promising method since it simplifies airflow distribution and provides fast cooling regulation. Real-time on-demand control of the cooling system is an effective way to improve its operational energy efficiency without sacrificing the thermal security of IT equipment. Accurate and fast temperature distribution prediction serves as one of the bases for ensuring the superior performance of advanced control algorithms. Thus, this study proposed a novel grey-box state-space model, to rapidly predict the dynamic temperature distribution for rack-based cooling DCs. Various heat transfer physics in the rack-based cooling DCs, including heat production caused by servers and heat movements by airflows, were modeled as a state-space structure using the zonal modeling approach. The coefficient matrices were identified through the prediction-error method (PEM), in order to avoid the extremely time-consuming process of obtaining accurate physical parameters regarding the system. This developed model was validated with an experimentally validated mechanistic model and computational fluid dynamics (CFD) simulations. Additionally, the impact of the prediction horizon's size and IT workload transient changes on the proposed model's prediction accuracy were investigated as well. Through simulation, the developed model achieves sufficient accuracy with an average root mean square error (RMSE) equal to 0.19 °C and less than 3% relative error for predicting 1-min dynamic temperature evolutions. Also, the developed model shows satisfying performances for the long-horizon prediction and transient scenario, which will facilitate advanced control techniques for DC cooling systems.

On the Real-Time Quantification of Flexibility Provided by District Heating Networks Considering Dynamic Temperature Distribution

District heating networks (DHNs) can provide electrical power systems (EPSs) flexibility by utilizing heat storage abilities. There are three challenges to the real-time quantification of the flexibility offered by the DHN. First, the balance between the heat supply and consumption in the long term must be considered. Second, the robustness and economy of the quantified result must be guaranteed. Third, the reliable dynamic temperature distribution in the DHN must be obtained. Thus, this paper designs a flexibility-providing mechanism for the DHN to interact with the EPS, consisting of the service period and the recovery period. The real-time quantification of the flexibility of the DHN is formulated as an optimization problem, where both flexibility revenue and the operation cost increase in the CHP units are considered. A novel thermal dynamic state estimation (TDSE) model with a bad data filtering process is proposed, which is used to provide the dynamic temperature distribution to support the real-time quantification of the flexibility. Numerical results verified that the proposed real-time quantification model has advantages in ensuring the benefits of the DHN. The performance of the TDSE in providing data support for flexibility quantification is also verified. Finally, key factors impacting the quantified flexibility are exploited by simulations.

Thermal FEM Analysis of Surge Arresters During HVdc Current Interruption Validated by Experiments

This paper deals with the development of an accurate finite-element model of an arrester to investigate the electrothermal and mechanical stress during dc current interruption. The comprehensive analysis performed on a ZnO surge arrester is supported by experiments during high-voltage dc circuit breaker current interruption. The performed experimental analysis comprises three sequential 26 kV/10 kA direct current interruption tests carried out within a period of one hour. The dynamic temperature and current distribution of the surge arrester columns during current interruption are measured. The finite-element simulation results are in good agreement with the test results. The influence of the surge arrester temperature on the current distribution among the surge arrester columns is analyzed. The impact of the surge arrester temperature on ZnO electrical characteristics and mechanical stress inside the surge arrester are also investigated. The surge arrester finite-element model can be used with full success for parameter optimization of the surge arresters to prevent possible failures when dc circuit breakers performed multiple interruptions in short period of time.

A high-fidelity real-time capable dynamic discretized model of proton exchange membrane fuel cells for the development of control strategies

Advanced modeling tools are widely applied to enhance control strategies of proton exchange membrane fuel cells. Existing dynamic fuel cell models, however, still face a conflict between computational speed and simulation accuracy. To overcome this contradiction, a novel discretized fuel cell model is presented based on existing segmented modeling approaches. With a new twofold subdivision method and a newly designed interpolation algorithm, this model allows accurate investigation of the dynamic temperature, two-phase water, and current density distribution. Compared to the finite element method, the simulation duration of hours is reduced to minutes, while the characteristic of high-fidelity is maintained. Model validations are performed on both cell-level and stack-level to demonstrate the model's high accuracy and real-time simulation capability, illustrating the high potential for applying the model in the development of fuel cell control strategies.

Parametric optimization of a cesaro fins employed latent heat storage system for melting performance enhancement

The large-scale use of sustainable energy necessitates the use of latent heat storage (LHS). This study aims to increase the melting performance of an LHS system by designing and optimizing the cesaro fins. To examine the melting behaviours of PCM (i.e., RT-82) in a finned LHS system, a 2-D melting heat transfer model is developed and numerically solved. As per literature, there is very literature including the NC (Natural Convection) and represents the coordination in fins, nanoparticles and metal foam. So, the impacts of natural convection, fin arrangement, nanoparticles and metal foam are explored for the Fourier number variations from 0.014 to 0.158 at a constant Stefan number of 0.20. The dynamic temperature distribution, as well as the natural convection and fin arrangement, are investigated to determine how phase change material melts over time by considering the non-thermal equilibrium between metal foam and PCM/nanoPCM. The findings illustrate that natural convection has a significant effect on melting behaviours in the LHTES (Latent Heat Thermal Energy Storage) system, with a 26.8% increase in melting/charging rate as compared to the situation without NC. The LHTES system with improved fin structure (i.e., Type-3) and increased number of fins (i.e., Type-5) configurations have more uniform temperature distribution and a higher melting rate by increasing the heat transfer cooperation between NC and thermal conduction. Low thermal conductivity causes poor performance in PCM (Phase Change Material) energy storage devices. The melting/charging of PCM in an LHTES system is greatly improved in this study by employing a porous metal foam (i.e., copper metal foam) or nanoparticle (i.e., Cu, CuO and Al2O3). The impact of nanoparticle volume fraction and metal foam porosity on the LHTES system's melting/charging performance is also investigated using the enthalpy-porosity technique. According to the results, PCM's melting/charging time is lowered by 63.4% when nanoparticles are mixed in and incorporated with metal foam for the Fourier number varies from 0.014 to 0.158 at a constant Stefan number of 0.20. PCM melting/charging time decreased when the metal foam porosity decreased or increased in the volume fraction of nanoparticles. High porosity metal foams with low volume fractions of nanoparticles can increase melting performance since it assures minimum PCM volume and increases natural convection.

Evaluating the thermal characteristics of laser powder bed fusion

• Dynamic temperature behaviour around the melt pool during the PBF-LB/M was evaluated. • Temperature distribution was influenced by morphological changes of metal powder. • Laser beam incident was characterized by direct heating and heat conduction. • Dynamic temperature behaviour provided a direct explanation for the cooling rate. This study investigates the dynamic temperature behaviour around a melt pool in metal-based powder bed fusion using a laser beam (PBF-LB/M) to clarify the influence of the associated morphological changes of the metal powder experimentally. Gas-atomized 18Ni (300-grade) maraging steel powders were processed by PBF-LB/M while high-speed photography with a two-colour radiometric thermal imaging system that was employed to correlate the temperature with melt pool behaviour. In addition, the cooling rate of the melt pool was measured directly using the dynamic temperature distribution. The temperature distribution of the melt pool was influenced by the morphological changes of the metal powder induced by physical and thermal interactions, and the melt pool exhibited an asymmetric temperature distribution in the direction parallel to the laser scan. The significant factors were droplet cohesion at low melt pool temperatures, remaining heat energy from previous laser beam irradiation, and the heat conduction inside the melt pool. The laser beam incident on the metal powder was primarily characterized by two modes: direct heating induced by laser beam irradiation and heat conduction through the single track, droplets, and substrate. In addition, the dynamic temperature behaviour provided a direct explanation for the cooling rate, the values of which ranged from 0.1 to 0.9 × 10 6 K/s owing to the self-cooling induced by PBF-LB/M.

Improved Temperature Monitoring and Protection Method of Three-Level NPC Application Based on Half-Bridge IGBT Modules

The overheating stress of power devices is the main cause of converter system failures, so real-time temperature acquisition, as well as timely over-temperature protection, under various operating conditions are key to achieving a more reliable energy conversion. In this paper, an improved three-dimensional (3-D) coupled thermal model for insulated gate bipolar transistors (IGBTs) and related heatsink is constructed by considering the thermal coupling effects at both the device- and the module level. Furthermore, the coupled thermal impedance parameters are thereby extracted based on the finite-element method (FEM), so that the dynamic loss and transient temperature distribution under actual working conditions can be obtained based on this electro-thermal coupling model. Meanwhile, to avoid the overheating damage of power devices due to untimely action of the negative temperature coefficient (NTC) thermistor over-temperature protection, this article aims at the idle IGBT device in the three-level neutral-point clamped (NPC) application using IGBT half-bridge modules, combined with the emitter terminal temperature, and finally realizes the multiple over-temperature protection settings. The improved temperature monitoring methods proposed in this paper are realized by the drive and control circuits and verified through simulation and experimental results.

Numerical modeling on thermal runaway triggered by local overheating for lithium iron phosphate battery

The popularization of electric vehicles is no longer a speculation thanks to the advancement of the lithium-ion battery. However, safety issues prevent a larger-scale application of lithium-ion batteries. In energy storage industry, thermal runaway is the direct cause of fire and explosions. The thermal runaway numerical model attempts to unveil the mechanism of thermal runaway triggered by local overheating in this work, which is distinct from accelerating rate calorimeter tests or oven tests. The model is based on the energy conservation equation, and the source term of which is modeled by the empirical Arrhenius equations. Further on, the results of the numerical model are validated by experiments. Based on the calculated analysis on the heat contribution to the thermal runaway, the heat contribution of anode-electrolyte reaction is up to 63.8% of the total heat, and internal short circuit is responsible for a little proportion of total heat generation during thermal runaway which is contrary to the general thoughts, yet still a crucial inducement of chain reactions. The results also confirm the conclusion that enhancement of heat dissipation condition of battery can retard thermal runaway progress and reduce the maximum thermal runaway temperature. In addition, interior dynamic temperature distribution of the battery can be achieved, and characteristics of local overheating are compared with those of oven tests, noteworthy temperature difference in local overheating occasion is highlighted.

A simulating method of moisture continuous diffusion under changing temperatures and analysis of moisture-induced stresses covering moisture desorption and reflow processes for the QFN

Moisture in the atmosphere could diffuse into plastic packaging devices, which would cause high hygro-mechanical stress and vapor stress and lead irreversible damages. In this paper, an analysis method for moisture continuous diffusion under changing temperature is proposed. Considering the vapor pressure induced strain as hygro-mechanical strain, a corresponding relationship between the parameters of hygro-mechanical stress and vapor stress is built. Subsequently, the formula of moisture-induced total stress is derived according to the superposition principle. A QFN (Quad Flat No‑lead) package model with gold wires inside is established to validate the effectiveness of the method. The moisture continuous diffusion from the absorption process to the desorption one is obtained, then the dynamic variations of the hygro-mechanical stress, vapor stress and moisture-induced total stress of the QFN during the desorption process and the entire reflow process are obtained for the first time ever. The results prove that the method in this paper can be employed for the analysis of moisture diffusion and moisture-induced stress under dynamic temperature and uneven humidity distribution.

Mathematical modeling of thermal processes in building enclosures

In the development of automatic control systems (ACS) for heat supply of buildings, mathematical models of their control objects play a significant role, one of the most important elements of which being the enclosing structures. The existing models of thermal processes indoors either do not take into account the dynamic effect of the enclosing structures, or are complex analytical expressions, the use of which in the engineering practice of the synthesis of ACS for the indoor thermal regime is difficult. The purpose of the study was to create a linear mathematical model of the thermal interaction of the external and internal air environment through the enclosing structures. To achieve this goal, a functional diagram of a heated room as a control object was developed, which describes the mechanism of thermal interaction of its elements. On its basis, design schemes of thermal processes for external and internal enclosing structures were developed. They were used to create systems of equations linking temperatures on the surfaces of layers of external enclosing structures and air temperatures outside and inside the premises. At the same time, the assumption was introduced that the dynamic temperature distribution over the thickness of the enclosing structure is linear. To find the temperatures at the boundaries of the layers, the average temperature value over the thickness of the structure was used. As a result of the study, equations for the dynamics of heat transfer processes were obtained, and a generalized structure of a model of multilayer external enclosing structures of a room as an element of a control object was developed. The resulting model can be used to develop a generalized mathematical model of the thermal regime of both individual rooms and the building as a whole.

Demonstration of treatment planning software for hyperthermic intraperitoneal chemotherapy in a rat model

Background Hyperthermic intraperitoneal chemotherapy (HIPEC) is administered to treat residual microscopic disease after cytoreductive surgery (CRS). During HIPEC, fluid (41–43 °C) is administered and drained through a limited number of catheters, risking thermal and drug heterogeneities within the abdominal cavity that might reduce effectiveness. Treatment planning software provides a unique tool for optimizing treatment delivery. This study aimed to investigate the influence of treatment-specific parameters on the thermal and drug homogeneity in the peritoneal cavity in a computed tomography based rat model. Method We developed computational fluid dynamics (CFD) software simulating the dynamic flow, temperature and drug distribution during oxaliplatin based HIPEC. The influence of location and number of catheters, flow alternations and flow rates on peritoneal temperature and drug distribution were determined. The software was validated using data from experimental rat HIPEC studies. Results The predicted core temperature and systemic oxaliplatin concentration were comparable to the values found in literature. Adequate placement of catheters, additional inflow catheters and higher flow rates reduced intraperitoneal temperature spatial variation by −1.4 °C, −2.3 °C and −1.2 °C, respectively. Flow alternations resulted in higher temperatures (up to +1.5 °C) over the peritoneal surface. Higher flow rates also reduced the spatial variation of chemotherapy concentration over the peritoneal surface resulting in a more homogeneous effective treatment dose. Conclusion The presented treatment planning software provides unique insights in the dynamics during HIPEC, which enables optimization of treatment-specific parameters and provides an excellent basis for HIPEC treatment planning in human applications.