CylinderAndInterfaceDomain Class Reference

Inheritance diagram for CylinderAndInterfaceDomain:

Public Member Functions
	CylinderAndInterfaceDomain (const double &Length, const double &Height)
	~CylinderAndInterfaceDomain ()
void	linear_interpolate (double Left[2], double Right[2], const double &s, Vector< double > &f)
void	macro_element_boundary (const unsigned &time, const unsigned &m, const unsigned &direction, const Vector< double > &s, Vector< double > &f)
Public Member Functions inherited from oomph::Domain
	Domain ()
	Constructor. More...
	Domain (const Domain &)=delete
	Broken copy constructor. More...
void	operator= (const Domain &)=delete
	Broken assignment operator. More...
virtual	~Domain ()
MacroElement *	macro_element_pt (const unsigned &i)
	Access to i-th macro element. More...
unsigned	nmacro_element ()
	Number of macro elements in domain. More...
void	output (const std::string &filename, const unsigned &nplot)
	Output macro elements. More...
void	output (std::ostream &outfile, const unsigned &nplot)
	Output macro elements. More...
virtual void	macro_element_boundary (const double &t, const unsigned &i_macro, const unsigned &i_direct, const Vector< double > &s, Vector< double > &f)
void	macro_element_boundary (const unsigned &i_macro, const unsigned &i_direct, const Vector< double > &s, Vector< double > &f)
void	output_macro_element_boundaries (const std::string &filename, const unsigned &nplot)
	Output all macro element boundaries as tecplot zones. More...
void	output_macro_element_boundaries (std::ostream &outfile, const unsigned &nplot)
	Output all macro element boundaries as tecplot zones. More...
virtual void	dmacro_element_boundary (const unsigned &t, const unsigned &i_macro, const unsigned &i_direct, const Vector< double > &s, Vector< double > &f)
virtual void	dmacro_element_boundary (const double &t, const unsigned &i_macro, const unsigned &i_direct, const Vector< double > &s, Vector< double > &f)
void	dmacro_element_boundary (const unsigned &i_macro, const unsigned &i_direct, const Vector< double > &s, Vector< double > &f)
virtual void	d2macro_element_boundary (const unsigned &t, const unsigned &i_macro, const unsigned &i_direct, const Vector< double > &s, Vector< double > &f)
virtual void	d2macro_element_boundary (const double &t, const unsigned &i_macro, const unsigned &i_direct, const Vector< double > &s, Vector< double > &f)
void	d2macro_element_boundary (const unsigned &i_macro, const unsigned &i_direct, const Vector< double > &s, Vector< double > &f)

Public Attributes
double	centre_x
double	centre_y

Private Attributes
double	Lower_left [2]
double	Lower_right [2]
double	Lower_mid_left [2]
double	Lower_mid_right [2]
double	Upper_left [2]
double	Upper_right [2]
double	Upper_mid_left [2]
double	Upper_mid_right [2]
double	Lower_centre_left [2]
double	Lower_centre_right [2]
double	Upper_centre_left [2]
double	Upper_centre_right [2]
GeomObject *	Cylinder_pt
	Geometric object that represents the rotating cylinder. More...

Additional Inherited Members
Protected Attributes inherited from oomph::Domain
Vector< MacroElement * >	Macro_element_pt
	Vector of pointers to macro elements. More...

Constructor & Destructor Documentation

CylinderAndInterfaceDomain()

CylinderAndInterfaceDomain::CylinderAndInterfaceDomain ( const double &  Length, const double &  Height  )

inline

  {   
 
   centre_x = Length/2.0;
   centre_y = Height/2.0; //3.0*Height/4.0;
  //Create a new ellipse object to represent the rotating cylinder
   Cylinder_pt = new GeneralEllipse(centre_x,centre_y,0.2*Height,0.2*Height);
 
   //Set some basic coordinates
 
   Lower_left[0] = 0.0;
   Lower_left[1] = 0.0;
   
   Upper_left[0] = 0.0;
   Upper_left[1] = Height;
 
   Lower_right[0] = Length;
   Lower_right[1] = 0.0;
   
   Upper_right[0] = Length;
   Upper_right[1] = Height;
 
 
   //Let's just do some mid coordinates
   Lower_mid_left[0] = Length/10.0;
   Lower_mid_left[1] = 0.0;
 
   Upper_mid_left[0] = Length/10.0;
   Upper_mid_left[1] = Height;
 
   Vector<double> xi(1), f(2);
   xi[0] = -3.0*atan(1.0);
   Cylinder_pt->position(xi,f);
 
   Lower_centre_left[0] = f[0];
   Lower_centre_left[1] = f[1];
   
   xi[0] = 3.0*atan(1.0);
   Cylinder_pt->position(xi,f);
 
   Upper_centre_left[0] = f[0];
   Upper_centre_left[1] = f[1];
 
 
   Lower_mid_right[0] = 9.0*Length/10.0;
   Lower_mid_right[1] = 0.0;
 
   Upper_mid_right[0] = 9.0*Length/10.0;
   Upper_mid_right[1] = Height;
 
 
   xi[0] = -1.0*atan(1.0);
   Cylinder_pt->position(xi,f);
 
   Lower_centre_right[0] = f[0];
   Lower_centre_right[1] = f[1];
   
   xi[0] = 1.0*atan(1.0);
   Cylinder_pt->position(xi,f);
 
   Upper_centre_right[0] = f[0];
   Upper_centre_right[1] = f[1];
   
   //There are six macro elements
   Macro_element_pt.resize(6); 
 
   // Build the macro elements
   for (unsigned i=0;i<6;i++)
    {Macro_element_pt[i]= new QMacroElement<2>(this,i);}
  }

