Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

collisions: Disable ion-neutral and electron-neutral by default #242

Merged
merged 2 commits into from
Aug 8, 2024
Merged
Show file tree
Hide file tree
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
9 changes: 8 additions & 1 deletion docs/sphinx/components.rst
Original file line number Diff line number Diff line change
Expand Up @@ -969,7 +969,8 @@ listed after all the species groups in the component list, so that all
the species are present in the state.

One of the most important is the `collisions`_ component. This sets collision
times for all species, which are then used
times for all species, which are then used in other components to calculate
quantities like heat diffusivities and viscosity closures.

.. _sound_speed:

Expand Down Expand Up @@ -1176,12 +1177,18 @@ species :math:`b` due to temperature differences, is given by:

- Ion-neutral and electron-neutral collisions

*Note*: These are disabled by default. If enabled, care is needed to
avoid double-counting collisions in atomic reactions e.g charge-exchange
reactions.

The cross-section for elastic collisions between charged and neutral
particles can vary significantly. Here for simplicity we just take
a value of :math:`5\times 10^{-19}m^2` from the NRL formulary.

- Neutral-neutral collisions

*Note* This is enabled by default.

The cross-section is given by

.. math::
Expand Down
4 changes: 2 additions & 2 deletions src/collisions.cxx
Original file line number Diff line number Diff line change
Expand Up @@ -33,13 +33,13 @@ Collisions::Collisions(std::string name, Options& alloptions, Solver*) {
.withDefault<bool>(true);
electron_neutral = options["electron_neutral"]
.doc("Include electron-neutral elastic collisions?")
.withDefault<bool>(true);
.withDefault<bool>(false);
ion_ion = options["ion_ion"]
.doc("Include ion-ion elastic collisions?")
.withDefault<bool>(true);
ion_neutral = options["ion_neutral"]
.doc("Include ion-neutral elastic collisions?")
.withDefault<bool>(true);
.withDefault<bool>(false);
neutral_neutral = options["neutral_neutral"]
.doc("Include neutral-neutral elastic collisions?")
.withDefault<bool>(true);
Expand Down
Loading