11.08.2013 Views

Untitled - Cdm.unimo.it

Untitled - Cdm.unimo.it

Untitled - Cdm.unimo.it

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

266 Polynomial Approximation of Differential Equations<br />

(12.1.5) Γ(0) = 0, λ Γ ′ (t) = ∂U<br />

(Γ(t),t), ∀t ∈]0,T].<br />

∂x<br />

In (12.1.5), λ is a pos<strong>it</strong>ive given constant.<br />

Let us provide a physical interpretation. We assume that the segment [−1,1]<br />

represents a one-dimensional bar of some prescribed material. In<strong>it</strong>ially, the interval<br />

[−1,0[ is in the solid state, while the interval [0,1] is in the liquid state. After appropriate<br />

scaling, the function U represents the temperature of the bar, in such a way that U = 0<br />

corresponds to the solidification temperature of the liquid. Below this level the material<br />

changes <strong>it</strong>s physical state. The in<strong>it</strong>ial configuration is expressed by cond<strong>it</strong>ions (12.1.4).<br />

Then, as the time passes, due to the lower temperature of the solid part, the liquid begins<br />

a solidification process and the boundary between the solid and the liquid regions is<br />

determined by the function Γ. For any t ∈]0,T], in the interval ] − 1,Γ(t)[ the<br />

temperature is governed by the heat equation (12.1.1) (see also section 10.2). In the<br />

remaining portion of the bar, the temperature is assumed to be zero (cond<strong>it</strong>ion (12.1.3)).<br />

Relations (12.1.2) are the boundary cond<strong>it</strong>ions of equation (12.1.1). In particular, at the<br />

endpoint x = −1 the heat flux is set to zero. Thus, the bar is thermally insulated, i.e.,<br />

there is no heat transfer across that point. Finally, the differential relation in (12.1.5)<br />

describes the evolution of Γ. This is known as the Stefan cond<strong>it</strong>ion and states that, for<br />

any t ∈]0,T], the moving point Γ(t) has a veloc<strong>it</strong>y proportional to the heat flux at<br />

(Γ(t),t). The parameter λ is the Stefan constant and <strong>it</strong> is proportional to the latent<br />

heat of the material. For high values of λ solidification is slow, and vice versa. For<br />

simplic<strong>it</strong>y, we do not take into account the effects of heat conduction inside the portion<br />

of liquid material.<br />

According to the l<strong>it</strong>erature, (12.1.1)-(12.1.5) is a single-phase Stefan problem and rep-<br />

resents an example of free boundary problem (see for instance friedman(1959), rubin-<br />

stein(1971), fasano and primicerio(1972), elliott and ockendon(1982), hill and<br />

dewynne(1987), chapter 7). A theoretical analysis shows that U and Γ are uniquely<br />

determined and Γ, as expected, is an increasing function.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!