FEATool  v1.7
Finite Element Analysis Toolbox
 All Files Functions Pages
ex_heattransfer2.m File Reference


EX_HEATTRANSFER2 1D Stationary heat transfer with radiation.

[ FEA, OUT ] = EX_HEATTRANSFER2( VARARGIN ) NAFEMS T2 benchmark example for heat transfer with radiation [1]. The left end of a 0.1 m rod is held at a temperature of 1000 K while the right end is radiating with an emissivity, em = 0.98, and Stefan-Bolzmann constant, sigma = 5.67e-8 Wm^2/K^4.

       +---------- L=0.1m ----------+ T(0.1)?
   T=1000K                   q_n = em*sigma*(T_amb^4-T^4)

The rod is made of iron with density 7850 kg/m^3, heat capacity 460 J/kgK, and thermal conductivity 55.563 W/mK. The steady state temperature at the left end is sought when the surrounding ambient temperature is 300 K.

  [1] The Standard NAFEMS Benchmarks,
      The National Agency for Finite Element Standards, UK, 1990.

Accepts the following property/value pairs.

Input       Value/{Default}        Description
hmax        scalar {0.02}          Grid cell size
sfun        string {sflag1}        Finite element shape function
istat       scalar {1}/0           Use stationary (=1), or time dependent solver
iplot       scalar {1}/0           Plot solution (=1)
Output      Value/(Size)           Description
fea         struct                 Problem definition struct
out         struct                 Output stuct
Generated fields of fea:
Generated fields of out: