Commit 069d2a7c authored by Elia Merzari's avatar Elia Merzari

Reverting a commit of April Novak that is not working properly with this...

Reverting a commit of April Novak that is not working properly with this version of Giraffe, to be investiagetd further
parent 6e5ca674
Pipeline #2753 failed with stages
in 314 minutes and 26 seconds
...@@ -76,10 +76,14 @@ NekSideIntegralVariableUserObject::computeQpIntegral() ...@@ -76,10 +76,14 @@ NekSideIntegralVariableUserObject::computeQpIntegral()
// There is an added correction factor that accounts for the surface area of the pin // 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 // The analytic expression is (2 / R / delta z) but we compute numerically as
// (2 / (Surf Area / (2*PI) ) ) = 4PI / Surf Area // (2 / (Surf Area / (2*PI) ) ) = 4PI / Surf Area
RealVectorValue coupled_grad_vec(_coupled_grad[_qp](0), _coupled_grad[_qp](1),
_coupled_grad[_qp](2)); // Implementation from April Novak, to be tested more carefully
RealVectorValue normals_vec(_normals[_qp](0), _normals[_qp](1), _normals[_qp](2)); // RealVectorValue coupled_grad_vec(_coupled_grad[_qp](0), _coupled_grad[_qp](1), _coupled_grad[_qp](2));
return -_diffusion_coefficient[_qp] * coupled_grad_vec * normals_vec * l_poly_val * f_poly_val * 4.0*M_PI / _surface_area_pp; // RealVectorValue normals_vec(_normals[_qp](0), _normals[_qp](1), _normals[_qp](2));
// return -_diffusion_coefficient[_qp] * coupled_grad_vec * normals_vec * l_poly_val * f_poly_val * 4.0*M_PI / _surface_area_pp;
// Original implementation from Matt Ellis
return -_diffusion_coefficient[_qp] * _coupled_grad[_qp].norm() * l_poly_val * f_poly_val * 4.0*M_PI / _surface_area_pp;
} }
void void
NekSideIntegralVariableUserObject::finalize() 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