Surface Temperature Distribution Research Articles

Study on Thermal Metamaterial for Uniform Surface Temperature Distribution

Abstract This paper proposes a thermal metamaterial based on the theory of thermodynamics transformation, which can achieve uniform temperature distribution. The required anisotropic thermal conductivity is theoretically derived. This thermal metamaterial structure is achieved by alternating aluminum and polyimide of different lengths. The temperature distribution under the condition of different convective heat transfer coefficients was simulated and compared with pure aluminum. The simulation and experimental results show that the metamaterial significantly reduced the non-uniformity of the temperature distribution.

Preparation and performance enhancement of MXene/Na2HPO4·12H2O@SiO2 phase change microcapsule

Currently, the leakage problem of phase change materials can be solved by the microencapsulation technology. Inorganic hydrated salt materials have the advantages of high latent heat density, high thermal conductivity and low cost. How to construct inorganic hydrated salt microcapsules with excellent heat storage capacity and photothermal conversion performance maintains a challenge. In this work, the MXene/Na2HPO4·12H2O@SiO2 (MXene/MPCMs) composite phase change microcapsules were designed and prepared. The Na2HPO4·12H2O was directly used as phase change material and TEOS was used as raw material to generate SiO2 wall material. The MXene played a key and important role in simultaneously enhancing their photothermal conversion efficiency while maintaining available latent heat. The prepared composite microcapsules (MPCM-1) showed more obvious heat storage and release characteristics, and the latent heat reached 169.62 J/g. Besides, when the hydrated salt was encapsulated in silica shell, the supercooling degree of hydrated salt was reduced from 14.2 °C to 10.1 °C. With the increase of MXene content, the maximum surface temperature of microcapsules significantly increased. Based on the results of surface temperature distribution, the photothermal conversion ability and heat transfer rate of microcapsules were effectively improved, which was mainly attributed to the remarkable light absorption conversion ability of MXene. When the MXene was added, the MXene/MPCMs still showed excellent heat storage properties (~ 156.64 J/g). At the same time, the MXene/MPCMs showed good chemical compatibility, crystallization characteristics and thermal stability. Considering available performance enhancement, the prepared microcapsules show good application prospect in solar energy utilization.

Imaging performance prediction of a hypersonic target in a geosynchronous orbit based on multi-dimensional information correlation of radiation, multi-spectral, and polarization

Hypersonic target detection based on infrared intensity characteristics is easily disturbed by sea surface and cloud flares when detected by space-based optical systems, which results in a low detection rate, high false alarm, and difficulty in stable detection. This paper explores a method to improve target detection performance based on the correlation of infrared radiation, multi-spectral and polarization. Firstly, the comprehensive factors that influence complex ambient illumination, atmospheric transmission, and clutter background on spectral-polarization characteristics of hypersonic targets are analyzed. Based on the global radiation scattering theory, the temperature distribution model of the hypersonic target is established by using FLUENT. The polarization emission and pBRDF model of the target is established, and the radiation polarization transfer model is generated. Secondly, the sea surface temperature distribution is obtained by inversion of Landsat8 remote sensing data. The radiation polarization transfer model of the sea surface is established based on the Cox-Munk model combined with pBRDF and the polarization emission model. Thirdly, the polarization scattering effect of atmospheric particles on the upward radiation of the interaction of the target with the sunlight is considered comprehensively, and the 6SV radiative transfer model is used to calculate the polarization effect of atmospheric particles on the upward radiation transmission of the target and the background. Then, combined with the point diffusion of the optical system and the photoelectric conversion of the detector, the multi-dimensional full-chain imaging prediction model of the hypersonic target-sea background-ambient atmosphere-optical system-detector is established. The imaging characteristics and detection performance of the target in different imaging dimensions are simulated and analyzed with the signal-to-clutter ratio (SCR). The research shows that in the direction of reflected sunlight from the sea surface, the sea surface glare is suppressed and the target is highlighted through a target detection method of multi-dimensional information. This method has better detection results than the infrared multi-spectral detection method.

Sizing horizontal metallic inclusions in insulators using lock-in inductive infrared thermography

