CShort Posted November 8, 2016 Share Posted November 8, 2016 I'm working with a relatively simple X1T model in which water 1 is injected into a sandstone containing water 2 for 22 years then chased through the system with water 1. During transport, the only attenuation process operating is an activity-KD retardation model. I'm trying to run this model using different porosity's, assuming that the rest of the volume is just non-reactive quartz. This should be a simple model. This model runs for certain values of concentrations and Kd's, but not others, although the model is fundamentally unchanged. The biggest problem is that the model is unable to initialize because the "Newton-Raphson did not converge...". I've been able to duplicate this issue in React. One issue I've noticed causes this problem is declaring a porosity. The modelled sandstone has a porosity of 0.25. When I run my react code using a blank porosity, blank bulk volume, and 0 inert volume (the default values), it's able to initialize. If I declare a porosity of 0.25, it is unable to initialize. Do you have any thoughts on what might be going on here or what I should look out for when trying to declare porosity values? I've referenced both the reaction guide and the reference guide and have not had much luck sorting this out. Thank you! Clark React_Issue.rea Kd_Ambrosia_Por=25_KD=0.029.sdat Quote Link to comment Share on other sites More sharing options...
Join the conversation
You can post now and register later. If you have an account, sign in now to post with your account.