Separable Equations — Question 7

PDF ↗

Question 7

A positive concentration obeys a power-law decay model dydt=−kyp,k>0,p∈ℝ.\frac{dy}{dt}=-k y^p,\qquad k>0,\quad p\in\mathbb R. Use fixed units for time and concentration. Exact observations are y(0)=16,y(1)=8,y(2)=4.y(0)=16,\qquad y(1)=8,\qquad y(2)=4. Assume one pair of constants (k,p)(k,p) describes the entire observation interval.

Tasks

  1. Express each of the two observed one-unit decay times as a separated definite integral, without assuming p=1p=1.

  2. Prove that these observations determine pp uniquely. Treat the case p=1p=1 without dividing by 1−p1-p.

  3. Find kk and the resulting solution, and predict when the concentration first reaches 11.

  4. Explain why the first two observations alone would not determine pp.

Original worksheet page 1: question and worked solution for 2-2-007
Show solutionHide solution

Question 7 – Solution

Strategy. Compare passage-time integrals on two concentration intervals related by scaling; the unknown constant kk then cancels.

Step 1: Write the two passage times. Since y>0y>0 and k>0k>0, the concentration decreases strictly. Separation gives 1=1k∫816s−pds,1=1k∫48s−pds.1=\frac 1k\int_8^{16}s^{-p}\,ds,\qquad 1=\frac 1k\int_4^8s^{-p}\,ds. These integrals are positive and finite for every real pp.

Step 2: Recover the exponent. In the first integral put s=2us=2u: ∫816s−pds=21−p∫48u−pdu.\int_8^{16}s^{-p}\,ds=2^{1-p}\int_4^8u^{-p}\,du. The observed equal passage times force 21−p=12^{1-p}=1, hence p=1\boxed{p=1}. Conversely, p=1p=1 makes the two integrals equal. No division by 1−p1-p was used, so the logarithmic case was retained throughout.

Step 3: Determine the rate and prediction. For p=1p=1, k=∫816dss=ln⁡2,ln⁡y16=−tln⁡2.k=\int_8^{16}\frac{ds}{s}=\ln 2,\qquad \ln\frac{y}{16}=-t\ln 2. Therefore y(t)=162−t\boxed{y(t)=16\,2^{-t}}. It gives exactly 16,8,416,8,4 at the measured times, and differentiating gives y′=−(ln⁡2)yy'=-(\ln 2)y. It is positive for every real tt. The value 11 is first reached at t=4\boxed{t=4} because the solution is strictly decreasing.

Step 4: Identify the information supplied by the third observation. With only y(0)=16y(0)=16 and y(1)=8y(1)=8, each real pp could be paired with k(p)=∫816s−pds>0.k(p)=\int_8^{16}s^{-p}\,ds>0. The corresponding separated solution decreases from 1616 to 88 in one unit of time and stays positive during that interval. Thus infinitely many exponents fit the first two data; the second measured halving time identifies the exponent within this model family.

Original worksheet page 2: question and worked solution for 2-2-007

Original worksheet layout. Use Enlarge or open the PDF for a closer view.