Related papers: A fully implicit, conservative, non-linear, electr…

The movement of ionic solutions is an essential part of biology and technology. Fluidics, from nano- to micro- to microfluidics, is a burgeoning area of technology which is all about the movement of ionic solutions, on various scales. Many…

This work focuses on the development of a self adjusting multirate strategy based on an implicit time discretization for the numerical solution of hyperbolic equations, that could benefit from different time steps in different areas of the…

Modeling multi-scale collisionless magnetized processes constitutes an important numerical challenge. By treating electrons as a fluid and ions kinetically, the so-called hybrid Particle-In-Cell (PIC) codes represent a promising…

The motion of a fully ionized plasma of electrons and ions is generally governed by the Vlasov-Maxwell-Landau system. We prove the global existence of solutions near Maxwellians to the Cauchy problem of the system for the long-range…

We present a conservative/dissipative time integration scheme for nonlinear mechanical systems. Starting from a weak form, we derive algorithmic forces and velocities that guarantee the desired conservation/dissipation properties. Our…

We construct a mean-field model that describes the nonlinear dynamics of a spin-polarized electron gas interacting with fixed, positively-charged ions possessing a magnetic moment that evolves in time. The mobile electrons are modeled by a…

Astrophysical fluids under the influence of magnetic fields are often subjected to single-fluid or two-fluid approximations. In the case of weakly ionized plasmas however, this can be inappropriate due to distinct responses from the…

We propose a model of ion-electron plasma (or nucleus-electron plasma) that accounts for the electronic structure around nuclei (i.e. ion structure) as well as for ion-ion correlations. The model equations are obtained through the…

A new one-dimensional fluid model for ions in weakly-ionized plasma is proposed. The model differs from the existing ones in two aspects. First, a more accurate approximation of the collision terms in the fluid equations is suggested. For…

A fundamental two-fluid model for describing dynamics of a plasma is the Euler-Poisson system, in which compressible ion and electron fluids interact with their self-consistent electrostatic force. Global smooth electron dynamics were…

In this paper, an implicit nonsymplectic exact energy-preserving integrator is specifically designed for a ten-dimensional phase-space conservative Hamiltonian system with five degrees of freedom. It is based on a suitable…

Ion transport, often described by the Poisson--Nernst--Planck (PNP) equations, is ubiquitous in electrochemical devices and many biological processes of significance. In this work, we develop conservative, positivity-preserving, energy…

Two-fluid ideal plasma equations are a generalized form of the ideal MHD equations in which electrons and ions are considered as separate species. The design of efficient numerical schemes for the these equations is complicated on account…

This study established a quantum-classical hybrid framework that integrates quantum computing paradigm with meshfree finite particle method. By harnessing quantum superposition and entanglement, it hybridized the critical computational…

High-resolution simulations of particle-based kinetic plasma models typically require a high number of particles and thus often become computationally intractable. This is exacerbated in multi-query simulations, where the problem depends on…

We present a broad family of high-order finite element algorithms for simulating the flow of electroneutral electrolytes. The governing partial differential equations that we solve are the electroneutral…

We propose an unconditionally energy-stable, orthonormality-preserving, component-wise splitting iterative scheme for the Kohn-Sham gradient flow based model in the electronic structure calculation. We first study the scheme discretized in…

Bounded plasmas are characterized by a rapid but smooth transition from quasi-neutrality in the volume to electron depletion close to the electrodes and chamber walls. The thin non-neutral region, the boundary sheath, comprises only a small…

Hybrid particle-field methods are computationally efficient approaches for modelling soft matter systems. So far applications of these methodologies have been limited to constant volume conditions. Here, we reformulate particle-field…

We propose Coadjoint Orbit FLIP (CO-FLIP), a high order accurate, structure preserving fluid simulation method in the hybrid Eulerian-Lagrangian framework. We start with a Hamiltonian formulation of the incompressible Euler Equations, and…