System of linear and nonlinear diff equations [duplicate]












-1












$begingroup$



This question already has an answer here:




  • System of differential equations - linear system, shooting method

    2 answers




I am looking for a help with solving this chemical-engineering problem. Is there anybody, who can help?



The task is about this:
Concentration profiles $cA,dA$ during extraction of A-component from solvent C into solvent B in countercurrent colony with axial mixing of continual phase of solvent B are described with following diff equations which I have to solve:
begin{align}
dA'' &= -(uB/E), dA'-kL ,(cA-dA/psi)\
cA' &= -(kL/uC), (cA-dA/psi),
end{align}



where $cA,dA$ represent concentrations of A-component in original and new solvent, $kL$,$uB$,$uC$,$E$ are constants, "psi" is partition coefficient. Then we have following initial and boundary conditions - z=0 m is inlet of feedstock, z=2 m is inlet of extraction agent.





  • $cA(0) = 100$


  • $uB* ,dA(2)+E,* dA'(2) = 0$


  • $dA'(0) = 0$


Firstly, I have to solve it for $psi = 100$ (at this moment system of diff equations is linear) and find value for $cA(2)$ in raffinate then. Then I have to calculate using shooting method and then use proper Runge-Kutta approximation (using $psi = 10*sqrt(dA)$, for non-physical option $dA$ $=<0$ I have to use $dA/psi = 0$).



Firstly I want understand the problem in general - the final numbers are not important. I want to try it using your advices (or some similar problems which I can use in my solution) and how to use the substition for the diff equation with rate 2, then how to calculate three diff equations with rate 1 with these written conditions for both cases - linear and nonlinear.



Thank you, J.










share|cite|improve this question









$endgroup$



marked as duplicate by Did, Nosrati, Cesareo, Shailesh, Leucippus Jan 21 at 3:57


This question has been asked before and already has an answer. If those answers do not fully address your question, please ask a new question.


















  • $begingroup$
    This, System of differential equations - linear system, shooting method, looks similar. Is this a duplicate?
    $endgroup$
    – LutzL
    Jan 14 at 20:09












  • $begingroup$
    Yes, its the same, but this version should be more clear for understanding. And I have not understand the comments below the task on that link..
    $endgroup$
    – user634771
    Jan 14 at 20:12












  • $begingroup$
    The comments were about a typographical problem. You should be able to understand the answers. What tools do you have available? You need a numerical ODE solver and a routine to find roots of scalar functions.
    $endgroup$
    – LutzL
    Jan 14 at 20:27










  • $begingroup$
    Maple. But for the analytic solution, I should do it on paper, without any software.
    $endgroup$
    – user634771
    Jan 14 at 20:29






  • 1




    $begingroup$
    The order-1 system is also nicely explained in the answer of Christoph in the other post.
    $endgroup$
    – LutzL
    Jan 14 at 20:47
















-1












$begingroup$



This question already has an answer here:




  • System of differential equations - linear system, shooting method

    2 answers




I am looking for a help with solving this chemical-engineering problem. Is there anybody, who can help?



The task is about this:
Concentration profiles $cA,dA$ during extraction of A-component from solvent C into solvent B in countercurrent colony with axial mixing of continual phase of solvent B are described with following diff equations which I have to solve:
begin{align}
dA'' &= -(uB/E), dA'-kL ,(cA-dA/psi)\
cA' &= -(kL/uC), (cA-dA/psi),
end{align}



where $cA,dA$ represent concentrations of A-component in original and new solvent, $kL$,$uB$,$uC$,$E$ are constants, "psi" is partition coefficient. Then we have following initial and boundary conditions - z=0 m is inlet of feedstock, z=2 m is inlet of extraction agent.





  • $cA(0) = 100$


  • $uB* ,dA(2)+E,* dA'(2) = 0$


  • $dA'(0) = 0$


Firstly, I have to solve it for $psi = 100$ (at this moment system of diff equations is linear) and find value for $cA(2)$ in raffinate then. Then I have to calculate using shooting method and then use proper Runge-Kutta approximation (using $psi = 10*sqrt(dA)$, for non-physical option $dA$ $=<0$ I have to use $dA/psi = 0$).



Firstly I want understand the problem in general - the final numbers are not important. I want to try it using your advices (or some similar problems which I can use in my solution) and how to use the substition for the diff equation with rate 2, then how to calculate three diff equations with rate 1 with these written conditions for both cases - linear and nonlinear.



