NavierStokesProblem.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 "FEProblem.h"
#include "libmesh/libmesh_config.h"
#include <petscsnes.h>

class NonlinearSystem;

class NavierStokesProblem : public FEProblem
{
public:
  static InputParameters validParams();

  NavierStokesProblem(const InputParameters & parameters);

#if PETSC_RELEASE_GREATER_EQUALS(3, 20, 0)

  TagID massMatrixTagID() const { return getMatrixTagID(_mass_matrix); }

  TagID LMatrixTagID() const { return getMatrixTagID(_L_matrix); }

  void clearIndexSets() { _index_sets.clear(); }

  /*
   * Given a \p KSP \p node and where we are in the field split tree, given by \p tree_position,
   * return the next \p KSP object in the tree on the way to the Schur complement \p KSP object.
   * Each invocation of this method moves through one level of our \p _index_sets data member. This
   * method will call itself recursively until it reaches the Schur complement \p KSP
   */
  KSP findSchurKSP(KSP node, unsigned int tree_position);

  void setupLSCMatrices(KSP schur_ksp);

  virtual ~NavierStokesProblem();

protected:
  virtual void initPetscOutputAndSomeSolverSettings() override;

private:
  const bool _commute_lsc;
  const TagName & _mass_matrix;
  const TagName & _L_matrix;
  const bool _have_mass_matrix;
  const bool _have_L_matrix;

  const bool _pressure_mass_matrix_as_pre;

  const std::vector<unsigned int> & _schur_fs_index;

  Mat _Q_scale = nullptr;
  Mat _L = nullptr;

  std::vector<IS> _index_sets;
#endif
};