Wall Heat Flux Research Articles

Numerical study on surface corrosion deposition of fuel elements and its influence on flow heat transfer

Corrosion of pressurized water reactors (PWR) in nuclear power plants can lead to serious safety hazards. This study aims to analyze the deposition of corrosion products using FLUENT software. Deposition models and thermal resistance models were developed, and the effects of deposits on the reactor’s thermal–hydraulic characteristics were evaluated. Additionally, the impact of various parameters on deposition and thermal–hydraulic characteristics was examined. Results show that deposits accumulate extensively in the inlet section of the fuel cladding, while appearing as spot deposits in the outlet section. For deposit thicknesses below 30 μm, the surface temperature of the cladding gradually increases. However, when the thickness exceeds 30 μm, the surface temperature rapidly rises. Furthermore, the study reveals that the deposition amount decreases with increasing inlet flow velocity, exhibits an upward trend with higher inlet temperature, and increases with a higher wall heat flux density. This research provides important insights for understanding core deposition and thermal–hydraulic characteristics in nuclear reactor systems. It offers valuable guidance for enhancing safety and operational efficiency in nuclear power plants.

Numerical study of film cooling in single-element injector gaseous CH4/O2 rocket engine with coupled wall function

The investigation of film cooling in CH4/O2 rocket engines holds paramount importance in the advancement of rocket propulsion. However, the wall heat flux is always overestimated in numerical simulation processes. Hence, this article proposes a numerical framework that employs the Reynolds averaged Navier–Stokes method to simulate the single-element gaseous CH4/gaseous O2 combustion chamber, which is the basis for rocket engine simulations. The coupled wall function that considers chemical reaction effects is introduced to enhance the accuracy of wall heat flux prediction. The impact of utilizing a coupled wall treatment on the prediction of wall heat flux and its fundamental parameters are examined. In addition, a single-element combustion chamber experiment is performed to validate the simulation. The results demonstrate that the implementation of the coupled wall function hardly influences the main flow characteristics, whereas the wall heat flux calculated with general wall functions shows an overestimation, which can be reduced by the coupled wall function. Finally, the case with a coupled wall function can improve the cooling efficiency with greater accuracy and the cooling systems with optimized design. The proposed numerical framework and the findings of the study provide profound insights that can improve the design and optimization of rocket engines.

Viscous dissipation effects on slip flow heat transfer in rhombic microchannels

This study aimed to numerically investigate the viscous dissipation effect on forced convection in rhombic microchannels for gases in a slip flow regime. The numerical analysis was carried out by assuming a 2D steady-state flow. The solution of governing equations was obtained by adopting the finite element method and assuming that the fluid is Newtonian, with constant thermophysical properties, and in a fully developed laminar flow regime. The solution of the momentum equation is obtained by considering a first-order boundary condition, while the solution of the energy balance equation is obtained by assuming a constant wall heat flux (H2 boundary condition) and taking into account the wall temperature jump. The validation of the numerical model was carried out using the data available in the scientific literature. The numerical outcomes obtained for several values of the acute angle of the rhombus, the Knudsen number (i.e., rarefaction effects), and the Brinkman number (i.e., viscous dissipation effects) reveal that viscous forces play an opposite role with respect to rarefaction and significantly affect the convective heat transfer coefficient, especially for low Knudsen numbers and for high values of the acute angle. In particular, it was observed that Nu is significantly affected by the Brinkman number for acute angles higher than 50°.

Flow boiling pressure drop correlation in small to micro passages

