Transient analysis of multilayered magneto-electro-thermoelastic strip due to nonuniform heat supply

This paper is concerned with the theoretical analysis of the multilayered magneto-electro-thermoelastic strip due to unsteady and nonuniform heat supply in the width direction. The transient two-dimensional temperature is analyzed by the methods of Laplace and finite sine transformations. We obtain the exact solution for the simply supported and multilayered magneto-electro-thermoelastic strip under a plane strain state. As an example, numerical calculations are carried out for a three-layered composite strip constructed of piezoelectric and magnetostrictive materials and the numerical results for temperature change, displacement, stress, electric potential, and magnetic potential distributions in a transient state are shown in figures and tables. The displacement and stress distributions are compared with those obtained for a similar thermoelastic strip. Furthermore the effects of the span-to-thickness ratio on the temperature change, displacements, stresses, electric potential, and magnetic potential are investigated.     

Thermal stresses in a thin circular disc of orthotropic material due to an instantaneous point heat source

Summary In this paper, we have analysed the transient plane thermal stress problem of a circular disc of orthotropic material with instantaneous point heat source. The variation of with time along different radius vectors is exhibited graphically and compared with that of the isotropic case.

---

Wärmespannungen in einer dünnen Kreisscheibe aus orthotropem Material zufolge einer punktförmigen instantanen Wärmequelle

Zusammenfassung In dieser Arbeit wurde das instationäre ebene Wärmespannungsproblem einer Kreisscheibe aus orthotropem Material zufolge einer punktförmigen instantanen Wärmequelle untersucht. Die Veränderung von über die Zeit für verschiedene Radien ist graphisch dargestellt und wird mit dem isotropen Fall verglichen.

---

Nomenclature r, polar coordinates - T temperature rise - ² ratio of conductivities - ² thermal diffusivity in -direction - J _n Bessel function ofn-th order - t time - p introduced in equation (2) - h heat transfer coefficient - a radius of circular disc - f(r, ) temperature distribution at initial state - r ₀, ₀ a point on the disc - T ₀ strength of point heat source - Dirac delta function - F stress function in two dimensions - ₁, ₂ coefficient of thermal expansion - a ₁₁,a ₁₂,a ₂₂,a ₆₆ elastic constants With 9 Figures     

A meshless analysis of three-dimensional transient heat conduction problems

In this paper, we consider a numerical modeling of a three-dimensional transient heat conduction problem. The modeling is carried out using a meshless reproducing kernel particle (RKPM) method. In the mathematical formulation, a variational method is employed to derive the discrete equations. The essential boundary conditions of the formulated problems are enforced by the penalty method. Compared with numerical methods based on meshes, the RKPM needs only scattered nodes, rather than having to mesh the domain of the problem. An error analysis of the RKPM for three-dimensional transient heat conduction problem is also presented in this paper. In order to demonstrate the applicability of the proposed solution procedures, numerical experiments are carried out for a few selected three-dimensional transient heat conduction problems.     

Multi-material topology optimization for the transient heat conduction problem using a sequential quadratic programming algorithm

Transient heat conduction analysis involves extensive computational cost. It becomes more serious for multi-material topology optimization, in which many design variables are involved and hundreds of iterations are usually required for convergence. This article aims to provide an efficient quadratic approximation for multi-material topology optimization of transient heat conduction problems. Reciprocal-type variables, instead of relative densities, are introduced as design variables. The sequential quadratic programming approach with explicit Hessians can be utilized as the optimizer for the computationally demanding optimization problem, by setting up a sequence of quadratic programs, in which the thermal compliance and weight can be explicitly approximated by the first and second order Taylor series expansion in terms of design variables. Numerical examples show clearly that the present approach can achieve better performance in terms of computational efficiency and iteration number than the solid isotropic material with penalization method solved by the commonly used method of moving asymptotes. In addition, a more lightweight design can be achieved by using multi-phase materials for the transient heat conductive problem, which demonstrates the necessity for multi-material topology optimization.     

Thermal performance of perforated plate matrix heat exchangers with effects from outer wall and flow channel geometry

The performance of high effectiveness (high NTU) perforated plate matrix heat exchangers (MHEs) is dependent on the geometry of the flow channels, as well as the longitudinal heat conduction through the outer walls. The effect of the above factors on the performance of MHEs is investigated in this paper numerically. The results obtained with the present model are validated with our own experimental results as well as those in the literature. The results show a strong influence of longitudinal heat conduction through the outer wall on the performance of MHEs. A parametric study has been carried out to arrive at the optimum flow channel geometry under given operating conditions.     

Real-time identification of a high-magnitude boundary heat flux on a plate

Identification of high-magnitude heat flux in real time is a challenging problem, since most of the currently available algorithms require large computation time in comparison with the time scale of the real physical problem. This paper presents a methodology that allows for quantifying the unsteady heat flux in real time by using the steady-state Kalman filter. Two different cases have been used to verify this algorithm and the estimates are in excellent agreement with the reference values.     

