International Journal of Thermal Sciences最新文献

{"title":"Multi-objective optimization of double-layer microchannel heat sink with trapezoidal cross-sections based on computational fluid dynamics","authors":"Rafat Mohammadi,&nbsp;Vahid Dadras","doi":"10.1016/j.ijthermalsci.2025.109879","DOIUrl":"10.1016/j.ijthermalsci.2025.109879","url":null,"abstract":"<div><div>This study presents a multi-objective optimization of a trapezoidal double-layer microchannel heat sink (TDL-MCHS) to minimize both thermal resistance and pumping power. A three-dimensional fluid-solid conjugate heat transfer model, coupled with a multi-objective genetic algorithm, was employed. The optimization considered five design variables, including the flow rate ratio between the upper and lower channels and four geometric parameters related to the channel cross-sections. A parametric study explored the design space, and response surface approximation was applied to improve computational efficiency. The results showed that the optimized TDL-MCHS achieved up to a 42 % reduction in thermal resistance, though at the cost of a significant increase in pumping power. Conversely, minimizing pumping power by 42 % led to a 4 % reduction in thermal resistance. The Pareto-optimal front highlights the trade-offs between thermal performance and energy consumption, providing valuable insights for the efficient design of TDL-MCHSs in electronic cooling applications.</div></div>","PeriodicalId":341,"journal":{"name":"International Journal of Thermal Sciences","volume":"214 ","pages":"Article 109879"},"PeriodicalIF":4.9,"publicationDate":"2025-03-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143645126","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Study on solar radiation transfer model of double skin façade with spray aerosol","authors":"Yanjin Wang,&nbsp;Fangfang Wang,&nbsp;Jintao Xiong","doi":"10.1016/j.ijthermalsci.2025.109875","DOIUrl":"10.1016/j.ijthermalsci.2025.109875","url":null,"abstract":"<div><div>Spraying droplets into the cavity of the double skin façade can improve its thermal performance. However, the spray system is opened, the mist droplets (aerosols) dispersed in the cavity alter the transmission mechanism of solar radiation through the double skin façade. This study develops a solar radiation transfer model for double skin facade with spray aerosol based on the Mie scattering theory and the net radiation method. The model calculates the transmittance, reflectance, and absorptance of the double skin façade with spray aerosol. Experimental validation shows good agreement with the model, with a maximum error of approximately 11.2 % for solar heat gain. Additionally, this study examines key factors that influence the optical properties of the double skin facade with spray aerosol, including aerosol particle number concentration, average radius, cavity distance, and incidence angle. The results indicate that total transmittance decreases as aerosol concentration and average radius increase. However, when the concentration exceeds 600 particles/cm³ and the average radius exceeds 15 μm, the reduction in transmittance becomes less pronounced. Changes in cavity distance and incidence angle have a minimal effect on transmittance at high aerosol concentrations. By controlling aerosol concentration and average radius, solar heat gain can be effectively reduced. The model accurately describes solar radiation transmission in real conditions, helping assess the optical and thermal properties of double skin façades with spray aerosol.</div></div>","PeriodicalId":341,"journal":{"name":"International Journal of Thermal Sciences","volume":"214 ","pages":"Article 109875"},"PeriodicalIF":4.9,"publicationDate":"2025-03-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143637794","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Numerical study of flow and heat transfer characteristics of the novel small-channel thermal protection component","authors":"Wanxiang Yao ,&nbsp;Xudong Zhang ,&nbsp;Tianqi Shao ,&nbsp;Yixuan Zhang ,&nbsp;Puyan Xu ,&nbsp;Yifan Li ,&nbsp;Man Fan ,&nbsp;Feng Shi ,&nbsp;Weixue Cao ,&nbsp;Bin Yang","doi":"10.1016/j.ijthermalsci.2025.109880","DOIUrl":"10.1016/j.ijthermalsci.2025.109880","url":null,"abstract":"<div><div>The thermal protection of high temperature as well as large heat flux surface was a key technological challenge in the development of thermal protection science. In this paper, a novel small-channel thermal protection component for curved surface cooling was proposed with reference to tree branching laws. The internal flow and heat transfer characteristics were investigated using numerical simulation techniques. Firstly, a variety of operating conditions were designed to explore the correlation between different material properties and thermal protection effect. Secondly, the heat flow state inside the thermal protection component was analyzed to obtain the optimal operating conditions. Finally, the results showed that the critical Reynolds number for flow within the thermal protection component was 4k and the coefficient of local resistance was 1.67. The thermal insulation coefficient and thermal resistance of the component were 81.12 % and 2.76E-4 at different operating conditions, respectively. The pressure difference between the different stages of flow channels were kept at a steady value during the boiling phase transition heat, respectively. This research was important for the development of electronic communication microelectronics, aerospace and solar energy.</div></div>","PeriodicalId":341,"journal":{"name":"International Journal of Thermal Sciences","volume":"214 ","pages":"Article 109880"},"PeriodicalIF":4.9,"publicationDate":"2025-03-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143637751","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A modeling method for the radiative characteristic parameters of a composite medium containing base fluid and randomly dispersed nanoparticles","authors":"Li Jiayu,&nbsp;Rong Teng","doi":"10.1016/j.ijthermalsci.2025.109862","DOIUrl":"10.1016/j.ijthermalsci.2025.109862","url":null,"abstract":"<div><div>Radiative heat transfer in a composite medium containing randomly dispersed nanoparticles exists widely in nature and industrial applications. The prediction of radiative characteristic parameters is a crucial issue for the simulation of radiative heat transfer in particulate composite media. The morphology and distribution of nanoparticles can affect the interaction between electromagnetic radiation and the nanoparticles, thereby influencing the radiative characteristic parameters of the composite medium. To solve this problem, an electromagnetic model is constructed for a composite medium containing base fluid and randomly distributed nanoparticles. The computational domain is divided into several nanoscale cubic grid cells. Then, the effective radiative characteristic parameters of a grid cell are simulated using finite-element method (FEM), incorporating the dependent scattering effects from nanoparticles in adjacent grid cells. FEM scattering models are established based on varying degrees of interparticle interaction, and the influence of these interaction degrees on the effective radiative characteristic parameters is analyzed. The multigrid Monte Carlo (MC) program is used to simulate the radiative transfer with the inputs of effective radiative characteristic parameters. Finally, the absorptivity of the composite medium containing base fluid and randomly dispersed nanoparticles is obtained. The simulation results presented in this study indicate that the influence of dependent scattering on the radiative characteristic parameters of a particulate composite medium increases with an increasing nanoparticle volume fraction (<span><math><mrow><msub><mi>f</mi><mi>v</mi></msub></mrow></math></span>). The absorptivity of the composite medium does not definitely increase with increasing <span><math><mrow><msub><mi>f</mi><mi>v</mi></msub></mrow></math></span>. The established method can be used to analyze the influences of morphology and the distribution of the nanoparticles in a particulate composite medium. Both the dependent scattering of nanoparticles and the interactions between nanoparticles and the base fluid are taken into account.</div></div>","PeriodicalId":341,"journal":{"name":"International Journal of Thermal Sciences","volume":"214 ","pages":"Article 109862"},"PeriodicalIF":4.9,"publicationDate":"2025-03-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143637753","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Accurate estimation of interfacial thermal conductance between silicon and diamond enabled by a machine learning interatomic potential","authors":"Ali Rajabpour ,&nbsp;Bohayra Mortazavi ,&nbsp;Pedram Mirchi ,&nbsp;Julien El Hajj ,&nbsp;Yangyu Guo ,&nbsp;Xiaoying Zhuang ,&nbsp;Samy Merabia","doi":"10.1016/j.ijthermalsci.2025.109876","DOIUrl":"10.1016/j.ijthermalsci.2025.109876","url":null,"abstract":"<div><div>Thermal management at silicon-diamond interface is critical for advancing high-performance electronic and optoelectronic devices. In this study, we calculate the interfacial thermal conductance between silicon and diamond using a computationally efficient machine learning (ML) interatomic potential trained on density functional theory (DFT) data. Using non-equilibrium molecular dynamics (NEMD) simulations, we compute the interfacial thermal conductance (ITC) for various system sizes. Our results reveal an extremely close agreement with experimental data than those obtained using traditional semi-empirical potentials such as Tersoff and Brenner which overestimate ITC. In addition, we analyze the frequency-dependent heat transfer spectrum, providing insights into the contributions of different phonon modes to the interfacial thermal conductance. The ML potential accurately captures the phonon dispersion relations and lifetimes, in good agreement with DFT calculations and experimental observations. It is shown that the Tersoff potential predicts higher phonon group velocities and phonon lifetimes compared to the DFT results. Furthermore, it predicts higher interfacial bonding strength, which is consistent with higher interfacial thermal conductance as compared to the ML potential. This study highlights the use of ML interatomic potentials to improve the accuracy and computational efficiency of thermal transport simulations of complex material interface systems.</div></div>","PeriodicalId":341,"journal":{"name":"International Journal of Thermal Sciences","volume":"214 ","pages":"Article 109876"},"PeriodicalIF":4.9,"publicationDate":"2025-03-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143645125","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Numerical study on the liquid phase structural evolution of high-temperature metal-oxide mixture in the magnetic field","authors":"Kerong Guo ,&nbsp;Houjun Gong ,&nbsp;Yang Li ,&nbsp;Yuanfeng Zan ,&nbsp;Zumao Yang ,&nbsp;Wenbin Zhuo","doi":"10.1016/j.ijthermalsci.2025.109866","DOIUrl":"10.1016/j.ijthermalsci.2025.109866","url":null,"abstract":"<div><div>In reactor safety analysis, the stratification phenomenon of the molten pool is crucial for the design of in-vessel retention techniques. During the course of the experimental study on the stratification of the molten pool using an electromagnetic cold crucible, the electromagnetic field affects the evolution of the structural morphology of the molten pool. This study constructs a multi-physics field model coupling electromagnetic field, flow field, temperature field, and two-phase flow to investigate the morphological structure, heat transfer, and fluid dynamics of immiscible two-phase liquids in the electromagnetic field. The model focuses on the effects of Lorentz force, buoyancy, surface tension, temperature gradients, and solidification on the two-phase liquid structure. The simulated result of the liquid-phases’ structure aligns well with experimental results. The computational results show that when subjected to an electromagnetic field, the metal surface undergoes a significant Lorentz force owing to the skin effect. Therefore, the metal was pushed towards the center of the molten pool. The buoyancy force causes the metal to reside above the molten pool. And under the combined effects of the Lorentz force and surface tension, the metal adopts a semi-spherical shape. In the absence of the Lorentz force, the buoyancy force predominates over the interaction forces between the two liquid phases, causing the metal to spread over the molten pool. In addition, natural convection due to temperature gradients affects the molten pool flow. During the solidification of the molten pool, the solidification of the oxide on the sidewalls restricts the flow and morphology of the metal. The study finds that the molten pool is primarily influenced by Lorentz force, followed by buoyancy force and natural convection, while surface tension has the least impact on the molten pool's morphology. These findings contribute to the understanding of the complex morphological evolution process of immiscible liquid phases with different conductivities in an electromagnetic field.</div></div>","PeriodicalId":341,"journal":{"name":"International Journal of Thermal Sciences","volume":"214 ","pages":"Article 109866"},"PeriodicalIF":4.9,"publicationDate":"2025-03-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143636435","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Icing characteristics of supercooled sessile water droplets on the top of cold micro-pillars","authors":"Ruoxiao Huang,&nbsp;Xuan Zhang,&nbsp;Shuang Zhao,&nbsp;Yubo Gao,&nbsp;Long Zhang,&nbsp;Mengjie Song","doi":"10.1016/j.ijthermalsci.2025.109871","DOIUrl":"10.1016/j.ijthermalsci.2025.109871","url":null,"abstract":"<div><div>Icing and frosting problems on cold surfaces affect the normal operation of equipment and optimizing the anti-icing and ice-phobic properties of structured surfaces needs exploration of the droplet icing process on typical micro-pillars. Based on the apparent heat capacity method, the icing characteristics of sessile water droplets on the top of cold micro-pillars are numerically studied with the supercooling degree considered. The effects of the micro-pillar diameter and height as well as the droplet volume and surface temperature are obtained. As the micro-pillar diameter becomes smaller, the icing rate of the droplet decreases and the freezing time increases. A higher micro-pillar enlarges the thermal resistance, slows down the movement of the freezing front, and results in an increase in the freezing time. The freezing time goes up as the droplet volume and the surface temperature increase. This changing trend becomes more conspicuous for a smaller micro-pillar diameter. Furthermore, the relationship between the freezing time and the micro-pillar diameter and height is derived from heat transfer analysis. The freezing time is negatively related to the square of the micro-pillar diameter. When the micro-pillar height increases one time, the droplet freezing time will increase by 3.42 %. The findings in this work give insights into the icing mechanism of supercooled sessile water droplets on the top of cold micro-pillars and provide references for the design and optimization of anti-icing and anti-frosting surfaces.</div></div>","PeriodicalId":341,"journal":{"name":"International Journal of Thermal Sciences","volume":"214 ","pages":"Article 109871"},"PeriodicalIF":4.9,"publicationDate":"2025-03-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143637279","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Numerical and experimental evaluation of temperature field and melt flow in keyhole laser welding of dissimilar duplex stainless steel and nickel base alloy","authors":"Xuefeng Li ,&nbsp;Awatif M.A. Elsiddieg ,&nbsp;Aisha M. Alqahtani ,&nbsp;Mohamed Ben Ammar ,&nbsp;Ali Alzahrani ,&nbsp;Mohamed Hussien ,&nbsp;Saipunidzam Mahamad","doi":"10.1016/j.ijthermalsci.2025.109858","DOIUrl":"10.1016/j.ijthermalsci.2025.109858","url":null,"abstract":"<div><div>To achieve high quality joint in keyhole laser welding of two dissimilar metals, phase transition behavior, the temperature and velocity field according to the variation of the process parameters were evaluated by utilizing both experimental and numerical approach. Due to the existing complex phenomena, the comprehensive analysis of the weld geometry and temperature field dependency in keyhole formation was performed either numerically or experimentally. An accurate numerical simulation of temperature and velocity fields, as well as material phase change at circular geometry path of laser beam movement were analyzed on dissimilar metals of duplex 2205 stainless steel and AISI 685 alloy metals to estimate such mentioned phenomena that could not be merely evaluated via experiments. A multi-physics numerical model that employed the finite volume method (FVM) and volume of fluid method (VOF) was utilized. The major novelty of dissimilar circular weld joint was simultaneous estimation the effect of different size and thereby volume of AISI 685 alloy and duplex 2205 alloy on the parts heat sink capacity, temperature gradient, melting ratio, fusion zone microstructure and fusion zone melt volume. The main reason for this is the asymmetric temperature distribution, resulting from the combined effects of material properties and the differing geometries and material volumes of the welded parts. To distinguish the laser process parameters, impact on the weld characterization according to the numerical simulation, the findings demonstrated that increasing the speed of the laser beam leads to the formation of bulge on the part's surface and around the keyhole while simultaneously diminishing the vapor volume. Furthermore, the laser beam's deviation from −0.25 mm at the AISI 685 alloy sheet to +0.25 at duplex 2205 led to the temperature reduction up to 300 °C at 1 mm distance from the joint centerline. Comparing the weld bead geometrical changes according to the variation of laser power and welding speed implies that the predicted temperature field of numerical simulation results is in good agreement with experimental results of weld bead geometry. The maximum error for experimental temperature measurement according to the variation of welding speed and laser power was less than 3 percent. By increasing laser power from 300 to 400 W, not only has the weld bead width become twofold, but also it penetrated toward the thickness completely, and the amount of weld bead overlap evidently increased more than 40 percent. The dissimilar joint fusion zone is mainly composed of cellular and columnar dendrite microstructure mainly created from nickel base alloy solidification according to the rapid heating followed by fast cooling induced by laser heating during welding.</div></div>","PeriodicalId":341,"journal":{"name":"International Journal of Thermal Sciences","volume":"214 ","pages":"Article 109858"},"PeriodicalIF":4.9,"publicationDate":"2025-03-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143637605","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Rotational flow and heat transfer in a serpentine cooling channel with realistic internal cooling schemes of a turbine blade","authors":"Jie Wen ,&nbsp;Chenghua Zhu ,&nbsp;Yanan Chen ,&nbsp;Guoqiang Xu ,&nbsp;Hao Li ,&nbsp;Jiale Wang","doi":"10.1016/j.ijthermalsci.2025.109863","DOIUrl":"10.1016/j.ijthermalsci.2025.109863","url":null,"abstract":"<div><div>Modern advanced turbine blade mid-chord cooling systems typically have three passages with different geometric shapes and cooling schemes. The current study conducts experimental and numerical analysis of the aerothermodynamic performance in a blade-shaped serpentine channel. The channel features asymmetric cross sections, 180-degree tip and hub turns, a minor secondary inlet, staggered ribs and bleed holes. The main inlet Reynolds number (Re) and rotation number (Ro) respectively vary between 17000 and 33000 and from 0 to 0.4, and the mass flow ratio of the minor secondary coolant to the main (MR) ranges from 0 to 0.2. It is revealed that the flow interactions between bleed holes and ribs significantly improve wall heat transfer. The rotation effect on heat transfer is less pronounced in a realistic channel than in a smooth one. The minor secondary stream can increase the channel heat transfer, and the ideal MR falls between 0.1 and 0.15. The proportion of the mass flow rate of each bleed hole to the total remains almost consistent regardless of the Re and Ro. Finally, the correlations of averaged heat transfer with high accuracy (≤10 %) are developed, which could interest turbine blade researchers and designers.</div></div>","PeriodicalId":341,"journal":{"name":"International Journal of Thermal Sciences","volume":"214 ","pages":"Article 109863"},"PeriodicalIF":4.9,"publicationDate":"2025-03-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143637604","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Numerical research of a new pipe network cooling scheme without film holes for the gas turbine blade mid-chord region","authors":"Yu Sun,&nbsp;Xiaojun Fan,&nbsp;Jiao Wang,&nbsp;Yijun Wang,&nbsp;Junlin Cheng,&nbsp;Lu Luo,&nbsp;Yueru Li","doi":"10.1016/j.ijthermalsci.2025.109860","DOIUrl":"10.1016/j.ijthermalsci.2025.109860","url":null,"abstract":"<div><div>To explore new efficient cooling technology for advanced gas turbine blades and reduce dependence on film cooling, this paper proposes a novel pipe network cooling structure. The design connects leading-edge impingement cooling holes to trailing-edge slits through lateral pipes and incorporates independent vertical pipes to form a network structure. This cooling structure can be applied to a complete blade cooling system, demonstrating strong cooling performance in the mid-chord region despite the absence of film holes, while achieving a more uniform overall temperature distribution, showing promising developmental potential. Through experimental and numerical simulations, comparisons were made with typical gas turbine blade cooling structures and double-wall cooling structures. The results indicate that this new pipes network cooling structure offers superior cooling performance and achieves a more uniform temperature distribution. In addition, the study investigated the impact of lateral pipes shapes and the distances between transverse and vertical pipes relative to the end wall on cooling performance. The results showed that, under the same boundary conditions, hexagonal pipes performed better. The relative positions of transverse and vertical pipes significantly affected blade cooling efficiency. P₁/P₂ = 0.5, the temperature distribution was the most uniform; P₁/P₂ = 1, heat transfer in the mid-chord region improved.</div></div>","PeriodicalId":341,"journal":{"name":"International Journal of Thermal Sciences","volume":"214 ","pages":"Article 109860"},"PeriodicalIF":4.9,"publicationDate":"2025-03-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143628142","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}