dma...@mst.edu schrieb am Dienstag, 20. März 2018 um 06:23:09 UTC+1:

On Thursday, October 14, 2010 at 4:05:06 AM UTC-5, Daniel wrote:

Torsten Hennig <Torsten...@umsicht.fhg.de> wrote in message <1098368673.97608.12869...@gallium.mathforum.org>...

I don't understand why you make life so difficult.
For each node j in the electrode, define a
solution vector
(y_1^(j),...,y_6^(j),y_7^(j),...,y_10^(j))
and for each node j in the seperator, define a
solution vector
(y_1^(j),...,y_6^(j)).
For your above choice of the nuumber of elements,
glue the solution vectors together to an y-vector
of the form
y = (y_1^(1),...,y_10^(1),y_1^(2),...,y_10^(2), ..,y_1^(6),...,y_10^(6),y_1^(7),...,y_6^(7),...,
y_1^(10),...,y_6^(10))
with 60 + 24 = 84 components and call ODE15s
with this vector.
In the function routine, you can write the full
y-vector in local vectors which reflect the underlying
problem (e.g. it is convenient to define subvectors (y_1^(1),y_1^(2),...,y_1^(10)),
(y_2^(1),y_2^(2),...,y_2^(10)),
..
(y_6^(1),y_6^(2),...,y_6^(10)),
(y_7^(1),y_7^(2),...,y_7^(6)),
..
(y_10^(1),y_10^(2),...,y_10^(6))
for teh different solution variables, I think.)

Best wishes
Torsten.

Thanks! Tried that.
But now the DAE index appears to be greater than 1:
??? Error using ==> daeic12 at 77
This DAE appears to be of index greater than 1.
Don't know why exactly because all algebraic states can be described by the algebraic equations...

So again the question can Matlab solve the original PDE/ODE system without discretization?

Hello Daniel,
I by chance came across this thread on index error associated with modelling lithium ion battery cell using ode15s in MATLAB. If it will not be too much of a ask, how did you eventually solve this problem because I am presently having this same problem

with my battery model using spectral methods.

Regards,
Damola M. Ajiboye.

Same problem, also modeling Li-Ion batteries. Any news?

--- SoupGate-Win32 v1.05
* Origin: fsxNet Usenet Gateway (21:1/5)