An accurate and acceptable correlation for the prediction of two-phase pressure drop is considered a crucial step in heat exchanger design. Although many existing models and correlations were developed and proposed in the past, their ability to predict within acceptable error bands when applied in general to flow boiling flows is limited, even within their originally recorded ranges. The discrepancies are worse when predicting two-phase pressure drop in small to micro-scale passages. Therefore, the aim of the work described in this paper is to assess the most well-known models and correlations using a large experimental data-bank. The data-bank includes four refrigerants, namely R134a, R245fa, HFE-7100 and HFE-7200, small to micro heat exchangers made of different metals and different channel configurations. In addition, the data cover a large range of operating conditions, which can allow generally applicability of new correlations developed. This range covers channel hydraulic diameter of 0.46−4.26 mm, heated length of 20−500 mm, system pressure of 1 − 14 bar, mass flux of 50−700 kg/m2 s, wall heat flux ranging from 2 to 234 kW/m2 and exit vapour quality up to one. Twenty six existing models and correlations that were developed and proposed for vertical/horizontal flow, single/multi-channels and circular/non-circular channels were evaluated. Moreover, the effect of using different equations for calculating the two-phase mixture viscosity, void fraction, Fanning friction factor and Lockhart–Martinelli parameter was assessed and discussed. The mean absolute error of the existing correlations when compared with our data-bank was more than 30 %. Therefore, a new correlation for calculating the two-phase multiplier, which was strongly dependent on the Boiling number and the Lockhart–Martinelli parameter, was developed and then the frictional component and total two-phase pressure drop relationships were completed.

Numerical investigation on boiling crisis in a full-length 5 × 5 fuel assembly under typical pressurized water reactor conditions

Critical heat flux (CHF) presents a research priority in the Pressurized Water Reactor (PWR) safety, as it can lead to the severe accidents and release radioactive materials. In this study, a wall heat flux partition model, proposed by Demarly [1], was employed within the two-fluid model by Computational Fluid Dynamics (CFD) methodology to investigate the characteristics of CHF under typical PWR conditions. A multi-scale interface model was coupled to consider the interfacial transfers in bulk flow during flow regime transition. The models were validated by using experimental data on subcooled flow boiling, boiling crisis in vertical tubes, and rod bundles. The PWR sub-channel and bundle test (PSBT) benchmark, comprising a full-length 5 × 5 fuel assembly with three different types of spacers, was simulated in this study. The CHF value was evaluated, demonstrating an acceptable accuracy compared to the experimental data, with a deviation of less than 10 % in the selected conditions. The CHF location in the cross section was successfully predicted by CFD approach, while the axial location of CHF was underestimated than experimental measurement. Besides, the impacts of mixing vanes with different deflection angles on boiling heat transfer characteristics were analyzed, revealing the crucial role of mixing vanes in the distribution of vapor and dry-patch area fraction in both the circumferential and flow directions.

Combustion Characteristics and Application Performance of JP-10/Nano-Sized Aluminum Gel-Fueled Scramjet

Using nano-sized energetic materials as fuel additives to hydrocarbon gel fuel has gained in popularity due to the advantages of higher energy density and shorter ignition delay times. Combustion characteristics and application performance of gel fuel containing 16% wt aluminum nanoparticles were investigated on a direct-connect scramjet platform with an inflow Mach number of 2, an inflow total temperature of 1700 K, and a fuel equivalence ratio ranging from 0.6 to 1.0. The gel fuel could be stably injected and atomized efficiently by a nitrogen-assisted injector. The combustion efficiency of JP-10+Al gel fuel was 4.27% higher than that of JP-10 gel fuel, while the addition of aluminum increased the density specific impulse by 6.34%, indicating that the addition of aluminum nanoparticles indeed presented a positive effect on combustion and engine performance. The wall deposition phenomenon of gel fuel was particularly slight as compared to that of slurry fuel cases. Unlike nanoaluminum-blended slurry fuel, nanoaluminum-blended gel fuel controlled the potential risk of increasing thermal load to the least degree within the scope of existing research, with wall heat flux increasing by 7.43% solely as a result of elevated airflow temperature, indicating its advantageous application prospect.

A new wall function method for hypersonic laminar boundary layers

