Register to reply 
Y' = ayb*y^2c*y*x; x'=d*y^2 
Share this thread: 
#1
Nov510, 05:26 AM

P: 2

Hi,
I am a chemist with unfortunately too little experience with differential equations, and now I have a problem I need to solve: I have a system of DE of the form: y' = ayb*y^2c*y*x x'=d*y^2 a, b, c, and d are constants (related to each other) I have been trying to read about this, but I got stuck after having read about first order systems. I would very much appreciate any help, and I hope that I am posting in the right forum. 


#2
Nov610, 05:39 AM

P: 11

Hi cesca,
We'll need a little more information: initial conditions, the relationships between a, b, c and d, the interval where you hope to find solutions. Out of curiosity, I'd be interested in knowing the system you are attempting to model. I'm assuming all derivatives are with respect to t? 


#3
Nov710, 04:10 AM

P: 759

y' = ayb*y^2c*y*x
x'=d*y^2 supposing that y = y(t) You can find a particular solution on the form : y = k/t x = h d*kČ/t k and h constants to be computed (related to a, b, c, d) 


#4
Nov710, 06:59 AM

P: 1,666

Y' = ayb*y^2c*y*x; x'=d*y^2



#5
Nov710, 01:25 PM

P: 759

Sure, Jackmell is right : numerical calculus is tne best method on a practical viewpoint.
Nevertheless, analytical solution is possible (attachment). The result is obtained on a parametric form but the formulas are rather complicated. 


#6
Nov710, 02:00 PM

P: 1,666

You have: Let: [tex]X=CxA;\quad dX=Cdx\rightarrow DC\frac{Dy}{dX}=\frac{X}{y}B[/tex] should that not be: [tex]X=CxA;\quad dX=Cdx\rightarrow DC\frac{Dy}{dX}=\frac{X}{y}B[/tex] 


#7
Nov710, 03:49 PM

P: 759

You are right, that's a careless mistake. May be not the only one !
The point is the method of resolution. Well, to be corrected and more important, to be verified. 


#8
Nov810, 02:25 AM

P: 2

Wow  a lot of answers. I will try to go trhough them and ask again if I get stuck (which I probably will). A numerical solution is perfectly acceptable, expecially if the analytical one is very complicated (or not possible to get).
Odysseus is right, it is a timederivative, and I am looking at how the concentration of some chemical species change over time. I have a suggestino for a reaction mechanism, and sould like to know if that suggestion fits with my experimental data. Relationships between a, b, c, and d are: a = y*k*OHs b = y*(k+l) c = y*k d= y*l Initial conditions: I am not entirely sure, [hc], should be small but not zero, since this reaction mechanism would not be the predominant one before some [hc] has been buildt up. Thank a lot! 


#9
Nov810, 05:39 PM

P: 1,666

Hi. I've compared Jacquelin's method against the numerical results although I avoided a value that would give X=0 since we're taking logs below. So I let:
[tex]a=1[/tex] [tex]b=1[/tex] [tex]c=2[/tex] [tex]d=1[/tex] [tex]x_0=1[/tex] [tex]y_0=2[/tex] I believe there is a mistake in the last line of Jac's analysis. Starting with the expression: [tex]\frac{dX}{X}=\frac{DCYdY}{1+BY+CDY^2}[/tex] then: [tex]\ln(x)\biggr_{X_0}^{X}=J(Y)J(Y_0)[/tex] where J(Y) is the antiderivative on the right side and [itex]X_0=Cx_0a[/tex] and [itex]Y_0=\frac{y_0}{X_0}[/tex] X(Y) and Y are as before but I think t is calculated this way unless I'm missing something: [tex]t=\frac{1}{D}\int \frac{dx}{y^2}=\frac{1}{CD}\int \frac{d(X(Y))}{X(Y)^2 Y^2}[/tex] This gives a positive value of t for Y in the range of (0,2]. Here's the code I used to compare the numerical results to the analytic results. The final plot superimposes the two methods and they agree very well in the range I computed them although there are some further considerations for Y very close to zero.



Register to reply 