[CIG-CS] Particle problems and Semi-Lagrangian Schemes for Gamr

[Matthew Knepley]

On Tue, May 31, 2011 at 3:09 PM, Walter Landry <walter at geodynamics.org>wrote:

> Back to the list.
>
> Matthew Knepley <knepley at mcs.anl.gov> wrote:
> > Of course, you know I like this method. Its the "right thing to
> > do"TM However, I think the question to ask is:
> >
> >   "How do you guys check the solution?"
> >
> > If its something like conservation of mass fraction, you can assure
> > people that a semi-lagrangian method with an accurate time-stepper
> > will obey that.
>
> I think the only way to solve the problems that I detailed in my
> previous email ("Gamr status") is to track interfaces (unless you have
> a better idea?).  If we track interfaces well, then we get mass
> conservation.  Particles or semi-lagrangian would be used to track
> material properties like integrated strain or elastic stress.  So, for
> example, to reproduce what CitcomS does, neither is required.
>
> However, we would still need particles for tracers.  But those are not
> used when computing solutions, so it is not necessary to have such a
> dense array.


I would not do this. As long as you are using SL, I would use it all the way
and eliminate particles, making the code much simpler, especially in
parallel.

   Matt


> Cheers,
> Walter Landry
> walter at geodynamics.org
> _______________________________________________
> CIG-CS mailing list
> CIG-CS at geodynamics.org
> http://geodynamics.org/cgi-bin/mailman/listinfo/cig-cs
>



-- 
What most experimenters take for granted before they begin their experiments
is infinitely more interesting than any results to which their experiments
lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://geodynamics.org/pipermail/cig-cs/attachments/20110531/c99a7864/attachment.htm