> I don't know why, when using a small time step, I get oscillations. The smaller the time step the larger the oscillations. I would have thought, due to the Courant conditions, the results would be more precise with a smaller time steps. I think it is because you are using centered in space. I'm not expert on error analysis of numerical methods for finite difference, but I think centered in space can lead to oscillations. I almost always use -space 0; -time 1 to avoid this possibility. If you want to know more, I can ask Ken Kipp to give us a thorough explanation. David David Parkhurst (dlpark@xxxxxxxx) U.S. Geological Survey Box 25046, MS 413 Denver Federal Center Denver, CO 80225 Project web page: https://wwwbrr.cr.usgs.gov/projects/GWC_coupled
Please note that some U.S. Geological Survey (USGS) information accessed through this page may be preliminary in nature and presented prior to final review and approval by the Director of the USGS. This information is provided with the understanding that it is not guaranteed to be correct or complete and conclusions drawn from such information are the sole responsibility of the user.
Any use of trade, product, or firm names in this publication is for descriptive purposes only and does not imply endorsement by the U.S. Government.
The URL of this page is:
https://wwwbrr.cr.usgs.gov/projects/GWC_coupled/phreeqc/mail/msg00487.html
Email:dlpark@usgs.gov
Last modified: $Date: 2005-09-13 21:04:21 -0600 (Tue, 13 Sep 2005) $
Visitor number [an error occurred while processing this directive] since Jan 22, 1998.