(2021). Industrial Applications of Smoothed Particle Hydrodynamics. International Journal of Computational Fluid Dynamics: Vol. 35, Industrial Applications of Smoothed Particle Hydrodynamics. Guest Editors: Nathan J. Quinlan and Mingming Tong, pp. 1-2.

The aim of this paper is the presentation of a new modified weighted essentially non-oscillatory scheme with a seventh-order of accuracy, and, a novel nonlinear weighting method. Local smoothness indicators of this scheme constructed using Lagrange’s interpolation polynomial. We developed a new high-order global smoothness indicator, which guarantees the optimal order of accuracy of the scheme at critical

ABSTRACT For simulation of hypersonic turbulent and transitional flows at high Reynolds numbers, a gas-kinetic scheme (GKS) strongly coupled with the turbulent kinetic energy equation in shear stress transport (SST) k−ω turbulence model is developed. To extend the current method to transitional flows, Langtry-Menter γ−Reθ transition model is added together with SST k−ω turbulence model. Comparisons

This paper presents a new hybrid scheme which combines the fifth-order Targeted Essentially Non-Oscillatory (TENO) reconstruction and the tangent of hyperbola for interface capturing (THINC) reconstruction by using a modified boundary variation diminishing (MBVD) algorithm to select one of the two candidate reconstructions. The THINC scheme has a jump-like distribution can approximate a discontinuity

Very low Earth orbits (VLEOs) have an altitude of 100–450 km and have become a popular topic in space technology owing to their high application potential. However, spacecraft experience stronger aerothermodynamic effects at VLEOs. In this study, we used a parallel direct simulation Monte Carlo (DSMC) code called PDSC++ and compared it with a traditional DSMC code. The results are in good agreement

A single-equation R=k/ω turbulence model is devised from the shear-stress transport (SST) k-ω closure without an assumption of equal-diffusion coefficients; evolving diffusion/destruction terms are habitually avoided owing to their numerical stiffness and complexity. Rearranging retained terms yields a production term naturally with an elliptic identity in the R-equation and hence, the production term

In the context of the development of HALO3D (High Altitude Low Orbit 3D), an edge-based Finite Element Method multidisciplinary solver for hypersonic flows, this paper presents the simulation of flows with an imposed magnetic field. Applying the low-magnetic Reynolds number approximation, a current-continuity equation replaces Maxwell's equations with the Hall effect introduced through an electrical

Large eddy simulation is performed in heat transfer of turbulent pipe flow with NEK5000 code based on spectral element method to study uncertainty quantification by Grid Convergence Index. The turbulent pipe flow with heat transfer is considered with Reynolds number Re = 19,000 and Prandtl number Pr = 0.71. Three meshes (fine, medium, and coarse) with the fixed interpolation polynomial order, eight

To solve the problem that the attenuation effect of the bladder pressure pulsation attenuator cannot be quantitatively analysed, the dynamic mesh technology was combined with the user-defined function to calculate the serial and parallel attenuators. Overcoming the divergence problem in finite-element fluid calculation and the negative volume problem in the dynamic mesh updating process, the bladder

ABSTRACT Computational fluid dynamics (CFD) modelling was performed to simulate spatial and temporal airborne pathogen concentrations during an observed COVID-19 outbreak in a restaurant in Guangzhou, China. The reported seating configuration, overlap durations, room ventilation, layout, and dimensions were modelled in the CFD simulations to determine relative exposures and probabilities of infection

A high-fidelity, low-Mach computational fluid dynamics simulation tool that includes evaporating droplets and variable-density turbulent flow coupling is well-suited to ascertain transmission probability and supports risk mitigation methods development for airborne infectious diseases such as COVID-19. A multi-physics large-eddy simulation-based paradigm is used to explore droplet and aerosol pathogen

