Related papers: A higher-order accurate operator splitting spectra…
We develop a spectral method to solve the heat equation in a closed cylinder, achieving a quasi-optimal $\mathcal{O}(N\log N)$ complexity and high-order, spectral accuracy. The algorithm relies on a Chebyshev--Chebyshev--Fourier (CCF)…
We extend our framework for 3D radiative transfer calculations with a non-local operator splitting methods along (full) characteristics to spherical and cylindrical coordinate systems. These coordinate systems are better suited to a number…
With the aim of efficiently simulating three-dimensional multiphase turbulent flows with a phase-field method, we propose a new discretization scheme for the biharmonic term (the 4th-order derivative term) of the Cahn-Hilliard equation.…
We present a comprehensive study of semiclassical phase-space propagation in the Wigner representation, emphasizing numerical applications, in particular as an initial-value representation. Two semiclassical approximation schemes are…
The connections between standard theoretical tools used to study open quantum systems can sometimes seem opaque. Whether it is a Lindblad master equation, the equation of motion for the Wigner function or a dissipative Keldysh action,…
Ideal gases can be modeled by the Boltzmann equation from statistical physics. Instead of trying to track the position and velocity of a large number of gas molecules, it is possible to describe the particles with a particle distribution…
In this paper, we present an efficient fully spectral approximation scheme for exploring the one-dimensional steady-state neutron transport equation. Our methodology integrates the spectral-(Petrov-)Galerkin scheme in the spatial dimension…
In this paper, we study the biharmonic equation with the Navier boundary conditions in a polygonal domain. In particular, we propose a method that effectively decouples the 4th-order problem into a system of Poisson equations. Different…
We present a method for constructing numerical schemes with up to 3rd strong convergence order for solution of a class of stochastic differential equations, including equations of the Langevin type. The construction proceeds in two stages.…
Simulation of conditional master equations is important to describe systems under continuous measurement and for the design of control strategies in quantum systems. For large bosonic systems, such as BEC and atom lasers, full quantum field…
We present a variationally separable splitting technique for the generalized-$\alpha$ method for solving parabolic partial differential equations. We develop a technique for a tensor-product mesh which results in a solver with a linear cost…
We report on a computational approach based on the self-consistent solution of the steady-state Boltzmann transport equation coupled with the Poisson equation for the study of inhomogeneous transport in deep submicron semiconductor…
We present an explicit solver of the three-dimensional screened and unscreened Poisson's equation which combines accuracy, computational efficiency and versatility. The solver, based on a mixed plane-wave / interpolating scaling function…
The present work proposes a second-order time splitting scheme for a linear dispersive equation with a variable advection coefficient subject to transparent boundary conditions. For its spatial discretization, a dual Petrov--Galerkin method…
Well-conditioned spectral collocation and spectral methods have recently been proposed to solve differential equations. In this paper, we revisit the well-conditioned spectral collocation methods proposed in [T.~A. Driscoll, {\it J. Comput.…
Time evolution of the expectation values of various dynamical operators of the harmonic oscillator with dissipation is analitically obtained within the framework of the Lindblad's theory for open quantum systems. We deduce the density…
We solve the two-point function of the lowest dimensional scalar operator in the critical $\phi^4$ theory on $4-\epsilon$ dimensional real projective space in three different methods. The first is to use the conventional perturbation…
Here we present exact, stationary, parametric solutions to the Schr\"odinger--Poisson system. We confront two images: on one hand, we draw on the homotopy analysis method which leads us to a nonlinear integral scheme. Indeed, this approach…
In approximating solutions of nonstationary problems, various approaches are used to compute the solution at a new time level from a number of simpler (sub-)problems. Among these approaches are splitting methods. Standard splitting schemes…
The discrete spectra of certain two-dimensional Schrodinger operators are numerically calculated. These operators have interesting spectral properties, i.e. their kernels are multi-dimensional and the deformations of potentials via the…