Related papers: ExoCubed: A Riemann-Solver based Cubed-Sphere Dyna…

We present a new general relativistic (GR) code for hydrodynamic supernova simulations with neutrino transport in spherical and azimuthal symmetry (1D/2D). The code is a combination of the CoCoNuT hydro module, which is a Riemann-solver…

Computational modeling and simulation of fluid-structure interactions constitute a fundamental cornerstone for advancing aerospace engineering endeavors. This paper addresses the notion and implementation of the immersed boundary method for…

We describe a method to extract resonant orbits from N-body simulations exploiting the fact that they close in a frame rotating with a constant pattern speed. Our method is applied to the N-body simulation of the Milky Way by Shen et al.…

Kilometer-scale Earth system models are essential for capturing local climate change. However, these models are computationally expensive and produce petabyte-scale outputs, which limits their utility for applications such as probabilistic…

An algorithm for simulating self-gravitating cosmological astrophysical fluids is presented. The advantages include a large dynamic range, parallelizability, high resolution per grid element and fast execution speed. The code is based on a…

We present a computational method to solve the magnetohydrodynamic equations in spherical geometry. The technique is fully nonlinear and wholly spectral, and uses an expansion basis that is adapted to the geometry: Chandrasekhar-Kendall…

We use a general circulation model to study the three-dimensional (3-D) flow and temperature distributions of atmospheres on tidally synchronized extrasolar planets. In this work, we focus on the sensitivity of the evolution to the initial…

Climate models are often affected by long-term drift that is revealed by the evolution of global variables such as the ocean temperature or the surface air temperature. This spurious trend reduces the fidelity to initial conditions and has…

We present a one-equation subgrid scale model that evolves the turbulence energy corresponding to unresolved velocity fluctuations in large eddy simulations. The model is derived in the context of the Germano consistent decomposition of the…

Salmon's nearly geostrophic model for rotating shallow-water flow is derived in full spherical geometry. The model, which results upon constraining the velocity field to the height field in Hamilton's principle for rotating shallow-water…

Shallow water surface flows commonly entrain sediments, resulting in scouring and/or deposition of the underlying substrate that may strongly influence the pattern of subsequent flow. These coupled phenomena, which can be investigated…

The inner edge of the classical habitable zone is often defined by the critical flux needed to trigger the runaway greenhouse instability. This 1D notion of a critical flux, however, may not be so relevant for inhomogeneously irradiated…

We present a new particle-based (discrete element) numerical method for the simulation of granular dynamics, with application to motions of particles on small solar system body and planetary surfaces. The method employs the parallel N-body…

In this work, we consider the generalized variable-coefficient nonlinear Schr\"{o}dinger equation with non-vanishing boundary conditions at infinity including the simple and double poles of the scattering coefficients. By introducing an…

In recent years, molecular dynamics (MD) simulations have emerged as a pivotal tool for understanding the structure, dynamics, and phase behavior in charged soft matter systems. To explore phenomena across greater length and time scales in…

Writing high performance solvers for engineering applications is a delicate task. These codes are often developed on an application to application basis, highly optimized to solve a certain problem. Here, we present our work on developing a…

Semi-implicit time-stepping schemes for atmosphere and ocean models require elliptic solvers that work efficiently on modern supercomputers. This paper reports our study of the potential computational savings when using mixed precision…

We describe an implicit general relativistic hydrodynamics code. The evolution equations are formulated in comoving coordinates. A conservative finite differencing of the Einstein equations is outlined, and artificial viscosity and…

In this paper, we present an analytical study, in the one space dimensional case, of the fluid dynamics system proposed in [4] to model the formation of biofilms. After showing the hyperbolicity of the system, we show that, in a open…

We investigate through numerical simulations the hydrodynamic interactions between two rigid spherical particles suspended on the axis of a cylindrical tube filled with an elastoviscoplastic fluid subjected to pressure-driven flow. The…