#include <vorticity.h>

Inheritance diagram for GRINS::Vorticity:

Collaboration diagram for GRINS::Vorticity:

Public Member Functions
	Vorticity (const std::string &qoi_name)
	Constructor. More...

virtual	~Vorticity ()

virtual QoIBase *	clone () const
	Required to provide clone (deep-copy) for adding QoI object to libMesh objects. More...

virtual bool	assemble_on_interior () const
	Does the QoI need an element interior assembly loop? More...

virtual bool	assemble_on_sides () const
	Does the QoI need a domain boundary assembly loop? More...

virtual void	init (const GetPot &input, const MultiphysicsSystem &system, unsigned int qoi_num)
	Initialize local variables. More...

virtual void	init_context (AssemblyContext &context)

virtual void	element_qoi (AssemblyContext &context, const unsigned int qoi_index)
	Compute the qoi value. More...

virtual void	element_qoi_derivative (AssemblyContext &context, const unsigned int qoi_index)
	Compute the qoi derivative with respect to the solution. More...

Public Member Functions inherited from GRINS::QoIBase
	QoIBase (const std::string &qoi_name)

virtual	~QoIBase ()

virtual void	reinit (MultiphysicsSystem &)
	Reinitialize QoI. More...

virtual void	side_qoi (AssemblyContext &, const unsigned int)
	Compute the qoi value on the domain boundary. More...

virtual void	side_qoi_derivative (AssemblyContext &, const unsigned int)
	Compute the qoi derivative with respect to the solution on the domain boundary. More...

virtual void	parallel_op (const libMesh::Parallel::Communicator &communicator, libMesh::Number &sys_qoi, libMesh::Number &local_qoi)
	Call the parallel operation for this QoI and cache the value. More...

virtual void	thread_join (libMesh::Number &qoi, const libMesh::Number &other_qoi)
	Call the operation to accumulate this QoI from multiple threads. More...

virtual void	output_qoi (std::ostream &out) const
	Basic output for computed QoI's. More...

libMesh::Number	value () const
	Returns the current QoI value. More...

const std::string &	name () const
	Returns the name of this QoI. More...

Public Member Functions inherited from GRINS::ParameterUser
	ParameterUser (const std::string &user_name)

virtual	~ParameterUser ()

virtual void	set_parameter (libMesh::Number &param_variable, const GetPot &input, const std::string &param_name, libMesh::Number param_default)
	Each subclass can simultaneously read a parameter value from. More...

virtual void	set_parameter (libMesh::ParsedFunction< libMesh::Number, libMesh::Gradient > &func, const GetPot &input, const std::string &func_param_name, const std::string &param_default)
	Each subclass can simultaneously read a parsed function from. More...

virtual void	set_parameter (libMesh::ParsedFEMFunction< libMesh::Number > &func, const GetPot &input, const std::string &func_param_name, const std::string &param_default)
	Each subclass can simultaneously read a parsed function from. More...

virtual void	move_parameter (const libMesh::Number &old_parameter, libMesh::Number &new_parameter)
	When cloning an object, we need to update parameter pointers. More...

virtual void	move_parameter (const libMesh::ParsedFunction< libMesh::Number, libMesh::Gradient > &old_func, libMesh::ParsedFunction< libMesh::Number, libMesh::Gradient > &new_func)
	When cloning an object, we need to update parameter pointers. More...

virtual void	move_parameter (const libMesh::ParsedFEMFunction< libMesh::Number > &old_func, libMesh::ParsedFEMFunction< libMesh::Number > &new_func)
	When cloning an object, we need to update parameter pointers. More...

virtual void	register_parameter (const std::string &param_name, libMesh::ParameterMultiAccessor< libMesh::Number > &param_pointer) const
	Each subclass will register its copy of an independent. More...

Protected Attributes
const VelocityVariable *	_flow_vars

std::set< libMesh::subdomain_id_type >	_subdomain_ids
	List of sumdomain ids for which we want to compute this QoI. More...

Protected Attributes inherited from GRINS::QoIBase
std::string	_qoi_name

libMesh::Number	_qoi_value

Private Member Functions
	Vorticity ()
	User never call default constructor. More...

