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 (L1, L2, r, α) and canard angle from horizontal axis (θ). The objective functions include magnitude of drag and lift coefficients

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

(2020). Special Issue on Reduced Order Models in CFD. International Journal of Computational Fluid Dynamics: Vol. 34, Advances in reduced order methods in CFD (GRSP). Guest Editors: Simona Perotto and Gianluigi Rozza, pp. 91-92.

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 predictions

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 different

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

A coupling technique for interface simulations applied to the fields of ship hydrodynamics and multiphase flow computations is presented. The technique takes advantages of mass conservation property of the volume-of-fluid method and sharpened interface computation of the Level-Set approach. An improved scheme for the gradient computation of the Level-Set function during re-initialisation is also introduced

The multi-mesh, multi-solver paradigm makes use of multiple Computational Fluid Dynamics (CFD) solvers in a single overset framework. A framework-level implementation of the Generalised Minimal Residual algorithm applied to the full implicit overset system is presented. The method requires only minimal changes to existing, Python-wrapped CFD solvers and demonstrates improved convergence compared to

When the Reynolds number exceeds approximately 10, drag from porous media becomes nonlinear and cannot be handled by Darcy's theory. The Darcy–Forchheimer law, which considers drag through porous media as a quadratic function, covers this region up to the Reynolds number of the order 102. In this research, we study the optimal shape of a porous unit cell based on this law and topology optimisation

A combination of the lattice-Boltzmann and immersed-boundary methods is used to simulate the behaviour of circular particles in fluids, including their behaviour in stationary and flowing fluids. This method simulates the fluid domain with a regular Eulerian grid, and the non-grid uniform boundaries with a Lagrangian grid. The no-slip boundary condition of the fluid and particle boundaries is enforced

Additive Manufacturing (AM) is a rapidly developing new technology which allows the manufacture of arbitrarily complex geometries, and which is likely to transform heat exchanger design. To drive this transformation we need to develop computer modelling techniques to model fluid flow, heat exchange and phase change in arbitrarily complex domains, such as can be manufactured using AM. The present work

A linear stability analysis is performed to explain the formation of coherent structures in the narrow regions of a rectangular channel containing a cylindrical rod and a channel formed by two rectangular ducts connected through a narrow slot. This stability analysis introduces small perturbations of the velocity and the pressure in the governing equations. The linearised equations are transformed

In the paper, an improved third-order finite difference weighted essentially non-oscillatory scheme is presented for achieving the optimal order near critical points. A new global smoothness indicator is obtained with the way of Taylor expansion for the local smoothness indicator. In the framework of the conventional WENO-Z scheme, we present an improved third-order finite difference WENO scheme. The

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

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 and

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 y 1 + . The sensitivity of present simulations to y 1 + is strongly dependent on the choice of scale determining variable; ε

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 companion

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 correction

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 work, a

A typical weighted compact nonlinear scheme (WCNS) uses a convex combination of several low-order polynomials approximated over selected candidate stencils of the same width, achieving non-oscillatory interpolation near discontinuities and high-order accuracy for smooth solutions. In this paper, we present a new multi-resolution fifth-order WCNS by making use of the information of polynomials on three

Obtaining accurate solutions of unsteady flows during the design process of an aircraft can be a highly demanding task. Reduced basis methods (RBM) are commonly used to reduce the number of degrees of freedom while preserving high accuracy. RBM based on Proper Othogonal Decomposition (POD) have been extensively used but limitations exist with unsteady problems, where the temporal nonlinear dynamics

Reduced order modelling (ROM) techniques allow to reduce the cost of shape optimisation problems. In the present work, the compressible turbulent flow around a gas turbine profile is studied by a Discontinuous Galerkin (DG) scheme. The simulations are accelerated by the combination of two existing ROM approaches: Domain Decomposition and Local Proper Orthogonal Decomposition in a DG (LPOD-DG) framework

In this contribution, we derive an a posteriori error estimator for the second-order wave equation motivated by energy-based a priori estimates by Bernardi and Süli [“Time and Space Adaptivity for the Second-order Wave Equation.” Mathematical Models and Methods in Applied Sciences 15 (2): 199–225]. This estimate (which is valid for general discretisations) is then used to derive a POD-Greedy reduced

Arteries exhibit a remarkable ability to adapt in response to sustained alterations in hemodynamic loading, to heal in response to injuries, and to compensate in response to diverse disease conditions. Nevertheless, such compensatory adaptations are limited and many vascular disorders, if untreated, lead to significant morbidity or mortality. Parallel advances in vascular biology, medical imaging,

We consider the Navier–Stokes equations in a channel with a narrowing and walls of varying curvature. By applying the empirical interpolation method to generate an affine parameter dependency, the offline-online procedure can be used to compute reduced order solutions for parameter variations. The reduced-order space is computed from the steady-state snapshot solutions by a standard POD procedure.