shamrock.external.coala.solver_DG.solver_coag_kdv

shamrock.external.coala.solver_DG.solver_coag_kdv(eps, nbins, kpol, massgrid, massbins, gij, tensor_tabflux_coag, tensor_tabintflux_coag, dv, dt)

Function to compute SSPRK order 3 time solver with piecewise constant approximation Function for ballistic kernel with differential velocities dv See Zhang & Shu 2010 and Lombart et al., 2021

Parameters

epsscalar, type -> float

minimum value for mass distribution approximation gij

nbinsscalar, type -> integer

number of dust bins

kpolscalar, type -> integer

degree of polynomials for approximation

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

grid of masses given borders value of mass bins

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

arithmetic mean value of massgrid for each mass bins

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

components of g on the polynomial basis

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

array to evaluate coagulation flux

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

array to evaluate the term including the integral of coagulation flux

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

array of the differential velocity between grains

dtscalar, type -> float

timestep

Returns

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

evolved components of g on the polynomial basis after 1 timestep