Related papers: A solver based on pseudo-spectral analytical time-…
In this work, we suggest an easy-to-code higher-order finite volume semi-discrete scheme to analyze the nonlinear behavior of the electron-plasma oscillations by solving electron fluid equations numerically. The present method employs a…
We report a numerical simulation for the phase diagram of a simple two dimensional model, similar to one proposed by Noro and Frenkel [J. Chem. Phys. \textbf{114}, 2477 (2001)] for membrane proteins, but one that includes the role of the…
The scattered field formalism is combined to the particle-in-cell method to model relativistic laser-plasma dynamics in complex field configurations. Despite the strong nonlinearity of the interactions, we demonstrate the validity of this…
An oblique pulsar magnetosphere with a plasma conductivity is studied by using a pseudo-spectral method. In the pseudo-spectral method, the time-dependent Maxwell equations are solved, both electric and magnetic fields are expanded in terms…
The Gas-Kinetic Scheme (GKS), widely used in computational fluid dynamics for simulating hypersonic and other complicated flow phenomena, is extended in this work to electromagnetic problems by solving Maxwell's equations. In contrast to…
We focus on non-stationary Maxwell equations defined on a regular patch of elements as considered in the isogeometric analysis (IGA). We apply the time-integration scheme following the ideas developed by the finite difference community [M.…
We present a 3D finite element solver for the nonlinear Poisson-Nernst-Planck (PNP) equations for electrodiffusion, coupled to the Stokes system of fluid dynamics. The model serves as a building block for the simulation of macromolecule…
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,…
An improved numerical solver for the unified solution of compressible and incompressible fluids involving interfaces is proposed. The present method is based on the CIP-CUP (Cubic Interpolated Propagation / Combined, Unified Procedure)…
We use our vector Maxwell's nonlinear eigenmode solver to study the stationary solutions in 2D cross-section waveguides with Kerr nonlinear cores. This solver is based on the fixed power algorithm within the finite element method. First,…
We derive and interpret solutions of time-harmonic Maxwell's equations with a vertical and a horizontal electric dipole near a planar, thin conducting film, e.g. graphene sheet, lying between two unbounded isotropic and non-magnetic media.…
Pseudospectral numerical schemes for solving the Dirac equation in general static curved space are derived using a pseudodifferential representation of the Dirac equation along with a simple Fourier-basis technique. Owing to the presence of…
Active Flux (AF) is a modified Finite Volume method that evolves additional Degrees of Freedom (DoF) located on the cell interfaces to compute high-order approximations to the numerical fluxes through the respective interface. We present an…
We derive novel algorithms for optimization problems constrained by partial differential equations describing multiscale particle dynamics, including non-local integral terms representing interactions between particles. In particular, we…
A long-standing challenge encountered in modeling plasma dynamics is achieving practical Vlasov equation simulation in multiple spatial dimensions over large length and time scales. While direct multi-dimension Vlasov simulation methods…
We develop a spectrally accurate numerical method to compute solutions of a model partial differential equation used in plasma physics to describe diffusion in velocity space due to Fokker-Planck collisions. The solution is represented as a…
Ultrahigh peak power femtosecond laser pulses create extreme states of matter that are currently being probed with great interest. Plasma optics have been proposed for shaping and amplifying high-power pulses, but they are subject to huge…
We investigate the performance of a hybrid plasma solver on the test problem of an ion beam. The parallel solver is based on cell centered finite differences in space, and a predictor-corrector leapfrog scheme in time. The implementation is…
We propose an adaptive numerical solver for the study of viscoelastic 2D two-phase flows using the volume-of-fluid method. The scheme uses the robust log conformation tensor technique of Fattal & Kupferman (2004,2005} combined with the…
We investigate the dynamics of localized solutions of the relativistic cold fluid plasma model in the small but finite amplitude limit, for slightly overcritical plasma density. Adopting a multiple scale analysis, we derive a perturbed…