Fourier-Bessel Series Model for the Stefan Problem With Internal Heat Generation in Cylindrical Coordinates

2018 ◽  
Author(s):  
Lyudmyla Barannyk ◽  
John Crepeau ◽  
Patrick Paulus ◽  
Ali Siahpush
Keyword(s):  
Differential Equation ◽  
Steady State ◽  
Heat Generation ◽  
Numerical Solutions ◽  
Internal Heat ◽  
Internal Heat Generation ◽  
Time Dependent ◽  
Cylindrical Coordinates ◽  
System Approaches ◽  
Melting And Solidification

A nonlinear, first-order ordinary differential equation that involves Fourier-Bessel series terms has been derived to model the time-dependent motion of the solid-liquid interface during melting and solidification of a material with constant internal heat generation in cylindrical coordinates. The model is valid for all Stefan numbers. One of the primary applications of this problem is for a nuclear fuel rod during meltdown. The numerical solutions to this differential equation are compared to the solutions of a previously derived model that was based on the quasi-steady approximation, which is valid only for Stefan numbers less than one. The model presented in this paper contains exponentially decaying terms in the form of Fourier-Bessel series for the temperature gradients in both the solid and liquid phases. The agreement between the two models is excellent in the low Stefan number regime. For higher Stefan numbers, where the quasi-steady model is not accurate, the new model differs from the approximate model since it incorporates the time-dependent terms for small times, and as the system approaches steady-state, the curves converge. At higher Stefan numbers, the system approaches steady-state faster than for lower Stefan numbers. During the transient process for both melting and solidification, the temperature profiles become parabolic.

Integral Solutions of Phase Change With Internal Heat Generation

2004 ◽  
Author(s):  
Ali Siahpush ◽  
John Crepeau
Keyword(s):  
Differential Equation ◽  
Phase Change ◽  
Heat Generation ◽  
Numerical Solutions ◽  
Internal Heat ◽  
Melting Rate ◽  
Internal Heat Generation ◽  
Stefan Number ◽  
Temperature Boundary ◽  
Solid Liquid

This paper presents solutions to a one-dimensional solid-liquid phase change problem using the integral method for a semi-infinite material that generates internal heat. The analysis assumed a quadratic temperature profile and a constant temperature boundary condition on the exposed surface. We derived a differential equation for the solidification thickness as a function of the internal heat generation (IHG) and the Stefan number, which includes the temperature of the boundary. Plots of the numerical solutions for various values of the IHG and Stefan number show the time-dependant behavior of both the melting and solidification distances and rates. The IHG of the material opposes solidification and enhances melting. The differential equation shows that in steady-state, the thickness of the solidification band is inversely related to the square root of the IHG. The model also shows that the melting rate initially decreases and reaches a local minimum, then increases to an asymptotic value.

Hysteresis Heating of Railroad Bearing Thermoplastic Elastomer Suspension Element

2017 ◽  
Author(s):  
Oscar O. Rodriguez ◽  
Arturo A. Fuentes ◽  
Constantine Tarawneh ◽  
Robert E. Jones
Keyword(s):  
Finite Element Analysis ◽  
Finite Element ◽  
Steady State ◽  
Temperature Distribution ◽  
Heat Generation ◽  
Internal Heat ◽  
Thermoplastic Elastomer ◽  
Internal Heat Generation ◽  
Element Analysis ◽  
Railroad Bearing

Thermoplastic elastomers (TPE’s) are increasingly being used in rail service in load damping applications. They are superior to traditional elastomers primarily in their ease of fabrication. Like traditional elastomers they offer benefits including reduction in noise emissions and improved wear resistance in metal components that are in contact with such parts in the railcar suspension system. However, viscoelastic materials, such as the railroad bearing thermoplastic elastomer suspension element (or elastomeric pad), are known to develop self-heating (hysteresis) under cyclic loading, which can lead to undesirable consequences. Quantifying the hysteresis heating of the pad during operation is therefore essential to predict its dynamic response and structural integrity, as well as, to predict and understand the heat transfer paths from bearings into the truck assembly and other contacting components. This study investigates the internal heat generation in the suspension pad and its impact on the complete bearing assembly dynamics and thermal profile. Specifically, this paper presents an experimentally validated finite element thermal model of the elastomeric pad and its internal heat generation. The steady-state and transient-state temperature profiles produced by hysteresis heating of the elastomer pad are developed through a series of experiments and finite element analysis. The hysteresis heating is induced by the internal heat generation, which is a function of the loss modulus, strain, and frequency. Based on previous experimental studies, estimations of internally generated heat were obtained. The calculations show that the internal heat generation is impacted by temperature and frequency. At higher frequencies, the internally generated heat is significantly greater compared to lower frequencies, and at higher temperatures, the internally generated heat is significantly less compared to lower temperatures. However, during service operation, exposure of the suspension pad to higher loading frequencies above 10 Hz is less likely to occur. Therefore, internal heat generation values that have a significant impact on the suspension pad steady-state temperature are less likely to be reached. The commercial software package ALGOR 20.3TM is used to conduct the thermal finite element analysis. Different internal heating scenarios are simulated with the purpose of obtaining the bearing suspension element temperature distribution during normal and abnormal conditions. The results presented in this paper can be used in the future to acquire temperature distribution maps of complete bearing assemblies in service conditions and enable a refined model for the evolution of bearing temperature during operation.

Finite element thermal analysis of a moving porous fin with temperature-variant thermal conductivity and internal heat generation

2020 ◽  
Vol 1 (1) ◽  
pp. 110
Author(s):  
Gbeminiyi Sobamowo ◽  
Keyword(s):  
Heat Transfer ◽  
Thermal Conductivity ◽  
Finite Element ◽  
Heat Generation ◽  
Peclet Number ◽  
Numerical Solutions ◽  
Internal Heat ◽  
Internal Heat Generation ◽  
Péclet Number ◽  
Porous Fin

