Andew,
It seems like sign of the diffusoin term is wrong. It can't possible be right.
transportEq2=(TransientTerm()==
-PowerLawConvectionTerm(coeff=vel())
-DiffusionTerm(coeff=Disp))
Disp has a value of 0.0002. Is this correct?
On Jan 24, 2008 11:17 AM, asreeve <[EMAIL PROTECTED]> wrote:
> Daniel,
>
> I've attached the file. Its a very simple mesh...just one cell in each
> layer, so I could easily keep track of what's going on and check it
> against a 1-D analytic solution. I've attached both the mesh and geometry
> files.
>
> Andy
>
> ----------------------
> Andrew Reeve
> Dept. of Earth Science
> University of Maine
> Orono, ME 04469
>
>
> On Wed, 23 Jan 2008, Daniel Wheeler wrote:
>
> > Date: Wed, 23 Jan 2008 16:59:23 -0500
> > From: Daniel Wheeler <[EMAIL PROTECTED]>
> > Reply-To: [email protected]
> > To: Multiple recipients of list <[email protected]>
> > Subject: Re: using convection term
>
> >
> >
> > Hi Andrew,
> >
> > Can you attach './HalfCirc.msh' so I can run the script.
> >
> > Cheers
> >
--
Daniel Wheeler
