28char star_bhns_hydro_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Star_bhns/star_bhns_hydro.C,v 1.3 2014/10/13 08:53:41 j_novak Exp $" ;
55#include "utilitaires.h"
78 double mass = ggrav * mass_bh ;
82 xx.std_spectral_base() ;
85 yy.std_spectral_base() ;
88 zz.std_spectral_base() ;
94 rbh.std_spectral_base() ;
105 msr.std_spectral_base() ;
119 dpsi2.std_spectral_base() ;
123 lldpsi.std_spectral_base() ;
131 tmp1.std_spectral_base() ;
139 for (
int i=1;
i<=3;
i++)
199 llbsn.std_spectral_base() ;
203 tmp2.std_spectral_base() ;
260 dpsi2.std_spectral_base() ;
266 for (
int i=1;
i<=3;
i++)
321 tmp3.std_spectral_base() ;
Time evolution with partial storage (*** under development ***).
const Base_vect_cart & get_bvect_cart() const
Returns the Cartesian basis associated with the coordinates (x,y,z) of the mapping,...
Coord y
y coordinate centered on the grid
double get_ori_y() const
Returns the y coordinate of the origin.
Coord z
z coordinate centered on the grid
Coord x
x coordinate centered on the grid
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
Tensor field of valence 0 (or component of a tensorial field).
virtual void std_spectral_base()
Sets the spectral bases of the Valeur va to the standard ones for a scalar field.
void set_dzpuis(int)
Modifies the dzpuis flag.
Vector wit_w
Spatial projection of the fluid 3-velocity with respect to the co-orbiting observer.
Vector bsn
3-vector shift, divided by N , of the rotating coordinates, .
virtual void del_deriv() const
Deletes all the derived quantities.
Scalar psi4
Fourth power of the total conformal factor.
Vector d_psi
Gradient of (in the irrotational case) (Spherical components with respect to the mapping of the star...
Scalar gam0
Lorentz factor between the co-orbiting observer and the Eulerian one.
void hydro_euler_bhns(bool kerrschild, const double &mass_bh, const double &sepa)
Computes the hydrodynamical quantities relative to the Eulerian observer from those in the fluid fram...
Scalar gam
Lorentz factor between the fluid and the co-orbiting observer.
bool irrotational
true for an irrotational star, false for a corotating one
Scalar loggam
Logarithm of the Lorentz factor between the fluid and the co-orbiting observer.
Scalar ener
Total energy density in the fluid frame.
Scalar ener_euler
Total energy density in the Eulerian frame.
Scalar gam_euler
Lorentz factor between the fluid and Eulerian observers.
Scalar s_euler
Trace of the stress scalar in the Eulerian frame.
Scalar press
Fluid pressure.
Vector u_euler
Fluid 3-velocity with respect to the Eulerian observer.
Map & mp
Mapping associated with the star.
Tensor field of valence 1.
virtual void std_spectral_base()
Sets the standard spectal bases of decomposition for each component.
Scalar & set(int)
Read/write access to a component.
Cmp sqrt(const Cmp &)
Square root.
Cmp exp(const Cmp &)
Exponential.
Cmp log(const Cmp &)
Neperian logarithm.
const Base_vect * get_triad() const
Returns the vectorial basis (triad) on which the components are defined.
void annule_domain(int l)
Sets the Tensor to zero in a given domain.
virtual void set_etat_qcq()
Sets the logical state of all components to ETATQCQ (ordinary state).
virtual void set_etat_zero()
Sets the logical state of all components to ETATZERO (zero state).
Standard units of space, time and mass.