INSADMomentumNoBCBC.C

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* 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

#include "INSADMomentumNoBCBC.h"
#include "MooseMesh.h"
#include "INSADObjectTracker.h"
#include "NS.h"

registerMooseObject("NavierStokesApp", INSADMomentumNoBCBC);

InputParameters
INSADMomentumNoBCBC::validParams()
{
  InputParameters params = ADVectorIntegratedBC::validParams();

  params.addClassDescription("This class implements the 'No BC' boundary condition based on the "
                             "'Laplace' form of the viscous stress tensor.");
  params.addRequiredCoupledVar(NS::pressure, "pressure");
  params.addParam<bool>("integrate_p_by_parts",
                        true,
                        "Allows simulations to be run with pressure BC if set to false");
  MooseEnum viscous_form("traction laplace", "laplace");
  params.addParam<MooseEnum>("viscous_form",
                             viscous_form,
                             "The form of the viscous term. Options are 'traction' or 'laplace'");

  // Optional parameters
  params.addParam<MaterialPropertyName>("mu_name", "mu", "The name of the dynamic viscosity");
  return params;
}

INSADMomentumNoBCBC::INSADMomentumNoBCBC(const InputParameters & parameters)
  : ADVectorIntegratedBC(parameters),
    _p(adCoupledValue(NS::pressure)),
    _integrate_p_by_parts(getParam<bool>("integrate_p_by_parts")),
    _mu(getADMaterialProperty<Real>("mu_name")),
    _form(getParam<MooseEnum>("viscous_form"))
{
  std::set<SubdomainID> connected_blocks;
  for (const auto bnd_id : boundaryIDs())
  {
    const auto & these_blocks = _mesh.getBoundaryConnectedBlocks(bnd_id);
    connected_blocks.insert(these_blocks.begin(), these_blocks.end());
  }
  auto & obj_tracker = const_cast<INSADObjectTracker &>(
      _fe_problem.getUserObject<INSADObjectTracker>("ins_ad_object_tracker"));
  for (const auto block_id : connected_blocks)
  {
    obj_tracker.set("viscous_form", _form, block_id);
    obj_tracker.set("integrate_p_by_parts", _integrate_p_by_parts, block_id);
  }
}

ADReal
INSADMomentumNoBCBC::computeQpResidual()
{
  // The viscous term
  ADReal residual;
  if (_form == "laplace")
    residual = -_mu[_qp] * (_grad_u[_qp] * _normals[_qp]) * _test[_i][_qp];
  else
    residual =
        -_mu[_qp] * ((_grad_u[_qp] + _grad_u[_qp].transpose()) * _normals[_qp]) * _test[_i][_qp];

  if (_integrate_p_by_parts)
    // pIn * test
    residual += _p[_qp] * _normals[_qp] * _test[_i][_qp];

  return residual;
}