[Date Index]
[Thread Index]
[Author Index]
Trouble solving system of first order, two variable PDEs in
*To*: mathgroup at smc.vnet.net
*Subject*: [mg93341] Trouble solving system of first order, two variable PDEs in
*From*: Tim Mc <tmccomb at creol.ucf.edu>
*Date*: Tue, 4 Nov 2008 06:18:52 -0500 (EST)
Hello All,
I am trying to solve what will become a large system of PDEs (20-30 equations) (though now it is only three representative equations) using the Mathematica NDSolve function. The equations all only involve first derivatives and are a function of one time and one space variable. I have inserted the code below. When I attempt to evaluate I get the error:
NDSolve::ivone: Boundary values may only be specified for one \
independent variable. Initial values may only be specified at one \
value of the other independent variable
It seems no matter what I do (even stripping off the numbers and breaking the equations down to their bare bones) I always get this error for these type of boundary conditions.
Can any one give me any advice on whether I am doing something wrong in the software or if my set of equations are simply incorrectly posed?
This is my Mathematica code:
\[Sigma]es = 3.4*^25;
\[Sigma]as = 5*^-27;
\[Sigma]ep = 5*^-26;
\[Sigma]ap = 8*^-25;
\[Tau] = 840*^-6;
\[Alpha] = 3*^-3;
Nn = 6*^25;
\[CapitalGamma] = 0.85;
\[CapitalGamma]p = 0.01;
\[Lambda]p = 0.915*^-6;
\[Lambda]s = 1.064*^-6;
A = 3*^-10;
h = 6.626*^-34;
c = 3*^8;
n = 1.5;
eq1 = D[N2[z, t],
t] == (\[CapitalGamma]p*\[Lambda]p/(h*c*A)) (\[Sigma]ap*(Nn -
N2[z, t]))*Pp[z, t] -
N2[z, t]/\[Tau] + (\[CapitalGamma]*\[Lambda]s/(h*c*
A)) (\[Sigma]as*(Nn - N2[z, t]) - \[Sigma]es*N2[z, t])*
Pp[z, t];
eq2 = D[Pp[z, t], z] + (n/c)*
D[Pp[z, t],
z] == -\[CapitalGamma]p (\[Sigma]ap*(Nn -
N2[z, t]) - \[Sigma]ep*N2[z, t])*Pp[z, t] - \[Alpha]*
Pp[z, t];
eq3 = D[Ps[z, t], z] + (n/c)*
D[Ps[z, t],
z] == \[CapitalGamma] (\[Sigma]es*
N2[z, t] - \[Sigma]as*(Nn - N2[z, t]))*Ps[z, t] - \[Alpha]*
Ps[z, t];
bc1 = Ps[0, t] == 1;
bc2 = Pp[0, t] == 20;
bc3 = N2[z, 0] == 0;
solution = {N2, Pp, Ps} /.
First[NDSolve[{eq1, eq2, eq3, bc1, bc2, bc3}, {N2, Pp, Ps}, {z, 0,
10}, {t, 0, 0.001}]]
The error I see:
NDSolve::ivone: Boundary values may only be specified for one \
independent variable. Initial values may only be specified at one \
value of the other independent variable. >>
Thank You in Advance,
Tim
Prev by Date:
**Re: Eliminating common factors?**
Next by Date:
**Re: Coefficient[ ] quirks**
Previous by thread:
**Re: Create a cumulative list**
Next by thread:
**strange effect: program works only when board clicked**
| |