Modeling Plankton Herbivory Dynamics

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

Discussion Overview

The discussion revolves around modeling plankton herbivory dynamics using a set of differential equations. Participants are exploring the process of nondimensionalizing the model, identifying dimensionless parameters, and analyzing the stability of steady states. The focus includes theoretical aspects of the model and mathematical reasoning related to dimensional analysis.

Discussion Character

  • Technical explanation
  • Mathematical reasoning
  • Debate/contested

Main Points Raised

  • One participant presents a dimensionalized model of plankton herbivory and seeks guidance on nondimensionalization.
  • Several participants suggest using substitutions such as $p = \frac{P}{C}$ and $k = \frac{K}{C}$ to create dimensionless variables.
  • There is a discussion about the dimensions of parameters, with some participants asserting that $P$ and $C$ share the same dimensions, while others question the clarity of these statements.
  • One participant expresses uncertainty about how to derive the parameters $r$, $H$, $B$, $A$, and $D$ in the context of the differential equations.
  • Another participant attempts to analyze the stability of the positive steady state by examining the signs of the derivatives of the nondimensionalized equations.
  • There is a mention of a necessary condition for periodic solutions in the parameter plane defined by $a$ and $k$, with a request for clarification on how to establish this condition.

Areas of Agreement / Disagreement

Participants express differing views on the dimensionality of parameters and the appropriateness of certain substitutions. There is no consensus on the best approach to nondimensionalization or the implications for stability analysis.

Contextual Notes

Some participants note that the dimensions of certain parameters are not explicitly defined, leading to uncertainty in their analysis. The discussion includes unresolved mathematical steps related to the derivation of dimensionless parameters.

Who May Find This Useful

This discussion may be useful for researchers and students interested in mathematical modeling of ecological systems, particularly those focusing on herbivory dynamics and dimensional analysis in differential equations.

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
1K
  • · 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