Re: Convective diffusion equation in 2D
- To: mathgroup at smc.vnet.net
- Subject: [mg79770] Re: Convective diffusion equation in 2D
- From: Jens-Peer Kuska <kuska at informatik.uni-leipzig.de>
- Date: Sat, 4 Aug 2007 05:49:07 -0400 (EDT)
- References: <f8v09f$cbk$1@smc.vnet.net>
Hi,
With[{\[ScriptCapitalD] = 1/8, v = 1/4},
sol = NDSolve[{
D[u[x, y, t],
t] == \[ScriptCapitalD] (D[u[x, y, t], {x, 2}] +
D[u[x, y, t], {y, 2}]) - v*D[u[x, y, t], x],
u[-1, y, t] == 0, u[1, y, t] == 0,
u[x, -1, t] == 0, u[x, 1, t] == 0,
u[x, y, 0] == Piecewise[{{1, Sqrt[x^2 + y^2] <= 0.5}}, 0]},
u[x, y, t], {x, -1, 1}, {y, -1, 1}, {t, 0, 2}
]
]
may do that .. Ignore the warnings you get with the
statement above, the reason are the initial conditions
and you may have other.
Regards
Jens
dantimatter wrote:
> Hello all,
>
> I'm trying to find a nice and neat way to numerically solve the
> convective diffusion equation
>
> da/dt = D (d^2/dx^2 + d^2/dy^2) a - v da/dx
>
> where a is the concentration of my solute, D is the diffusion
> constant, and v is the surrounding fluid velocity in the x direction.
> I thought that there was a small chance that maybe someone else here
> has attempted something similar.
> Is it even possible to solve this equation? As always, any suggestions
> would be much appreciated.
>
> Cheers,
> Dan
>
>