MathGroup Archive 2003

[Date Index] [Thread Index] [Author Index]

Search the Archive

Trouble with iterative function


Good day everyone,

I have problem trying to implement an iterative function that attempts 
to estimate the temperature of a wood particle subjected to variable air 
temperature (air temperature is function of time, e.g. 
Tair[t_]=600+t*5).
The function is the following (hc is the heat transfer coeficient):
Fignition[t_, hc_] :=
  Block[{ Tfuel = Tini, Tini = 300, Fcp =
      2100, Fd = 0.005, Fl = 0.1, ?f = 520,
      FMC = 0.7, Wcp = 2417, Whv = 2257000},
    {While[Tfuel < 373, Print[Tfuel]; Save["test01.dat", Tfuel];{S = 
?f*Fv[
      Fd, Fl]*(Fcp + (Wcp*FMC));
             Q = hc*Fa[Fd, Fl]*(Tair[t] - Tfuel);
             Xtp = Tfuel + (Q/S);
             Tfuel = Xtp}]}
      {While[Tfuel > 373, Print[Tfuel]; Save["test01.dat", Tfuel]; {
          S = ?f*Fv[Fd, Fl]*(Fcp);
             Q = hc*Fa[Fd, Fl]*(Tair[t] - Tfuel);
             Xtp = Tfuel + (Q/S);
             Tfuel = Xtp}; If[Tfuel > 600, Break[]]
          ]
        }
    ]

My problem is that I want the Fignition[t_, hc_]  function to be 
computed for each time step by reading the Tair at that time step. the 
time step should be the iterator. Currently the function does its thing 
but uses a fixed Tair.
Hope I was clear. Thanks
Miguel


********************************************
Miguel Cruz
460 "The Bullpen" Science Complex
School of Forestry
University of Montana
Missoula, MT 59812

Phone: (406) 329 4748 (@ FiSL)
Fax:   (406) 329 4825 (@ FiSL)
********************************************




  • Prev by Date: Re: general solution for element of series
  • Next by Date: What's legit here?
  • Previous by thread: Re: Calling Install[] from packages
  • Next by thread: What's legit here?