MHB Peak in single ODE within a system

  • Thread starter Thread starter Carla1985
  • Start date Start date
  • Tags Tags
    Ode Peak System
Click For Summary
The discussion centers on a system of nonlinear ordinary differential equations (ODEs) and the conditions under which the variable Z exhibits a peak in its time course. It is noted that a peak occurs when the derivative of Z changes from positive to negative, leading to the condition where αk₊XY must exceed βZ. The conservation law D = X + Y + 2Z is confirmed to be a constant, allowing for simplifications in the system. The constraints on the parameters α, β, k₊, and k₋ are established as non-negative values. The conversation concludes with an acknowledgment of the need to explore these conditions further to understand the behavior of Z.
Carla1985
Messages
91
Reaction score
0
Hi all, I have the system of nonlinear ODEs:

$$
\begin{align}
\frac{dX}{dt}=&-k_+ A X+k_-Y \\
\frac{dY}{dt}=&\ k_+ A X-k_-Y-\alpha k_+ X Y +\beta Z \\
\frac{dZ}{dt}=&\ \alpha k_+ X Y -\beta Z
\end{align}
$$

I also have a conservation law that says $D=X+Y+2Z$. Obviously it is not possible to find exact solutions to the system as it is nonlinear, but I noticed when running an ODE solver that the time course curve of $Z$ sometimes has a peak and then fall, dependant on the parameters. I was wondering if there is any way to be able to find a condition on whether I get a peak or not?

Regards
Carla
 
Physics news on Phys.org
At a base level, you could say that a peak like that means there is an interval $(a,b)$ containing a point $c$ such that $\dfrac{dZ}{dt}>0 \; \forall \, t\in(a,c)$ and $\dfrac{dZ}{dt}<0 \; \forall \, t\in(c,b)$. So, based on your system of ODE's, you could compare $\alpha k_+ X Y -\beta Z$ with $0$. Put another way, when is $\alpha k_+ X Y > \beta Z?$

Do you have any constraints (other than your conservation law, which looks like it could reduce the order of the system by 1. I assume $D$ is a constant?) on the constants $\alpha, \beta, k_+,$ and $k_-?$
 
Does your "conservation law" mean that D= X+ Y+ 2Z is a constant? If so then, taking the derivative of both sides with respect to t, $\frac{dX}{dt}+ \frac{dX}{dt}+ 2\frac{dZ}{dt}= 0$. Using the derivative formulas given, it follows that $\alpha k_+XY- \beta Z= 0$ so you can replace the XY term by $\beta Z$ to get a linear system.
 
Hi,

Yes, sorry, I should have explained better, D is a constant value. The only other constraints are $\alpha, \beta, k_+, k_-\geq 0$. I will work through your suggestions. Thank you very much for the help.

Regards
Carla
 
Carla1985 said:
D is a constant value. The only other constraints are $\alpha, \beta, k_+, k_-\geq 0$.
Is $A$ a constant as well?
 
Krylov said:
Is $A$ a constant as well?

Yes $A$ is constant. The only variables are $X,Y$ and $Z$
 

Similar threads

  • · Replies 2 ·
Replies
2
Views
833
  • · Replies 6 ·
Replies
6
Views
2K
  • · Replies 13 ·
Replies
13
Views
2K
  • · Replies 2 ·
Replies
2
Views
3K
  • · Replies 6 ·
Replies
6
Views
2K
  • · Replies 2 ·
Replies
2
Views
3K
  • · Replies 16 ·
Replies
16
Views
3K
  • · Replies 1 ·
Replies
1
Views
2K
Replies
4
Views
2K
  • · Replies 2 ·
Replies
2
Views
2K