11.10.2013 Views

Conference, Proceedings

Conference, Proceedings

Conference, Proceedings

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

- from numeric solution of Eqs (9) and (10) expressed as [12]<br />

256<br />

2 ⎛ π f ⎞ ⎛ π f ⎞<br />

= cosh⎜<br />

⎟ + ⎜ ⎟<br />

⎜<br />

2 ⋅ l<br />

⎟<br />

cos<br />

⎜<br />

2 ⋅ l<br />

ψ<br />

⎟<br />

⎝ a ⎠ ⎝ a ⎠<br />

⎛ π f ⎞ ⎛ π f ⎞<br />

tg ϕ = tg⎜<br />

⎟ ⎜ ⎟<br />

⎜<br />

⋅ l<br />

⎟<br />

tanh<br />

⎜<br />

⋅ l<br />

⎟<br />

(21)<br />

⎝ a ⎠ ⎝ a ⎠<br />

4. Numerical modelling<br />

Because of a complex character of the thermal diffusivity identification problem direct analyses<br />

of any negative effects on the measurement procedure are difficult. For this reason a procedure<br />

of numerical validation of the measurement procedure has been started. The numerical analysis<br />

discussed here is of a preliminary character. It refers to real experiments performed while the<br />

temperature oscillation procedure has been validated [12]. This is why certain thermophysical<br />

parameters have been assumed for the modelled media (Table 1). Despite the fact that the<br />

analysis has been focussed mostly on the effect of a scanning mode procedure it was started<br />

from consideration of the heat loses and finite specimen lateral dimension on the identification<br />

result (Fig. 2). The analysis has been performed applying finite element modelling (FEM)<br />

software Comsol. Prior to main calculations all the numerical modelling procedures have been<br />

tested and optimal parameter values for the time step, finite element type and dimension etc.<br />

have been established. The numerical details have not been discussed here as we will rather<br />

focus on the obtained results.<br />

At first effects of convective heat transfer loses and a thermocouple (TC) location on the thermal<br />

diffusivity identification have been studied. The calculations have been performed applying<br />

thermal properties of the ice‐type solid with the temperature oscillation applied to the bottom<br />

specimen surface (Fig. 3; z=0 mm) while the upper and the side surfaces have been exposed to<br />

convection into the ambient of a certain temperature. Direct effects of the heat loses on the<br />

temperature recordings are shown in Fig. 4. Results of the thermal diffusivity identification<br />

listed in Table 2 prove that these effects are of a minor importance regarding thermal insulation<br />

applied in real experiment [12]. It should be underlined that the expected real value of the heat<br />

transfer coefficient is close to 4 W⋅m‐ 2 ⋅K ‐1 . The observed effects of TC location are non<br />

significant.<br />

Material<br />

Table 1: Thermal properties of the modelled media<br />

Density Heat capacity Thermal conductivity<br />

kg⋅m 3 J⋅kg ‐1 ⋅K W⋅m ‐1 ⋅K ‐1<br />

model ice‐like solid 925 2200 2<br />

model H2O<br />

(Comsol Script notation)<br />

925+75*flc2hs(T‐<br />

273,4)+(273‐T)<br />

2200+83250*(flc2hs(T‐<br />

271,2)‐flc2hs(T‐<br />

275,2))+2000*flc2hs(T‐<br />

273,4)<br />

2.2‐1.7*flc1hs(273,10)<br />

(20)

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

Saved successfully!

Ooh no, something went wrong!