A new wall function method for hypersonic laminar boundary layers (HLBLs) is proposed to reduce the near-wall grid dependence of skin friction $c_f$ and wall heat flux $q_w$ in numerical simulations, aiming for fast and accurate predictions. First, an analytic laminar velocity law of the wall is derived, which achieves a universal scaling of the near-wall velocity of HLBLs. Then an accurate temperature–velocity relation is deduced by introducing the general recovery factor to address the invalidation of the Walz relation under the cold wall effect. Based on the laminar laws of the wall, a new wall function method for HLBLs is proposed. To avoid introducing the boundary layer edge quantities, the laminar laws of the wall are reformed by modifying the outer boundary conditions of the differential equation in deriving the temperature–velocity relation. Unlike the wall function method in turbulence, the new wall function obtains directly the accurate $c_f$ and $q_w$ by post-processing without being involved in the simulation iteration. The numerical experiments of a Mach 8 HLBL over the flat plate show that effectively, the new wall function can enlarge the distance of the first grid point off the wall $\Delta y_1$ from $10^{-6}$ m to $10^{-3}$ m, which brings a 50 times enhancement of the simulation efficiency. Meanwhile, the simulation errors of $c_f$ and $q_w$ of the mesh with $\Delta y_1=10^{-3}$ m are reduced significantly from 24.2 % and 18.5 % to 0.5 % and 0.1 %, respectively. Due to the new wall function removing the boundary layer edge quantities, success is also achieved under the curved walls.

Analysis of the relationship between near-wall velocity distribution and wall heat flux under diesel spray flame impingement in an RCEM

To clarify the relationship between the near-wall flow and the heat flux through the wall, the velocity distribution of the diesel spray flame and the heat flux at the combustion chamber wall were measured in a rapid compression and expansion machine (RCEM), which has a twodimensional combustion chamber. The velocity distribution of the diesel spray flame near the wall was measured by the particle image velocimetry (PIV), and the heat flux through the wall where the diesel spray flame impinges was measured using a newly developed multi-point heat flux sensor. Furthermore, the velocity gradient tensor was calculated based on the measured velocity distribution results, and the effect of the shear flow characteristics on wall heat transfer was analyzed. As a result, the higher the injection pressure leads the higher the peak value of the heat flux. A fluctuation of the heat flux exhibits the same trend as the fluctuations of the velocity, and a high heat flux appeared in the high-velocity gradient tensor region, which suggests that the characteristics of shear flow near the wall can affect the heat flux.

Numerical Model of Simultaneous Multi-Regime Boiling Quenching of Metals

This work presents a heat transfer and boiling model that computes the evolution of the temperature field in a representative steel workpiece quenched from 850 or 930 °C by immersion in water flowing at average velocities of 0.2 or 0.6 m/s, respectively. Under these conditions, all three boiling regimes were present during cooling: stable vapor film, nucleate boiling, and single-phase convection. The model was based on the numerical solution of the heat conduction equation coupled to the solution of the energy and momentum equations for water. The mixture phase approach was adopted using the Lee model to compute the rates of water evaporation–condensation. Heat flux at the wall was calculated for all regimes using a single semi-mechanistic model. Therefore, the evolution of boiling regimes at every position on the wall surface was automatically determined. Predictions were validated using laboratory results, namely: (a) videorecording the upward motion of the wetting front along the workpiece wall surface; and (b) cooling curves obtained with embedded thermocouples in the steel probe. Wall heat flux calculations were used to determine the importance of the simultaneous presence of all three boiling regimes on the heat flux distribution. It was found that this simultaneous presence leads to high heat flux variations that should be avoided in production lines. In addition, it was determined that the corresponding inverse heat conduction problem to estimate the active heat transfer boundary condition must be set-up for 2D heat flow.

Controlling gas–liquid flow and enhancing heat transfer in a T-junction microchannel by wettability-engineered walls

