[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: Phreeqci



M Parkhurst,

I'm not running a transport problem.
In order to avoit convergence problems I use very small time stpes (-steps 86400 in 1500 steps for example). This leads to a quit big outpout file because for each time step many informatin is printed by default.
In order to avoid this huge amount of information, avoid convergence problems and go farther on time, I am trying to use the INCREMENTAL_REACTIONS true keyword after the KINETIKS 1 block but it does'nt work at all. The first time step is ok but it never converges at the Reaction step 2 (even for very small time steps). 
Thanks again.
Salvador


At 18:35 27/11/01 -0700, you wrote:

>I am unsing Phreeqci 2.5.0.2021 under Windows98, PIII/1000MH, 128 M RAM,
>and I am having a couple of problems:
>> 1) After running a case, when the last time step is reached, the
>calculations stops and then it tries automatically to open the newly
>generated output file, but it last a lot, a very long period of time
>(several minutes trying to open the file) and often it hungs the program
>(not responding...).
>
>Are you running a transport problem? It could be that the output file is
>very large and your system is having trouble with it. You can limit the
>data that is written to the output file with -print_cells and
>-print_frequency.
>
>> 2) In the input file I have defined a Kinetics bloc and a Rates bloc in a
>batch simulation. A very single 1st order kinetic law is defined but often
>the program hungs because it does not converge. It is very strange because
>the final concentration is only the half of the initial, the time steps are
>very small (1 second),....I am not able to get than the simulation go
>farther (seversl days) and the final concentration several orders of
>magnitude smaller than the initial.
>
>The time step for kinetics is defined either in the KINETICS data block or
>the TRANSPORT (or ADVECTION) data block, depending on the type of
>simulation you are running. The value is defined in seconds. The default
>value is 1 second. Sounds like you may be using the default value and
>should define your own time step.
>
>If these suggestions don't help, please send the input file (and the
>database file if it is different from the database files distributed with
>PhreeqcI). I'll take a look at it.
>
>David
>
>David Parkhurst (dlpark@xxxxxxxx)
>U.S. Geological Survey
>Box 25046, MS 413
>Denver Federal Center
>Denver, CO 80225 

----------------------------------
Salvador JORDANA
ENVIROSQUANTISCI
Av. Univ. Autònoma,3
Parc Tecnològic del Vallès
08290 CERDANYOLA DEL VALLÈS
BARCELONA-SPAIN
T+34 935 824 410
F+34 935 824 412
mailto:sjordana@xxxxxxxxxxxx
http://www.quantisci.es
----------------------------------



Project Home Page
Complete Water Resources Division Software
USGS Home Page
Water Resources Division Home Page
NRP Home Page
Help Page
USGS Privacy Statement       

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/msg00356.html
Email:dlpark@usgs.gov
Last modified: $Date: 2005-09-13 21:04:21 -0600 (Tue, 13 Sep 2005) $
Visitor number 1663 since Jan 22, 1998.