Source author record

Piotr K. Smolarkiewicz

Piotr K. Smolarkiewicz appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

3works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

3 published item(s)

preprint2021arXiv

Stratified flow past a sphere at moderate Reynolds numbers

A numerical study of stably stratified flows past spheres at Reynolds numbers $Re=200$ and $Re=300$ is reported. In these flow regimes, a neutrally stratified laminar flow induces distinctly different near-wake features. However, the flow behaviour changes significantly as the stratification increases and suppresses the scale of vertical displacements of fluid parcels. Computations for a range of Froude numbers $Fr\in [0.1,\infty]$ show that as Froude number decreases, the flow patterns for both Reynolds numbers become similar. The representative simulations of the lee-wave instability at $Fr=0.625$ and the two-dimensional vortex shedding at $Fr=0.25$ regimes are illustrated for flows past single and tandem spheres, thereby providing further insight into the dynamics of stratified flows past bluff bodies. In particular, the reported study examines the relative influence of viscosity and stratification on the dividing streamline elevation, wake structure and flow separation. The solutions of the Navier-Stokes equations in the incompressible Boussinesq limit are obtained on unstructured meshes suitable for simulations involving multiple bodies. Computations are accomplished using the finite volume, non-oscillatory forward-in-time (NFT) Multidimensional Positive Definite Transport Algorithm (MPDATA) based solver. The impact and validity of the numerical approximations, especially for the cases exhibiting strong stratification, are also discussed. Qualitative and quantitative comparisons with available laboratory experiments and prior numerical studies confirm the validity of the numerical approach.

preprint2020arXiv

Evolution of three-dimensional coherent structures in Hall magnetohydrodynamics

The work extends the computational model EULAG-MHD to include Hall magnetohydrodynamics (HMHD)---important to explore physical systems undergoing fast magnetic reconnection at the order of the ion inertial length scale. Examples include solar transients along with reconnections in magnetosphere, magnetotail and laboratory plasmas. The paper documents the results of two distinct sets of implicit large-eddy simulations in the presence and absence of the Hall forcing term, initiated with an unidirectional sinusoidal magnetic field. The HMHD simulation while benchmarking the code also emphasizes the complexity of three dimensional (3D) evolution over its two dimensional (2D) counterpart. The magnetic reconnections onset significantly earlier in HMHD. Importantly, the magnetic field generated by the Hall term breaks any inherent symmetry, ultimately making the evolution 3D. The resulting 3D reconnections develop magnetic flux ropes and magnetic flux tubes. Projected on the reconnection plane, the ropes and tubes appear as magnetic islands, which later break into secondary islands, and finally coalesce to generate an X-type neutral point. These findings are in agreement with the theory and contemporary simulations of HMHD, and thus verify our extension of the EULAG-MHD model. The second set explores the influence of the Hall forcing on generation and ascend of a magnetic flux rope from sheared magnetic arcades---a novel scenario instructive in understanding the coronal transients. The rope evolves through intermediate complex structures, ultimately breaking locally because of reconnections. Interestingly, the breakage occurs earlier in the presence of the Hall term, signifying faster dynamics leading to magnetic topology favorable for reconnections.

preprint2014arXiv

libmpdata++ 0.1: a library of parallel MPDATA solvers for systems of generalised transport equations

This paper accompanies first release of libmpdata++, a C++ library implementing the Multidimensional Positive-Definite Advection Transport Algorithm (MPDATA). The library offers basic numerical solvers for systems of generalised transport equations. The solvers are forward-in-time, conservative and non-linearly stable. The libmpdata++ library covers the basic second-order-accurate formulation of MPDATA, its third-order variant, the infinite-gauge option for variable-sign fields and a flux-corrected transport extension to guarantee non-oscillatory solutions. The library is equipped with a non-symmetric variational elliptic solver for implicit evaluation of pressure gradient terms. All solvers offer parallelisation through domain decomposition using shared-memory parallelisation. The paper describes the library programming interface, and serves as a user guide. Supported options are illustrated with benchmarks discussed in the MPDATA literature. Benchmark descriptions include code snippets as well as quantitative representations of simulation results. Examples of applications include: homogeneous transport in one, two and three dimensions in Cartesian and spherical domains; shallow-water system compared with analytical solution (originally derived for a 2D case); and a buoyant convection problem in an incompressible Boussinesq fluid with interfacial instability. All the examples are implemented out of the library tree. Regardless of the differences in the problem dimensionality, right-hand-side terms, boundary conditions and parallelisation approach, all the examples use the same unmodified library, which is a key goal of libmpdata++ design. The libmpdata++ library is implemented in C++, making use of the Blitz++ multi-dimensioanl array containers, and is released as free/libre and open-source software.