Related papers: Relativistic Magnetohydrodynamics: Renormalized ei…
The Velocity-Vorticity (VV) formulation of the incompressible Navier-Stokes equations has become popular in recent years, especially in numerical studies, due to its structural advantages. Recently, with L. Rebholz, we introduced a Voigt…
The equations of general relativistic magnetohydrodynamics (GRMHD) have become the standard mathematical framework for modeling high-energy plasmas in curved spacetimes. However, the fragility of the primitive variable reconstruction…
The paper describes a new upwind conservative numerical scheme for special relativistic resistive magnetohydrodynamics with scalar resistivity. The magnetic field is kept approximately divergence free and the divergence of the electric…
In this work, we develop efficient solvers for linear inverse problems based on randomized singular value decomposition (RSVD). This is achieved by combining RSVD with classical regularization methods, e.g., truncated singular value…
Generalising the Elsasser variables, we introduce the Q-variables. These are more flexible than the Elsasser variables, because they also allow to track waves with phase speeds different than the Alfven speed. We rewrite the MHD equations…
The evolution of electromagnetic and thermodynamic fields in a non-ideal fluid are studied in the framework of ultrarelativistic transverse magnetohydrodynamics (MHD), which is essentially characterized by electric and magnetic fields being…
Magnetohydrodynamics in divergence form describes a hyperbolic system of covariant and constraint-free equations. It comprises a linear combination of an algebraic constraint and Faraday's equations. Here, we study the problem of…
Ideal systems like MHD and Euler flow may develop singularities in vorticity (w = curl v). Viscosity and resistivity are dissipative regularizations. We propose a minimal, local, conservative, nonlinear, dispersive regularization of…
We report on recent extensions and improvements to the Legolas code, which is an open-source, finite element-based numerical framework to solve the linearised (magneto)hydrodynamic equations for a three-dimensional force- and thermally…
This paper explores the problem of generalized phase retrieval, which involves reconstructing a length-$n$ signal $\bm{x}$ from its $m$ phaseless samples $y_k = \left|\langle \bm{a}_k,\bm{x}\rangle\right|^2$, where $k = 1,2,...,m$, and…
In this paper, we are concerned with the two-dimensional (2D) incompressible magnetohydrodynamic (MHD) equations with velocity dissipation given by $(-\Delta)^{\alpha}$ and magnetic diffusion given by reducing about logarithmic diffusion…
We consider a general family of regularized Navier-Stokes and Magnetohydrodynamics (MHD) models on n-dimensional smooth compact Riemannian manifolds with or without boundary, with n greater than or equal to 2. This family captures most of…
We prove existence, uniqueness, and higher-order global regularity of strong solutions to a particular Voigt-regularization of the three-dimensional inviscid resistive Magnetohydrodynamic (MHD) equations. Specifically, the coupling of a…
We derive generic relativistic hydrodynamical equations with dissipative effects from the underlying Boltzmann equation in a mechanical and systematic way on the basis of so called the renormalization-group (RG) method. A macroscopic frame…
A high-order flux reconstruction solver has been developed and validated to perform implicit large-eddy simulations of industrially representative turbomachinery flows. The T106c low-pressure turbine and VKI LS89 high-pressure turbine cases…
This work presents the first numerical investigation of using Voigt regularization as a method for obtaining magnetohydrodynamic (MHD) equilibria without the assumption of nested magnetic flux surfaces. Voigt regularization modifies the MHD…
We present a numerical method to solve the equations of general relativistic hydrodynamics in a given external gravitational field. The method is based on a generalization of Roe's approximate Riemann solver for the non relativistic Euler…
We address the estimation of seismic wavefields by means of Multidimensional Deconvolution (MDD) for various redatuming applications. While offering more accuracy than conventional correlation-based redatuming methods, MDD faces challenges…
In this chapter, we aim at presenting the basic techniques necessary to go beyond the widely accepted paradigm of second-order numerics. We specifically focus on finite-volume schemes for hyperbolic conservation laws occuring in fluid…
A deflated restarted Lanczos algorithm is given for both solving symmetric linear equations and computing eigenvalues and eigenvectors. The restarting limits the storage so that finding eigenvectors is practical. Meanwhile, the deflating…