Related papers: A solver based on pseudo-spectral analytical time-…
A novel time domain solver of Maxwell's equations in passive (dispersive and absorbing) media is proposed. The method is based on the path integral formalism of quantum theory and entails the use of ({\it i}) the Hamiltonian formalism and…
In this work, we implement a new London equation module for superconductivity in the GPU-enabled ARTEMIS framework, and couple it to a finite-difference time-domain solver for Maxwell's equations. We apply this two-fluid approach to model a…
Relativistic plasmas are central to the study of black hole accretion, jet physics, neutron star mergers, and compact object magnetospheres. Despite the need to accurately capture the dynamics of these plasmas and the implications for…
This work presents a unified fluid modeling framework for reacting flows coupled with nonthermal plasmas (NTPs). Building upon the gas-plasma kinetics solver, ChemPlasKin, and the CFD library, OpenFOAM, the integrated solver, reactPlasFOAM,…
We determine conditions allowing to simplify the description of the impact of a short and arbitrarily intense laser pulse onto a cold plasma at rest. If both the initial plasma density and pulse profile have plane simmetry, then suitable…
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…
In this article, we propose high-order finite-difference entropy stable schemes for the two-fluid relativistic plasma flow equations. This is achieved by exploiting the structure of the equations, which consists of three independent flux…
In this paper, we present the numerical solution of two-phase flow problems of engineering significance with a space-time finite element method that allows for local temporal refinement. Our basis is the method presented in [3], which…
Collisionless plasmas, mostly present in astrophysical and space environments, often require a kinetic treatment as given by the Vlasov equation. Unfortunately, the six-dimensional Vlasov equation can only be solved on very small parts of…
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…
The self-consistent modeling of vacuum polarization due to virtual electron-positron fluctuations is of relevance for many near term experiments associated with high intensity radiation sources and represents a milestone in describing…
Resolvent analysis is a powerful tool for studying coherent structures in turbulent flows. However, its application beyond canonical flows with symmetries that can be used to simplify the problem to inherently three-dimensional flows and…
The time domain linear sampling method (TD-LSM) solves inverse scattering problems using time domain data by creating an indicator function for the support of the unknown scatterer. It involves only solving a linear integral equation called…
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 present a new and efficient phase-field solver for viscoelastic fluids with moving contact line based on a dual-resolution strategy. The interface between two immiscible fluids is tracked by using the Cahn-Hilliard phase-field model, and…
This paper is concerned with the large-time behavior of solutions to the Cauchy problem on the two-fluid Euler-Maxwell system with collisions when initial data are around a constant equilibrium state. The main goal is the rigorous…
Validity of fluid models breaks down for non-thermal or weakly collisional plasmas which often occur e.g. in the solar wind. In these regimes one has to resort to modelling through the first-principle Vlasov-Maxwell system, but its…
Previously developed method for finding asymptotic solutions of Vlasov equations using two-dimensional (in coordinate x and time t) Laplace transform is applied to low-collision electron-ion plasmas. Taking into account Coulomb collisions…
We present a novel multiscale numerical approach that combines parallel-in-time computation with hybrid domain adaptation for linear collisional kinetic equations in the diffusive regime. The method addresses the computational challenges of…
An application of approximate transformation groups to study dynamics of a system with distinct time scales is discussed. The utilization of the Krylov-Bogoliubov-Mitropolsky method of averaging to find solutions of the Lie equations is…