Xiangda Cui, Ahmed Bakkar and Wagdi George Habashi
This paper aims to introduce a three-dimensional smoothed particle hydrodynamics (SPH) framework for simulating supercooled large droplets (SLD) dynamics at aeronautical speeds.
Abstract
Purpose
This paper aims to introduce a three-dimensional smoothed particle hydrodynamics (SPH) framework for simulating supercooled large droplets (SLD) dynamics at aeronautical speeds.
Design/methodology/approach
To include the effects of the surrounding air, a multiphase model capable of handling high density-ratio problems is adopted. A diffusive term is incorporated to smooth the density field and avoid numerical instabilities. Additionally, a particle shifting technique is used to eliminate anisotropic particle distributions.
Findings
The framework is validated against low-speed droplet impingement experimental results and then applied to the droplet impingement at high speeds typical of SLD conditions. Preliminary parametric studies are conducted to investigate the post-impact splashing. It is observed that a thicker water film can decrease the crown diameter and a smaller impact angle can suppress upward and forward splashing.
Originality/value
A three-dimensional multiphase SPH framework for SLD dynamics at a wide range of impact speed is developed and validated. The effects of particle resolution, water film thickness and impact angle on the post-impact crown evolution are investigated.
Details
Keywords
Adam Targui and Wagdi George Habashi
Responsible for lift generation, the helicopter rotor is an essential component to protect against ice accretion. As rotorcraft present a smaller wing cross-section and a lower…
Abstract
Purpose
Responsible for lift generation, the helicopter rotor is an essential component to protect against ice accretion. As rotorcraft present a smaller wing cross-section and a lower available onboard power compared to aircraft, electro-thermal heating pads are favored as they conform to the blades’ slender profile and limited volume. Their optimization is carried out here taking into account, for the first time, the highly three-dimensional (3D) nature of the flow and ice accretion, in contrast to the current state-of-the-art that is limited to two-dimensional (2D) airfoils.
Design/methodology/approach
Conjugate heat transfer simulation results are provided by the truly 3D finite element Navier–Stokes analysis package-ICE code, embedded in a proprietary rotorcraft simulation toolkit, with reduced-order modeling providing a time-efficient evaluation of the objective and constraint functions at every iteration. The proposed methodology optimizes heating pads extent and power usage and is versatile enough to address in a computationally efficient manner a wide variety of optimization formulations.
Findings
Low-error reduced-order modeling strategies are introduced to make the tackling of complex 3D geometries feasible in todays’ computers, with the developed framework applied to four problem formulations, demonstrating marked reductions to power consumption along with improved aerodynamics.
Originality/value
The present paper proposes a 3D framework for the optimization of electro-thermal rotorcraft ice protection systems, in hover and forward flight. The current state-of-the-art is limited to 2D airfoils.
Details
Keywords
Jory Seguin, Song Gao, Wagdi George Habashi, Dario Isola and Guido Baruzzi
This paper aims to describe the physical and numerical modeling of a new computational fluid dynamics solver for hypersonic flows in thermo-chemical non-equilibrium. The code uses…
Abstract
Purpose
This paper aims to describe the physical and numerical modeling of a new computational fluid dynamics solver for hypersonic flows in thermo-chemical non-equilibrium. The code uses a blend of numerical techniques to ensure accuracy and robustness and to provide scalability for advanced hypersonic physics and complex three-dimensional (3D) flows.
Design/methodology/approach
The solver is based on an edge-based stabilized finite element method (FEM). The chemical and thermal non-equilibrium systems are loosely-coupled to provide flexibility and ease of implementation. Chemical non-equilibrium is modeled using a laminar finite-rate chemical kinetics model while a two-temperature model is used to account for thermodynamic non-equilibrium. The systems are solved implicitly in time to relax numerical stiffness. Investigations are performed on various canonical hypersonic geometries in two-dimensional and 3D.
Findings
The comparisons with numerical and experimental results demonstrate the suitability of the code for hypersonic non-equilibrium flows. Although convergence is shown to suffer to some extent from the loosely-coupled implementation, trading a fully-coupled system for a number of smaller ones improves computational time. Furthermore, the specialized numerical discretization offers a great deal of flexibility in the implementation of numerical flux functions and boundary conditions.
Originality/value
The FEM is often disregarded in hypersonics. This paper demonstrates that this method can be used successfully for these types of flows. The present findings will be built upon in a later paper to demonstrate the powerful numerical ability of this type of solver, particularly with respect to robustness on highly stretched unstructured anisotropic grids.