Calculation of conjugate heat transfer problem with volumetric heat generation using the Galerkin method (original) (raw)

The Galerkin method solution of the conjugate heat transfer problems for the cross-flow conditions

A conjugate heat transfer model of fluid flow across a solid heat conducting structure has been built. Two examples are presented: a.) air-stream cooling of the solid structure and b.) flow across rods with volumetric heat generation. To construct the model, a Volume Average Technique (VAT) has been applied to the momentum and the energy transport equations for a fluid and a solid phase to develop a specific form of porous media flow equations. The model equations have been solved with the semi-analytical Galerkin method. The detailed velocity and temperature fields in the fluid flow and the solid structure have been obtained. Using the solution fields, the whole-section drag coefficient Cd and the whole-section Nusselt number Nu have been also calculated. To validate the developed solution procedure, the results have been compared to the results of the finite volume method and to the experimental data. The comparison demonstrates an excellent agreement.

APPLICATION OF GALERKIN METHOD TO CONJUGATE HEAT TRANSFER CALCULATION

Numerical Heat Transfer, Part B: Fundamentals, 2003

A fast-running computational algorithm based on the volume averaging technique (VAT) is developed and solutions are obtained using the Galerkin method (GM). The goal is to extend applicability of the GM to the area of heat exchangers in order to provide a reliable benchmark for numerical calculations of conjugate heat transfer problems. Using the VAT, the computational algorithm is fast-running, but still able to present a detailed picture of temperature fields in air flow as well as in the solid structure of the heat sink. The calculated whole-section drag coefficient C d and Nusselt number Nu were compared with finite-volume method (FVM) results and with experimental data to verify the computational model. The comparison shows good agreement. The present results demonstrate that the selected Galerkin approach is capable to perform heat exchanger calculations where the thermal conductivity of the solid structure has to be taken into account.

Modeling of conjugate heat transfer using Galerkin approach

An algorithm for simulation of conjugate heat transfer in an electronic chip heat sink is described. Applying Volume Averaging Theory (VAT) to a system of transport equations, a heat exchanger structure is modeled as a homogeneous porous media. The interaction between the fluid and the structure, the VAT equation closure requirement, is accomplished with drag and heat transfer coefficients taken from the available literature and inserted into a computer code. The system of partial differential equations is solved using the Galerkin method to decompose the temperature field into a series of eigenfunctions. An example calculation is performed for an aluminum heat sink exposed to force convection airflow. The geometry of the simulation domain and boundary conditions follow the geometry of the experimental test section used in the Morrin-Martinelli-Gier Memorial Heat Transfer Laboratory at University of California, Los Angeles. A comparison of the whole-section drag coefficient and Nusselt number as functions of Reynolds number shows good agreement with finite volume method results as well as with experimental data. The calculated temperature fields reveal the local heat flow distribution and enable optimization of the surface geometry.

Finite Element Method for Analysis of Conjugate Heat Transfer between Solid and Unsteady Viscous Flow

Engineering Journal, 2009

A fractional four-step finite element method for analyzing conjugate heat transfer between solid and unsteady viscous flow is presented. The secondorder semi-implicit Crank-Nicolson scheme is used for time integration and the resulting nonlinear equations are linearized without losing the overall time accuracy. The streamline upwind Petrov-Galerkin method (SUPG) is applied for the weighted formulation of the Navier-Stokes equations. The method uses a three-node triangular element with equal-order interpolation functions for all the variables of the velocity components, the pressure and the temperature. The main advantage of the method presented is to consistently couple heat transfer along the fluid-solid interface. Four test cases, which are the lid-driven cavity flow, natural convection in a square cavity, transient flow over a heated circular cylinder and forced convection cooling across rectangular blocks, are selected to evaluate the efficiency of the method presented.

An immersed-boundary method for conjugate heat transfer analysis

Journal of Mechanical Science and Technology, 2017

An immersed-boundary method is proposed for the analysis of conjugate problems of convective heat transfer in conducting solids. Inside the solid body, momentum forcing is applied to set the velocity to zero. A thermal conductivity ratio and a heat capacity ratio, between the solid body and the fluid, are introduced so that the energy equation is reduced to the heat diffusion equation. At the solid fluid interface, an effective conductivity is introduced to satisfy the heat flux continuity. The effective thermal conductivity is obtained by considering the heat balance at the interface or by using a harmonic mean formulation. The method is first validated against the analytic solution to the heat transfer problem in a fully developed laminar channel flow with conducting solid walls. Then it is applied to a laminar channel flow with a heated, block-shaped obstacle to show its validity for geometry with sharp edges. Finally the validation for a curvilinear solid body is accomplished with a laminar flow through arrayed cylinders.

