Line data    Source code

       1             : //* This file is part of the MOOSE framework
       2             : //* https://mooseframework.inl.gov
       3             : //*
       4             : //* All rights reserved, see COPYRIGHT for full restrictions
       5             : //* https://github.com/idaholab/moose/blob/master/COPYRIGHT
       6             : //*
       7             : //* Licensed under LGPL 2.1, please see LICENSE for details
       8             : //* https://www.gnu.org/licenses/lgpl-2.1.html
       9             : 
      10             : #include "INSFEFluidKernelStabilization.h"
      11             : #include "libmesh/quadrature.h"
      12             : 
      13             : InputParameters
      14         861 : INSFEFluidKernelStabilization::validParams()
      15             : {
      16         861 :   InputParameters params = INSFEFluidKernelBase::validParams();
      17         861 :   return params;
      18             : }
      19             : 
      20         462 : INSFEFluidKernelStabilization::INSFEFluidKernelStabilization(const InputParameters & parameters)
      21             :   : INSFEFluidKernelBase(parameters),
      22         462 :     _u_dot(_bTransient ? _var.uDot() : _zero),
      23         462 :     _du_dot_du(_bTransient ? _var.duDotDu() : _zero),
      24         924 :     _tauc(getMaterialProperty<Real>("tauc")),
      25         924 :     _taum(getMaterialProperty<Real>("taum")),
      26        1386 :     _taue(getMaterialProperty<Real>("taue"))
      27             : {
      28         462 : }
      29             : 
      30             : void
      31    12954520 : INSFEFluidKernelStabilization::precalculateResidual()
      32             : {
      33    12954520 :   _vel_elem = RealVectorValue(0, 0, 0);
      34             :   // calculating element average velocity
      35    64772600 :   for (_qp = 0; _qp < _qrule->n_points(); _qp++)
      36             :   {
      37    51818080 :     _vel_elem(0) += _u_vel[_qp];
      38    51818080 :     _vel_elem(1) += _v_vel[_qp];
      39    51818080 :     _vel_elem(2) += _w_vel[_qp];
      40             :   }
      41    12954520 :   _vel_elem(0) = _vel_elem(0) / _qp;
      42    12954520 :   _vel_elem(1) = _vel_elem(1) / _qp;
      43    12954520 :   _vel_elem(2) = _vel_elem(2) / _qp;
      44    12954520 : }