We present a methodology to size the area and depth of horizontal metallic inclusions of unknown geometry, embedded in insulators, using lock-in inductive thermography. The method is based on averaging the amplitude and phase thermograms along circles concentric with the center of the heated region. We present an analytical calculation of the surface temperature distribution produced by a horizontal circular heat source, taking into account conductive coupling with the air, and we propose to fit the averaged data to this model. We present lock-in thermography data on resin samples containing embedded calibrated circular and rectangular Cu slabs that we excite inductively. The results indicate that it is possible to determine the area and depth of rectangular heat sources with aspect ratio below 1.5 with accuracy better than 10 %.

Preflight and Postflight Thermal/Structural Analysis of BOLT-2: The Holden Mission

BOLT-2: The Holden Mission, a follow-on to the BOLT experiment, was designed to measure the characteristics of turbulent flow in the hypersonic flight regime. This paper summarizes preflight transient thermal and structural analyses for the predicted flight trajectories that led to the selection of suitable materials for the flight experiment. Following the successful flight in March 2022, the preflight analysis methods, along with temperature measurements from the flight, were used to conduct a postflight thermal analysis based on the as-flown trajectory that attempted to match the flight data. Conservatism in the preflight design analysis was compared to the flight measurements. Databases of computed laminar and turbulent heating were generated as part of this effort that can aid in future analysis of the experimental data. The effects of modeling transitional heating during ascent were quantified when compared to a fully turbulent heating assumption. A set of estimated three-dimensional surface temperature distributions have been computed and were found to be a good match to the available flight data. Postflight structural modeling was used to estimate the differential thermal expansion at the nosetip/frustum joint interface.

Numerical Simulation of Water Film Flow and Breakup on Anti-Icing Surface

The flow and morphological characteristics of liquid water on the icing and anti-icing surfaces of aircraft are closely related to the icing characteristics and anti-icing surface temperature distribution. To predict the flow and breakup characteristics of a water film, a 3D model of continuous water film flow and a model of water film breakup into rivulets on an anti-icing surface were constructed based on the icing model, and the corresponding methods for solving the models were developed. Using the NACA0012 airfoil as a simulation object, the changing characteristics of height and velocity for a continuous water film with time and the morphological characteristics of rivulets formed from the breakup of a continuous water film were simulated numerically. The results indicate that, with an increase in inflow velocity, the time required for the water film to completely cover the surface and reach stability decreases. Downstream in the water droplet impact zone, the calculated values of continuous water film height align well with experiments, as well as the stream height at the continuous water film rupture location with the experimental values. With the reasonable contact angle, the calculation error of the stream width is within 10%.

Geographic analysis of the distribution of surface temperature, vegetation cover and its relationship with socioeconomic indicators – Cuiabá/MT

In Brazil and around the world, the problems resulting from urban planning, changes in land cover and the increase in urban infrastructure influence temperature change. The aim of this research is to understand the difference in temperature over the last 36 years (1985-2021) in the city of Cuiabá/MT, and its relationship with vegetation and socioeconomic and racial-ethnic factors. Annual images (36 images) from the Landsat 5 and 8 satellites between 1985 and 2021 were used, and Google Earth Engine was used to obtain this data for the urban area of the city of Cuiabá. Mapping of the spatial distribution of surface temperature, vegetation and changes in land use and cover was carried out using data from the MapBiomas project and socio-economic data from the Instituto Brasileiro de Geografia e Estatística (IBGE, Brazilian Institute of Geography and Statistics in English) 2010 census. The results indicated that over the last 36 years Cuiabá's surface temperature has increased by an average of 13°C, as part of the result of the growth of urban infrastructure and loss of vegetation cover, specifically in the peripheral areas. There is an inequality between the high-income population living in neighborhoods with milder temperatures and the low-income population living in neighborhoods with higher temperatures. In the temperature and skin color data, the white population generally lives in neighborhoods with milder temperatures, unlike the black population, which mostly suffers from high temperatures.

Effect of the Backboard on the Pyrolysis Front Morphology of Downward Flame Spread: A Simulation Study

In order to investigate the influence of the backboard on the morphology of the pyrolysis front, this study employed the fire dynamics simulator (FDS) to conduct simulation research, aiming to quantitatively analyze the reasons for the differences in the morphology of the pyrolysis front under two conditions: with and without the backboard. By comparing the surface temperature distribution, flow field, and air entrainment velocity at the bottom of the pyrolysis front made of polymethyl methacrylate (PMMA) under different conditions, it was found that the presence of the backboard not only reduced the surface temperature of PMMA at the bottom of the pyrolysis front and the temperature difference between the edge and middle, but also changed the flow field at the edge and suppressed the horizontal air entrainment velocity at that position, especially for thinner PMMA. These factors ultimately result in a larger angle of the pyrolysis front during the burning stable stage under the condition with the backboard compared to the condition without the backboard. In addition, the trend of horizontal air entrainment velocity at the edge position under different conditions also provide data support for the qualitative observation and analysis of previous experiments phenomena.

