Vitamin D has been proven to be a strong stimulator of mechanisms associated with the elimination of pathogens. Because of its recognized effectiveness against viral infections, during SARS-CoV-2 infection, the effects of Vitamin D supplementation have been the object of debate. This study aims to contribute to this debate by the means of a qualitative phenomenological mathematical model in which the role of Vitamin D and its interactions with the innate immune system are explicitly considered.

The present work extends a previous paper where an agent-based and two-dimensional partial differential diffusion model was introduced for describing immune cell dynamics (leukocytes) in cancer-on-chip experiments. In the present work, new features are introduced for the dynamics of leukocytes and for their interactions with tumor cells, improving the adherence of the model to what is observed in laboratory experiments. Each system's solution realization is a family of biased random walk trajectories, affected by the chemotactic gradients and in turn affecting them.

The simulations for the inverse problem of radiative transfer, even if built with a correct Bayesian approach, do not represent the full source of errors present in the experimental data. We point out two categories of errors (atmospheric model errors and non-Gaussian instrumental errors due to the optics and hardware, that are not considered by standard methods. Moreover, we show cases taken from FORUM simulated radiances using an End to End simulator, where se show how the instrument reacts to a non homogeneousneous filed of view.

To evaluate changes in the Soil Organic Carbon (SOC) index, one of the key indicators of land degradation neutrality, soil carbon modeling is of primary importance. In litera-ture, the analysis has been focused on the stability characterization of soil carbon steady states and in the calculation of the resilience of the stable equilibria. Neither stability nor resilience, however, provide any information about transient dynamics, and models with highly resilient equilibria can exhibit dramatic transient responses to perturbations.

In this paper, an algorithm for the numerical evaluation of hypersingular finite-part integrals with rapidly oscillating kernels is proposed. The method is based on an interpolatory procedure at zeros of the orthogonal polynomials with respect to the first kind Chebyshev weight. Bounds of the error and of the amplification factor are also provided. Numerically stable procedure are obtained and the corresponding algorithms can be implemented in a fast way.

During melting under gravity in the presence of a horizontal thermal gradient, buoyancy-driven convection in the liquid phase affects significantly the evolution of the liquid-solid interface. Due to the obvious engineering interest in understanding and controlling melting processes, fluid dynamicists and applied mathematicians have spent many efforts to model and simulate them numerically. Their endeavors concentrated in the twenty-five years period between the publication of the paper by Brent, Voller & Reid (1988) and that by Mansutti & Bucchignani (2011).

Recent literature confirms the crucial influence of non-viscous deformations together with temperature impact on glacier and rock glacier flow numerical simulation. Along this line, supported by the successful test on a one-dimensional set-up developed by two of the author, we propose the numerical solution of a two-dimensional rock-glacier flow model based on an ice constitutive law of second grade differential type .