Additional Inherited Members
Static Public Attributes inherited from GRINS::ParameterUser
static std::string	zero_vector_function = std::string("{0}")
	A parseable function string with LIBMESH_DIM components, all 0. More...

Detailed Description

Vorticity QoI.

This class implement a vorticity QoI that can be used to both compute QoI values and drive QoI-based adaptive refinement. Currently, this QoI is only implemented in 2D and will error if it detects a three-dimensional problem.

Definition at line 44 of file vorticity.h.

Constructor & Destructor Documentation

GRINS::Vorticity::Vorticity ( const std::string & qoi_name )

Constructor.

Constructor takes GetPot object to read any input options associated with this QoI

Definition at line 45 of file vorticity.C.

     : QoIBase(qoi_name)
   {
     return;
   }

GRINS::Vorticity::~Vorticity ( )

virtual

Definition at line 51 of file vorticity.C.

   {
     return;
   }

GRINS::Vorticity::Vorticity ( )

private

User never call default constructor.

Referenced by clone().

Member Function Documentation

bool GRINS::Vorticity::assemble_on_interior ( ) const

inlinevirtual

Does the QoI need an element interior assembly loop?

This is pure virtual to force to user to specify.

Implements GRINS::QoIBase.

Definition at line 98 of file vorticity.h.

   {
     return true;
   }

bool GRINS::Vorticity::assemble_on_sides ( ) const

inlinevirtual

Does the QoI need a domain boundary assembly loop?

This is pure virtual to force to user to specify.

Implements GRINS::QoIBase.

Definition at line 104 of file vorticity.h.

   {
     return false;
   }

QoIBase * GRINS::Vorticity::clone ( ) const

virtual

Required to provide clone (deep-copy) for adding QoI object to libMesh objects.

Implements GRINS::QoIBase.

Definition at line 56 of file vorticity.C.

References Vorticity().

   {
     return new Vorticity( *this );
   }

void GRINS::Vorticity::element_qoi	(	AssemblyContext &	context,
		const unsigned int	qoi_index
	)

virtual

Compute the qoi value.

Currently, only implemented for 2D. Assumes that the vorticity will be computed over area of input subdomain id. Vorticity computed as $\int_{\Omega} \nabla \times \mathbf{u} \; d\mathbf{x}$

Todo:: Need to generalize this to the multiple QoI case

Reimplemented from GRINS::QoIBase.

Definition at line 100 of file vorticity.C.

References _flow_vars, _subdomain_ids, GRINS::MultcomponentVectorVariable::u(), and GRINS::MultcomponentVectorVariable::v().

   {
 
     if( _subdomain_ids.find( (&context.get_elem())->subdomain_id() ) != _subdomain_ids.end() )
       {
         libMesh::FEBase* element_fe;
         context.get_element_fe<libMesh::Real>(_flow_vars->u(), element_fe);
         const std::vector<libMesh::Real> &JxW = element_fe->get_JxW();
 
         unsigned int n_qpoints = context.get_element_qrule().n_points();
 
         libMesh::Number& qoi = context.get_qois()[qoi_index];
 
         for( unsigned int qp = 0; qp != n_qpoints; qp++ )
           {
             libMesh::Gradient grad_u = 0.;
             libMesh::Gradient grad_v = 0.;
             context.interior_gradient( _flow_vars->u(), qp, grad_u );
             context.interior_gradient( _flow_vars->v(), qp, grad_v );
             qoi += (grad_v(0) - grad_u(1)) * JxW[qp];
           }
       }
 
     return;
   }

void GRINS::Vorticity::element_qoi_derivative	(	AssemblyContext &	context,
		const unsigned int	qoi_index
	)

virtual

Compute the qoi derivative with respect to the solution.

Currently, only implemented for 2D. Assumes that the vorticity will be computed over area of input subdomain id.

Todo:: Need to generalize this to the multiple QoI case

Reimplemented from GRINS::QoIBase.

Definition at line 128 of file vorticity.C.