A machine learning-driven approach to predicting thermo-elasto-hydrodynamic lubrication in journal bearings

Traditional methods of evaluating the performance of journal bearings, for example thermal-elastic-hydrodynamic- lubrication theory, are limited to simplified conditions that often fail to accurately model real-world components. Numerical models that include additional phenomena such as cavitation and fully coupled effects like deformation, temperature, pressure and viscosity can be more accurate but require a large amount of computational overhead, making analysis slower and more costly. To address this limitation, a novel machine learning-driven approach is developed to predict the 2D distribution of surface deformation, film thickness, temperature, and pressure across the bearing surface as a function of design variables such as load and speed. The training dataset, generated using a fully coupled Reynolds’ Equation solver implemented in OpenFOAM, contains a significantly extended range of conditions than in previous studies with approximately 39000000 points encompassing 4925 different test cases. Modelled bearing speeds range from 2000 to 10000rpm, while load values are varied between 1 and 30kN. Predicting surface deformation, film thickness, temperature and pressure across the bearing surface results in a mean absolute percentage error below 0.4% or better. The work also demonstrates that the trained models have a strong ability to generalise the prediction beyond the original training data range with only a 1% error at up to 200% of the highest trained speed. This work also demonstrates that machine learning-based processes are a practical alternative to physics-based numerical modelling, especially in cases where rapid performance evaluation is desired as real-time calculation is possible with significantly reduced computational cost. This has the potential to enable development of rapid design optimisation tools and real-time performance monitoring at high resolution and with low latency. Using consumer hardware, it is found that the neural network-based approach is faster than the existing numerical modelling technique by a factor of over 10000, enabling real-time predictions of lubrication systems.

Numerical simulation of an industrial radiant tube burner using OpenFOAM

This study presents a numerical investigation into the effects of physical models on the prediction accuracy of the wall temperature distribution in an industrial radiant tube burner. Utilizing a reacting flow solver based on OpenFOAM, we explored the effects of various physical models, including those for chemistry, combustion, heat transfer, and radiation properties. The choice of combustion model significantly influences prediction accuracy, playing a more dominant role than the chemistry mechanism. Moreover, the simulations captured a distinctive triple flame structure inside the burner, representing the coexistence of rich premixed, non-premixed, and lean premixed flame structures. Conditional scatter plots displayed the development of both premixed and non-premixed flame structures, converging on the fuel-lean side. Notably, accurate prediction of wall temperature distribution depends on the incorporation of a precise heat transfer model, coupled with a detailed radiation property model. Regarding the distribution of tube surface temperature in the main radiation zone (a distance from the burner nozzle greater than 1 m), the most accurate prediction exhibits a maximum deviation of less than 56 K and an average deviation of 24 K compared to experimental results. The simulation closely matched experimental data for exhaust concentration of NO within an error margin of 20 ppm. However, discrepancy was observed in the CO concentration, which was attributed to the simplified representations of fuel chemistry and composition, as well as the difficulties in accurately capturing the unsteady flame dynamics near the wall.

Heat-transfer characteristics of the two-phase-change process between a microencapsulated phase-change-material-slurry droplet and air: A numerical study

A concept is proposed for spray-heat exchangers with the aim of addressing the mismatch of heat supply-demand within and outside a water droplet. The concept involves using a microencapsulated phase-change-material (MPCM) slurry as the working medium, enabling two-phase-change heat transfer with air. The temperature of the water droplet rapidly decreases owing to continuous evaporation, leading to a rapid decrease in its heat-transfer capacity. The phase change of the MPCM particles not only provides sufficient heat for evaporation but also maintains a high surface temperature, thus sustaining effective evaporation. A numerical model was developed to theoretically investigate the effect of the MPCM phase change on external convection. The heat-mass transfer of a droplet with a radius of 0.25 mm composed of water and 20 vol% MPCM was simulated. The results indicate that the MPCM phase change effectively reduces the temperature decrease rate. However, when the overall surface temperature of the MPCM droplet exceeds that of water, the total heat-transfer coefficient is not necessarily higher. Further local analysis reveals that maintaining both an overall high surface temperature and a reasonable local temperature distribution is crucial to achieve a higher total heat-transfer coefficient. The MPCM phase change significantly influences the local surface-temperature distribution.

