Module for photoionization with the Helmholtz approximation. More...
Functions/Subroutines | |
| subroutine, public | photoi_helmh_initialize (tree, cfg, is_used, eta) |
| Initialize options for Helmholtz photoionization. | |
| subroutine, public | photoi_helmh_compute (tree, i_photo) |
| subroutine, public | photoi_helmh_bc (box, nb, iv, coords, bc_val, bc_type) |
Detailed Description
Module for photoionization with the Helmholtz approximation.
The equations solved are nabla^2 phi_i - lambda_i^2 * phi_i = f, where f is the photon production term, and there can be modes i = 1, ..., N. The photon absorption term is then given by sum(c_i * phi_i).
For the case N=2, the following coefficients can be used: lambdas = 4425.36_dp, 750.06_dp 1/(m bar) coeffs = 337555.5_dp, 19972.0_dp 1/(m bar)^2
TODO: look up values for the case N=3
Function/Subroutine Documentation
◆ photoi_helmh_bc()
| subroutine, public m_photoi_helmh::photoi_helmh_bc | ( | type(box_t), intent(in) | box, |
| integer, intent(in) | nb, | ||
| integer, intent(in) | iv, | ||
| real(dp), dimension(ndim, box%n_cell**(ndim-1)), intent(in) | coords, | ||
| real(dp), dimension(box%n_cell**(ndim-1)), intent(out) | bc_val, | ||
| integer, intent(out) | bc_type ) |
- Todo
- Think about good (and efficient) boundary conditions for Helmholtz equations, see also Bourdon et al. PSST 2017. Setting a Dirichlet zero b.c. is inaccurate if the streamer gets close to that boundary, otherwise it should be quite reasonable.
Definition at line 210 of file m_photoi_helmh.f90.
◆ photoi_helmh_compute()
| subroutine, public m_photoi_helmh::photoi_helmh_compute | ( | type(af_t), intent(inout) | tree, |
| integer, intent(in) | i_photo ) |
Definition at line 162 of file m_photoi_helmh.f90.
◆ photoi_helmh_initialize()
| subroutine, public m_photoi_helmh::photoi_helmh_initialize | ( | type(af_t), intent(inout) | tree, |
| type(cfg_t), intent(inout) | cfg, | ||
| logical, intent(in) | is_used, | ||
| real(dp), intent(in) | eta ) |
Initialize options for Helmholtz photoionization.
- Parameters
-
[in] is_used Whether the module is used [in] eta Efficiency
Bourdon two terms lambdas in SI units are: [7305.62,44081.25] Bourdon two terms coeffs in SI units are: [11814508.38,998607256]
Bourdon three terms lambdas in SI units are: [4147.85 10950.93 66755.67] bourdon three terms coeffs in SI units are: [ 1117314.935 28692377.5 2748842283 ]
Definition at line 44 of file m_photoi_helmh.f90.
