Travis Leadbetter, Prashant K. Purohit, Celia Reina
Given a particle system obeying overdamped Langevin dynamics, we demonstrate that it is always possible to construct a thermodynamically consistent macroscopic model which obeys a gradient flow with respect to its non-equilibrium free energy. To do so, we significantly extend the recent Stochastic Thermodynamics with Internal Variables (STIV) framework, a method for producing macroscopic thermodynamic models far-from-equilibrium from the underlying mesoscopic dynamics and an approximate probability density of states parameterized with so-called internal variables. Though originally explored for Gaussian probability distributions, we here allow for an arbitrary choice of the approximate probability density while retaining a gradient flow dynamics. This greatly extends its range of applicability and automatically ensures consistency with the second law of thermodynamics, without the need for secondary verification. We demonstrate numerical convergence, in the limit of increasing internal variables, to the true probability density of states for both a multi-modal relaxation problem, a protein diffusing on a strand of DNA, and for an externally driven particle in a periodic landscape. Finally, we provide a reformulation of STIV with the quasi-equilibrium approximations in terms of the averages of observables of the mesostate, and show that these, too, obey a gradient flow.
{"title":"From Langevin dynamics to macroscopic thermodynamic models: a general framework valid far from equilibrium","authors":"Travis Leadbetter, Prashant K. Purohit, Celia Reina","doi":"10.1515/jnet-2025-0071","DOIUrl":"https://doi.org/10.1515/jnet-2025-0071","url":null,"abstract":"Given a particle system obeying overdamped Langevin dynamics, we demonstrate that it is always possible to construct a thermodynamically consistent macroscopic model which obeys a gradient flow with respect to its non-equilibrium free energy. To do so, we significantly extend the recent Stochastic Thermodynamics with Internal Variables (STIV) framework, a method for producing macroscopic thermodynamic models far-from-equilibrium from the underlying mesoscopic dynamics and an approximate probability density of states parameterized with so-called internal variables. Though originally explored for Gaussian probability distributions, we here allow for an arbitrary choice of the approximate probability density while retaining a gradient flow dynamics. This greatly extends its range of applicability and automatically ensures consistency with the second law of thermodynamics, without the need for secondary verification. We demonstrate numerical convergence, in the limit of increasing internal variables, to the true probability density of states for both a multi-modal relaxation problem, a protein diffusing on a strand of DNA, and for an externally driven particle in a periodic landscape. Finally, we provide a reformulation of STIV with the quasi-equilibrium approximations in terms of the averages of observables of the mesostate, and show that these, too, obey a gradient flow.","PeriodicalId":16428,"journal":{"name":"Journal of Non-Equilibrium Thermodynamics","volume":"95 1","pages":""},"PeriodicalIF":6.6,"publicationDate":"2025-10-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145311559","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
In this work, a theoretical study is carried out on the effects of the thickness of a semiconductor thin film on the transport of heat and particles under the action of an external temperature difference. The dependence of the Seebeck effect on the thickness is considered. The thickness film is introduced through the transport coefficients of the material, namely, the thermal and electrical conductivities and the intrinsic Seebeck coefficient, by resorting to known results of irreversible thermodynamics. Graphs of the generated electric potential difference by an external temperature difference versus film thickness are obtained in the range 0.1 µm–1 µm. We compare the results of two slightly doped materials, namely, silicon and bismuth telluride of the n- and p- types. The n-type silicon shows an optimal thickness where the generated electric potential difference is maximum, while the electric potential difference in the n-type bismuth telluride decreases with decreasing thickness. The electric response of p-type silicon and p-type bismuth telluride also worsens as the thickness decreases. The results presented may be useful in the design of thermoelectric devices on the sub-micrometer length scale.
{"title":"Size scaling effects on heat and particle transport in semiconductor thin films: a near-equilibrium thermodynamic approach","authors":"Ruth Estephania Gonzalez-Narvaez, Iván Rivera, Víctor Hernández, Aldo Figueroa, Federico Vázquez","doi":"10.1515/jnet-2025-0056","DOIUrl":"https://doi.org/10.1515/jnet-2025-0056","url":null,"abstract":"In this work, a theoretical study is carried out on the effects of the thickness of a semiconductor thin film on the transport of heat and particles under the action of an external temperature difference. The dependence of the Seebeck effect on the thickness is considered. The thickness film is introduced through the transport coefficients of the material, namely, the thermal and electrical conductivities and the intrinsic Seebeck coefficient, by resorting to known results of irreversible thermodynamics. Graphs of the generated electric potential difference by an external temperature difference versus film thickness are obtained in the range 0.1 µm–1 µm. We compare the results of two slightly doped materials, namely, silicon and bismuth telluride of the <jats:italic>n</jats:italic>- and <jats:italic>p</jats:italic>- types. The <jats:italic>n</jats:italic>-type silicon shows an optimal thickness where the generated electric potential difference is maximum, while the electric potential difference in the <jats:italic>n</jats:italic>-type bismuth telluride decreases with decreasing thickness. The electric response of <jats:italic>p</jats:italic>-type silicon and <jats:italic>p</jats:italic>-type bismuth telluride also worsens as the thickness decreases. The results presented may be useful in the design of thermoelectric devices on the sub-micrometer length scale.","PeriodicalId":16428,"journal":{"name":"Journal of Non-Equilibrium Thermodynamics","volume":"117 1","pages":""},"PeriodicalIF":6.6,"publicationDate":"2025-10-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145282849","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Konrad Giżyński, Karol Makuch, Jan Paczesny, Paweł Żuk, Anna Maciołek, Robert Hołyst
We studied planar compressible Poiseuille flows of an ideal gas, both in steady and unsteady states, to identify the minimal number of state parameters required to describe changes in internal energy. In previous work (Phys. Rev. E 104, 055107 (2021)), five parameters were needed for steady flows. Here, using global non-equilibrium thermodynamics, we reduce this number to three: non-equilibrium entropy S*, volume V, and number of particles N. The internal energy U(S*, V, N) of such systems in stationary and non-stationary states is the function of non-equilibrium entropy S*, volume V and number of particles N in the system irrespective of any processes, number of boundary conditions or imposed constraints. We tested this by placing a cylinder inside the channel, finding that U depends on the cylinder’s location yc only via the state parameters S*(yc) and N(yc) for V = const. Moreover, in cases where the flow becomes unstable and parameters such as velocity and pressure oscillate, U depends on time t only through S*(t) and N(t) for V = const. These results demonstrate that this formulation of internal energy remains robust and consistent, even in unsteady flows with varying boundary conditions.
{"title":"The internal energy as a function of state parameters in steady and unsteady Poiseuille flows","authors":"Konrad Giżyński, Karol Makuch, Jan Paczesny, Paweł Żuk, Anna Maciołek, Robert Hołyst","doi":"10.1515/jnet-2025-0003","DOIUrl":"https://doi.org/10.1515/jnet-2025-0003","url":null,"abstract":"We studied planar compressible Poiseuille flows of an ideal gas, both in steady and unsteady states, to identify the minimal number of state parameters required to describe changes in internal energy. In previous work (Phys. Rev. E 104, 055107 (2021)), five parameters were needed for steady flows. Here, using global non-equilibrium thermodynamics, we reduce this number to three: non-equilibrium entropy <jats:italic>S</jats:italic> <jats:sup>*</jats:sup>, volume <jats:italic>V</jats:italic>, and number of particles <jats:italic>N</jats:italic>. The internal energy <jats:italic>U</jats:italic>(<jats:italic>S</jats:italic> <jats:sup>*</jats:sup>, <jats:italic>V</jats:italic>, <jats:italic>N</jats:italic>) of such systems in stationary and non-stationary states is the function of non-equilibrium entropy <jats:italic>S</jats:italic> <jats:sup>*</jats:sup>, volume <jats:italic>V</jats:italic> and number of particles <jats:italic>N</jats:italic> in the system irrespective of any processes, number of boundary conditions or imposed constraints. We tested this by placing a cylinder inside the channel, finding that <jats:italic>U</jats:italic> depends on the cylinder’s location <jats:italic>y</jats:italic> <jats:sub> <jats:italic>c</jats:italic> </jats:sub> only via the state parameters <jats:italic>S</jats:italic> <jats:sup>*</jats:sup>(<jats:italic>y</jats:italic> <jats:sub> <jats:italic>c</jats:italic> </jats:sub>) and <jats:italic>N</jats:italic>(<jats:italic>y</jats:italic> <jats:sub> <jats:italic>c</jats:italic> </jats:sub>) for <jats:italic>V</jats:italic> = const. Moreover, in cases where the flow becomes unstable and parameters such as velocity and pressure oscillate, <jats:italic>U</jats:italic> depends on time <jats:italic>t</jats:italic> only through <jats:italic>S</jats:italic> <jats:sup>*</jats:sup>(<jats:italic>t</jats:italic>) and <jats:italic>N</jats:italic>(<jats:italic>t</jats:italic>) for <jats:italic>V</jats:italic> = const. These results demonstrate that this formulation of internal energy remains robust and consistent, even in unsteady flows with varying boundary conditions.","PeriodicalId":16428,"journal":{"name":"Journal of Non-Equilibrium Thermodynamics","volume":"29 1","pages":""},"PeriodicalIF":6.6,"publicationDate":"2025-09-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145078044","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Extended irreversible thermodynamics is a theory that expands the classical framework of nonequilibrium thermodynamics by going beyond the local-equilibrium assumption. A notable example of this is the Maxwell–Cattaneo heat flux model, which introduces a time lag in the heat flux response to temperature gradients. In this paper, we develop a variational formulation of the equations of extended irreversible thermodynamics by introducing an action principle for a nonequilibrium Lagrangian that treats thermodynamic fluxes as independent variables. A key feature of this approach is that it naturally extends both Hamilton’s principle of reversible continuum mechanics and the earlier variational formulation of classical irreversible thermodynamics. The variational principle is initially formulated in the material (Lagrangian) description, from which the Eulerian form is derived using material covariance (or relabeling symmetries). The tensorial structure of the thermodynamic fluxes dictates the choice of objective rate in the Eulerian description, and plays a central role in the emergence of nonequilibrium stresses – arising from both viscous and thermal effects – that are essential to ensure thermodynamic consistency. This framework naturally results in the Cattaneo–Christov model for heat flux. We also investigate the extension of the approach to accommodate higher-order fluxes and the general form of entropy fluxes. The variational framework presented in this paper has promising applications in the development of structure-preserving and thermodynamically consistent numerical methods. It is particularly relevant for modeling systems where entropy production is a delicate issue that requires careful treatment to ensure consistency with the laws of thermodynamics.
{"title":"A variational principle for extended irreversible thermodynamics: heat conducting viscous fluids","authors":"François Gay-Balmaz","doi":"10.1515/jnet-2025-0022","DOIUrl":"https://doi.org/10.1515/jnet-2025-0022","url":null,"abstract":"Extended irreversible thermodynamics is a theory that expands the classical framework of nonequilibrium thermodynamics by going beyond the local-equilibrium assumption. A notable example of this is the Maxwell–Cattaneo heat flux model, which introduces a time lag in the heat flux response to temperature gradients. In this paper, we develop a variational formulation of the equations of extended irreversible thermodynamics by introducing an action principle for a nonequilibrium Lagrangian that treats thermodynamic fluxes as independent variables. A key feature of this approach is that it naturally extends both Hamilton’s principle of reversible continuum mechanics and the earlier variational formulation of classical irreversible thermodynamics. The variational principle is initially formulated in the material (Lagrangian) description, from which the Eulerian form is derived using material covariance (or relabeling symmetries). The tensorial structure of the thermodynamic fluxes dictates the choice of objective rate in the Eulerian description, and plays a central role in the emergence of nonequilibrium stresses – arising from both viscous and thermal effects – that are essential to ensure thermodynamic consistency. This framework naturally results in the Cattaneo–Christov model for heat flux. We also investigate the extension of the approach to accommodate higher-order fluxes and the general form of entropy fluxes. The variational framework presented in this paper has promising applications in the development of structure-preserving and thermodynamically consistent numerical methods. It is particularly relevant for modeling systems where entropy production is a delicate issue that requires careful treatment to ensure consistency with the laws of thermodynamics.","PeriodicalId":16428,"journal":{"name":"Journal of Non-Equilibrium Thermodynamics","volume":"78 1","pages":""},"PeriodicalIF":6.6,"publicationDate":"2025-09-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144928613","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Umar Farooq, Hafiz Hamza Riaz, Tauqir Muhammad, Samar Ali, Tzu Chi Chan
Despite advancements in cooling solutions for electronic devices, heat dissipation remains the primary challenge in optimizing heat sink performance in a competitive industry. The current study numerically investigates the performance of plate-fin heat sink (PHS) and pin-fin heat sink (PnHS) using a hybrid nanofluid (GnP-MWCNT/Water) as the working fluid. Key performance parameters, including pressure drop, thermal resistance, effectivity, and turbulent kinetic energy, are analyzed across different Re and nanofluid concentrations. Findings indicate that while the PnHS exhibits higher convective heat transfer due to increased flow disturbances, it also suffers from greater thermal resistance and pressure drop. In contrast, the PHS offers an adequate balance between heat dissipation and flow efficiency, leading to higher overall effectiveness. At Re = 5,334, the pressure drop for 0.20 % nanofluid is 88.5 Pa in the PnHS and 73 Pa in the PHS. Additionally, at Re = 1,333, the effectiveness values for PHS are 0.333 (water) and 0.354 (nanofluid), while for PnHS, they are lower at 0.186 and 0.195, respectively. The current study highlights the interplay between enhanced heat transfer and increased flow resistance, emphasizing the importance of optimizing fin design and nanofluid concentration for efficient thermal management.
{"title":"Heat transfer analysis of plate versus pin fin heat sinks with GnP-MWCNT/water hybrid nanofluid","authors":"Umar Farooq, Hafiz Hamza Riaz, Tauqir Muhammad, Samar Ali, Tzu Chi Chan","doi":"10.1515/jnet-2025-0033","DOIUrl":"https://doi.org/10.1515/jnet-2025-0033","url":null,"abstract":"Despite advancements in cooling solutions for electronic devices, heat dissipation remains the primary challenge in optimizing heat sink performance in a competitive industry. The current study numerically investigates the performance of plate-fin heat sink (PHS) and pin-fin heat sink (PnHS) using a hybrid nanofluid (GnP-MWCNT/Water) as the working fluid. Key performance parameters, including pressure drop, thermal resistance, effectivity, and turbulent kinetic energy, are analyzed across different Re and nanofluid concentrations. Findings indicate that while the PnHS exhibits higher convective heat transfer due to increased flow disturbances, it also suffers from greater thermal resistance and pressure drop. In contrast, the PHS offers an adequate balance between heat dissipation and flow efficiency, leading to higher overall effectiveness. At Re = 5,334, the pressure drop for 0.20 % nanofluid is 88.5 Pa in the PnHS and 73 Pa in the PHS. Additionally, at Re = 1,333, the effectiveness values for PHS are 0.333 (water) and 0.354 (nanofluid), while for PnHS, they are lower at 0.186 and 0.195, respectively. The current study highlights the interplay between enhanced heat transfer and increased flow resistance, emphasizing the importance of optimizing fin design and nanofluid concentration for efficient thermal management.","PeriodicalId":16428,"journal":{"name":"Journal of Non-Equilibrium Thermodynamics","volume":"31 1","pages":""},"PeriodicalIF":6.6,"publicationDate":"2025-08-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144899995","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
The current trend of reducing the size of electronic devices in the industry has extensively increased the demand for effective heat dissipation, thereby intensifying the need for high-performance heat-dissipating devices. A promising approach to solve this challenge is the use of single-phase (SP), two-phase (TP), and supercritical fluids in micro-channels (MCs). Two-phase cooling is applicable only to those devices in which the tip temperature is high enough to allow the cooling fluid to convert into a two-phase state. In all other cases, only single-phase cooling can be utilized. In this work, numerical and experimental investigations on MC have been performed using water as the working fluid to predict TP behavior and heat dissipation from electronic devices using SP and TP flow. A numerical model of flow boiling heat transfer was developed based on conservation equations, which is solved to identify the existence of single and two-phase regions in the MC and to study the variation of pressure along its length at different heating powers. Further, experiments were performed in both SP and TP conditions to observe the nature of flow regimes and the impact of various parameters on effective heat dissipation through MCs well as temperature distribution. Numerical results were validated with experimental results, which showed good agreement. Several experiments were also carried out to develop an empirical correlation between mass flow rate and heat power to maintain the electronic device temperature below 40 °C. The developed correlation is experimentally validated at three different heat powers 6 W, 8 W and 10 W.
{"title":"Numerical and experimental heat transfer analysis of two-phase flow through microchannel for development of heat dissipation correlation","authors":"Santosh Kumar Rai, Vikas Goyat, Mahesh Kumar Gupta, Gyander Ghangas, Dhowmya Bhatt, Arun Uniyal, Pardeep Kumar, Nikhil Vivek Shrivas","doi":"10.1515/jnet-2025-0044","DOIUrl":"https://doi.org/10.1515/jnet-2025-0044","url":null,"abstract":"The current trend of reducing the size of electronic devices in the industry has extensively increased the demand for effective heat dissipation, thereby intensifying the need for high-performance heat-dissipating devices. A promising approach to solve this challenge is the use of single-phase (SP), two-phase (TP), and supercritical fluids in micro-channels (MCs). Two-phase cooling is applicable only to those devices in which the tip temperature is high enough to allow the cooling fluid to convert into a two-phase state. In all other cases, only single-phase cooling can be utilized. In this work, numerical and experimental investigations on MC have been performed using water as the working fluid to predict TP behavior and heat dissipation from electronic devices using SP and TP flow. A numerical model of flow boiling heat transfer was developed based on conservation equations, which is solved to identify the existence of single and two-phase regions in the MC and to study the variation of pressure along its length at different heating powers. Further, experiments were performed in both SP and TP conditions to observe the nature of flow regimes and the impact of various parameters on effective heat dissipation through MCs well as temperature distribution. Numerical results were validated with experimental results, which showed good agreement. Several experiments were also carried out to develop an empirical correlation between mass flow rate and heat power to maintain the electronic device temperature below 40 °C. The developed correlation is experimentally validated at three different heat powers 6 W, 8 W and 10 W.","PeriodicalId":16428,"journal":{"name":"Journal of Non-Equilibrium Thermodynamics","volume":"69 1","pages":""},"PeriodicalIF":6.6,"publicationDate":"2025-08-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144778550","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Lactose permease, a secondary active transporter from Escherichia coli, facilitates the co-transport of protons and lactose across the cytoplasmic membrane. Unlike passive diffusion mechanisms, lactose permease operates via conformational switching that alternately exposes the binding pocket to either membrane side. In this study, we present a theoretical treatment combining irreversible thermodynamic principles and kinetic modeling to quantify its operation. Onsager’s reciprocity is applied to analyze proton-lactose coupling, and a bisubstrate kinetic scheme is employed to simulate system behavior under various proton gradients and lactose concentrations. The complete catalytic cycle is characterized by associated rate constants and energetic transitions, highlighting that lactose permease exhibits a dissipative nature as a hallmark of secondary active transport. Altogether, this study provides a novel thermodynamic perspective on lactose permease, aiming to bridge molecular transport kinetics with the formalism of irreversible processes. This work is the first to integrate Onsager relations with the Michaelis-Menten kinetic model to quantify the energetic efficiency of lactose permease.
{"title":"Applying irreversible thermodynamics to the paradigmatic secondary transporter: lactose permease (LacY)","authors":"Jordi H. Borrell","doi":"10.1515/jnet-2025-0054","DOIUrl":"https://doi.org/10.1515/jnet-2025-0054","url":null,"abstract":"Lactose permease, a secondary active transporter from <jats:italic>Escherichia coli</jats:italic>, facilitates the co-transport of protons and lactose across the cytoplasmic membrane. Unlike passive diffusion mechanisms, lactose permease operates via conformational switching that alternately exposes the binding pocket to either membrane side. In this study, we present a theoretical treatment combining irreversible thermodynamic principles and kinetic modeling to quantify its operation. Onsager’s reciprocity is applied to analyze proton-lactose coupling, and a bisubstrate kinetic scheme is employed to simulate system behavior under various proton gradients and lactose concentrations. The complete catalytic cycle is characterized by associated rate constants and energetic transitions, highlighting that lactose permease exhibits a dissipative nature as a hallmark of secondary active transport. Altogether, this study provides a novel thermodynamic perspective on lactose permease, aiming to bridge molecular transport kinetics with the formalism of irreversible processes. This work is the first to integrate Onsager relations with the Michaelis-Menten kinetic model to quantify the energetic efficiency of lactose permease.","PeriodicalId":16428,"journal":{"name":"Journal of Non-Equilibrium Thermodynamics","volume":"159 1","pages":""},"PeriodicalIF":6.6,"publicationDate":"2025-08-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144786588","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
High-temperature proton exchange membrane fuel cells (HT-PEMFCs) inherently produce waste heat, leading to component degradation, increased cooling demands, and reduced efficiency and longevity. To mitigate these challenges, this study introduces isopropanol-acetone-hydrogen chemical heat pumps (IAH-CHPs), selected for their proven ability to efficiently upgrade and store the waste heat from HT-PEMFCs in a high-value form. Grounded in thermodynamic and electrochemical principles, a comprehensive mathematical model, incorporating key irreversible losses, is developed to evaluate the potential. Numerical calculations predict a 29 % increase in the hybrid system’s maximum power density compared to a standalone HT-PEMFC operating at 443 K, along with corresponding enhancements of 14.17 % in energy efficiency and 14.16 % in exergy efficiency. Preliminary predictions confirm the feasibility of this approach, and the optimal operating ranges for maximizing power density are identified. Additionally, exhaustive parametric studies reveal the impacts of various structural and operational parameters – such as leakage current density, phosphoric acid doping, relative humidity, operating temperatures, and critical factors within the heat pump cycle – on the system’s thermodynamic performance and key current density indicators. Local sensitivity analyses highlight effective performance regulation strategies. These results provide essential insights for mitigating waste heat challenges, enhancing system efficiency, and extending the operational lifespan for HT-PEMFCs.
{"title":"A Hydrogen-fueled hybrid system based on HT-PEMFCs for simultaneous electrical power generation and high-value heat storage","authors":"Houcheng Zhang, Han Wang, Min Kuang, Yejian Xue","doi":"10.1515/jnet-2024-0122","DOIUrl":"https://doi.org/10.1515/jnet-2024-0122","url":null,"abstract":"High-temperature proton exchange membrane fuel cells (HT-PEMFCs) inherently produce waste heat, leading to component degradation, increased cooling demands, and reduced efficiency and longevity. To mitigate these challenges, this study introduces isopropanol-acetone-hydrogen chemical heat pumps (IAH-CHPs), selected for their proven ability to efficiently upgrade and store the waste heat from HT-PEMFCs in a high-value form. Grounded in thermodynamic and electrochemical principles, a comprehensive mathematical model, incorporating key irreversible losses, is developed to evaluate the potential. Numerical calculations predict a 29 % increase in the hybrid system’s maximum power density compared to a standalone HT-PEMFC operating at 443 K, along with corresponding enhancements of 14.17 % in energy efficiency and 14.16 % in exergy efficiency. Preliminary predictions confirm the feasibility of this approach, and the optimal operating ranges for maximizing power density are identified. Additionally, exhaustive parametric studies reveal the impacts of various structural and operational parameters – such as leakage current density, phosphoric acid doping, relative humidity, operating temperatures, and critical factors within the heat pump cycle – on the system’s thermodynamic performance and key current density indicators. Local sensitivity analyses highlight effective performance regulation strategies. These results provide essential insights for mitigating waste heat challenges, enhancing system efficiency, and extending the operational lifespan for HT-PEMFCs.","PeriodicalId":16428,"journal":{"name":"Journal of Non-Equilibrium Thermodynamics","volume":"1 1","pages":""},"PeriodicalIF":6.6,"publicationDate":"2025-07-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144747576","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Fuzhang Wang, Isaac Lare Animasaun, Taseer Muhammad
Accurately predicting turbulent water flow in duct systems remains a challenging problem, particularly when anisotropic turbulence effects are significant. Bridging the gap between industrial applications and academic research requires a deeper understanding of such complex flows. This study investigates a less commonly analyzed configuration involving a horizontal aluminum duct transitioning into a converging wavy duct. The wavy section consists of 2.5 full sinusoidal periods, ending in a reduced outlet diameter. In addition, the effect of incorporating four minor/secondary inlets, arranged as branches at different angles, was examined and presented herein. Aluminum was selected for its low density and corrosion resistance, which are beneficial in experimental and industrial setups. Initially, the duct was analyzed in an unbranched configuration. The study then progressed to include the four secondary/minor branch inlets at various angles. The simulation results were validated by comparison with a solution for a simple flow in a 70 mm duct. Additional verification was provided by employing other CFD codes, along with grid convergence index and mesh sensitivity analyses, improving the confidence in the simulation results. Branch angles influences turbulence intensity depending on flow conditions and angle magnitude. Sharper branch angles are particularly effective, inducing greater turbulence at the converged outlet. Higher inlet temperatures and velocities lead to increased Reynolds stress due to enhanced energy transfer and elevated turbulent kinetic energy. Specifically, an increase in inlet velocity at a 45° branch angle further augments turbulent momentum transfer, resulting in more controlled mixing along the duct.
{"title":"Anisotropic turbulent flow of water through converging wavy-aluminum-circular pipe with five half-cycles: insight into the significance of four-branch minor-inlet angle","authors":"Fuzhang Wang, Isaac Lare Animasaun, Taseer Muhammad","doi":"10.1515/jnet-2025-0046","DOIUrl":"https://doi.org/10.1515/jnet-2025-0046","url":null,"abstract":"Accurately predicting turbulent water flow in duct systems remains a challenging problem, particularly when anisotropic turbulence effects are significant. Bridging the gap between industrial applications and academic research requires a deeper understanding of such complex flows. This study investigates a less commonly analyzed configuration involving a horizontal aluminum duct transitioning into a converging wavy duct. The wavy section consists of 2.5 full sinusoidal periods, ending in a reduced outlet diameter. In addition, the effect of incorporating four minor/secondary inlets, arranged as branches at different angles, was examined and presented herein. Aluminum was selected for its low density and corrosion resistance, which are beneficial in experimental and industrial setups. Initially, the duct was analyzed in an unbranched configuration. The study then progressed to include the four secondary/minor branch inlets at various angles. The simulation results were validated by comparison with a solution for a simple flow in a 70 mm duct. Additional verification was provided by employing other CFD codes, along with grid convergence index and mesh sensitivity analyses, improving the confidence in the simulation results. Branch angles influences turbulence intensity depending on flow conditions and angle magnitude. Sharper branch angles are particularly effective, inducing greater turbulence at the converged outlet. Higher inlet temperatures and velocities lead to increased Reynolds stress due to enhanced energy transfer and elevated turbulent kinetic energy. Specifically, an increase in inlet velocity at a 45<jats:italic>°</jats:italic> branch angle further augments turbulent momentum transfer, resulting in more controlled mixing along the duct.","PeriodicalId":16428,"journal":{"name":"Journal of Non-Equilibrium Thermodynamics","volume":"69 1","pages":""},"PeriodicalIF":6.6,"publicationDate":"2025-07-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144533270","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Yasir Ul Umair Bin Turabi, Zeshan Faiz, Shahzad Munir, Shafee Ahmad, M.M. Alam, Hakim AL Garalleh
Enclosure design is essential for thermal engineering technology and applications, including electronics, heat transfer equipment, power reactors, cooling mechanisms, solar energy systems, and nuclear power plants. This study aims to analyze the numerical and Response Surface Methodology (RSM) optimization for natural convection and entropy generation in a wavy triangular cavity with Casson fluid under inclined magnetohydrodynamic and radiation influences. The finite element approach (FEM) is utilized to compute the numerical solution for the simulation framework, while RSM is applied to determine the optimal heat transfer rate among four different parameters. The study presents streamlines, velocity profiles, isothermal lines, total entropy generation, and average Nusselt number in graphically and tabularly. The results show that an increase in the number of undulations and the Casson parameter leads to an increase in the thermal transfer rate and total entropy generation, whereas the Hartmann number has a decreasing effect on both. The Nusselt number rises with the rising number of undulations and the radiation parameter. The peak stream function is observed at an inclination angle of 60°. The significant R2 value of 0.9967 shows a good agreement between the expected and actual values.
{"title":"Numerical and optimization analysis of natural convection and entropy-generation in wavy triangular cavity with Casson fluid under magnetohydrodynamics and radiation","authors":"Yasir Ul Umair Bin Turabi, Zeshan Faiz, Shahzad Munir, Shafee Ahmad, M.M. Alam, Hakim AL Garalleh","doi":"10.1515/jnet-2024-0101","DOIUrl":"https://doi.org/10.1515/jnet-2024-0101","url":null,"abstract":"Enclosure design is essential for thermal engineering technology and applications, including electronics, heat transfer equipment, power reactors, cooling mechanisms, solar energy systems, and nuclear power plants. This study aims to analyze the numerical and Response Surface Methodology (RSM) optimization for natural convection and entropy generation in a wavy triangular cavity with Casson fluid under inclined magnetohydrodynamic and radiation influences. The finite element approach (FEM) is utilized to compute the numerical solution for the simulation framework, while RSM is applied to determine the optimal heat transfer rate among four different parameters. The study presents streamlines, velocity profiles, isothermal lines, total entropy generation, and average Nusselt number in graphically and tabularly. The results show that an increase in the number of undulations and the Casson parameter leads to an increase in the thermal transfer rate and total entropy generation, whereas the Hartmann number has a decreasing effect on both. The Nusselt number rises with the rising number of undulations and the radiation parameter. The peak stream function is observed at an inclination angle of 60°. The significant <jats:italic>R</jats:italic> <jats:sup>2</jats:sup> value of 0.9967 shows a good agreement between the expected and actual values.","PeriodicalId":16428,"journal":{"name":"Journal of Non-Equilibrium Thermodynamics","volume":"26 1","pages":""},"PeriodicalIF":6.6,"publicationDate":"2025-06-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144488465","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}