Next Article in Journal
Construction of Disease-Symptom Knowledge Graph from Web-Board Documents
Previous Article in Journal
Visual Demands of Walking Are Reflected in Eye-Blink-Evoked EEG-Activity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analytical Solutions to Temperature Field in Various Relative-Scale Media Subjected to a Reciprocating Motion Point Heat Source

1
School of Mechanical Engineering, Southeast University, Nanjing 211189, China
2
School of Mechanical Engineering, Yangzhou University, Yangzhou 225127, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2022, 12(13), 6612; https://doi.org/10.3390/app12136612
Submission received: 19 May 2022 / Revised: 26 June 2022 / Accepted: 27 June 2022 / Published: 29 June 2022
(This article belongs to the Section Applied Thermal Engineering)

Abstract

:
To reveal the temperature rise evolution mechanism of isotropic media subjected to reciprocating motion constant-strength point heat source, various forms of analytical solutions are derived on the basis of differentiated relative scales, and non-dimensionalized parameters are designed to characterize the thermal distribution regularities by utilizing numerical calculations. Temperature rise curves of media subjected to a reciprocating motion point heat source allow similar quasi-steady-state characteristics to appear, which finally achieve a stable state, so that the values of surplus temperature oscillate around the constant time-average quantity. The time to reach quasi-steady state, the time-averaged quantity and the fluctuation amplitude of surplus temperature are comprehensively impacted by the dimensionless distance parameter γ, the convective heat transfer parameter ω and the velocity and travel parameter β. This work discusses influence rules of temperature evolution in various relative-scale media and further enriches the moving heat source theory.

1. Introduction

Studies on heat transfer in media subjected to moving heat sources have been motivated by a broad range of applications. The research directions mainly focus on the heat transfer in various relative-scale media subjected to heat source with different shapes and moving paths. Heat transfer characteristics may appear differentiated among media with the irregular shapes and discrete heating area [1]. Generally, the studies can be classified by relative scales including the infinite [2], the semi-infinite [3] and the finite [4] media. Zhao et al. [5] presented a mapped transient infinite element in the global coordinate system, of which heat transfer functions were derived to solve thermal problems in the infinite media. Sheng et al. [6] revealed the temperature evolution mechanism of infinite/finite-length cylindrical solids and helped successfully predict the temperature rise in the ball screw for engineering applications. Winczek et al. [7] proposed an analytical semi-infinite body model for temperature field description of multi-pass arc weld surfacing. Discussed in a rectangular coordinate system, when the sizes are negligible compared with these of the heat transfer media in all three directions, the heat source can be regarded as a point heat source [8]. Division of the line and the surface heat source [9] obeys the similar rule. In addition, studies defined the spherical [10], the circular [11] and the cylindrical [12] heat sources in their respective coordinate systems. Wang et al. [13] established the common point and linear heat source model to simulate the temperature fields of electron beam welding and predict the weld shapes. Azar et al. [14] utilized analytical and numerical approaches to model a gas metal arc welding process with the discretely distributed point heat source. According to the various heat generation laws, the temperature solutions to a moving point heat source can also be divided into periodic-type [15], gamma-type [16] variable power and so on. The temperature field distribution in media subjected to straight [17], alterable-direction [18] and specific included angle motion heat sources have also been investigated; however, few studies considering the reciprocating characteristics of heat source can be seen in our visible references.
For the specific media subjected to an excitation heat source with constant unidirectional velocity, the transient heat conduction equation is commonly utilized to derive the corresponding temperature-response function, of which the simplified solution has been demonstrated effective by transforming into a moving coordinate whose origin moving at the coincident velocity with the heat source [16]. Green’s function method is one of the most significant tools for applying field theory in different branches of physic systems [19]. In the light of heat sources with the reciprocating characteristics, Green’s functions are no doubt able to construct a suitable mathematical description for the temperature field distribution.
Green’s functions can be utilized to deal with all kinds of unsteady-state heat conduction equations with heat source, non-homogeneous boundary and specific initial conditions. Green’s function method is also known as pulse decomposition or point source superposition method, which decomposes the physical quantities into pulses according to the time and space instead of the spectrum. The decomposed active factors including sources, boundary and initial disturbances are uniformly regarded as point sources. Sources generate fields and the heat conduction equations formed as partial differential are used to describe the distribution and variation of temperature fields.
The early investigations mainly focused on heat generation laws, shapes of heat sources, boundary conditions, etc., barely touching on the reciprocating motion of heat sources. In this work, periodically reciprocating characteristics of point heat source has been taken into prior consideration. Based on Green’s function method, analytical solutions in the form of generalized incomplete Gamma functions to temperature field of various relative-scale media are derived. Non-dimensionalized parameters are put forward to characterize the temperature field distribution. Influences of the variables ξ, γ, κ, ω and the reduced parameter β on the temperature field are discussed.

2. Problem Description and Mathematical Modelling

2.1. Path Simplification

