Heather's Mathematical Model Question on Herbicide Absorption

Click For Summary
SUMMARY

The discussion centers on a mathematical model for herbicide absorption in plant leaves, detailing the differential equations governing the process. The initial equation for concentration, C(t), is derived as $$C(t)=C_a(1-e^{-\alpha t})$$ when no herbicide is present initially. When accounting for the plant's circulatory system, the modified equation becomes $$C(t)=\frac{\alpha C_a}{\alpha+\beta}(1-e^{-(\alpha+\beta)t})$$, illustrating how the herbicide disperses throughout the plant. The key difference between the two models is the limiting concentration, with the first model approaching Ca and the second approaching $$\frac{\alpha C_a}{\alpha+\beta}$$.

PREREQUISITES
  • Understanding of ordinary differential equations (ODEs)
  • Familiarity with initial value problems (IVPs)
  • Knowledge of integrating factors in solving linear ODEs
  • Basic concepts of herbicide absorption dynamics
NEXT STEPS
  • Study the application of separation of variables in ODEs
  • Explore the use of integrating factors in linear differential equations
  • Research the dynamics of chemical absorption in biological systems
  • Learn about the impact of rate constants on concentration over time
USEFUL FOR

Mathematicians, biologists, and environmental scientists interested in modeling chemical absorption processes, as well as students studying differential equations and their applications in real-world scenarios.

MarkFL
Gold Member
MHB
Messages
13,284
Reaction score
12
Here is the question:

Mathematical Model Question?

A simple model for the absorption of a herbicide placed on the surface of a leaf is that it will be absorbed at a rate proportional to the difference in concentration between that on the surface and that on the interior. Assume the rate constant is α.

(a) If a layer with constant concentration Ca is placed on the leaf, write an equation for the concentration of herbicide in the leaf, C(t) as time progresses.

(b) Solve the equation for concentration C(t) if there is no herbicide in the leaf to begin.

(c) Now suppose the plant circulatory system (xylem and phloem) disperses the herbicide throughout the plant with rate proportional to concentration in the leaf and rate constant β. Modify your equation accordingly and hence your solution.

(d) What is the difference in outcome between cases in (b) and (c)?

I have posted a link there to this topic so the OP can see my work.
 
Physics news on Phys.org
Hello Heather,

(a) The time rate of change of herbicide in the leaf is proportional to the difference between the constant concentration on the surface and the variable concentration in the leaf, so we may state:

$$\frac{dC}{dt}=\alpha\left(C_a-C \right)$$ where $$0<\alpha$$

(b) We are given the initial value $$C(0)=0$$ and asked to solve the resulting initial value problem (IVP).

To solve the ordinary differential equation (ODE), we may either separate variables or write in standard linear form, and use in integration factor. I will demonstrate both methods.

i) Separate variables:

We may thus write the ODE as:

$$\frac{1}{C-C_a}\,dC=-\alpha\,dt$$

Integrate, switching dummy variables and using the boundaries as the limits of integration:

$$\int_0^{C}\frac{1}{u-C_a}\,du=-\alpha\int_0^t\,dv$$

Applying the fundamental theorem of calculus (FTOC), we obtain:

$$\left[\ln\left|u-C_a \right| \right]_0^C=-\alpha\left[v \right]_0^t$$

$$\ln\left|\frac{C-C_a}{-C_a} \right|=-\alpha t$$

Convert from logarithmic to exponential form:

$$\frac{C-C_a}{-C_a}=e^{-\alpha t}$$

Solving for $C(t)$, we obtain:

$$C(t)=C_a\left(1-e^{-\alpha t} \right)$$

ii) Express as linear ODE

$$\frac{dC}{dt}+\alpha C=\alpha C_a$$

Compute the integrating factor:

$$\mu(t)=e^{\alpha\int\,dt}=e^{\alpha t}$$

Multiply the ODE by this factor:

$$e^{\alpha t}\frac{dC}{dt}+\alpha e^{\alpha t}C=\alpha C_ae^{\alpha t}$$

