Heat Transfer Module Research Articles

Thermal propagation analysis of 2G HTS stacked wires based on a dimensional coupling method

Stacked structures of high-temperature superconducting (HTS) tapes are usually employed to enhance current-carrying capacity in cable or magnet applications, but their multilayer characteristics may result in local hot spots due to uneven cooling. Besides, the inconsistencies along the length of tapes introduced during manufacturing process and the significant transverse Lorentz forces experienced in magnet operations can lead to weak parts, bringing uncertainties to thermal stability. Finite-element methods are widely used to predict thermal propagations, but 3D models encounter challenges such as distorted mesh elements and convergence issues, particularly when combining electromagnetic and heat transfer modules for long-distance wires. In this work, we have developed a Dimensional Coupling Method (DCM) to assess the thermal impact of weak parts in stacked wires applying coupled 1D and 2D models. The 2D model analyzes the electromagnetic and heat characteristics of stacked surfaces, and provides an initial heat source for the 1D model, which evaluates thermal propagation longitudinally. Simulation results of the 1D module are then transferred back to update the 2D outcomes. Models of distinct dimensions are coupled sequentially in physical steps but simultaneously in the time domain. Our approach is verified by 3D model benchmarks and offers a computational cost reduction of approximately 60 % compared to the benchmarks, making it more suitable for applications with large Iop/Ic ratios. Specially, multi-layer stacked wires under low ratios cases are also been analyzed. What’s more, two influencing factors of heat propagation, the weak-part length and position, are also investigated.

Effect of insulation materials and power on microwave heating characteristics of SiC susceptor under material–specific parametric conditions

High-quality microwave absorber materials play a crucial role in the efficient and uniform processing of variety of materials through microwave hybrid heating. Silicon carbide (SiC) is an excellent microwave absorber at room temperature, which enhances the heating rate of the target materials significantly. In the present study, effects of heat insulation materials and their configurations on microwave heating characteristics of SiC were investigated at different power levels. The responses were assessed in terms of the heating capability (maximum temperature achieved), radiation heat loss, temperature gradient, and heating uniformity under different parametric conditions. A 3D COMSOL Multiphysics simulation model, coupled with microwave heating and heat transfer module, was developed with surface-to-ambient radiation boundary conditions. Electromagnetic and thermal properties of the SiC utilized in the model were selected based on physical properties, such as porosity/density, particle size, and chemical composition of SiC. The model outputs were validated with experimental findings. Results showed that microwave heating characteristics were highly sensitive to the susceptor (SiC) properties, insulation materials and power levels. Insulation materials prevented heat losses significantly; at the same time, different insulation configurations altered the electrical field intensity distribution within the SiC sample, which affected the responses. Microwave power levels also influenced the responses and quantity of insulation required.

Thermal modeling, simulation, and experimental validation of micro electric discharge machining of tungsten carbide (WC)

Tungsten carbide is a material that is widely used in tools and dies applications due to its high hardness, excellent toughness, and wear resistance. However, machining this material using conventional methods can be very challenging due to its exceptional properties. In such cases, electrical discharge machining can be used as an alternative method for machining this material. The aim of this study is to examine the thermal properties of EDM of tungsten carbide, both through finite element modeling and experimental conditions. The study aims to evaluate the temperature distribution in tungsten carbide during EDM die sinking. This is accomplished by drilling 0.5 mm diameter and 1 mm deep micro holes using micro-drilling processes. Seven experiments were performed using an EDM die-sinking machine, and a 3D axis symmetrical model was created and simulated using the COMSOL Multiphysics heat transfer module. The temperature profile of tungsten carbide material during a single spark machining was obtained by considering only 30% of energy in the form of Gaussian distribution that gets transferred to the workpiece. The temperature profile of this model was then used to estimate the material removal rate. By comparing the numerical and experimental results, it was found that there was only a 4.8% average percentage error, indicating very close agreement between the numerical and experimental results. The findings suggest that the finite element method of the COMSOL Multiphysics heat transfer module can accurately predict and simulate the real-time results of EDM machining. These results could be useful for EDM machine operators for pre-estimating MRR and maximum temperatures before proceeding with the actual machining.

Analytical and Experimental Behaviour of GFRP-Reinforced Concrete Columns under Fire Loading

