Related papers: A simple solver for the two-fluid plasma model bas…
We present and analyze a parallel solver for the solution of fluid structure interaction problems described by a fictitious domain approach. In particular, the fluid is modeled by the non-stationary incompressible Navier-Stokes equations,…
Existence and admissibility of $\delta$-shock type solution is discussed for the following nonconvex strictly hyperbolic system arising in studues of plasmas: \pa_t u + \pa_x \big(\Sfrac{u^2+v^2}{2} \big) &=0 \pa_t v +\pa_x(v(u-1))&=0. The…
Three pseudospectral algorithms are described (Euler, leapfrog and trapez) for solving numerically the time dependent nonlinear Schroedinger equation in one, two or three dimensions. Numerical stability regions in the parameter space are…
We present a way to combine Vlasov and two-fluid codes for the simulation of a collisionless plasma in large domains while keeping full information of the velocity distribution in localized areas of interest. This is made possible by…
A collection of codes (in MATLAB & Fortran 77), and examples, for solving reaction-diffusion equations in one and two space dimensions is presented. In areas of the mathematical community spectral methods are used to remove the stiffness…
We develop a linearly-scaling variant of the Force Coupling Method [K. Yeo and M. R. Maxey, J. Fluid Mech. 649, 205-231 (2010)] for computing hydrodynamic interactions among particles confined to a doubly-periodic geometry with either a…
We propose a numerical algorithm that integrates quantum two-level systems (TLSs) into the finite-difference time-domain (FDTD) framework for simulating quantum emitters in arbitrary 3D photonic environments. Conventional methods struggle…
Quasi-periodic trajectories with two or more incommensurate frequencies are ubiquitous in nonlinear dynamics, yet the classical Fourier-based time-spectral method is tied to strictly periodic responses. We introduce a torus time-spectral…
We study a nonlinear coupled fluid-structure system modelling the blood flow through arteries. The fluid is described by the incompressible Navier-Stokes equations in a 2D rectangular domain where the upper part depends on a structure…
High-intensity laser plasma interactions create complex computational problems because they involve both fluid and kinetic regimes, which need models that maintain physical precision while keeping computational speed. The research…
We study the existence and multiplicity of solutions and the global solution curve of the following free boundary value problem, arising in plasma physics, see R. Temam [18], and H. Berestycki and H. Brezis [3]: find a function $u(x)$ and a…
In this work, the Navier-Stokes (NS) solver is combined with the Direct simulation Monte Carlo (DSMC) solver in a direct way, under the wave-particle formulation [J. Comput. Phys. 401, 108977 (2020)]. Different from the classical domain…
The Vlasov--Maxwell equations are used for the kinetic description of magnetized plasmas. As they are posed in an up to 3+3 dimensional phase space, solving this problem is extremely expensive from a computational point of view. In this…
A solver for the Poisson equation for 1D, 2D and 3D regular grids is presented. The solver applies the convolution theorem in order to efficiently solve the Poisson equation in spectral space over a rectangular computational domain.…
Two fluid simulations using local Landau-fluid closures derived from linear theory provide an efficient computational framework for plasma modelling, since they bridge the gap between computationally intensive kinetic simulations and fluid…
We present a six-moment multi-fluid model, which solves the governing equations for both ions and electrons, with pressure anisotropy along and perpendicular to the magnetic field direction, as well as the complete set of Maxwell equations.…
In this paper, we study the mathematical structure and numerical approximation of elliptic problems posed in a (3D) domain $\Omega$ when the right-hand side is a (1D) line source $\Lambda$. The analysis and approximation of such problems is…
The propagation and absorption of electromagnetic waves in plasma is one of the fundamental issues in plasma physics. The electromagnetic particle-in-cell method with the finite-difference time-domain solver plus Monte Carlo collision model…
We study and implement a simple method, based on the Perfectly Matched Layer approach, to treat non reflecting boundary conditions with the Smoothed Particles Hydrodynamics numerical algorithm. The method is based on the concept of physical…
The capability to accurately predict flood flows via numerical simulations is a key component of contemporary flood risk management practice. However, modern flood models lack the capacity to accurately model flow interactions with linear…