val46 
April 4, 2012 05:53 
plot residuals in multiregion solver
Hello,
I would like to plot the residuals from a logfile with gnuplot.
A nice tutorial for that can be found here.
With single region solvers it works nice but not with multiregion solvers.
As an example, here is a timestep of a tutorial case:
Code:
Time = 0.0469758
Solving for fluid region bottomAir
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Ux, Initial residual = 0.05933919, Final residual = 1.875811e09, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.07216778, Final residual = 2.727341e09, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 0.06497582, Final residual = 2.468975e09, No Iterations 2
DILUPBiCG: Solving for h, Initial residual = 0.1391416, Final residual = 5.374545e08, No Iterations 2
Min/max T:300 319.7154
GAMG: Solving for p_rgh, Initial residual = 0.1779339, Final residual = 0.0008159564, No Iterations 3
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors (bottomAir): sum local = 5.094676e07, global = 1.166084e10, cumulative = 1.841265e07
GAMG: Solving for p_rgh, Initial residual = 0.008835792, Final residual = 2.173498e07, No Iterations 1000
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors (bottomAir): sum local = 1.498523e10, global = 2.259186e12, cumulative = 1.841288e07
Solving for fluid region topAir
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Ux, Initial residual = 0.08491977, Final residual = 4.188711e09, No Iterations 4
DILUPBiCG: Solving for Uy, Initial residual = 0.1317485, Final residual = 9.808667e08, No Iterations 4
DILUPBiCG: Solving for Uz, Initial residual = 0.094053, Final residual = 3.943191e08, No Iterations 4
DILUPBiCG: Solving for h, Initial residual = 0.00873557, Final residual = 9.642108e08, No Iterations 3
Min/max T:299.9996 300
GAMG: Solving for p_rgh, Initial residual = 0.03020841, Final residual = 0.0001659505, No Iterations 2
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors (topAir): sum local = 3.574824e07, global = 2.221026e07, cumulative = 0.0002244136
GAMG: Solving for p_rgh, Initial residual = 0.003304726, Final residual = 8.552481e08, No Iterations 13
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors (topAir): sum local = 2.001775e10, global = 5.684051e11, cumulative = 0.0002244136
Solving for solid region heater
DICPCG: Solving for T, Initial residual = 0.05991691, Final residual = 1.506036e08, No Iterations 2
DICPCG: Solving for T, Initial residual = 2.892845e06, Final residual = 8.470107e10, No Iterations 1
Min/max T:min(T) [0 0 0 1 0 0 0] 300 max(T) [0 0 0 1 0 0 0] 500
Solving for solid region leftSolid
DICPCG: Solving for T, Initial residual = 0.9645239, Final residual = 1.318008e07, No Iterations 2
DICPCG: Solving for T, Initial residual = 1.649023e05, Final residual = 5.392089e09, No Iterations 1
Min/max T:min(T) [0 0 0 1 0 0 0] 300 max(T) [0 0 0 1 0 0 0] 300
Solving for solid region rightSolid
DICPCG: Solving for T, Initial residual = 0.6954516, Final residual = 1.022314e07, No Iterations 2
DICPCG: Solving for T, Initial residual = 1.153665e05, Final residual = 3.904825e09, No Iterations 1
Min/max T:min(T) [0 0 0 1 0 0 0] 300 max(T) [0 0 0 1 0 0 0] 300
ExecutionTime = 3.04 s ClockTime = 3 s
I'm struckling to plot for example Tmax of the heater region.
Any hint is appreciated.
Best Regards,
Toni