Thank you, J.










share|cite|improve this question









$endgroup$



marked as duplicate by Did, Nosrati, Cesareo, Shailesh, Leucippus Jan 21 at 3:57


This question has been asked before and already has an answer. If those answers do not fully address your question, please ask a new question.


















  • $begingroup$
    This, System of differential equations - linear system, shooting method, looks similar. Is this a duplicate?
    $endgroup$
    – LutzL
    Jan 14 at 20:09












  • $begingroup$
    Yes, its the same, but this version should be more clear for understanding. And I have not understand the comments below the task on that link..
    $endgroup$
    – user634771
    Jan 14 at 20:12












  • $begingroup$
    The comments were about a typographical problem. You should be able to understand the answers. What tools do you have available? You need a numerical ODE solver and a routine to find roots of scalar functions.
    $endgroup$
    – LutzL
    Jan 14 at 20:27










  • $begingroup$
    Maple. But for the analytic solution, I should do it on paper, without any software.
    $endgroup$
    – user634771
    Jan 14 at 20:29






  • 1




    $begingroup$
    The order-1 system is also nicely explained in the answer of Christoph in the other post.
    $endgroup$
    – LutzL
    Jan 14 at 20:47














-1












-1








-1





$begingroup$



This question already has an answer here:




  • System of differential equations - linear system, shooting method

    2 answers




I am looking for a help with solving this chemical-engineering problem. Is there anybody, who can help?



The task is about this:
Concentration profiles $cA,dA$ during extraction of A-component from solvent C into solvent B in countercurrent colony with axial mixing of continual phase of solvent B are described with following diff equations which I have to solve:
begin{align}
dA'' &= -(uB/E), dA'-kL ,(cA-dA/psi)\
cA' &= -(kL/uC), (cA-dA/psi),
end{align}



where $cA,dA$ represent concentrations of A-component in original and new solvent, $kL$,$uB$,$uC$,$E$ are constants, "psi" is partition coefficient. Then we have following initial and boundary conditions - z=0 m is inlet of feedstock, z=2 m is inlet of extraction agent.





  • $cA(0) = 100$


  • $uB* ,dA(2)+E,* dA'(2) = 0$


  • $dA'(0) = 0$


Firstly, I have to solve it for $psi = 100$ (at this moment system of diff equations is linear) and find value for $cA(2)$ in raffinate then. Then I have to calculate using shooting method and then use proper Runge-Kutta approximation (using $psi = 10*sqrt(dA)$, for non-physical option $dA$ $=<0$ I have to use $dA/psi = 0$).



Firstly I want understand the problem in general - the final numbers are not important. I want to try it using your advices (or some similar problems which I can use in my solution) and how to use the substition for the diff equation with rate 2, then how to calculate three diff equations with rate 1 with these written conditions for both cases - linear and nonlinear.



Thank you, J.










share|cite|improve this question









$endgroup$





This question already has an answer here:




  • System of differential equations - linear system, shooting method

    2 answers




I am looking for a help with solving this chemical-engineering problem. Is there anybody, who can help?



The task is about this:
Concentration profiles $cA,dA$ during extraction of A-component from solvent C into solvent B in countercurrent colony with axial mixing of continual phase of solvent B are described with following diff equations which I have to solve:
begin{align}
dA'' &= -(uB/E), dA'-kL ,(cA-dA/psi)\
cA' &= -(kL/uC), (cA-dA/psi),
end{align}



where $cA,dA$ represent concentrations of A-component in original and new solvent, $kL$,$uB$,$uC$,$E$ are constants, "psi" is partition coefficient. Then we have following initial and boundary conditions - z=0 m is inlet of feedstock, z=2 m is inlet of extraction agent.





  • $cA(0) = 100$


  • $uB* ,dA(2)+E,* dA'(2) = 0$


  • $dA'(0) = 0$


Firstly, I have to solve it for $psi = 100$ (at this moment system of diff equations is linear) and find value for $cA(2)$ in raffinate then. Then I have to calculate using shooting method and then use proper Runge-Kutta approximation (using $psi = 10*sqrt(dA)$, for non-physical option $dA$ $=<0$ I have to use $dA/psi = 0$).