Boundary element method analyses of transient heat conduction in an unbounded solid layer containing inclusions

The boundary element method (BEM) is used to compute the three-dimensional transient heat conduction through an unbounded solid layer that may contain heterogeneities, when a pointwise heat source placed at some point in the media is excited. Analytical solutions for the steady-state response of this solid layer when subjected to a spatially sinusoidal harmonic heat line source are presented when the solid layer has no inclusions. These solutions are incorporated into a BEM formulation as Greens functions to avoid the discretization of flat media interfaces. The solution is obtained in the frequency domain, and time responses are computed by applying inverse (Fast) Fourier Transforms. Complex frequencies are used to prevent the aliasing phenomena. The results provided by the proposed Greens functions and BEM formulation are implemented and compared with those computed by a BEM code that uses the Greens functions for an unbounded media which requires the discretization of all solid interfaces with boundary elements. The proposed BEM model is then used to evaluate the temperature field evolution through an unbounded solid layer that contains cylindrical inclusions with different thermal properties, when illuminated by a plane heat source. In this model zero initial conditions are assumed. Different simulation analyses using this model are then performed to evaluate the importance of the thermal properties of the inclusions on transient heat conduction through the solid layer.     

Thermal Transport Properties of Layered Materials: Identification by a New Numerical Algorithm for Transient Measurements

Transient methods are widely used to determine thermal transport properties. In some situations they can be used for homogeneous media to measure several properties either simultaneously or separately. In this context an analytic model is available and a well-posed inverse problem of parameter identification has to be solved. The examination of composite media is more complicated. The algorithm proposed here allows simultaneous determination of the thermal conductivity and thermal diffusivity of layered dielectrics by transient measurements. It is based on a plane source that acts both as a resistive heater and temperature sensor. For the technique to be successful two essential aspects have to be considered: firstly, the mathematical modeling of the measured data (the forward problem) and secondly, the problem of ill-posedness of the inverse problem. For the proposed measurement configuration, a new fast data analysis algorithm based on an analytic solution for the forward problem is presented. In principle, a numerical solution such as an FEM solution of the heat conduction equation can be used instead of the analytical one, but the computational effort is much greater. The inverse problem is formulated as an output-least-squares problem, which leads to a transcendent algebraic system of equations. The method was successfully tested for different situations.Paper presented at the Fifteenth Symposium on Thermophysical Properties, June 22--27, 2003, Boulder, Colorado, U.S.A.     

The Dual-Phase-Lag Heat Conduction Model in Thin Slabs Under a Fluctuating Volumetric Thermal Disturbance

In the present work, the thermal behavior of a thin slab, under the effect of a fluctuating volumetric thermal disturbance described by the dual-phase-lag heat conduction model is investigated. It is found that the use of the dual-phase-lag heat conduction model is essential at large frequencies of the volumetric disturbance. It is found that the hyperbolic wave model deviates from the diffusion model when and the dual-phase-lag model deviates from the diffusion model when . where is the angular velocity of the fluctuating wall temperature, is the phase-lag in the heat flux vector and is the phase-lag in the temperature gradient vector.     

Thermal Behavior of a Stagnant Gas Confined in a Horizontal Microchannel as Described by the Dual-Phase-Lag Heat Conduction Model

The transient thermal behavior of a stagnant gas confined in a horizontal microchannel is investigated analytically under the effect of the dual-phase-lag heat conduction model. The microchannel is formed from two infinite horizontal parallel plates where the upper plate is heated isothermally and the lower one is kept adiabatic. The model that combines both the continuum approach and the possibility of slip at the boundary is adopted in this study. The effects of the Knudsen number Kn, the thermal relaxation time _q, and the thermal retardation time _T on the microchannel thermal behavior are investigated using three heat conduction models. It is found that the deviations between the predictions of the parabolic and the hyperbolic models are insignificant. On the other hand, the deviations between the parabolic and dual-phase-lag models are significant under the same operating conditions.     

Wave characteristic of femtosecond heat conduction in thin films

The hyperbolic heat conduction equation (HHC) is solved for submicrometer gold film irradiated with a short-pulse laser. The transient temperature profiles are calculated. It is shown that the solutions of HHC and standard heat diffusion equation are significantly different for submicrometer films.Paper presented at the Third Workshop on Subsecond Thermophysics, September 17–18, 1992, Graz, Austria.     

Model of Fractional Heat Conduction in a Thermoelastic Thin Slim Strip under Thermal Shock and Temperature-Dependent Thermal Conductivity

The present paper paper, we estimate the theory of thermoelasticity a thin slim strip under the variable thermal conductivity in the fractional-order form is solved. Thermal stress theory considering the equation of heat conduction based on the time-fractional derivative of Caputo of order α is applied to obtain a solution. We assumed that the strip surface is to be free from traction and impacted by a thermal shock. The transform of Laplace (LT) and numerical inversion techniques of Laplace were considered for solving the governing basic equations. The inverse of the LT was applied in a numerical manner considering the Fourier expansion technique. The numerical results for the physical variables were calculated numerically and displayed via graphs. The parameter of fractional order effect and variation of thermal conductivity on the displacement, stress, and temperature were investigated and compared with the results of previous studies. The results indicated the strong effect of the external parameters, especially the time-fractional derivative parameter on a thermoelastic thin slim strip phenomenon.     

