shamrock.external.coala.interface_coala_shamrock.coala_source_term_k0

shamrock.external.coala.interface_coala_shamrock.coala_source_term_k0(nbins, massgrid, rhodust, rhodust_eps, tensor_tabflux_coag, v_dust)

Function to compute the source for coagulation and fragmentation in continuity equation for piecewise constant approximation (see Lombart et al., 2021) Function for ballistic kernel with differential velocities dv Used to evaluate the source term, then hydro code applies time solver

/!Only coagulation so far

Parameters

nbinsscalar, type -> integer

number of dust bins

massgrid1D array (dim = nbins+1), type -> float

grid of masses given borders value of mass bins

rhodust1D array (dim = nbins), type -> float

dust density for each grain size

rhodust_epsscalar, type -> float

threshold value for rhodust

tensor_tabflux_coag3D array (dim = (nbins,nbins,nbins)), type -> float

array to evaluate coagulation flux

v_dust1D array (dim = (nbins)), type -> float

array of the dust velocities (could also be delta_v in monofluid since it is a delta)

Returns

S_coag1D array (dim = nbins), type -> float

Source term for dust coagulation in continuity equation DG operator for piecewise constant approximation in each binls