References Eigen::bfloat16_impl::atan(), GlobalParameters::Cylinder_pt, f(), GlobalParameters::Height, i, Global_Physical_Variables::Length, and OscillatingCylinder::position().

~CylinderAndInterfaceDomain()

CylinderAndInterfaceDomain::~CylinderAndInterfaceDomain ( )

inline

{}

Member Function Documentation

linear_interpolate()

void CylinderAndInterfaceDomain::linear_interpolate ( double  Left[2], double  Right[2], const double &  s, Vector< double > &  f  )

inline

  {
   for(unsigned i=0;i<2;i++)
    {
     f[i] = Left[i] + (Right[i] - Left[i])*0.5*(s+1.0);
    }
  }

References f(), i, and s.

macro_element_boundary()

void CylinderAndInterfaceDomain::macro_element_boundary ( const unsigned &  t, const unsigned &  i_macro, const unsigned &  i_direct, const Vector< double > &  s, Vector< double > &  f  )

inlinevirtual

Vector representation of the i_macro-th macro element boundary i_direct (e.g. N/S/W/E in 2D) at time level t (t=0: present; t>0: previous): f(s)

Implements oomph::Domain.

 {
 
  using namespace QuadTreeNames;
 
#ifdef WARN_ABOUT_SUBTLY_CHANGED_OOMPH_INTERFACES
   // Warn about time argument being moved to the front
   OomphLibWarning(
    "Order of function arguments has changed between versions 0.8 and 0.85",
    "CylinderAndInterfaceDomain::macro_element_boundary(...)",
    OOMPH_EXCEPTION_LOCATION);
#endif
 
  Vector<double> xi(1);  
 
  //Switch on the macro element
  switch(m)
   {
    //Macro element 0, is the left-hand film
   case 0:
    
    switch(direction)
     {
     case N:
      linear_interpolate(Upper_left,Upper_mid_left,s[0],f);
      break;
 
     case S:  
       linear_interpolate(Lower_left,Lower_mid_left,s[0],f);
      break;
 
     case W:
       linear_interpolate(Lower_left,Upper_left,s[0],f);
      break;
 
     case E:
       linear_interpolate(Lower_mid_left,Upper_mid_left,s[0],f);
      break;
 
     default:
      std::ostringstream error_stream;
      error_stream << "Direction is incorrect: " << direction << std::endl;
 
      throw OomphLibError(error_stream.str(),
                          OOMPH_CURRENT_FUNCTION,
                          OOMPH_EXCEPTION_LOCATION);
     }
    
    break;
    
    //Macro element 1, is immediately left of the cylinder
   case 1:
    
    switch(direction)
     {
     case N:
       linear_interpolate(Upper_mid_left,Upper_centre_left,s[0],f);
      break;
 
     case S:  
       linear_interpolate(Lower_mid_left,Lower_centre_left,s[0],f);
      break;
 
     case W:
       linear_interpolate(Lower_mid_left,Upper_mid_left,s[0],f);
      break;
 
     case E:
      xi[0] = 5.0*atan(1.0) - 2.0*atan(1.0)*0.5*(1.0+s[0]);
      Cylinder_pt->position(xi,f);
      break;
 
     default:
      std::ostringstream error_stream;
      error_stream << "Direction is incorrect: " << direction << std::endl;
      
      throw OomphLibError(error_stream.str(),
                          OOMPH_CURRENT_FUNCTION,
                          OOMPH_EXCEPTION_LOCATION);
     }
    
    break;
 
   //Macro element 2, is immediately above the cylinder
   case 2:
    
    switch(direction)
     {
     case N:
       linear_interpolate(Upper_mid_left,Upper_mid_right,s[0],f);
      break;
      
     case S:  
      xi[0] = 3.0*atan(1.0) - 2.0*atan(1.0)*0.5*(1.0+s[0]);
      Cylinder_pt->position(xi,f);
      break;
 
     case W:
       linear_interpolate(Upper_centre_left,Upper_mid_left,s[0],f);
      break;
 
     case E:
       linear_interpolate(Upper_centre_right,Upper_mid_right,s[0],f);
      break;
 
     default:
      std::ostringstream error_stream;
      error_stream << "Direction is incorrect: " << direction << std::endl;
      
      throw OomphLibError(error_stream.str(),
                          OOMPH_CURRENT_FUNCTION,
                          OOMPH_EXCEPTION_LOCATION);
     }
    
    break;
 
    //Macro element 3, is immediately right of the cylinder
   case 3:
 
    switch(direction)
     {
     case N:
       linear_interpolate(Upper_centre_right,Upper_mid_right,s[0],f);
      break;
      
     case S:  
       linear_interpolate(Lower_centre_right,Lower_mid_right,s[0],f);
      break;
 
     case W:
      xi[0] = -atan(1.0) + 2.0*atan(1.0)*0.5*(1.0+s[0]);
      Cylinder_pt->position(xi,f);
      break;
 
     case E:
       linear_interpolate(Lower_mid_right,Upper_mid_right,s[0],f);
      break;
 
     default:
      std::ostringstream error_stream;
      error_stream << "Direction is incorrect: " << direction << std::endl;
 
      throw OomphLibError(error_stream.str(),
                          OOMPH_CURRENT_FUNCTION,
                          OOMPH_EXCEPTION_LOCATION);
     }
 
    break;
    
    //Macro element 4, is immediately below cylinder
   case 4:
    
    switch(direction)
     {
     case N:
      //linear_interpolate(Lower_centre_left,Lower_centre_right,s[0],f);
      xi[0] = -3.0*atan(1.0) + 2.0*atan(1.0)*0.5*(1.0+s[0]);
      Cylinder_pt->position(xi,f);
      break;
      
     case S:  
       linear_interpolate(Lower_mid_left,Lower_mid_right,s[0],f);
      break;
 
     case W:
       linear_interpolate(Lower_mid_left,Lower_centre_left,s[0],f);
      break;
 
     case E:
       linear_interpolate(Lower_mid_right,Lower_centre_right,s[0],f);
      break;
 
     default:
      std::ostringstream error_stream;
      error_stream << "Direction is incorrect: " << direction << std::endl;
 
      throw OomphLibError(error_stream.str(),
                          OOMPH_CURRENT_FUNCTION,
                          OOMPH_EXCEPTION_LOCATION);
     }  
 
    break;
 
    //Macro element 5, is right film
   case 5:
    
    switch(direction)
     {
     case N:
       linear_interpolate(Upper_mid_right,Upper_right,s[0],f);
      break;
      
     case S:  
       linear_interpolate(Lower_mid_right,Lower_right,s[0],f);
      break;
 
     case W:
       linear_interpolate(Lower_mid_right,Upper_mid_right,s[0],f);
      break;
 
     case E:
       linear_interpolate(Lower_right,Upper_right,s[0],f);
      break;
 
     default:
      std::ostringstream error_stream;
      error_stream << "Direction is incorrect: " << direction << std::endl;
 
      throw OomphLibError(error_stream.str(),
                          OOMPH_CURRENT_FUNCTION,
                          OOMPH_EXCEPTION_LOCATION);
     }
 
    break;
    
   default:
    std::ostringstream error_stream;
      error_stream << "Wrong domain number: " << m<< std::endl;
 
      throw OomphLibError(error_stream.str(),
                          OOMPH_CURRENT_FUNCTION,
                          OOMPH_EXCEPTION_LOCATION);
   }
 
 }

