A model is proposed to describe the spike-frequency adaptation observed in many neuronal systems. We assume that adaptation is mainly due to a calcium-activated potassium current, and we consider two coupled stochastic differential equations for which an analytical approach combined with simulation techniques and numerical methods allow to obtain both qualitative and quantitative results about asymptotic mean firing rate, mean calcium concentration and the firing probability density. A related algorithm, based on the Hazard Rate Method, is also devised and described.

Direct numerical approximation of a continuous-time infinite horizon control problem, requires to recast the model as a discrete-time, finite-horizon control model. The quality of the optimization results can be heavily degraded if the discretization process does not take into account features of the original model to be preserved.

Numerical results on the translocation of long biopolymers through mid-sized and wide pores are presented. The simulations are based on a novel methodology which couples molecular motion to a mesoscopic fluid solvent. Thousands of events of long polymers (up to 8000 monomers) are monitored as they pass through nanopores. Comparison between the different pore sizes shows that wide pores can host a larger number of multiple biopolymer segments, as compared to smaller pores.

In this paper we deal with high oscillatory systems and numerical methods for the approximation of their solutions. Some classical schemes developed in the literature are recalled and a recent approach based on the expression of the oscillatory solution by means of the exponential map is considered. Moreover we introduce a new method based on the Cayley map and provide some numerical tests in order to compare the different approaches

We provide numerical evidence that electronic preturbulent phenomena in graphene could be observed, under current experimental conditions, through current fluctuations, echoing the detachment of vortices past localized micron-sized impurities. Vortex generation, due to micron-sized constriction, is also explored with special focus on the effects of relativistic corrections to the normal Navier-Stokes equations. These corrections are found to cause a delay in the stability breakout of the fluid as well as a small shift in the vortex shedding frequency.

We present a highly optimized implementation of a Monte Carlo (MC) simulator for the three-dimensional Ising spin-glass model with bimodal disorder, i.e.; the 3D Edwards-Anderson model running on CUDA enabled GPUs. Multi-GPU systems exchange data by means of the Message Passing Interface (MPI). The chosen MC dynamics is the classic Metropolis one, which is purely dissipative, since the aim was the study of the critical off-equilibrium relaxation of the system.

We present a lattice Boltzmann realization of Grad's extended hydrodynamic approach to nonequilibrium flows. This is achieved by using higher-order isotropic lattices coupled with a higher-order regularization procedure. The method is assessed for flow across parallel plates and three-dimensional flows in porous media, showing excellent agreement of the mass flow with analytical and numerical solutions of the Boltzmann equation across the full range of Knudsen numbers, from the hydrodynamic regime to ballistic motion.