linearised_navier_stokes_eigenvalue_elements.cc

Go to the documentation of this file.

//LIC// ====================================================================
//LIC// This file forms part of oomph-lib, the object-oriented, 
//LIC// multi-physics finite-element library, available 
//LIC// at http://www.oomph-lib.org.
//LIC// 
//LIC//    Version 1.0; svn revision $LastChangedRevision: 1097 $
//LIC//
//LIC// $LastChangedDate: 2015-12-17 11:53:17 +0000 (Thu, 17 Dec 2015) $
//LIC// 
//LIC// Copyright (C) 2006-2016 Matthias Heil and Andrew Hazel
//LIC// 
//LIC// This library is free software; you can redistribute it and/or
//LIC// modify it under the terms of the GNU Lesser General Public
//LIC// License as published by the Free Software Foundation; either
//LIC// version 2.1 of the License, or (at your option) any later version.
//LIC// 
//LIC// This library is distributed in the hope that it will be useful,
//LIC// but WITHOUT ANY WARRANTY; without even the implied warranty of
//LIC// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
//LIC// Lesser General Public License for more details.
//LIC// 
//LIC// You should have received a copy of the GNU Lesser General Public
//LIC// License along with this library; if not, write to the Free Software
//LIC// Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA
//LIC// 02110-1301  USA.
//LIC// 
//LIC// The authors may be contacted at oomph-lib@maths.man.ac.uk.
//LIC// 
//LIC//====================================================================
// The element-independent guts for imposition of "constant volume"
// constraints in free surface/interface problems.


#include "linearised_navier_stokes_eigenvalue_elements.h"

namespace oomph
{
 
//=====================================================================
/// \short Fill in the residuals for the volume constraint
//====================================================================
 void LinearisedNavierStokesEigenfunctionNormalisationElement::
 fill_in_generic_contribution_to_residuals_normalisation(
  Vector<double> &residuals)
 {
  //Note that this element can only be used with our linearised navier
  //stokes elements
  //Read part
  {
   const int local_eqn = this->eigenvalue_local_eqn(0);
   if(local_eqn >= 0)
    {
     residuals[local_eqn] -= (*Normalisation_pt).real();
    }
  }
  //Imaginary part
  {
   const int local_eqn = this->eigenvalue_local_eqn(1);
   if(local_eqn >= 0)
    {
     residuals[local_eqn] -= (*Normalisation_pt).imag();
    }
  }
  //Bifurcation constraint
  const int local_eqn = this->eigenvalue_local_eqn(2);
  if(local_eqn >= 0)
   {
    residuals[local_eqn] += this->eigenvalue(0);
   }
  
 }
 
//===========================================================================
/// \short Constructor: Pass pointer to target volume. "Pressure" value that
/// "traded" for the volume contraint is created internally (as a Data
/// item with a single pressure value)
//===========================================================================
 LinearisedNavierStokesEigenfunctionNormalisationElement::
 LinearisedNavierStokesEigenfunctionNormalisationElement(
  std::complex<double>* const &normalisation_pt)
 {
  // Store pointer to normalisation
  Normalisation_pt = normalisation_pt;
  
  // Create data, add as internal data and record the index
  // (gets deleted automatically in destructor of GeneralisedElement)
  External_or_internal_data_index_of_eigenvalue=
   add_internal_data(new Data(3)); 
  
  // ...and stored the "traded pressure" value as first value
  Index_of_eigenvalue=0; 
 } 
  
//======================================================================
/// \short Constructor: Pass pointer to target volume, pointer to Data
/// item whose value specified by index_of_traded_pressure represents
/// the "Pressure" value that "traded" for the volume contraint.
/// The Data is stored as external Data for this element.
//======================================================================
/* LinearisedNavierStokesEigenfunctionNormalisationElement::LinearisedNavierStokesEigenfunctionNormalisationElement(
  double* prescribed_volume_pt,
  Data* p_traded_data_pt,
  const unsigned& index_of_traded_pressure)
 {
  // Store pointer to prescribed volume
  Prescribed_volume_pt = prescribed_volume_pt;
  
  // Add as external data and record the index
  External_or_internal_data_index_of_traded_pressure=
   add_external_data(p_traded_data_pt);
  
  // Record that it is external data
  Traded_pressure_stored_as_internal_data=false;
  
  // Record index
  Index_of_traded_pressure_value=index_of_traded_pressure;
  } */



}