Source author record

Giacomo Dimarco

Giacomo Dimarco appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

19works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

19 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.

preprint2016arXiv

A new class of high order semi-Lagrangian schemes for rarefied gas dynamics

In this paper we genealize the fast semi-Lagrangian scheme developed in [J. Comput. Phys., Vol. 255, 2013, pp 680-698] to the case of high order reconstructions of the distribution function. The original first order accurate semi-Lagrangian scheme is supplemented with polynomial reconstructions of the distribution function and of the collisional operator leading to an effective high order accurate numerical scheme for all regimes, from extremely rarefied gas to highly collisional siuation. The main idea relies on updating at each time step the extreme points of the distribution function for each velocity of the lattice instead of updating the solution in the cell centers, these extremes points being located at different positions for any fixed velocity of the lattice. The result is a class of scheme which permits to preserve the structure of the solution over very long times compared to existing schemes from the literature. We propose a proof of concept of this new approach along with numerical tests and comparisons with classical numerical methods.

preprint2016arXiv

An efficient numerical method for solving the Boltzmann equation in multidimensions

In this paper we deal with the extension of the Fast Kinetic Scheme (FKS) [J. Comput. Phys., Vol. 255, 2013, pp 680-698] originally constructed for solving the BGK equation, to the more challenging case of the Boltzmann equation. The scheme combines a robust and fast method for treating the transport part based on an innovative Lagrangian technique supplemented with fast spectral schemes to treat the collisional operator by means of an operator splitting approach. This approach along with several implementation features related to the parallelization of the algorithm permits to construct an efficient simulation tool which is numerically tested against exact and reference solutions on classical problems arising in rarefied gas dynamic. We present results up to the $3$D$\times 3$D case for unsteady flows for the Variable Hard Sphere model which may serve as benchmark for future comparisons between different numerical methods for solving the multidimensional Boltzmann equation. For this reason, we also provide for each problem studied details on the computational cost and memory consumption as well as comparisons with the BGK model or the limit model of compressible Euler equations.

preprint2016arXiv

Implicit-explicit linear multistep methods for stiff kinetic equations

We consider the development of high order asymptotic-preserving linear multistep methods for kinetic equations and related problems. The methods are first developed for BGK-like kinetic models and then extended to the case of the full Boltzmann equation. The behavior of the schemes in the Navier-Stokes regime is also studied and compatibility conditions derived. We show that, compared to IMEX Runge-Kutta methods, the IMEX multistep schemes have several advantages due to the absence of coupling conditions and to the greater computational efficiency. The latter is of paramount importance when dealing with the time discretization of multidimensional kinetic equations.

preprint2014arXiv

Macroscopic models of collective motion with repulsion

We study a system of self-propelled particles which interact with their neighbors via alignment and repulsion. The particle velocities result from self-propulsion and repulsion by close neighbors. The direction of self-propulsion is continuously aligned to that of the neighbors, up to some noise. A continuum model is derived starting from a mean-field kinetic description of the particle system. It leads to a set of non conservative hydrodynamic equations. We provide a numerical validation of the continuum model by comparison with the particle model. We also provide comparisons with other self-propelled particle models with alignment and repulsion.

preprint2013arXiv

Hydrodynamics of the Kuramoto-Vicsek model of rotating self-propelled particles

We consider an Individual-Based Model for self-rotating particles interacting through local alignment and investigate its macroscopic limit. This model describes self-propelled particles moving in the plane and trying to synchronize their rotation motion with their neighbors. It combines the Kuramoto model of synchronization and the Vicsek model of swarm formation. We study the mean-field kinetic and hydrodynamic limits of this system within two different scalings. In the small angular velocity regime, the resulting model is a slight modification of the 'Self-Organized Hydrodynamic' model which has been previously introduced by the first author. In the large angular velocity case, a new type of hydrodynamic model is obtained. A preliminary study of the linearized stability is proposed.

preprint2013arXiv

