Related papers: Convection in a Rapidly Rotating Spherical Shell: …
We are interested here in describing the linear response of the ocean to some wind forcing, which admits fast time oscillations and may be resonant with the Coriolis force. In addition to the usual Ekman layer, we exhibit another - much…
In this study, we use numerical simulations to investigate the flow field induced by a single magnetic microrobot rotating with a constant angular speed about an axis perpendicular to an underlying surface. A parallel solver for steady…
A weakly magnetized neutron star (NS) undergoing disk accretion should release about a half of its power in a compact region known as the accretion boundary layer. Latitudinal spread of the accreted matter and efficient radiative cooling…
Turbulent convection in the interiors of the Sun and the Earth occurs at high Rayleigh numbers $Ra$, low Prandtl numbers $Pr$, and different levels of rotation rates. To understand the combined effects better, we study rotating turbulent…
In planetary fluid cores, the density depends on temperature and chemical composition, which diffuse at very different rates. This leads to various instabilities, bearing the name of double-diffusive convection. We investigate rotating…
Stability is an important aspect of numerical methods for hyperbolic conservation laws and has received much interest. However, continuity in time is often assumed and only semidiscrete stability is studied. Thus, it is interesting to…
A conservative primitive variable discrete exterior calculus (DEC) discretization of the Navier-Stokes equations is performed. An existing DEC method (Mohamed, M. S., Hirani, A. N., Samtaney, R. (2016). Discrete exterior calculus…
We propose a semismooth Newton method for non-Newtonian models of incompressible flow where the constitutive relation between the shear stress and the symmetric velocity gradient is given implicitly; this class of constitutive relations…
We investigate the properties of forced inertial modes of a rotating fluid inside a spherical shell. Our forcing is tidal like, but its main property is that it is on the large scales. Our solutions first confirm some analytical results…
Accurate simulations of ice sheet dynamics, mantle convection, lava flow, and other highly viscous free-surface flows involve solving the coupled Stokes/free-surface equations. In this paper, we theoretically analyze the stability and…
We describe a new algorithm to solve the time dependent, frequency integrated radiation transport (RT) equation implicitly, which is coupled to an explicit solver for equations of magnetohydrodynamics (MHD) using {\sf Athena++}. The…
We investigate the linear onset of thermal convection in rotating spherical shells with a focus on the influence of mechanical boundary conditions and thermal driving modes. Using a spectral method, we determine critical Rayleigh numbers,…
In close exoplanetary systems, tidal interactions drive orbital and spin evolution of planets and stars over long timescales. Tidally-forced inertial waves (restored by the Coriolis acceleration) in the convective envelopes of low-mass…
We present a simple modification of the direct-forcing immersed boundary method (IBM) proposed by Uhlmann [J. Comput. Phys, 2005] in order to enable it to be applied to particulate flows with solid-to-fluid density ratios around unity. The…
Simulating the interaction of fluids with immersed moving solids is playing an important role for gaining a better quantitative understanding of how fluid dynamics is altered by the presence of obstacles and which forces are exerted on the…
Kraichnan seminal ideas on inverse cascades yielded new tools to study common phenomena in geophysical turbulent flows. In the atmosphere and the oceans, rotation and stratification result in a flow that can be approximated as…
We present a generalized form of open boundary conditions, and an associated numerical algorithm, for simulating incompressible flows involving open or outflow boundaries. The generalized form represents a family of open boundary…
Numerical schemes for the solution of the Euler equations have recently been developed, which involve the discretisation of the internal energy equation, with corrective terms to ensure the correct capture of shocks, and, more generally,…
An efficient, iterative semi-implicit (SI) numerical method for the time integration of stiff wave systems is presented. Physics-based assumptions are used to derive a convergent iterative formulation of the SI scheme which enables the…
This paper proposes a simple new closure principle for turbulent shear flows. The turbulent flow field is divided into an outer and an inner region. The inner region is made up of a log-law region and a wall layer. The wall layer is viewed…