Rewrite the left side as the differentiation of a product:

$$\frac{d}{dt}\left(e^{\alpha t}C \right)=\alpha C_ae^{\alpha t}$$

Integrate with respect to $t$

$$\int\,d\left(e^{\alpha t}C \right)=\alpha C_a\int e^{\alpha t}\,dt$$

$$e^{\alpha t}C=C_ae^{\alpha t}+c_1$$

Solve for $C(t)$:

$$C(t)=C_a+c_1e^{-\alpha t}$$

Use the initial values to determine the parameter $c_1$ (the constant of integration):

$$C(0)=C_a+c_1=0\,\therefore\,c_1=-C_a$$

And so we find the solution satisfying the given conditions is:

$$C(t)=C_a\left(1-e^{-\alpha t} \right)$$

(c) Here we are told herbicide is also leaving the leaf at a rate proportional to $C(t)$, and so the ODE becomes:

$$\frac{dC}{dt}=\alpha\left(C_a-C \right)-\beta C$$ where $$C(0)=0$$ and $$0<\alpha,\beta$$

While we could still use both methods to solve the resulting IVP, let's just use the linear method. So, let's express the ODE in standard linear form:

$$\frac{dC}{dt}+(\alpha+\beta)C=\alpha C_a$$

Computing the integrating factor, we find:

$$\mu(t)=e^{(\alpha+\beta)\int\,dt}=e^{(\alpha+\beta)t}$$

And so the ODE becomes:

$$e^{(\alpha+\beta)t}\frac{dC}{dt}+(\alpha+\beta)e^{(\alpha+\beta)t}C=\alpha C_ae^{(\alpha+\beta)t}$$

Expressing the left side as the differentiation of a product, we obtain:

$$\frac{d}{dt}\left(e^{(\alpha+\beta)t}C \right)=\alpha C_ae^{(\alpha+\beta)t}$$

Integrate with respect to $t$:

$$\int\,d\left(e^{(\alpha+\beta)t}C \right)=\alpha C_a\int e^{(\alpha+\beta)t}\,dt$$

$$e^{(\alpha+\beta)t}C=\frac{\alpha C_a}{\alpha+\beta}e^{(\alpha+\beta)t}+c_1$$

Solving for $C(t)$, we get:

$$C(t)=\frac{\alpha C_a}{\alpha+\beta}+c_1e^{-(\alpha+\beta)t}$$

Use the initial values to determine the parameter $c_1$ (the constant of integration):

$$C(0)=\frac{\alpha C_a}{\alpha+\beta}+c_1=0\,\therefore\,c_1=-\frac{\alpha C_a}{\alpha+\beta}$$

And so we find the solution satisfying the given conditions is:

$$C(t)=\frac{\alpha C_a}{\alpha+\beta}\left(1-e^{-(\alpha+\beta)t} \right)$$

(d) The main difference I see is in the limiting concentration in the leaf, that is the concentration as $t\to\infty$:

For the model in part (b), we have:

$$\lim_{t\to\infty}C(t)=C_a$$

And for the model in part (c), we have:

$$\lim_{t\to\infty}C(t)=\frac{\alpha C_a}{\alpha+\beta}$$

Notice if $$\beta=0$$, then the two models are equivalent, as we should expect.

Another difference is that the second model approaches its limiting value more rapidly, since the herbicide is being dispersed, or removed from the leaf, resulting in a slower rate of change. Here is a plot of two possible curves, letting $C_a=\alpha=\beta=1$:

View attachment 1215
 

Attachments

  • heather.jpg
    heather.jpg
    5.5 KB · Views: 114

Similar threads

  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 4 ·
Replies
4
Views
2K
Replies
1
Views
6K
  • · Replies 3 ·
Replies
3
Views
4K
  • · Replies 2 ·
Replies
2
Views
4K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 3 ·
Replies
3
Views
10K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 2 ·
Replies
2
Views
2K
Replies
2
Views
3K