Finite element approaches to the PTC thermistor problem
Introduction
In this paper we consider an initial boundary value problem of coupled non-linear partial differential equations which consists of a heat-flow equation with the joule heating as a source and a current flow equation with a temperature dependent electrical conductivity. This problem is a thermo-electric problem called as a PTC thermistor problem which is mathematically modeled as follows:
Heat-flow problem:
Current-flow problem:where U(x,t) is the temperature, φ(x,t) is the electrical potential, β is a surface heat transfer coefficient, α is the ratio of electrical conductivity and σ is a temperature dependent electrical conductivity which decreases as the temperature increases.
The existence, uniqueness and behaviour of solutions to this problem were given in [1], [2], [3], [4], [5], [6]. Numerical solutions of the PTC thermistor problem have been undertaken by employing various forms of the finite element [7], [8], [9], finite difference [10], [11], [12], [13] and heat balance integral methods [14], [15].
In this paper, SC and PG methods have been applied to the PTC thermistor problem with an electrical conductivity proposed in the next section to obtain its approximate steady-state solutions. The numerical results obtained by both methods are compared with exact ones and also earlier work.
Section snippets
Electrical conductivities for the problem
Many authors [7], [11], [15] considered the above given thermistor problem with a step function electrical conductivity of the form
In a postscript, Howison [16] makes the point that a step function is a crude representation for σ(U). Recently, Kutluay and Wood [12] introduced a slightly more realistic model for σ(U) of the formwhich is mathematically equivalent to a ramp function. They have applied a standard explicit finite
Warm/cold solution
We obtain the exact steady-state solution to this phase by using a standard analytical method aswhereInserting into Eq. (13) leads to a non-linear equation for the free interface s1:which can be solved numerically. If s1 becomes equal to 1, Eq. (14) givesand so the problem enters a warm phase. Since 0<s1⩽1 in this phase, we obtain an
The finite element solutions
The region [0,1] is partitioned into N finite elements of equal length h by the knots xi such that 0=x0<x1<⋯<xN−1<xN=1. Over the element [xm,xm+1], an approximation to the exact solution U(x,t) in terms of quadratic B-splines Qm with the required properties can be given as [17]where the δm are time dependent element parameters to be determined from the boundary and weighted residual conditions. In terms of a local coordinate system ξ given by ξ=x−xm (0⩽ξ⩽h) the base
Numerical results and conclusion
In order to obtain the steady-state temperature distribution U(x,t) and the locations of the free boundaries s1 and s2, respectively, in the warm/cold and hot/warm phases constructed in the previous sections we used criteriawith k=0.01 for both the SC and PG solutions.
In order to show how good the numerical solutions of the problem in comparison with the exact ones we used the error norm L2 defined by
To compare our results
References (19)
- et al.
Stationary solutions to the thermistor problem
J. Math. Anal. Appl.
(1993) - et al.
Numerical solutions of the thermistor equations
J. Comput. Appl. Math.
(1997) - et al.
A variety of finite difference methods to the thermistor with a new modified electrical conductivity
Appl. Math. Comput.
(1999) - et al.
Numerical solutions of the thermistor problem with a ramp electrical conductivity
Appl. Math. Comput.
(2004) - et al.
A heat balance integral model of the thermistor
Int. J. Heat Mass Transfer
(1995) Remark on existence and uniqueness for the thermistor problem under mixed boundary conditions
Q. Appl. Math.
(1989)- et al.
Existence results for a non-linear elliptic system modelling a temperature dependent electrical resistor
Ann. Mat. Pura Appl.
(1992) - et al.
Temperature surges in thermistors
SIAM J. Appl. Math.
(1992) The thermistor problem
Nonlinear Differ. Equations Appl.
(1997)