Enhancing Parabolic Trough Solar Collector Efficiency through Absorber Tube Rotation: A Computational Study

This work uses CFD simulations to examine how absorber tube rotation affects parabolic trough solar collector performance, including output temperature, thermal efficiency, and surface temperature distribution. Simulation results are closely matched to experimental data, with variances ranging from 0.54% to 4.63%. It examines how mass flow rates (0.01 kg/s, 0.015 kg/s, and 0.02 kg/s) and input temperatures (32°C to 38°C) affect system performance in stationary and rotational modes at 1 and 2 rpm. Tube rotation considerably improves thermal mixing in the heat transfer fluid (HTF), resulting in more equal temperature distribution and possibly improved heat transfer efficiency. Increased rotational speed causes complicated dynamics that slightly impair thermal efficiency due to centrifugal forces affecting fluid flow patterns. The study shows that thermal efficiency peaks in the morning, emphasizing the importance of climatic and operational factors on collector performance. DOI: https://doi.org/10.52783/pst.276

Exploring the body surface temperature of the lumbosacral region and relevant back-shu points in patients with lumbar disc herniation induced low back pain based on infrared thermography.

To observe the body surface temperature of the lumbosacral region and relevant back-shu points in patients with lumbar disc herniation (LDH) induced low back pain utilizing infrared thermography, and to explore the functional attribute changes of acupoints under pathological conditions. A total of 50 patients with LDH induced low back pain were included as the observation group, and 45 healthy subjects were included as the control group. Using infrared thermography, the body surface temperature of the lumbosacral region and bilateral Sanjiaoshu (BL 22), Shenshu (BL 23), Qihaishu (BL 24), Dachangshu (BL 25), Guanyuanshu (BL 26), Xiaochangshu (BL 27), and Pangguangshu (BL 28) was measured in both groups. The temperature difference values between the bilateral lumbosacral regions and back-shu points of the two groups were calculated. Additionally, the body surface temperature of the affected and healthy sides of the lumbosacral region and relevant back-shu points was compared in the observation group. Compared with the control group, the body surface temperature of the lumbosacral region and the bilateral temperature difference values of the lumbosacral regions were increased in the observation group (P<0.001). The body surface temperature difference values of bilateral Shenshu (BL 23), Qihaishu (BL 24), Dachangshu (BL 25), Guanyuanshu (BL 26) and Xiaochangshu (BL 27) in the observation group were higher than those in the control group (P<0.05, P<0.01, P<0.001). In the observation group, the body surface temperature of the affected side of the lumbosacral region as well as Shenshu (BL 23) and Dachangshu (BL 25) was elevated compared with that of healthy side (P<0.001). The patients with LDH induced low back pain have imbalanced and asymmetrical distribution of body surface temperature in the lumbosacral region and related back-shu points, Shenshu (BL 23) and Dachangshu (BL 25) have the relative specificity.

Flow structure and cooling effect of seeping gas film in supersonic turbulent boundary layer

Seeping gas film cooling is an active cooling method that involves injecting a cooling gas film into the boundary layer through porous media. In this study, a sintered wire mesh porous material is used to generate the gas film through seepage flow. The flow of the gas film in the supersonic turbulent boundary layer is observed using nanometer particle plane laser scattering technology, while the temperature distribution of the porous wall surface is measured by infrared thermometry. It is observed that as the injection rate of the gas film exceeds 0.2 %, the supersonic turbulent boundary layer thickens significantly. The rate of thickness growth is promoted from 5.5 % to 17.1 % as the injection rate increases from 0.3 % to 0.8 %. Additionally, the injection of the gas film induces an oblique shock wave, and the shock angle increases with the enhancement of the injection rate, resulting in an uneven cooling effect along the injector length. The transient cooling efficiency of the seeping gas film under thermal non-equilibrium conditions is analyzed. The results indicate that the temperature amplitude on the porous wall surface decreases significantly when the injection rate is promoted from 0.2 % to 0.3 %. However, continuing to increase the coolant injection rate does not lead to a more significant improvement in cooling efficiency. When the injection rate reaches 0.5 %, the relative cooling efficiency reaches 50 %–60 %, which is much lower than the predicted values of the thin film theory.

