Application of Fractional Derivatives in Thermal Analysis of Disk Brakes

This paper presents a Fractional Derivative Approach for thermal analysis of disk brakes. In this research, the problem is idealized as one-dimensional. The formulation developed contains fractional semi integral and derivative expressions, which provide an easy approach to compute friction surface temperature and heat flux as functions of time. Given the heat flux, the formulation provides a means to compute the surface temperature, and given the surface temperature, it provides a means to compute surface heat flux. A least square method is presented to smooth out the temperature curve and eliminate/reduce the effect of statistical variations in temperature due to measurement errors. It is shown that the integer power series approach to consider simple polynomials for least square purposes can lead to significant error. In contrast, the polynomials considered here contain fractional power terms. The formulation is extended to account for convective heat loss from the side surfaces. Using a simulated experiment, it is also shown that the present formulation predicts accurate values for the surface heat flux. Results of this study compare well with analytical and experimental results.     

Analytical solution of time fractional Cattaneo heat equation for finite slab under pulse heat flux

A fractional Cattaneo model is derived for studying the heat transfer in a finite slab irradiated by a short pulse laser. The analytical solutions for the fractional Cattaneo model, the classical Cattaneo-Vernotte model, and the Fourier model are obtained with finite Fourier and Laplace transforms. The effects of the fractional order parameter and the relaxation time on the temperature fields in the finite slab are investigated. The results show that the larger the fractional order parameter, the slower the thermal wave. Moreover, the higher the relaxation time, the slower the heat flux propagates. By comparing the fractional order Cattaneo model with the classical Cattaneo-Vernotte and Fourier models, it can be found that the heat flux predicted using the fractional Cattaneo model always transports from the high temperature to the low one, which is in accord with the second law of thermodynamics. However, the classical Cattaneo-Vernotte model shows that the unphysical heat flux sometimes transports from the low temperature to the high one.     

Fractional four-step finite element method for analysis of thermally coupled fluid-solid interaction problems

An integrated fluid-thermal-structural analysis approach is presented. In this approach, the heat conduction in a solid is coupled with the heat convection in the viscous flow of the fluid resulting in the thermal stress in the solid. The fractional four-step finite element method and the streamline upwind Petrov-Galerkin (SUPG) method are used to analyze the viscous thermal flow in the fluid. Analyses of the heat transfer and the thermal stress in the solid are performed by the Galerkin method. The second-order semiimplicit Crank-Nicolson scheme is used for the time integration. The resulting nonlinear equations are linearized to improve the computational efficiency. The integrated analysis method uses a three-node triangular element with equal-order interpolation functions for the fluid velocity components, the pressure, the temperature, and the solid displacements to simplify the overall finite element formulation. The main advantage of the present method is to consistently couple the heat transfer along the fluid-solid interface. Results of several tested problems show effiectiveness of the present finite element method, which provides insight into the integrated fluid-thermal-structural interaction phenomena.     

Approximate solutions to the Stefan problem with internal heat generation

Using a quasi-static approach valid for Stefan numbers less than one, we derive approximate equations governing the movement of a phase change front for materials which generate internal heat. These models are applied for both constant surface temperature and constant surface heat flux boundary conditions, in cylindrical, spherical, plane wall and semi-infinite geometries. Exact solutions with the constant surface temperature condition are obtained for the steady-state solidification thickness using the cylinder, sphere, and plane wall geometries which show that the thickness depends on the inverse square root of the internal heat generation. Under constant surface heat flux conditions, closed form equations can be obtained for the three geometries. In the case of the semi-infinite wall, we show that for constant temperature and constant heat flux out of the wall conditions, the solidification layer grows then remelts.     

ELASTO－PLASTICITY ANALYSIS BASED ON COLLOCATION WITH THE MOVING LEAST SQUARE METHOD

