Eric Sonnenthal's Profile - TOUGH User Community

Skip to main content

Latest Posts

Re: FLOW VECTOR CANCELLING EACH OTHER FOR SOURCE/SINK TERMS
Hi Pierre, I'm not exactly sure what you did here, but will take a closer look. The area-weighting is the dvecta factor, calculated in t2f in the following loop: c.....…
- TOUGHREACT Developer
- Eric_Sonnenthal
- 1 yr ago
- TOUGHREACT
- Reported - view
Re: FLOW VECTOR CANCELLING EACH OTHER FOR SOURCE/SINK TERMS
Kenny- Writing the fluxes are a special feature in TOUGHREACT/TReactMech. Yes, there are some issues on the fluxes from GENER which are not known in the fluxes at connection when the fluxes/area are…
- TOUGHREACT Developer
- Eric_Sonnenthal
- 1 yr ago
- TOUGHREACT
- Reported - view
Re: modify output flowdata.tect
Hi Pierre, It should be pretty straightforward since the specific energy is the heat capacity of the phase x the fraction of the phase (solid, liquid, or gas) x density of the phase.…
- TOUGHREACT Developer
- Eric_Sonnenthal
- 1 yr ago
- TOUGHREACT
- Reported - view
Re: Confusion about CTOT and CGUESS
Hi Majid, Every concentration can be entered in different units depending on the ICON value. It is described in the manual. If ICON = 0, then it is pH or ppm for other species.…
- TOUGHREACT Developer
- Eric_Sonnenthal
- 1 yr ago
- TOUGHREACT
- Reported - view
Re: TOUGHREACT/ECO2N:CANNOT FIND PARAMETERS AT ELEMENT * 1*
Hi, The salt mass fraction is above halite saturation. Try starting with something much lower< 0.3.. Also start with a small amont of dissolved CO2, say 0.0001. Eric
- TOUGHREACT Developer
- Eric_Sonnenthal
- 1 yr ago
- TOUGHREACT
- Reported - view
Re: Mesh Creation
Hi Uzoma, We have a new online MESH and input file generator toughio-dash you can try out: https://toughio-dash.lbl.gov I have only tried it a few times, and person who wrote it has now left,…
- TOUGHREACT Developer
- Eric_Sonnenthal
- 1 yr ago
- 1
- Pre & Post Processors
- Reported - view
Re: How do I conduct a batch modeling of NaCl-water with primary minerals
Hi Efengwenbe, The last part of your question was "is there another way?" The way I solve these problems is to set up single grid block or 1-D problem with the observed mineralogy,…
- TOUGHREACT Developer
- Eric_Sonnenthal
- 1 yr ago
- 1
- TOUGHREACT
- Reported - view
Re: How do I conduct a batch modeling of NaCl-water with primary minerals
Hi Efengwengbe, Yes, you could put a small concentration (eg., 1e-10) for the initial and the guess concentration and then use icon=2, followed by the mineral (or gas) you want the solution to…
- TOUGHREACT Developer
- Eric_Sonnenthal
- 1 yr ago
- 1
- TOUGHREACT
- Reported - view
Re: How do I conduct a batch modeling of NaCl-water with primary minerals
Hi Efenwengbe, Pitzer sample problems 1 & 2 in the distribution package (Problem1-BlackSea-evap-eos1 & Problem2-Williston12-evap-eos1) in the Sample_Problems_Pitzer folder are single grid block…
- TOUGHREACT Developer
- Eric_Sonnenthal
- 1 yr ago
- TOUGHREACT
- Reported - view
Re: How to calculate what portion of a constituent adsorbs.
Hi Mikey, It looks like it is calculated in each time step, but I don't see where the cumulative amount is saved. I'll look into it. Eric
- TOUGHREACT Developer
- Eric_Sonnenthal
- 1 yr ago
- TOUGHREACT
- Reported - view
Re: SMCO2 is not increased
Hi Bobby, Allison, SmCO2 may not have worked in some of the older versions. You can always sum up the CO2 added to minerals and in solution from the tecplot output files.…
- TOUGHREACT Developer
- Eric_Sonnenthal
- 1 yr ago
- 1
- TOUGHREACT
- Reported - view
Re: How to monitor the temperature change of the sink with time?
Hi, I just saw this so maybe you figured it out already. The temperature of the fluid in the grid block you are taking fluid out of (with the MASS option) can be written out in FOFT or time.…
- TOUGHREACT Developer
- Eric_Sonnenthal
- 2 yrs ago
- TOUGHREACT
- Reported - view
Re: Multiple questions regarding data entry for a surface complexation problem
Hi Bobby, Regarding questions 3 onwards: 3. As for carbonate species, B(OH)3 is entered as a total concentration, so you just need to convert B concentration to B(OH3) concentration.…
- TOUGHREACT Developer
- Eric_Sonnenthal
- 2 yrs ago
- TOUGHREACT
- Reported - view
Re: Multiple questions regarding data entry for a surface complexation problem
Hi Babak, I'll try to answer some of these now, and get back to you on the rest. 1. When setting up the chemistry, the concentration you input is the total concentration,…
- TOUGHREACT Developer
- Eric_Sonnenthal
- 2 yrs ago
- TOUGHREACT
- Reported - view
Re: Restart time in inchem file is not the same as starting time in INCON - 1D Radial CO2 injection - Problem 5
Hi Majid, Your second question was about the "ramifications of running it this way". This general approach of running a steady-state, or "pseudo-steady-state" and then injection is the correct way…
- TOUGHREACT Developer
- Eric_Sonnenthal
- 2 yrs ago
- 1
- TOUGHREACT
- Reported - view
Re: Restart time in inchem file is not the same as starting time in INCON - 1D Radial CO2 injection - Problem 5
Hi Majid, The end time in flow.inp was only 1000 years, rather than 10,000 years. That fixed this problem. However, then the GENER block inside flow.inp also did not extend to the end time,…
- TOUGHREACT Developer
- Eric_Sonnenthal
- 2 yrs ago
- 1
- TOUGHREACT
- Reported - view
Re: Error in reading INCON - Problem 5, eco2n_1D Radial
Hi Majid, If you remove the last two lines, then the starting time is back to 0, but the conditions are whatever is in the INCON. TOUGHREACT looks at the starting simulation time and then looks in…
- TOUGHREACT Developer
- Eric_Sonnenthal
- 2 yrs ago
- 1
- TOUGHREACT
- Reported - view
Re: Error in reading INCON - Problem 5, eco2n_1D Radial
Hi Majid, At the bottom of SAVE, there are *** because the number of time steps and iterations exceeded the fortran format. To restart a file from time 0, delete these two lines,…
- TOUGHREACT Developer
- Eric_Sonnenthal
- 2 yrs ago
- 1
- TOUGHREACT
- Reported - view
Re: Error in reading INCON - Problem 5, eco2n_1D Radial
Hi Majid, Make sure the last time in GENER is 3000 years, and is equal to or greater than the final time in flow.inp. If you are still having troubles, can you send the modified input files,…
- TOUGHREACT Developer
- Eric_Sonnenthal
- 2 yrs ago
- TOUGHREACT
- Reported - view
Re: Permeability-Porosity relation
Hi Shijun, Thanks for pointing this out, I will look into and make the correction if necessary. Maybe this was supposed to be the half aperture, I don't remember.…
- TOUGHREACT Developer
- Eric_Sonnenthal
- 2 yrs ago
- TOUGHREACT
- Reported - view
Re: Error when use MESHMAKER to generate XYZ mesh in TOUGHREACT ECO2N sample problem
Hi Huanwei, Yes, a radial grid should be created with the RZ2D option. Aside from the TOUGHREACT Sample problems, there is a good discussion of the options for MESHMAKER in the TOUGH2 Manual.…
- TOUGHREACT Developer
- Eric_Sonnenthal
- 2 yrs ago
- TOUGHREACT
- Reported - view
Re: Problem 5 –1-D Radial Model for CO2 Sequestration in a Deep Saline Formation (ECO2N
Hi Alonso, Just change the end time in the PARAM block in flow.inp from 3.15569E08 to 3.15569E09 (seconds). The sample problems are usually set so they run in a short time just so the user can check…
- TOUGHREACT Developer
- Eric_Sonnenthal
- 2 yrs ago
- TOUGHREACT
- Reported - view
Re: Output flux of aqueous species
Hi Heather, There isn't an option to write out the chemical fluxes. However, if you multiply the fluid fluxes in flowvector.tec by the concentrations in the aqueous concentration tecplot file,…
- TOUGHREACT Developer
- Eric_Sonnenthal
- 2 yrs ago
- TOUGHREACT
- Reported - view
Re: TOUGHREACT_Debugging
Hi Shijun, Everything should be done in a stepwise approach, first flow, then injection, chemistry, perm change. etc., until you are sure the files and the results are correct. Look at the .…
Attachment
- TOUGHREACT Developer
- Eric_Sonnenthal
- 2 yrs ago
- TOUGHREACT
- Reported - view
Re: Tuning diffusion coefficients for water isotopes
Hi Tyler, I haven't had a chance to get back to this issue, but it is really important to figure out. I started putting in some code specific for water isotopes,…
- TOUGHREACT Developer
- Eric_Sonnenthal
- 2 yrs ago
- TOUGHREACT
- Reported - view

1 2 3 4 5 6

Back to top

Home

View all topics

TOUGH Website