shamrock.model_gsph.GSPHModel_f64_3_M6_SolverConfig#
- class shamrock.model_gsph.GSPHModel_f64_3_M6_SolverConfig#
Bases:
pybind11_objectMethods
Set adiabatic equation of state: P = (gamma-1) rho u
Set isothermal equation of state: P = cs^2 rho
Set first-order piecewise constant reconstruction.
Set HLLC approximate Riemann solver.
Set iterative Riemann solver (van Leer 1997).
- add_ext_force_point_mass(self: shamrock.pyshamrock.model_gsph.GSPHModel_f64_3_M6_SolverConfig, *, central_mass: float, Racc: float) None#
- set_cfl_cour(self: shamrock.pyshamrock.model_gsph.GSPHModel_f64_3_M6_SolverConfig, arg0: float) None#
- set_cfl_force(self: shamrock.pyshamrock.model_gsph.GSPHModel_f64_3_M6_SolverConfig, arg0: float) None#
- set_eos_adiabatic(self: shamrock.pyshamrock.model_gsph.GSPHModel_f64_3_M6_SolverConfig, gamma: float) None#
Set adiabatic equation of state: P = (gamma-1) rho u
Parameters#
- gammafloat
Adiabatic index (e.g., 5/3 for monatomic gas, 7/5 for diatomic)
- set_eos_isothermal(self: shamrock.pyshamrock.model_gsph.GSPHModel_f64_3_M6_SolverConfig, cs: float) None#
Set isothermal equation of state: P = cs^2 rho
Parameters#
- csfloat
Sound speed
- set_particle_mass(self: shamrock.pyshamrock.model_gsph.GSPHModel_f64_3_M6_SolverConfig, arg0: float) None#
- set_reconstruct_piecewise_constant(self: shamrock.pyshamrock.model_gsph.GSPHModel_f64_3_M6_SolverConfig) None#
Set first-order piecewise constant reconstruction.
Sets all gradients to zero. Most diffusive but most stable. Good for very strong shocks or initial testing.
- set_riemann_hllc(self: shamrock.pyshamrock.model_gsph.GSPHModel_f64_3_M6_SolverConfig) None#
Set HLLC approximate Riemann solver.
Fast approximate Riemann solver that captures contact discontinuities. Recommended for general use - good balance of accuracy and speed.
- set_riemann_iterative(self: shamrock.pyshamrock.model_gsph.GSPHModel_f64_3_M6_SolverConfig, *, tolerance: float = 1e-06, max_iter: int = 20) None#
Set iterative Riemann solver (van Leer 1997).
This is the most accurate but slower Riemann solver. Uses Newton-Raphson iteration to find the pressure in the star region.
Parameters#
- tolerancefloat
Convergence tolerance for Newton-Raphson iteration (default: 1e-6)
- max_iterint
Maximum number of iterations (default: 20)
- set_tree_reduction_level(self: shamrock.pyshamrock.model_gsph.GSPHModel_f64_3_M6_SolverConfig, arg0: int) None#
- set_two_stage_search(self: shamrock.pyshamrock.model_gsph.GSPHModel_f64_3_M6_SolverConfig, arg0: bool) None#