Heat Transfer Prediction With Unknown Thermal Boundary Conditions Using Physics-Informed Neural Networks

Author(s):  
Shengze Cai ◽  
Zhicheng Wang ◽  
Chryssostomos Chryssostomidis ◽  
George Em Karniadakis

Abstract Simulating convective heat transfer using traditional numerical methods requires explicit definition of thermal boundary conditions on all boundaries of the domain, which is almost impossible to fulfill in real applications. Here, we address this ill-posed problem using machine learning techniques by assuming that we have some extra measurements of the temperature at a few locations in the domain, not necessarily located on the boundaries with the unknown thermal boundary condition. In particular, we employ physics-informed neural networks (PINNs) to represent the velocity and temperature fields while simultaneously enforce the Navier-Stokes and energy equations at random points in the domain. In PINNs, all differential operators are computed using automatic differentiation, hence avoiding discretization in either space or time. The loss function is composed of multiple terms, including the mismatch in the velocity and temperature data, the boundary and initial conditions, as well as the residuals of the Navier-Stokes and energy equations. Here, we develop a data-driven strategy based on PINNs to infer the temperature field in the prototypical problem of convective heat transfer in flow past a cylinder. We assume that we have just a couple of temperature measurements on the cylinder surface and a couple more temperature measurements in the wake region, but the thermal boundary condition on the cylinder surface is totally unknown. Upon training the PINN, we can discover the unknown boundary condition while simultaneously infer the temperature field everywhere in the domain with less than 5% error in the Nusselt number prediction. In order to assess the performance of PINN, we carried out a high fidelity simulation of the same heat transfer problem (with known thermal boundary conditions) by using the high-order spectral/hp-element method (SEM), and quantitatively evaluated the accuracy of PINN’s prediction with respect to SEM. We also propose a method to adaptively select the location of sensors in order to minimize the number of required temperature measurements while increasing the accuracy of the inference in heat transfer.

Author(s):  
H. D. Madhawa Hettiarachchi ◽  
Mihajlo Golubovic ◽  
William M. Worek

Slip-flow and heat transfer in rectangular microchannels are studied numerically for constant wall temperature (T) and constant wall heat flux (H2) boundary conditions under thermally developing flow. Navier-Stokes and energy equations with velocity slip and temperature jump at the boundary are solved using finite volume method in a three dimensional cartesian coordinate system. A modified convection-diffusion coefficient at the wall-fluid interface is defined to incorporate the temperature-jump boundary condition. Validity of the numerical simulation procedure is stabilized. The effect of rarefaction on heat transfer in the entrance region is analyzed in detail. The velocity slip has an increasing effect on the Nusselt (Nu) number whereas temperature jump has a decreasing effect, and the combined effect could result increase or decrease in the Nu number. For the range of parameters considered, there could be high as 15% increase or low as 50% decrease in fully developed Nu is plausible for T thermal boundary condition while it could be high as 20% or low as 35% for H2 thermal boundary condition.


2021 ◽  
Author(s):  
Janet Peifer ◽  
Onno Bokhove ◽  
Steve Tobias

<p>Rayleigh-Bénard convection (RBC) is a fluid phenomenon that has been studied for over a century because of its utility in simplifying very complex physical systems. Many geophysical and astrophysical systems, including planetary core dynamics and components of weather prediction, are modeled by including rotational forcing in classic RBC. Our understanding of these systems is confined by experimental and numerical limits, as well as theoretical assumptions. </p><p>The role of thermal boundary condition choice on experimental studies of geophysical and astrophysical systems has been often been overlooked, which could account for some lack of agreement between experimental and numerical models as well as the actual flows. The typical thermal boundary conditions prescribed at the top and the bottom of a convection system are fixed temperature conditions, despite few real geophysical systems being bounded with a fixed temperature. A constant heat flux is generally more applicable for real large-scale geophysical systems. However, when this condition is applied in numerical systems, the lack of fixed temperature can cause a temperature drift. In this study, we seek to minimize temperature drifting by applying a fixed temperature condition on one boundary and a fixed thermal flux on the other.</p><p>Experimental boundary conditions are also often assumed to be a fixed temperature. However, the actual condition is determined by the ratio of the height and thermal conductivity of the boundary material to that of the contained fluid, known as the Biot number. The relationship between the Biot number and thermal boundary condition behavior is defined by the Robin, or 'thin-lid', boundary condition such that low Biot number boundaries are essentially fixed thermal flux and high Biot number boundaries are essentially fixed temperature. </p><p>This study seeks to strengthen the link between numerical and experimental models and geophysical flows by investigating the effects of thermal boundary conditions and their relationship to real-world processes. Both fixed temperature and fixed flux boundary conditions are considered. In addition, the Robin boundary condition is studied at a range of Biot numbers spanning from fixed temperature to fixed flux, allowing intermediate conditions to be investigated. Each system is studied at increasingly rapid rotation rates, corresponding to decreasing Ekman numbers as low as Ek=10<sup>-5</sup> Heat transport is analyzed using the Nusselt number, Nu, and the form of the solution is described by the number of convection rolls and time-dependency. Further investigations will analyze Nu and fluid movement within a system with heterogeneous heat flux condition on the  sidewall boundary conditions, which is useful in the study of planetary core dynamics. The results of this study have implications for improvements in modeling geophysical systems both experimentally and numerically. </p>


