grid_movement_structure.hpp

Go to the documentation of this file.

#pragma once

#include "./mpi_structure.hpp"

#include <iostream>
#include <cstdlib>
#include <fstream>
#include <cmath>
#include <ctime>

#include "geometry_structure.hpp"
#include "config_structure.hpp"
#include "matrix_structure.hpp"
#include "vector_structure.hpp"
#include "linear_solvers_structure.hpp"

using namespace std;

class CGridMovement {
public:

    CGridMovement(void);

         virtual ~CGridMovement(void);
  
    virtual void SetSurface_Deformation(CGeometry *geometry, CConfig *config);
  
};

class CFreeFormBlending {

protected:
  unsigned short Order, 
   Degree,              
   nControl;            
public:

  CFreeFormBlending();

  virtual ~CFreeFormBlending();

  virtual su2double GetBasis(short val_i, su2double val_t);

  virtual su2double GetDerivative(short val_i, su2double val_t, short val_order);

  virtual void SetOrder(short val_order, short n_controlpoints);

  su2double GetOrder();

  su2double GetDegree();
};

class CBSplineBlending : public CFreeFormBlending{

private:
  vector<su2double>          U;  
  vector<vector<su2double> > N;  
  unsigned short KnotSize;       
public:

  CBSplineBlending(short val_order, short n_controlpoints);

  ~CBSplineBlending();

  su2double GetBasis(short val_i, su2double val_t);

  su2double GetDerivative(short val_i, su2double val_t, short val_order_der);

  void SetOrder(short val_order, short n_controlpoints);

};

class CBezierBlending : public CFreeFormBlending{

private:

  vector<su2double> binomial; 
  su2double GetBernstein(short val_n, short val_i, su2double val_t);

  su2double GetBernsteinDerivative(short val_n, short val_i, su2double val_t, short val_order_der);

  su2double Binomial(unsigned short n, unsigned short m);

public:

  CBezierBlending(short val_order, short n_controlpoints);

  ~CBezierBlending();

  su2double GetBasis(short val_i, su2double val_t);

  su2double GetDerivative(short val_i, su2double val_t, short val_order_der);

  void SetOrder(short val_order, short n_controlpoints);

};

class CFreeFormDefBox : public CGridMovement {
public:
    unsigned short nDim;                  
    unsigned short nCornerPoints,         
    nControlPoints, nControlPoints_Copy;  
    su2double **Coord_Corner_Points,        
    ****Coord_Control_Points,               
    ****ParCoord_Control_Points,        
    ****Coord_Control_Points_Copy,  
    ****Coord_SupportCP;                        
    unsigned short lOrder, lOrder_Copy, 
    mOrder, mOrder_Copy,                                
    nOrder, nOrder_Copy;                                    
    unsigned short lDegree, lDegree_Copy, 
    mDegree, mDegree_Copy,                              
    nDegree, nDegree_Copy;                              
    su2double *ParamCoord, *ParamCoord_,    
    *cart_coord, *cart_coord_;          
  su2double ObjFunc;            
    su2double *Gradient;            
  su2double **Hessian;    
    su2double MaxCoord[3];      
    su2double MinCoord[3];      
    string Tag;                     
    unsigned short Level;                               
    vector<su2double> CartesianCoord[3];        
    vector<su2double> ParametricCoord[3];   
    vector<unsigned short> MarkerIndex; 
    vector<unsigned long> VertexIndex;  
    vector<unsigned long> PointIndex;       
    unsigned long nSurfacePoint;                
    vector<string> ParentFFDBox;                    
    vector<string> ChildFFDBox;                 
  vector<unsigned short> Fix_IPlane;  
  vector<unsigned short> Fix_JPlane;  
  vector<unsigned short> Fix_KPlane;  
  CFreeFormBlending** BlendingFunction;


public:
    
    CFreeFormDefBox(void);
    
  CFreeFormDefBox(unsigned short Degree[], unsigned short BSplineOrder[], unsigned short kind_blending);

