Related papers: Arbitrarily high-order energy-conserving methods f…
We consider the problem of ''energy conserving'' lattice Boltzmann models. A major difficulty observed in previous studies is the coupling between the viscous and thermal waves even at moderate wave numbers. We propose a theoretical…
We propose and study a class of arbitrarily high-order numerical discretizations that preserve multiple invariants and are essentially explicit (they do not require the solution of any large systems of algebraic equations). In space, we use…
We consider a class of multiscale parabolic problems with diffusion coefficients oscillating in space at a possibly small scale $\varepsilon$. Numerical homogenization methods are popular for such problems, because they capture efficiently…
We study relaxation-based approaches for conserving mass and energy in the numerical solution of Schr\"odinger-Poisson (SP) type systems. Relaxation-based methods offer a general approach that can be applied as post-time step processing to…
I recently proposed a method of bosonization based on the use of coherent states of fermion composites, whose validity was restricted to smooth structure functions. In the present paper I remove this limitation and derive results which hold…
A high-order accurate adjoint-based optimization framework is presented for unsteady multiphysics problems. The fully discrete adjoint solver relies on the high-order, linearly stable, partitioned solver introduced in [1], where different…
In this article, using a known method, a computation is performed of the derivatives of the microcanonical entropy, with respect to the energy up to the 4-th order, using a Laplace transform technique, and adapted it to the case where the…
Using a recent characterization of energy-preserving B-series, we derive the explicit conditions on the coefficients of a Runge-Kutta method that ensure energy preservation (for Hamiltonian systems) up to a given order in the step size,…
This work is devoted to the numerical simulation of a Vlasov-Poisson model describing a charged particle beam under the action of a rapidly oscillating external electric field. We construct an Asymptotic Preserving numerical scheme for this…
We report a few sumerical tests comparing some newly defined energy-preserving integrators and symplectic methods, using either constant and variable stepsize.
We investigate the existence of solutions of reversible and irreversible port-Hamiltonian systems. To this end, we utilize the associated exergy, a function that is composed of the system's Hamiltonian and entropy, to prove global existence…
We consider stopping criteria that balance algebraic and discretization errors for the conjugate gradient algorithm applied to high-order finite element discretizations of Poisson problems. Firstly, we introduce a new stopping criterion…
By using the variational minimizing method with a special constraint and the direct variational minimizing method without constraint, we study second order Hamiltonian systems with a singular potential $V\in C^2(R^n\backslash O,R)$ and…
We construct a finite element discretization and time-stepping scheme for the incompressible Euler equations with variable density that exactly preserves total mass, total squared density, total energy, and pointwise incompressibility. The…
In this work we design and analyze a free energy satisfying finite difference method for solving Poisson-Nernst-Planck equations in a bounded domain. The algorithm is of second order in space, with numerical solutions satisfying all three…
In this paper, high order well-balanced finite difference weighted essentially non-oscillatory methods to solve general systems of balance laws are presented. Two different families are introduced: while the methods in the first one…
This work presents a novel formulation and numerical strategy for the simulation of geometrically nonlinear structures. First, a non-canonical Hamiltonian (Poisson) formulation is introduced by including the dynamics of the stress tensor.…
In this work, we introduce a novel algorithm for the Biot problem based on a Hybrid High-Order discretization of the mechanics and a Symmetric Weighted Interior Penalty discretization of the flow. The method has several assets, including,…
A recent reformulation [1] of the problem of Coulomb gases in the presence of a dynamical dielectric medium showed that finite temperature simulations of such systems can be accomplished on the basis of completely local Hamiltonians on a…
A new class of high-order maximum principle preserving numerical methods is proposed for solving parabolic equations, with application to the semilinear Allen--Cahn equation. The proposed method consists of a $k$th-order multistep…