Characteristics of gas–liquid flow and heat transfer in a cross-flow T-junction microchannel with wettability-engineered walls are numerically investigated in this paper. The validated diffuse interface method is adopted for interface capture. First, the effects of wall wettability on bubble formation and transportation are studied. Three flow patterns are observed due to different combinations of the bottom and the top wall contact angles. On this basis, two methods are proposed to enhance the heat transfer. One is to increase the two-phase interfacial contact area by dividing the microchannel into three functional regions, which can promote the heat exchange at the two-phase interface. The other is to increase the velocity fluctuation intensity by alternating the contact angle along the channel, which can enhance mixing between the hot liquid layer adjacent to the wall and the cool liquid core. These two methods are applicative for steady and unsteady problems, respectively. The flow states, velocity vectors, and streamlines are used to analyze the fluid and thermal mixing mechanism. Meanwhile, a quantitative comparison of the wall temperature is made at a given wall heat flux. The obtained results can provide fresh insights into the gas–liquid flow control and the heat transfer enhancement in a microchannel, which are valuable for the design of microreactors and radiators.

Experimental study of the effects of isolated roughness elements on the stability and transition of a hypersonic boundary layer on a flat plate

Surface roughness elements on hypersonic vehicles can cause early boundary layer transition, increasing wall skin friction and heat flux and affecting aircraft range and thermal protection systems. Accurate prediction of the transition caused by these roughness elements is crucial for the design of hypersonic vehicles. In this work, wind tunnel experiments on isolated roughness-induced boundary layer transition at Ma = 6 are conducted. Infrared thermography and high-frequency pressure sensors are utilized to investigate the effects of different roughness element configurations (cylindrical, diamond, ramp) on the hypersonic boundary layer instability and transition. The experimental results show that all three roughness elements can effectively enhance the generation of second mode waves and promote boundary layer transition. Compared to smooth surfaces, they exhibit similar frequency band range, faster growth, and earlier saturation. Among them, the ramp roughness element most effectively triggers the boundary layer transition, with a relatively small heat flux increase. Furthermore, bispectral analysis illustrates that all three roughness elements undergo self-interactions that lead to spectral broadening, ultimately resulting in boundary layer transitions.

Infinitely Fast Heterogeneous Catalysis Model for Premixed Hydrogen Flame-Wall Interaction

In most computations of reacting flows, walls are considered as inert. This approach can produce non-physical results that lead to very large values of heat release at the wall for hydrogen Flame-Wall Interaction (FWI). This paper confirms this unexpected behavior at walls for H2 flames. It also explores the possibility that the problem does not come from a deficiency of the chemical scheme at low temperatures but from the wall boundary condition: due to the large amount of radical species produced by H2 flames, the “inert” assumption may not be valid anymore. The purpose of this study is twofold: (1) to propose a simplified model named Infinitely Fast Heterogeneous Catalysis (IFHC) that mimics surface reactions and (2) to test it on the canonical Head-On Quenching (HOQ) configuration. A large set of parameters is investigated: a skeletal (San Diego) and a detailed (Burke) chemical schemes are tested, the grid-dependence and the influence of operating conditions (ϕ, Tw and P) are assessed. Results show that the IFHC model can inhibit the chemical pathways that are responsible for the heat generation peak observed on the wall in hydrogen FWI with inert surfaces. The IFHC model also allows computations to converge when the mesh is refined, a property that cannot be reached for inert wall treatments. Using IFHC, neither the predicted quenching distance nor the flame structure are altered. However, the maximum wall heat flux is significantly decreased when compared to inert walls. Over the range of operating conditions studied, the maximum wall heat flux is reached at high pressure, low wall-temperature and equivalence ratios close to the maximum of reactivity of H2 flames (ϕ≈1.7). This study shows that, for H2 FWI, walls should be at least modeled with simple surface reactions – such as IFHC – and that assuming inert walls is not possible.Novelty and significance statementThe work presented in this paper is new, original and of interest at the present time as many groups simulate near-walls H2/Air flames. The IFHC model is relevant for all premixed hydrogen combustion close to walls and is, thus, relevant for most realistic premixed hydrogen combustion scenarios. Besides being conservative in terms of mass, momentum, energy and atom balance, IFHC allows to reach mesh convergence during quenching contrary to simplistic inert wall treatments. It also allows to describe walls at an affordable computational cost compared to a full set of heterogeneous reactions. This work also opens the discussion with experimentalists since, to the best of the authors’ knowledge, the dimensional maximum wall heat fluxes are not yet available in the literature and are, thus, provided here for future comparison.

