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();
};