    ~CFreeFormDefBox(void);
    
  void Set_Fix_IPlane(unsigned short val_plane);

  void Set_Fix_JPlane(unsigned short val_plane);
  
  void Set_Fix_KPlane(unsigned short val_plane);
  
  unsigned short Get_Fix_IPlane(unsigned short val_index);
  
  unsigned short Get_Fix_JPlane(unsigned short val_index);
  
  unsigned short Get_Fix_KPlane(unsigned short val_index);
  
  unsigned short Get_nFix_IPlane(void);
  
  unsigned short Get_nFix_JPlane(void);
  
  unsigned short Get_nFix_KPlane(void);
  
    void Set_MarkerIndex(unsigned short val_iMarker);
    
    void Set_VertexIndex(unsigned long val_iVertex);
    
    void Set_PointIndex(unsigned long val_iPoint);
    
    void Set_CartesianCoord(su2double *val_coord);
    
    void Set_ParametricCoord(su2double *val_coord);
    
    void SetParentFFDBox(string val_iParentFFDBox);
    
    void SetChildFFDBox(string val_iChildFFDBox);
    
    void Set_CartesianCoord(su2double *val_coord, unsigned long val_iSurfacePoints);
    
    void Set_ParametricCoord(su2double *val_coord, unsigned long val_iSurfacePoints);

    unsigned short Get_MarkerIndex(unsigned long val_iSurfacePoints);
    
    unsigned long Get_VertexIndex(unsigned long val_iSurfacePoints);
    
    unsigned long Get_PointIndex(unsigned long val_iSurfacePoints);
    
    su2double *Get_CartesianCoord(unsigned long val_iSurfacePoints);
    
    su2double *Get_ParametricCoord(unsigned long val_iSurfacePoints);
    
    unsigned long GetnSurfacePoint(void);
    
    unsigned short GetnParentFFDBox(void);
    
    unsigned short GetnChildFFDBox(void);
    
    string GetParentFFDBoxTag(unsigned short val_ParentFFDBox);
    
    string GetChildFFDBoxTag(unsigned short val_ChildFFDBox);
    
    void SetSupportCPChange(CFreeFormDefBox *FFDBox);
    
    void SetnCornerPoints(unsigned short val_ncornerpoints);
    
    unsigned short GetnCornerPoints(void);
    
    unsigned short GetnControlPoints(void);
  
    void SetnControlPoints(void);
    
    unsigned long GetnSurfacePoints(void);
    
    void SetUnitCornerPoints(void);
    
    void SetCoordCornerPoints(su2double *val_coord, unsigned short val_icornerpoints);

    void SetCoordCornerPoints(su2double val_xcoord, su2double val_ycoord, su2double val_zcoord, unsigned short val_icornerpoints);
    
    void SetCoordControlPoints(su2double *val_coord, unsigned short iDegree, unsigned short jDegree, unsigned short kDegree);   

    void SetCoordControlPoints_Copy(su2double *val_coord, unsigned short iDegree, unsigned short jDegree, unsigned short kDegree);

    void SetParCoordControlPoints(su2double *val_coord, unsigned short iDegree, unsigned short jDegree, unsigned short kDegree);    
    
    su2double GetCoordCornerPoints(unsigned short val_dim, unsigned short val_icornerpoints);
    
    su2double *GetCoordCornerPoints(unsigned short val_icornerpoints);
    
    su2double *GetCoordControlPoints(unsigned short val_iindex, unsigned short val_jindex, unsigned short val_kindex);
    
    su2double *GetParCoordControlPoints(unsigned short val_iindex, unsigned short val_jindex, unsigned short val_kindex);
    
    void SetControlPoints_Parallelepiped(void);
    
    void SetSupportCP(CFreeFormDefBox *FFDBox);

    void SetControlPoints(unsigned short *val_index, su2double *movement);
    
    void SetOriginalControlPoints(void);
    
    void SetTecplot(CGeometry *geometry, unsigned short iFFDBox, bool original);
        
