PartMC  2.6.1
Functions/Subroutines | Variables
pmc_coag_kernel_constant Module Reference

Constant coagulation kernel. More...

Functions/Subroutines

subroutine kernel_constant (aero_particle_1, aero_particle_2, aero_data, env_state, k)
 Constant coagulation kernel. More...
 
subroutine kernel_constant_minmax (v1, v2, aero_data, env_state, k_min, k_max)
 Minimum and maximum values of the constant coagulation kernel. More...
 
subroutine soln_constant_exp (bin_grid, aero_data, time, num_conc, radius_at_mean_vol, env_state, aero_binned)
 Exact solution with a constant coagulation kernel and an exponential initial condition. More...
 

Variables

real(kind=dp), parameter beta_0 = 0.25d0 / (60d0 * 2d8)
 Coefficient for constant kernel. More...
 

Detailed Description

Constant coagulation kernel.

Function/Subroutine Documentation

◆ kernel_constant()

subroutine pmc_coag_kernel_constant::kernel_constant ( type(aero_particle_t), intent(in)  aero_particle_1,
type(aero_particle_t), intent(in)  aero_particle_2,
type(aero_data_t), intent(in)  aero_data,
type(env_state_t), intent(in)  env_state,
real(kind=dp), intent(out)  k 
)

Constant coagulation kernel.

Parameters
[in]aero_particle_1First particle.
[in]aero_particle_2Second particle.
[in]aero_dataAerosol data.
[in]env_stateEnvironment state.
[out]kCoagulation kernel.

Definition at line 29 of file coag_kernel_constant.F90.

◆ kernel_constant_minmax()

subroutine pmc_coag_kernel_constant::kernel_constant_minmax ( real(kind=dp), intent(in)  v1,
real(kind=dp), intent(in)  v2,
type(aero_data_t), intent(in)  aero_data,
type(env_state_t), intent(in)  env_state,
real(kind=dp), intent(out)  k_min,
real(kind=dp), intent(out)  k_max 
)

Minimum and maximum values of the constant coagulation kernel.

Parameters
[in]v1Volume of first particle.
[in]v2Volume of second particle.
[in]aero_dataAerosol data.
[in]env_stateEnvironment state.
[out]k_minCoagulation kernel minimum value.
[out]k_maxCoagulation kernel maximum value.

Definition at line 50 of file coag_kernel_constant.F90.

◆ soln_constant_exp()

subroutine pmc_coag_kernel_constant::soln_constant_exp ( type(bin_grid_t), intent(in)  bin_grid,
type(aero_data_t), intent(in)  aero_data,
real(kind=dp), intent(in)  time,
real(kind=dp), intent(in)  num_conc,
real(kind=dp), intent(in)  radius_at_mean_vol,
type(env_state_t), intent(in)  env_state,
type(aero_binned_t), intent(inout)  aero_binned 
)

Exact solution with a constant coagulation kernel and an exponential initial condition.

Given input paramaters $R$ and $N_0$ we let the mean volume be $v_\mu = \frac{4\pi}{3} R^3$ and define the rescaled time $\tau = N_0 \beta_0 t$, where $\beta_0$ is the fixed constant kernel value. We also set the parameter $\lambda = 1$. Then the solution is

\[ n(D,t) \ {\rm d}\ln D = \frac{\pi}{2} D^3 \frac{N_0}{v_\mu} \frac{4}{(\tau + 2)^2} \exp\left(-\frac{v}{v_\mu} \frac{2}{\tau + 2} \exp(-\lambda \tau) - \lambda \tau\right) {\rm d}\ln D \]

This thus has initial condition

\[ n(D,t) \ {\rm d}\ln D = \frac{\pi}{2} D^3 \frac{N_0}{v_\mu} \exp\left(-\frac{v}{v_\mu}\right) {\rm d}\ln D \]

Parameters
[in]bin_gridBin grid.
[in]aero_dataAerosol data.
[in]timeCurrent time.
[in]num_concParticle number concentration (#/m^3).
[in]radius_at_mean_volMean init radius (m).
[in]env_stateEnvironment state.
[in,out]aero_binnedOutput state.

Definition at line 92 of file coag_kernel_constant.F90.

Variable Documentation

◆ beta_0

real(kind=dp), parameter pmc_coag_kernel_constant::beta_0 = 0.25d0 / (60d0 * 2d8)

Coefficient for constant kernel.

Definition at line 22 of file coag_kernel_constant.F90.