Related papers: pizza: An open-source pseudo-spectral code for sph…
We compare 1D nonlocal turbulent convection models with 3D hydrodynamic numerical simulations. We study the validity of closure models and turbulent coefficients by varying the Prandtl number, the P$\acute{e}$clet number, and the depth of…
This work introduces a new higher-order super-compact (HOSC) implicit finite difference scheme for analyzing three-dimensional (3D) natural convection and entropy generation in non-Newtonian fluids. The proposed scheme achieves fourth-order…
Reservoir computing is applied to model the large-scale evolution and the resulting low-order turbulence statistics of a two-dimensional turbulent Rayleigh-B\'{e}nard convection flow at a Rayleigh number ${\rm Ra}=10^7$ and a Prandtl number…
In this paper, a novel stochastic extra-step quasi-Newton method is developed to solve a class of nonsmooth nonconvex composite optimization problems. We assume that the gradient of the smooth part of the objective function can only be…
We describe a newly developed hydrodynamic code for studying accretion disk processes. The numerical method uses a finite volume, nonlinear, Total Variation Diminishing (TVD) scheme to capture shocks and control spurious oscillations. It is…
We describe a spectral method for the numerical solution of the Vlasov-Poisson system where the velocity space is decomposed by means of an Hermite basis, and the configuration space is discretized via a Fourier decomposition. The novelty…
In this work, we consider a Shallow-Water Quasi Geostrophic equation on the sphere, as a model for global large-scale atmospheric dynamics. This equation, previously studied by Verkley (2009) and Schubert et al. (2009), possesses a rich…
The objective of this study is to develop a fully compressible magnetohydrodynamic solver for fast simulations of the global dynamo of the Sun using unstructured grids and GPUs. Accurate modeling of the Sun's convective layers is vital to…
In this work we propose a new algorithm for the computation of statistical equilibrium quantities on a cubic lattice when both an energy and a statistical temperature are involved. We demonstrate that the pivot algorithm used in situations…
Many time-dependent partial differential equations (PDEs) can be transformed into an ordinary differential equations (ODEs) containing moderately stiff and non-stiff terms after spatial semi-discretization. In the present paper, we…
This paper presents a unified approach to the modeling and computation of the Kubo conductivity of incommensurate bilayer heterostructures at finite temperature. Firstly, we derive an expression for the large-body limit of Kubo-Greenwood…
In this thesis we study the evolution of systems of concentric shells interacting gravitationally and in the process (1) propose and implement a nearly energy-conserving numerical integration scheme for evolving the concentric spherical…
As we approach solar convection simulations that seek to model the interaction of small-scale granulation and supergranulation and even larger scales of convection within the near-surface shear layer (NSSL), the treatment of the boundary…
This paper continues our study of quasicrystals initiated in Part I. We propose a general mechanism for constructing quasicrystals, existing globally in time, in spatially-extended systems (partial differential equations with Euclidean…
The hard-sphere model is one of the most extensively studied models in statistical physics. It describes the continuous distribution of spherical particles, governed by hard-core interactions. An important quantity of this model is the…
This paper studies the dissipative generalized surface quasi-geostrophic equations in a supercritical regime where the order of the dissipation is small relative to order of the velocity, and the velocities are less regular than the…
Double-diffusive convection, often referred to as semi-convection in astrophysics, occurs in thermally and compositionally stratified systems which are stable according to the Ledoux-criterion but unstable according to the Schwarzchild…
The quasistatic normal-compliance contact problem of isotropic homogeneous linear visco-elastic bodies with Coulomb friction at small strains in Kelvin-Voigt rheology is considered. The discretization is made by a semi-implicit formula in…
(ABRIDGED) We present a new Schwarzschild orbit-superposition code designed to model discrete datasets composed of velocities of individual kinematic tracers in a dynamical system. This constitutes an extension of previous implementations…
We study the implementation of a Chebyshev spectral method with forward Euler integrator to investigate a peridynamic nonlocal formulation of Richards' equation. We prove the convergence of the fully-discretization of the model showing the…