Evaluation of Two Finite Element Schemes for Conjugate Heat Transfer Problems

This paper presents two finite element schemes for solving conjugate heat transfer problems, where heat conduction in a solid is coupled with heat convection in viscous fluid flow. For solving viscous incompressible thermal flow in fluid region, the Streamline Upwind Finite Element method and the Streamline Upwind Petrov-Galerkin method are selected, while heat conduction in solid region is solved using the standard Galerkin method. The methods use the three-node triangular element with equal-order interpolation functions for all the variables of the velocity components, the pressure and the temperature. The main advantage of the presented approach is to consistently couple heat transfer along the solid-fluid interface. Three test cases, conjugate Couette flow problem in parallel plate channel, counter-flow in heat exchanger, and conjugate natural convection in a square cavity with a conducting wall, are selected to evaluate the presented algorithms.

Numerical Study of Heat Transfer Coefficient in Porous Media

2009

⎯ In this study, convective heat transfer in a porous flat plate channel flow is simulated by a direct numerical method. The solid materials consist of uniform distributed blocks, which resemble the porous medium within the channel. The solid materials are assumed to be isothermal and the channel walls are under adiabatic conditions. The Navier–Stokes equations are solved directly in the fluid region without the assumption of volume averaging. The two energy transport equations are solved for the solid and fluid flow separately. The results indicate that the mean bulk temperature across the channel develops faster if the channel aspect ratio gets smaller. On the other hand, the Nusselt number has the highest value at the channel inlet and gradually approaches to a minimum or developing condition at a distance which depends on the value of the aspect ratio.

Numerical Investigation of Conjugate Heat Transfer Problems

2005

Conjugate heat transfer (CHT) analysis has been carried out for laminar flow past flat plate and turbulent flow between parallel plates using a commercial CFD software CFXTASCFlow. Navier-Stokes equations alongwith k − ε turbulence model in the fluid and conduction equation in the solid have been solved simultaneously to obtain the flow features. The computed temperature distribution of the flow past flat plat matches very well with analytical and other numerical results. For the turbulent flow between parallel plates, nondimensional temperature distribution and Nusselt number distribution matches very well with the analytical and experimental results for different values of Reynolds number and thickness of the heat transfer plate.

Conjugate Heat Transfer Model Based on SIMPLE and Coupled Energy and Heat Equations

The International Journal of Multiphysics, 2021

In this study, a numerical weak coupling strategy for the modeling of a conjugate heat transfer phenomenon is considered. Where the incompressible Navier-Stokes equations are solved using the Semi-Implicit Method for Pressure Linked Equations (SIMPLE) as a first step, and then the heat conduction equation for solid is solved in a second step considering the convective velocity field resulting from the first step. A finite-difference approach is used for both discretized time and spatial operators. In this paper, a two-dimensional simulation case study of a steady uniform streamwise flow around heated rectangular and triangle solids is presented. The simulation is forward in time until the steady-state regime is reached as the residuals converge and tend to zero. The spatial analysis of the temperature is obtained through the numerical resolution of the incompressible Navier-Stokes energy equation and the heat diffusion equation for the fluid and solid media, respectively. The result...

Conjugate heat transfer inside a porous channel

Heat and Mass Transfer, 2005

Analytical and numerical analyses have been performed for fully developed forced convection in a fluid-saturated porous medium channel bounded by two parallel plates. The channel walls are assumed to be finite in thickness. Conduction heat transfer inside the channel wall is also accounted and the full problem is treated as a conjugate heat transfer problem. The flow in the porous material is described by the Darcy-Brinkman momentum equation. The outer surfaces of the solid walls are treated as isothermal. A temperature dependent volumetric heat generation is considered inside the solid wall only. Analytical expressions for velocity, temperature, and Nusselt number are obtained after simplifying and solving the governing differential equations with reasonable approximations. Subsequent results obtained by numerical calculations show an excellent agreement with the analytical results.