Firstly I want understand the problem in general - the final numbers are not important. I want to try it using your advices (or some similar problems which I can use in my solution) and how to use the substition for the diff equation with rate 2, then how to calculate three diff equations with rate 1 with these written conditions for both cases - linear and nonlinear.



Thank you, J.





This question already has an answer here:




  • System of differential equations - linear system, shooting method

    2 answers








ordinary-differential-equations






share|cite|improve this question













share|cite|improve this question











share|cite|improve this question




share|cite|improve this question










asked Jan 14 at 19:33









user634771user634771

1




1




marked as duplicate by Did, Nosrati, Cesareo, Shailesh, Leucippus Jan 21 at 3:57


This question has been asked before and already has an answer. If those answers do not fully address your question, please ask a new question.









marked as duplicate by Did, Nosrati, Cesareo, Shailesh, Leucippus Jan 21 at 3:57


This question has been asked before and already has an answer. If those answers do not fully address your question, please ask a new question.














  • $begingroup$
    This, System of differential equations - linear system, shooting method, looks similar. Is this a duplicate?
    $endgroup$
    – LutzL
    Jan 14 at 20:09












  • $begingroup$
    Yes, its the same, but this version should be more clear for understanding. And I have not understand the comments below the task on that link..
    $endgroup$
    – user634771
    Jan 14 at 20:12












  • $begingroup$
    The comments were about a typographical problem. You should be able to understand the answers. What tools do you have available? You need a numerical ODE solver and a routine to find roots of scalar functions.
    $endgroup$
    – LutzL
    Jan 14 at 20:27










  • $begingroup$
    Maple. But for the analytic solution, I should do it on paper, without any software.
    $endgroup$
    – user634771
    Jan 14 at 20:29






  • 1




    $begingroup$
    The order-1 system is also nicely explained in the answer of Christoph in the other post.
    $endgroup$
    – LutzL
    Jan 14 at 20:47


















  • $begingroup$
    This, System of differential equations - linear system, shooting method, looks similar. Is this a duplicate?
    $endgroup$
    – LutzL
    Jan 14 at 20:09












  • $begingroup$
    Yes, its the same, but this version should be more clear for understanding. And I have not understand the comments below the task on that link..
    $endgroup$
    – user634771
    Jan 14 at 20:12












  • $begingroup$
    The comments were about a typographical problem. You should be able to understand the answers. What tools do you have available? You need a numerical ODE solver and a routine to find roots of scalar functions.
    $endgroup$
    – LutzL
    Jan 14 at 20:27










  • $begingroup$
    Maple. But for the analytic solution, I should do it on paper, without any software.
    $endgroup$
    – user634771
    Jan 14 at 20:29






  • 1




    $begingroup$
    The order-1 system is also nicely explained in the answer of Christoph in the other post.
    $endgroup$
    – LutzL
    Jan 14 at 20:47
















$begingroup$
This, System of differential equations - linear system, shooting method, looks similar. Is this a duplicate?
$endgroup$
– LutzL
Jan 14 at 20:09






$begingroup$
This, System of differential equations - linear system, shooting method, looks similar. Is this a duplicate?
$endgroup$
– LutzL
Jan 14 at 20:09














$begingroup$
Yes, its the same, but this version should be more clear for understanding. And I have not understand the comments below the task on that link..
$endgroup$
– user634771
Jan 14 at 20:12






$begingroup$
Yes, its the same, but this version should be more clear for understanding. And I have not understand the comments below the task on that link..
$endgroup$
– user634771
Jan 14 at 20:12














$begingroup$
The comments were about a typographical problem. You should be able to understand the answers. What tools do you have available? You need a numerical ODE solver and a routine to find roots of scalar functions.
$endgroup$
– LutzL
Jan 14 at 20:27




$begingroup$
The comments were about a typographical problem. You should be able to understand the answers. What tools do you have available? You need a numerical ODE solver and a routine to find roots of scalar functions.
$endgroup$
– LutzL
Jan 14 at 20:27












$begingroup$
Maple. But for the analytic solution, I should do it on paper, without any software.
$endgroup$
– user634771
Jan 14 at 20:29




$begingroup$
Maple. But for the analytic solution, I should do it on paper, without any software.
$endgroup$
– user634771
Jan 14 at 20:29




1




1




$begingroup$
The order-1 system is also nicely explained in the answer of Christoph in the other post.
$endgroup$
– LutzL
Jan 14 at 20:47




