Source giving the acceleration due to gravity in the spherical, Newtonian Lane-Emden star solution. More...

#include <LaneEmdenGravitationalField.hpp>

Public Types
using	sourced_variables = tmpl::list< Tags::MomentumDensity< 3 >, Tags::EnergyDensity >

using	argument_tags = tmpl::list< Tags::MassDensityCons, Tags::MomentumDensity< 3 >, ::Tags::AnalyticSolution< NewtonianEuler::Solutions::LaneEmdenStar >, domain::Tags::Coordinates< 3, Frame::Inertial > >

Public Member Functions
	LaneEmdenGravitationalField (const LaneEmdenGravitationalField &)=default

LaneEmdenGravitationalField &	operator= (const LaneEmdenGravitationalField &)=default

	LaneEmdenGravitationalField (LaneEmdenGravitationalField &&) noexcept=default

LaneEmdenGravitationalField &	operator= (LaneEmdenGravitationalField &&) noexcept=default

void	pup (PUP::er &) noexcept

Static Public Member Functions
static void	apply (gsl::not_null< tnsr::I< DataVector, 3 > * > source_momentum_density, gsl::not_null< Scalar< DataVector > * > source_energy_density, const Scalar< DataVector > &mass_density_cons, const tnsr::I< DataVector, 3 > &momentum_density, const NewtonianEuler::Solutions::LaneEmdenStar &star, const tnsr::I< DataVector, 3 > &x) noexcept

Detailed Description

Source giving the acceleration due to gravity in the spherical, Newtonian Lane-Emden star solution.

The gravitational field \(g^i\) enters the NewtonianEuler system as source terms for the conserved momentum and energy:

\begin{align*} \partial_t S^i + \partial_j F^{j}(S^i) &= S(S^i) = \rho g^i \partial_t e + \partial_j F^{j}(e) &= S(e) = S_i g^i, \end{align*}

where \(S^i\) is the conserved momentum density, \(e\) is the conserved energy, \(F^{j}(u)\) is the flux of the conserved quantity \(u\), and \(\rho\) is the fluid mass density.

Note: This source is specialized to the Lane-Emden solution because it queries a LaneEmdenStar analytic solution for the gravitational field that generates the fluid acceleration. This source does not integrate the fluid density to compute a self-consistent gravitational field (i.e., as if one were solving a coupled Euler + Poisson system).

Public Types

Public Member Functions

Static Public Member Functions

Detailed Description