120 #include "et_rot_mag.h" 121 #include "utilitaires.h" 211 if ( withbphi == 0 ) {
346 ost <<
"Electromagnetic quantities" << endl ;
347 ost <<
"----------------------" << endl ;
350 ost <<
"***************************" << endl ;
351 ost <<
"**** perfect conductor ****" << endl ;
352 ost <<
"***************************" << endl ;
353 ost <<
"Prescribed charge : " <<
Q*j_unit*
pow(r_unit,3)/v_unit
357 ost <<
"***************************" << endl ;
358 ost <<
"**** isolator ****" << endl ;
359 ost <<
"***************************" << endl ;
361 ost <<
"Prescribed current amplitude : " <<
a_j*j_unit
362 <<
" [A/m2]" << endl ;
363 ost <<
"Magnetic Momentum : " <<
MagMom()/1.e32
364 <<
" [10^32 Am^2]" << endl ;
365 ost <<
"Radial magnetic field polar value : " <<
369 ost <<
"Tangent magnetic field equatorial value : " <<
374 ost <<
"Central magnetic field values : " <<
375 Magn()(0)(0,0,0,0)*mag_unit/1.e9
378 ost <<
"Radial electric field polar value : " <<
383 ost <<
"Radial electric field equatorial value : " <<
388 ost <<
"Magnetic/fluid pressure : " 389 << 1/(2*mu_si)*(
pow(
Magn()(0)(0,0,0,0),2) +
391 pow(
Magn()(2)(0,0,0,0),2))*mag_unit*mag_unit
392 / (
press()(0,0,0,0)*rho_unit*
pow(v_unit,2)) << endl ;
393 ost <<
"Computed charge : " <<
Q_comput() <<
" [C]" << endl ;
394 ost <<
"Interior charge : " <<
Q_int() <<
" [C]" << endl ;
397 ost <<
"Gyromagnetic ratio : " <<
GyroMag() << endl ;
const Tbl & xi_surf() const
Description of the stellar surface: returns a 2-D Tbl containing the values of the radial coordinate...
Tenseur E_em
electromagnetic energy density in the Eulerian frame
Component of a tensorial field *** DEPRECATED : use class Scalar instead ***.
double GyroMag() const
Gyromagnetic ratio .
double Q
In the case of a perfect conductor, the requated baryonic charge.
Equation of state base class.
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
Base class for coordinate mappings.
Class for isolated rotating stars *** DEPRECATED : use class Star_rot instead ***.
virtual double mass_g() const
Gravitational mass.
Tenseur Magn() const
Computes the magnetic field spherical components in Lorene's units.
void operator=(const Etoile_rot &)
Assignment to another Etoile_rot.
virtual ~Et_rot_mag()
Destructor.
Tenseur press
Fluid pressure.
virtual const Itbl & l_surf() const
Description of the stellar surface: returns a 2-D Itbl containing the values of the domain index l on...
virtual void del_hydro_euler()
Sets to ETATNONDEF (undefined state) the hydrodynamical quantities relative to the Eulerian observer...
Cmp A_phi
-component of the electromagnetic potential 1-form divided by .
virtual void del_deriv() const
Deletes all the derived quantities.
Tenseur Srr_em
rr component of the electromagnetic stress 3-tensor, as measured in the Eulerian frame. (not used and set to 0, should be supressed)
virtual ostream & operator>>(ostream &) const
Operator >> (virtual function called by the operator <<).
void operator=(const Et_rot_mag &)
Assignment to another Et_rot_mag.
double a_j
Amplitude of the curent/charge function.
void sauve(FILE *) const
Save in a file.
bool is_conduct() const
Tells if the star is made of conducting or isolating material.
Class for magnetized (isolator or perfect conductor), rigidly rotating stars.
virtual ostream & operator>>(ostream &) const
Operator >> (virtual function called by the operator <<).
virtual void sauve(FILE *) const
Save in a file.
Map & mp
Mapping associated with the star.
Et_rot_mag(Map &mp_i, int nzet_i, bool relat, const Eos &eos_i, const int cond)
Standard constructor.
virtual void del_deriv() const
Deletes all the derived quantities.
int fwrite_be(const int *aa, int size, int nb, FILE *fich)
Writes integer(s) into a binary file according to the big endian convention.
Cmp pow(const Cmp &, int)
Power .
Cmp j_phi
-component of the current 4-vector
int nzet
Number of domains of *mp occupied by the star.
int fread_be(int *aa, int size, int nb, FILE *fich)
Reads integer(s) from a binary file according to the big endian convention.
Cmp j_t
t-component of the current 4-vector
double Q_int() const
Computed charge from the integration of charge density over the star (i.e.
virtual void sauve(FILE *) const
Save in a file.
virtual void set_der_0x0() const
Sets to 0x0 all the pointers on derived quantities.
Tenseur Elec() const
Computes the electric field spherical components in Lorene's units.
Cmp A_t
t-component of the elecctromagnetic potential 1-form, divided by .
virtual void set_der_0x0() const
Sets to 0x0 all the pointers on derived quantities.
int get_nt(int l) const
Returns the number of points in the co-latitude direction ( ) in domain no. l.
int conduc
Flag: conduc=0->isolator, 1->perfect conductor.
void sauve(FILE *) const
Save in a file.
double Q_comput() const
Computed charge deduced from the asymptotic behaviour of At [SI units].
virtual void del_hydro_euler()
Sets to ETATNONDEF (undefined state) the hydrodynamical quantities relative to the Eulerian observer...
double MagMom() const
Magnetic Momentum in SI units.
Cmp B_phi
-component of the magnetic field
Standard electro-magnetic units.
Tensor handling *** DEPRECATED : use class Tensor instead ***.
Tenseur Jp_em
component of the electromagnetic momentum density 3-vector, as measured in the Eulerian frame...
Tenseur Spp_em
component of the electromagnetic stress 3-tensor, as measured in the Eulerian frame.