shamrock.external.coala.solver_DG.L_func_coag_kdv

shamrock.external.coala.solver_DG.L_func_coag_kdv(nbins, kpol, massgrid, gij, tensor_tabflux_coag, tensor_tabintflux_coag, dv)

Function to compute the DG operator L for piecewise polynomial approximation (see Lombart et al., 2021) Function for ballistic kernel with differential velocities dv It is used for the time solver

Parameters

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

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

Returns

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

DG operator for piecewise polynomials approximation in each bin