Merging for RTES

include/dirackernels/EnthalpySink.h

@@ -0,0 +1,36 @@
+//* This file is part of the MOOSE framework
+//* https://www.mooseframework.org
+//*
+//* 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 "DiracKernel.h"
+#include "PorousFlowPolyLineSink.h"
+
+class SinglePhaseFluidProperties;
+
+/**
+ * Point source (or sink) that adds (removes) fluid at a mass flux rate specified by a postprocessor
+ */
+class EnthalpySink : public PorousFlowPolyLineSink
+{
+public:
+  static InputParameters validParams();
+
+  EnthalpySink(const InputParameters & parameters);
+
+  virtual Real computeQpBaseOutflow(unsigned current_dirac_ptid) const override;
+
+protected:
+  /// Pressure
+  const VariableValue & _pressure;
+  /// Specified inlet temperature
+  const Real & _T_in;
+  /// Fluid properties UserObject
+  const SinglePhaseFluidProperties & _fp;
+};

include/dirackernels/FunctionEnthalpySink.h

@@ -0,0 +1,35 @@
+//* This file is part of the MOOSE framework
+//* https://www.mooseframework.org
+//*
+//* 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 "PorousFlowPolyLineSink.h"
+
+class SinglePhaseFluidProperties;
+class Function;
+/**
+ * Point source (or sink) that adds (removes) fluid at a mass flux rate specified by a postprocessor
+ */
+class FunctionEnthalpySink : public PorousFlowPolyLineSink
+{
+public:
+  static InputParameters validParams();
+
+  FunctionEnthalpySink(const InputParameters & parameters);
+
+  virtual Real computeQpBaseOutflow(unsigned current_dirac_ptid) const override;
+
+protected:
+  /// Pressure
+  const VariableValue & _pressure;
+   /// The function being used for evaluation
+  const Function & _func;
+  /// Fluid properties UserObject
+  const SinglePhaseFluidProperties & _fp;
+};

include/dirackernels/PointEnthalpySink.h

@@ -0,0 +1,46 @@
+//* This file is part of the MOOSE framework
+//* https://www.mooseframework.org
+//*
+//* 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 "DiracKernel.h"
+
+class SinglePhaseFluidProperties;
+
+/**
+ * Point source that adds heat energy corresponding to adding a fluid with flux rate specified by a
+ * postprocessor at given temperature (specified by a postprocessor).
+ */
+class PointEnthalpySink : public DiracKernel
+{
+public:
+  static InputParameters validParams();
+
+  PointEnthalpySink(const InputParameters & parameters);
+
+  virtual void addPoints() override;
+  virtual Real computeQpResidual() override;
+  virtual Real computeQpJacobian() override;
+  virtual Real computeQpOffDiagJacobian(unsigned int jvar) override;
+
+protected:
+  /// The constant mass flux (kg/s)
+  const PostprocessorValue & _mass_flux;
+  /// Pressure
+  const VariableValue & _pressure;
+  /// Quadpoint temperature
+  const MaterialProperty<Real> * const _temperature;
+
+  /// Fluid properties UserObject
+  const SinglePhaseFluidProperties & _fp;
+  /// The location of the point source
+  const Point _p;
+  /// Presure variable number
+  unsigned int _p_var_num;
+};

include/dirackernels/PointEnthalpySourceFromFunction.h

@@ -0,0 +1,44 @@
+//* This file is part of the MOOSE framework
+//* https://www.mooseframework.org
+//*
+//* 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 "DiracKernel.h"
+
+class SinglePhaseFluidProperties;
+
+/**
+ * Point source that adds heat energy corresponding to adding a fluid with flux rate specified by a
+ * postprocessor at given temperature (specified by a postprocessor).
+ */
+class PointEnthalpySourceFromFunction : public DiracKernel
+{
+public:
+  static InputParameters validParams();
+
+  PointEnthalpySourceFromFunction(const InputParameters & parameters);
+
+  virtual void addPoints() override;
+  virtual Real computeQpResidual() override;
+  virtual Real computeQpJacobian() override;
+  virtual Real computeQpOffDiagJacobian(unsigned int jvar) override;
+protected:
+  /// The constant mass flux (kg/s)
+  const PostprocessorValue & _mass_flux;
+  /// Pressure
+  const VariableValue & _pressure;
+   /// The function being used for inlet temperature (K)
+  const Function & _func;
+  /// Fluid properties UserObject
+  const SinglePhaseFluidProperties & _fp;
+  /// The location of the point source
+  const Point _p;
+  /// Presure variable number
+  unsigned int _p_var_num;
+};

include/dirackernels/PointEnthalpySourceSinkFromFunction.h

