[CIG-SHORT] Uniform time-stepping for long runs on viscoelastic domain (Brad Aagaard)

Matthew Knepley knepley at mcs.anl.gov
Mon Jun 8 04:16:19 PDT 2015


On Mon, Jun 8, 2015 at 1:23 AM, Michal Dichter <mdichter at unm.edu> wrote:

>  Dear Brad,
>
> Please disregard my last response. It seems that the cause of the problem
> that I have been describing has been with the nondimensionalization of
> time. (I stumbled upon a PyLith FAQ sort-of page where the importance of
> nondimensionalization was emphasized.) I changed my "relaxation_time"
> parameter to be equal to dt rather than to the minimum relaxation time.
> This enabled my code to run to completion and output the expected number of
> files. Thank you again for your help in this matter.
>
> PyLith developers,

I think some sort of diagnostic thing, even if only in the log, for looking
at the non-dimensionalization would be useful. We have
been struggling with this in PETSc as well. What if we flag things that are
1e6 times more or less than the non-dimensionalization?

  Thanks,

     Matt

> All the best,
>
>
>  Michal
>
>
>   ____________
> Michal Dichter
> Department of Physics and Astronomy
> University of New Mexico
> (860) 459-0956
>
> _______________________________________________
> CIG-SHORT mailing list
> CIG-SHORT at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/cig-short
>



-- 
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://lists.geodynamics.org/pipermail/cig-short/attachments/20150608/b0770d26/attachment.html>


More information about the CIG-SHORT mailing list