39 #include "libmesh/quadrature.h"
40 #include "libmesh/elem.h"
41 #include "libmesh/unstructured_mesh.h"
42 #include "libmesh/fem_system.h"
54 this->
_dim = system->get_mesh().mesh_dimension();
59 libMesh::Gradient grad_u, grad_v;
60 grad_u = context.interior_gradient(this->
_flow_vars.
u(), qp);
61 grad_v = context.interior_gradient(this->
_flow_vars.
v(), qp);
63 libMesh::Real vorticity_value;
64 vorticity_value = fabs(grad_v(0) - grad_u(1));
68 libMesh::Gradient grad_w;
69 grad_w = context.interior_gradient(this->
_flow_vars.
w(), qp);
71 libMesh::Real vorticity_component_0 = grad_w(1) - grad_v(2);
72 libMesh::Real vorticity_component_1 = grad_u(2) - grad_v(0);
74 libMesh::Real term = vorticity_component_0*vorticity_component_0
75 + vorticity_component_1*vorticity_component_1
76 + vorticity_value*vorticity_value;
78 vorticity_value += std::sqrt(term);
81 return vorticity_value;
const VelocityFEVariables & _flow_vars
libMesh::Real vorticity(AssemblyContext &context, unsigned int qp) const
void init_variables(libMesh::FEMSystem *system)