  void SetParaview(CGeometry *geometry, unsigned short iFFDBox, bool original);

  void SetCyl2Cart_ControlPoints(CConfig *config);
  
  void SetCart2Cyl_ControlPoints(CConfig *config);
  
  void SetCyl2Cart_CornerPoints(CConfig *config);
  
  void SetCart2Cyl_CornerPoints(CConfig *config);
  
  void SetSphe2Cart_ControlPoints(CConfig *config);
  
  void SetCart2Sphe_ControlPoints(CConfig *config);
  
  void SetSphe2Cart_CornerPoints(CConfig *config);
  
  void SetCart2Sphe_CornerPoints(CConfig *config);

    su2double *GetParametricCoord_Analytical(su2double *cart_coord);
    
    su2double *GetParametricCoord_Iterative(unsigned long iPoint, su2double *xyz, su2double *guess, CConfig *config);
    
    void CrossProduct(su2double *v1, su2double *v2, su2double *v3);
    
    su2double DotProduct(su2double *v1, su2double *v2);
    
    su2double *EvalCartesianCoord(su2double *ParamCoord);
    
    unsigned short GetlOrder(void);
    
    unsigned short GetmOrder(void);
    
    unsigned short GetnOrder(void);
  
    void SetlOrder(unsigned short val_lOrder);
    
    void SetmOrder(unsigned short val_mOrder);
    
    void SetnOrder(unsigned short val_nOrder);
    
    bool GetPointFFD(CGeometry *geometry, CConfig *config, unsigned long iPoint);
    
    void SetDeformationZone(CGeometry *geometry, CConfig *config, unsigned short iFFDBox);
    
    su2double *GetFFDGradient(su2double *val_coord, su2double *xyz);
    
  void GetFFDHessian(su2double *uvw, su2double *xyz, su2double **val_Hessian);
  
  su2double GetDerivative1(su2double *uvw, unsigned short val_diff, unsigned short *ijk, unsigned short *lmn);
    
    su2double GetDerivative2(su2double *uvw, unsigned short dim, su2double *xyz, unsigned short *lmn);
    
    su2double GetDerivative3(su2double *uvw, unsigned short dim, unsigned short diff_this,
                          unsigned short *lmn);
    
  su2double GetDerivative4(su2double *uvw, unsigned short val_diff, unsigned short val_diff2,
               unsigned short *ijk, unsigned short *lmn);
    
  su2double GetDerivative5(su2double *uvw, unsigned short dim, unsigned short diff_this, unsigned short diff_this_also,
              unsigned short *lmn);
    
    su2double GetNorm(su2double *a);
    
    void SetTag(string val_tag);
    
    string GetTag(void);
    
    void SetLevel(unsigned short val_level);
    
    unsigned short GetLevel(void);
  
    su2double Determinant_3x3(su2double A00, su2double A01, su2double A02, su2double A10, su2double A11,
                         su2double A12, su2double A20, su2double A21, su2double A22);
  
};

class CVolumetricMovement : public CGridMovement {
protected:

    unsigned short nDim;        
    unsigned short nVar;        
    unsigned long nPoint;       
    unsigned long nPointDomain;     
    unsigned long nIterMesh;    
  CSysMatrix StiffMatrix; 
  CSysVector LinSysSol;
  CSysVector LinSysRes;

public:

    CVolumetricMovement(CGeometry *geometry, CConfig *config);
    
    ~CVolumetricMovement(void);
  
    void UpdateGridCoord(CGeometry *geometry, CConfig *config);
  
    void UpdateDualGrid(CGeometry *geometry, CConfig *config);
  
    void UpdateMultiGrid(CGeometry **geometry, CConfig *config);
  
    su2double SetFEAMethodContributions_Elem(CGeometry *geometry, CConfig *config);
  
  void SetFEA_StiffMatrix3D(CGeometry *geometry, CConfig *config, su2double **StiffMatrix_Elem, unsigned long PointCorners[8], su2double CoordCorners[8][3],
                            unsigned short nNodes, su2double ElemVolume, su2double ElemDistance);
  