A meshless approach based on the moving least square method is developed for elasto-plasticity analysis, in which the incremental formulation is used. In this approach, the displacement shape functions are constructed by using the moving least square approximation, and the discrete governing equations for elasto-plastic material are constructed with the direct collocation method. The boundary conditions are also imposed by collocation. The method established is a truly meshless one, as it does not need any mesh, either for the purpose of interpolation of the solution variables, or for the purpose of construction of the discrete equations. It is simply formulated and very efficient, and no post-processing procedure is required to compute the derivatives of the unknown variables, since the solution from this method based on the moving least square approximation is already smooth enough. Numerical examples are given to verify the accuracy of the meshless method proposed for elasto-rdasticity analysis.     

Exact solution of heat transfer from a stretching surface with variable heat flux

An exact expression of the temperature distribution is constructed for the heat transfer from a stretching surface with prescribed power law heat flux. The stretching velocity is inversely proportional to the one third power of the distance measured along the surface from a thin slit. The final result is expressed in terms of hypergeometric functions. Although the exact solution is accomplished, some physically unrealistic phenomena are encounters for specific conditions. The temperature parameter which prescribe the surface heat flux, strongly affects those situations. Two types of temperature distribution are discussed: dimensionless temperatures with and without scaling to the dimensionless surface temperature. The expression of the temperature distribution without scaling is lucid to understand the heat transfer characteristics. Received on 23 July 1997     

Similarity analysis of axisymmetric free convection on a horizontal infinite plate subject to a mixed thermal boundary condition in a micropolar fluid

Similarity analysis of the problem of axisymmetric free convection on a horizontal infinite plate of a micropolar fluid is considered assuming that the plate is subjected to a mixed thermal boundary condition. It is shown that the thermal boundary condition is characterized by a positive parameter m and the two cases of m = 0 and m = 1 correspond to prescribed plate temperature and prescribed surface heat flux respectively. If one has to compute the heat transfer coefficient for various values of m, there is no need to solve the boundary value problem every time; it is enough to solve a certain polynomial equation provided the solution is known for any particular value of m. Received on 3 November 1998     

The effect of variable viscosity on the entropy generation and pumping power in a laminar fluid flow through a duct subjected to constant heat flux

Entropy generation and pumping power required for a laminar viscous flow in a duct subjected to constant heat flux has been investigated. The temperature dependence of the viscosity is taken into consideration. The ratio of pumping power to total heat flux decreases considerably and entropy generation increases along the duct length for viscous fluids. Therefore, it is shown that an optimum duct length may be obtained which minimizes total energy losses due to both entropy generation and pumping power. For low heat-flux conditions, entropy generation due to viscous friction becomes dominant and the dependence of viscosity on temperature must be considered in order to determine entropy generation accurately. Received on 17 May 1999     

单层微通道热沉拓扑结构对芯片中心热点的降温效能研究

通过热流固耦合模拟分析得到了不同微通道结构热沉基底的温度场及微通道内速度场,研究了相同入流功率下不同单层微通道拓扑结构对中心有高热流密度热点芯片的散热能力。结果表明:相同入流功率(0.05W)下,不同结构的散热能力排序由高到低为Y分形、弯曲散射、直散射(双侧出流)、直螺旋、直散射(单侧出流)、圆螺旋、树状分形、直槽结构;采用中心入流可有效降低芯片中心热点附近的温度,对于中心入流的散射结构,采用对称出流结构可提升其流动传热性能;Y分形结构具有良好的流动传热特性,对于热源面和中心热点均具有良好的散热效果。     

Thin Layer Thermography - a new heat transfer measurement technique

 In order to overcome the shortcomings of conventional temperature instrumentation to compute heat transfer rates, a novel technique for heat transfer measurement is presented. Stemming from infrared thermography, the potential for further development of the ‘Thin Layer Thermography’ is demonstrated. This new measurement technique is based on the wavelength dependent transmissivity of thin layers. It captures temperature distributions on the wall surface and simultaneously in a prescribed and well defined depth of the wall of a given object of interest. This enables the calculation of a temperature gradient normal to the surface and therefore the determination of the wall heat flux. Received: 8 April 1998/Accepted: 12 August 1998     

Performance enhancement of a DC-operated micropump with electroosmosis in a hybrid nanofluid: fractional Cattaneo heat flux problem

