binaire_ana_shift.C
/*
* Method of class Binaire to set some analytical form to the shift vector.
*
* (see file binaire.h for documentation).
*/
/*
* Copyright (c) 2000-2001 Eric Gourgoulhon
* Copyright (c) 2000-2001 Keisuke Taniguchi
*
* This file is part of LORENE.
*
* LORENE is free software; you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by
* the Free Software Foundation; either version 2 of the License, or
* (at your option) any later version.
*
* LORENE is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with LORENE; if not, write to the Free Software
* Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
*
*/
char binaire_ana_shift_C[] = "$Header: /cvsroot/Lorene/C++/Source/Binaire/binaire_ana_shift.C,v 1.3 2014/10/13 08:52:44 j_novak Exp $" ;
/*
* $Id: binaire_ana_shift.C,v 1.3 2014/10/13 08:52:44 j_novak Exp $
* $Log: binaire_ana_shift.C,v $
* Revision 1.3 2014/10/13 08:52:44 j_novak
* Lorene classes and functions now belong to the namespace Lorene.
*
* Revision 1.2 2004/03/25 10:28:59 j_novak
* All LORENE's units are now defined in the namespace Unites (in file unites.h).
*
* Revision 1.1.1.1 2001/11/20 15:19:30 e_gourgoulhon
* LORENE
*
* Revision 2.3 2000/03/17 15:36:26 eric
* Suppression de l'appel a analytical_omega().
*
* Revision 2.2 2000/03/17 15:27:11 eric
* Appel de la fonction analytical_omega() pour fixer la valeur de omega.
*
* Revision 2.1 2000/03/16 09:37:28 eric
* Utilisation du cas incompressible plutot que n=1.
*
* Revision 2.0 2000/03/15 16:43:35 eric
* *** empty log message ***
*
*
* $Header: /cvsroot/Lorene/C++/Source/Binaire/binaire_ana_shift.C,v 1.3 2014/10/13 08:52:44 j_novak Exp $
*
*/
// Headers C
#include "math.h"
// Headers Lorene
#include "binaire.h"
#include "unites.h"
namespace Lorene {
void Binaire::analytical_shift(){
// Does nothing for a Newtonian star
// ---------------------------------
if ( !star1.is_relativistic() ){
assert( !star2.is_relativistic() ) ;
return ;
}
using namespace Unites ;
for (int i=0; i<2; i++) {
// Radius of the star:
double a0 = et[i]->ray_eq() ;
// Mass ratio
double p_mass = et[i]->mass_g() / et[1-i]->mass_g() ;
// G M Omega R / (1+p)
double www = ggrav * et[i]->mass_g() * omega
* separation() / (1. + p_mass) ;
const Map& mp = et[i]->get_mp() ;
Cmp tmp(mp) ;
Cmp tmp_ext(mp) ;
int nzet = et[i]->get_nzet() ;
int nzm1 = mp.get_mg()->get_nzone() - 1 ;
// Computation of w_shift
// ----------------------
et[i]->set_w_shift().set_etat_qcq() ;
// X component
// -----------
et[i]->set_w_shift().set(0) = 0 ;
// Y component
// -----------
// For the incompressible case :
tmp = - 6 * www / a0 * ( 1 - (mp.r)*(mp.r) / (3*a0*a0) ) ;
// For the compressible (n=1) case :
// Mtbl xi = M_PI * mp.r / a0 ;
// Mtbl sinc = sin(xi) / xi ;
// The value of sinc is set to 1 at the origin
// for (int k=0; k<mp.get_mg()->get_np(0); k++) {
// for (int j=0; j<mp.get_mg()->get_nt(0); j++) {
// sinc.set(0, k, j, 0) = 1 ;
// }
// }
// tmp = - 4 * www / a0 * ( 1 + sinc ) ;
tmp.annule(nzet, nzm1) ;
tmp_ext = - 4 * www / mp.r ;
tmp_ext.annule(0, nzet-1) ;
et[i]->set_w_shift().set(1) = tmp + tmp_ext ;
// Z component
// -----------
et[i]->set_w_shift().set(2) = 0 ;
// Sets the standard spectral bases for Cartesian components
et[i]->set_w_shift().set_std_base() ;
// Computation of khi_shift
// ------------------------
tmp = 2 * www / a0 * (mp.y) * ( 1 - 3 * (mp.r)*(mp.r) / (5*a0*a0) ) ;
tmp.annule(nzet, nzm1) ;
tmp_ext = 0.8 * www * a0*a0 * (mp.sint) * (mp.sinp)
/ (mp.r * mp.r) ;
tmp_ext.annule(0, nzet-1) ;
et[i]->set_khi_shift() = tmp + tmp_ext ;
// Sets the standard spectral bases for a scalar field
et[i]->set_khi_shift().set_std_base() ;
}
}
}
Lorene::Binaire::analytical_shift
void analytical_shift()
Sets some analytical template for the shift vector (via the members {\tt w_shift} and {\tt khi_shift}...
Definition binaire_ana_shift.C:69
Lorene::Binaire::et
Etoile_bin * et[2]
Array of the two stars (to perform loops on the stars): {\tt et[0]} contains the address of {\tt star...
Definition binaire.h:124
Lorene::Binaire::separation
double separation() const
Returns the coordinate separation of the two stellar centers [{\tt r_unit}].
Definition binaire.C:457
Lorene::Binaire::star2
Etoile_bin star2
Second star of the system.
Definition binaire.h:118
Lorene::Binaire::omega
double omega
Angular velocity with respect to an asymptotically inertial observer.
Definition binaire.h:129
Lorene::Binaire::star1
Etoile_bin star1
First star of the system.
Definition binaire.h:115
Lorene::Cmp
Component of a tensorial field *** DEPRECATED : use class Scalar instead ***.
Definition cmp.h:446
Lorene::Cmp::annule
void annule(int l)
Sets the Cmp to zero in a given domain.
Definition cmp.C:348
Lorene
Lorene prototypes.
Definition app_hor.h:64
Lorene::Map
Map(const Mg3d &)
Constructor from a multi-domain 3D grid.
Definition map.C:139
Unites
Standard units of space, time and mass.