NekBoundaryFlux.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 "ConservativeFieldTransfer.h"
 
class NekBoundaryFlux : public ConservativeFieldTransfer
{
public:
  static InputParameters validParams();
 
  NekBoundaryFlux(const InputParameters & parameters);
 
  virtual void sendDataToNek() override;
 
  virtual void readDataFromNek() override;
 
protected:
  void checkInitialFluxValues(const Real & nek_flux, const Real & moose_flux) const;
 
  bool
  normalizeFlux(const double moose_integral, double nek_integral, double & normalized_nek_integral);
 
  bool normalizeFluxBySideset(const std::vector<double> & moose_integral,
                              std::vector<double> & nek_integral,
                              double & normalized_nek_integral);
 
  const bool & _conserve_flux_by_sideset;
 
  const Real & _initial_flux_integral;
 
  const std::vector<int> * _boundary;
 
  const Real _reference_flux_integral;
 
  const PostprocessorValue * _flux_integral = nullptr;
 
  const VectorPostprocessorValue * _flux_integral_vpp = nullptr;
};