calculation interrupted

All about WUFI 2D
Post Reply
avp081090
WUFI User
WUFI User
Posts: 38
Joined: Sun May 12, 2013 8:18 pm -1100

calculation interrupted

Post by avp081090 » Sun Mar 23, 2014 6:09 am -1100

My calculation was interupted at 96 % and it said that Maxit was to small..What does this mean? can i use the calculation even though it interupted it?
should i just increase maxit to solve this problem?

Best regards

veitner
WUFI SupportTeam IBP
WUFI SupportTeam IBP
Posts: 135
Joined: Tue Mar 08, 2005 11:14 pm -1100

Re: calculation interrupted

Post by veitner » Sun Mar 23, 2014 8:46 pm -1100

Please see explanation in following post

http://www.wufi-forum.com/viewtopic.php ... axit#p2100

Regards,

avp081090
WUFI User
WUFI User
Posts: 38
Joined: Sun May 12, 2013 8:18 pm -1100

Re: calculation interrupted

Post by avp081090 » Sun Apr 06, 2014 9:52 pm -1100

veitner wrote:Please see explanation in following post

http://www.wufi-forum.com/viewtopic.php ... axit#p2100

Regards,
Okay i have tried increasing my maxit and I have experienced the following:
- Calculation interrupted - maxit to small (ended at 97 %)
- Calculation interrupted - maxit to small (ended at 100 %)

What does this percentage mean? I understand that it is a measure of how far the calculations but fx when it interrupts at 97 % - is it then missing calculations for 3 % of the model geometry ?

Best regards

veitner
WUFI SupportTeam IBP
WUFI SupportTeam IBP
Posts: 135
Joined: Tue Mar 08, 2005 11:14 pm -1100

Re: calculation interrupted

Post by veitner » Sun Apr 06, 2014 10:12 pm -1100

Calculation interrupted - maxit to small (ended at 97 %)
I'm sorry but I do know knothing about this percentage. Where did you find that?
I only know
maxit_msg.png
maxit_msg.png (622 Bytes) Viewed 10013 times
.

avp081090
WUFI User
WUFI User
Posts: 38
Joined: Sun May 12, 2013 8:18 pm -1100

Re: calculation interrupted

Post by avp081090 » Sun Apr 06, 2014 10:19 pm -1100

veitner wrote:
Calculation interrupted - maxit to small (ended at 97 %)
I'm sorry but I do know knothing about this percentage. Where did you find that?
I only know
maxit_msg.png
.
Oh i can see my post could be misunderstood - its not that it displays "(ended at 97 %)" im talking about the large blue line at the bottom of the screen. Sorry for the misconception.

Best regards

veitner
WUFI SupportTeam IBP
WUFI SupportTeam IBP
Posts: 135
Joined: Tue Mar 08, 2005 11:14 pm -1100

Re: calculation interrupted

Post by veitner » Sun Apr 06, 2014 11:28 pm -1100

I see - the blue bar is an overall progress indicator e.g. 97% means current_time_step/total_time_steps=0.97
And *** interrupted - maxit too small *** means the convergence criterion was not reached within maxit iterations - which does not necessarily mean you got a "wrong solution" but there is a chance to get a slightly better solution for that specific timestep by increasing maxit. Sometimes one can get some kind of oscillation: current temperature field results in temperature dependent moisture transport coefficients leading to a moisture field yielding to different moisture dependent heat transport coefficients resulting in different temperature field and so on - without finding a "converged solution". WUFI assumes a solution of both temperatur and r.h. as converged as soon as two subsequent solutions of temperature and r.h. do not differ more than "convergence crit." (maximum norm, residual). The independent solutions of both temperature and r.h. are exact for each single inner iteration - but the coupled solution has not underrun the user definable convergence criterion perhaps because of huge nonlinearity.

Following chart (hypothetical values) might illustrate the iterative process and meaning of the convergence criterion.
konvergenz.png
konvergenz.png (12.32 KiB) Viewed 10010 times
iterative Naeherungsloesung: iterative approximative solution
unbekannte Loesung: unknown solution (the exact or analytical solution)
Unterschied zum vorhergehenden Loesungsschritt: difference to last solution (residual)
Konvergenzkriterium: convergence criterion
Abweichung: difference between "converged solution" and "exact solution"

The left axis is for the red and the black line
The right axis for the blue and dark green line

The convergence criterion has been reached at about iteration 800 (light green line) and at this point the difference between the exact and iterative solution is about 0.0056 degree celsius.

But remember: the exact solution is unknown - therefor WUFI :lol:

veitner
WUFI SupportTeam IBP
WUFI SupportTeam IBP
Posts: 135
Joined: Tue Mar 08, 2005 11:14 pm -1100

Re: calculation interrupted

Post by veitner » Mon Apr 07, 2014 12:01 am -1100

Here comes another chart illustrating the residuals comparing two different approaches of source term linearization.
iter_small.png
iter_small.png (4.18 KiB) Viewed 10010 times
This shows the real course of the residuals during the iterative process for a single time step.

You can create a similiar chart by setting "Print residuals at iteration modulo" (computational settings/enhanced) to 1. Then redo the calculation and open the "result analyzer" (shows up without logarithmic axis scaling). Do not forget to reset the "print residuals modulo" before saving your data because a value of 1 leads to a really huge "iter.txt" file.

avp081090
WUFI User
WUFI User
Posts: 38
Joined: Sun May 12, 2013 8:18 pm -1100

Re: calculation interrupted

Post by avp081090 » Mon Apr 07, 2014 12:35 am -1100

Tanks, makes a lot mre sense now :)

Best Regards

avp081090
WUFI User
WUFI User
Posts: 38
Joined: Sun May 12, 2013 8:18 pm -1100

Re: calculation interrupted

Post by avp081090 » Mon Apr 07, 2014 1:48 am -1100

But just to clarify one last thing: If it interrupts the calculations at 97 % saying that maxit is to small. Does it mean that 97 % of the time interval have converged?

Best Regards

veitner
WUFI SupportTeam IBP
WUFI SupportTeam IBP
Posts: 135
Joined: Tue Mar 08, 2005 11:14 pm -1100

Re: calculation interrupted

Post by veitner » Mon Apr 07, 2014 1:56 am -1100

Hmhh - it does not interrupt the calculation at 97%

It does not reach the convergence criterion at timestep n/8760=0.97->n=0.97*8760=8498 (assuming you intend to simulate 1 year=8760 timesteps of 1 h each).
After that it continues with timestep 8499 using the "unconverged" solution as a starting point. But unconverged does not mean bad - it just could be better and has not reached the user configurable convergence criterion.

Post Reply