add a test for coupled elemental sln idaholab#9690

YaqiWang · Aug 30, 2017 · bd128b3 · bd128b3
1 parent 6d12923
commit bd128b3
Show file tree

Hide file tree

Showing 6 changed files with 165 additions and 0 deletions.
diff --git a/test/include/postprocessors/ElementMomentSum.h b/test/include/postprocessors/ElementMomentSum.h
@@ -0,0 +1,37 @@
+/****************************************************************/
+/*               DO NOT MODIFY THIS HEADER                      */
+/* MOOSE - Multiphysics Object Oriented Simulation Environment  */
+/*                                                              */
+/*           (c) 2010 Battelle Energy Alliance, LLC             */
+/*                   ALL RIGHTS RESERVED                        */
+/*                                                              */
+/*          Prepared by Battelle Energy Alliance, LLC           */
+/*            Under Contract No. DE-AC07-05ID14517              */
+/*            With the U. S. Department of Energy               */
+/*                                                              */
+/*            See COPYRIGHT for full restrictions               */
+/****************************************************************/
+
+#ifndef ELEMENTMOMENTSUM_H
+#define ELEMENTMOMENTSUM_H
+
+#include "ElementIntegralVariablePostprocessor.h"
+
+// Forward Declarations
+class ElementMomentSum;
+
+template <>
+InputParameters validParams<ElementMomentSum>();
+
+class ElementMomentSum : public ElementIntegralVariablePostprocessor
+{
+public:
+  ElementMomentSum(const InputParameters & parameters);
+
+protected:
+  virtual void execute() override;
+
+  const DenseVector<Number> & _elemental_sln;
+};
+
+#endif // ELEMENTMOMENTSUM_H
diff --git a/test/src/base/MooseTestApp.C b/test/src/base/MooseTestApp.C
@@ -221,6 +221,7 @@
 #include "NumAdaptivityCycles.h"
 #include "TestDiscontinuousValuePP.h"
 #include "RandomPostprocessor.h"
+#include "ElementMomentSum.h"
 #include "ChannelGradientVectorPostprocessor.h"
 
 // Functions
@@ -561,6 +562,7 @@ MooseTestApp::registerObjects(Factory & factory)
   registerPostprocessor(NumAdaptivityCycles);
   registerPostprocessor(TestDiscontinuousValuePP);
   registerPostprocessor(RandomPostprocessor);
+  registerPostprocessor(ElementMomentSum);
 
   registerVectorPostprocessor(LateDeclarationVectorPostprocessor);
   registerVectorPostprocessor(ChannelGradientVectorPostprocessor);

diff --git a/test/src/postprocessors/ElementMomentSum.C b/test/src/postprocessors/ElementMomentSum.C
@@ -0,0 +1,44 @@
+/****************************************************************/
+/*               DO NOT MODIFY THIS HEADER                      */
+/* MOOSE - Multiphysics Object Oriented Simulation Environment  */
+/*                                                              */
+/*           (c) 2010 Battelle Energy Alliance, LLC             */
+/*                   ALL RIGHTS RESERVED                        */
+/*                                                              */
+/*          Prepared by Battelle Energy Alliance, LLC           */
+/*            Under Contract No. DE-AC07-05ID14517              */
+/*            With the U. S. Department of Energy               */
+/*                                                              */
+/*            See COPYRIGHT for full restrictions               */
+/****************************************************************/
+
+#include "ElementMomentSum.h"
+
+#include "MooseVariable.h"
+
+template <>
+InputParameters
+validParams<ElementMomentSum>()
+{
+  InputParameters params = validParams<ElementIntegralVariablePostprocessor>();
+  return params;
+}
+
+ElementMomentSum::ElementMomentSum(const InputParameters & parameters)
+  : ElementIntegralVariablePostprocessor(parameters),
+    _elemental_sln(coupledElementalSln("variable"))
+{
+}
+
+void
+ElementMomentSum::execute()
+{
+  unsigned int ndofs = _elemental_sln.size();
+  Real v = 0;
+  for (unsigned int i = 0; i < ndofs; ++i)
+    v += _elemental_sln(i);
+  v /= ndofs;
+  v *= _current_elem_volume;
+
+  _integral_value += v;
+}
diff --git a/test/tests/postprocessors/coupled_elemental_sln/coupled_elemental_sln.i b/test/tests/postprocessors/coupled_elemental_sln/coupled_elemental_sln.i
@@ -0,0 +1,67 @@
+[Mesh]
+  type = GeneratedMesh
+  dim = 2
+  nx = 3
+  ny = 3
+[]
+
+[Variables]
+  [./u]
+  [../]
+[]
+
+[Kernels]
+  [./diff]
+    type = Diffusion
+    variable = u
+  [../]
+  [./time]
+    type = TimeDerivative
+    variable = u
+  [../]
+[]
+
+[BCs]
+  [./left]
+    type = DirichletBC
+    variable = u
+    boundary = left
+    value = 0
+  [../]
+  [./right]
+    type = DirichletBC
+    variable = u
+    boundary = right
+    value = 1
+  [../]
+[]
+
+[Postprocessors]
+  [./L2_norm]
+    type = ElementL2Norm
+    variable = u
+  [../]
+  [./integral]
+    type = ElementIntegralVariablePostprocessor
+    variable = u
+  [../]
+  [./direct_sum]
+    type = ElementMomentSum
+    variable = u
+  [../]
+  [./direct_sum_old]
+    type = ElementMomentSum
+    variable = u
+    implicit = false
+  [../]
+[]
+
+[Executioner]
+  type = Transient
+  num_steps = 3
+  nl_abs_tol = 1e-12
+[]
+
+[Outputs]
+  csv = true
+[]
diff --git a/test/tests/postprocessors/coupled_elemental_sln/gold/coupled_elemental_sln_out.csv b/test/tests/postprocessors/coupled_elemental_sln/gold/coupled_elemental_sln_out.csv
@@ -0,0 +1,6 @@
+time,L2_norm,direct_sum,direct_sum_old,integral
+0,0,0,0,0
+1,0.54729958835168,0.46610169491525,0,0.46610169491525
+2,0.57484469915105,0.49712726228095,0.46610169491525,0.49712726228095
+3,0.5771390732329,0.49975654765402,0.49712726228095,0.49975654765402
+
diff --git a/test/tests/postprocessors/coupled_elemental_sln/tests b/test/tests/postprocessors/coupled_elemental_sln/tests
@@ -0,0 +1,9 @@
+[Tests]
+  [./test]
+    type = 'CSVDiff'
+    input = 'coupled_elemental_sln.i'
+    csvdiff = 'coupled_elemental_sln_out.csv'
+    max_parallel = 3
+    max_threads = 3
+  [../]
+[]