The purpose of this investigation is to theoretically shed some light on the effect of the unsteady electroosmotic flow (EOF) of an incompressible fractional second-grade fluid with low-dense mixtures of two spherical nanoparticles, copper, and titanium. The flow of the hybrid nanofluid takes place through a vertical micro-channel. A fractional Cattaneo model with heat conduction is considered. For the DC-operated micropump, the Lorentz force is responsible for the pressure difference through the microchannel. The Debye-Hükel approximation is utilized to linearize the charge density. The semi-analytical solutions for the velocity and heat equations are obtained with the Laplace and finite Fourier sine transforms and their numerical inverses. In addition to the analytical procedures, a numerical algorithm based on the finite difference method is introduced for the given domain. A comparison between the two solutions is presented. The variations of the velocity heat transfer against the enhancements in the pertinent parameters are thoroughly investigated graphically. It is noticed that the fractional-order parameter provides a crucial memory effect on the fluid and temperature fields. The present work has theoretical implications for biofluid-based microfluidic transport systems.

     

Planar solidification of a warm flowing-liquid under different boundary conditions

Planar solidification of a warm flowing liquid with the convective heat transfer to the growing solid layer, has been analysed for the boundary conditions of constant temperature, constant heat flux and convective heat flux at the surface respectively. The mathematical formulation of the problem resulted in a coupled set of two differential equations in temperature and solid thickness as function of position, time and the problem parameters. Analytical expressions for the temperature distribution within the growing solid layer, the rate of solidification and the solidification time are obtained. The perturbation techniques employed here is simple and straight forward in contrast with the earlier techniques. Good agreement between the experimental results and the present solutions is obtained for the convective heat flux boundary condition. The results of this analysis are useful in the design and analysis of experiments dealing with freezing/melting in one dimension. The role of the parameter Stefan number which is small for phase change materials, is discussed in context with the storage of thermal energy.     

Flow and heat transfer of an electrically conducting fluid of second grade in a porous medium over a stretching sheet subject to a transverse magnetic field

An analysis is performed for flow and heat transfer of a steady laminar boundary layer flow of an electrically conducting fluid of second grade in a porous medium subject to a transverse uniform magnetic field past a semi-infinite stretching sheet with power-law surface temperature or power-law surface heat flux. The effects of viscous dissipation, internal heat generation of absorption and work done due to deformation are considered in the energy equation. The variations of surface temperature gradient for the prescribed surface temperature case (PST) and surface temperature for the prescribed heat flux case (PHF) with various parameters are tabulated. The asymptotic expansions of the solutions for large Prandtl number are also given for the two heating conditions. It is shown that, when the Eckert number is large enough, the heat flow may transfer from the fluid to the wall rather than from the wall to the fluid when Eckert number is small. A physical explanation is given for this phenomenon.     

Study of the analytical solution to the heat transfer problem and surface temperature in a semi-infinite body with a constant heat flux at the surface and an initial temperature distribution

In many practical cases, one heats a semi-infinite solid with a constant heat flux source. For such an unsteady heat transfer problem, if the body has a uniform initial temperature, the analytical solution has been given by Carslaw and Jaeger. The surface temperature of the semi-infinite body follows the $\sqrt t $ -rule, that is, the surface temperature changes in proportion to square root of heating time. But if, instead of the uniform initial temperature, the body has a temperature distribution at the beginning of heating, the analytical solution has not yet been developed. Analytical solutions to the same problem with an exponential or a linear initial temperature distribution are obtained in this paper. It is shown, that in the case of a linear initial temperature distribution the surface temperature also changes according to $\sqrt t $ -rule Approximating the initial temperature distribution near the surface by its tangent at the surface, it is found that the surface temperature within a short time after the start of heating should also satisfy the $\sqrt t $ -rule, in spite of an arbitrary initial temperature distribution. The experimental data support this argument. Furthermore, the constant heat flux can be calculated after relationship between the surface temperature and heating time according to the equation derived in this paper, if the initial temperature distribution or its first-order derivative at the surface is known.     

Infrared temperature sensing of mechanically loaded specimens: Thermal analysis

