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); }
 
  virtual ~NavierStokesProblem();
 
  virtual void initialSetup() override;
 
protected:
  virtual void initPetscOutputAndSomeSolverSettings() override;
 
private:
  struct FieldSplitPostSetUpContext
  {
    NavierStokesProblem * problem;
    std::size_t tree_position;
  };
 
  static PetscErrorCode fieldSplitPostSetUpCallback(PC pc);
 
  void fieldSplitPostSetUp(PC pc, std::size_t tree_position);
 
  void setFieldSplitPostSetUp(PC pc, std::size_t tree_position);
 
  void setupLSCMatrices(PC schur_pc);
 
  const bool _commute_lsc;
  const TagName & _mass_matrix;
  const TagName & _L_matrix;
  const bool _have_mass_matrix;
  const bool _have_L_matrix;
 
  const enum class SetSchurPreType { FALSE, MASS, A11_AND_MASS } _set_schur_pre;
 
  const std::vector<unsigned int> & _schur_fs_index;
 
  Mat _Q_scale = nullptr;
  Mat _L = nullptr;
 
  std::vector<IS> _index_sets;
 
  std::vector<FieldSplitPostSetUpContext> _field_split_post_setup_contexts;
#endif
};