hot_strot_diff_cfc_hydro.C

/*
 * Function Hot_star_rot_diff_CFC::hydro_euler
 *
 * (see file hot_star_rot_diff_cfc.h for documentation)
 *
 */

/*
 *   Copyright (c) 2025 Santiago Jaraba
 *
 *   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
 *
 */


// C headers
#include <cstdlib>

// Lorene headers
#include "hot_star_rot_diff_cfc.h"
#include "utilitaires.h"

namespace Lorene {
void Hot_star_rot_diff_CFC::hydro_euler(){
  
    // u_euler (fluid 3-velocity w.r.t. the Eulerian frame)
    // -----------------------------------------------------
    

    u_euler.set(1).set_etat_zero() ;
    u_euler.set(2).set_etat_zero() ;
    
    u_euler.set(3) = omega_field ;
    u_euler.set(3).std_spectral_base() ;
    u_euler.set(3).mult_rsint() ;
    u_euler.set(3) += beta(3) ;
    
    u_euler = u_euler / nn ;

    // v2 (square of the norm of u_euler)
    // ----------------------------------
    
    v2 = contract(contract(gamma.cov(), 0, u_euler, 0), 0, u_euler, 0) ;
    
    
    // gam_euler (Lorentz factor between the fluid and Eulerian observers)
    // -------------------------------------------------------------------
    
    gam_euler = 1. / sqrt(1. - v2) ; 
    
    gam_euler.std_spectral_base() ;
    
    
    // ener_euler (energy density w.r.t. the Eulerian observer)
    // ------------------------------------------------------
    
    ener_euler = gam_euler*gam_euler*(ener + press) - press ;
    
    ener_euler.std_spectral_base() ;
    if (spectral_filter > 0) {
      ener_euler.exponential_filter_r(1, nzet-1, spectral_filter) ;
    }
    // j_euler (momentum density 3-vector w.r.t. the Eulerian observer)
    // ----------------------------------------------------------------
    
    j_euler = (ener_euler + press)*u_euler ;
    j_euler.std_spectral_base() ;
    
    if (spectral_filter > 0) {
      j_euler.exponential_filter_r(1, nzet-1, spectral_filter) ;
    }
    
    // s_euler (trace of the stress tensor w.r.t. the Eulerian observer)
    // ----------------------------------------------------------------
    
    s_euler = (ener_euler + press)*v2 + 3*press ;
    s_euler.std_spectral_base() ;
    if (spectral_filter > 0) {
      s_euler.exponential_filter_r(1, nzet-1, spectral_filter) ;
    }
    // stress_euler (stress tensor w.r.t. the Eulerian observer)
    // ---------------------------------------------------------
    
    
    stress_euler = (ener_euler + press)*u_euler*u_euler + press*gamma.con() ;
    stress_euler.std_spectral_base() ;
    if (spectral_filter > 0) {
      stress_euler.exponential_filter_r(1, nzet-1, spectral_filter) ;
    }
    
    // The derived quantities are obsolete
    // -----------------------------------
    
    del_deriv() ;                

}
}