$begingroup$
The order-1 system is also nicely explained in the answer of Christoph in the other post.
$endgroup$
– LutzL
Jan 14 at 20:47










1 Answer
1






active

oldest

votes


















0












$begingroup$

Hint.



For $psi = C^{te}$. The linear DE system has the structure



$$
dA''+a* dA'+b* cA-b *dA/psi = 0\
cA' +c* cA-c*dA/psi=0
$$



or



$$
left(partial^2+apartial-frac{b}{psi}right)dA =-b*cA\
left(partial +cright)cA=c*dA/psi
$$



or



$$
left(partial +cright)left(partial^2+apartial-frac{b}{psi}right)dA +b*c*dA/psi = 0
$$



or



$$
left(left(partial +cright)left(partial^2+apartial-frac{b}{psi}right)+frac{b*c}{psi}right)dA=0
$$



so



$$
dA = C_1 e^{r_1 t}+C_2 e^{r_2 t}+C_3 e^{r_3 t}
$$



where $r_k, k = 1,2,3$ are the differential operator characteristic values (roots)



and then



$$
cA = -frac 1bleft(partial^2+apartial-frac{b}{psi}right)dA
$$



so now we have $dA, cA$ with three constants to determine, ($C_k$) using the boundary conditions. For the linear DE system no shooting method is needed.



NOTE



Applying the differential operator on $dA$ we get at



$$
dA = C_1 e^{r_1 t}+C_2 e^{r_2 t}+C_3 e^{r_3 t}\
cA = frac{left(b-psi r_1 left(a+r_1right)right)C_1}{b psi }e^{r_1 t}+frac{left(b-psi r_2 left(a+r_2right)right)C_2}{b psi }e^{r_2 t}+frac{left(b-psi r_3 left(a+r_3right)right)C_3}{b psi }e^{r_3 t}
$$



now with the boundary conditions we have



$$
cA(0) = frac{left(b-psi r_1 left(a+r_1right)right)C_1}{b psi }+frac{left(b-psi r_2 left(a+r_2right)right)C_2}{b psi }+frac{left(b-psi r_3 left(a+r_3right)right)C_3}{b psi }=100\
uB dA(2)+E dA'(2) = uB left(C_1 e^{2 r_1}+C_2 e^{2 r_2}+C_3 e^{2 r_3}right)+E left(C_1 e^{2 r_1} r_1+C_2 e^{2 r_2} r_2+C_3
e^{2 r_3} r_3right) = 0\
dA'(0) = C_1r_1+C_2r_2+C_3r_3 = 0
$$



so we have a linear system to solve for $(C_1,C_2,C_3)$






share|cite|improve this answer











$endgroup$













  • $begingroup$
    I understood the first part of yours, but still I have not deal with founding of the boundary condition that I need...
    $endgroup$
    – user634771
    Jan 20 at 8:15










  • $begingroup$
    @user634771 Please. See note attached.
    $endgroup$
    – Cesareo
    Jan 20 at 10:52










  • $begingroup$
    In general, I think I get it, but can you tell me what does it mean using this symbol "ϕ". Thanks
    $endgroup$
    – user634771
    Jan 20 at 11:19










  • $begingroup$
    @user634771 Sorry. Is is a typo now corrected.
    $endgroup$
    – Cesareo
    Jan 20 at 13:12


















1 Answer
1






active

oldest

votes








1 Answer
1






active

oldest

votes









active

oldest

votes






active

oldest

votes









0












$begingroup$

Hint.



For $psi = C^{te}$. The linear DE system has the structure



$$
dA''+a* dA'+b* cA-b *dA/psi = 0\
cA' +c* cA-c*dA/psi=0
$$



or



$$
left(partial^2+apartial-frac{b}{psi}right)dA =-b*cA\
left(partial +cright)cA=c*dA/psi
$$



or



$$
left(partial +cright)left(partial^2+apartial-frac{b}{psi}right)dA +b*c*dA/psi = 0
$$



or



$$
left(left(partial +cright)left(partial^2+apartial-frac{b}{psi}right)+frac{b*c}{psi}right)dA=0
$$



so



$$
dA = C_1 e^{r_1 t}+C_2 e^{r_2 t}+C_3 e^{r_3 t}
$$



where $r_k, k = 1,2,3$ are the differential operator characteristic values (roots)



and then



$$
cA = -frac 1bleft(partial^2+apartial-frac{b}{psi}right)dA
$$