Gradient method for inverse heat conduction problem in nanoscale

An inverse heat conduction problem for nanoscale structures was studied. The conduction phenomenon is modelled using the Boltzmann transport equation. Phonon‐mediated heat conduction in one dimension is considered. One boundary, where temperature observation takes place, is subject to a known boundary condition and the other boundary is exposed to an unknown temperature. The gradient method is employed to solve the described inverse problem. The sensitivity, adjoint and gradient equations are derived. Sample results are presented and discussed. Copyright © 2004 John Wiley & Sons, Ltd.     

Effects of wall thickness and material property on inverse heat conduction analysis of a hollow cylindrical tube

In this study, we examined the effects of a hollow cylindrical tube’s thickness and material properties on estimated time delay and waveform distortion in a one-dimensional inverse heat transfer analysis model using the thermal resistance method and an input estimation algorithm. Results indicated a persistent time delay for various heat flux amounts applied to different tube thicknesses. As the tube thickness increased, the numerically determined temperature data also experienced a time delay, which affected the inverse heat transfer response curve. Results also indicated that the transient heat flux waveform estimated for different material properties showed higher levels of distortion for materials having relatively low thermal conductivity. These materials also exhibited greater time delays. To address these issues, we applied a Fourier number (a dimensionless number representing the tube’s thickness and material properties) and proposed an equation to calculate sharpness, which can subsequently be used to predict probable time delays and heat flux waveform distortion. In conclusion, a correction is required when a low Fourier number is used in inverse heat transfer analysis.     

Deformation due to Mechanical and Thermal Sources in a Thermoelastic Body with Voids under Axi-symmetric Distributions

The Laplace and Hankel transforms have been employed to find the general solution of a homogeneous, isotropic, thermoelastic half-space with voids for a plane axi-symmetric problem. The application of a thermoelastic half-space with voids subjected to a normal force and a thermal source acting at the origin has been considered to show the utility of the solution obtained. To obtain the solution in a physical form, a numerical inversion technique has been applied. The results in the form of displacements, stresses, temperature distribution, and change in volume fraction field are computed numerically and illustrated graphically for a magnesium crystal-like material to depict the effects of voids in the theory of coupled thermoelasticity (CT) and uncoupled thermoelasticity (UCT) for an insulated boundary and a temperature gradient boundary.     

Flexural response to impulsive loads of a fluid-saturated thin poroelastic circular plate

Biot's poroelastic theory is used with classic plate theory and plane stress theory to determine the constitutive relationships for a thin poroelastic plate. The dynamic equations for the thin poroelastic plate are derived from the extended Hamilton's principle. The dynamic equations are then transformed to frequency domain and Galerkin's finite element method is used to derive the stiffness matrix of a triangular plate element. When impulsive loads and elastic boundary conditions are applied, the finite element frequency domain analysis for the thin poroelastic plates is achieved. Vibration behavior of thin elastic and poroelastic circular plates is accurately predicted.     

Wave propagation and transient response of a FGM plate under a point impact load based on higher-order shear deformation theory

In this paper, the wave propagation and transient response of an infinite functionally graded plate under a point impact load are presented. The effective material properties of functionally graded materials (FGMs) for the plate are assumed to vary continuously through the plate thickness and be distributed according to a volume fraction power law along the plate thickness. Based on the higher-order shear deformation theory and considering the effect of the rotary inertia, the governing equations of the wave propagation in the functionally graded plate are derived by using the Hamilton’s principle. The analytic dispersion relation of the functionally graded plate is obtained by means of integral transforms and a complete discussion of dispersion for the functionally graded plate is given. Then, using the dispersion relation and integral transforms, exact integral solutions for the functionally graded plate under a point impact load are obtained. The transient response curves of the functionally graded plates are plotted and the influence of volume fraction distributions on transient response of functionally graded plates is analyzed. Finally, the solutions of the higher-order shear deformation theory and the first-order shear deformation theory are studied.     

An infinitely large plate weakened by a circular-arc crack subjected to partially distributed loads

On the basis of the complex-variable approach for the first boundary condition problems, a mapping function is proposed to transform the contour surface of a circular arc crack into a unit circle. By this mapping, direct stress integration along the contour surface can be performed for the case when uniform tractions are applied on part of the crack edge. General complex stress functions are obtained by evaluating the Cauchy integral for the governing boundary equation. After the obtained stress functions are differentiated with respect to a reference angle in the mapped plane, the general complex stress functions for the circular-arc crack problem, when concentrated loads are applied on the crack surface, can be obtained. The importance of this solution lies in its general applicability to crack problems with arbitrary loading.