Infrared temperature-sensing techniques have the major advantages of virtually no interference of the sensor with the sensed phenomenon and fast inherent response. On the other hand, infrared temperature sensing, as a superficial measurement technique, does not indicate the specimen's core temperatures, and hence a complementary thermal analysis is required. A thermal analysis of surface temperature measurements of a mechanically loaded cylindrical specimen is presented. The specimen is modeled as an infinite cylinder, suddenly exposed to a uniformly distributed volumetric heat source. This heat source results from the conversion of mechanical energy into thermal energy. A closed-form solution is obtained and numerical examples are given for metallic and polymeric specimens. The current analysis provides the upper boundaries for temperature differences between the core and the surface temperatures when compared with the actual problem of a finite specimen. It is shown that surface temperature measurement is a good indication of the core temperature for metallic specimens but may lead to some poor results in the case of polymeric specimens. It is found that the transient thermal response of the infinite cylinder to sudden heating behaves like a first-order process. In the case of cyclic loading, the typical time scale of loading is found to be at least two orders of magnitude shorter than the typical time scale of heat transfer. Hence, the specimen is affected by the average power of heat generation, and not by the instantaneous effect of heating within a single loading cycle.     

Simulation of liquid sloshing in curved‐wall containers with arbitrary Lagrangian–Eulerian method

There are many challenges in the numerical simulation of liquid sloshing in horizontal cylinders and spherical containers using the finite element method of arbitrary Lagrangian–Eulerian (ALE) formulation: tracking the motion of the free surface with the contact points, defining the mesh velocity on the curved wall boundary and updating the computational mesh. In order to keep the contact points slipping along the curved side wall, the shape vector in each time advancement is defined to modify the kinematical boundary conditions on the free surface. A special function is introduced to automatically smooth the nodal velocities on the curved wall boundary based on the liquid nodal velocities. The elliptic partial differential equation with Dirichlet boundary conditions can directly rezone the inner nodal velocities in more than a single freedom. The incremental fractional step method is introduced to solve the finite element liquid equations. The numerical results that stemmed from the algorithm show good agreement with experimental phenomena, which demonstrates that the ALE method provides an efficient computing scheme in moving curved wall boundaries. This method can be extended to 3D cases by improving the technique to compute the shape vector. Copyright © 2007 John Wiley & Sons, Ltd.     

Numerical investigation of fully developed turbulent fluid flow and heat transfer in a square duct

Three-dimensional fully developed turbulent fluid flow and heat transfer in a square duct are numerically investigated with the author's anisotropic low-Reynolds-number k-ε turbulence model. Special attenton has been given to the regions close to the wall and the corner, which are known to influence the characteristics of secondary flow a great deal. Hence, instead of the common wall function approach, the no-slip boundary condition at the wall is directly used. Velocity and temperature profiles are predicted for fully developed turbulent flows with constant wall temperature. The predicted variations of both local wall shear stress and local wall heat flux are shown to be in close agreement with available experimental data. The present paper also presents the budget of turbulent kinetic energy equation and the systematic evaluation for existing wall function forms. The commonly adopted wall function forms that are valid for two-dimensional flows are found to be inadequate for three-dimensional turbulent flows in a square duct.     

Melting or sublimation of solids with volume heat sources

Summary The problem of the phase change produced in semi-infinite media by volume heat sources is considered for quasi steady state conditions and where the molten or sublimated material is removed as soon as it is formed. A general formal solution for the temperature distribution and the velocity of the phase change surface is presented, and it is shown that for source functions of physical interest the velocity of the surface is the ratio of the total heat flux to the energy required to produce the phase change of a unit volume of material. Several examples are considered.     

Solution of non-linear inverse heat conduction problems using the method of lines

 Two space marching methods for solving the one-dimensional nonlinear inverse heat conduction problems are presented. The temperature-dependent thermal properties and the boundary condition on the accessible part of the boundary of the body are known. Additional temperature measurements in time are taken with a sensor located in an arbitrary position within the solid, and the objective is to determine the surface temperature and heat flux on the remaining part of the unspecified boundary. The methods have the advantage that time derivatives are not replaced by finite differences and the good accuracy of the method results from an appropriate approximation of the first time derivative using smoothing polynomials. The extension of the first method presented in this study to higher dimensions inverse heat conduction problems is straightforward. Received on 3 May 1999