References Eigen::bfloat16_impl::atan(), GlobalParameters::Cylinder_pt, Global_Physical_Variables::E, f(), m, N, OOMPH_CURRENT_FUNCTION, OOMPH_EXCEPTION_LOCATION, OscillatingCylinder::position(), s, oomph::QuadTreeNames::S, and oomph::QuadTreeNames::W.

Member Data Documentation

centre_x

double CylinderAndInterfaceDomain::centre_x

centre_y

double CylinderAndInterfaceDomain::centre_y

Cylinder_pt

GeomObject* CylinderAndInterfaceDomain::Cylinder_pt

private

Geometric object that represents the rotating cylinder.

Lower_centre_left

double CylinderAndInterfaceDomain::Lower_centre_left[2]

private

Lower_centre_right

double CylinderAndInterfaceDomain::Lower_centre_right[2]

private

Lower_left

double CylinderAndInterfaceDomain::Lower_left[2]

private

Lower_mid_left

double CylinderAndInterfaceDomain::Lower_mid_left[2]

private

Lower_mid_right

double CylinderAndInterfaceDomain::Lower_mid_right[2]

private

Lower_right

double CylinderAndInterfaceDomain::Lower_right[2]

private

Upper_centre_left

double CylinderAndInterfaceDomain::Upper_centre_left[2]

private

Upper_centre_right

double CylinderAndInterfaceDomain::Upper_centre_right[2]

private

Upper_left

double CylinderAndInterfaceDomain::Upper_left[2]

private

Upper_mid_left

double CylinderAndInterfaceDomain::Upper_mid_left[2]

private

Upper_mid_right

double CylinderAndInterfaceDomain::Upper_mid_right[2]

private

Upper_right

double CylinderAndInterfaceDomain::Upper_right[2]

private

---

The documentation for this class was generated from the following file:

• adaptive_interface.cc