SIMPLESolveNonlinearAssembly.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

// Moose includes
#include "RhieChowMassFlux.h"
#include "SIMPLESolveBase.h"
#include "INSFVRhieChowInterpolatorSegregated.h"

class SIMPLESolveNonlinearAssembly : public SIMPLESolveBase
{
public:
  SIMPLESolveNonlinearAssembly(Executioner & ex);

  static InputParameters validParams();

  virtual void linkRhieChowUserObject() override;

  virtual void checkIntegrity() override;

  virtual bool solve() override;

protected:
  virtual std::vector<std::pair<unsigned int, Real>> solveMomentumPredictor() override;
  virtual std::pair<unsigned int, Real> solvePressureCorrector() override;

  virtual void checkTimeKernels(NonlinearSystemBase & system);

  std::pair<unsigned int, Real> solveAdvectedSystem(const unsigned int system_num,
                                                    NonlinearSystemBase & system,
                                                    const Real relaxation_factor,
                                                    libMesh::SolverConfiguration & solver_config,
                                                    const Real abs_tol);

  std::pair<unsigned int, Real> solveSolidEnergySystem();

  std::vector<unsigned int> _momentum_system_numbers;

  std::vector<NonlinearSystemBase *> _momentum_systems;

  const unsigned int _pressure_sys_number;

  NonlinearSystemBase & _pressure_system;

  // ************************ For quick equation check ********************* //

  const bool _has_turbulence_systems;

  // ************************ Energy Eq Variables ************************** //

  const unsigned int _energy_sys_number;

  NonlinearSystemBase * _energy_system;

  // ********************* Solid Energy Eq Variables *********************** //

  const unsigned int _solid_energy_sys_number;

  NonlinearSystemBase * _solid_energy_system;

  // ********************* Passive Scalar Eq. Variables *********************** //

  std::vector<NonlinearSystemBase *> _passive_scalar_systems;

  // ********************* Turbulence Eq Variables ************************* //

  const std::vector<SolverSystemName> & _turbulence_system_names;

  // The number(s) of the system(s) corresponding to the turbulence equation(s)
  std::vector<unsigned int> _turbulence_system_numbers;

  std::vector<NonlinearSystemBase *> _turbulence_systems;

  const std::vector<Real> _turbulence_equation_relaxation;

  std::vector<Real> _turbulence_field_min_limit;

  Moose::PetscSupport::PetscOptions _turbulence_petsc_options;

  SIMPLESolverConfiguration _turbulence_linear_control;

  const Real _turbulence_l_abs_tol;

  // ********************* SIMPLE iteration variables ********************** //

  const std::vector<Real> _turbulence_absolute_tolerance;

  // *************************** Other variables *************************** //

  INSFVRhieChowInterpolatorSegregated * _rc_uo;

  const TagName _pressure_tag_name;

  const TagID _pressure_tag_id;
};