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