Addressing the challenges of implementation of high-order finite volume schemes for atmospheric dynamics of unstructured meshes

The solution of the non-hydrostatic compressible Euler equations using Weighted Essentially Non-Oscillatory (WENO) schemes in two and three-dimensional unstructured meshes, is presented. Their key characteristics are their simplicity; accuracy; robustness; non-oscillatory properties; versatility in handling any type of grid topology; computational and parallel efficiency. Their defining characteristic is a non-linear combination of a series of high-order reconstruction polynomials arising from a series of reconstruction stencils. In the present study an explicit TVD Runge-Kutta 3rd -order method is employed due to its lower computational resources requirement compared to implicit type time advancement methods. The WENO schemes (up to 5th -order) are applied to the two dimensional and three dimensional test cases: a 2D rising

Higher-order CFD and Interface Tracking Methods on Highly-Parallel MPI and GPU systems

A computational investigation of the effects on parallel performance of higher-order accurate schemes was carried out on two different computational systems: a traditional CPU based MPI cluster and a system of four Graphics Processing Units (GPUs) controlled by a single quad-core CPU. The investigation was based on the solution of the level set equations for interface tracking using a High-Order Upstream Central (HOUC) scheme. Different variants of the HOUC scheme were employed together with a 3rd-order TVD Runge-Kutta time integration. An increase in performance of two orders of magnitude was seen when comparing a single CPU core to a single GPU with a greater increase at higher orders of accuracy and at lower precision

Validation of a magneto- and ferro-hydrodynamic model for non-isothermal flows in conjunction with Newtonian and non-Newtonian fluids

This work focuses on the validation of a magnetohydrodynamic (MHD) and ferrohydrodynamic (FHD) model for non-isothermal flows in conjunction with Newtonian and non- Newtonian fluids. The importance of this research field is to gain insight into the interaction of non-linear viscous behaviour of blood flow in the presence of MHD and FHD effects, because its biomedical application such as magneto resonance imaging (MRI) is in the centre of research interest. For incompressible flows coupled with MHD and FHD models, the Lorentz force and a Joule heating term appear due to the MHD effects and the magnetization and magnetocaloric terms appear due to the FHD effects in the non-linear momentum and temperature equations, respectively. Tzirtzilakis and Loukopoulos [1] investigated the effects of MHD and FHD for incompressible non-isothermal flows in conjunction with Newtonian fluids in a small rectangular channel. Their model excluded the non-linear viscous behaviour of blood flows considering blood as a Newtonian biofluid. Tzirakis et al. [2, 3] modelled the effects of MHD and FHD for incompressible isothermal flows in a circular duct and through a stenosis in conjunction with both Newtonian and non-Newtonian fluids, although their approach neglects the non-isothermal magnetocaloric FHD effects. Due to the fact that there is a lack of experimental data available for non-isothermal and non-Newtonian blood flows in the presence of MHD and FHD effects, therefore the objective of this study is to establish adequate validation test cases in order to assess the reliability of the implemented non-isothermal and non-Newtonian MHD-FHD models. The non-isothermal Hartmann flow has been chosen as a benchmark physical problem to study velocity and temperature distributions for Newtonian fluids and non-Newtonian blood flows in a planar microfluidic channel. In addition to this, the numerical behaviour of an incompressible and non-isothermal non-Newtonian blood flow has been investigated from computational aspects when a dipole-like rotational magnetic field generated by infinite conducting wires. The numerical results are compared to available computational data taken from literature

WENO schemes on arbitrary mixed-element unstructured meshes in three space dimensions

The paper extends weighted essentially non-oscillatory (WENO) methods to three dimensional mixed-element unstructured meshes, comprising tetrahedral, hexahedral, prismatic and pyramidal elements. Numerical results illustrate the convergence rates and non-oscillatory properties of the schemes for various smooth and discontinuous solutions test cases and the compressible Euler equations on various types of grids. Schemes of up to fifth order of spatial accuracy are considered

Growth of a Richtmyer-Meshkov turbulent layer after reshock

This paper presents a numerical study of a reshocked turbulent mixing layer using high-order accurate Implicit Large-Eddy-Simulations (ILES). Existing theoretical approaches are discussed, and the theory of Youngs (detailed in Ref. 1) is extended to predict the behaviour of a reshocked mixing layer formed initially from a shock interacting with a broadband instability. The theory of Mikaelian2 is also extended to account for molecular mixing in the single-shocked layer prior to reshock. Simulations are conducted for broadband and narrowband initial perturbations and results for the growth rate of the reshocked layer and the decay rate of turbulent kinetic energy show excellent agreement with the extended theoretical approach. Reshock causes a marginal decrease in mixing parameters for the narrowband layer, but a significant increase for the broadband initial perturbation. The layer properties are observed to be very similar post-reshock, however, the growth rate exponent for the mixing layer width is higher in the broadband case, indicating that the reshocked layer still has a dependence (although weakened) on the initial conditions. These results have important implications for Unsteady Reynolds Averaged Navier Stokes modelling of such instabilities

