Related papers: Application of an Upwind Integration Method to Pla…
We propose an efficient algorithm for the evaluation of the potential and its gradient of gravitational/electrostatic $N$-body systems, which we call particle mesh multipole method (PMMM or PM$^3$). PMMM can be understood both as an…
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…
Fundamental solution of a space fractional convection equation of order $\alpha$ is the probability density function of L\'{e}vy flights with long-tailed $\alpha$-stable jump length distribution. By studying an upwind second-order implicit…
We study the role of the Hall current and electron inertia in collisionless magnetic reconnection within the framework of full two-fluid MHD. At spatial scales smaller than the electron inertial length, a topological change of magnetic…
Arc plasmas have promising applications in many fields. To explore their property is of interest. This paper presents detailed pressure-based finite volume simulation of argon arc. In the modeling, the whole cathode region is coupled to…
Much of the progress in our understanding of dynamo mechanisms has been made within the theoretical framework of magnetohydrodynamics (MHD). However, for sufficiently diffuse media, the Hall effect eventually becomes non-negligible. We…
During relativistic magnetic reconnection, antiparallel magnetic fields undergo a rapid change in topology, releasing a large amount of energy in the form of non-thermal particle acceleration. This work explores the application of mesh…
Colliding hypersonic flows play a decisive role in many astrophysical objects. In this paper, we look at the idealized model of a 2D plane parallel isothermal slab (CDL) and at symmetric settings, where both flows have equal parameters. We…
In an idealized system where four current channels interact in a two-dimensional periodic setting, we follow the detailed evolution of current sheets (CSs) forming in between the channels, as a result of a large-scale merging. A central…
Numerical geodynamo simulations with parameters close to an Earth-like regime would be of great interest for understanding the dynamics of the Earth's liquid outer core and the associated geomagnetic field. Such simulations are far too…
Weakly ionised plasmas, formed in high enthalpy hypersonic flows, can be actively manipulated via imposed magnetic fields - a concept termed magnetohydrodynamic (MHD) flow control. Imposed MHD effects, within flows which exhibit multiple…
In the finite difference method which is commonly used in computational micromagnetics, the demagnetizing field is usually computed as a convolution of the magnetization vector field with the demagnetizing tensor that describes the…
In this work we explore the evolution of magnetic fields inside strongly magnetized neutron stars in axisymmetry. We model numerically the coupled field evolution in the core and the crust. Our code models the Hall drift and Ohmic effects…
Steady state simulations} of magnetized electron fluid equations with strong anisotropic diffusion based on the first-order hyperbolic approach is carried out using cell-centered higher order upwind schemes, linear and weighted essentially…
Strongly correlated materials are a natural target for fault-tolerant quantum computers, but they require tools beyond those developed for molecules. Electronic wavefunctions vary rapidly near nuclei yet remain delocalized across many unit…
We present a simulation capability for micro-scale light-emitting diodes (uLEDs) that achieves comparable accuracy to CPU-based finite-difference time-domain simulation but is more than 10^7 times faster. Our approach is based on the…
We present new simulations of the evolution of axially symmetric magnetic fields in neutron star crusts under the influence of the Hall effect and subdominant Ohmic dissipation. In the Hall effect, differential rotation of the electron…
Finite element methods usually construct basis functions and quadrature rules for multidimensional domains via tensor products of one-dimensional counterparts. While straightforward, this approach results in integration spaces larger than…
Numerical methods to improve the treatment of magnetic fields in smoothed field magnetohydrodynamics (SPMHD) are developed and tested. Chapter 2 is a review of SPMHD. In Chapter 3, a mixed hyperbolic/parabolic scheme is developed which…
Accurate prediction of rarefied gas flows is important for space vehicle design, particularly in rarefied regimes where the Navier-Stokes equations are no more valid. While the direct simulation Monte Carlo (DSMC) method acts as a numerical…