Related papers: A pencil distributed finite difference code for st…
The use of reduced and mixed precision computing has gained increasing attention in high-performance computing (HPC) as a means to improve computational efficiency, particularly on modern hardware architectures like GPUs. In this work, we…
An efficient algorithm for computing the terms appearing in the Generalised Kolmogorov Equation (GKE) written for the indefinite plane channel flow is presented. The algorithm, which features three distinct strategies for parallel…
A high-order scheme for direct numerical simulations of turbulent combustion is discussed. Its implementation in the massively parallel and publicly available Pencil Code is validated with the focus on hydrogen combustion. Ignition delay…
In order to prevent velocity, pressure, and temperature spikes at material discontinuities occurring when the interface-capturing schemes inconsistently simulate compressible multi-material flows(when the specific heats ratio is…
Preserving scalar boundedness is important for numerical schemes used in turbulent compressible multi-component flow simulations to prevent unphysical results and unstable simulations. However, ensuring scalar boundedness for high-order,…
A new scheme that tightly couples kinetic turbulence codes across a spatial interface is introduced. This scheme evolves from considerations of competing strategies and down-selection. It is found that the use of a composite kinetic…
Intermittent turbulent-laminar patterns characterize the transition to turbulence in pipe, plane Couette and plane channel flows. The time evolution of turbulent-laminar bands in plane channel flow is studied via direct numerical…
A new slow growth formulation for DNS of wall-bounded turbulent flow is developed and demonstrated to enable extension of slow growth modeling concepts to complex boundary layer flows. As in previous slow growth approaches, the formulation…
A low-Mach-number flow, in the laminar regime, has intrinsically two characteristic spatial scales for a given time scale, or two characteristic temporal scales for a given spatial scale, and these dual scales are very different due to the…
In wall-bounded flows, the laminar regime remain linearly stable up to large values of the Reynolds number while competing with nonlinear turbulent solutions issued from finite amplitude perturbations. The transition to turbulence of plane…
In this paper, we present the results of a numerical study of air-water turbulent bubbly flow in a periodic vertical square duct. The study is conducted using a novel numerical technique which leverages Volume of Fluid method for interface…
We have conducted a direct numerical simulation (DNS) study of dilute turbulent particulate flow in a vertical plane channel, considering thousands of finite-size rigid particles with resolved phase interfaces. The particle diameter…
Numerical schemes used for the integration of complex flow simulations should provide accurate solutions for the long time integrations these flows require. To this end, the performance of various high-order accurate numerical schemes is…
We focus on implementing and optimizing a sixth-order finite-difference solver for simulating compressible fluids on a GPU using third-order Runge-Kutta integration. Since graphics processing units perform well in data-parallel tasks, this…
In particle-laden turbulent flows the turbulence in carrier fluid phase gets affected by the dispersed particle phase for volume fraction above $10^{-4}$ and hence reverse coupling or two-way coupling becomes relevant in that volume…
More and more massive parallel codes running on several hundreds of thousands of cores enter the computational science and engineering domain, allowing high-fidelity computations on up to trillions of unknowns for very detailed analyses of…
Wall-pressure fluctuations are a practically robust input for real-time control systems aimed at modifying wall-bounded turbulence. The scaling behaviour of the wall-pressure--velocity coupling requires investigation to properly design a…
We present an improved method for computing incompressible viscous flow around suspended rigid particles using a fixed and uniform computational grid. The main idea is to incorporate Peskin's regularized delta function approach [Acta…
Turbulent wall flows offer the most direct means for understanding the effects of boundaries and viscosity on turbulent fluctuations. Available data on mean-square fluctuations in these flows show apparent contradiction with classical…
A low diffusive flux difference splitting based kinetic scheme is developed based on a discrete velocity Boltzmann equation, with a novel three velocity model. While two discrete velocities are used for upwinding, the third discrete…