相关论文: A comparison of spectral element and finite differ…
The incompressible Hall-magnetohydrodynamics (Hall--MHD) system presents substantial analytical and computational challenges due to its stiff, highly nonlinear Hall term and the strict requirement that the magnetic field remains solenoidal.…
Recently, secondary island formation due to the tearing instability of the Sweet-Parker current sheet was identified as a possible mechanism that can lead to fast reconnection (less sensitive dependence on Lundquist number $S$) both in…
We present an adaptive scheme for isogeometric phase-field modeling, to perform suitably graded hierarchical refinement and coarsening on both single- and multi-patch geometries by considering truncated hierarchical spline constructions…
We study dendritic microstructure evolution using an adaptive grid, finite element method applied to a phase-field model. The computational complexity of our algorithm, per unit time, scales linearly with system size, rather than the…
I provide a pedagogic review of adaptive mesh refinement (AMR) radiation hydrodynamics (RHD) methods and codes used in simulations of star formation, at a level suitable for researchers who are not computational experts. I begin with a…
Many problems in physics are inherently of multi-scale nature. The issues of MHD turbulence or magnetic reconnection, namely in the hot and sparse, almost collision-less astrophysical plasmas, can stand as clear examples. The Finite Element…
An adaptive finite element method is presented for the elastic scattering of a time-harmonic plane wave by a periodic surface. First, the unbounded physical domain is truncated into a bounded computational domain by introducing the…
Three algebraically stabilized finite element schemes for discretizing convection-diffusion-reaction equations are studied on adaptively refined grids. These schemes are the algebraic flux correction (AFC) scheme with Kuzmin limiter, the…
We report a new finite-difference scheme for two-dimensional magnetohydrodynamics (MHD) simulations, with and without rotation, in unstructured grids with quadrilateral cells. The new scheme is implemented within the code VULCAN/2D, which…
In this work we present a proof of concept of CUDA-capable, resistive, multi-fluid models of relativistic magnetohydrodynamics (RMHD). Resistive and multi-fluid codes for simulating models of RMHD suffer from stiff source terms, so it is…
Fully adaptive computations of the resistive magnetohydrodynamic (MHD) equations are presented in two and three space dimensions using a finite volume discretization on locally refined dyadic grids. Divergence cleaning is used to control…
We develop a micromorphic-based approach for finite element stabilization of reaction-convection-diffusion equations, by gradient enhancement of the field of interest via introducing an auxiliary variable. The well-posedness of the…
We report the development of a discontinuous spectral element flow solver that includes the implementation of both spectral difference and flux reconstruction formulations. With this high order framework, we have constructed a foundation…
A 2D incompressible viscoresistive-MHD model [Mahapatra et. al., Phys. Plasmas 28, 072103 (2021)] is used to study the scaling of reconnection parameters in magnetic island coalescence problem under two interesting scenarios. Firstly, the…
Spectral element methods (SEM), which are extensions of finite element methods (FEM), are important emerging techniques for solving partial differential equations in physics and engineering. SEM can potentially deliver better accuracy due…
A mixed mimetic spectral element method is applied to solve the rotating shallow water equations. The mixed method uses the recently developed spectral element histopolation functions, which exactly satisfy the fundamental theorem of…
Melting and solidification processes are often affected by natural convection of the liquid, posing a multi-physics problem involving fluid flow, convective and diffusive heat transfer, and phase-change reactions. Enthalpy methods formulate…
The purpose of this research is to describe an efficient iterative method suitable for obtaining high accuracy solutions to high frequency time-harmonic scattering problems. The method allows for both refinement of local polynomial degree…
A two-grid scheme based on mixed finite-element approximations to the incompressible Navier-Stokes equations is introduced and analyzed. In the first level the standard mixed finite-element approximation over a coarse mesh is computed. In…
When a numerical simulation has to handle a physics problem with a wide range of time-dependent length scales, dynamically adaptive discretizations can be the method of choice. We present a major upgrade to the numerical relativity code…