References _flow_vars, _subdomain_ids, GRINS::MultcomponentVectorVariable::u(), and GRINS::MultcomponentVectorVariable::v().

   {
 
     if( _subdomain_ids.find( (&context.get_elem())->subdomain_id() ) != _subdomain_ids.end() )
       {
         // Element
         libMesh::FEBase* element_fe;
         context.get_element_fe<libMesh::Real>(_flow_vars->u(), element_fe);
 
         // Jacobian times integration weights
         const std::vector<libMesh::Real> &JxW = element_fe->get_JxW();
 
         // Grad of basis functions
         const std::vector<std::vector<libMesh::RealGradient> >& du_phi =
           context.get_element_fe(_flow_vars->u())->get_dphi();
         const std::vector<std::vector<libMesh::RealGradient> >& dv_phi =
           context.get_element_fe(_flow_vars->v())->get_dphi();
 
         // Local DOF count and quadrature point count
         const unsigned int n_T_dofs = context.get_dof_indices(0).size();
         unsigned int n_qpoints = context.get_element_qrule().n_points();
 
         // Warning: we assume here that vorticity is the only QoI!
         // This should be consistent with the assertion in grins_mesh_adaptive_solver.C
         libMesh::DenseSubVector<libMesh::Number> &Qu = context.get_qoi_derivatives(qoi_index, _flow_vars->u());
         libMesh::DenseSubVector<libMesh::Number> &Qv = context.get_qoi_derivatives(qoi_index, _flow_vars->v());
 
         // Integration loop
         for( unsigned int qp = 0; qp != n_qpoints; qp++ )
           {
             for( unsigned int i = 0; i != n_T_dofs; i++ )
               {
                 Qu(i) += - dv_phi[i][qp](1) * JxW[qp];
                 Qv(i) += du_phi[i][qp](0) * JxW[qp];
               }
           }
       }
 
     return;
   }

void GRINS::Vorticity::init	(	const GetPot &	input,
		const MultiphysicsSystem &	system,
		unsigned int	qoi_num
	)

virtual

Initialize local variables.

Any local variables that need information from libMesh get initialized here. For example, variable indices.

Reimplemented from GRINS::QoIBase.

Definition at line 62 of file vorticity.C.

References GRINS::VariablesParsing::QOI, and GRINS::VariablesParsing::velocity_variable_name().

   {
     // Extract subdomain on which to compute to qoi
     int num_ids = input.vector_variable_size( "QoI/Vorticity/enabled_subdomains" );
 
     if( num_ids == 0 )
       {
         std::cerr << "Error: Must specify at least one subdomain id on which to compute vorticity." << std::endl;
         libmesh_error();
       }
 
     for( int i = 0; i < num_ids; i++ )
       {
         libMesh::subdomain_id_type s_id = input( "QoI/Vorticity/enabled_subdomains", -1, i );
         _subdomain_ids.insert( s_id );
       }
 
     _flow_vars = &(GRINSPrivate::VariableWarehouse::get_variable_subclass<VelocityVariable>(VariablesParsing::velocity_variable_name(input,"Vorticity",VariablesParsing::QOI)));
   }

void GRINS::Vorticity::init_context ( AssemblyContext & context )

virtual

Reimplemented from GRINS::QoIBase.

Definition at line 84 of file vorticity.C.

References _flow_vars, GRINS::MultcomponentVectorVariable::u(), and GRINS::MultcomponentVectorVariable::v().

   {
     libMesh::FEBase* u_fe = NULL;
     libMesh::FEBase* v_fe = NULL;
 
     context.get_element_fe<libMesh::Real>(_flow_vars->u(), u_fe);
     context.get_element_fe<libMesh::Real>(_flow_vars->v(), v_fe);
 
     u_fe->get_dphi();
     u_fe->get_JxW();
 
     v_fe->get_dphi();
 
     return;
   }

Member Data Documentation

const VelocityVariable* GRINS::Vorticity::_flow_vars

protected

Definition at line 86 of file vorticity.h.

Referenced by element_qoi(), element_qoi_derivative(), and init_context().

std::set<libMesh::subdomain_id_type> GRINS::Vorticity::_subdomain_ids

protected

List of sumdomain ids for which we want to compute this QoI.

Definition at line 89 of file vorticity.h.

Referenced by element_qoi(), and element_qoi_derivative().

The documentation for this class was generated from the following files:

src/qoi/include/grins/vorticity.h
src/qoi/src/vorticity.C

Public Member Functions

Protected Attributes

Private Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation