Abstract: Freeze and thaw occurs seasonally in river banks in northern countries. Little is known on how the riverbank soil temperature responds to air temperature changes and how freeze and thaw develops in a river bank seasonally. This study presents a two-dimensional heat conduction model for numerical investigations of seasonal freeze and thaw processes in an idealized river bank. The model uses the finite difference method and it is convenient for applications. The model is validated with an analytical solution and a field case with soil temperature distributions. It is then applied to the idealized river bank in terms of partially and fully saturated conditions with or without ice cover influence. Simulated results illustrate the response processes of the river bank to seasonal air temperature variations. It promotes the understanding of freeze and thaw processes in river banks and prepares for further investigation of frost and thaw impacts on riverbank stability.
Abstract: This article presents the application of the semi-analytic method (SAM) in the thermal management solution (TMS) of the energy storage system (ESS). The TMS studied in this work is fluid cooling. In fluid cooling, both effective heat conduction and heat convection are indispensable due to the heat transfer from solid to fluid. Correspondingly, an efficient TMS requires a design investigation of the following parameters: fluid inlet temperature, ESS initial temperature, fluid flow rate, working c rate, continuous working time, and materials properties. Their variation induces a change of thermal performance in the battery module, which is usually evaluated by numerical simulation. Compared to complicated computation resources and long computation time in simulation, the SAM is developed in this article to predict the thermal influence within a few seconds. In SAM, a fast prediction model is reckoned by combining numerical simulation with theoretical/empirical equations. The SAM can explore the thermal effect of boundary parameters in both steady-state and transient heat transfer scenarios within a short time. Therefore, the SAM developed in this work can simplify the design cycle of TMS and inspire more possibilities in TMS design.
Abstract: A two-step multigrid approach is proposed to solve the inverse heat conduction problem in a 3-D object under laser irradiation. In the first step, the location of the laser center is estimated using a coarse and uniform grid system. In the second step, the front-surface temperature is recovered in good accuracy using a multiple grid system in which fine mesh is used at laser spot center to capture the drastic temperature rise in this region but coarse mesh is employed in the peripheral region to reduce the total number of sensors required. The effectiveness of the two-step approach and the multiple grid system are demonstrated by the illustrative inverse solutions. If the measurement data for the temperature and heat flux on the back surface do not contain random error, the proposed multigrid approach can yield more accurate inverse solutions. When the back-surface measurement data contain random noise, accurate inverse solutions cannot be obtained if both temperature and heat flux are measured on the back surface.
Abstract: This paper discusses the implementation of the boundary element method (BEM) on an Excel spreadsheet and how it can be used in teaching vector calculus and simulation. There are two separate spreadheets, within which Laplace equation is solved by the BEM in two dimensions (LIBEM2) and axisymmetric three dimensions (LBEMA). The main algorithms are implemented in the associated programming language within Excel, Visual Basic for Applications (VBA). The BEM only requires a boundary mesh and hence it is a relatively accessible method. The BEM in the open spreadsheet environment is demonstrated as being useful as an aid to teaching and learning. The application of the BEM implemented on a spreadsheet for educational purposes in introductory vector calculus and simulation is explored. The development of assignment work is discussed, and sample results from student work are given. The spreadsheets were found to be useful tools in developing the students’ understanding of vector calculus and in simulating heat conduction.
Abstract: In the paper introduced the diagnostic technique making possible the research of internal structures in composite materials reinforced fibres using in different applications. The main reason of damages in structures of these materials is the changing distribution of load in constructions in the lifetime. Appearing defect is largely complicated because of the appearance of disturbing of continuity of reinforced fibres, binder cracks and loss of fibres adhesiveness from binders. Defect in composite materials is usually more complicated than in metals. At present, infrared thermography is the most effective method in non-destructive testing composite. One of IR thermography methods used in non-destructive evaluation is vibrothermography. The vibrothermography is not a new non-destructive method, but the new solution in this test is use ultrasonic waves to thermal stimulation of materials. In this paper, both modelling and experimental results which illustrate the advantages and limitations of ultrasonic IR thermography in inspecting composite materials will be presented. The ThermoSon computer program for computing 3D dynamic temperature distribuions in anisotropic layered solids with subsurface defects subject to ulrasonic stimulation was used to optimise heating parameters in the detection of subsurface defects in composite materials. The program allows for the analysis of transient heat conduction and ultrasonic wave propagation phenomena in solids. The experiments at MIAT were fulfilled by means of FLIR SC 7600 IR camera. Ultrasonic stimulation was performed with the frequency from 15 kHz to 30 kHz with maximum power up to 2 kW.
Abstract: In the current study, two-dimensional unsteady heat conduction in a functionally graded cylinder is studied analytically. The temperature distribution is in radial and longitudinal directions. Heat conduction coefficients are considered a power function of radius both in radial and longitudinal directions. The proposed solution can exactly satisfy the boundary conditions. Analytical unsteady temperature distribution for different parameters of functionally graded cylinder is investigated. The achieved exact solution is useful for thermal stress analysis of functionally graded cylinders. Regarding the analytical approach, this solution can be used to understand the concepts of heat conduction in functionally graded materials.
Abstract: This article presents a numerical method to find the
heat flux in an inhomogeneous inverse heat conduction problem with
linear boundary conditions and an extra specification at the terminal.
The method is based upon applying the satisfier function along with
the Ritz-Galerkin technique to reduce the approximate solution of the
inverse problem to the solution of a system of algebraic equations.
The instability of the problem is resolved by taking advantage of
the Landweber’s iterations as an admissible regularization strategy.
In computations, we find the stable and low-cost results which
demonstrate the efficiency of the technique.
Abstract: The present work numerically analyzes the transient heat transfer in the absorber plates of a flat-plate solar collector based on the dual-phase-lag (DPL) heat conduction model. An efficient numerical scheme involving the hybrid application of the Laplace transform and control volume methods is used to solve the linear hyperbolic heat conduction equation. This work also examines the effect of different medium parameters on the behavior of heat transfer. Results show that, while the heat-flux phase lag induces thermal waves in the medium, the temperature-gradient phase lag smoothens the thermal waves by promoting non-Fourier diffusion-like conduction into the medium.
Abstract: Bio-based carbon nanotubes (CNTs) have received considerable research attention due to their comparative advantages of high level stability, simplistic use, low toxicity and overall environmental friendliness. New potentials for improvement in heat transfer applications are presented due to their high aspect ratio, high thermal conductivity and special surface area. Phonons have been identified as being responsible for thermal conductivities in carbon nanotubes. Therefore, understanding the mechanism of heat conduction in CNTs involves investigating the difference between the varieties of phonon modes and knowing the kinds of phonon modes that play the dominant role. In this review, a reference to a different number of studies is made and in addition, the role of phonon relaxation rate mainly controlled by boundary scattering and three-phonon Umklapp scattering process was investigated. Results show that the phonon modes are sensitive to a number of nanotube conditions such as: diameter, length, temperature, defects and axial strain. At a low temperature (
Abstract: Numerical approach based on the electrical simulation method is proposed to solve a nonlinear transient heat conduction problem with nonlinear boundary for a spherical body. This problem represents a strong nonlinearity in both the governing equation for temperature dependent thermal property and the boundary condition for combined convective and radiative cooling. By analysing the equivalent electrical model using the electrical circuit simulation program HSPICE, transient temperature and heat flux distributions at sphere can be obtained easily and fast. The solutions clearly illustrate the effect of the radiation-conduction parameter Nrc, the Biot number and the linear coefficient of temperature dependent conductivity and heat capacity. On comparing the results with corresponding numerical solutions, the accuracy and efficiency of this computational method is found to be good.
Abstract: Transient simulation of the hydrogen-assisted self-ignition of propane-air mixtures were carried out in platinum-coated micro-channels from ambient cold-start conditions, using a two-dimensional model with reduced-order reaction schemes, heat conduction in the solid walls, convection and surface radiation heat transfer. The self-ignition behavior of hydrogen-propane mixed fuel is analyzed and compared with the heated feed case. Simulations indicate that hydrogen can successfully cause self-ignition of propane-air mixtures in catalytic micro-channels with a 0.2 mm gap size, eliminating the need for startup devices. The minimum hydrogen composition for propane self-ignition is found to be in the range of 0.8-2.8% (on a molar basis), and increases with increasing wall thermal conductivity, and decreasing inlet velocity or propane composition. Higher propane-air ratio results in earlier ignition. The ignition characteristics of hydrogen-assisted propane qualitatively resemble the selectively inlet feed preheating mode. Transient response of the mixed hydrogen- propane fuel reveals sequential ignition of propane followed by hydrogen. Front-end propane ignition is observed in all cases. Low wall thermal conductivities cause earlier ignition of the mixed hydrogen-propane fuel, subsequently resulting in low exit temperatures. The transient-state behavior of this micro-scale system is described, and the startup time and minimization of hydrogen usage are discussed.
Abstract: A solar receiver is designed for operation under
extremely uneven heat flux distribution, cyclic weather, and cloud
transient cycle conditions, which can include large thermal stress and
even receiver failure. In this study, the effect of different oil velocity
on convection coefficient factor and impact of wind velocity on local
Nusselt number by Finite Volume Method will be analyzed. This
study is organized to give an overview of the numerical modeling
using a MATLAB software, as an accurate, time efficient and
economical way of analyzing the heat transfer trends over stationary
receiver tube for different Reynolds number. The results reveal when
oil velocity is below 0.33m/s, the value of convection coefficient is
negligible at low temperature. The numerical graphs indicate that
when oil velocity increases up to 1.2 m/s, heat convection coefficient
increases significantly. In fact, a reduction in oil velocity causes a
reduction in heat conduction through the glass envelope. In addition,
the different local Nusselt number is reduced when the wind blows
toward the concave side of the collector and it has a significant effect
on heat losses reduction through the glass envelope.
Abstract: In this article, we used the residual correction method
to deal with transient thermoelastic problems with a hollow spherical
region when the continuum medium possesses spherically isotropic
thermoelastic properties. Based on linear thermoelastic theory, the
equations of hyperbolic heat conduction and thermoelastic motion
were combined to establish the thermoelastic dynamic model with
consideration of the deformation acceleration effect and non-Fourier
effect under the condition of transient thermal shock. The approximate
solutions of temperature and displacement distributions are obtained
using the residual correction method based on the maximum principle
in combination with the finite difference method, making it easier and
faster to obtain upper and lower approximations of exact solutions.
The proposed method is found to be an effective numerical method
with satisfactory accuracy. Moreover, the result shows that the effect
of transient thermal shock induced by deformation acceleration is
enhanced by non-Fourier heat conduction with increased peak stress.
The influence on the stress increases with the thermal relaxation time.
Abstract: The reheating furnace is used to reheat the steel slabs
before the hot-rolling process. The supported system includes the
stationary/moving beams, and the skid buttons which block some
thermal radiation transmitted to the bottom of the slabs. Therefore, it is
important to analyze the steel slab temperature distribution during the
heating period. A three-dimensional mathematical transient heat
transfer model for the prediction of temperature distribution within the
slab has been developed. The effects of different skid button height
(H=60mm, 90mm, and 120mm) and different gap distance between
two slabs (S=50mm, 75mm, and 100mm) on the slab skid mark
formation and temperature profiles are investigated. Comparison with
the in-situ experimental data from Steel Company in Taiwan shows
that the present heat transfer model works well for the prediction of
thermal behavior of the slab in the reheating furnace. It is found that
the skid mark severity decreases with an increase in the skid button
height. The effect of gap distance is important only for the slab edge
planes, while it is insignificant for the slab central planes.
Abstract: A solution methodology without using integral
transformation is proposed to develop analytical solutions for
transient heat conduction in nonuniform hollow cylinders with
time-dependent boundary condition at the outer surface. It is shown
that if the thermal conductivity and the specific heat of the medium
are in arbitrary polynomial function forms, the closed solutions of the
system can be developed. The influence of physical properties on the
temperature distribution of the system is studied. A numerical
example is given to illustrate the efficiency and the accuracy of the
solution methodology.
Abstract: In this paper, we introduced a gradient-based inverse
solver to obtain the missing boundary conditions based on the
readings of internal thermocouples. The results show that the method
is very sensitive to measurement errors, and becomes unstable when
small time steps are used. The artificial neural networks are shown to
be capable of capturing the whole thermal history on the run-out
table, but are not very effective in restoring the detailed behavior of
the boundary conditions. Also, they behave poorly in nonlinear cases
and where the boundary condition profile is different.
GA and PSO are more effective in finding a detailed
representation of the time-varying boundary conditions, as well as in
nonlinear cases. However, their convergence takes longer. A
variation of the basic PSO, called CRPSO, showed the best
performance among the three versions. Also, PSO proved to be
effective in handling noisy data, especially when its performance
parameters were tuned. An increase in the self-confidence parameter
was also found to be effective, as it increased the global search
capabilities of the algorithm. RPSO was the most effective variation
in dealing with noise, closely followed by CRPSO. The latter
variation is recommended for inverse heat conduction problems, as it
combines the efficiency and effectiveness required by these
problems.
Abstract: An alternative approach is proposed to develop the analytic solution for one dimensional heat conduction with one mixed type boundary condition and general time-dependent heat transfer coefficient. In this study, the physic meaning of the solution procedure is revealed. It is shown that the shifting function takes the physic meaning of the reciprocal of Biot function in the initial time. Numerical results show the accuracy of this study. Comparing with those given in the existing literature, the difference is less than 0.3%.
Abstract: This paper presents a complete dynamic modeling
of a membrane distillation process. The model contains two
consistent dynamic models. A 2D advection-diffusion equation
for modeling the whole process and a modified heat equation
for modeling the membrane itself. The complete model describes
the temperature diffusion phenomenon across the feed, membrane,
permeate containers and boundary layers of the membrane. It gives
an online and complete temperature profile for each point in the
domain. It explains heat conduction and convection mechanisms that
take place inside the process in terms of mathematical parameters, and
justify process behavior during transient and steady state phases. The
process is monitored for any sudden change in the performance at any
instance of time. In addition, it assists maintaining production rates
as desired, and gives recommendations during membrane fabrication
stages. System performance and parameters can be optimized
and controlled using this complete dynamic model. Evolution of
membrane boundary temperature with time, vapor mass transfer along
the process, and temperature difference between membrane boundary
layers are depicted and included. Simulations were performed over
the complete model with real membrane specifications. The plots
show consistency between 2D advection-diffusion model and the
expected behavior of the systems as well as literature. Evolution
of heat inside the membrane starting from transient response till
reaching steady state response for fixed and varying times is
illustrated.
Abstract: Owing to the lack of synchronization between the solar energy availability and the heat demands in a specific application, the energy storing sub-system is necessary to maintain the continuity of thermal process. The present work is dealing with an active solar heating storing system in which an air solar collector is connected to storing unit where this energy is distributed and provided to the heated space in a controlled manner. The solar collector is a box type absorber where the air flows between a number of vanes attached between the collector absorber and the bottom plate. This design can improve the efficiency due to increasing the heat transfer area exposed to the flowing air, as well as the heat conduction through the metal vanes from the top absorbing surface. The storing unit is a packed bed type where the air is coming from the air collector and circulated through the bed in order to add/remove the energy through the charging / discharging processes, respectively. The major advantage of the packed bed storage is its high degree of thermal stratification. Numerical solution of the packed bed energy storage is considered through dividing the bed into a number of equal segments for the bed particles and solved the energy equation for each segment depending on the neighbor ones. The studied design and performance parameters in the developed simulation model including, particle size, void fraction, etc. The final results showed that the collector efficiency was fluctuated between 55%-61% in winter season (January) under the climatic conditions of Misurata in Libya. Maximum temperature of 52ºC is attained at the top of the bed while the lower one is 25ºC at the end of the charging process of hot air into the bed. This distribution can satisfy the required load for the most house heating in Libya.
Abstract: A mathematical heat transfer model for the prediction of transient heating of the slab in a direct-fired walking beam type reheating furnace has been developed by considering the nongray thermal radiation with given furnace environments. The furnace is modeled as radiating nongray medium with carbon dioxide and water with five-zoned gas temperature and the furnace wall is considered as a constant temperature lower than furnace gas one. The slabs are moving with constant velocity depending on the residence time through the non-firing, charging, preheating, heating, and final soaking zones. Radiative heat flux obtained by considering the radiative heat exchange inside the furnace as well as convective one from the surrounding hot gases are introduced as boundary condition of the transient heat conduction within the slab. After validating thermal radiation model adopted in this work, thermal fields in both model and real reheating furnace are investigated in terms of radiative heat flux in the furnace and temperature inside the slab. The results show that the slab in the furnace can be more heated with higher slab emissivity and residence time.