Researcher profile

Giacomo Dimarco

Giacomo Dimarco contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
6works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

6 published item(s)

preprint2025arXiv

Micro-Macro Tensor Neural Surrogates for Uncertainty Quantification in Collisional Plasma

Plasma kinetic equations exhibit pronounced sensitivity to microscopic perturbations in model parameters and data, making reliable and efficient uncertainty quantification (UQ) essential for predictive simulations. However, the cost of uncertainty sampling, the high-dimensional phase space, and multiscale stiffness pose severe challenges to both computational efficiency and error control in traditional numerical methods. These aspects are further emphasized in presence of collisions where the high-dimensional nonlocal collision integrations and conservation properties pose severe constraints. To overcome this, we present a variance-reduced Monte Carlo framework for UQ in the Vlasov--Poisson--Landau (VPL) system, in which neural network surrogates replace the multiple costly evaluations of the Landau collision term. The method couples a high-fidelity, asymptotic-preserving VPL solver with inexpensive, strongly correlated surrogates based on the Vlasov--Poisson--Fokker--Planck (VPFP) and Euler--Poisson (EP) equations. For the surrogate models, we introduce a generalization of the separable physics-informed neural network (SPINN), developing a class of tensor neural networks based on an anisotropic micro-macro decomposition, to reduce velocity-moment costs, model complexity, and the curse of dimensionality. To further increase correlation with VPL, we calibrate the VPFP model and design an asymptotic-preserving SPINN whose small- and large-Knudsen limits recover the EP and VP systems, respectively. Numerical experiments show substantial variance reduction over standard Monte Carlo, accurate statistics with far fewer high-fidelity samples, and lower wall-clock time, while maintaining robustness to stochastic dimension.

preprint2020arXiv

Implicit-Explicit multistep methods for hyperbolic systems with multiscale relaxation

We consider the development of high order space and time numerical methods based on Implicit-Explicit (IMEX) multistep time integrators for hyperbolic systems with relaxation. More specifically, we consider hyperbolic balance laws in which the convection and the source term may have very different time and space scales. As a consequence the nature of the asymptotic limit changes completely, passing from a hyperbolic to a parabolic system. From the computational point of view, standard numerical methods designed for the fluid-dynamic scaling of hyperbolic systems with relaxation present several drawbacks and typically lose efficiency in describing the parabolic limit regime. In this work, in the context of Implicit-Explicit linear multistep methods we construct high order space-time discretizations which are able to handle all the different scales and to capture the correct asymptotic behavior, independently from its nature, without time step restrictions imposed by the fast scales. Several numerical examples confirm the theoretical analysis.

preprint2020arXiv

Social climbing and Amoroso distribution

We introduce a class of one-dimensional linear kinetic equations of Boltzmann and Fokker--Planck type, describing the dynamics of individuals of a multi-agent society questing for high status in the social hierarchy. At the Boltzmann level, the microscopic variation of the status of agents around a universal desired target, is built up introducing as main criterion for the change of status a suitable value function in the spirit of the prospect theory of Kahneman and Twersky. In the asymptotics of grazing interactions, the solution density of the Boltzmann type kinetic equation is shown to converge towards the solution of a Fokker--Planck type equation with variable coefficients of diffusion and drift, characterized by the mathematical properties of the value function. The steady states of the statistical distribution of the social status predicted by the Fokker--Planck equations belong to the class of Amoroso distributions with Pareto tails, which correspond to the emergence of a \emph{social elite}. The details of the microscopic kinetic interaction allow to clarify the meaning of the various parameters characterizing the resulting equilibrium. Numerical results then show that the steady state of the underlying kinetic equation is close to Amoroso distribution even in an intermediate regime in which interactions are not grazing.

preprint2019arXiv

Kinetic modeling of alcohol consumption

In most countries, alcohol consumption distributions have been shown to possess universal features. Their unimodal right-skewed shape is usually modeled in terms of the Lognormal distribution, which is easy to fit, test, and modify. However, empirical distributions often deviate considerably from the Lognormal model, and both Gamma and Weibull distributions appear to better describe the survey data. In this paper we explain the appearance of these distributions by means of classical methods of kinetic theory of multi-agent systems. The microscopic variation of alcohol consumption of agents around a universal \emph{social} accepted value of consumption, is built up introducing as main criterion for consumption a suitable value function in the spirit of the prospect theory of Kahneman and Twersky. The mathematical properties of the value function then determine the unique macroscopic equilibrium which results to be a generalized Gamma distribution. The modeling of the microscopic kinetic interaction allows to clarify the meaning of the various parameters characterizing the generalized Gamma equilibrium.

preprint2019arXiv

The Aw-Rascle traffic model: Enskog-type kinetic derivation and generalisations

We study the derivation of second order macroscopic traffic models from kinetic descriptions. In particular, we recover the celebrated Aw-Rascle model as the hydrodynamic limit of an Enskog-type kinetic equation out of a precise characterisation of the microscopic binary interactions among the vehicles. Unlike other derivations available in the literature, our approach unveils the multiscale physics behind the Aw-Rascle model. This further allows us to generalise it to a new class of second order macroscopic models complying with the Aw-Rascle consistency condition, namely the fact that no wave should travel faster than the mean traffic flow.

preprint2010arXiv

Exponential Runge-Kutta methods for stiff kinetic equations

We introduce a class of exponential Runge-Kutta integration methods for kinetic equations. The methods are based on a decomposition of the collision operator into an equilibrium and a non equilibrium part and are exact for relaxation operators of BGK type. For Boltzmann type kinetic equations they work uniformly for a wide range of relaxation times and avoid the solution of nonlinear systems of equations even in stiff regimes. We give sufficient conditions in order that such methods are unconditionally asymptotically stable and asymptotic preserving. Such stability properties are essential to guarantee the correct asymptotic behavior for small relaxation times. The methods also offer favorable properties such as nonnegativity of the solution and entropy inequality. For this reason, as we will show, the methods are suitable both for deterministic as well as probabilistic numerical techniques.