  void SetFEA_StiffMatrix2D(CGeometry *geometry, CConfig *config, su2double **StiffMatrix_Elem, unsigned long PointCorners[8], su2double CoordCorners[8][3],
                            unsigned short nNodes, su2double ElemVolume, su2double ElemDistance);
    
  su2double ShapeFunc_Hexa(su2double Xi, su2double Eta, su2double Zeta, su2double CoordCorners[8][3], su2double DShapeFunction[8][4]);
  
  su2double ShapeFunc_Tetra(su2double Xi, su2double Eta, su2double Zeta, su2double CoordCorners[8][3], su2double DShapeFunction[8][4]);
  
  su2double ShapeFunc_Pyram(su2double Xi, su2double Eta, su2double Zeta, su2double CoordCorners[8][3], su2double DShapeFunction[8][4]);
  
  su2double ShapeFunc_Prism(su2double Xi, su2double Eta, su2double Zeta, su2double CoordCorners[8][3], su2double DShapeFunction[8][4]);
  
  su2double ShapeFunc_Triangle(su2double Xi, su2double Eta, su2double CoordCorners[8][3], su2double DShapeFunction[8][4]);
  
  su2double ShapeFunc_Quadrilateral(su2double Xi, su2double Eta, su2double CoordCorners[8][3], su2double DShapeFunction[8][4]);
  
  su2double GetHexa_Volume(su2double CoordCorners[8][3]);
  
  su2double GetTetra_Volume(su2double CoordCorners[8][3]);
  
  su2double GetPrism_Volume(su2double CoordCorners[8][3]);
  
  su2double GetPyram_Volume(su2double CoordCorners[8][3]);
  
  su2double GetTriangle_Area(su2double CoordCorners[8][3]);
  
  su2double GetQuadrilateral_Area(su2double CoordCorners[8][3]);
    
  void AddFEA_StiffMatrix(CGeometry *geometry, su2double **StiffMatrix_Elem, unsigned long PointCorners[8], unsigned short nNodes);
  
  void ComputeDeforming_Element_Volume(CGeometry *geometry, su2double &MinVolume, su2double &MaxVolume);
  
  
  void ComputeDeforming_Wall_Distance(CGeometry *geometry, CConfig *config, su2double &MinDistance, su2double &MaxDistance);
    
    void SetBoundaryDisplacements(CGeometry *geometry, CConfig *config);
    
    void SetDomainDisplacements(CGeometry *geometry, CConfig *config);
  
    void Rigid_Rotation(CGeometry *geometry, CConfig *config, unsigned short iZone, unsigned long iter);
    
    void Rigid_Pitching(CGeometry *geometry, CConfig *config, unsigned short iZone, unsigned long iter);
  
    void Rigid_Plunging(CGeometry *geometry, CConfig *config, unsigned short iZone, unsigned long iter);
  
    void Rigid_Translation(CGeometry *geometry, CConfig *config, unsigned short iZone, unsigned long iter);
  
  void SetVolume_Scaling(CGeometry *geometry, CConfig *config, bool UpdateGeo);
  
  void SetVolume_Translation(CGeometry *geometry, CConfig *config, bool UpdateGeo);
  
  void SetVolume_Rotation(CGeometry *geometry, CConfig *config, bool UpdateGeo);
  
  void SetVolume_Deformation(CGeometry *geometry, CConfig *config, bool UpdateGeo, bool Derivative = false);

  void SetBoundaryDerivatives(CGeometry *geometry, CConfig *config);

  void UpdateGridCoord_Derivatives(CGeometry *geometry, CConfig *config);

    su2double Determinant_3x3(su2double A00, su2double A01, su2double A02, su2double A10, su2double A11, su2double A12, su2double A20, su2double A21, su2double A22);


    void Set_nIterMesh(unsigned long val_nIterMesh);

    unsigned long Get_nIterMesh(void);
};

