Related papers: A simple solver for the two-fluid plasma model bas…
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…
Implicit particle-in-cell codes offer advantages over their explicit counterparts in that they suffer weaker stability constraints on the need to resolve the higher frequency modes of the system. This feature may prove particularly valuable…
In this paper, a perfectly matched layer (PML) method is proposed to solve the time-domain electromagnetic scattering problems in 3D effectively. The PML problem is defined in a spherical layer and derived by using the Laplace transform and…
Rigorous computer simulations of propagating electromagnetic fields have become an important tool for optical metrology and optics design of nanostructured components. As has been shown in previous benchmarks some of the presently used…
Simulation of unsteady creeping flows in complex geometries has traditionally required the use of a time-stepping procedure, which is typically costly and unscalable. To reduce the cost and allow for computations at much larger scales, we…
With the stagnation of processor core performance, further reductions in the time-to-solution for geophysical fluid problems are becoming increasingly difficult with standard time integrators. Parallel-in-time exposes and exploits…
We present an explicit multiscale algorithm for solving differential equations for problems with high-frequency modes that can be averaged over by separating and scaling the fast and slow dynamics within a single equation. We introduce a…
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…
This paper is concerned with boundary integral equation methods for solving the two-dimensional fluid-solid interaction problem. We reduce the problem to three differential systems of boundary integral equations via direct and indirect…
This paper is concerned with global-in-time, nonoverlapping domain decomposition methods for the mixed formulation of the diffusion problem. Two approaches are considered: one uses the time-dependent Steklov-Poincar\'e operator and the…
In this work, a one-dimensional simulation code was developed for both single-phase and two-phase systems, focusing on time-dependent Euler equations for gas and particles. These equations, non-linear hyperbolic conservation laws, describe…
Electromagnetic particle-in-cell (PIC) codes are widely used to perform computer simulations of a variety of physical systems, including fusion plasmas, astrophysical plasmas, plasma wakefield particle accelerators, and secondary photon…
We present a parareal in time algorithm for the simulation of neutron diffusion transient model. The method is made efficient by means of a coarse solver defined with large time steps and steady control rods model. Using finite element for…
Efficient simulation of the Navier-Stokes equations for fluid flow is a long standing problem in applied mathematics, for which state-of-the-art methods require large compute resources. In this work, we propose a data-driven approach that…
In this paper we present a technique for constructing robust solvers for stiff algebraic source terms, such as those typically used for modelling relaxation processes in hyperbolic systems of partial differential equations describing…
In this paper, we present a model based on a local thermodynamic equilibrium, weakly ionized plasma-mixture model used for medical and technical applications in etching processes. We consider a simplified model based on the Maxwell-Stefan…
Within OpenFOAM, we develop a pressure-based solver for the Euler equations written in conservative form using density, momentum, and total energy as variables. Under simplifying assumptions, these equations are used to describe…
We construct smooth, non-symmetric plasma equilibria which possess closed, nested flux surfaces and solve the magnetohydrostatic (steady three-dimensional incompressible Euler) equations with a small force. The solutions are also `nearly'…
Pseudospectral time domain (PSTD) methods are widely used in many branches of acoustics for the numerical solution of the wave equation, including biomedical ultrasound and seismology. The use of the Fourier collocation spectral method in…
In this work, we introduce PHOENIX, a highly optimized explicit open-source solver for two-dimensional nonlinear Schr\"odinger equations with extensions. The nonlinear Schr\"odinger equation and its extensions (Gross-Pitaevskii equation)…