INSFEFluidKernelBase.h

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

#pragma once

#include "Kernel.h"
#include "SinglePhaseFluidProperties.h"
#include "MooseVariable.h"

class INSFEFluidKernelBase : public Kernel
{
public:
  static InputParameters validParams();

  INSFEFluidKernelBase(const InputParameters & parameters);
  virtual ~INSFEFluidKernelBase() {}

protected:
  Real velocityDiv() { return _grad_u_vel[_qp](0) + _grad_v_vel[_qp](1) + _grad_w_vel[_qp](2); }
  RealVectorValue velocityDot() const;

  const VariableSecond & _second_u;
  // Coupled variables
  MooseVariable & _u_var;
  MooseVariable & _v_var;
  MooseVariable & _w_var;

  const VariableValue & _u_vel;
  const VariableValue & _v_vel;
  const VariableValue & _w_vel;
  const VariableValue & _pressure;
  const VariableValue & _temperature;

  const MaterialProperty<Real> & _rho;
  bool _has_porosity;
  const VariableValue & _porosity;

  bool _bTransient;
  const VariableValue & _u_vel_dot;
  const VariableValue & _v_vel_dot;
  const VariableValue & _w_vel_dot;

  // Gradients
  const VariableGradient & _grad_u_vel;
  const VariableGradient & _grad_v_vel;
  const VariableGradient & _grad_w_vel;
  const VariableGradient & _grad_pressure;
  const VariableGradient & _grad_temperature;

  // Variable numberings
  unsigned _u_vel_var_number;
  unsigned _v_vel_var_number;
  unsigned _w_vel_var_number;
  unsigned _pressure_var_number;
  unsigned _temperature_var_number;

  // Material properties
  const MaterialProperty<RealTensorValue> & _viscous_stress_tensor;
  const MaterialProperty<Real> & _dynamic_viscosity;
  const MaterialProperty<Real> & _turbulence_viscosity;
  const MaterialProperty<RealTensorValue> & _inertia_resistance_coeff;
  const MaterialProperty<RealTensorValue> & _viscous_resistance_coeff;

  unsigned int mapVarNumber(unsigned int var) const;

  const SinglePhaseFluidProperties & _eos;

  RealVectorValue _vec_g;
};

inline RealVectorValue
INSFEFluidKernelBase::velocityDot() const
{
  RealVectorValue vec_vel_dot(_u_vel_dot[_qp], _v_vel_dot[_qp], _w_vel_dot[_qp]);
  return vec_vel_dot;
}