Dear getfem users,I am using getfem for solving an axisymmetric problem and I have to integrate a term of the following form :
\int_\Omega 1/r u v dr dzwhere u and v are respectively solution and test function, r is the radial distance and z is the height. Of course u and v are enforced to be zero on the axis. Consequently the integrand should be regular even for the elements touching the axis. I wonder how to perform an accurate quadrature of this term particularly for the elements touching the axis because I do not know how to describe the function 1/r properly by a mesh_fem object (only in the area close to the axis). I could do something by hand for the elements touching the axis but it would be helpful if someone had a better proposition.
Thanks in advance for your proposition(s). Best regards, Ronan
<<attachment: ronan_perrussel.vcf>>
_______________________________________________ Getfem-users mailing list [email protected] https://mail.gna.org/listinfo/getfem-users
