Solving Equations: From t=0 to 1d

Join the discussion
Registration is free. Start your own thread to ask a follow-up.
4 replies · 3K views · last reply
jegues
Messages
1,085
Reaction score
3

Homework Statement



See figure attached for problem statement.

Printed textbook page, problem 1.6, on a radioactive contaminant in a closed reactor: decay rate = -kc, mass balance dc/dt = -kc, and parts (a) Euler's method from t = 0 to 1 d with k = 0.2 d^-1 and step 0.1, and (b) a semilog plot.


Homework Equations





The Attempt at a Solution



Here's as far as I got,

Handwritten numerical problem on grid paper: decay rate -Ky with K in day^-1, ODE y' = -Ky, Euler steps y1 = y0 + 0.1 F(x0, y0) with y(0) = 10 Bq/L, giving y1 = 10 - K.


The part that confuses me is the range we should solve this equation. It says,

[tex]\text{from } t=0 \text{ to } 1d[/tex]

I put the [tex]^{-1}[/tex] in there with pencil because I thought it was a typo. Is it?

If I can figure out where to stop I really just have to keep repeating the formula,

[tex]y_{n} = y_{n-1} + \Delta t F(x_{n-1}, y_{n-1})[/tex]

Right?

Thanks again!
 
Attachments
  • Printed textbook page, problem 1.6: radioactive contaminant decay in a closed reactor, with the ODE dc/dt = -kc and parts (a) Euler's method from t = 0 to 1 d and (b) semilog plot of ln c versus t.
    Q1.6.JPG
    61.1 KB · Views: 615
  • Handwritten numerical methods solution on grid paper: problem 1.6 with decay rate -Ky, y' = -Ky, part a from t=0 to 1d with k=0.2 d^-1 and step 0.1, y(0)=10, and y1 = 10 + (0.1)(-10K) = 10 - K.
    A1.6.jpg
    34.1 KB · Views: 616
Physics news on Phys.org
jegues said:
[tex]\text{from } t=0 \text{ to } 1d[/tex]

I put the [tex]^{-1}[/tex] in there with pencil because I thought it was a typo. Is it?
No, t is time and it's units are days or d, not d-1.

If I can figure out where to stop I really just have to keep repeating the formula,

[tex]y_{n} = y_{n-1} + \Delta t F(x_{n-1}, y_{n-1})[/tex]

Right?

Thanks again!
Yes, that is the correct approach. You stop once you have reached t = 1 day.
 
Redbelly98 said:
No, t is time and it's units are days or d, not d-1.


Yes, that is the correct approach. You stop once you have reached t = 1 day.

But how will I reach 1 day?

I have,

[tex]10-k[/tex] where k has units [tex]day^{-1}[/tex].

Can you explain please?
 
Well, if you're going to be careful about the units -- and it's a good thing if you are -- then you need to include the units correctly on all quantities. Note that Δt should really be 0.1d, not simply 0.1 as stated in the problem statement.

So you really have

y1 = y0 - k·y0·Δt
. . .= 10 Bq/L - (0.2 d-1)*(10 Bq/L)*(0.1 d)
. . .= 10 Bq/L - 0.2 Bq/L = 9.8 Bq/L

That is y1, so that is the concentration at 1·Δt or 0.1 days.
When you calculate yn, you will have the concentration after n·Δt or n·0.1 days.