constvelocitymodule.cpp

Go to the documentation of this file.

/*
 * constvelocitymodule.cpp
 *
 *  Created on: Sep 15, 2008
 *      Author: marcosm
 */
#include "constvelocitymodule.h"

ConstVelocityModule::ConstVelocityModule(OrthoMesh &mesh, Function3D & f)
  :m_mesh(mesh),m_velFunction(f),m_vPressure(m_mesh.numCells())
{
  assert(f.n_components() == 3);
  m_vPressure=0.0;
}

ConstVelocityModule::ConstVelocityModule(OrthoMesh &mesh, Function3D & f,Function3D &fPressure)
  :m_mesh(mesh),m_velFunction(f),m_vPressure(m_mesh.numCells())
{
  assert(f.n_components() == 3);
  mesh.setCentralValuesFromFunction(fPressure,m_vPressure);
}




ConstVelocityModule::~ConstVelocityModule()
{


}

void ConstVelocityModule::getVelocitiesAtFaces(Matrix & vel)
{
  m_mesh.setFacesValuesFromFunction(m_velFunction,vel);
}




void ConstVelocityModule::iterate(TransportBase &trans)
{
  //This is a dummy dynamic module just to test the transport.
  //So, dont do anything.
  
  
  return;

}








 void ConstVelocityModule::getNormalVelocityAtFaces(VecDouble &vNormal) 
{
  assert(vNormal.size()==m_mesh.numFaces());
  OrthoMesh::Face_It face    = m_mesh.begin_face();
  OrthoMesh::Face_It endFace = m_mesh.end_face();
  
  for (;face!=endFace;face++)
  {

    Point<3> P = face->barycenter();
    unsigned faceIndex = face->index();
    for (unsigned deg=0;deg<m_velFunction.n_components();deg++)
    {
      vNormal(faceIndex) = m_velFunction(P,face->getNormalOrientation());
    }
  }
}


 const VecDouble& ConstVelocityModule::getPressureAtCells() 
{
  return m_vPressure;
}


 void ConstVelocityModule::setDt(double dt) 
{

}

---