Multiscale simulation strategies and mesoscale modelling of gas and liquid flows

This paper was presented at the 2nd Micro and Nano Flows Conference (MNF2009), which was held at Brunel University, West London, UK. The conference was organised by Brunel University and supported by the Institution of Mechanical Engineers, IPEM, the Italian Union of Thermofluid dynamics, the Process Intensification Network, HEXAG - the Heat Exchange Action Group and the Institute of Mathematics and its Applications.This paper presents a review of multiscale simulation strategies for the modelling of micro- and nanoscale flows. These have been developed in the last two decades in an attempt to bridge the application gap between molecular and continuum simulation methods preventing the simulation of many micro- and nanofluidic devices. The paper is focused on hybrid molecular-continuum methods and reviews different coupling strategies, including geometrical decomposition in conjunction with state- and flux coupling, pointwise coupling, the heterogeneous multiscale method and the equation free approach. The different applications of these methods are briefly discussed

Advances and challenges in computational research of micro and nano flows

This paper presents an overview of past and current research in computational modelling of micro- and nanofluidic systems with particular focus on recent advances in multiscale modelling. Different mesoscale and hybrid molecular-continuum methods are presented. The contributions of these methods to a broad range of applications, as well as the physical and computational modelling challenges associated with the development of these methods, are also discussed

Thermodynamics at solid-liquid interfaces

The variation of the liquid properties in the vicinity of a solid surface complicates the description of heat transfer along solid-liquid interfaces. Using Molecular Dynamics simulations, this investigation aims to understand how the material properties, particularly the strength of the solid-liquid interaction, affect the thermal conductivity of the liquid at the interface. The molecular model consists of liquid argon confined by two parallel, smooth, solid walls, separated by a distance of 6.58σ. We find that the component of the thermal conductivity parallel to the surface increases with the affinity of the solid and liquid

Computational uncertainty in hybrid atomistic-continuum frameworks

This paper was presented at the 3rd Micro and Nano Flows Conference (MNF2011), which was held at the Makedonia Palace Hotel, Thessaloniki in Greece. The conference was organised by Brunel University and supported by the Italian Union of Thermofluiddynamics, Aristotle University of Thessaloniki, University of Thessaly, IPEM, the Process Intensification Network, the Institution of Mechanical Engineers, the Heat Transfer Society, HEXAG - the Heat Exchange Action Group, and the Energy Institute.Over the past decade micro and nanofluidics emerged as vital tools in the ongoing drive towards the development of nano-scale analysis and manufacturing systems. Accurate numerical modelling of the phenomena involved at these scales is ssential in order to speed up the industrial design process for nanotechnology. However a parameter often ignored in hybrid simulations is the uncertainty level introduced in the numerical modelling of phenomena taking place at micro and nanoscales. The main interest of the present study is the propagation of the inherent atomistic fluctuations to the continuum solver in the case of multiscale modelling and hybrid solvers

Implicit large eddy simulation of weakly-compressible turbulent channel flow

This paper concerns the accuracy of several high-resolution and high-order finite volume schemes in Implicit Large Eddy Simulation of weakly-compressible turbulent channel flow. The main objective is to investigate the properties of numerical schemes, originally designed for compressible flows, in low Mach compressible, near-wall turbulent flows. Variants of the Monotone Upstream-centred Scheme for Conservation Laws and Weighted Essentially Non-Oscillatory schemes for orders of accuracy ranging from second to ninth order, as well as with and without low Mach corrections, have been investigated. The performance of the schemes has been assessed against incompressible Direct Numerical Simulations. Detailed comparisons of the velocity profiles, turbulent shear stresses and higher-order turbulent statistics reveal that the low Mach correction can significantly reduce the numerical dissipation of the methods in low Mach boundary layer flows. The effects of the low Mach correction have more profound impact on second and third-order schemes, but they also improve the accuracy of fifth order schemes. The ninth-order Weighted Essentially Non-Oscillatory scheme is the least dissipative scheme and it is shown that the implementation of the low Mach correction in conjunction with this scheme has a significant anti-dissipative effect that adversely affects the accuracy. Finally, the computational cost required for obtaining the improved accuracy using increasingly higher order schemes is also discussed