so now we have $dA, cA$ with three constants to determine, ($C_k$) using the boundary conditions. For the linear DE system no shooting method is needed.



NOTE



Applying the differential operator on $dA$ we get at



$$
dA = C_1 e^{r_1 t}+C_2 e^{r_2 t}+C_3 e^{r_3 t}\
cA = frac{left(b-psi r_1 left(a+r_1right)right)C_1}{b psi }e^{r_1 t}+frac{left(b-psi r_2 left(a+r_2right)right)C_2}{b psi }e^{r_2 t}+frac{left(b-psi r_3 left(a+r_3right)right)C_3}{b psi }e^{r_3 t}
$$



now with the boundary conditions we have



$$
cA(0) = frac{left(b-psi r_1 left(a+r_1right)right)C_1}{b psi }+frac{left(b-psi r_2 left(a+r_2right)right)C_2}{b psi }+frac{left(b-psi r_3 left(a+r_3right)right)C_3}{b psi }=100\
uB dA(2)+E dA'(2) = uB left(C_1 e^{2 r_1}+C_2 e^{2 r_2}+C_3 e^{2 r_3}right)+E left(C_1 e^{2 r_1} r_1+C_2 e^{2 r_2} r_2+C_3
e^{2 r_3} r_3right) = 0\
dA'(0) = C_1r_1+C_2r_2+C_3r_3 = 0
$$



so we have a linear system to solve for $(C_1,C_2,C_3)$






share|cite|improve this answer











$endgroup$













  • $begingroup$
    I understood the first part of yours, but still I have not deal with founding of the boundary condition that I need...
    $endgroup$
    – user634771
    Jan 20 at 8:15










  • $begingroup$
    @user634771 Please. See note attached.
    $endgroup$
    – Cesareo
    Jan 20 at 10:52










  • $begingroup$
    In general, I think I get it, but can you tell me what does it mean using this symbol "ϕ". Thanks
    $endgroup$
    – user634771
    Jan 20 at 11:19










  • $begingroup$
    @user634771 Sorry. Is is a typo now corrected.
    $endgroup$
    – Cesareo
    Jan 20 at 13:12
















0












$begingroup$

Hint.



For $psi = C^{te}$. The linear DE system has the structure



$$
dA''+a* dA'+b* cA-b *dA/psi = 0\
cA' +c* cA-c*dA/psi=0
$$



or



$$
left(partial^2+apartial-frac{b}{psi}right)dA =-b*cA\
left(partial +cright)cA=c*dA/psi
$$



or



$$
left(partial +cright)left(partial^2+apartial-frac{b}{psi}right)dA +b*c*dA/psi = 0
$$



or



$$
left(left(partial +cright)left(partial^2+apartial-frac{b}{psi}right)+frac{b*c}{psi}right)dA=0
$$



so



$$
dA = C_1 e^{r_1 t}+C_2 e^{r_2 t}+C_3 e^{r_3 t}
$$



where $r_k, k = 1,2,3$ are the differential operator characteristic values (roots)



and then



$$
cA = -frac 1bleft(partial^2+apartial-frac{b}{psi}right)dA
$$



so now we have $dA, cA$ with three constants to determine, ($C_k$) using the boundary conditions. For the linear DE system no shooting method is needed.



NOTE



Applying the differential operator on $dA$ we get at



$$
dA = C_1 e^{r_1 t}+C_2 e^{r_2 t}+C_3 e^{r_3 t}\
cA = frac{left(b-psi r_1 left(a+r_1right)right)C_1}{b psi }e^{r_1 t}+frac{left(b-psi r_2 left(a+r_2right)right)C_2}{b psi }e^{r_2 t}+frac{left(b-psi r_3 left(a+r_3right)right)C_3}{b psi }e^{r_3 t}
$$



now with the boundary conditions we have



$$
cA(0) = frac{left(b-psi r_1 left(a+r_1right)right)C_1}{b psi }+frac{left(b-psi r_2 left(a+r_2right)right)C_2}{b psi }+frac{left(b-psi r_3 left(a+r_3right)right)C_3}{b psi }=100\
uB dA(2)+E dA'(2) = uB left(C_1 e^{2 r_1}+C_2 e^{2 r_2}+C_3 e^{2 r_3}right)+E left(C_1 e^{2 r_1} r_1+C_2 e^{2 r_2} r_2+C_3
e^{2 r_3} r_3right) = 0\
dA'(0) = C_1r_1+C_2r_2+C_3r_3 = 0
$$



