Solving Diffusion Equations: F'(t) / [1 - F(t)] = p + q F(t)

  • Context: Graduate 
  • Thread starter Thread starter Enuma_Elish
  • Start date Start date
  • Tags Tags
    Diffusion
Click For Summary
SUMMARY

The differential equation F'(t) / [1 - F(t)] = p + q F(t) has the solution F(t) = [1 - exp(-bt)]/[1 + a exp(-bt)], where a = q / p and b = p + q. The discussion highlights the process of deriving a system of equations with solutions F1(t) and F2(t) based on boundary conditions. The solution is confirmed to hold under the condition F(0)=0, leading to the conclusion that K must equal 1 when the ratio of p and q is not -1. The necessity of setting up a system of differential equations for two particular solutions is questioned.

PREREQUISITES
  • Understanding of differential equations, specifically first-order equations.
  • Familiarity with boundary conditions in mathematical modeling.
  • Knowledge of exponential functions and their properties.
  • Ability to manipulate algebraic expressions involving ratios and limits.
NEXT STEPS
  • Study the derivation of solutions for first-order differential equations.
  • Learn about boundary value problems and their implications in differential equations.
  • Explore the application of exponential decay in modeling real-world phenomena.
  • Investigate the implications of varying parameters p and q in differential equations.
USEFUL FOR

Mathematicians, physicists, and engineers involved in modeling dynamic systems, particularly those working with differential equations and boundary conditions.

Enuma_Elish
Messages
78
Reaction score
0
The diff. eq. F'(t) / [1 - F(t)] = p + q F(t) or F'(t) = [p + q F(t)][1 - F(t)] has the solution F(t) = [1 - exp(-bt)]/[1 + a exp(-bt)] where a = q / p and b = p + q.

How can I write the system of equations with the solution {F1(t), F2(t)} = { [1 - exp(-b*t)]/[1 + a exp(-b*t)] , [1 - exp(-b*(t-T))]/[1 + a2 exp(-b*(t-T))] } in a similar way, with F1' and F2' on the left-hand side, and p and q (or equivalently a and b) on the right-hand side?

I tried defining L1(t) = F1(t)[1 - F2(t - T)], L2(t) = F2(t - T)[1 + F1(t)] and differentiating, but that did not seem to get me to L1' / (1 - L1) = p + q*L1, L2' / (1 - L2) = p + q*L2.

Thanks for any thoughts and ideas.
 
Last edited:
Physics news on Phys.org
Hello Enuma_Elish,

I looked into your equation and came up with the following conclusions. The solution to the DE is given by you without any boundary condition. Now solving the DE gives as solution the following with K as integration constant:

[tex]F(t)=\frac{K-e^{-(p+q)t}}{K+\frac{p}{q}e^{-(p+q)t}}[/tex]

Assuming now the following boundary condition, [itex]F(0)=0[/itex] holds, it is easily written as your solution. Indeed, we have:

[tex]F(0)=0=\frac{K-1}{K+\frac{p}{q}}[/tex]

From which, [itex]K=1[/itex] and [itex]\frac{p}{q} \neq -1[/itex]. This gives thus:

[tex]F(t)=\frac{1-e^{-(p+q)t}}{1+\frac{p}{q}e^{-(p+q)t}}[/tex]

However notice the assumption on the ratio of p and q, if this is not fullfilled you will have to integrate the DE again and you will obtain a different solution. Try it.
The second part is to write it as a system of equations, but I don't think it is necessary, you can obtain the other solution by setting the boundary condition to [itex]F(T)=0[/itex]. This gives then [itex]K=e^{-(p+q)T}[/itex], from which:

[tex]F(t)=\frac{1-e^{-(p+q)(t-T)}}{1+\frac{p}{q}e^{-(p+q)(t-T)}}[/tex]

I don't see any reason for setting a system of DEs up for two different particular solutions of the same DE. Hope this helps.

coomast
 

Similar threads

  • · Replies 2 ·
Replies
2
Views
1K
  • · Replies 4 ·
Replies
4
Views
3K
  • · Replies 0 ·
Replies
0
Views
4K
  • · Replies 2 ·
Replies
2
Views
3K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 1 ·
Replies
1
Views
4K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 8 ·
Replies
8
Views
4K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 23 ·
Replies
23
Views
7K