Inevitable acceleration and deceleration exists in the actual reciprocating motion process of a heat source. The accompanying quadratic function of time presents great difficulties in theoretical derivation [20]. The Laplace transform is one of the most effective means to eliminate the partial derivatives of time variables in differential equations [21]. Nevertheless, neglecting the stages of variable speed motion is another more simple and convenient method, and displacement curves of the practical and the simplified motion are exhibited in Figure 1. Relationship between practical and simplified displacements can be deducted as
S 1 S 2 S 2 = 1 1 + τ u / τ v
The two displacement curves trends to be close enough if the motion time τv of variable speed is short enough compared with that τu of uniform speed, so that the motion of point heat source is assumed as periodically reciprocating motion at the uniform velocity of U from the outset O. The period of motion is T. Then, the displacement function S(τ) can be expressed as
S ( τ ) = { U ( τ K T ) ,     K T τ ( K + 1 / 2 ) T U [ ( K + 1 ) T τ ] , ( K + 1 / 2 ) T < τ < ( K + 1 ) T
where K is the integer of the ratio of time to period, K = [τ/T]. Representation method of time-displacement in the following article obeys the calculation rule of S(τ).

2.2. Mathematical Modelling

The reciprocating motion means that the heat source starts to make straight movement from some point in the medium and return to the outset after a full motion cycle. Then, the heat source repeats this kind of reciprocating process. As the heat transfer medium is infinite, we assume the outset of heat source as origin point and the move path direction as x axis. Then, we construct the corresponding rectangular coordinate system to study the heat transfer problem. The schematic diagram of the physical model is exhibited as Figure 2.
Heat conduction happens in a homogeneous and isotropic medium of the infinite extension in all directions, of which the initial temperature is t. The point heat source starts periodically reciprocating motion along the x axis in a rectangular coordinate region (x, y, z). Introducing the surplus temperature θ = t(x, y, z, t) − t and heating rate operator gp of moving heat source, the general mathematical model of three-dimensional heat transfer problem subjected to a point heat source can be expressed as
2 θ x 2 + 2 θ y 2 + 2 θ z 2 + 1 k g p ( t ) δ ( x x ) = 1 α θ t ,             < x , y , z < +
with initial condition
θ = I ( x , y , z , t ) ,           t = 0 ,           < x , y , z < +
where δ function is the density representation of point source as a centrally distributed physical quantity, δ(x − x′) represents the space distribution of point heat resource in the x axis [22], m, k is the thermal conductivity, and W/m⋅K, k = αρc. α is the thermal diffusivity, m2/s.
Due to the periodical characteristics of point heat source, departure and back travels consume the same time, so that the integral limit 0~t can be split by the interval T/2. Solutions to temperature field can be calculated as
θ = ( 0 T 2 + T 2 T + T 3 T 2 + ( K 1 ) T ( K 1 2 ) T + ( K 1 2 ) T K T + K T t ) G g p ( τ ) d τ
Tallying up the reciprocating regular pattern, the integral form of surplus temperature can be generally expressed as
θ = ( N = 0 K 1 N T ( N + 1 2 ) T + N = 0 K 1 ( N + 1 2 ) T ( N + 1 ) T + K T t ) G g p ( τ ) d τ
where K = 1, 2, ··· and NT~(N + 1/2)T represents the time period of departure while (N + 1/2)T~(N + 1)T represents that of back tracking. KT~t represents the rest time after deducting all the complete cycles. G is Green’s function of the corresponding heat transfer problem.

2.3. Integral Form of Surplus Temperature

Substituting (tτ) for t, separated variable method has been utilized to work out the one-dimensional Green’s function in x direction formed as [23]
G ( x , t | x , τ ) = α k 1 4 π α ( t τ ) exp [ ( x x ) 2 4 α ( t τ ) ]
Green’s functions in the other two directions have the similar forms so that corresponding Green’s function of the multidimensional heat transfer problem in an infinite medium can be expressed as the product form of three one-dimensional problem solutions. It is worth noting that the coefficient (α/k) can only appear once, no matter whether the heat transfer problem is one or multidimensional [22]. So, the Green’s function of above 3D heat transfer problem can be expressed as
G ( x , y , z , t | x , y , z , τ ) = α k 1 [ 4 π α ( t τ ) ] 3 / 2 exp [ [ x S ( τ ) ] 2 + y 2 + z 2 4 α ( t τ ) ]
In our study, the heating power is a constant unchangeable with time so that operator gp(t) = Q. For a reciprocating motion heat source, x′ = S(τ) is introduced in calculation. Then, the integral form of surplus temperature can be calculated as
θ = Q ( 4 π α ) 3 / 2 α k ( N = 0 K 1 N T ( N + 1 2 ) T + N = 0 K 1 ( N + 1 2 ) T ( N + 1 ) T + K T t ) 1 ( t τ ) 3 / 2 exp [ [ x S ( τ ) ] 2 + y 2 + z 2 4 α ( t τ ) ] d τ
where Q is the rate of heat transfer from a point source [24], W. Locating point heat source at x = x′, the temperature distribution produced by point heat source can be represented as T(x) = (Q/ρc) × δ(xx′). This formula constructs a bridge between temperature and heat conversion.
The formula can be solved in three sections, introducing
θ 1 N = Q / ( ρ c ) ( 4 π α ) 3 / 2 N T ( N + 1 2 ) T 1 ( t τ ) 3 / 2 exp [ [ x S ( τ ) ] 2 + y 2 + z 2 4 α ( t τ ) ] d τ
Bringing in the displacement calculation formula of the first half cycle, the position variable (xx′) over time can be noted: xU(τNT) = x + UNTUt + U(tτ). Substituting variables as XN = x + UNTUt; RN2 = XN2 + y2 + z2 and φ = RN2/[4α(tτ)], then Equation (10) can be rewritten in the following form
θ 1 N = Q 4 k π 3 / 2 R N exp ( X N U 2 α ) × R N 2 / [ 4 α ( t N T ) ] R N 2 / { 4 α [ t ( N + 1 / 2 ) T ] } ϕ 1 2 exp ( ϕ U 2 R N 2 16 α 2 ϕ 1 ) d ϕ
The expression can be denoted by the generalized incomplete Gamma functions [12] in the form as
θ 1 N = Q 4 k π 3 / 2 R N exp ( X N U 2 α ) × [ Γ ( 1 2 , R N 2 4 α ( t N T ) ; U 2 R N 2 16 α 2 ) Γ ( 1 2 , R N 2 4 α [ t ( N + 1 / 2 ) T ] ; U 2 R N 2 16 α 2 ) ]
The generalized incomplete Gamma function is proved useful in the analytical study of heat conduction problems in an infinite medium [16], of which the general form can be expressed as
Γ ( α , t ; b ) = t x α 1 exp ( x b x 1 ) d x
Similarly, substituting variables as YN = x + UtU(N + 1)T; SN2 = YN2 + y2 + z2 and η = RN2/[4α(tτ)]. The second half cycle part can be expressed as
θ 2 N = Q 4 k π 3 / 2 S N exp ( Y N U 2 α ) × [ Γ ( 1 2 , S N 2 4 α [ t ( N + 1 / 2 ) T ] ; U 2 S N 2 16 α 2 ) Γ ( 1 2 , S N 2 4 α [ t ( N + 1 ) T ] ; U 2 S N 2 16 α 2 ) ]
In particular, the role of rest time should be discussed separately according to whether the half cycle is full or not.
If tKTT/2,
θ r = Q 4 k π 3 / 2 R K exp ( X K U 2 α ) × Γ ( 1 2 , R K 2 4 α ( t K T ) ; U 2 R K 2 16 α 2 )
If tKT > T/2,
θ r = θ 1 K + Q 4 k π 3 / 2 S K exp ( Y K U 2 α ) × Γ ( 1 2 , S K 2 4 α [ t ( K + 1 / 2 ) T ] ; U 2 S K 2 16 α 2 )
where θ 1 K , XK, YK are obtained from θ 1 N , XN, YN by changing N to K. Then, the complete integral form of temperature field can be expressed as
θ = N = 0 K 1 ( θ 1 N + θ 2 N ) + θ r
stipulating when K = 0, θ = θr.

3. Non-Dimensionalization and Parameter Discussions

3.1. Dimensionless Parameters Design and Non-Dimensionalization

The essence of dimension is expansion of the quantity coefficient according to basic quantity class. The different selection of basic quantity class, or even of families the unit system belongs to, may cause the expression form of dimensional equations to be changed. Non-dimensionalization means all the dimensions are constructed into one through appropriate variable substitutions. Thus, some or all units involved in the physical quantity equations are removed to ensure the form of the equations remains independent from unit system families. For convenient analysis, dimensionless parameters are designed by introducing the variables group listed as
  • Θ = 8 α k θ   π 3 / 2 Q U : dimensionless surplus temperature;
  • ξ = U x 2 α ; μ = U y 2 α ; ς = U z 2 α : dimensionless coordinate;
  • γ 2 = μ 2 + ς 2 : dimensionless distance away from the x axis;
  • κ = t T : dimensionless time;
  • β = U 2 T 2 α : dimensionless velocity and travel parameter.
Then, the non-dimensionalized forms of each parts can be expressed as
Θ 1 N = e [ ξ + β ( N κ ) ] [ ξ + β ( N κ ) ] 2 + γ 2 × [ Γ ( 1 2 , [ ξ + β ( N κ ) ] 2 + γ 2 2 β ( κ N ) ; [ ξ + β ( N κ ) ] 2 + γ 2 4 ) Γ ( 1 2 , [ ξ + β ( N κ ) ] 2 + γ 2 2 β ( κ N 1 / 2 ) ; [ ξ + β ( N κ ) ] 2 + γ 2 4 ) ]
Θ 2 N = e ξ + β ( κ N 1 ) [ ξ + β ( κ N 1 ) ] 2 + γ 2 × [ Γ ( 1 2 , [ ξ + β ( κ N 1 ) ] 2 + γ 2 2 β ( κ N 1 / 2 ) ; [ ξ + β ( κ N 1 ) ] 2 + γ 2 4 ) Γ ( 1 2 , [ ξ + β ( κ N 1 ) ] 2 + γ 2 2 β ( κ N 1 ) ; [ ξ + β ( κ N 1 ) ] 2 + γ 2 4 ) ]
If κ − K ≤ 1/2,
Θ r = e [ ξ + β ( K κ ) ] [ ξ + β ( K κ ) ] 2 + γ 2 [ Γ ( 1 2 , [ ξ + β ( K κ ) ] 2 + γ 2 2 β ( κ K ) ; [ ξ + β ( K κ ) ] 2 + γ 2 4 ) ]
If κ − K > 1/2,
Θ r = Θ 1 K + e ξ + β ( κ K 1 ) [ ξ + β ( κ K 1 ) ] 2 + γ 2 × [ Γ ( 1 2 , [ ξ + β ( κ K 1 ) ] 2 + γ 2 2 β ( κ K 1 / 2 ) ; [ ξ + β ( κ K 1 ) ] 2 + γ 2 4 ) ]
specifying when K = 0, Θ = Θr.

3.2. Discussions on Dimensionless Parameters

Non-dimensionalization helps appropriately reduce the number of variables in equations, reveal the mathematical essence of physical formulas, and facilitate the relative value comparisons of various items, so as to simplify the calculation and facilitate the analysis. This is guided by the generalized incomplete Gamma formula [13]
Γ ( 1 2 , x ; b ) = π 2 [ e 2 b erfc ( x b x ) + e 2 b erfc ( x + b x ) ]
where erfc is complementary error function. Numerical results with various values of dimensionless parameters are graphically discussed in this section.
To reveal the temperature evolution regularities over time, Figure 3a shows that the curves experience a rapid ascending process at the initial period of heat source effect and then gradually flatten out. The dimensionless surplus temperature finally rises to a stable fluctuation state over time. The curves perform the characteristics of almost unchangeable time-averaged quantity but periodically fluctuated instantaneous quantity. The surplus temperature fluctuation mainly results from periodic motion of the point heat source. Dimensionless parameter γ represents the distance away from x axis. Meanwhile, the smaller the γ values, the higher the time-averaged quantity and the more violent fluctuation the dimensionless surplus temperature develops, the less time used to reach a stable state. As moving along the x axis, the point heat source impacts much more on the regions nearby (smaller γ values) rather than those far away. The fluctuation amplitude A is defined as the maximum offset quantity compared to the time-average value during oscillation process. Figure 3b clarifies more clearly the relationship between the fluctuation amplitude and the distance parameter. With fixed β values, the smaller the distance away from the x axis, the larger fluctuation amplitude. The fluctuation amplitudes A over dimensionless velocity and travel parameter β with various γ values show a similar rule, which performs a rapid increase at first and then reaches a stable state when parameter β starts to increase from zero. However, β only affects the fluctuation amplitude within a certain range, in which the larger β values aggravate the oscillation of surplus temperature With smaller γ values, fluctuation amplitudes reach a stable state more rapidly and remain unchanged as parameter β increases. Moreover, the values of A are larger at a stable state, which indicates that the distance parameter is still the main influence factor of fluctuation amplitude compared with parameter β.
Figure 3c reveals the influence of reduced parameter β, the smaller β values, the greater time consumption for surplus temperature to reach a stable state and the lower final time-average value it ascends to. For a reciprocating moving heat source with a fixed travel, the lower velocity means larger time intervals between twice heating effects, and untimely heat compensation for that transferred to low-temperature regions causes the difficulty in temperature rise. Additionally, the value of β is inversely proportional to that of the thermal diffusivity. For regions in the infinite medium with inapparent temperature fluctuation, a lower β value means a faster heat transfer of media, so that the surplus temperature increases more slowly due to the increased heat dissipation. Moreover, more time is needed to reach a steady state, and its steady-state value is higher.
Figure 3d displays the evolution process of dimensionless surplus temperature. For suitable β and γ values, the received heat is more than the dissipated. Because of the reciprocating motion, for each value the point heat source passes by, it will cause the accumulation of residual heat, so that the surplus temperature continuously increases over time. The ξ-Θ curves have a peak value and symmetry about some point on the horizontal coordinate ξ. As the point heat source periodically moves between (0, 0, 0) and (UT/2, 0, 0) on the x axis, the most affected point is located at (UT/2), where the influence of the heat source continues for the longest time, and surplus temperature values are at the maximum. The point location at the converted horizontal coordinate is ξ = β/4. What needs to be particularly indicated is that the heating conditions are not symmetrical about ξ = β/4 at the early stage of the heat source motion, which is also the case for the temperature distribution. Only when the reciprocating time is enough, ξ = β/4 is flanked by almost the same heating conditions, and the Θ curve is considered symmetrical.
Figure 3e further confirms the symmetry as peak points of the fluctuation amplitude a and b, d and e are symmetrical about ξ = β/4, respectively. The most violent oscillation appears in the adjacent regions outside the abscissa interval [0, β/2] and then decreases quickly to almost zero. Due to shorter heating intervals, the fluctuation amplitude curves of surplus temperature reach a local trough of the wave at the symmetrical points, which indicates that, on the contrary, oscillation at the position of the heat source motion center is not violent. Fluctuation curves in the regions too far away from the symmetric center show more smooth features for the opposite reason. The regions of maximum surplus temperature and most violent fluctuation amplitude are abhorrent.

4. Analytical Solutions of Semi-Infinite Media

Assuming the boundary surface xoy of a specific semi-infinite medium, when t > 0, the medium conducts convective heat transfer with ambiance on the boundary z = 0. The initial temperature is set ti, while the convective heat transfer coefficient is h. Rate of heat transfer from a point heat source valued Q starts from the original point o to perform a periodically reciprocating motion along the positive direction of the x axis with the velocity U and the period T. A graphical description of the heat transfer problem is shown in Figure 4.
The mathematical model of the heat conduction problem in semi-infinite media can be expressed as
2 θ x 2 + 2 θ y 2 + 2 θ z 2 + 1 k g p ( t ) δ ( x x ) = 1 α θ t ,       < x , y < + ,   0 z < +
The boundary condition is
k θ z + h θ = 0 ,       z = 0 ,       t > 0
The initial condition is
θ = I ( z , t ) ,           t = 0 ,           0 z < +

4.1. Solutions of Green’s Function

The solution process of the Green’s function is similar with that of the heat transfer problem in an infinite medium. In fact, the three-dimensional unsteady-state heat transfer process, in a semi-infinite medium with the third kind of boundary condition, can been seen as the superposition of three one-dimensional problems. The Green’s function finally appears in the form of the product of three single functions.
G ( x , t | x , τ ) = α k 1 4 π α ( t τ ) exp [ ( x x ) 2 4 α ( t τ ) ]
G ( y , t | y , τ ) = 1 4 π α ( t τ ) exp [ ( y y ) 2 4 α ( t τ ) ]
G ( z , t | z , τ ) = [ 1 4 π α ( t τ ) { exp [ ( z z ) 2 4 α ( t τ ) ] + exp [ ( z + z ) 2 4 α ( t τ ) ] } H e H ( z + z ) + H 2 α ( t τ ) erfc ( z + z 4 α ( t τ ) + H α ( t τ ) ) ]
where H = h/k. So, for the r = (x, y, z) in semi-infinite media, Green’s function can be calculated as
G ( r , t | r , τ ) = G ( x , t | x , τ ) G ( y , t | y , τ ) G ( z , t | z , τ )

4.2. Integral Form of Surplus Temperature

For a moving inner point heat source, the integral form of surplus temperature θ = t(r,τ) − t should be calculated as [22,23]
θ ( r , t ) = k α R G | τ = 0 F ( r ) d v + 0 τ G ( r , t | r , τ ) g p ( r , τ )   d   τ
where F(r′) is the initial surplus temperature distribution and F(r′) = tit. R is the whole integral domain. In particular, the ambient surplus temperature can be expressed as [22]
θ amb = θ 0 R G ( x , y , z , t | x , y , z , 0 )   d x   d y   d z                               = θ 0 [ erf     ( z 4 α t ) + e H z + H 2 α t erfc     ( z 4 α t + H α t ) ]
where erf is the error function.
The calculation of the second half obeys the rule of Equation (4), and the resulting expressions are given directly to avoid repeats and save article space. So,
θ 12 N = Q H e H z X N U 2 α 2 k π [ χ ( 0 , z 4 α ( t N T ) ; X N 2 + y 2 x 2 , z 2 ( U 2 4 α 2 H 2 ) 16 α 2 , H z 2 ) χ ( 0 , z 4 α [ t ( N + 1 / 2 ) T ] ; X N 2 + y 2 x 2 , z 2 ( U 2 4 α 2 H 2 ) 16 α 2 , H z 2 ) ]
where χ(α, t; a, b, c) is the custom expression which cannot be defined by the ready-made functions. The form of χ can be expressed as
χ ( α , t ; a , b , c ) = t x α 1 exp ( a x 2 b x 2 ) erfc ( x + c x 1 ) d x
Similarly,
θ 22 N = Q H e H z Y N U 2 α 2 k π [ χ ( 0 , z 4 α [ t ( N + 1 / 2 ) T ] ; Y N 2 + y 2 x 2 , z 2 ( U 2 4 α 2 H 2 ) 16 α 2 , H z 2 ) χ ( 0 , z 4 α [ t ( N + 1 ) T ] ; Y N 2 + y 2 x 2 , z 2 ( U 2 4 α 2 H 2 ) 16 α 2 , H z 2 ) ]
If t − KTT/2
θ r 2 = Q H e H z X K U 2 α 2 k π [ χ ( 0 , z 4 α ( t K T ) ; X K 2 + y 2 x 2 , z 2 ( U 2 4 α 2 H 2 ) 16 α 2 , H z 2 ) ]
If tKT > T/2,
θ r 2 = θ 12 K + Q H e H z Y K U 2 α 2 k π [ χ ( 0 , z 4 α [ t ( N + 1 / 2 )     T ] ; Y K 2 + y 2 x 2 , z 2 ( U 2 4 α 2 H 2 ) 16 α 2 , H z 2 ) ]
where θ 12 K , XK, YK are obtained from θ 1 N , XN, YN by changing N to K. Then, the complete integral form of temperature field can be expressed as
θ ( x , y , z , t ) = θ amb + N = 0 K 1 [ ( θ 11 N θ 12 N ) ( θ 21 N θ 22 N ) ] + ( θ r 1 θ r 2 )
where θ 11 N , θ 21 N and θ r 1 are the calculation expressions of the parts which satisfied ready-made functions in the integral solutions and they have the same form with θ 1 N , θ 2 N and θ r in Equation (17).

4.3. Non-Dimensionalization and Parameter Discussion

By introducing one more dimensionless parameter ω = (2αH)/U, the non-dimensionalized forms of each parts can be expressed as
Θ amb = Θ 0 [ erf ( ς 2 κ β ) + e ( ω ς + κ β ς 2 2 ) erfc ( ς 2 2 κ β + ω κ β 2 ) ]
Θ 12 N = π ω e ω ς ξ β ( N κ ) × [ χ ( 0 , ς 2 β ( κ N ) ; [ ξ + β ( κ N ) ] 2 + μ 2 ς 2 , ( 1 ω 2 ) ς 2 4 , ω ς 2 ) χ ( 0 , ς 2 β [ κ ( N + 1 / 2 ) ] ; [ ξ + β ( κ N ) ] 2 + μ 2 ς 2 , ( 1 ω 2 ) ς 2 4 , ω ς 2 ) ]
Θ 22 N = π ω e ω ς ξ β ( N κ ) × [ χ ( 0 , ς 2 β [ κ ( N + 1 / 2 ) ] ; [ ξ + β [ κ ( N + 1 ) ] ] 2 + μ 2 ς 2 , ( 1 ω 2 ) ς 2 4 , ω ς 2 ) χ ( 0 , ς 2 β [ κ ( N + 1 ) ] ; [ ξ + [ κ ( N + 1 ) ] ] 2 + μ 2 ς 2 , ( 1 ω 2 ) ς 2 4 , ω ς 2 ) ]
If κ − K ≤ 1/2,
Θ r 2 = π ω e ω ς ξ β ( K κ ) [ χ ( 0 , ς 2 β ( κ K ) ; [ ξ + β ( κ K ) ] 2 + μ 2 ς 2 , ( 1 ω 2 ) ς 2 4 , ω ς 2 ) ]
If κ − K > 1/2,
Θ r 2 = Θ 12 K + π ω e ω ς ξ β ( K κ ) × [ χ ( 0 , ς 2 β [ κ ( K + 1 / 2 ) ] ; [ ξ + β [ κ ( K + 1 ) ] ] 2 + μ 2 ς 2 , ( 1 ω 2 ) ς 2 4 , ω ς 2 ) ]
specifying when K = 0, Θ = Θamb + (Θr1 − Θr2).
It is easy to see that β characterizes the relative speed and motion range of the heat source, ω characterizes the relative intensity of boundary convective heat transfer, and Θ0 characterizes the difference between initial and ambient temperature. These three comprehensively describe the influence factors of temperature field from various angles.
Finally, two points are discussed about the temperature field.
(1)
If the initial temperature is equal to the ambient temperature, which means Θamb = 0, then
Θ = N = 0 K 1 [ ( Θ 11 N Θ 12 N ) + ( Θ 21 N Θ 22 N ) ] + ( Θ r 1 Θ r 2 )
Hence, the initial temperature only has an influence on the component Θamb of Θ but has nothing to do with the other components.
(2)
The situation of ω = 0.
With the two conditions listed as
h = 0, which means the adiabatic boundary. ω = 0 leads Θ 12 N = Θ 22 N = Θ r 2 = 0.
② Neglecting the influence of initial temperature, Θamb = 0.
Then, the expression of dimensionless surplus temperature degenerates into the form as
Θ = N = 0 K 1 ( Θ 11 N + Θ 21 N ) + Θ r 1
The equation is the surplus temperature solution in the region 0 < z < +∞ of an infinite medium induced by a 2Q-strength point heat source moving along x axis. Therefore, if Θamb = 0, it can be inferred that the temperature field of a semi-infinite medium subjected to a periodic-motion point heat source is similar to that in an infinite medium. However, the former temperature values are lower than the latter, due to surface cooling.

5. Conclusions

In this study, analytical solutions to the various relative-scale media are derived, including infinite and semi-infinite ones. Dimensionless parameters are designed to reveal the temperature evolution regularities from the heating mechanism. Some conclusions are made below:
  • The temperature field evolution is the comprehensive result of dimensionless parameters γ, β, ω and Θ0. In general, the smaller the distance away from the heat source, the shorter the time intervals brought by larger velocity, the lower the relative intensity of boundary-convective heat transfer and the higher the initial temperature, the greater the benefit will be to heat accumulation and temperature rise. Among which, the most important influence factor can be γ, as the regions far away from the point heat source are little impacted.
  • The periodically reciprocating motion of point heat source results in surplus temperature appearing as a special feature, oscillating around the stable time-averaged quantity. For an infinite medium, the surplus is symmetrical about the horizontal coordinate at ξ = β/4. The maximum temperature value is obtained at the same location, but the peak value of fluctuation amplitude appears outside the abscissa interval [0, β/2], the regions of maximum surplus temperature and most violent fluctuation amplitude are abhorrent.
  • The reduced parameter β has a critical influence on temperature distribution, the amplitude of temperature fluctuation, the time to reach steady-state and the stable time-averaged quantity. In the regions of an infinite medium where temperature fluctuation is not obvious, surplus temperature reaches a steady state more quickly, and the stable time-averaged quantity is larger when the reduced parameter β increases.
  • Only if some specific conditions are satisfied, such as adiabatic boundary and neglecting the influence of the initial temperature, the analytical solutions to the temperature field of the semi-infinite media degenerate into a similar form to these infinite media.
Studies on the temperature evolution mechanism help further enrich the moving heat source theory and have great potential to provide theoretical guidance for engineering practice, such as in the fields of welding, heat exchangers and so on.

Author Contributions

J.Z. and X.L. proposed the idea of this work. X.S. made main contributions to the derivation of specific formulas and he wrote numerical algorithm programs in Matlab with the help of Y.X., D.H. and Y.L. are responsible for the edition and analysis of the numerical results. X.L. finally reviewed and checked the whole article. All authors have read and agreed to the published version of the manuscript.

Funding

This research was founded by National Key Research and Development Plan, grant number No. 2019YFB2006402.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

X.S. wants to express his thanks to his wife and lovely daughter for their financial and spiritual support during the study for his doctorate. Your selfless help is his motivation for constantly researching in this field.

Conflicts of Interest

The author(s) declared no potential conflict of interest with respect to the research, authorship, and/or publication of this article. The authors declare that they have no conflict of interest.

Glossary

Nomenclature
Adimensionless fluctuation amplitude, (-)
a, b, cgeneral expression of variables, (-)
cspecific heat capacity of mass, (J/kg⋅K)
gheating rate operator, (-)
Hprescribed form, H = h/k, (-)
hconvective heat transfer coefficient, (W/m2⋅°C)
I(*)initial condition, (-)
Kinteger of the ratio of time to period, (-)
kthermal conductivity coefficient, (W/m⋅K)
Nsequence number of complete cycles, (-)
Oorigin coordinate, (-)
QHeat power, (W)
rposition vector, (-)
Sdisplacement, (m)
Tmotion period, (s);temperature distribution, (°C)
ttemperature, (°C); time, (s)
Uuniform motion velocity, (m/s)
x:y,zcoordinate position, (m)
Greek symbol
αthermal diffusivity, (m2/s)
θsurplus temperature, (°C)
ρmedium density, (kg/m3)
τtime, (s)
δ(xx′)distribution of physical quantities in space, (m)
Superscripts
Kthe former k complete periods
Nsequence number of complete cycles
Subscripts
1simplified situation;former half of the kth complete period
2practical situation;latter half of the kth complete period
ppoint heat source
uuniform velocity motion
vvariable velocity motion
environment

References

  1. Cui, C.; Huang, X.Y.; Liu, C.Y. The Green function and its application to heat transfer in a low permeability porous channel. J. Electron. Packag. 2000, 122, 274–278. [Google Scholar] [CrossRef]
  2. Povstenko, Y.; Kyrylych, T. Time-fractional heat conduction in an infinite plane containing an external crack under heat flux loading. Comput. Math. Appl. 2019, 78, 1386–1395. [Google Scholar] [CrossRef]
  3. Berrone, L.R.; Tarzia, D.A.; Villa, L.T. Asymptotic behaviour of a non-classical heat conduction problem for a semi-infinite material. Math. Method Appl. Sci. 2000, 23, 1161–1177. [Google Scholar] [CrossRef]
  4. Gembarovic, J.; Gembarovic, J., Jr. A comparison of DHW algorithm for temperature distribution calculation with Fourier’s algorithm for transmission of heat between discrete bodies. Int. J. Thermophys. 2007, 28, 891–905. [Google Scholar] [CrossRef]
  5. Zhao, C.B.; Valliappan, S. Mapped transient infinite elements for heat transfer problems in infinite media. Comput. Methods Appl. Mech. Eng. 1993, 108, 119–131. [Google Scholar] [CrossRef]
  6. Sheng, X.; Zhang, J.R.; Lu, Y.Q.; Lu, X. Temperature evolution on infinite/finite- length cylindrical solids subjected to reciprocating motion heat source. Case Stud. Therm. Eng. 2021, 28, 101559. [Google Scholar] [CrossRef]
  7. Winczek, J. Modeling of Temperature field during multi-pass GMAW surfacing or rebuilding of steel elements taking into account the heat of the deposit metal. Appl. Sci. 2017, 7, 6. [Google Scholar] [CrossRef] [Green Version]
  8. Winczek, J.; Gawronska, E.; Gucwa, M.; Sczygiol, N. Theoretical and experimental investigation of temperature and phase transformation during SAW overlaying. Appl. Sci. 2019, 9, 1472. [Google Scholar] [CrossRef] [Green Version]
  9. Levin, P. A general solution of 3-D quasi-steady-state problem of a moving heat source on a semi-infinite solid. Mech. Res. Commun. 2008, 35, 151–157. [Google Scholar] [CrossRef]
  10. Zubair, S.M.; Chaudhry, M.A. Temperature solutions due to continuously operating spherical surface-heat sources in an infinite medium. Int. Commun. Heat Mass Transf. 1991, 18, 805–811. [Google Scholar] [CrossRef]
  11. Mebarek-Oudina, F. Numerical modeling of the hydrodynamic stability in vertical annulus with heat source of different lengths. Eng. Sci. Technol. 2017, 20, 1324–1333. [Google Scholar] [CrossRef]
  12. BniLam, N.; Al-Khoury, R. Transient heat conduction in an infinite medium subjected to multiple cylindrical heat sources: An application to shallow geothermal systems. Renew. Energ. 2016, 97, 145–154. [Google Scholar] [CrossRef]
  13. Wang, Y.J.; Fu, P.F.; Guan, Y.J.; Lu, Z.J.; Wei, Y.T. Research on modeling of heat source for electron beam welding fusion-solidification zone. Chin. J. Aeronaut. 2013, 26, 217–223. [Google Scholar] [CrossRef] [Green Version]
  14. Azar, A.S.; As, S.K.; Akselsen, O.M. Analytical modeling of weld bead shape in dry hyperbaric GMAW using Ar-He chamber gas mixtures. J. Mater. Eng. Perform. 2013, 22, 673–680. [Google Scholar] [CrossRef]
  15. Zubair, S.M.; Chaudhry, M.A. Temperature solutions due to steady, periodic-type, moving-point-heat sources in an infinite medium. Int. Commun. Heat Mass Transf. 1994, 21, 207–215. [Google Scholar] [CrossRef]
  16. Chaudhry, M.A.; Zubair, S.M. Analytic study of temperature solutions due to gamma-type moving point-heat source. Int. J. Heat Mass Transf. 1993, 36, 1633–1637. [Google Scholar] [CrossRef]
  17. Mendez, P.F.; Lu, Y.; Wang, Y. Scaling analysis of a moving point heat source in steady-state on a semi-infinite soild. J. Heat Trans.-ASME 2018, 140, 081301. [Google Scholar] [CrossRef]
  18. Winczek, J. Analytical solution to transient temperature field in a half-infinite body caused by moving volumetric heat source. Int. J. Heat Mass Transf. 2010, 53, 5774–5781. [Google Scholar] [CrossRef]
  19. Honkonen, J. Contour-ordered Green’s functions in stochastic field theory. Theor. Math. Phys. 2013, 175, 827–834. [Google Scholar] [CrossRef]
  20. Sheng, X.; Lu, X.; Zhang, J.R.; Lu, Y.Q. An analytical solution to temperature field distribution in a thick rod subjected to periodic-motion heat sources and application in ball screws. Eng. Optim. 2021, 53, 2144–2163. [Google Scholar] [CrossRef]
  21. Arif, M.; Kuman, P.; Kuman, W.; Mostafa, Z. Heat transfer analysis of radiator using different shaped nanoparticles water-based ternary hybrid nanofluid with applications: A fractional model. Case Stud. Therm. Eng. 2022, 31, 101837. [Google Scholar] [CrossRef]
  22. Hu, H.P. Theory of Heat Conduction; USTC Press: Hefei, China, 2010. [Google Scholar]
  23. Özisik, M.N. Heat Conduction; John Wiley and Sons: New York, NY, USA, 1980. [Google Scholar]
  24. Fujii, F. Theory of the steady laminar natural convection above a horizontal line heat source and a point heat source. Int. J. Heat Mass Transf. 1963, 6, 597–606. [Google Scholar] [CrossRef]
Figure 1. Displacement curves of practical and simplified motion.
Figure 1. Displacement curves of practical and simplified motion.
Applsci 12 06612 g001
Figure 2. Schematic diagram of moving heat source in an infinite medium.
Figure 2. Schematic diagram of moving heat source in an infinite medium.
Applsci 12 06612 g002
Figure 3. Graphical exhibition of numerical results with various dimensionless parameter values: (a) Dimensionless surplus temperature profiles over time parameter for various values of γ with fixed ξ = 0 and β = 2; (b) Fluctuation amplitude profiles over heat conduction parameter for various values of γ with fixed ξ = 0; (c) Dimensionless surplus temperature profiles over time parameter for various values of β with fixed ξ = 0 and γ = 2; (d) Dimensionless surplus temperature profiles over location parameter for various values of κ with fixed β = 2 and γ = 1; (e) Fluctuation amplitude profiles over location parameter for various values of β with fixed γ = 1.
Figure 3. Graphical exhibition of numerical results with various dimensionless parameter values: (a) Dimensionless surplus temperature profiles over time parameter for various values of γ with fixed ξ = 0 and β = 2; (b) Fluctuation amplitude profiles over heat conduction parameter for various values of γ with fixed ξ = 0; (c) Dimensionless surplus temperature profiles over time parameter for various values of β with fixed ξ = 0 and γ = 2; (d) Dimensionless surplus temperature profiles over location parameter for various values of κ with fixed β = 2 and γ = 1; (e) Fluctuation amplitude profiles over location parameter for various values of β with fixed γ = 1.
Applsci 12 06612 g003
Figure 4. Semi-infinite medium subjected to reciprocating a point heat source.
Figure 4. Semi-infinite medium subjected to reciprocating a point heat source.
Applsci 12 06612 g004
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sheng, X.; Xu, Y.; Huang, D.; Zhang, J.; Lu, Y.; Lu, X. Analytical Solutions to Temperature Field in Various Relative-Scale Media Subjected to a Reciprocating Motion Point Heat Source. Appl. Sci. 2022, 12, 6612. https://doi.org/10.3390/app12136612

AMA Style

Sheng X, Xu Y, Huang D, Zhang J, Lu Y, Lu X. Analytical Solutions to Temperature Field in Various Relative-Scale Media Subjected to a Reciprocating Motion Point Heat Source. Applied Sciences. 2022; 12(13):6612. https://doi.org/10.3390/app12136612

Chicago/Turabian Style

Sheng, Xin, Yadong Xu, Dacheng Huang, Jianrun Zhang, Yunqiao Lu, and Xi Lu. 2022. "Analytical Solutions to Temperature Field in Various Relative-Scale Media Subjected to a Reciprocating Motion Point Heat Source" Applied Sciences 12, no. 13: 6612. https://doi.org/10.3390/app12136612

APA Style

Sheng, X., Xu, Y., Huang, D., Zhang, J., Lu, Y., & Lu, X. (2022). Analytical Solutions to Temperature Field in Various Relative-Scale Media Subjected to a Reciprocating Motion Point Heat Source. Applied Sciences, 12(13), 6612. https://doi.org/10.3390/app12136612

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop