ADGateValve1PhaseUserObject.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html

#pragma once

#include "ADFlowJunctionUserObject.h"

class ADNumericalFlux3EqnBase;

class ADGateValve1PhaseUserObject : public ADFlowJunctionUserObject
{
public:
  ADGateValve1PhaseUserObject(const InputParameters & params);

  virtual void initialize() override;
  virtual void execute() override;
  virtual void threadJoin(const UserObject & uo) override;
  virtual void finalize() override;

  const std::vector<ADReal> & getFlux(const unsigned int & connection_index) const override;

protected:
  const Real & _f_open;
  const Real & _f_open_min;

  const ADVariableValue & _A;
  const ADVariableValue & _rhoA;
  const ADVariableValue & _rhouA;
  const ADVariableValue & _rhoEA;

  const unsigned int _rhoA_jvar;
  const unsigned int _rhouA_jvar;
  const unsigned int _rhoEA_jvar;

  const ADMaterialProperty<Real> & _p;

  const ADNumericalFlux3EqnBase & _numerical_flux;

  const std::string & _component_name;

  const MaterialProperty<RealVectorValue> & _dir;

  std::vector<std::vector<ADReal>> _solutions;
  std::vector<std::vector<ADReal>> _fluxes;
  std::vector<std::vector<dof_id_type>> _dof_indices;

  std::vector<ADReal> _stored_p;

  std::vector<unsigned int> _elem_ids;
  std::vector<unsigned int> _local_side_ids;

  std::vector<ADReal> _areas;
  std::vector<RealVectorValue> _directions;

  std::vector<unsigned int> _connection_indices;

  static const std::vector<std::pair<std::string, unsigned int>> _varname_eq_index_pairs;

public:
  static InputParameters validParams();
};