SPIbaseflow.hpp

Go to the documentation of this file.

#ifndef SPIBASEFLOW_H
#define SPIBASEFLOW_H
#include <iostream>
#include <vector>
#include <petscksp.h>
#include <slepceps.h>
#include <slepcpep.h>
#include <string>
#include <tuple>
#include "SPIgrid.hpp"
#include "SPIparams.hpp"
#include "SPIMat.hpp"
#include "SPIVec.hpp"
#include "SPIprint.hpp"
 
namespace SPI{
    struct SPIbaseflow{
        SPIbaseflow(std::string _name="baseflow");                 // constructor with no arguments
        SPIbaseflow(
                SPIVec U,
                SPIVec V,
                SPIVec Ux,
                SPIVec Uxx, // added for streamwise viscous marching BL
                SPIVec Uy,
                SPIVec Uxy,
                SPIVec Vx, // added for full terms
                SPIVec Vxx, // added for streamwise viscous marching BL
                SPIVec Vy,
                SPIVec W,
                SPIVec Wx,
                SPIVec Wy,
                SPIVec Wxy, 
                SPIVec P,
                std::string _name="baseflow");                  // constructor with baseflow values
        std::string name;                                       
        PetscErrorCode ierr;                                    
        PetscBool flag_init=PETSC_FALSE;                        
        PetscInt print();                                       // print baseflow to screen
        SPIVec U,       // streamwise baseflow
               V,       // wall-normal baseflow
               Ux,      // streamwise baseflow derivative with respect to streamwise
               Uxx,     // streamwise baseflow second derivative with respect to streamwise
               Uy,      // streamwise baseflow derivative with respect to wall-normal
               Uxy,     // streamwise baseflow mixed derivative
               Vx,      // wall-normal baseflow derivative with respect to streamwise
               Vxx,     // wall-normal baseflow second derivative with respect to streamwise
               Vy,      // wall-normal baseflow derivative with respect to wall-normal
               W,       // spanwise baseflow
               Wx,      // spanwise baseflow derivative with respect to streamwise
               Wy,      // spanwise baseflow derivative with respect to wall-normal
               Wxy,     // spanwise baseflow mixed derivative
               P;       //  pressure baseflow
        ~SPIbaseflow(); // destructor to delete memory
    };
    SPIbaseflow blasius(SPIparams &params, SPIgrid1D &grid);
    int _bblf( const PetscScalar input[3], PetscScalar output[3]); // Blasius boundary layer flow ODE
    SPIbaseflow channel(SPIparams &params, SPIgrid1D &grid);
}
#endif // SPIBASEFLOW_H