Modeling Plankton Herbivory Dynamics

  • Context: MHB 
  • Thread starter Thread starter Dustinsfl
  • Start date Start date
  • Tags Tags
    Dynamics Modeling
Click For Summary
SUMMARY

The discussion centers on modeling plankton herbivory dynamics using a system of differential equations. The dimensionalized model is represented as $\frac{dP}{dt} = rP\left[(K-P)-\frac{BH}{C+P}\right]$ and $\frac{dH}{dt} = DH\left[\frac{P}{C+P} - AH\right]$, where $r$, $K$, $A$, $B$, $C$, and $H$ are positive constants. Participants discuss the process of nondimensionalization, specifically using substitutions such as $p = \frac{P}{C}$ and $k = \frac{K}{C}$. The conversation also touches on stability analysis in the context of the parameters $a$ and $k$ and their implications for periodic solutions.

PREREQUISITES
  • Understanding of differential equations and their applications in population dynamics.
  • Familiarity with dimensional analysis and nondimensionalization techniques.
  • Knowledge of stability analysis in dynamical systems.
  • Basic concepts of ecological modeling, particularly in herbivory dynamics.
NEXT STEPS
  • Explore the process of nondimensionalization in differential equations.
  • Study stability analysis methods for nonlinear systems.
  • Investigate the implications of parameter planes in dynamical systems.
  • Research periodic solutions in ecological models and their significance.
USEFUL FOR

Mathematicians, ecologists, and researchers involved in modeling population dynamics, particularly those focusing on herbivory interactions in aquatic ecosystems.

Dustinsfl
Messages
2,217
Reaction score
5
This a plankton herbivore model.

The dimensionalized model is

$\displaystyle
\frac{dP}{dt} = rP\left[(K-P)-\frac{BH}{C+P}\right], \quad \frac{dH}{dt} = DH\left[\frac{P}{C+P} - AH\right]
$

where $r$, $K$, $A$, $B$, $C$, and $H$ are positive constants.

The dimensions of K, P, B, H, C have to be population (that is the only way I can see it to make since) then we have pop^2 - pop^2.

Then D or A has to be (pop)^{-1}.

I am trying to nondimensionalize to

$\displaystyle
\frac{dp}{d\tau} = p\left[(k-p) - \frac{h}{1+p}\right], \quad \frac{dh}{d\tau} = dh\left[\frac{p}{1+p} -ah\right]$

I am not sure what is a good starting point. I need a hint on one dimensionless unit.
 
Physics news on Phys.org
Start with $p = \dfrac{P}{C}.$
 
Danny said:
Start with $p = \dfrac{P}{C}.$

How were you able to identify that as a substitution?

$k = \frac{K}{C}$ then correct?
 
Last edited:
Well, $P$ and $C$ have the same dimenson, right? So let $P = Cp$ so that $p$ is dimensionless.
 
Last edited:
Danny said:
Well, $p$ and $C$ have the same dimenson, right? So let $P = pC$ so that $p$ is dimensionless.

How I see the problem I believe they do unless I am wrong. I wasn't told any dimensions but I know P and H are population and t is time.
 
dwsmith said:
How I see the problem I believe they do unless I am wrong.

You might want to rephrase this statement - it makes no sense!
 
Danny said:
You might want to rephrase this statement - it makes no sense!

How I see the problem, I believe P and C have the same dimensions unless I am wrong.
 
You said they were both "population". What different dimensions could there be? I suppose P could be measured in "people" and C in "thousands of people" but the ratio would still be dimensionless- if P is 120000 people and C is 400 "thousands of people", P/C= 12/3 (people/1000 people)= 4/1000.
 
Are these substitutions correct?

$k=\dfrac{K}{C}$ and $\tau=tr$.

I can't figure out $r,H,B$ in the first DE as well and $A$ and $D$ in the second.
 
  • #10
So skipping the nondimensionalizing, next I am trying to show that for 0 < k < 1 the positive steady state is stable by noting the signs $\frac{dp}{d\tau}$ and $\frac{dh}{d\tau}$.

So the derivatives are

$$
k - 2p - \dfrac{h - p}{(1 + p)^2} \ \text{and} \ \dfrac{dp}{1 + p} - 2dah
$$

I don't get how the derivatives are going to help answer this question. Could it make $k - 2p - \dfrac{h - p}{(1 + p)^2} < 0$?
 
  • #11
The dimensions I know are:

$$
p =\frac{P}{C}, \ k=\frac{K}{C}, \ h=\frac{HB}{C^2}, \ \tau=Crt
$$

That means
$$
a = \frac{AC^2}{B} \ \text{and} \ d = \frac{D}{Cr}
$$

But those wouldn't be dimensionless.

I found the problem. A has to be 1 over population. I am still not sure of D though.
 
Last edited:
  • #12
dwsmith said:
This a plankton herbivore model.

The dimensionalized model is

$\displaystyle
\frac{dP}{dt} = rP\left[(K-P)-\frac{BH}{C+P}\right], \quad \frac{dH}{dt} = DH\left[\frac{P}{C+P} - AH\right]
$

where $r$, $K$, $A$, $B$, $C$, and $H$ are positive constants.

The dimensions of K, P, B, H, C have to be population (that is the only way I can see it to make since) then we have pop^2 - pop^2.

Then D or A has to be (pop)^{-1}.

I am trying to nondimensionalize to

$\displaystyle
\frac{dp}{d\tau} = p\left[(k-p) - \frac{h}{1+p}\right], \quad \frac{dh}{d\tau} = dh\left[\frac{p}{1+p} -ah\right]$

I am not sure what is a good starting point. I need a hint on one dimensionless unit.

How can I find the a,k parameter plane so I can answer:
Hence show that in the a,k parameter plane a necessary condition for a periodic solution to exist is that a,k lie in the domain bounded by a = 0 and $a=4(k-1)/(k+1)^3$. Hence show that if a < 4/27 there is a window of values of k where periodic solutions are possible.
 

Similar threads

  • · Replies 2 ·
Replies
2
Views
965
  • · Replies 11 ·
Replies
11
Views
2K
  • · Replies 4 ·
Replies
4
Views
11K
  • · Replies 2 ·
Replies
2
Views
10K
  • · Replies 8 ·
Replies
8
Views
4K
  • · Replies 1 ·
Replies
1
Views
10K
  • · Replies 3 ·
Replies
3
Views
1K
Replies
2
Views
2K
  • · Replies 9 ·
Replies
9
Views
3K
  • · Replies 32 ·
2
Replies
32
Views
3K