Opposing Mixed Convection Heat Transfer for Turbulent Single-Phase Flows

Convection, wherein forced and natural convections are prominent, is known as mixed convection. Specifically, when a forced convection flow is downward, this flow is called opposing flow. The objectives of this study are to gain a comprehensive understanding of opposing flow mixed convection heat transfer and to establish the prediction methodology by evaluating existing correlations and models. Several heat transfer correlations have been reported related to single-phase opposing flow; however, these correlations are based on experiments conducted in various channel geometries, working fluids, and thermal flow parameter ranges. Because the definition of nondimensional parameters and their validated range confirmed by experiments differ for each correlation reported in previous studies, establishing a guideline for deciding which correlation should be selected based on its range of applicability and extrapolation performance is important. This study reviewed the existing heat transfer correlations for turbulent opposing flow mixed convection and the single-phase heat transfer correlations implemented in the thermal–hydraulic system codes. Furthermore, the authors evaluated the predictive performance of each correlation by comparing them with the experimental data obtained under various experimental conditions. The Jackson and Fewster, Churchill, and Swanson and Catton correlations can accurately predict all the experimental data. The effect of the difference in the thermal boundary conditions, i.e., uniform heat flux and uniform wall temperature, on the turbulent mixed convection heat transfer coefficient is not substantial. The authors confirmed that heat transfer correlations using the hydraulic-equivalent diameter as a characteristic length can be used for predictions regardless of channel-geometry differences. Furthermore, correlations described based on nondimensional dominant parameters can be used for predictions regardless of the differences in working fluids. The authors investigated the extrapolation performance of the mixed convection heat transfer correlations for a wide range of nondimensional parameters and observed that the Jackson and Fewster, Churchill, and Aicher and Martin correlations exhibit excellent extrapolation performance with respect to natural and forced convection flows, indicating that they can be applied beyond the parameter range validated experimentally.

Computational Fluid Dynamics and Its Potential Applications for the ENT Clinician.

This article is an examination of computational fluid dynamics in the field of otolaryngology, specifically rhinology. The historical development and subsequent application of computational fluid dynamics continues to enhance our understanding of various sinonasal conditions and surgical planning in the field today. This article aims to provide a description of computational fluid dynamics, the methods for its application, and the clinical relevance of its results. Consideration of recent research and data in computational fluid dynamics demonstrates its use in nonhistological disease pathology exploration, accompanied by a large potential for surgical guidance applications. Additionally, this article defines in lay terms the variables analyzed in the computational fluid dynamic process, including velocity, wall shear stress, area, resistance, and heat flux.

Steady and Unsteady Complex Heat Transfer in Optically Thick Medium During Film Boiling

Abstract The paper presents the results of a study of radiative-convective heat transfer at film boiling of a liquid on a vertical heated plate. Both a steady-state problem of heat transfer and a transient problem were considered. The latter describes the instantaneous (flash) boiling up of a liquid on a heated surface. The novelty of the present study is the use of the optically thick medium approximation in a mathematical model when studying the process of radiation-convective heat transfer in the film boiling regime. For the first time, radiation heat transfer was considered for an optically thick medium. An analytical solution of the steady-state problems is obtained for boundary conditions involving a constant wall temperature and a constant wall heat flux. The effect of radiation and the temperature difference between the wall and liquid on the temperature profiles in the vapor phase is shown. The effect of radiation becomes more pronounced with an increase in the temperature difference between the wall and the liquid. As a result of solving the transient problem, the variation in time of the temperature profile and the heat transfer coefficient in the vapor film were obtained. The effect of radiation (Stark number) on the heat transfer coefficients is elucidated. An increase in the radiative heat flux leads to an increase in the Nusselt number, as well as the time it takes for the heat transfer process to reach a steady-state regime.