Although weighted compact nonlinear scheme (WCNS) has been successfully applied to the simulation of complex configurations, due to the excessive dissipation, the scheme is still not suitable for many circumstances. In this article, an optimised five-point-stencil WCNS is proposed. On the one hand, by assembling two fourth-order sub-schemes, the background upwind scheme is optimised for low dissipation

ABSTRACT This work presents a computational model for a discontinuous fibre composite manufacturing process. The alignment mechanism of this novel process, called the High Performance Discontinuous Fibre (HiPerDiF) method, involves highly coupled fluid–structure interactions. Fibres with a length on the order of a few millimetres are placed in a water suspension, sprayed between two parallel plates

ABSTRACT This paper addresses the influence of the scheme used in the convective terms of turbulence and transition transport equations on the numerical accuracy of transitional flow simulations. Three transition models are combined with the k − ω Shear Stress Transport turbulence model: the γ − R e θ model, the γ model and the Amplification Factor Transport model. Two airfoils and a 6:1 prolate spheroid

The vorticity confinement (VC) method was applied to tip vortices shed by edges of wings in order to predict induced drag using far-field integration. The optimal VC parameter was determined by its...

When the tire of aircraft landing gear passes through the contaminated runway, besides the side water spray, the water spray between the multi-wheels will collide and form the pattern similar to the tail of the rooster. This kind of water spray is called the cock tail flow. In this paper, numerical analysis is used to study the formation process of the cock tail flow caused by the main aircraft tires

A series of challenging fluid flow applications are used to demonstrate the powerful capabilities of the SPH method. The applications are classified according to whether they are industrial, geophysical or biophysical in nature. The versatility and flexibility of SPH allows it to be used to predict wide ranges of flow types with diverse coupled secondary physics and chemistry. The demonstration examples

One of the main issues with naive Smoothed Particle Hydrodynamics (SPH) implementations is the lack of uniform accuracy in the computational domain. If not mastered correctly, this leads to non-physical predictions when dealing with large-domain hydraulic problems or with very fine resolutions. The present article addresses two recommended methodologies to achieve the best numerical accuracy with single-precision

This paper aims to introduce the core features of the Fuzzy Logic (FL) approach in simulations of turbulent Supersonic flows as a cost-effective methodology. Computational Fluid Dynamic (CFD) analysis is conducted to study the effect of Mach number (M e ), nozzle pressure ratio (NPR) and expansion corners (EC) on primary recirculation region length (Xr1). An Open FOAM numerical solver was used for

The objective of this paper is to check the ability of the high-order OSMP scheme [Daru, V., and C. Tenaud. 2004. “HighOrder One-Step Monotonicity-Preserving Schemes for Unsteady Compressible Flow Calculations.” Journal of Computational Physics 193 (2): 563–594] to accurately compute turbulent compressible flows with a special focus on the effect of the MP constraints on solutions of wall bounded turbulent

This paper develops a new hybrid finite difference scheme which consists of a nonlinear WENOCU4 scheme to maintain the numerical stability near discontinuities and a fourth-order linear central scheme elsewhere to accurately resolve smooth fluctuations and to speed up the computation. The central scheme is constructed in a robust skew-symmetric form which satisfies energy conservation property. A new

A compact high-order gas-kinetic scheme (GKS) with spectral resolution will be presented and used in the simulation of acoustic and shock waves. Based on high-order gas evolution model, the GKS provides a time accurate solution at a cell interface, from which both time accurate flux function and flow variables can be obtained for updating cell averaged conservative flow variables and gradients. Based

A mathematical model is developed to study on the influence of magnetic field on the flow of biofluids. This paper describes the laminar, viscous incompressible flow, which is a fully developed flow of a conducting biomagnetic fluid in a rectangular duct with square cross-section, is numerically studied under the influence of magnetic field. Therefore, numerical simulation was carried out by using