Fire engineering endeavours to mitigate injury or the loss of life in the event of a fire. This is achieved primarily through fire prevention, containment, and extinguishment measures. Should prevention fail, the structural integrity of buildings, coupled with effective evacuation strategies, becomes paramount. While glass fibre-reinforced polymer (GFRP) materials have demonstrated efficacy in reinforcing concrete elements, their performance under fire conditions, notably in comparison to steel, necessitates a deeper understanding for structural applications. This study experimentally and numerically investigates the fire performance of GFRP-reinforced concrete (RC) columns subjected to only fire load without additional external loads. The research aims to ascertain the fire resistance based on the thickness of the concrete coating and the ultimate tensile strength of GFRP rebars after 90 min of fire exposure. Four GFRP-RC columns were subjected to a standardized fire curve on all sides in the experimental program. In the analytical program, a theoretical model was developed using the heat transfer module of the COMSOL software. The results of both analyses were very close, indicating the reliability of the procedure used. Based on the findings, recommendations regarding the fire resistance of GFRP-RC columns were formulated for structural applications. Results from this research provide the civil engineering community with data that will help them continue using FRP materials as internal reinforcement for concrete.

Numerical simulation of temperature field of a new type of freezing device under seepage effect.

In order to investigate the development of the temperature field of a new type of freezing reinforcement under seepage conditions, in this paper, COMSOL finite element software was used to simplify the model and simulate the effect of groundwater seepage on the development of the temperature field of frozen pipes by coupling the Darcy's law module and the heat transfer module for porous media. The heads of water were also varied to simulate the change in seepage velocity to further investigate the effect of seepage velocity on the temperature field. The results of the study show that the freezing wall formed in the high head region was thinner than that in the low head region due to the effect of seepage, and this phenomenon was aggravated with the increase of seepage rate; The effect of seepage action on the temperature field had a hysteresis along the seepage direction; When the seepage rate was greater than 1.65 m/d, the soil in the center of the device feezed better and could form a tight and dense freezing wall comparable to the size of the freezing device; When the seepage rate was greater than 5.78 m/d, the temperature of the center soil body gradually increased, and eventually the freezing curtain cannot be formed.

Numerical Simulation and Modeling of Mechano–Electro–Thermal Behavior of Electrical Contact Using COMSOL Multiphysics

Electrical contacts are important circuit components with diverse industrial applications, and their failure can lead to multiple unwanted effects. Hence, the behavior of electrical contacts is a widely studied topic in the scientific literature based on various approaches, tools, and techniques. The present study proposes a new approach to numerical modeling and simulation based on the Holm contact theory, aiming to study the dependence between the electric potential and the temperature within an electrical contact. Structured in five sections, the research was conducted using COMSOL Multiphysics software and its solid-state mechanics, electric current, and heat transfer modules in order to highlight contact behavior from mechanical, electrical and thermal points of view: the von Mises stress, contact force, electric field amplitude, variation of the electrical potential along the current path, temperature gradient, and dependence of temperature along the contact elements edges were obtained by simulation, and are graphically represented. The results show that the temperature increase follows a parabolic curve, and that for values higher than 4 mV of voltage drop, the temperature of the contact increases to 79.25 degrees (and up to 123.81 degrees for 5 mV) over the ambient temperature, thus the integrity of insulation can be compromised. These values are close (10–12%) to the analytically calculated ones, and also in line with research assessed in the literature review.

Chiral photothermoelectric effect driven by high-Q quasi-bound states in the continuum

Bound states in continuum (BIC) have been proposed as a means to efficiently improve light-matter interactions in metasurfaces. While breaking the mirror symmetry of structure and developing BIC into a reachable and observable quasi-BIC, it is usually accompanied by the chiral phenomenon with high quality (Q) factor. Here, we report a spin-sensitive photodetector in the infrared (NIR) region comprising a silicon metasurface with chiral quasi-BIC, a silver layer, and a thermoelectric layer. A chiral quasi-BIC supported by a silicon metasurface can be realized under normal incidence. Based on finite element method simulations, a silicon metasurface with a silver layer shows a high Q-factor of 958.6 with a giant absorption circular dichroism of 0.83. Subsequently, we study the thermal performance of the chiral absorbers by using the heat transfer module of COMSOL. Combined with the thermoelectric material bismuth telluride, we calculate the differential photothermoelectric effects of the system under circular-polarized light irradiation. When the incident flux is 100 W cm−2, the output voltage under right-circular-polarization (left-circular-polarization) light reaches 0.59 mV (0.08 mV), which can be used for polarization detection. Therefore, our designed structure incorporating thermoelectricity broadens the applications of chiral BIC in sensors and detectors.