Author(s):  
Masoud Darbandi ◽  
Shidvash Vakilipour

In this work, we present a more realistic inlet boundary condition to simulate compressible and incompressible flows through micro and nano channels considering consistent momentum and heat transfer specifications there. At solid walls, a constant wall temperature with suitable jump is applied as the wall thermal boundary condition; however, two types of thermal inlet boundary conditions are investigated at the inlet. We firstly examine the classical inlet boundary condition, which specifies a uniform temperature distribution right at the real inlet. Alternatively, we apply the same boundary condition but at a fictitious place far upstream of the real channel inlet. To validate our results, the results obtained after employing the former boundary conditions are evaluated against other available numerical results and Lattice Boltzmann solutions. In this validation, the friction factor and Nusselt number are treated as the most important hydrodynamics and heat transfer factors to be appraised. Next, we present the results after applying the second type of inlet boundary conditions and compare them with those of applying the first type of boundary condition. This study suggests an innovation in micro/nano scale heat transfer treatment close to the channel inlets.


Energies ◽  
2021 ◽  
Vol 14 (20) ◽  
pp. 6448
Author(s):  
Raoudha Chaabane ◽  
Annunziata D’Orazio ◽  
Abdelmajid Jemni ◽  
Arash Karimipour ◽  
Ramin Ranjbarzadeh

In recent decades, research utilizing numerical schemes dealing with fluid and nanoparticle interaction has been relatively intensive. It is known that CuO nanofluid with a volume fraction of 0.1 and a special thermal boundary condition with heat supplied to part of the wall increases the average Nusselt number for different aspect ratios ranges and for high Rayleigh numbers. Due to its simplicity, stability, accuracy, efficiency, and ease of parallelization, we use the thermal single relaxation time Bhatnagar-Gross-Krook (SRT BGK) mesoscopic approach D2Q9 scheme lattice Boltzmann method in order to solve the coupled Navier–Stokes equations. Convection of CuO nanofluid in a square enclosure with a moderate Rayleigh number of 105 and with new boundary conditions is highlighted. After a successful validation with a simple partial Dirichlet boundary condition, this paper extends the study to deal with linear and sinusoidal thermal boundary conditions applied to part of the wall.


2014 ◽  
Vol 522-524 ◽  
pp. 1707-1712 ◽  
Author(s):  
Qing Long Peng ◽  
Zhao Hui Qi ◽  
Xia Gan ◽  
Chao Li

How to use numerical simulation method to analyze human body heat transfer and human thermal comfort is introduced in this paper systematically. Under the same working conditions, numerical simulation of human body heat transfer has been finished based on three thermal boundary conditions, and then the results are compared. The results show that the third thermal boundary condition is better than the first and the second one, which have some problems in simulation and are not good at reflecting the fact on thermal comfort of human body. The third thermal boundary condition which is made to adapt the surrounding flow field automatically can get a more accurate result on calculating the heat transfer of different parts on human body and reflect hot or cool feeling preferably, which proves that the method put forward in this article to research the human body comfort is feasible.


1992 ◽  
Vol 15 (4) ◽  
pp. 789-794 ◽  
Author(s):  
G. Ramanaiah ◽  
V. Kumaran

The Darcy-Brinkman free convection near a wedge and a cone in a porous medium with high porosity has been considered. The surfaces are subjected to a mixed thermal boundary condition characterized by a parameterm;m=0,1,∞correspond to the cases of prescribed temperature, prescribed heat flux and prescribed heat transfer coefficient respectively. It is shown that the solutions for differentmare dependent and a transformation group has been found, through which one can get solution for anymprovided solution for a particular value ofmis known. The effects of Darcy number on skin friction and rate of heat transfer are analyzed.


Sign in / Sign up

Export Citation Format

Share Document