decoupledibpm.h

Go to the documentation of this file.

 
#pragma once
 
#include <petibm/bodypack.h>
 
#include "../navierstokes/navierstokes.h"
 
class DecoupledIBPMSolver : protected NavierStokesSolver
{
public:
    DecoupledIBPMSolver() = default;
 
    DecoupledIBPMSolver(const MPI_Comm &world, const YAML::Node &node);
 
    ~DecoupledIBPMSolver();
 
    PetscErrorCode destroy();
 
    PetscErrorCode init(const MPI_Comm &world, const YAML::Node &node);
 
    using NavierStokesSolver::ioInitialData;
 
    PetscErrorCode advance();
 
    PetscErrorCode write();
 
    using NavierStokesSolver::finished;
 
protected:
    petibm::type::BodyPack bodies;
 
    petibm::type::LinSolver fSolver;
 
    Mat H;
 
    Mat E;
 
    Mat EBNH;
 
    Mat BNH;
 
    Vec f;
 
    Vec df;
 
    Vec rhsf;
 
    PetscLogStage stageRHSForces;
 
    PetscLogStage stageSolveForces;
 
    PetscLogStage stageIntegrateForces;
 
    PetscViewer forcesViewer;
 
    virtual PetscErrorCode assembleRHSVelocity();
 
    virtual PetscErrorCode assembleRHSForces();
 
    virtual PetscErrorCode solveForces();
 
    virtual PetscErrorCode applyNoSlip();
 
    virtual PetscErrorCode updateForces();
 
    virtual PetscErrorCode createExtraOperators();
 
    virtual PetscErrorCode createExtraVectors();
 
    virtual PetscErrorCode writeRestartDataHDF5(const std::string &filePath);
 
    virtual PetscErrorCode readRestartDataHDF5(const std::string &filePath);
 
    virtual PetscErrorCode writeLinSolversInfo();
 
    virtual PetscErrorCode writeForcesASCII();
 
};  // DecoupledIBPMSolver