Research on the influence and optimization of sunshade effect on radiative cooling performance

Intense solar radiation is the main reason why radiative cooling is difficult to achieve during the daytime, as the radiative heat flux from the sun will offset the heat emitted by the cooler. Most previous studies have focused on providing coolers with high solar reflectance and high mid-infrared spectral emissivity by designing photonic structures. In this study, a macroscopic sunshade device based on the solar visual motion trajectory is proposed to reduces the direct solar radiation received by the cooler surface through setting a blocking on the solar path, while hardly affecting the radiative capacity of the cooler. Here, we built a three-dimensional sunshade device simulation model based on the surface-to-surface radiative heat transfer module in COMSOL Multiphysics® software, which can effectively improve the cooling performance of the cooler. Then, using polished aluminum plates, black paint-coated cooler, and reflective radi-cool film to test the impact of macroscopic sunshade device based on the solar apparent motion trajectory on improving the cooling performance of difference coolers. The results verified the effectiveness of the proposed sunshade devices in improving radiative cooling performance, and the results showed that the maximum temperature reduction could reach 3.35 °C, 26.96 °C, and 3.42 °C, respectively.

Geothermal Energy Extraction Using a Novel Combined Coaxial and U-Shaped Closed-Loop System

In the present study, a novel combined coaxial and U-shaped closed-loop system was proposed to enhance the thermal performance of the U-shaped geothermal system. The proposed geothermal system consists of two separate energy extraction units, called primary and secondary systems. The primary system consists of an insulated U-shaped tube heat exchanger. In order to extract more thermal energy, a coaxial double-pipe heat exchanger was coupled to the primary system as the secondary system. A conjugate heat transfer (CHT) module incorporating a k-ε turbulent model was employed to simulate the heat extraction and fluid flow within the proposed system. Further, the study extended to analyse the effect of various individual parameters on the thermal performance of the proposed closed-loop system. The results depicted that the proposed system can extract up to 42% higher thermal power than that of an insulated U-shaped tube heat exchanger. Further, the coefficient of performance (COP) of the combined system ranges between 1.03 to 1.22 depending upon the operational parameters. Overall, the proposed geothermal system effectively utilized both the heat loss from the hot primary fluid as well as the thermal energy available in the rock surrounding the insulation of U-shaped geothermal heat exchangers. The outcomes of the present study may be productive for researchers and industries interested in the extraction of clean and renewable energy.

Assessment of thermal and radiation induced creep in the dual cooled lead lithium blanket

The creep in recently designed blanket for the fusion nuclear science facility is described by decomposing it into thermal, irradiation and cavity swelling creep. Initial estimates of thermal creep using elastic high temperature rules showed that the recent heat transfer improvements at the first wall (FW) is feasible, and the FW can withstand heat flux of 0.35 MW/m2 at 3000 hrs without rupture. A viscoplastic model that is based on Norton's power law and relies on Multiphysics coupling of solid mechanics and heat transfer modules is used to capture the inelastic deformation in the blanket. Results showed that the design peak heat flux of 0.25 MW/m2 produced maximum thermal creep strain of 0.45 % and the relaxation of the thermal stress at the FW. Irradiation creep is prescribed to be proportional to the displacement damage dose and applied stress. The displacements from irradiation creep radially decrease from the FW to the back wall, but the maximum deformation is found at the back wall that is connected to the Helium manifold due to the high stresses at the region. Cavity swelling creep is dominant at steady state and combined with radiation creep to produce displacement of about 8 mm at the FW.

Energy, economic and environmental (3E) analysis of residential building walls enhanced with phase change materials

The building sector has recently been responsible for the world's energy consumption growth. Heating and cooling accounts for more than 30 % of building energy consumption. A phase change storage wall heat transfer module considering phase change material hysteresis in TRNSYS was developed in this article. A binary variable is defined to determine the enthalpy curve in heat transfer calculation, which solves the problems of overestimate the latent heat and difficult temperature convergence in the previous hysteresis model for non-complete transformation processes. An extensive study, including energy, economic and environmental analyses, has been proposed to assess the benefits of phase change materials which integrated into building walls in Xi'an. In the study of phase change storage wall building, the energy efficiency, heating and cooling loads of phase change walls with different structures were calculated. It was found that the optimal melting temperature of the phase change layer was close to the summer indoor setting temperature. EPS-Brick-PCM is the optimal phase change insulation composite exterior wall structure. A life cycle economic analysis and carbon emission study were conducted on the phase change wall. It was found that the phase change exterior wall without insulation has better economic efficiency and more carbon emission reduction. The phase change exterior wall with insulation is less economic efficient and has less carbon emission reduction. The thicker the insulation layer of the phase change insulation composite wall, the worse the economic and environmental benefits.