Heat transfer measurements on a flat plate and cylinders at high velocity conditions: Comparison of calculation models with infrared thermography

In this study, various heat transfer measurement models using infrared (IR) thermography were applied to a flat plate and cylinders to examine the accuracy and suitability of the measurements in a transonic blow-down test facility recently developed at Korea Aerospace University. Experiments were performed under three inlet velocity conditions for the flat plate, and for the cylinders, two diameters, two turbulence intensity levels, and five inlet velocity conditions. The surface temperature distribution during the experiment was acquired using an IR camera, and the surface heat transfer coefficient (HTC) was calculated using four HTC calculation models based on the one-dimensional semi-infinite (1D SI) solid assumption: step change (STEP), Duhamel's superposition (DUH), 1D finite volumetric model (1D FVM), and linear regression model (LRM). The LRM showed minimum deviation with reference correlation for the flat plate heat transfer measurements, while STEP, DUH, and 1D FVM showed increased deviations at higher mainstream velocity conditions. The measurements by LRM at the cylinder stagnation location showed the most comparable value with the reference correlation, whereas STEP and FVM showed about 10% deviation from the reference. It is concluded that the LRM is the most suitable heat transfer calculation model for high velocity tests because the model excludes the ramping data during the transient tests, and it can also obtain the local adiabatic wall temperature experimentally.

Seasonal variations in chlorophyll–a and sea surface temperature in the exclusive economic zone of Sri Lanka

Oceanic phytoplankton plays a crucial role in regulating the biogeochemical cycle in marine ecosystems. Light harvesting pigments such as chlorophyll–a (chl–a) can be used to understand marine primary productivity, especially phytoplankton abundance. Sea surface temperature is closely linked with phytoplankton abundance and climatic changes. The objective of the current study is to determine spatial and temporal variations of chl–a and sea surface temperature and their co-variation, over 12 years using statistical methodologies and Remote Sensing, Geographic Information System (GIS) and Cloud Computing Approach. In this study, the available sea surface temperature and chl–a data acquired from a Moderate Resolution Imaging Spectroradiometer (MODIS) satellite sensor were proceeded through Google Earth Engine and Earth Science Data System. Data visualization was carried out using ArcMap™ software. Results show a weak negative linear relationship (−0.40, P < 0.05) between sea surface temperature and chl–a. Therefore, additional favorable environmental factors such as coastal upwelling, and nutrient discharge through river run-off can control the primary productivity of phytoplankton. We also illustrated spatial and temporal maps for chl–a and sea surface temperature distribution in the exclusive economic zone of Sri Lanka. Chlorophyll–a concentration is enhanced during the southwest monsoon, and the lowest concentration is observed during the first inter-monsoon period. Next, the highest sea surface temperature is recorded during the first inter-monsoon period, and the lowest sea surface temperature is observed during the northeast monsoon. Seasonal behavior of chl–a concentration and sea surface temperature are changed as chl–a (SWM) > chl–a (NEM) > chl–a (SIM) > chl–a (FIM) and SST (FIM) > SST (SIM) > SST (SWM) > SST (NEM), respectively. In addition, both trend lines for each variable show a positive gradient, which may be attributed to global warming. Finally, we observed the alteration of chl–a and sea surface temperature that linked to changes in ocean chemistry during the global COVID–19 pandemic.

Crossflow Instability on a Swept Fin–Cone with Variable Nose Bluntness in Mach 6 Flow

Experiments and computational analysis have been conducted to study instabilities in the boundary layer over a swept fin–cone geometry in hypersonic flow. Experiments were carried out at the United States Air Force Academy’s Mach 6 Ludwieg Tube. Infrared thermography was employed to assess the surface temperature distribution, particularly the thermal striations on the fin that indicated a crossflow-dominated transition. Increasing the unit Reynolds number increased heating and moved the crossflow-transition front closer to the fin leading edge. Stability analysis based on the linear parabolized stability equations identified unstable stationary crossflow modes across the range of tested unit Reynolds numbers, and the most amplified disturbances were found to have wavelengths of 2–4 mm. The experimentally observed thermal striations were analyzed and found to have initial wavelengths of 2–5 mm but eventually formed higher-amplitude structures with a wavelength of approximately 10 mm. Discrete roughness elements of different diameters and spacings were implemented near the fin leading edge. Boundary-layer transition onset was found to be delayed with a subsequent lowering of the fin surface heating, with discrete roughness elements having a wavelength of 13.33 mm.