Solidification is an important issue that needs to be addressed when modelling supercooled large droplets (SLD) impacting aircraft surfaces. This work develops a 3D triple-phase smoothed particle hydrodynamics (SPH) solver that can simultaneously simulate SLD impingement and solidification on cold surfaces by solving the Euler equations for fluid dynamics, the energy equation for heat transfer, and

CFD is probably amongst the numerical methods fields that benefited most from advances in parallelisation techniques. This is due to the ever increasing size and complexity of CFD applications, as ...

ABSTRACT The present study performs direct numerical simulations of turbulent channel flows using a spectral method in a large computational domain. Because of applying Fourier discretisation in the spectral method, parallelisation of the method may incur heavy communication overhead, thereby resulting in poor scalability. We design and improve the spectral code by exploring parallel techniques, including

Computational fluid dynamics (CFD) computations often require coupling of some values on one computational domain or portion thereof to another. This is the case both for multi-physics couplings, and for the implementation of various features, such as restarting on a different mesh, boundary condition couplings, and others. This usually involves locating points or elements relative to other elements

The Smoothed Particle Hydrodynamics code SPHERA (RSE SpA) is validated on a laboratory urban dam-break flood. Comparisons with a porous Shallow Water Equations - Finite Volume Method 2D model are also reported. SPHERA provides a performance comparable with the state-of-the-art code, with results closer to the measured values during the most risky flood stage. Some improvements of SPHERA predictions

(2020). Correction. International Journal of Computational Fluid Dynamics: Vol. 34, No. 6, pp. I-I.

We have examined the effects on fluid and particle motion due to solid particles passing around circular obstacles in particulate flows. Particle interaction with internal obstacles, outer boundary and with the fluid is inspected. Eulerian approach using a fixed computational mesh is used across which the solid particles move freely in fluid. Treatment of fluid and particle interaction inside the whole

The vorticity confinement (VC) method was used with total variation diminishing (TVD) schemes to reduce possible over-confinement and applied to convected and stationary vortices. The optimal VC parameter was determined by application to 2-D vortices. The conservativity of the scheme is investigated in terms of kinetic energy and momentum. Grid convergence study was conducted for the range of grid

This paper discusses selected features of the Smoothed Particle Hydrodynamics (SPH) method when used in combination with a Finite Element (FE) software tool. Special attention is paid to the implementation of wall boundary conditions in view of their relevance to the simulation of fluid-structure interaction (FSI). Several industrial applications of the hybrid SPH-FE method are reviewed in detail.

ABSTRACT For complex engineering and scientific applications, Computational Fluid Dynamics (CFD) simulations require a huge amount of computational power. As such, it is of paramount importance to carefully assess the performance of CFD codes and to study them in depth for enabling optimisation and portability. In this paper, we study three complex CFD codes, OpenFOAM, Alya and CHORUS representing

The design process is a critical issue in order to improve rocket engine performance. Indeed, the prediction of the aerodynamic forces acting on the nozzle, in particular the off-axis loads, is challenging even for modern methods. In this study, the flow of an exit Mach number Truncated Ideal Contour (TIC) nozzle is numerically investigated by means of Large-Eddy Simulation (LES), using the unstructured

Abstract This paper facilitates the integration of TRIZ and CFD techniques to resolve design challenges and drive innovation in flow handling equipment new product development processes. A new methodology is proposed based on an integrated product lifecycle suited for flow handling equipment product development processes. The methodology accomplishes novel design solutions by designing geometrical

In this paper a characteristics-based open boundary condition (CBC) is proposed for the magnetohydrodynamic (MHD) system of equations. The algorithm is carefully designed and implemented in the context of a high-order flux reconstruction (FR) scheme under the Generalised Lagrange Multiplier (GLM)-MHD system of equations. It is implemented by adding the contribution of the characteristic equation directly

This paper presents results of Artificial Neural Networks (ANN) applications to sub-grid Large Eddy Simulation (LES) model. The training data for the ANN is provided by simulation of Homogeneous Isotropic Turbulence at different Reynolds numbers. The results show that the correlation coefficients are superior to other sub-grid models, using a similar set of input variables. As the ANN model extrapolates