Computational and experimental study on the resistance welding process of a glass fiber‐reinforced epoxy‐based composite with thermoplastic interlayer adherent

AbstractIn this work, resistance welding of a glass fiber‐reinforced epoxy composite (GFRC) was studied with numerical optimization and experimental validation. A steel mesh and polymethyl methacrylate (PMMA) films were used as the heating element and adherent interlayers, respectively. A transient heat transfer module was implemented to conduct the parametric optimization study, with variables of electricity power, clamping distance and weld time. The optimal welding condition was then confirmed as 20 W, 0.4 mm and 30 s, with a melting degree of 95.2%. A thermal meter and a thermal camera validated the simulated temperature results. Welding quality was experimentally characterized by single lap shear tests and scanning electron microscopy (SEM). The highest lap shear strength of 3.8 ± 0.3 MPa was captured on the specimen welded with the optimized condition. This was 76% that of the benchmark made with the adhesive bonding method but it was over 200 times faster.Highlights Resistance welding of GFRC with PMMA films and a steel mesh is studied with FEA. Simulation results are quantitatively validated with experimental methods. Optimal welding conditions are confirmed in association with welding tests.

Numerical Simulation for Cooling of Integrated Toroidal Octagonal Inductor Using Nanofluid in a Microchannel Heat Sink

This paper presents a comprehensive numerical simulation study focused on the cooling of integrated toroidal octagonal inductor using nanofluids within a microchannel heat sink. The investigation utilizes COMSOL Multiphysics 6.0 integrated with the Fluid Flow and Conjugate Heat Transfer Module. The primary objective is to explore and understand fluid flow and heat transfer characteristics within the integrated inductor. The study involves testing three distinct fluids, water, CuO-water nanofluid, and Al2O3-water nanofluid, under laminar flow conditions within microchannels. The choice of fluid plays a significant role in heat transfer, interacting with the microchannel geometry to optimize performance. Three-dimensional computational fluid dynamics (CFD) models are meticulously developed; focusing on toroidal inductors equipped with micro pin fins heat sinks. The study commences by detailing the geometry of the micro coil and the integrated heat sink. The simulation encompasses a mathematical model that captures the intricate interplay between the governing Navier-Stokes equations for fluid dynamics and the heat transfer equations within the integrated inductor. As φ increases, temperature, viscosity, and pressure decrease. CuO-water and Al2O3-water nanofluids play a significant role in influencing laminar flow and key thermal parameters in the toroidal inductor. These nanofluids, which consist of base fluids (water) with dispersed nanoparticles (CuO or Al2O3), are employed as cooling agents to enhance heat transfer. The presence of nanoparticles in the fluid alters its thermal properties, leading to changes in the flow dynamics and overall heat dissipation within the toroidal inductor.The laminar flow characteristics are affected by the nanofluid's viscosity, density, and thermal conductivity. Additionally, the Nusselt number, Reynolds number, and thermal resistance are key thermal parameters that reflect the performance of the cooling system. The nanofluid's influence on these parameters is crucial for understanding and optimizing the thermal management of the integrated toroidal inductor. The enhancement of heat dissipation in the toroidal inductor is achieved through improved thermal properties of the nanofluid. Higher nanoparticle concentrations result in better heat transfer rates, leading to lower temperatures in the toroidal inductor. This, in turn, improves the overall efficiency and performance of the cooling system. The viscosity of the nanofluid is influenced by the presence of nanoparticles. The pressure within the microchannels is also affected by the nanoparticle concentration. An increase in φ can lead to changes in pressure drop along the microchannels. Understanding these variations is crucial for designing an effective cooling system.

Experimental and numerical analysis of an asphalt solar collector with a conductive asphalt mixture

