reconstruction.hpp File Reference

Effective face (volume element) interpolation for the Inutsuka (2002) GSPH formulation. More...

#include "shambackends/math.hpp"

Include dependency graph for reconstruction.hpp:

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Classes
struct	shammodels::gsph::EffectiveFace< Tscal >
	Result of the effective face interpolation. More...

Namespaces
namespace	shammodels
	namespace for models

Functions
template<class Tscal>
EffectiveFace< Tscal >	shammodels::gsph::lin_v2_sast_ij (Tscal vol_a, Tscal vol_b, Tscal h_a, Tscal h_b, Tscal rab_inv)
	Linear (1st order) interpolation of the effective face between a pair.

Detailed Description

Effective face (volume element) interpolation for the Inutsuka (2002) GSPH formulation.

Author

Guo Yansong (guo.y.nosp@m.anso.nosp@m.ng.ng.nosp@m.y@gm.nosp@m.ail.c.nosp@m.om)

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)

Inutsuka's GSPH momentum equation replaces the standard SPH kernel-gradient/rho^2 weighting by an effective face between each particle pair, described by an effective squared volume element V2_ij and an effective face location s*. These quantities are obtained by interpolating each particle's specific volume (1/rho for equal-mass particles) along the line joining the pair.

This file implements the linear (1st order) interpolation, which only needs the volume at each particle (no gradients).

Reference:

• Inutsuka, S. (2002) "Reformulation of Smoothed Particle Hydrodynamics with Riemann Solver"

Definition in file reconstruction.hpp.

Function Documentation

lin_v2_sast_ij()

template<class Tscal>

EffectiveFace< Tscal > shammodels::gsph::lin_v2_sast_ij ( Tscal vol_a, Tscal vol_b, Tscal h_a, Tscal h_b, Tscal rab_inv )

inline

Linear (1st order) interpolation of the effective face between a pair.

Reconstructs the volume function V(s) linearly between the two particles using only their volume elements (no gradients needed), following Inutsuka (2002).

Template Parameters

Tscal

Scalar type

Parameters

vol_a	Specific volume of particle a (1/rho_a for equal-mass particles)
vol_b	Specific volume of particle b (1/rho_b for equal-mass particles)
h_a	Smoothing length of particle a
h_b	Smoothing length of particle b
rab_inv	Inverse of the pair separation (1/|r_a - r_b|)

Returns

EffectiveFace with V2 and s_ast

Definition at line 63 of file reconstruction.hpp.

Here is the call graph for this function: