But then again, I could be wrong.

  • Thread starter Thread starter masterchiefo
  • Start date Start date
  • Tags Tags
    Differential
masterchiefo
Messages
211
Reaction score
2

Homework Statement


Chemical reactions being studied in which a body A undergoes transformations
according to the following scheme:
http://prntscr.com/8shuvb

k1, k2, k3 , k4 are the rate constants .
We denote x (t ), y ( t) , z (t) the respective concentrations of the products A, B, C at a given time t
( t expressed in minutes).
The initial conditions x (0) = 1, y (0) = 0 and z ( 0) = 0 .
Is arranged above the vessel where the reaction takes place by a burette which is poured
product A at a constant speed in the tank. Under these experimental conditions,
functions x , y, z defined on the interval [0 ; + infinite [ check the following differential system :

dx/dt (1-2x +y+ z)
dy/dt (x - y)
dz/dt (x - z)

Question 1:
Calculate d/dt( x + y + z) and , using the initial conditions, deduce that :
y(t) + z(t) = 1 + t - x(t)
Question 2:
Demonstrate that x is a solution of the differential equation (E) : dx/dt + 3x = 2 + t, then resolve
Equation ( E) knowing that it validates the initial condition x (0) = 1 .

Homework Equations

The Attempt at a Solution

Question 1:
dx/dt+dy/dt+dz/dt= 1-2x+y+z+x-y+x-z
dx/dt+dy/dt+dz/dt= 1

integral dx/dt+ integral dy/dt+ integral dz/dt= integral ( 1 dt)
x(t)+y(t)+z(t)=t*c
y(t)+z(t)=t*c-x(t)

Question 2:
dx/dt+3x=2+t
dx/dt =2+t-3x

dx/dt(x) =1
1= 2+x-3x
1=2+-2x
... does not work. :(
I am not sure what to do here, some tips would help.
 
Physics news on Phys.org
masterchiefo said:

Homework Statement


Chemical reactions being studied in which a body A undergoes transformations
according to the following scheme:
http://prntscr.com/8shuvb

k1, k2, k3 , k4 are the rate constants .
We denote x (t ), y ( t) , z (t) the respective concentrations of the products A, B, C at a given time t
( t expressed in minutes).
The initial conditions x (0) = 1, y (0) = 0 and z ( 0) = 0 .
Is arranged above the vessel where the reaction takes place by a burette which is poured
product A at a constant speed in the tank. Under these experimental conditions,
functions x , y, z defined on the interval [0 ; + infinite [ check the following differential system :

dx/dt (1-2x +y+ z)
dy/dt (x - y)
dz/dt (x - z)

Question 1:
Calculate d/dt( x + y + z) and , using the initial conditions, deduce that :
y(t) + z(t) = 1 + t - x(t)
Question 2:
Demonstrate that x is a solution of the differential equation (E) : dx/dt + 3x = 2 + t, then resolve
Equation ( E) knowing that it validates the initial condition x (0) = 1 .

Homework Equations

The Attempt at a Solution

Question 1:
dx/dt+dy/dt+dz/dt= 1-2x+y+z+x-y+x-z
dx/dt+dy/dt+dz/dt= 1

integral dx/dt+ integral dy/dt+ integral dz/dt= integral ( 1 dt)
x(t)+y(t)+z(t)=t*c

You need an additive constant when integrating, not a multiplicative one!
 
Rather a weird question so I wonder if you have reproduced it properly.

First substances are called A, B, C then the same substances are called x , y, z.

Then four rate constants k1, 2, 3, 4, are mentioned, afterwards they have all become 1.

The rate of addition is also 1.

OK if the model has been so simplified, you can simplify it a bit further yourself: since y and z are formed and decompose at exactly the same rate and are initially at the same concentration, y = z at all times it is sufficient to formulate e.g.

Constant supply rate 1 ⇒ x ⇔ w

where the forwards and backward rate constants are both 2; at the end of the calculation get y and z back from y = z = w/2.
____________________

(Corrected)
 
Last edited:
ehild said:
You need an additive constant when integrating, not a multiplicative one!
oops, made a mistake, but is that a good answer ?
y(t)+z(t)=t+c-x(t)

and for question 2 how do I do it?
 
masterchiefo said:
oops, made a mistake, but is that a good answer ?
y(t)+z(t)=t+c-x(t)

and for question 2 how do I do it?
It is correct. Find the value of c from the initial conditions.
Substitute y(t)+z(t) back into the first equation, you get the desired differential equation for x(t).
It is a first-order linear equation, solve with one of the standard methods you certainly have learnt.
 
  • Like
Likes masterchiefo
ehild said:
It is correct. Find the value of c from the initial conditions.
Substitute y(t)+z(t) back into the first equation, you get the desired differential equation for x(t).
It is a first-order linear equation, solve with one of the standard methods you certainly have learnt.
this is the first equation ?
dx/dt + 3x = 2 + thow am I supposed to substitute ?
 
masterchiefo said:
this is the first equation ?
dx/dt + 3x = 2 + thow am I supposed to substitute ?
I meant y(t)+z(t) substituted into dx/dt = (1-2x +y+ z).
You should get the equation dx/dt + 3x = 2 + t, but you have to prove it.
The next task is to solve this differential equation and give x(t) using the initial condition x(0)=1.
 
  • Like
Likes masterchiefo
ehild said:
I meant y(t)+z(t) substituted into dx/dt = (1-2x +y+ z).
You should get the equation dx/dt + 3x = 2 + t, but you have to prove it.
The next task is to solve this differential equation and give x(t) using the initial condition x(0)=1.
dx/dt=1-2x+(y+z)
dx/dt=1-2x+(1+t-x(t))
dx/dt=2-2x+t-x(t)
dx/dt=2-3x+t
dx/dt+3x=2+t

by doing this I prove that X is a solution ?
and now I just resolve this as a normal linear differential equation ?
 
You derived the de for x(t). Resolve it as a normal linear first order differential equation.
 
  • Like
Likes masterchiefo
  • #10
masterchiefo said:
dx/dt (1-2x +y+ z)
dy/dt (x - y)
dz/dt (x - z)
A differential system is a system of equations, each of which should have = in it somewhere.
Is this the system?
dx/dt = 1 - 2x + y + z
dy/dt = x - y
dz/dt = x - z
masterchiefo said:
Question 1:
Calculate d/dt( x + y + z) and , using the initial conditions, deduce that :
y(t) + z(t) = 1 + t - x(t)
Question 2:
Demonstrate that x is a solution of the differential equation (E) : dx/dt + 3x = 2 + t
This doesn't make any sense to me. You need to have x as some function of t. The differential equation you show above (equation E) is a first-order, nonhomogeneous equation. One technique for solution is to find an integrating factor.
masterchiefo said:
, then resolve
Equation ( E) knowing that it validates the initial condition x (0) = 1 .
 
  • Like
Likes masterchiefo
  • #11
ehild said:
You derived the de for x(t). Resolve it as a normal linear first order differential equation.
okay, 1 and 2 is done thank you very much. This problem confusing me a lot.

Now for Question 3:
Prove that d/dt(y - z) + y - z = 0 and deduce that y = z
d/dt(y - z) = 0 so,.
y-z=0
y = z

Not sure if I am doing this right.
 
  • #12
masterchiefo said:
okay, 1 and 2 is done thank you very much. This problem confusing me a lot.

Now for Question 3:
Prove that d/dt(y - z) + y - z = 0 and deduce that y = z
d/dt(y - z) = 0 so,.
y-z=0
y = z

Not sure if I am doing this right.
How did you get the equation in red?

Start with the original equations in the problem.
dy/dt (x - y)
dz/dt (x - z)
Subtract them. What do you get?
 
  • #13
ehild said:
How did you get the equation in red?

Start with the original equations in the problem.
dy/dt (x - y)
dz/dt (x - z)
Subtract them. What do you get?
(x-y)-(x-z)= -y+z

d/dt(y - z) + (y - z) = 0
d/dt(-y+z) + (-y+z) =0
 
  • #14
masterchiefo said:
(x-y)-(x-z)= -y+z

d/dt(y - z) + (y - z) = 0
d/dt(-y+z) + (-y+z) =0
Sorry, the original equation missed the equal signs.

So it should have been
dy/dt= (x - y)
dz/dt =(x - z)
When you subtract them, it becomes d/dt(y-z)+(y-z)=0
You can consider y-z as a new function v = y-z, and you have the differential equation dv/dt+v=0. What is the solution for v(t) and what is the initial condition?
 
  • #15
ehild said:
Sorry, the original equation missed the equal signs.

So it should have been
dy/dt= (x - y)
dz/dt =(x - z)
When you subtract them, it becomes d/dt(y-z)+(y-z)=0
You can consider y-z as a new function v = y-z, and you have the differential equation dv/dt+v=0. What is the solution for v(t) and what is the initial condition?
I don't understand how it becomes d/dt(y-z)+(y-z)=0
when I subtract them I get (x - y)-(x - z) = z-y
 
  • #16
masterchiefo said:
I don't understand how it becomes d/dt(y-z)+(y-z)=0
when I subtract them I get (x - y)-(x - z) = z-y
Subtract both sides. On the left sides, you have the derivatives of y and z : dy/dt and dz/dt. What is their difference?
 
  • Like
Likes masterchiefo
  • #17
ehild said:
Subtract both sides. On the left sides, you have the derivatives of y and z : dy/dt and dz/dt. What is their difference?
ah Okay,

so its like this: dy/dt - dz/dt = x-y-(x-z)

dy/dt - dz/dt =-y+z
dy/dt - dz/dt +y-z =0
d/dt(y-z) + ( y-z) =0

dv/dt +v =0

I find v(t)= C * e-t
 
  • #18
masterchiefo said:
ah Okay,

so its like this: dy/dt - dz/dt = x-y-(x-z)

dy/dt - dz/dt =-y+z
dy/dt - dz/dt +y-z =0
d/dt(y-z) + ( y-z) =0

dv/dt +v =0

I find v(t)= C * e-t
Good. What were the initial conditions for y and z?
 
  • Like
Likes masterchiefo
  • #19
ehild said:
Good. What were the initial conditions for y and z?
y (0) = 0 and z ( 0) = 0
 
  • #20
masterchiefo said:
y (0) = 0 and z ( 0) = 0
What initial condition does it mean for v?
 
  • Like
Likes masterchiefo
  • #21
ehild said:
What initial condition does it mean for v?
v(t)=y(t)-z(t)

v(0)=0
 
  • #22
masterchiefo said:
v(t)=y(t)-z(t)

v(0)=0
So what is v(t)?
 
  • Like
Likes masterchiefo
  • #23
ehild said:
So what is v(t)?
v(t)= 0
y(t)-z(t)=0
y(t)=z(t)

You are awesome, thank you so much for the help.
Last thing, I want to push this problem a little bit.
now since we know that y(t) is equal to z(t), if I want to find the equation for them.

y(t)= 1 + t - x(t) - z(t)
z(t)= 1 + t - x(t) - y(t)

is this how I do it?
 
  • #24
masterchiefo said:
v(t)= 0
y(t)-z(t)=0
y(t)=z(t)

You are awesome, thank you so much for the help.
Last thing, I want to push this problem a little bit.
now since we know that y(t) is equal to z(t), if I want to find the equation for them.

y(t)= 1 + t - x(t) - z(t)
z(t)= 1 + t - x(t) - y(t)

is this how I do it?
Not quite.
You know that y(t)=z(t), and y(t)+z(t)=2y(t) =1+t-x(t). You also know x(t). What is y(t) then?
 
  • Like
Likes masterchiefo
  • #25
ehild said:
Not quite.
You know that y(t)=z(t), and y(t)+z(t)=2y(t) =1+t-x(t). You also know x(t). What is y(t) then?
x(t)= 1/9 *(4/e3t+3t+5)
2y(t) = 1+t- x(t)

2y(t) = 1+t-1/9 *(4/e3t+3t+5)

y(t) = (1+t-1/9 *(4/e3t+3t+5))/2
and
z(t) = (1+t-1/9 *(4/e3t+3t+5))/2

both same equation.
 
  • #26
masterchiefo said:
okay, 1 and 2 is done thank you very much. This problem confusing me a lot.

Now for Question 3:
.

There is some confusion here, so if anything is clear we should see it stated. If the differential equation for x is solved please tell us the solution.

Also because the rest is easy - once we know x the rest is just algebra. We have an algebraic relationship between x, y, z and t and also y=z as pointed out previously.
 
  • #27
masterchiefo said:
x(t)= 1/9 *(4/e3t+3t+5)
2y(t) = 1+t- x(t)

2y(t) = 1+t-1/9 *(4/e3t+3t+5)

y(t) = (1+t-1/9 *(4/e3t+3t+5))/2
and
z(t) = (1+t-1/9 *(4/e3t+3t+5))/2

both same equation.
It looks correct, but can be simplified a bit.
 
  • #28
epenguin said:
There is some confusion here, so if anything is clear we should see it stated. If the differential equation for x is solved please tell us the solution.

Also because the rest is easy - once we know x the rest is just algebra. We have an algebraic relationship between x, y, z and t and also y=z as pointed out previously.
x(t)= 1/9 *(4/e3t+3t+5)
 
  • #29
Is that same thing as x(t) = (1/9)* (4e-3t + 3t + 5) and have you checked it works?
 
  • #30
epenguin said:
Is that same thing as x(t) = (1/9)* (4e-3t + 3t + 5) and have you checked it works?
See post #25. It is the same.
 
  • Like
Likes masterchiefo
  • #31
Sorry.

Have you checked it? :oldbiggrin:

If it is right then the problem is solved for y and z too because

y(t) + z(t) = 1 + t - x(t)

and y = z at all times.
 
Last edited:
  • Like
Likes masterchiefo
  • #32
epenguin said:
Sorry.

Have you checked it? :oldbiggrin:

If it is right then the problem is solved for y and z too because

y(t) + z(t) = 1 + t - x(t)

and y = z at all times.
yes it works.
=)
 
  • #33
epenguin said:
Sorry.

Have you checked it? :oldbiggrin:
Of course, I always do.
epenguin said:
If it is right then the problem is solved for y and z too because

y(t) + z(t) = 1 + t - x(t)

and y = z at all times.
It has been solved by the OP, see Post #27
 
  • #34
ehild said:
It has been solved by the OP, see Post #27

Sorry I got lost on the way there.

Bit of a pity in such a long thread that the working of the essential bit was not shown. However I was trying to work it out not using integrating factors but a couple of other methods. Late because it was one of those days you may have had where it came out something like the OP's (correct) solution but not exactly right, t n finding the error the next attack came nearly but not quite right ands so on. In the end it is quite simple.

The first was the way it's usually taught (mostly for second order) for linear non-homogeneous d.e.'s with constant coefficients.

Denoting differentiation with ' to the equation

x' + 3x =. t. +. 2

following teaching, we assume a particular solution of form x = at + b. Then putting this into the above equation we get

a + 3(at + b) = 2 + t

From which we find we must have

a = 1/3 , b = 5/9

This is to be added to the general solution to

x' + 3x = 0

Which is x = Ae-3t

Into the solution

x = Ae-3t + t/3 + 5/9

We put x. =. x0 = 1, and find A = 4/9 so the solution for this initial condition is x = (4/9)e-3t + t/3 + 5/9

in accord with the OP. Slightly more suggestive is when initial condition is x0 = 0, and then we get

x = (5/9)e-3t + t/3I'll write up another method tomorrow.
 
Last edited:
  • #35
epenguin said:
Slightly more suggestive is when initial condition is x0 = 0, and then we get

x = (5/9)e-3t + t/3I'll write up another method tomorrow.
It is nice that you wrote one method to solve a linear first order differential equation in detail, and promised to give the other ones. These methods (mainly the integrating factor method) are taught during the Calculus courses everywhere. The OP was familiar with one of them. He did not write the process as he had no problem with it. I also think it is useful for the visitors of PF when they can see full solutions, but full solutions are forbidden even then, when the OP has solved the problem already. I have got warnings and infraction for less.
By the way, the equation to be solved was dx/dt + 3x = 2 + t, but you wrote
x' + x =. t. +. 2
, and your last equation does not correspond to x(0)=0.
 
Last edited:
  • #36
ehild said:
It is nice that you wrote one method to solve a linear first order differential equation in detail, and promised to give the other ones. These methods (mainly the integrating factor method) are taught during the Calculus courses everywhere. The OP was familiar with one of them. He did not write the process as he had no problem with it. I also think it is useful for the visitors of PF when they can see full solutions, but full solutions are forbidden even then, when the OP has solved the problem already. I have got warnings and infraction for less.
If the OP has already done the problem and gotten a solution, I don't have a problem with giving a full solution in that case, and I don't believe the other mentors who are currently active would, either.
 
  • #37
Mark44 said:
If the OP has already done the problem and gotten a solution, I don't have a problem with giving a full solution in that case, and I don't believe the other mentors who are currently active would, either.
It is a pleasure if that is allowed. Is it really? I am often inclined to show a simpler or nicer solution after the OP gave his one, but I have got warnings and even infraction for "almost" full solutions in the past.
 
  • #38
ehild said:
It is a pleasure if that is allowed. Is it really? I am often inclined to show a simpler or nicer solution after the OP gave his one, but I have got warnings and even infraction for "almost" full solutions in the past.

Here's the PF rule that applies:
Giving Full Answers: On helping with questions: Any and all assistance given to homework assignments or textbook style exercises should be given only after the questioner has shown some effort in solving the problem. If no attempt is made then the questioner should be asked to provide one before any assistance is given. Under no circumstances should complete solutions be provided to a questioner, whether or not an attempt has been made.
I believe the intent of the last sentence is that if the OP has provided at most only an attempt, a full solution should not be given. On the other hand, if the OP has shown a solution, it's OK for a poster to provide an alternative solution.

That's my take. I will bring it up with the mentors and see if they with me (or not). If they do, I'll look into changing that rule.
 
  • #39
Mark44 said:
Here's the PF rule that applies:

I believe the intent of the last sentence is that if the OP has provided at most only an attempt, a full solution should not be given. On the other hand, if the OP has shown a solution, it's OK for a poster to provide an alternative solution.

That's my take. I will bring it up with the mentors and see if they with me (or not). If they do, I'll look into changing that rule.
Thank you Mark. It would be nice if a thread could be finished with a comprehensive full solution. As it is now, the threads usually end when the OP got enough hint to solve the problem, and the process of solution is not shown, even by the OP.
 
  • #40
Yes it seems the rules allow detailed solution and discussion of a solved problem.

Especially in this case where the boot was rather on the other foot - if the OP had not posted the correct solution I would probably have posted one of my incorrect ones!

I say there is a sound educational or learning principle involved.
epenguin said:
And because of the Polya principle - when you've got an answer it isn't finished! They will miss concluding comments we might make. This solution is just an example of a wider principle, that technique can be applied to other things, there have been recent advances around the theme being treated... This question or subject matter is related to another.

This last is related thing could do a little to combat what is a vice made too prevalent just by our focussed question-and-answer thing that pervades education which misleads many students into seeing Science as just a fearful ritual :bow: in which you have to be able to recite the responses or quote the texts. Or see the world as divided into chapters and subjects, self-sufficient and unrelated.
...

Then we mustn't let happen what I have unfortunately seen a lot of in other worlds and situations - rules rigidly applied to the detriment of the forgotten objective they were invented for in the first place! :oldbiggrin:
 
Last edited:
  • #41
Firstly probably the most important thing to for anyone to realize straightaway about this problem is that without solving the d.e. It should be obvious enough that long-term, that is at high t, the solution has to be close to x = t/3 , in fact x = y = z = t/3 .

Then about the solution, I invite comments about this not-quite-the-same reasoning or 'method'.

Denoting differentiation w.r.t. t by ' sign, the equation is

x' + 3x = 2 + t

Change variable to X = x - 2/3 , equation is

X' + 3X = t

Differentiating,

X'' + 3X' = 1

Changing variable again to Y = X' - 1/3, this is the homogeneous linear d.e.

Y' + 3Y = 0

Y'/Y = 3

Solution

Y = Ae-3t

X' = Ae-3t + 1/3

X = Ae-3t + t/3 + K

x =. Ae-3t + t/3 + K'

and the arbitrary constants are worked out as before.

Comments:

  • in solving d.e.'s I've often enough found it useful to differentiate them (though usually thee are other methods too).
  • This gives nicely the most important part of the solution - the long term t/3 term
  • The virtue I see is that I/anyone does it from a starting point of greater ignorance, no need to know that if the 'forcing function' (RHS) is a polynomial you have to know or assume then'particular' solution is another polynomial
  • I think it works for any polynomial - just have to differentiate it enough times (n) to make a constant.
  • (I guess the first substitution was superfluous.)
 
Back
Top