Asphalt solar collectors (ASCs) have emerged as a promising renewable energy source. They offer a practical solution for generating renewable energy and thereby reduce the reliance on fossil fuels. In this study, a computational analysis and experimental validation of an ASC with a conductive asphalt mixture were conducted. The ASC that was investigated had 3000 mm long serpentine copper tubes positioned at the center of a 1000 mm ×60 mm ×8 mm sample of hot-mix asphalt. Ten 250 W infrared bulbs installed in two rows served as the simulated heat source of the ASC. The ASC had a water flow in the range of 14–97 L⁄h and solar intensities of 600, 800, and 1000W/m2. The thermal response of the ASC was modeled using the COMSOL Multiphysics 5.6.0.280 heat-transfer module. The experimental results indicted that the efficiency of the ASC was 53.56% or 64.1% with either reference or conductive mixtures, respectively. When a conductive mixture was used in the ASC, its surface temperature decreased by 11.37% when compared to that obtained using the reference mixture. The numerical simulation revealed that the average deviation between the numerical and experimental temperature data of the water that flowed within the ASC tubes was between 6% and 14%. The use of conductive asphalt mix would enhance the thermal performance of ASCs.

Geometric modification and placement of high heat flux ic chips on substrates of different materials for enhanced heat transfer

This article presents the significance of using different sizes and positions of IC (Integrated Circuit) chips that are mounted on substrates of different materials. The ICs are cooled by forced convection in a horizontal wind tunnel. COMSOL Multiphysics 5.4 solve the IC chip cooling problem by selecting a conjugate heat transfer module with laminar flow. FR4, ba-kelite, single and multi-layer copper clad boards are used as substrate materials. Numerical simulations are performed for FR4, bakelite with a constant heat flux of 5000 W/m2, at 2.5 m/s air velocity. In contrast, single and multi-layer copper clad boards are studied for 10000 W/m2 with 1.5 m/s air velocity. The prime objective of this research is to use adequate size and placement of chips generating high heat fluxes for enhanced heat transfer. Results showed that larger chips placed at bottom rows and sequentially decreasing sizes in the subsequent rows for the same overall input for high substrate thermal conductivity give more heat dissipation. Among all configurations A0 – F considered in the study, case E provides minimum tempera-ture using single and multi-layer copper clad boards. In addition to modification of chip sizes, geometric spacing of X1/X2 = 1.8 results in lower maximum temperature.

Computational study of a microwave plasma reactor based on the TM112 mode for diamond deposition

One of the main features of microwave plasma reactors is the electric field structure in the resonant cavity, which must be both intense and uniform in front of the substrate. For this reason, transverse magnetic modes are often used, especially axisymmetric modes because they produce an axisymmetric plasma. Microwave plasma reactors can be differentiated according to the chosen mode, because this has a direct influence on the diamond film growth process, among other features such as the coupling technique and the used quartz window. Another attractive characteristic of said reactors is obtaining large activation areas of the plasma. In this paper, we propose a microwave plasma reactor based on the TM112 cylindrical mode, which is subject to a computational study. Unlike axisymmetric modes, which activate the plasma on the cavity axis, the TM112 cylindrical mode presents two activation plasma areas. The reactor was designed following the methodology described by Silva et al., and using the Plasma, Radiofrequency (RF), and Heat transfer modules of the software COMSOL Multiphysics. The obtained results are presented in two stages. The first one is related to the initial electric field distribution of the TM112 mode. Next, the generation of the hydrogen plasma was simulated from the interaction of H2 gas with the TM112 microwave field. The plasma activation process is described in detail from graphics of the time evolution of the electron density, hydrogen density, and their respective temperatures until a steady state is reached. Additionally, the influence of the pressure on the concentration and the temperature of both electrons and gas in a steady state is analyzed. The presented results can be useful for the design of plasma reactors for diamond deposition.

Optimization of rifled tube by parametric changes using CFD for goodness factor enhancement

AbstractIn this study, extensive research was conducted using computational fluid dynamics (CFD) ANSYS®–FLUENT, with the k–ε turbulence heat transfer module to optimize the rifled‐tube heat exchanger through parametric modifications for heat transfer and Goodness factor enhancement. This is based on predicting the performance of various heat transfer surfaces such as rectangular and two opposite right‐angled triangular ribs. The objective of this research was to enhance the heat transfer efficiency of a rifled tube with rectangular ribs, a height of 0.775 mm, and helix angles of 30° and 58°. The research findings demonstrated that the heat transfer rate improved owing to the establishment of a secondary helical or swirl flow near the wall region, leading to an enhanced heat transfer. The modification of the rib geometry from rectangular to two opposite right‐angled triangles increased the heat transfer by 17%, whereas the pressure drop remained relatively constant at 0.62 mm of rib height. This benefit is directly translated into heat transfer without any change in pressure drop. The values of the friction factor (f), Nusselt number (Nu), and goodness factor are appropriate for the new geometry in relation to the reference geometry. Rifled‐tube designers can use CFD heat transfer analysis to optimize the design, reduce the need for physical prototypes, and test similar applications and rib configurations.

