Related papers: Kinetic Scheme for Solving M1 Model of Radiative T…
We present a hierarchical approach for enhancing the robustness of numerical solvers for modelling radiative MHD flows in multi-dimensions. This approach is based on clustering the entries of the global Jacobian in a hierarchical manner…
Aims. Numerical test-particle simulations are a reliable and frequently used tool to test analytical transport theories and to predict mean-free paths. The comparison between solutions of the diffusion equation and the particle flux is used…
The ideal MHD equations are a central model in astrophysics, and their solution relies upon stable numerical schemes. We present an implementation of a new method, which possesses excellent stability properties. Numerical tests demonstrate…
Frequency-dependent/hybrid approaches for stellar irradiation are of primary importance in numerical simulations of massive star formation. We seek to compare outflow and accretion mechanisms in star formation simulations. We investigate…
In this article we study the numerical solution of the $L^1$-Optimal Transport Problem on 2D surfaces embedded in $R^3$, via the DMK formulation introduced in [FaccaCardinPutti:2018]. We extend from the Euclidean into the Riemannian setting…
We present a new algorithm based on a Cartesian mesh for the numerical approximation of kinetic models for chemosensitive movements set in an arbitrary geometry. We investigate the influence of the geometry on the collective behavior of…
A direct numerical solution of the radiative transfer equation or any kinetic equation is typically expensive, since the radiative intensity depends on time, space and direction. An expansion in the direction variables yields an equivalent…
This article presents a high order conservative flux optimization (CFO) finite element method for the elliptic diffusion equations. The numerical scheme is based on the classical Galerkin finite element method enhanced by a flux…
For an accurate treatment of the shock wave propagation in high-energy astrophysical phenomena, such as supernova shock breakouts, gamma-ray bursts and accretion disks, knowledge of radiative transfer plays a crucial role. In this paper we…
We study an inverse problem for Light Sheet Fluorescence Microscopy (LSFM), where the density of fluorescent molecules needs to be reconstructed. Our first step is to present a mathematical model to describe the measurements obtained by an…
We introduce a new code for computing time-dependent continuum radiative transfer and non-equilibrium ionization states in static density fields with periodic boundaries. Our code solves the moments of the radiative transfer equation,…
We develop, and implement in a Finite Volume environment, a density-based approach for the Euler equations written in conservative form using density, momentum, and total energy as variables. Under simplifying assumptions, these equations…
In this paper, we present a kinetic model with flexible velocities that satisfy positivity preservation conditions for the Euler equations. Our 1D kinetic model consists of two velocities and employs both the asymmetrical and symmetrical…
This paper aims to develop and analyze a numerical scheme for solving the backward problem of semilinear subdiffusion equations. We establish the existence, uniqueness, and conditional stability of the solution to the inverse problem by…
Due to the nonlocal feature of fractional differential operators, the numerical solution to fractional partial differential equations usually requires expensive memory and computation costs. This paper develops a fast scheme for fractional…
In our previous work, numerical schemes for a simplified version of 3-wave kinetic equations, in which only the simple forward-cascade terms of the collision operators are kept, have been successfully designed, especially to capture the…
Two-stream radiative transfer is a workhorse in the Earth, planetary and exoplanetary sciences due to its simplicity and ease of implementation. However, a longstanding limitation of the two-stream approximation is its inaccuracy in the…
We present a novel numerical implementation of radiative transfer in the cosmological smoothed particle hydrodynamics (SPH) simulation code {\small GADGET}. It is based on a fast, robust and photon-conserving integration scheme where the…
In this work, we prove rigorous error estimates for a hybrid method introduced in [15] for solving the time-dependent radiation transport equation (RTE). The method relies on a splitting of the kinetic distribution function for the…
We discuss the problem of ultra high energy particles propagation in astrophysical backgrounds. We present two different computational schemes based on both kinetic and Monte Carlo approaches. The kinetic approach is an analytical…