A new computational methodology is developed for large eddy simulation (LES) of turbulent reacting flows using graphic processing units (GPUs). The LES is based on the filtered density function (FDF) of the scalar-composition in conjunction with a discontinuous Galerkin (DG) discretisation scheme on a structured rectangular mesh. This hybrid solver is developed in a manner suitable for GPU computing

Multiphysics problems involve the couplings of different sets of partial differential equations. Partitioned methods consider the individual solutions of each set, which upon iterating, converge to the monolithic solution. The main drawback of partitioned methods is the additional iterative loop, which can be done a la Jacobi (parallel) or a la Gauss–Seidel (sequential). The latter method has worse

This paper addresses the MPI parallelisation of a 3D multiphase smoothed particle hydrodynamics (SPH) solver, focusing on supercooled large droplets (SLD) impingement. SPH uses moving particles to represent fluid flows, demanding particular parallelisation strategies. A cell system provides a spatial reference to generate dynamical neighbour lists, and to pack, send and receive particles. A re-indexing

ABSTRACT In this work, we propose an efficient ‘r-to-h’ adaptation algorithm for moving boundary problems using only vertex displacements and some local h-adaptation operations. Moving boundary algorithms usually require the regeneration of the mesh to avoid his invalidation in particular when the motion is large and/or several bodies are involved. The proposed approach aims to move the mesh nodes

Aerodynamic characteristics of tilt-rotor UAVs are usually simulated by computational fluid dynamics (CFD) only, and the fuselage/wing is fixed in simulations. It is difficult for this method to consider drastic changes of the aerodynamic forces and flowfield originating from conversion of flight modes. And during flight-mode conversions, the rapid changes of postures of tilt-rotors will unavoidably

Fault-tolerant computing options based on the use of restart information stored on and off node and the use of reserve processes have been developed, implemented and tested in a large-scale, production field solver taken from the domain of computational fluid dynamics. The tests conducted to date have shown good results, with recovery rates approaching 100% under realistic node failure scenarios. Even

Aerodynamic characteristics of tilt-rotor UAVs are usually simulated by computational fluid dynamics (CFD) only, and the fuselage/wing is fixed in simulations. It is difficult for this method to consider drastic changes of the aerodynamic forces and flowfield originating from conversion of flight modes. And during flight-mode conversions, the rapid changes of postures of tilt-rotors will unavoidably

High fidelity multiphysics simulations are ubiquitous in science and engineering but still face many challenges to run efficiently in today's supercomputers. This work reports advanced technologies present in EdgeCFD, a hybrid parallel variational multiscale multiphysics finite element software capable of running industrial problems involving turbulent incompressible, compressible, free-surface flows

Canard is one of the aerodynamic add-on devices which can reduce drag coefficient of the car. In this paper, different parameters of the canard geometry are determined using a multi-objective optimisation. Design variables are entrance velocity (U), geometrical parameters of canard (L 1, L 2, r, α) and canard angle from horizontal axis (θ). The objective functions include magnitude of drag and lift

The shark-skin-inspired riblets could effectively reduce the turbulence wall-skin friction compared with the smooth surface. Nevertheless, computations of the flow field over riblets in the practical configurations (e.g., airplanes and wind turbines) remain unfeasible for its huge computational cost. In this study, a general strategy for modelling the effects of riblets in turbulent flows is proposed

We analyse and study instability problems related to the solution of the advection–diffusion–reaction equation (ADR) using a standard finite element scheme. With this aim, this work has been carried out in the following way: first, three weak formulations are obtained from the general problem. In specific, we study the existence and uniqueness of the solution for each of the aforementioned formulations

