Related papers: Legolas 2.0: Improvements and extensions to an MHD…
Scarcity of hydrocarbon resources and high exploration risks motivate the development of high fidelity algorithms and computationally viable approaches to exploratory geophysics. Whereas early approaches considered least-squares…
An optimisation scheme is developed to accurately represent the sub-grid scale forcing of a high dimensional chaotic ocean system. Using a simple parameterisation scheme, the velocity components of a 30km resolution shallow water ocean…
We present a new 1-D multi-physics simulation code with use cases intended for, but not limited to, hydrodynamic escapeproblems of planetary atmospheres and planetary accretion models. Our formulation treats an arbitrary number of species…
Two-dimensional free-surface flow over localised topography is examined with the emphasis on the stability of hydraulic-fall solutions. A Gaussian topography profile is assumed with a positive or negative amplitude modelling a bump or a…
We modify an existing magnetohydrodynamics algorithm to make it more compatible with a dimensionally-split (DS) framework. It is based on the standard reconstruct-solve-average strategy (using a Riemann solver), and relies on constrained…
Axisymmetric magnetohydrodynamics (MHD) can be invoked for describing astrophysical magnetized flows and formulated to model stellar magnetospheres including main sequence stars (e.g. the Sun), compact stellar objects [e.g. magnetic white…
We present a multi-dimensional numerical code to solve isothermal magnetohydrodynamic (IMHD) equations for use in modeling astrophysical flows. First, we have built a one-dimensional code which is based on an explicit finite-difference…
We present the extension of the differentiable hydrodynamics code, diffhydro, enabling scalable PDE-constrained inference and integrated hybrid physics-ML models for a wide range of astrophysical applications. New physics additions include…
We have developed a three-dimensional (3D) spectral hydrodynamic code to study vortex dynamics in rotating, shearing, stratified systems (e.g. the atmosphere of gas giant planets, protoplanetary disks around newly forming protostars). The…
We propose several continuous data assimilation (downscaling) algorithms based on feedback control for the 2D magnetohydrodynamic (MHD) equations. We show that for sufficiently large choices of the control parameter and resolution and…
We describe a numerical method to solve the magnetohydrodynamic (MHD) equations. The fluid variables are updated along each direction using the flux conservative, 2nd order, total variation diminishing (TVD), upwind scheme of Jin and Xin.…
A new approach is presented to compute the seismic normal modes of a fully heterogeneous, rotating planet. Special care is taken to separate out the essential spectrum in the presence of a fluid outer core. The relevant…
It is well known that the linear stability of solutions of partial differential equations which are integrable can be very efficiently investigated by means of spectral methods. We present here a direct construction of the eigenmodes of the…
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…
We present an open-source Matlab framework, titled iFluid, for simulating the dynamics of integrable models using the theory of generalized hydrodynamics (GHD). The framework provides an intuitive interface, enabling users to define and…
We present the derivation of hydrodynamical equations for a perfect fluid in General Relativity, within the 3+1 decomposition of spacetime framework, using only primitive variables. Primitive variables are opposed to conserved variables, as…
We have written and tested a new general relativistic magnetohydrodynamics (GRMHD) code, capable of evolving MHD fluids in dynamical spacetimes with adaptive-mesh refinement (AMR). Our code solves the Einstein-Maxwell-MHD system of coupled…
A hybrid spectral/finite-element code is developed to numerically solve the resistive finite-pressure magnetohydrodynamic equilibria without the necessity of postulating nested magnetic flux surfaces in the non-axisymmetric toroidal…
Numerical solutions of the incompressible magnetohydrodynamic (MHD) equations are reported for the interior of a rotating, perfectly-conducting, rigid spherical shell that is insulator-coated on the inside. A previously-reported spectral…
We have developed a new global eigenvalue code, Multiscale Analysis for plasma Stabilities (MAS), for studying plasma problems with wave toroidal mode number n and frequency omega in a broad range of interest in general tokamak geometry,…