36 #include "libmesh/quadrature.h"
37 #include "libmesh/elem.h"
38 #include "libmesh/unstructured_mesh.h"
39 #include "libmesh/fem_system.h"
50 this->
_dim = system->get_mesh().mesh_dimension();
59 libMesh::Gradient grad_u, grad_v;
63 libMesh::Real vorticity_value;
64 vorticity_value = fabs(grad_v(0) - grad_u(1));
68 libMesh::Gradient grad_w;
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;
VariableIndex w_var() const
PrimitiveFlowVariables _flow_vars
virtual void init(libMesh::FEMSystem *system)
VariableIndex u_var() const
libMesh::Real vorticity(AssemblyContext &context, unsigned int qp) const
VariableIndex v_var() const
void init_variables(libMesh::FEMSystem *system)