The Moment Guided Monte Carlo method for the Boltzmann equation

In this work we propose a generalization of the Moment Guided Monte Carlo method developed in [11]. This approach permits to reduce the variance of the particle methods through a matching with a set of suitable macroscopic moment equations. In order to guarantee that the moment equations provide the correct solutions, they are coupled to the kinetic equation through a non equilibrium term. Here, at the contrary to the previous work in which we considered the simplified BGK operator, we deal with the full Boltzmann operator. Moreover, we introduce an hybrid setting which permits to entirely remove the resolution of the kinetic equation in the limit of infinite number of collisions and to consider only the solution of the compressible Euler equation. This modification additionally reduce the statistical error with respect to our previous work and permits to perform simulations of non equilibrium gases using only a few number of particles. We show at the end of the paper several numerical tests which prove the efficiency and the low level of numerical noise of the method.

preprint2012arXiv

High order asymptotic-preserving schemes for the Boltzmann equation

In this note we discuss the construction of high order asymptotic preserving numerical schemes for the Boltzmann equation. The methods are based on the use of Implicit-Explicit (IMEX) Runge-Kutta methods combined with a penalization technique recently introduced in [F. Filbet, S. Jin: A class of asymptotic preserving schemes for kinetic equations and related problems with stiff sources,J. Comp. Phys. 229, (2010), pp. 7625-7648.].

preprint2012arXiv

Towards an ultra efficient kinetic scheme. Part I: basics on the BGK equation

In this paper we present a new ultra efficient numerical method for solving kinetic equations. In this preliminary work, we present the scheme in the case of the BGK relaxation operator. The scheme, being based on a splitting technique between transport and collision, can be easily extended to other collisional operators as the Boltzmann collision integral or to other kinetic equations such as the Vlasov equation. The key idea, on which the method relies, is to solve the collision part on a grid and then to solve exactly the transport linear part by following the characteristics backward in time. The main difference between the method proposed and semi-Lagrangian methods is that here we do not need to reconstruct the distribution function at each time step. This allows to tremendously reduce the computational cost of the method and it permits for the first time, to the author's knowledge, to compute solutions of full six dimensional kinetic equations on a single processor laptop machine. Numerical examples, up to the full three dimensional case, are presented which validate the method and assess its efficiency in 1D, 2D and 3D.

preprint2010arXiv

Direct simulation Monte Carlo schemes for Coulomb interactions in plasmas

We consider the development of Monte Carlo schemes for molecules with Coulomb interactions. We generalize the classic algorithms of Bird and Nanbu-Babovsky for rarefied gas dynamics to the Coulomb case thanks to the approximation introduced by Bobylev and Nanbu (Theory of collision algorithms for gases and plasmas based on the Boltzmann equation and the Landau-Fokker-Planck equation, Physical Review E, Vol. 61, 2000). Thus, instead of considering the original Boltzmann collision operator, the schemes are constructed through the use of an approximated Boltzmann operator. With the above choice larger time steps are possible in simulations; moreover the expensive acceptance-rejection procedure for collisions is avoided and every particle collides. Error analysis and comparisons with the original Bobylev-Nanbu (BN) scheme are performed. The numerical results show agreement with the theoretical convergence rate of the approximated Boltzmann operator and the better performance of Bird-type schemes with respect to the original scheme.

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.

preprint2010arXiv

Fluid Simulations with Localized Boltzmann Upscaling by Direct Simulation Monte-Carlo

In the present work, we present a novel numerical algorithm to couple the Direct Simulation Monte Carlo method (DSMC) for the solution of the Boltzmann equation with a finite volume like method for the solution of the Euler equations. Recently we presented in [14],[16],[17] different methodologies which permit to solve fluid dynamics problems with localized regions of departure from thermodynamical equilibrium. The methods rely on the introduction of buffer zones which realize a smooth transition between the kinetic and the fluid regions. In this paper we extend the idea of buffer zones and dynamic coupling to the case of the Monte Carlo methods. To facilitate the coupling and avoid the onset of spurious oscillations in the fluid regions which are consequences of the coupling with a stochastic numerical scheme, we use a new technique which permits to reduce the variance of the particle methods [11]. In addition, the use of this method permits to obtain estimations of the breakdowns of the fluid models less affected by fluctuations and consequently to reduce the kinetic regions and optimize the coupling. In the last part of the paper several numerical examples are presented to validate the method and measure its computational performances.

