ibpm.h

Go to the documentation of this file.

 
#pragma once
 
#include <petibm/bodypack.h>
 
#include "../navierstokes/navierstokes.h"
 
class IBPMSolver : protected NavierStokesSolver
{
public:
    IBPMSolver() = default;
 
    IBPMSolver(const MPI_Comm &world, const YAML::Node &node);
 
    ~IBPMSolver();
 
    PetscErrorCode destroy();
 
    PetscErrorCode init(const MPI_Comm &world, const YAML::Node &node);
 
    using NavierStokesSolver::ioInitialData;
 
    using NavierStokesSolver::advance;
 
    PetscErrorCode write();
 
    using NavierStokesSolver::finished;
 
protected:
    petibm::type::BodyPack bodies;
 
    Vec P;
 
    IS isDE[2];
 
    PetscLogStage stageIntegrateForces;
 
    PetscViewer forcesViewer;
 
    virtual PetscErrorCode assembleRHSPoisson();
 
    virtual PetscErrorCode createOperators();
 
    virtual PetscErrorCode createVectors();
 
    virtual PetscErrorCode setNullSpace();
 
    virtual PetscErrorCode writeSolutionHDF5(const std::string &filePath);
 
    virtual PetscErrorCode writeRestartDataHDF5(const std::string &filePath);
 
    virtual PetscErrorCode readRestartDataHDF5(const std::string &filePath);
 
    virtual PetscErrorCode writeForcesASCII();
 
};  // IBPMSolver