@@ -0,0 +1,55 @@
+//* This file is part of the MOOSE framework
+//* https://www.mooseframework.org
+//*
+//* 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 "DiracKernel.h"
+#include "PorousFlowSumQuantity.h"
+class Function;
+class SinglePhaseFluidProperties;
+/**
+ * Point source (or sink) that adds (removes) fluid at a constant mass flux rate for times
+ * between the specified start and end times. If no start and end times are specified,
+ * the source (sink) starts at the start of the simulation and continues to act indefinitely
+ */
+class PointEnthalpySourceSinkFromFunction : public DiracKernel
+{
+public:
+  static InputParameters validParams();
+
+  PointEnthalpySourceSinkFromFunction(const InputParameters & parameters);
+
+  virtual void addPoints() override;
+  virtual Real computeQpResidual() override;
+  virtual Real computeQpJacobian() override;
+  virtual Real computeQpOffDiagJacobian(unsigned int jvar) override;
+
+protected:
+  /**
+   * This is used to hold the total fluid flowing into the line sink for each time step.
+   * Hence, it is positive for production wells where fluid is flowing
+   * from porespace into the line sink (and hence removed from the model)
+   */
+  PorousFlowSumQuantity & _total_outflow_enthalpy;
+  /// The function being used for mass flux (kg/s) function
+  const Function & _func_mass_flux;
+  /// The function being used for temperature (K) function
+  const Function & _func_temperature;
+  /// Pressure
+  const VariableValue & _pressure;
+  /// Quadpoint temperature
+  const MaterialProperty<Real> * const _temperature;
+  /// Fluid properties UserObject
+  const SinglePhaseFluidProperties & _fp;
+  /// The location of the point source (sink)
+  const Point _p;
+  /// Presure variable number
+  unsigned int _p_var_num;
+
+};

include/dirackernels/PointSourceSinkFromFunction.h

@@ -0,0 +1,41 @@
+//* This file is part of the MOOSE framework
+//* https://www.mooseframework.org
+//*
+//* 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 "DiracKernel.h"
+#include "PorousFlowSumQuantity.h"
+class Function;
+/**
+ * Point source (or sink) that adds (removes) fluid at a mass flux rate specified by a postprocessor
+ */
+class PointSourceSinkFromFunction : public DiracKernel
+{
+public:
+  static InputParameters validParams();
+
+  PointSourceSinkFromFunction(const InputParameters & parameters);
+
+  virtual void addPoints() override;
+  virtual Real computeQpResidual() override;
+  // virtual Real computeQpJacobian() override;
+  // virtual Real computeQpOffDiagJacobian(unsigned int jvar) override;
+
+protected:
+  /// The function being used for evaluation mass flux
+  const Function & _func;
+  /// The location of the point source (sink)
+  const Point _p;
+  /**
+   * This is used to hold the total fluid flowing into the line sink for each time step.
+   * Hence, it is positive for production wells where fluid is flowing
+   * from porespace into the line sink (and hence removed from the model)
+   */
+  PorousFlowSumQuantity & _total_outflow_mass;
+};

include/dirackernels/PorousFlowSquarePulsePointEnthalpySink.h

@@ -0,0 +1,51 @@
+//* This file is part of the MOOSE framework
+//* https://www.mooseframework.org
+//*
+//* 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 "DiracKernel.h"
+
+class SinglePhaseFluidProperties;
+
+/**
+ * Point source that adds heat energy corresponding to adding a fluid with flux rate specified by a
+ * postprocessor at given temperature (specified by a postprocessor).
+ */
+class PorousFlowSquarePulsePointEnthalpySink : public DiracKernel
+{
+public:
+  static InputParameters validParams();
+
+  PorousFlowSquarePulsePointEnthalpySink(const InputParameters & parameters);
+
+  virtual void addPoints() override;
+  virtual Real computeQpResidual() override;
+  virtual Real computeQpJacobian() override;
+  virtual Real computeQpOffDiagJacobian(unsigned int jvar) override;
+
+protected:
+  /// The constant mass flux (kg/s)
+  const Real _mass_flux;
+  /// Pressure
+  const VariableValue & _pressure;
+  /// Quadpoint temperature
+  const MaterialProperty<Real> * const _temperature;
+
+  /// Fluid properties UserObject
+  const SinglePhaseFluidProperties & _fp;
+  /// The location of the point source
+  const Point _p;
+  /// The time at which the point source (sink) starts operating
+  const Real _start_time;
+
+  /// The time at which the point source (sink) stops operating
+  const Real _end_time;
+  /// Presure variable number
+  unsigned int _p_var_num;
+};

include/dirackernels/PorousFlowSquarePulsePointEnthalpySource.h

@@ -0,0 +1,52 @@
+//* This file is part of the MOOSE framework
+//* https://www.mooseframework.org
+//*
+//* 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 "DiracKernel.h"
+class SinglePhaseFluidProperties;
+/**
+ * Point source (or sink) that adds (removes) fluid at a constant mass flux rate for times
+ * between the specified start and end times. If no start and end times are specified,
+ * the source (sink) starts at the start of the simulation and continues to act indefinitely
+ */
+class PorousFlowSquarePulsePointEnthalpySource : public DiracKernel
+{
+public:
+  static InputParameters validParams();
+
+  PorousFlowSquarePulsePointEnthalpySource(const InputParameters & parameters);
+
+  virtual void addPoints() override;
+  virtual Real computeQpResidual() override;
+  virtual Real computeQpJacobian() override;
+  virtual Real computeQpOffDiagJacobian(unsigned int jvar) override;
+
+protected:
+  /// The constant mass flux (kg/s)
+  const Real _mass_flux;
+  /// Pressure
+  const VariableValue & _pressure;
+  /// The location of the point source (sink)
+  const Point _p;
+  /// Input temperature
+  const Real _temperature;
+
+  /// Fluid properties UserObject
+  const SinglePhaseFluidProperties & _fp;
+
+  /// The time at which the point source (sink) starts operating
+  const Real _start_time;
+
+  /// The time at which the point source (sink) stops operating
+  const Real _end_time;
+
+  /// Pressure variable number
+  unsigned int _p_var_num;
+};