so we have a linear system to solve for $(C_1,C_2,C_3)$






share|cite|improve this answer











$endgroup$













  • $begingroup$
    I understood the first part of yours, but still I have not deal with founding of the boundary condition that I need...
    $endgroup$
    – user634771
    Jan 20 at 8:15










  • $begingroup$
    @user634771 Please. See note attached.
    $endgroup$
    – Cesareo
    Jan 20 at 10:52










  • $begingroup$
    In general, I think I get it, but can you tell me what does it mean using this symbol "ϕ". Thanks
    $endgroup$
    – user634771
    Jan 20 at 11:19










  • $begingroup$
    @user634771 Sorry. Is is a typo now corrected.
    $endgroup$
    – Cesareo
    Jan 20 at 13:12














0












0








0





$begingroup$

Hint.



For $psi = C^{te}$. The linear DE system has the structure



$$
dA''+a* dA'+b* cA-b *dA/psi = 0\
cA' +c* cA-c*dA/psi=0
$$



or



$$
left(partial^2+apartial-frac{b}{psi}right)dA =-b*cA\
left(partial +cright)cA=c*dA/psi
$$



or



$$
left(partial +cright)left(partial^2+apartial-frac{b}{psi}right)dA +b*c*dA/psi = 0
$$



or



$$
left(left(partial +cright)left(partial^2+apartial-frac{b}{psi}right)+frac{b*c}{psi}right)dA=0
$$



so



$$
dA = C_1 e^{r_1 t}+C_2 e^{r_2 t}+C_3 e^{r_3 t}
$$



where $r_k, k = 1,2,3$ are the differential operator characteristic values (roots)



and then



$$
cA = -frac 1bleft(partial^2+apartial-frac{b}{psi}right)dA
$$



so now we have $dA, cA$ with three constants to determine, ($C_k$) using the boundary conditions. For the linear DE system no shooting method is needed.



NOTE



Applying the differential operator on $dA$ we get at



$$
dA = C_1 e^{r_1 t}+C_2 e^{r_2 t}+C_3 e^{r_3 t}\
cA = frac{left(b-psi r_1 left(a+r_1right)right)C_1}{b psi }e^{r_1 t}+frac{left(b-psi r_2 left(a+r_2right)right)C_2}{b psi }e^{r_2 t}+frac{left(b-psi r_3 left(a+r_3right)right)C_3}{b psi }e^{r_3 t}
$$



now with the boundary conditions we have



$$
cA(0) = frac{left(b-psi r_1 left(a+r_1right)right)C_1}{b psi }+frac{left(b-psi r_2 left(a+r_2right)right)C_2}{b psi }+frac{left(b-psi r_3 left(a+r_3right)right)C_3}{b psi }=100\
uB dA(2)+E dA'(2) = uB left(C_1 e^{2 r_1}+C_2 e^{2 r_2}+C_3 e^{2 r_3}right)+E left(C_1 e^{2 r_1} r_1+C_2 e^{2 r_2} r_2+C_3
e^{2 r_3} r_3right) = 0\
dA'(0) = C_1r_1+C_2r_2+C_3r_3 = 0
$$



so we have a linear system to solve for $(C_1,C_2,C_3)$






share|cite|improve this answer











$endgroup$



Hint.



For $psi = C^{te}$. The linear DE system has the structure



$$
dA''+a* dA'+b* cA-b *dA/psi = 0\
cA' +c* cA-c*dA/psi=0
$$



or



$$
left(partial^2+apartial-frac{b}{psi}right)dA =-b*cA\
left(partial +cright)cA=c*dA/psi
$$



or



$$
left(partial +cright)left(partial^2+apartial-frac{b}{psi}right)dA +b*c*dA/psi = 0
$$



or



$$
left(left(partial +cright)left(partial^2+apartial-frac{b}{psi}right)+frac{b*c}{psi}right)dA=0
$$



so



$$
dA = C_1 e^{r_1 t}+C_2 e^{r_2 t}+C_3 e^{r_3 t}
$$



where $r_k, k = 1,2,3$ are the differential operator characteristic values (roots)



and then



$$
cA = -frac 1bleft(partial^2+apartial-frac{b}{psi}right)dA
$$