ABSTRACT In-flight ice accretion code predictions depend on heat loss over rough surfaces. The equivalent sand grain roughness models the friction coefficient, but an additional model is needed for heat transfer predictions. In this paper, a two parameters model based on a turbulent Prandtl number correction is derived from the sublayer Stanton-based model. For flow over rough surfaces, the new model

The resolution and robustness properties of a numerical scheme are two mutually restricted aspects for the detailed simulation of complex flows. In the present paper, an improved seventh-order finite difference weighted essentially non-oscillatory scheme (WENO-PR) is developed by distributing larger weights to the less smooth stencils. To meet the optimization requirements, the high-order smoothness

Since the discovery of vortex structures in turbulence, research in the area of vortex visualisation methods has been of interest. The prevalent Galilean invariant visualisation methods, like Q, Lambda2, and Omega methods, are based on the local pointwise analysis of invariants of the velocity gradient tensor. When plotting the vortex iso-surfaces, different methods will introduce different thresholds

Accurate modelling of fluid and particle flow is essential for reliable predictions of gas–solid reactions in the multiphase reactor. Due to the promising features such as short contact time, less back mixing of solids and low residence time, downer column has been used in industries than conventional riser configuration. However, the effective simulation and validation of the hydrodynamic characteristics

Multi-physics simulations are at the heart of today's engineering applications. The trend is towards more realistic and detailed simulations, which demand highly resolved spatial and temporal scales of various physical mechanisms to solve engineering problems in a reasonable amount of time. As a consequence, numerical codes need to run efficiently on high-performance computers. Therefore, the framework

ABSTRACT This paper presents an isogeometric (IGA) solver for time-dependent, incompressible magnetohydrodynamics (MHD). In this paper, a combination of inf-sup stable mixed discretisations is considered to discretise the hydrodynamic pair (i.e. velocity and pressure) and magnetic pair (i.e. magnetic field and magnetic pressure). The one step θ-method is used for the temporal discretisation. Manufactured

ABSTRACT The sensitivity of a large range of turbulence closures to the near-wall grid resolution is assessed through analytical and numerical techniques. We determine model boundedness, limiting behaviour, and solution dependence of channel flow simulations to the near-wall cell size . The sensitivity of present simulations to is strongly dependent on the choice of scale determining variable; and

Reduced Order Models (ROMs), also known as Reduced Basis Methods (RBMs), have received considerable attention in recent years for their ability to drastically reduce CFD cost, particularly when dealing with parametrised problems in a multi-query setting. This Special Issue gathers recent advances in ROM/RBM techniques for complex flow problems relevant to applications in mechanical and aerospace engineering

ABSTRACT In this contribution we consider cardiovascular hemodynamic modelling in patient-specific artery branches. To this aim, we first propose a procedure based on non-uniform rational basis splines (NURBS) to parametrise the artery volume which identifies the computational domain. Then, we adopt an isogeometric hierarchically reduced model which suitably combines separation of variables with a

ABSTRACT This paper demonstrates how, by means of Reduced Order Modelling (ROM), CFD cost can be drastically reduced and, in addition, a more complete investigation of a continuous design space obtained by adding experimental fluid dynamics (EFD) and flight fluid dynamics (FFD) data. The ‘Gappy’ Proper Orthogonal Decomposition (GPOD) method is used to enrich the EFD and FFD data to the size of their

ABSTRACT This paper investigates the recently introduced data-driven correction reduced order model (DDC-ROM) in the numerical simulation of the quasi-geostrophic equations. The DDC-ROM uses available data to model the correction term that is generally used to represent the missing information in low-dimensional ROMs. Physical constraints are added to the DDC-ROM to create the constrained data-driven

ABSTRACT Reversible axial flow fans, used in tunnel ventilation systems are designed in a way that as the direction of impeller rotation changes, suction and discharge directions change without any sensible change in flow rate and pressure. Also, their reversibility must be higher than 90%. Therefore, the blade profiles of these fans (S-shaped or elliptic) are designed symmetrically. In the current