Investigation of forced convection of a novel hybrid nanofluid containing NEPCM in a square microchannel: Application of single-phase and Eulerian-Eulerian two-phase models

In the present study, laminar forced convection flow and heat transfer of novel hybrid nanofluid of Al 2 O 3 -n-octadecane/water through a horizontal microchannel under a constant wall heat flux condition are investigated. N-octadecane is regarded as a nano-encapsulated phase change material (NEPCM). Simulations are performed utilizing single-phase and Eulerian-Eulerian two-phase models. Governing equations are discretized by finite volume method and then solved using SIMPLE algorithm. Effects of nanoparticles volume fraction ( φ ), Reynolds number ( Re ), and melting range ( T mr ) on thermal and hydrodynamic parameters of hybrid nanofluid are evaluated. It is shown that predictions by Eulerian-Eulerian approach are more accurate than those by single-phase approach when compared to available experimental data. Addition of Al 2 O 3 and n-octadecane nanoparticles to water-based fluid leads to a reduction in fluid and wall temperatures, which is intensified at higher φ and Re . Also, these nanoparticles do not have a significant effect on hydrodynamic development of flow, while they delay its thermal development. NEPCM with a higher T mr only slightly improves thermal performance of hybrid nanofluid, which is more noticeable in single-phase procedure. Total efficiency of hybrid nanofluid ( η hnf ) is amplified by raising φ for both approaches. This augmentation is more evident in Eulerian-Eulerian model due to better thermal and hydrodynamic performances. Maximum increases in η hnf concerning single-phase and Eulerian-Eulerian models are 27% and 60%, respectively, occurring at φ Al 2 O 3 = 1 % and φ n ‐ octadecane = 10 % . Therefore, this type of hybrid nanofluids, because of its enhanced thermo-physical properties, can be considered as a suitable candidate for thermal energy storage, heat transfer, and renewable energy applications.

Experimental study on shock interaction control of double wedge in high-enthalpy hypersonic flow subject to plasma synthetic jet

The hypersonic shock-shock interaction flow field at double-wedge geometries controlled by plasma synthetic jet actuator is experimentally studied in a Ma = 8 high-enthalpy shock tunnel with the purpose of exploring a novel technique for reducing surface heat flux in a real flight environment. The results demonstrate that increasing the discharge energy is advantageous in eliminating the shock wave, shifting the shock wave interaction point, and shortening the control response time. The oblique shock wave can be completely removed when the actuator's discharge energy grows from 0.4 J to 11.5 J, and the displacement of the shock wave interaction point increases by 124.56%, while the controlled response time is shortened by 30 μs. Besides, the reduction in diameter of the jet exit is firstly proved to have a negative impact on energy deposition in a working environment with incoming flow, which reduces the discharge energy and hence decreases the control effect. The shock wave control response time lengthens when the jet exits away from the second wedge. Along with comparing the change in wall heat flux at the second wedge over time, the control effect of plasma synthetic jet actuator with and without inflation is also analyzed. When plasma synthetic jet works in inflatable mode, both the ability to eliminate shock waves and the shifting effect of the shock wave interaction point are increased significantly, and the wall heat flux is also reduced.

Numerical Investigations of the Thermal-Hydraulic Characteristics of Microchannel Heat Sinks Inspired by Leaf Veins