so now we have $dA, cA$ with three constants to determine, ($C_k$) using the boundary conditions. For the linear DE system no shooting method is needed.



NOTE



Applying the differential operator on $dA$ we get at



$$
dA = C_1 e^{r_1 t}+C_2 e^{r_2 t}+C_3 e^{r_3 t}\
cA = frac{left(b-psi r_1 left(a+r_1right)right)C_1}{b psi }e^{r_1 t}+frac{left(b-psi r_2 left(a+r_2right)right)C_2}{b psi }e^{r_2 t}+frac{left(b-psi r_3 left(a+r_3right)right)C_3}{b psi }e^{r_3 t}
$$



now with the boundary conditions we have



$$
cA(0) = frac{left(b-psi r_1 left(a+r_1right)right)C_1}{b psi }+frac{left(b-psi r_2 left(a+r_2right)right)C_2}{b psi }+frac{left(b-psi r_3 left(a+r_3right)right)C_3}{b psi }=100\
uB dA(2)+E dA'(2) = uB left(C_1 e^{2 r_1}+C_2 e^{2 r_2}+C_3 e^{2 r_3}right)+E left(C_1 e^{2 r_1} r_1+C_2 e^{2 r_2} r_2+C_3
e^{2 r_3} r_3right) = 0\
dA'(0) = C_1r_1+C_2r_2+C_3r_3 = 0
$$



so we have a linear system to solve for $(C_1,C_2,C_3)$







share|cite|improve this answer














share|cite|improve this answer



share|cite|improve this answer








edited Jan 20 at 13:11

























answered Jan 14 at 21:20









CesareoCesareo

8,8993516




8,8993516












  • $begingroup$
    I understood the first part of yours, but still I have not deal with founding of the boundary condition that I need...
    $endgroup$
    – user634771
    Jan 20 at 8:15










  • $begingroup$
    @user634771 Please. See note attached.
    $endgroup$
    – Cesareo
    Jan 20 at 10:52










  • $begingroup$
    In general, I think I get it, but can you tell me what does it mean using this symbol "ϕ". Thanks
    $endgroup$
    – user634771
    Jan 20 at 11:19










  • $begingroup$
    @user634771 Sorry. Is is a typo now corrected.
    $endgroup$
    – Cesareo
    Jan 20 at 13:12


















  • $begingroup$
    I understood the first part of yours, but still I have not deal with founding of the boundary condition that I need...
    $endgroup$
    – user634771
    Jan 20 at 8:15










  • $begingroup$
    @user634771 Please. See note attached.
    $endgroup$
    – Cesareo
    Jan 20 at 10:52










  • $begingroup$
    In general, I think I get it, but can you tell me what does it mean using this symbol "ϕ". Thanks
    $endgroup$
    – user634771
    Jan 20 at 11:19










  • $begingroup$
    @user634771 Sorry. Is is a typo now corrected.
    $endgroup$
    – Cesareo
    Jan 20 at 13:12
















$begingroup$
I understood the first part of yours, but still I have not deal with founding of the boundary condition that I need...
$endgroup$
– user634771
Jan 20 at 8:15




$begingroup$
I understood the first part of yours, but still I have not deal with founding of the boundary condition that I need...
$endgroup$
– user634771
Jan 20 at 8:15












$begingroup$
@user634771 Please. See note attached.
$endgroup$
– Cesareo
Jan 20 at 10:52




$begingroup$
@user634771 Please. See note attached.
$endgroup$
– Cesareo
Jan 20 at 10:52












$begingroup$
In general, I think I get it, but can you tell me what does it mean using this symbol "ϕ". Thanks
$endgroup$
– user634771
Jan 20 at 11:19




$begingroup$
In general, I think I get it, but can you tell me what does it mean using this symbol "ϕ". Thanks
$endgroup$
– user634771
Jan 20 at 11:19












$begingroup$
@user634771 Sorry. Is is a typo now corrected.
$endgroup$
– Cesareo
Jan 20 at 13:12




$begingroup$
@user634771 Sorry. Is is a typo now corrected.
$endgroup$
– Cesareo
Jan 20 at 13:12



Popular posts from this blog

'app-layout' is not a known element: how to share Component with different Modules

android studio warns about leanback feature tag usage required on manifest while using Unity exported app?

WPF add header to Image with URL pettitions [duplicate]