Shamrock 2025.10.0
Astrophysical Code
Loading...
Searching...
No Matches
Namespaces | Functions
forces.hpp File Reference

file containing formulas for sph forces More...

#include "shambase/numeric_limits.hpp"
#include "shambackends/math.hpp"
+ Include dependency graph for forces.hpp:
+ This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Namespaces

namespace  shamrock
 namespace for the main framework
 

Functions

template<class Tvec , class Tscal >
Tvec shamrock::sph::sph_pressure_symetric (const Tscal &m_b, const Tscal &rho_a_sq, const Tscal &rho_b_sq, const Tscal &P_a, const Tscal &P_b, const Tscal &omega_a, const Tscal &omega_b, const Tvec &nabla_Wab_ha, const Tvec &nabla_Wab_hb)
 [2] eq.34, with \(q^a_{ab} = q^b_{ab} = 0\)
 
template<class Tvec , class Tscal >
Tvec shamrock::sph::sph_pressure_symetric_av (const Tscal &m_b, const Tscal &rho_a_sq, const Tscal &rho_b_sq, const Tscal &P_a, const Tscal &P_b, const Tscal &omega_a, const Tscal &omega_b, const Tscal &qa_ab, const Tscal &qb_ab, const Tvec &nabla_Wab_ha, const Tvec &nabla_Wab_hb)
 [2] eq.34
 
template<class Tvec , class Tscal >
Tscal shamrock::sph::duint_dt_pressure (const Tscal &pmass, const Tscal &P_a, const Tscal &inv_omega_a_2_rho_a, const Tvec &v_ab, const Tvec &grad_W_ab)
 [2] eq.35
 
template<class Tscal >
Tscal shamrock::sph::lambda_shock_conductivity (const Tscal &pmass, const Tscal &alpha_u, const Tscal &vsig_u, const Tscal &u_ab, const Tscal &Fab_inv_omega_a_rho_a, const Tscal &Fab_inv_omega_b_rho_b)
 [2] eq.42
 
template<class Tvec , class Tscal >
void shamrock::sph::add_to_derivs_sph_artif_visco_cond (const Tscal &pmass, const Tscal &rho_a_sq, const Tscal &omega_a_rho_a_inv, const Tscal &rho_a_inv, const Tscal &rho_b, const Tscal &omega_a, const Tscal &omega_b, const Tscal &Fab_a, const Tscal &Fab_b, const Tscal &u_a, const Tscal &u_b, const Tscal &P_a, const Tscal &P_b, const Tscal &alpha_u, const Tvec &v_ab, const Tvec &r_ab_unit, const Tscal &vsig_u, const Tscal &qa_ab, const Tscal &qb_ab, Tvec &dv_dt, Tscal &du_dt)
 

Detailed Description

file containing formulas for sph forces

Author
Timothée David–Cléris (tim.s.nosp@m.hamr.nosp@m.ock@p.nosp@m.roto.nosp@m.n.me)
Yona Lapeyre (yona..nosp@m.lape.nosp@m.yre@e.nosp@m.ns-l.nosp@m.yon.f.nosp@m.r)

Definition in file forces.hpp.

Function Documentation

◆ add_to_derivs_sph_artif_visco_cond()

template<class Tvec , class Tscal >
void shamrock::sph::add_to_derivs_sph_artif_visco_cond ( const Tscal &  pmass,
const Tscal &  rho_a_sq,
const Tscal &  omega_a_rho_a_inv,
const Tscal &  rho_a_inv,
const Tscal &  rho_b,
const Tscal &  omega_a,
const Tscal &  omega_b,
const Tscal &  Fab_a,
const Tscal &  Fab_b,
const Tscal &  u_a,
const Tscal &  u_b,
const Tscal &  P_a,
const Tscal &  P_b,
const Tscal &  alpha_u,
const Tvec &  v_ab,
const Tvec &  r_ab_unit,
const Tscal &  vsig_u,
const Tscal &  qa_ab,
const Tscal &  qb_ab,
Tvec &  dv_dt,
Tscal &  du_dt 
)
inline

Definition at line 173 of file forces.hpp.

◆ duint_dt_pressure()

template<class Tvec , class Tscal >
Tscal shamrock::sph::duint_dt_pressure ( const Tscal &  pmass,
const Tscal &  P_a,
const Tscal &  inv_omega_a_2_rho_a,
const Tvec &  v_ab,
const Tvec &  grad_W_ab 
)
inline

[2] eq.35

Template Parameters
Tvec
Tscal
Parameters
P_a
omega_a_rho_a_inv
rho_a_inv
pmass
v_ab
grad_W_ab
Returns
Tscal

Definition at line 138 of file forces.hpp.

+ Here is the call graph for this function:

◆ lambda_shock_conductivity()

template<class Tscal >
Tscal shamrock::sph::lambda_shock_conductivity ( const Tscal &  pmass,
const Tscal &  alpha_u,
const Tscal &  vsig_u,
const Tscal &  u_ab,
const Tscal &  Fab_inv_omega_a_rho_a,
const Tscal &  Fab_inv_omega_b_rho_b 
)
inline

[2] eq.42

Template Parameters
Tvec
Tscal
Parameters
pmass
alpha_u
vsig_u
u_abdefined as : \(u_a - u_b\)
Fab_inv_omega_a_rho_a
Fab_inv_omega_b_rho_b
Returns
Tscal

Definition at line 161 of file forces.hpp.

+ Here is the call graph for this function:

◆ sph_pressure_symetric()

template<class Tvec , class Tscal >
Tvec shamrock::sph::sph_pressure_symetric ( const Tscal &  m_b,
const Tscal &  rho_a_sq,
const Tscal &  rho_b_sq,
const Tscal &  P_a,
const Tscal &  P_b,
const Tscal &  omega_a,
const Tscal &  omega_b,
const Tvec &  nabla_Wab_ha,
const Tvec &  nabla_Wab_hb 
)
inline

[2] eq.34, with \(q^a_{ab} = q^b_{ab} = 0\)

Template Parameters
Tvec
Tscal
Parameters
m_b
rho_a_sq
rho_b_sq
P_a
P_b
omega_a
omega_b
nabla_Wab_ha
nabla_Wab_hb
Returns
Tvec

Definition at line 60 of file forces.hpp.

+ Here is the call graph for this function:

◆ sph_pressure_symetric_av()

template<class Tvec , class Tscal >
Tvec shamrock::sph::sph_pressure_symetric_av ( const Tscal &  m_b,
const Tscal &  rho_a_sq,
const Tscal &  rho_b_sq,
const Tscal &  P_a,
const Tscal &  P_b,
const Tscal &  omega_a,
const Tscal &  omega_b,
const Tscal &  qa_ab,
const Tscal &  qb_ab,
const Tvec &  nabla_Wab_ha,
const Tvec &  nabla_Wab_hb 
)
inline

[2] eq.34

Template Parameters
Tvec
Tscal
Parameters
m_b
rho_a_sq
rho_b_sq
P_a
P_b
omega_a
omega_b
qa_ab
qb_ab
nabla_Wab_ha
nabla_Wab_hb
Returns
Tvec

Definition at line 100 of file forces.hpp.

+ Here is the call graph for this function: