Services & Resources / Wolfram Forums
-----
 /
MathGroup Archive
2003
*January
*February
*March
*April
*May
*June
*July
*August
*September
*October
*November
*December
*Archive Index
*Ask about this page
*Print this page
*Give us feedback
*Sign up for the Wolfram Insider

MathGroup Archive 2003

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

Search the Archive

Trouble with iterative function

  • To: mathgroup at smc.vnet.net
  • Subject: [mg39937] Trouble with iterative function
  • From: "Miguel Cruz" <pyro_gest at hotmail.com>
  • Date: Wed, 12 Mar 2003 02:31:25 -0500 (EST)
  • Sender: owner-wri-mathgroup at wolfram.com

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?