Drexel dragonThe Math ForumDonate to the Math Forum



Search All of the Math Forum:

Views expressed in these public forums are not endorsed by Drexel University or The Math Forum.


Math Forum » Discussions » Software » comp.soft-sys.matlab

Topic: Pdepe
Replies: 5   Last Post: Mar 26, 2013 3:55 AM

Advanced Search

Back to Topic List Back to Topic List Jump to Tree View Jump to Tree View   Messages: [ Previous | Next ]
Torsten

Posts: 1,454
Registered: 11/8/10
Re: Pdepe
Posted: Mar 25, 2013 11:05 AM
  Click to see the message monospaced in plain text Plain Text   Click to reply to this topic Reply

"Igor" wrote in message <kipm7d$61d$1@newscl01ah.mathworks.com>...
> I cannot provide the whole code, just key components:
>
> I have system of 16 eq like this(only first rows) for steady state:
> .....
> function [c,f,s] = pdex14pde....
> c = ones(NumVars,1);%c
> f = [(D_bio+D_O2)*DuDx(1)-w*u(1);....
> s = [(BC0_O2-u(1))*alfax + (-0.25*R8-2* 1.5*R21 - 0.5*R22) - R1*F;... %w,R,alfax-some constant coef.
> ..........
> function u0 = pdex1ic(~)
> global NumVars;
> u0 = zeros(NumVars,1);%u0
> ...........
> function [pl,ql,pr,qr] = pdex1bc(~,ul,~,ur,~)
> pl = [ul(1)-BC0_O2;...
> ql = [0;1;....
> pr = [w*ur(1);...
> qr = ones....
>
> This code takes 300 sec for calculation.
>
> But then I make some BC0_OM1=(360*cos(t*(2*pi))+360) it takes hours.
>
> Maybe it will help,
> Igor


I doubt that a steady-state solution for your problem exists if a boundary condition depends on time.
Or did I misunderstand something ?

Best wishes
Torsten.


Date Subject Author
3/23/13
Read Re: Pdepe
Nasser Abbasi
3/25/13
Read Re: Pdepe
Torsten
3/26/13
Read Re: Pdepe
Torsten

Point your RSS reader here for a feed of the latest messages in this topic.

[Privacy Policy] [Terms of Use]

© Drexel University 1994-2014. All Rights Reserved.
The Math Forum is a research and educational enterprise of the Drexel University School of Education.