preprint2009arXiv

A Multiscale Kinetic-Fluid Solver with Dynamic Localization of Kinetic Effects

This paper collects the efforts done in our previous works [P. Degond, S. Jin, L. Mieussens, A Smooth Transition Between Kinetic and Hydrodynamic Equations, J. Comp. Phys., 209 (2005) 665--694.],[P.Degond, G. Dimarco, L. Mieussens, A Moving Interface Method for Dynamic Kinetic-fluid Coupling, J. Comp. Phys., Vol. 227, pp. 1176-1208, (2007).],[P. Degond, J.G. Liu, L. Mieussens, Macroscopic Fluid Model with Localized Kinetic Upscaling Effects, SIAM Multi. Model. Sim. 5(3), 940--979 (2006)] to build a robust multiscale kinetic-fluid solver. Our scope is to efficiently solve fluid dynamic problems which present non equilibrium localized regions that can move, merge, appear or disappear in time. The main ingredients of the present work are the followings ones: a fluid model is solved in the whole domain together with a localized kinetic upscaling term that corrects the fluid model wherever it is necessary; this multiscale description of the flow is obtained by using a micro-macro decomposition of the distribution function [P. Degond, J.G. Liu, L. Mieussens, Macroscopic Fluid Model with Localized Kinetic Upscaling Effects, SIAM Multi. Model. Sim. 5(3), 940--979 (2006)]; the dynamic transition between fluid and kinetic descriptions is obtained by using a time and space dependent transition function; to efficiently define the breakdown conditions of fluid models we propose a new criterion based on the distribution function itself. Several numerical examples are presented to validate the method and measure its computational efficiency.

preprint2009arXiv

Fluid Solver Independent Hybrid Methods for Multiscale Kinetic equations

In some recent works [G. Dimarco, L. Pareschi, Hybrid multiscale methods I. Hyperbolic Relaxation Problems, Comm. Math. Sci., 1, (2006), pp. 155-177], [G. Dimarco, L. Pareschi, Hybrid multiscale methods II. Kinetic equations, SIAM Multiscale Modeling and Simulation Vol 6., No 4,pp. 1169-1197, (2008)] we developed a general framework for the construction of hybrid algorithms which are able to face efficiently the multiscale nature of some hyperbolic and kinetic problems. Here, at variance with respect to the previous methods, we construct a method form-fitting to any type of finite volume or finite difference scheme for the reduced equilibrium system. Thanks to the coupling of Monte Carlo techniques for the solution of the kinetic equations with macroscopic methods for the limiting fluid equations, we show how it is possible to solve multiscale fluid dynamic phenomena faster with respect to traditional deterministic/stochastic methods for the full kinetic equations. In addition, due to the hybrid nature of the schemes, the numerical solution is affected by less fluctuations when compared to standard Monte Carlo schemes. Applications to the Boltzmann-BGK equation are presented to show the performance of the new methods in comparison with classical approaches used in the simulation of kinetic equations.

preprint2009arXiv

The Moment Guided Monte Carlo Method

In this work we propose a new approach for the numerical simulation of kinetic equations through Monte Carlo schemes. We introduce a new technique which permits to reduce the variance of particle methods through a matching with a set of suitable macroscopic moment equations. In order to guarantee that the moment equations provide the correct solutions, they are coupled to the kinetic equation through a non equilibrium term. The basic idea, on which the method relies, consists in guiding the particle positions and velocities through moment equations so that the concurrent solution of the moment and kinetic models furnishes the same macroscopic quantities.