* TDiff 6.0 Script (Field Precision)
* File: ThermalFlux.TIN
* Date: 10/08/2008
* Time: 09:13:50

* ---- CONTROL ----

Mesh = ThermalFlux
Mode = TVar
Geometry = Cylin
DUnit =   1.0000E+02
TMax =   1.0100E+01
DtMin =   5.0000E-05
DtMax =   1.0000E-02

* ---- REGIONS ----

* Region 1: FLUXSURFACE
Cond(1) =   2.1000E+02
Cp(1) =   9.0000E+02
Dens(1) =   2.7000E+03
Source(1) =   2.4300E+07

* Region 2: ALUMINUM
Cond(2) =   2.1000E+02
Cp(2) =   9.0000E+02
Dens(2) =   2.7000E+03
Source(2) =   2.4300E+07

* ---- DIAGNOSTICS ----
SetTime =   1.0000E+01
History =   0.0000E+00  0.0000E+00
History =   0.0000E+00  2.5000E+00

EndFile


----- ANALYSIS ----

Calculations: aluminum sphere of radius 5.0 cm
  k = 210 W/m2-degK
  rho = 2700 kg/m3
  Cp = 900 J/kg-degK
  V (volume) = 5.236E-4 m3
  S (surface area) = 3.1416E-2 m2
  M (mass) = 1.4137 kg
  Q (total energy added)
  M*Cp*(dT/dt) = Q
  
For 100 degC temperature rise in 10 s, need Q = 1.2723E4 W

1) Volumetric source: Q/V = 2.430E7 W/m3
2) External thermal flux: Q/S = 4.0498 W/m2
3) Equivalent volumetric source in a layer of thickness 0.001 m: 4.0489E8 W/m3

-----------------------------------------------

Setup for a uniform volumetric source

* Region 1: FLUXSURFACE
Cond(1) =   2.1000E+02
Cp(1) =   9.0000E+02
Dens(1) =   2.7000E+03
Source(1) =   2.4300E+07

* Region 2: ALUMINUM
Cond(2) =   2.1000E+02
Cp(2) =   9.0000E+02
Dens(2) =   2.7000E+03
Source(2) =   2.4300E+07

-----------------------------------------------

Setup for a specified thermal flux

* Region 1: FLUXSURFACE
Cond(1) =   2.1000E+02
Cp(1) =   9.0000E+02
Dens(1) =   2.7000E+03
Source(1) =   4.0498E8

* Region 2: ALUMINUM
Cond(2) =   2.1000E+02
Cp(2) =   9.0000E+02
Dens(2) =   2.7000E+03