相关论文: GAMERA-OP: A three-dimensional finite-volume MHD s…

This paper develops the genuinely multidimensional HLL Riemann solver for the two-dimensional special relativistic hydrodynamic equations on Cartesian meshes and studies its physical-constraint preserving (PCP) property. Based on the…

Many astrophysical and terrestrial scenarios involving magnetic fields can be approached in axial geometry. Although the smoothed particle hydrodynamics (SPH) technique has been successfully extended to magneto-hydrodynamics (MHD), a…

This paper puts forth a coarse grid projection (CGP) multiscale method to accelerate computations of quasigeostrophic (QG) models for large scale ocean circulation. These models require solving an elliptic sub-problem at each time step,…

A novel orthogonalization-free method together with two specific algorithms are proposed to solve extreme eigenvalue problems. On top of gradient-based algorithms, the proposed algorithms modify the multi-column gradient such that earlier…

In axisymmetric fusion reactors, the equilibrium magnetic configuration can be expressed in terms of the solution to a semi-linear elliptic equation known as the Grad-Shafranov equation, the solution of which determines the poloidal…

In this work, we take a modern high-resolution finite-volume scheme for solving the rotational shallow-water equations and extend it with features required to run real-world ocean simulations. Our contributions include a spatially varying…

This paper explores the effects of numerical algorithms on global magnetohydrodynamics (MHD) simulations of solar wind (SW) in the inner heliosphere. To do so, we use sunRunner3D, a 3-D MHD model that employs the boundary conditions…

Three dimensional (3D) topology optimization problems always involve huge numbers of Degrees of Freedom (DOFs) in finite element analysis (FEA) and design variables in numerical optimization, respectively. This will inevitably lead to large…

We extend the general-relativistic magnetohydrodynamics (GRMHD) capabilities of Athena++ to incorporate radiation. The intensity field in each finite-volume cell is discretized in angle, with explicit transport in both space and angle…

We consider deterministic Markov decision processes (MDPs) and apply max-plus algebra tools to approximate the value iteration algorithm by a smaller-dimensional iteration based on a representation on dictionaries of value functions. The…

Due to the prevalence of magnetic fields in astrophysical environments, magnetohydrodynamic (MHD) simulation has become a basic tool for studying astrophysical fluid dynamics. To further advance the precision of MHD simulations, we have…

We introduce a max-plus analogue of the Petrov-Galerkin finite element method, to solve finite horizon deterministic optimal control problems. The method relies on a max-plus variational formulation, and exploits the properties of…

This paper presents the implementation of a HLLC finite volume solver using GPU technology for the solution of shallow water problems in two dimensions. It compares both CPU and GPU approaches for implementing all the solver's steps. The…

Magnetohydrodynamics (MHD) describes the interaction between electrically conducting fluids and electromagnetic fields. We propose and analyze a symplectic, second-order algorithm for the evolutionary MHD system in Els\"asser variables. We…

We describe a magnetohydrodynamic (MHD) constrained energy functional for equilibrium calculations that combines the topological constraints of ideal MHD with elements of Taylor relaxation. Extremizing states allow for partially chaotic…

Much recent research effort has been directed to the development of efficient algorithms for solving minimax problems with theoretical convergence guarantees due to the relevance of these problems to a few emergent applications. In this…

The metriplectic framework, which permits to formulate an algebraic structure for dissipative systems, is applied to visco-resistive Magneto-Hydrodynamics (MHD), adapting what had already been done for non-ideal Hydrodynamics (HD). The…

Complex-variable matrix optimization problems (CMOPs) in Frobenius norm emerge in many areas of applied mathematics and engineering applications. In this letter, we focus on solving CMOPs by iterative methods. For unconstrained CMOPs, we…

A new magnetohydrodynamics (MHD) code based on initial value approach, GMEC_I, has been developed for simulating various MHD physics in tokamak plasmas, as the MHD foundation of the gyrokinetic-MHD energetic particle simulation code (GMEC)…

Polyhedral Omega is a new algorithm for solving linear Diophantine systems (LDS), i.e., for computing a multivariate rational function representation of the set of all non-negative integer solutions to a system of linear equations and…