class CSurfaceMovement : public CGridMovement {
protected:
  CFreeFormDefBox** FFDBox; 
    unsigned short nFFDBox; 
    unsigned short nLevel;  
    bool FFDBoxDefinition;  
  vector<su2double> GlobalCoordX[MAX_NUMBER_FFD];
  vector<su2double> GlobalCoordY[MAX_NUMBER_FFD];
  vector<su2double> GlobalCoordZ[MAX_NUMBER_FFD];
  vector<string> GlobalTag[MAX_NUMBER_FFD];
  vector<unsigned long> GlobalPoint[MAX_NUMBER_FFD];

public:
    
    CSurfaceMovement(void);
    
    ~CSurfaceMovement(void);
  
    void SetHicksHenne(CGeometry *boundary, CConfig *config, unsigned short iDV, bool ResetDef);
  
    void SetSurface_Bump(CGeometry *boundary, CConfig *config, unsigned short iDV, bool ResetDef);

  void SetAngleOfAttack(CGeometry *boundary, CConfig *config, unsigned short iDV, bool ResetDef);

    void SetCST(CGeometry *boundary, CConfig *config, unsigned short iDV, bool ResetDef);

    void SetNACA_4Digits(CGeometry *boundary, CConfig *config);
    
    void SetParabolic(CGeometry *boundary, CConfig *config);
    
    void SetAirfoil(CGeometry *boundary, CConfig *config);
  
    void SetRotation(CGeometry *boundary, CConfig *config, unsigned short iDV, bool ResetDef);

    void Moving_Walls(CGeometry *geometry, CConfig *config, unsigned short iZone, unsigned long iter);
  
    void Surface_Translating(CGeometry *geometry, CConfig *config,
                        unsigned long iter, unsigned short iZone);
  
    void Surface_Plunging(CGeometry *geometry, CConfig *config,
                           unsigned long iter, unsigned short iZone);
  
    void Surface_Pitching(CGeometry *geometry, CConfig *config,
                             unsigned long iter, unsigned short iZone);
  
    void Surface_Rotating(CGeometry *geometry, CConfig *config,
                        unsigned long iter, unsigned short iZone);
  
  void HTP_Rotation(CGeometry *geometry, CConfig *config,
                    unsigned long iter, unsigned short iZone);

    void AeroelasticDeform(CGeometry *geometry, CConfig *config, unsigned long ExtIter, unsigned short iMarker, unsigned short iMarker_Monitoring, vector<su2double>& displacements);
    
    void SetBoundary_Flutter3D(CGeometry *geometry, CConfig *config, 
                             CFreeFormDefBox **FFDBox, unsigned long iter, unsigned short iZone);
    
  void SetCollective_Pitch(CGeometry *geometry, CConfig *config);
  
  void SetExternal_Deformation(CGeometry *geometry, CConfig *config, unsigned short iZone, unsigned long iter);
  
    void SetTranslation(CGeometry *boundary, CConfig *config, unsigned short iDV, bool ResetDef);

  void SetScale(CGeometry *boundary, CConfig *config, unsigned short iDV, bool ResetDef);

    void CopyBoundary(CGeometry *geometry, CConfig *config);
  
    void SetSurface_Deformation(CGeometry *geometry, CConfig *config);
    
    void SetParametricCoord(CGeometry *geometry, CConfig *config, CFreeFormDefBox *FFDBox, unsigned short iFFDBox);
    
    void UpdateParametricCoord(CGeometry *geometry, CConfig *config, CFreeFormDefBox *FFDBox, unsigned short iFFDBox);
    
  void CheckFFDIntersections(CGeometry *geometry, CConfig *config, CFreeFormDefBox *FFDBox, unsigned short iFFDBox);
  
  void CheckFFDDimension(CGeometry *geometry, CConfig *config, CFreeFormDefBox *FFDBox, unsigned short iFFDBox);

    void SetParametricCoordCP(CGeometry *geometry, CConfig *config, CFreeFormDefBox *FFDBoxParent, CFreeFormDefBox *FFDBoxChild);
    
