Related papers: A Numerical Technique for Coupling the Momentum an…
We present an exponentially convergent semi-implicit meshless algorithm for the solution of Navier-Stokes equations in complex domains. The algorithm discretizes partial derivatives at scattered points using radial basis functions as…
In this paper, we propose and analyze a fully discrete finite element projection method for the magnetohydrodynamic (MHD) equations. A modified Crank--Nicolson method and the Galerkin finite element method are used to discretize the model…
Dissipation can be represented in Hamiltonian mechanics in an extended phase space as a symplectic process. The method uses an auxiliary variable which represents the excitation of unresolved dynamics and a Hamiltonian for the interaction…
Numerical simulation of compressible fluid flows is performed using the Euler equations. They include the scalar advection equation for the density, the vector advection equation for the velocity and a given pressure dependence on the…
In this work I develop a new framework for anisotropic hydrodynamics that generalizes the leading order of the hydrodynamic expansion to the full (3+1)-dimensional anisotropic massive case. Following previous works, my considerations are…
We revise the structure-preserving finite element method in [K. Hu, Y. MA and J. Xu. (2017) Stable finite element methods preserving $\nabla \cdot \mathbf{B}=0$ exactly for MHD models. Numer. Math.,135, 371-396]. The revised method is…
In this paper, we present and analyze a new finite difference method for computing three dimensional wave maps into spheres. By introducing the angular momentum as an auxiliary variable, we recast the governing equation as a first order…
We propose a new unstructured numerical subgrid method for solving the shallow water equations using a finite volume method with enhanced bathymetry resolution. The method employs an unstructured triangular mesh with support for…
A numerical method is presented to compute the eigenmodes supported by three dimensional (3D) metamaterials using the Method of Moments (MoM). The method relies on interstitial equivalent currents between layers. First, a parabolic…
To simulate the dynamics of fluid with polydisperse particles on macroscale level, one has to solve hydrodynamic equations with several relaxation terms, representing momentum transfer from fluid to particles and vice versa. For small…
In this paper, we study a nonlinear system of first order partial differential equations describing the macroscopic behavior of an ensemble of interacting self-propelled rigid bodies. Such system may be relevant for the modelling of bird…
A fully implicit finite difference scheme has been developed to solve the hydrodynamic equations coupled with radiation transport. Solution of the time dependent radiation transport equation is obtained using the discrete ordinates method…
In the framework of a mixed finite element method, a structure-preserving formulation for incompressible magnetohydrodynamic (MHD) equations with general boundary conditions is proposed. A leapfrog-type temporal scheme fully decouples the…
In multi-phase fluid flow, fluid-structure interaction, and other applications, partial differential equations (PDEs) often arise with discontinuous coefficients and singular sources (e.g., Dirac delta functions). These complexities arise…
In this paper we consider the spatial semi-discretization of conservative PDEs. Such finite dimensional approximations of infinite dimensional dynamical systems can be described as flows in suitable matrix spaces, which in turn leads to the…
Active Flux is a third order accurate numerical method which evolves cell averages and point values at cell interfaces independently. It naturally uses a continuous reconstruction, but is stable when applied to hyperbolic problems. In this…
One dimensional (1D) simulations of the flow and flooding of open channels are known to be inaccurate as the flow is multi-dimensional in nature, especially at the flooded regions. However, multi-dimensional simulations, even in two…
We consider the design of structure-preserving discretization methods for the solution of systems of boundary controlled Partial Differential Equations (PDEs) thanks to the port-Hamiltonian formalism. We first provide a novel general…
We study three different time integration methods for a dynamic pore network model for immiscible two-phase flow in porous media. Considered are two explicit methods, the forward Euler and midpoint methods, and a new semi-implicit method…
In this paper we propose a stable and robust strategy to approximate the 3d incompressible hydrostatic Euler and Navier-Stokes systems with free surface. Compared to shallow water approximation of the Navier-Stokes system, the idea is to…