ADJunctionOneToOne1PhaseUserObject.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"
#include "SlopeReconstruction1DInterface.h"
 
class ADNumericalFlux3EqnBase;
class SinglePhaseFluidProperties;
 
class ADJunctionOneToOne1PhaseUserObject : public ADFlowJunctionUserObject,
                                           public SlopeReconstruction1DInterface<true>
{
public:
  ADJunctionOneToOne1PhaseUserObject(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;
 
  virtual std::vector<ADReal> computeElementPrimitiveVariables(const Elem * elem) const override;
 
protected:
  const ADVariableValue & _A_linear;
 
  // piecewise constant conserved variables
  MooseVariable * _A_var;
  MooseVariable * _rhoA_var;
  MooseVariable * _rhouA_var;
  MooseVariable * _rhoEA_var;
  std::vector<MooseVariable *> _U_vars;
 
  const SinglePhaseFluidProperties & _fp;
 
  const ADNumericalFlux3EqnBase & _numerical_flux;
 
  const std::string & _junction_name;
 
  const MaterialProperty<RealVectorValue> & _dir;
 
  std::vector<std::vector<ADReal>> _primitive_solutions;
  std::vector<std::vector<ADReal>> _neighbor_primitive_solutions;
 
  std::vector<std::vector<ADReal>> _fluxes;
  std::vector<std::vector<dof_id_type>> _dof_indices;
 
  std::vector<unsigned int> _elem_ids;
  std::vector<unsigned int> _local_side_ids;
 
  std::vector<ADReal> _areas_linear;
  std::vector<RealVectorValue> _directions;
  std::vector<Point> _positions;
  std::vector<Point> _neighbor_positions;
  std::vector<Real> _delta_x;
  std::vector<bool> _has_neighbor;
 
  std::vector<unsigned int> _connection_indices;
 
  static const std::vector<std::pair<std::string, unsigned int>> _varname_eq_index_pairs;
 
  static Threads::spin_mutex _spin_mutex;
 
public:
  static InputParameters validParams();
};