Dear all,
I am trying to carry out the convergence analysis on the 2D axisymmetric
mesh. The problem is a simple Stokes flow. Before comparing the result
to the reference mesh, I compared the solution to itself.
// Compare the computed solution to the reference solution
ExactSolution ExSol(equation_systems);
ExSol.attach_reference_solution(&equation_systems);
string sys_name("Cell");
string var_name;
var_name = "vr";
ExSol.compute_error(sys_name, var_name);
Real h1_vr = ExSol.h1_error(sys_name, var_name);
var_name = "vz";
ExSol.compute_error(sys_name, var_name);
Real h1_vz = ExSol.h1_error(sys_name, var_name);
var_name = "p";
ExSol.compute_error(sys_name, var_name);
Real L2_p = ExSol.l2_error(sys_name, var_name);
The error should be zero, but the result is not very good...
H1 vr: 0.0189777
H1 vz: 1.74623e-09
L2 p: 2.2159e-08
Is this because internally, libmesh differentiate and integrate the
error over normal 2D domain? If so, which part should I modify?
Thanks,
Yusuke
------------------------------------------------------------------------------
Live Security Virtual Conference
Exclusive live event will cover all the ways today's security and
threat landscape has changed and how IT managers can respond. Discussions
will include endpoint security, mobile security and the latest in malware
threats. http://www.accelacomm.com/jaw/sfrnl04242012/114/50122263/
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
