MathGroup Archive 2009

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

Search the Archive

interpolating function as initial condition

  • To: mathgroup at smc.vnet.net
  • Subject: [mg102804] interpolating function as initial condition
  • From: venkat reddy <vemkatma at gmail.com>
  • Date: Thu, 27 Aug 2009 06:37:58 -0400 (EDT)

* hi all,

I am relatively new to mathematica. I am trying to solve  4 couples pdes in
which one is second order pde. What I like to do is solve the equation
for some time span(like t=0 to t=tfinal) with some arbitray initial
conditions(e.g. cf[0,z]=tf[0,z]=something). I get solution for my 4
dependent variables(cf,tf,cc,tc) in terms of interpolated functions. I want
to use these interpolated functions as initial condition(like
cf[o,z]=InterpolatingFunction for  cf at tfinal) for further calculations=
. I
have trouble in using the interpolating function as initial condition.

Please can some body help me out.

Thank you.

here is the code for your reference



sol=NDSolve[

{

D[cf[t,z],t]+v D[cf[t,z],z]-a1 (cc[t,z]-cf[t,z])=C5 0,

D[tf[t,z],t]+v D[tf[t,z],z]-a2 (tc[t,z]-tf[t,z])-a3 (tw-tf[t,z])=C5 0,

D[cc[t,z],t]-a4 (cf[t,z]-cc[t,z])-a5 r[cc[t,z],tc[t,z]]=C5 0,

D[tc[t,z],t]-a6 (tf[t,z]-tc[t,z])-a7 r[cc[t,z],tc[t,z]]-a8 D[tc[t,z],z,z]=
=C5 0,

cf[t,0]=C5 c0,tf[t,0]=C5 t0,(D[tc[t,z],z]/.z=C2=AE0)=C5 0,(D[tc[t,z],=
z]/.z=C2=AElength)=C5 0,

cf[0,z]=C5 c0,tf[0,z]=C5 t0,cc[0,z]=C5 c0,tc[0,z]=C5 t0

},

{cf,tf,cc,tc},{t,0,tfinal}, {z,0,length}
]

{{cf=C2=AEInterpolatingFunction[{{0.,10000.},{0.,1.}},<>],

tf=C2=AEInterpolatingFunction[{{0.,10000.},{0.,1.}},<>],

cc=C2=AEInterpolatingFunction[{{0.,10000.},{0.,1.}},<>],

tc=C2=AEInterpolatingFunction[{{0.,10000.},{0.,1.}},<>]}}
*



  • Prev by Date: Re: ClearAll ?? or what
  • Next by Date: Re: image from graphics3d - rid black line
  • Previous by thread: Re: Problems with functions of a Complex Variable
  • Next by thread: accumulate coefficients of a polynomial