    void GetCartesianCoordCP(CGeometry *geometry, CConfig *config, CFreeFormDefBox *FFDBoxParent, CFreeFormDefBox *FFDBoxChild);

  void SetCartesianCoord(CGeometry *geometry, CConfig *config, CFreeFormDefBox *FFDBox, unsigned short iFFDBox, bool ResetDef);
    
  bool SetFFDCPChange_2D(CGeometry *geometry, CConfig *config, CFreeFormDefBox *FFDBox, CFreeFormDefBox **ResetFFDBox, unsigned short iDV, bool ResetDef);
  
  bool SetFFDCPChange(CGeometry *geometry, CConfig *config, CFreeFormDefBox *FFDBox, CFreeFormDefBox **ResetFFDBox, unsigned short iDV, bool ResetDef);
  
  bool SetFFDGull(CGeometry *geometry, CConfig *config, CFreeFormDefBox *FFDBox, CFreeFormDefBox **ResetFFDBox, unsigned short iDV, bool ResetDef);
  
  bool SetFFDNacelle(CGeometry *geometry, CConfig *config, CFreeFormDefBox *FFDBox, CFreeFormDefBox **ResetFFDBox, unsigned short iDV, bool ResetDef);
  
  bool SetFFDCamber_2D(CGeometry *geometry, CConfig *config, CFreeFormDefBox *FFDBox, CFreeFormDefBox **ResetFFDBox, unsigned short iDV, bool ResetDef);
  
  bool SetFFDTwist_2D(CGeometry *geometry, CConfig *config, CFreeFormDefBox *FFDBox, CFreeFormDefBox **ResetFFDBox, unsigned short iDV, bool ResetDef);
  
  bool SetFFDThickness_2D(CGeometry *geometry, CConfig *config, CFreeFormDefBox *FFDBox, CFreeFormDefBox **ResetFFDBox, unsigned short iDV, bool ResetDef);
  
  bool SetFFDCamber(CGeometry *geometry, CConfig *config, CFreeFormDefBox *FFDBox, CFreeFormDefBox **ResetFFDBox, unsigned short iDV, bool ResetDef);
  
  bool SetFFDThickness(CGeometry *geometry, CConfig *config, CFreeFormDefBox *FFDBox, CFreeFormDefBox **ResetFFDBox, unsigned short iDV, bool ResetDef);
  
  void SetFFDAngleOfAttack(CGeometry *geometry, CConfig *config, CFreeFormDefBox *FFDBox, CFreeFormDefBox **ResetFFDBox, unsigned short iDV, bool ResetDef);
  
  bool SetFFDTwist(CGeometry *geometry, CConfig *config, CFreeFormDefBox *FFDBox, CFreeFormDefBox **ResetFFDBox, unsigned short iDV, bool ResetDef);
  
  bool SetFFDRotation(CGeometry *geometry, CConfig *config, CFreeFormDefBox *FFDBox,CFreeFormDefBox **ResetFFDBox, unsigned short iDV, bool ResetDef);
  
  bool SetFFDControl_Surface(CGeometry *geometry, CConfig *config, CFreeFormDefBox *FFDBox, CFreeFormDefBox **ResetFFDBox, unsigned short iDV, bool ResetDef);
    
    void ReadFFDInfo(CGeometry *geometry, CConfig *config, CFreeFormDefBox **FFDBox, string val_mesh_filename);
    
  void ReadFFDInfo(CGeometry *geometry, CConfig *config, CFreeFormDefBox **FFDBox);
  
  void MergeFFDInfo(CGeometry *geometry, CConfig *config);
  
    void WriteFFDInfo(CGeometry *geometry, CConfig *config);
    
    bool GetFFDBoxDefinition(void);
    
  bool CheckFFDBoxDefinition(CConfig* config, unsigned short iDV);

    unsigned short GetnFFDBox(void);
    
    unsigned short GetnLevel(void);

  void SetSurface_Derivative(CGeometry *geometry, CConfig *config);
};

#include "grid_movement_structure.inl"