|
|
The work-in-progress implementation (on the `planetary` branch) of allowing different "gas" particles to have different equations of state, as well as the addition of more complex equations of state like iron and rock for planetary simulations such as giant impacts.
|
|
|
|
|
|
Status
|
|
|
------
|
|
|
* General structure implemented but incomplete and untested.
|
|
|
|
|
|
New Hydro Scheme
|
|
|
----------------
|
|
|
|
... | ... | @@ -37,3 +41,11 @@ Add a new integer flag particle property to identify the material of every parti |
|
|
Use the material ID flag with the EoS:
|
|
|
* `src/hydro/MinimalMultiMat/hydro.h`: Add `p->mat_id` as a third argument (after e.g. `denisty` and `u`) to all the equation of state functions, and to the functions themselves in `src/equation_of_state/tillotson/equation_of_state.h`. Here set the type `int` rather than `material_id` for generality.
|
|
|
* `src/equation_of_state/tillotson/equation_of_state.h`: In the EoS functions, declare a new `Til_params *mat` struct pointer then use a switch case to select which material's parameters to point it at and use in the equations.
|
|
|
|
|
|
To Do (Other than all the testing and fixing etc.)
|
|
|
--------------------------------------------------
|
|
|
|
|
|
* Add more equations of state (most of which will be read from tables).
|
|
|
* Integrate with more/all hydro schemes (careful with EoS that e.g. don't include entropy).
|
|
|
* (Not any time soon) Try SPH-modifying schemes to improve interactions of different-material particles.
|
|
|
|