Commit cb870daf authored by Elia Merzari's avatar Elia Merzari
Browse files

Fixing bug in surface correction

parent d8184540
......@@ -74,7 +74,7 @@ NekSideIntegralVariableUserObject::computeQpIntegral()
// There is an added correction factor that accounts for the surface area of the pin
// The analytic expression is (2 / R / delta z) but we compute numerically as
// (2 / (Surf Area / (2*PI) ) ) = 4PI / Surf Area
return -_diffusion_coefficient[_qp] * _coupled_grad[_qp].norm() * l_poly_val * f_poly_val * 4.0 / M_PI / _surface_area_pp;
return -_diffusion_coefficient[_qp] * _coupled_grad[_qp].norm() * l_poly_val * f_poly_val * 4.0*M_PI / _surface_area_pp;
}
void
NekSideIntegralVariableUserObject::finalize()
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment