NekRSProblem.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"

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

  static InputParameters validParams();

  ~NekRSProblem();

  virtual void initialSetup() override;

  void sendBoundaryHeatFluxToNek();

  void sendBoundaryDeformationToNek();

  void sendVolumeDeformationToNek();

  void sendVolumeHeatSourceToNek();

  void getBoundaryTemperatureFromNek();

  void getVolumeTemperatureFromNek();

  virtual void adjustNekSolution() override;

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

  virtual void addExternalVariables() override;

  virtual double maxInterpolatedTemperature() const;

  virtual double minInterpolatedTemperature() const;

  virtual const bool hasMovingNekMesh() const override { return nekrs::hasMovingMesh(); }

protected:
  bool normalizeHeatSource(const double moose_integral,
                           const double nek_integral,
                           double & normalized_nek_integral);

  void checkInitialFluxValues(const Real & nek_flux, const Real & moose_flux) const;

  virtual void addTemperatureVariable() override { return; }

  void calculateMeshVelocity(int e, const field::NekWriteEnum & field);

  const bool & _has_heat_source;

  const bool & _conserve_flux_by_sideset;

  const Real & _abs_tol;

  const Real & _rel_tol;

  const PostprocessorValue * _flux_integral = nullptr;

  const VectorPostprocessorValue * _flux_integral_vpp = nullptr;

  const PostprocessorValue * _source_integral = nullptr;

  const PostprocessorValue * _min_T = nullptr;

  const PostprocessorValue * _max_T = nullptr;

  double * _T = nullptr;

  double * _flux_face = nullptr;

  double * _flux_elem = nullptr;

  double * _source_elem = nullptr;

  double * _displacement_x = nullptr;

  double * _displacement_y = nullptr;

  double * _displacement_z = nullptr;

  double * _mesh_velocity_elem = nullptr;

  unsigned int _temp_var;

  unsigned int _avg_flux_var;

  unsigned int _disp_x_var;

  unsigned int _disp_y_var;

  unsigned int _disp_z_var;

  unsigned int _heat_source_var;
};