[Date Index]
[Thread Index]
[Author Index]
using interpolating function as initial condition to Ndsolve
*To*: mathgroup at smc.vnet.net
*Subject*: [mg102766] using interpolating function as initial condition to Ndsolve
*From*: Santhosh <santhosh833 at gmail.com>
*Date*: Wed, 26 Aug 2009 07:45:23 -0400 (EDT)
hi all,
I am solving a set of four pdes (first order in time and second order
in space) using NDsolve. I am able to solve for arbitrary initial
conditions. What I am interested is fist solve with some random
initial conditions and get the steady state solution. I want to use
this solution as initial condition(e.g cf[0,z]=f(z)....) for the next
computations. I could not figure out how to use interpolating function
as initial condition.
can anybody have any idea?
Thank you.
here is my code
sol = NDSolve[
{
D[cf[t, z], t] + v D[cf[t, z], z] - a1 (cc[t, z] - cf[t, z]) ==
0,
D[tf[t, z], t] + v D[tf[t, z], z] - a2 (tc[t, z] - tf[t, z]) -
a3 (tw - tf[t, z]) == 0,
D[cc[t, z], t] - a4 (cf[t, z] - cc[t, z]) -
a5 r[cc[t, z], tc[t, z]] == 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] == 0,
cf[t, 0] == c0,
tf[t, 0] == t0, (D[tc[t, z], z] /. z -> 0) ==
0, (D[tc[t, z], z] /. z -> length) == 0,
cf[0, z] == c0, tf[0, z] == t0, cc[0, z] == c0, tc[0, z] == t0
},
{cf, tf, cc, tc}, {t, 0, tfinal}, {z, 0, length}
]
Prev by Date:
**Re: Select text with keyboard?**
Next by Date:
**Re: Viewing packages in mathematica**
Previous by thread:
**Re: image from graphics3d - rid black line**
Next by thread:
**Problems with Evaluate[Symbol[xxx]]**
| |