Thermal response of an asphalt solar collector: Numerical and experimental validation

This study conducts a numerical analysis and experimental validation of the thermal performance of an asphalt solar collector. Both parallel and serpentine copper tubes with a total length of 3000 mm were embedded in the middle of a 1000 × 600 × 80 mm hot mix asphalt model. The sides and bottom of the test sections were insulated to prevent heat loss to the surrounding environment. Using four values of water mass flow rate and three values of artificial solar radiation, the thermal performance of the system was examined. A mass flow inlet of 0.0038–0.035 kg/s is specified for the water inside the embedded tubes, while 600, 800, and 1000 W/m2 solar radiation is applied on the upper section of the asphalt mix during this study. The thermal response was simulated using the COMSOL Multiphysics 5.6.0.280 heat transfer module. Using a reference asphalt solar collector, the desired parameter effects were investigated independently. The experimental results showed that the efficiencies of the asphalt solar collector with serpentine bare tubes and parallel bare tubes are 53.56 % and 52.85 %, respectively. The numerical simulation reveals that there is no significant thermal effect in the area under the embedded tubes when compared to the area above the tubes that are exposed to solar radiation from above. The average error found between the numerical and experimental results for the water temperature difference inside the embedded tubes of the asphalt solar collector ranged between 6 % and 14 %. This good agreement indicated that the computational fluid dynamics model is reliable and capable of simulating the thermal response of an asphalt solar collector.

Heating and evaporation of a mono-component spheroidal droplet with non-uniform surface temperature

A new mathematical model for spheroidal droplet heating and evaporation is proposed. This model takes into account the effect of liquid finite thermal conductivity and is based on the previously obtained analytical solution for the vapour mass fraction at the droplet surface and a new correlation for the convective heat transfer coefficient incorporated into the numerical code. The heat transfer equation in the liquid phase is solved numerically using the finite-element heat transfer module of COMSOL Multiphysics. It is shown that the lifetime of spheroidal (prolate and oblate) droplets is shorter than that of spherical droplets of the same volume. The difference in the lifetimes of spheroidal and spherical droplets, predicted by the new model, is shown to increase with increasing aspect ratios for prolate droplets and decreasing aspect ratios for oblate droplets. As in the case of stationary spherical droplets, the d2-law is shown to be valid for spheroidal droplets after the completion of the heat-up period. The predictions of this model agree with experimental observations. The duration of the heat-up period is shown to decrease with increasing aspect ratios for prolate droplets and decreasing aspect ratios for oblate droplets. The maximal surface temperatures are predicted near the regions where the surface curvature is maximal. The aspect ratios are shown to be weak functions of time, in agreement with experimental observations.

R&D of a Hydraulic Hydrogen Compression System for Refuelling Stations

Abstract The article presents a hydraulic hydrogen compression solution designed to serve as a booster compressor. It can be adapted to changing parameters of the inlet pressure of hydrogen and allows stabilising the hydrogen accumulation process in the high-pressure storage. The main results of this study were obtained using a numerical model developed to explore the thermodynamic processes that occur during the hydraulic compression of hydrogen. The modelling was carried out using COMSOL Multiphysics® 6.0 software with the CFD and heat transfer modules. The compression chamber in the form of a cylinder with a volume of 1.14 l and wall thickness of 5 mm was used in the computational model. The aim of these simulations was to investigate the temperature change limits of hydrogen, cylinder walls and working fluid, as well as to estimate the actual value of pressure inside the cylinder. The considered process of pressure increase in the cylinder chamber was modelled as a continuous change of volume filled with working fluid with discrete time step of 0.01 s, taking into account the increase of temperature inside the cylinder. The derived modelling results for different durations of compression stroke ts from 0.5 to 20.0 s were presented. The curves of energy consumption and temperature rise during the compression process were calculated for initial hydrogen pressures P1 = 3.0, 10.0, 15.0 and 20.0 MPa and compression ratio Kc = 5.0. The results of simulation of thermodynamic processes and their analysis allowed estimating energy consumption in the system of hydraulic compression and determining conditions which would lead to the increase in efficiency of hydrogen compression operation systems under consideration.