Services & Resources / Wolfram Forums / MathGroup Archive
-----

MathGroup Archive 2008

[Date Index] [Thread Index] [Author Index]

Search the Archive

Partial differential equation with evolving boundary conditions

  • To: mathgroup at smc.vnet.net
  • Subject: [mg91477] Partial differential equation with evolving boundary conditions
  • From: "Ingolf Dahl" <ingolf.dahl at telia.com>
  • Date: Sat, 23 Aug 2008 01:41:41 -0400 (EDT)
  • Organization: Goteborg University
  • Reply-to: <ingolf.dahl at telia.com>

Best friends!
I am trying to solve a partial differential equation (in principle the heat
equation) with boundary conditions that also evolve with time. For instance
this works for me:
 
Manipulate[s=NDSolve[{\[Lambda];emax;dutycycle;\!\(
\*SubscriptBox[\(\[PartialD]\), \(y, y\)]\ \(\[Theta][y,
t]\)\)+elecy[t]*(5-\[Theta][y,t])*y*(1-y)==\[Lambda]*\!\(
\*SubscriptBox[\(\[PartialD]\), \(t\)]\(\[Theta][y, t]\)\),
\[Theta][y,0]==0,
-(\!\(
\*SubscriptBox[\(\[PartialD]\), \(y1\)]\(\[Theta][y1,
t]\)\)/.{y1->0})*0.1+\[Theta][0,t]==0,+( \!\(
\*SubscriptBox[\(\[PartialD]\), \(y1\)]\(\[Theta][y1,
t]\)\)/.{y1->1})*0.1+\[Theta][1,t]==
0},\[Theta],{y,0,1},{t,0,2}];Plot3D[\[Theta][y,t]/.s,{y,0,1},{t,0,2},PlotSty
le->Automatic,PlotRange->{0,5}],{{emax,25.,"emax"},0,100,Appearance->"Labele
d"},{{dutycycle,0.25,"dutycycle"},0,1,Appearance->"Labeled"},
{{\[Lambda],1.,"\[Lambda]"},0.001,10,Appearance->"Labeled"},Initialization:>
(elecy[t_]:=Which[0<=Mod[t,1]<=dutycycle,emax,dutycycle<Mod[t,1]<=1,0.]),
ContinuousAction->False,ControlPlacement->Top]

but this does not work
 
Manipulate[s=NDSolve[{\[Lambda];\[Lambda]s;emax;dutycycle;\!\(
\*SubscriptBox[\(\[PartialD]\), \(y, y\)]\ \(\[Theta][y,
t]\)\)+elecy[t]*(5-\[Theta][y,t])*y*(1-y)==\[Lambda]*\!\(
\*SubscriptBox[\(\[PartialD]\), \(t\)]\(\[Theta][y, t]\)\),
\[Theta][0,t]==-\[Lambda]s*\!\(
\*SubscriptBox[\(\[PartialD]\), \(t\)]\(\[Theta][0,
t]\)\),\[Theta][1,t]==-\[Lambda]s*\!\(
\*SubscriptBox[\(\[PartialD]\), \(t\)]\(\[Theta][1, t]\)\),
\[Theta][y,0]==0},\[Theta],{y,0,1},{t,0,2}];Plot3D[\[Theta][y,t]/.s,{y,0,1},
{t,0,2},PlotStyle->Automatic,PlotRange->{0,5}],{{emax,25.,"emax"},0,100,Appe
arance->"Labeled"},{{dutycycle,0.25,"dutycycle"},0,1,Appearance->"Labeled"},
{{\[Lambda],1.,"\[Lambda]"},0.001,10,Appearance->"Labeled"},{{\[Lambda]s,1.,
"\[Lambda]s"},0.001,10,Appearance->"Labeled"},Initialization:>(elecy[t_]:=Wh
ich[0<=Mod[t,1]<=dutycycle,emax,dutycycle<Mod[t,1]<=1,0.]),
ContinuousAction->False,ControlPlacement->Top]
 
I obtain the error message
 
NDSolve::bdord: Boundary condition \[Theta][0,t]+1. (\[Theta]^(0,1))[0,t]
should have derivatives of order lower than the differential order of the
partial differential equation. >>
 
plus a lot more. As I understand it, the second case should be solvable in
principle. Is this error a deficiency of NDSolve, or have I made some
mistake? This problem can be seen as a partial differential equation (in y
and t) coupled to two ordinary differential equations (for the boundary
conditions, only in t). What is the best way to solve such problems?
 
This is the simplification of a more involved case, where the boundary
conditions also are coupled to the partial differential equation, but I have
tried to boil down the problem here.
 
Best regards
 
Ingolf Dahl
Sweden




  • Prev by Date: Re: Integral of radial solution (hydrogen atom) is not evaluated
  • Next by Date: Re: Non Default stylesheet format and syntax coloring
  • Previous by thread: Re: Re: Help to remove equivalent (redundant) solutions
  • Next by thread: Re: Partial differential equation with evolving boundary conditions