This paper focuses on finite element analysis of the thermal behaviour of a moving porous fin with temperature-variant thermal conductivity and internal heat generation. The numerical solutions are used to investigate the effects of Peclet number, Hartmann number, porous and convective parameters on the temperature distribution, heat transfer and efficiency of the moving fin. The results show that when the convective and porous parameters increase, the adimensional fin temperature decreases. However, the value of the fin temperature is amplified as the value Peclet number is enlarged. Also, an increase in the thermal conductivity and the internal heat generation cause the fin temperature to fall and the rate of heat transfer from the fin to decrease. Therefore, the operational parameters of the fin must be carefully selected to avoid thermal instability in the fin.

Analysis of Nonlinear Heat Transfer in Solids of Various Geometries with Variable Thermal Conductivity and Heat Source

2017 ◽  
Vol 377 ◽  
pp. 1-16
Author(s):  
Raseelo Joel Moitsheki ◽  
Oluwole Daniel Makinde
Keyword(s):  
Heat Transfer ◽  
Thermal Conductivity ◽  
Exact Solutions ◽  
Heat Generation ◽  
Numerical Solutions ◽  
Internal Heat ◽  
Variable Thermal Conductivity ◽  
Internal Heat Generation ◽  
Temperature Dependent ◽  
Power Law Function

In this paper we consider heat transfer in a hot body with different geometries. Here, the thermal conductivity and internal heat generation are both temperature-dependent. This assumption rendered the model considered to be nonlinear. We assume that thermal conductivity is given by a power law function. We employ the preliminary group classification to determine the cases of internal heat generation for which the principal Lie algebra extends by one. Exact solutions are constructed for the case when thermal conductivity is a differential consequence of internal heat generation term. We derive the approximate numerical solutions for the cases where exact solutions are difficult to construct or are nonexistent. The effects of parameters appearing in the model on temperature profile are studied.

On the Stefan Problem With Internal Heat Generation and Prescribed Heat Flux Conditions at the Boundary

2019 ◽  
Author(s):  
Lyudmyla Barannyk ◽  
Sidney D. V. Williams ◽  
Olufolahan Irene Ogidan ◽  
John C. Crepeau ◽  
Alexey Sakhnov
Keyword(s):  
Heat Flux ◽  
Steady State ◽  
Phase Change ◽  
Heat Generation ◽  
Separation Of Variables ◽  
Outer Boundary ◽  
Internal Heat ◽  
Internal Heat Generation ◽  
The Difference ◽  
Solid Liquid Interface

Abstract We study the evolution of the solid-liquid interface during melting and solidification of a material with constant internal heat generation and prescribed heat flux at the boundary for a plane wall and a cylinder. The equations are solved by splitting them into transient and steady-state components and then using separation of variables. This results in an ordinary differential equation for the interface that involves infinite series. The initial value problem is solved numerically, and solutions are compared to the previously published quasi-static solutions. We show that when the internal heat generation and the heat flux at the boundary are close in value to each other, the motion of the phase change front takes longer to reach steady-state than when the values are farther apart. As the difference between the internal heat generation and the heat flux increases, the transient solutions become more dominant and the numerical solution of the phase change front does not reach steady-state before the outer boundary or centerline is reached. The difference between the internal heat generation and the heat flux at the boundary can be used to control the motion and speed of the interface. The problem has applications for a nuclear fuel rod during meltdown.

Numerical Treatment of the MHD Convective Heat and Mass Transfer in an Electrically Conducting Fluid over an Infinite Solid Surface in Presence of the Internal Heat Generation

2003 ◽  
Vol 58 (11) ◽  
pp. 601-611 ◽  
Author(s):  
N. T. Eldabe ◽  
A. G. El-Sakka ◽  
Ashraf Fouad
Keyword(s):  
Mass Transfer ◽  
Heat Generation ◽  
Solid Surface ◽  
Numerical Solutions ◽  
Internal Heat ◽  
Internal Heat Generation ◽  
Conducting Fluid ◽  
Electrically Conducting ◽  
Linear Partial Differential Equations ◽  
Electrically Conducting Fluid

Numerical solutions of a set of non-linear partial differential equations are investigated. We obtained the velocity distribution of a conducting fluid flowing over an infinite solid surface in the presence of an uniform magnetic field and internal heat generation. The temperature and concentration distributions of the fluid are studied as well as the skin-friction, rate of mass transfer and local wall heat flux. The effect of the parameters of the problem on these distributions is illustrated graphically.

Analysis of Creep in Cylinders, Spheres and Thin Discs

1965 ◽  
Vol 7 (1) ◽  
pp. 82-92 ◽  
Author(s):  
E. M. Smith
Keyword(s):  
Steady State ◽  
Heat Generation ◽  
Time Variation ◽  
Creep Behaviour ◽  
Internal Heat ◽  
Internal Heat Generation ◽  
Steady State Creep ◽  
Time Interval ◽  
Simultaneous Solution ◽  
Short Period

A method of analysing the creep behaviour of cylinders, spheres and thin discs is described, in which the effects of time variation of the environmental conditions of temperature, pressure, rotational speed and internal heat generation may be accommodated. Deformation is considered as a series of steps each consisting of a short period of steady state creep followed by instantaneous stress readjustment to resatisfy the equations of equilibrium and compatibility. The relationships for stress readjustment can be written down once for all, since they are independent of the type of creep law used and of the length of time interval taken. The simultaneous solution of these relationships is presented in a form suitable for incorporation in computer programmes.

Sign in / Sign up

Export Citation Format

Share Document