Related papers: Multigrid iterative algorithm using pseudo-compres…
High order algorithms have emerged in numerical astrophysics as a promising avenue to reduce truncation error (proportional to a power of the linear resolution $\Delta x$) with only a moderate increase to computational expense. Significant…
We have devised a variational sinc collocation method (VSCM) which can be used to obtain accurate numerical solutions to many strong-coupling problems. Sinc functions with an optimal grid spacing are used to solve the linear and non-linear…
Two methods for solid body representation in flow simulations available in the Pencil Code are the immersed boundary method and overset grids. These methods are quite different in terms of computational cost, flexibility and numerical…
This paper solves the discretised multiphase flow equations using tools and methods from machine-learning libraries. The idea comes from the observation that convolutional layers can be used to express a discretisation as a neural network…
We present a multigrid algorithm for the solution of the linear systems of equations stemming from the $p-$version of the Virtual Element discretization of a two-dimensional Poisson problem. The sequence of coarse spaces are constructed…
We present a three-dimensional (3D) common-refinement method for non-matching meshes between discrete non-overlapping subdomains of incompressible fluid and nonlinear hyperelastic structure. To begin, we first investigate the accuracy of…
We present a novel high order semi-implicit hybrid finite volume/virtual element numerical scheme for the solution of compressible flows on Voronoi tessellations. The method relies on the flux splitting of the compressible Navier-Stokes…
To design a method to solve the issues of handling 'dirty' and highly complex geometries, the topology-free method combined with the immersed boundary method is presented for viscous and incompressible flows at a high Reynolds number. The…
This paper introduces a novel approach to compute the numerical fluxes at the cell boundaries in the finite volume approach. Explicit gradients used in deriving the reconstruction polynomials are replaced by high-order gradients computed by…
Preserving scalar boundedness is important for numerical schemes used in turbulent compressible multi-component flow simulations to prevent unphysical results and unstable simulations. However, ensuring scalar boundedness for high-order,…
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…
Conventional mathematical models for simulating incompressible fluid flow problems are based on the Navier-Stokes equations expressed in terms of pressure and velocity. In this context, pressure-velocity coupling is a key issue, and…
The multigrid algorithm is a multilevel approach to accelerate the numerical solution of discretized differential equations in physical problems involving long-range interactions. Multiresolution analysis of wavelet theory provides an…
In this paper, a two-dimensional incompressible miscible displacement model is considered, and a novel decoupled and linearized high-order finite difference scheme is developed, by utilizing the multi-time-step strategy to treat the…
We propose and analyze a new method for the unsteady incompressible magnetohydrodynamics equations on convex domains with hybrid approximations of both vector-valued and scalar-valued fields. The proposed method is convection-semirobust,…
As a phase space language for quantum mechanics, the Wigner function approach bears a close analogy to classical mechanics and has been drawing growing attention, especially in simulating quantum many-body systems. However, deterministic…
The use of multigrid and related preconditioners with the finite element method is often limited by the difficulty of applying the algorithm effectively to a problem, especially when the domain has a complex shape or adaptive refinement. We…
In this work, we introduce a Variational Multi-Scale (VMS) method for the numerical approximation of parabolic problems, where sub-grid scales are approximated from the eigenpairs of associated elliptic operator. The abstract method is…
This article performs a unified convergence analysis of a variety of numerical methods for a model of the miscible displacement of one incompressible fluid by another through a porous medium. The unified analysis is enabled through the…
Well designed lattice-Boltzmann codes exploit the essentially embarrassingly parallel features of the algorithm and so can be run with considerable efficiency on modern supercomputers. Such scalable codes permit us to simulate the behaviour…