A microchannel heat sink (MCHS) is a potential solution for chip and battery thermal management. The new microchannel structure is beneficial for further improving the thermal-hydraulic performance of MCHSs. Inspired by leaf veins, six new channel structures were designed, and the effects of the channel structures (three parallel structures named PAR I, II, and III and three pinnate structures named PIN I, II, and III), channel depths (0.4, 0.8, and 1.6 mm), and heat fluxes (20, 50, and 80 kW/m2) were investigated via numerical simulation. The cooling medium was water, and the heating area was 40 × 40 mm2. Both PAR II and PIN III exhibit superior overall performance, characterized by the highest Nusselt number and the lowest heating wall temperature. Moreover, PIN III demonstrates the lowest standard deviation in heating wall temperature, while PAR II exhibits the lowest friction factor. The greater the channel depth is, the larger the solid–liquid contact area is, leading to a reduced wall temperature at the interface under identical conditions of inlet Reynolds number and heating wall heat flux. Consequently, an increase in the Nusselt number corresponds to an increase in the friction factor. The maximum value and standard deviation of the heating wall temperature increase with increasing heat flux, while the Nusselt number and friction factor remain unaffected. The overheating near the two right angles of the outlet should be carefully considered for an MCHS with a single inlet–outlet configuration.

Adjoint-based linear sensitivity of a supersonic boundary layer to steady wall blowing–suction/heating–cooling

For a Mach $4.5$ flat-plate adiabatic boundary layer, we study the sensitivity of the first, second Mack modes and streaks to steady wall-normal blowing/suction and wall heat flux. The global instabilities are characterised in frequency space with resolvent gains and their gradients with respect to wall-boundary conditions are derived through a Lagrangian-based method. The implementation is performed in the open-source high-order finite-volume code BROADCAST and algorithmic differentiation is used to access the high-order state derivatives of the discretised governing equations. For the second Mack mode, the resolvent optimal gain decreases when suction is applied upstream of Fedorov's mode $S$ /mode $F$ synchronisation point, leading to stabilisation, and the converse when applied downstream. The largest suction gradient is in the region of branch I of mode $S$ neutral curve. For heat-flux control, strong heating at the leading edge stabilises both the first and second Mack modes, the former being more sensitive to wall-temperature control. Streaks are less sensitive to any boundary control in comparison with the Mack modes. Eventually, we show that an optimal actuator consisting of a single steady heating strip located close to the leading edge manages to damp the linear growth of all three instability mechanisms.

Water drops on nanofiber-coated substrates: Influence of wall temperature and coating thickness on hydrodynamics and wall heat flux distribution

Heat and mass transport during drop impact or a gentle deposition onto a heated substrate depends on the surface chemistry, morphology, thermal and mechanical properties of the substrate, as well as the substrate temperature. In particular, if a substrate is coated with a porous layer of a wettable material, the drop spreading is accompanied by the imbibition of a liquid into the layer. The wetted area is significantly enhanced in comparison with the area which can be covered by a drop spreading over a bare substrate. As a result, the liquid evaporation rate increases and the evaporation time decreases. Understanding the interaction between the liquid spreading, imbibition and evaporation is important, both for enhancement of heat transfer in cooling applications and for functionalizing of porous media. In our experimental studies, a disk of infrared-transparent calcium fluoride is used as the base substrate. The disk is coated with a submicrometer layer of black chromium nitride and an electrically conductive chromium layer, used for heating of the substrate by electrical current. Nanofiber coating layers are applied on top of the chromium layer using electrospinning. The dynamics of drop spreading is captured by a high-speed camera in a side view, the imbibition is observed with a top view camera, and the temperature distribution at the substrate-coating interface is captured by an infrared camera. Based on the transient temperature field, corresponding heat flux distributions are determined. For the detailed analysis of the interplay between hydrodynamics and heat transfer, the nanofiber coating thickness and substrate temperature have been varied. It has been found that the drop drying time and maximal imbibed area depend non-monotonously on the coating thickness. Changing the initial wall temperature leads to qualitatively different distribution of the heat flux at the substrate surface.