NekRSSeparateDomainProblem.h

Go to the documentation of this file.

/********************************************************************/
/*                  SOFTWARE COPYRIGHT NOTIFICATION                 */
/*                             Cardinal                             */
/*                                                                  */
/*                  (c) 2021 UChicago Argonne, LLC                  */
/*                        ALL RIGHTS RESERVED                       */
/*                                                                  */
/*                 Prepared by UChicago Argonne, LLC                */
/*               Under Contract No. DE-AC02-06CH11357               */
/*                With the U. S. Department of Energy               */
/*                                                                  */
/*             Prepared by Battelle Energy Alliance, LLC            */
/*               Under Contract No. DE-AC07-05ID14517               */
/*                With the U. S. Department of Energy               */
/*                                                                  */
/*                 See LICENSE for full restrictions                */
/********************************************************************/

#pragma once

#include "NekRSProblemBase.h"
#include "NekTimeStepper.h"
#include "NekRSMesh.h"
#include "Transient.h"

#include <memory>

class NekRSSeparateDomainProblem : public NekRSProblemBase
{
public:
  NekRSSeparateDomainProblem(const InputParameters & params);

  static InputParameters validParams();

  ~NekRSSeparateDomainProblem();

  virtual void initialSetup() override;

  void sendBoundaryVelocityToNek(const nek_mesh::NekMeshEnum pp_mesh);

  void sendBoundaryTemperatureToNek(const nek_mesh::NekMeshEnum pp_mesh);

  void sendBoundaryScalarToNek(const nek_mesh::NekMeshEnum pp_mesh,
                               const int scalarId,
                               const double scalarValue);

  virtual void syncSolutions(ExternalProblem::Direction direction) override;

  virtual void addExternalVariables() override;

protected:
  void velocity(const nek_mesh::NekMeshEnum pp_mesh, const int elem_id, const double velocity);

  void temperature(const nek_mesh::NekMeshEnum pp_mesh, const int elem_id, const double temperature);

  void scalar(const nek_mesh::NekMeshEnum pp_mesh, const int elem_id, const int scalarId, const double scalarValue);

  const MultiMooseEnum _coupling_type;

  const std::vector<int> & _outlet_boundary;

  const std::vector<int> & _inlet_boundary;

  bool _inlet_coupling;

  bool _outlet_coupling;

  const MultiMooseEnum _coupled_scalars;

  bool _scalar01_coupling;
  bool _scalar02_coupling;
  bool _scalar03_coupling;

  const PostprocessorValue * _toNekRS_velocity = nullptr;

  const PostprocessorValue * _toNekRS_temp = nullptr;

  const PostprocessorValue * _toNekRS_scalar01 = nullptr;
  const PostprocessorValue * _toNekRS_scalar02 = nullptr;
  const PostprocessorValue * _toNekRS_scalar03 = nullptr;

  const PostprocessorValue * _transfer_in = nullptr;

  const nek_mesh::NekMeshEnum _pp_mesh = nek_mesh::fluid;
};