Numerical convergence problem
Numerical convergence problem
Dear Sir,
Sometimes HYDRUS does not stop on its own after showing the error i.e. numerical solution not converged! even if I have not checked the option Hit enter at the end. Please help me regarding this issue.
Sometimes HYDRUS does not stop on its own after showing the error i.e. numerical solution not converged! even if I have not checked the option Hit enter at the end. Please help me regarding this issue.
Re: Numerical convergence problem
Since the Richards equation is a highly nonlinear equation, there is always a chance that the numerical solution will not converge. However, for physically realistic problems, this should be rather rate. J.
Re: Numerical convergence problem
Dear Sir,
My issue is not with the nonconvergence of the simulation, the problem is sometimes it does not stops on its own after showing the error, I need to forced stop it. I am attaching the image which shows the problem I am facing.
Thank you
K Rajsekhar
My issue is not with the nonconvergence of the simulation, the problem is sometimes it does not stops on its own after showing the error, I need to forced stop it. I am attaching the image which shows the problem I am facing.
Thank you
K Rajsekhar
 Attachments

 HYDRUS image.JPG (67.38 KiB) Viewed 1344 times
Re: Numerical convergence problem
I would need to see your project. However, I would not have time to look at it until the end of September as I’m currently traveling in Himachal Pradesh, India. J.
Re: Numerical convergence problem
Dear Sir,
I have attached my project, please check it. And have a happy stay in INDIA. If you have time do please visit IISc Bangalore. I am doing PhD under prof. Sekhar Muddu in IISc. I guess you visited IISc long back for a short course on HYDRUS.
I have attached my project, please check it. And have a happy stay in INDIA. If you have time do please visit IISc Bangalore. I am doing PhD under prof. Sekhar Muddu in IISc. I guess you visited IISc long back for a short course on HYDRUS.
 Attachments

 daily.rar
 (513 Bytes) Downloaded 16 times
Re: Numerical convergence problem
Dear Sir,
Please check this attachment, in the earlier attachment forgot to add some files.
Please check this attachment, in the earlier attachment forgot to add some files.
 Attachments

 daily.rar
 (16.45 KiB) Downloaded 13 times
Re: Numerical convergence problem
Dear Sir,
what I realized is when the head is less than 10000000 the above problem arises i.e. hydrus doesnot stop after showing the error. Don't know how to tackle this issue.
what I realized is when the head is less than 10000000 the above problem arises i.e. hydrus doesnot stop after showing the error. Don't know how to tackle this issue.
Re: Numerical convergence problem
Dear Sir,
Issue is resolved. Actually the problem was with hCritA, as I fixed the hCritA to 15000 and for n=1.05 h values are around 10e+20.
Issue is resolved. Actually the problem was with hCritA, as I fixed the hCritA to 15000 and for n=1.05 h values are around 10e+20.
Re: Numerical convergence problem
Dear Sir,
The problem still persists. The issue is, if the solution is not converged it should stop, but in my case it is going into infinite loop with all NaN values. Please help me fixing the issue. In the attachments, I have enclosed the project with name daily and also my executable. As an example I have taken n=1.05, which is creating the issue, but in my Monte Carlo simulations different combination of parameters create the same issue. So, what I want is when the solution is not converged, it should stop and not go to infinite loop.
Thank you
Rajsekhar
The problem still persists. The issue is, if the solution is not converged it should stop, but in my case it is going into infinite loop with all NaN values. Please help me fixing the issue. In the attachments, I have enclosed the project with name daily and also my executable. As an example I have taken n=1.05, which is creating the issue, but in my Monte Carlo simulations different combination of parameters create the same issue. So, what I want is when the solution is not converged, it should stop and not go to infinite loop.
Thank you
Rajsekhar
 Attachments

 daily.rar
 (528.26 KiB) Downloaded 13 times
Re: Numerical convergence problem
Dear Sir,
I have attached the screen shot of the simulation runs of two settings. One with max no of iterations=10 and other with max no of iterations=3. when I choose max no of iterations=10, the simulation is not stopping after convergence issue, but with the option max no of iterations=3, after convergence issue the simulation stops. I want my simulation to stop after the convergence issue with max no of iterations=10.
Thank you.
I have attached the screen shot of the simulation runs of two settings. One with max no of iterations=10 and other with max no of iterations=3. when I choose max no of iterations=10, the simulation is not stopping after convergence issue, but with the option max no of iterations=3, after convergence issue the simulation stops. I want my simulation to stop after the convergence issue with max no of iterations=10.
Thank you.
 Attachments

 max_no_of_iterations=10.JPG (94.66 KiB) Viewed 587 times

 max_no_of_iterations=3.JPG (89.47 KiB) Viewed 587 times
Re: Numerical convergence problem
Dear Sir,
It is happening for energy balance boundary condition only. For PenmannMonteith equation it is working fine, that is after nonconvergence it is automatically stopping. But in my work I am using energy balance BC, so want the simulation to stop after nonconvergence.
It is happening for energy balance boundary condition only. For PenmannMonteith equation it is working fine, that is after nonconvergence it is automatically stopping. But in my work I am using energy balance BC, so want the simulation to stop after nonconvergence.
Re: Numerical convergence problem
You like use some unrealistic combination of parameters, which somewhere apparently leads to a division by zero, which leads to NaN (standing for not a number, is a member of a numeric data type that can be interpreted as a value that is undefined or unrepresentable, especially in floatingpoint arithmetic). I do not know how to help with that. J.
Re: Numerical convergence problem
Dear Sir,
Sorry for bothering you. But, is there a way to stop the simulation automatically when there is a NaN value. I mean in the code itself, if there is a NaN value, the code will detect the NaN and stops automatically.
Sorry for bothering you. But, is there a way to stop the simulation automatically when there is a NaN value. I mean in the code itself, if there is a NaN value, the code will detect the NaN and stops automatically.
Re: Numerical convergence problem
I wish I knew how to do it in the Intel compiler. J.