Investigations on Hot Air Anti-Icing Characteristics with Internal Jet-Induced Swirling Flow

The tangential jet-induced swirling flow is a highly efficient technology for enhancing heat transfer. This paper explores the application of swirling flow of an airfoil/aero-engine in a hot air anti-icing chamber, aiming to improve the anti-icing performance and achieve a more uniform temperature on the surface. A series of numerical computations adopting the SST k − ω turbulent model was carried out to obtain the internal flow and heat transfer characteristics, as well as the surface temperature distributions, considering water evaporation and solid heat conduction. Three jet arrangements, including impingement jets, offset jets, and swirl jets, were studied and compared, which evidently showed that the swirling effect was helpful to elevate the internal heat transfer. Compared to the impingement jets at the Reynolds number of 40,000, the Nusselt number with the offset jets is increased by 19.5%, while the corresponding Nusselt number of the swirl jets is augmented by 44.3%. The swirling flow significantly elevates the swirl number within the internal chamber, intensifying the vortex strength near the wall and increasing the circumferential velocity, which also results in an augmentation of internal pressure loss. By adopting the swirling internal flow, the temperature distribution on the anti-icing surface is more uniform and is increased by up to about 4.1 K in the leading edge when the internal-to-external temperature difference is 80 K. Simultaneously, the heat absorption of water evaporation and the matches between the internal heat transfer and external icing load are of particular importance to determine the anti-icing performance, and this has been discussed in this paper.

Intra-Operative Infrared Thermography of Ocular Surface During Cryodestruction of Epibulbar Tumours

The paper presents the possibilities of intraoperative use of infrared (IR) thermography to image the temperature distribution of the outer ocular surface during cryodestruction of conjunctival tumours with localization in the ciliary body projection to determine an individual freezing regimen. IR thermography was carried out using a portable device for infrared real-time thermography. Thermal fields were visualized, and the duration of cooling and warming of the structures of the outer ocular surface around the area of freezing, including the cornea, were recorded. It was found that intraoperative IR thermography during cryodestruction of epibulbar tumours in the projection of the ciliary body allows real-time noninvasive imaging of the temperature fi eld dynamics around the freezing zone, which can be used to select individual cooling exposure to reduce the risk of complications by avoiding excessive cooling of intraocular structures. It was also determined that IR thermography makes it possible to assess the individual dynamics of rewarming of the ocular tissues to the initial level if a repeated cycle of tumour cryodestruction is necessary.

A bibliometric and performance evaluation of nano-PCM-integrated photovoltaic panels: Energy, exergy, environmental and sustainability perspectives

One of the major problems regarding PV panels is the decline in power output and efficiency whilst exposed to temperatures surpassing their operating temperature. In order to preclude such undesirable situation, it is imperative to cool PV panels and provide a uniform distribution of surface temperatures during the implementation of the cooling method. Thermal management can be achieved at the surface temperatures of PV panels by utilizing phase change materials (PCMs). In this study, along with PCM, the potential of enhancing output parameters by decreasing the surface temperature of PV panels with the addition of nanoparticles (Al2O3) at different concentrations (0.05%, 0.1%, and 0.15% w/v) to PCM (RT35) is examined. The study compared five systems: a reference PV panel (PV), PV panel cooled with PCM without nanoparticles (PVPCM-0), and PV panels with PCM containing different concentrations of nanoparticles (PVPCM-0.05, PVPCM-0.1, and PVPCM-0.15). Among the five different systems, the PV panel containing 0.15% w/v nanoparticles (referred to as PVPCM-0.15) demonstrated the most effective cooling capability. Moreover, the PVPCM-0.15 system provided the highest performance with a 19.49% increase in panel power output. PV systems have average energy and exergy efficiency values of 9.06% and 3.79% for PV panel, 9.60% and 5.15% for PVPCM-0, 9.70% and 5.12% for PVPCM-0.05, 10.28% and 6.01% for PVPCM-0.1, and 10.44% and 7.29% for PVPCM-0.15. Upon analyzing the sustainability metrics, it was determined that the PVPCM-0.15 system was more energy and environmentally sustainable than the others.