Search results

1 – 10 of 30
Article
Publication date: 6 July 2015

He-yong Xu, Shi-long Xing and Zheng-yin Ye

The purpose of this paper is to investigate and improve a new method of unstructured rotational dynamic overset grids, which can be used to simulate the unsteady flows…

Abstract

Purpose

The purpose of this paper is to investigate and improve a new method of unstructured rotational dynamic overset grids, which can be used to simulate the unsteady flows around rotational parts of aircraft.

Design/methodology/approach

The computational domain is decomposed into two sub-domains, namely, the rotational sub-domain which contains the rotational boundaries, and the stationary sub-domain which contains the remainder flow field including the stationary boundaries. The artificial boundaries and restriction boundaries are used as the restriction condition to generate the entire computational grid, and then the overset grids are established according to the radius parameters of artificial boundaries set previously. The deformation of rotational boundary is treated by using the linear spring analogy method which is suitable for the dynamic unstructured grid. The unsteady Navier-Stokes/Euler equations are solved separately in the rotational sub-domain and stationary sub-domain, and data coupling is accomplished through the overlapping area. The least squares method is used to interpolate the flow variables for the artificial boundary points with a higher calculating precision. Implicit lower-upper symmetric-Gauss-Seidel (LU-SGS) time stepping scheme is implemented to accelerate the inner iteration during the unsteady simulation.

Findings

The airfoil steady flow, airfoil pitching unsteady flow, three-dimensional (3-D) rotor flow field, rotor-fuselage interaction unsteady flow field and the flutter exciting system unsteady flow field are numerically simulated, and the results have good agreements with the experimental data. It is shown that the present method is valid and efficient for the prediction of complicated unsteady problems which contain rotational dynamic boundaries.

Research limitations/implications

The results are entirely based on computational fluid dynamics (CFD), and the 3D simulations are based on the Euler equations in which the viscous effect is ignored. The current work shows further applicable potential to simulate unsteady flow around rotational parts of aircraft.

Practical implications

The current study can be used to simulate the two-dimensional airfoil pitching, 3-D rotor flow field, rotor-fuselage interaction and the flutter exciting system unsteady flow. The work will help the aircraft designer to get the unsteady flow character around rotational parts of aircraft.

Originality/value

A new type of rotational dynamic overset grids is presented and validated, and the current work has a significant contribution to the development of unstructured rotational dynamic overset grids.

Details

Aircraft Engineering and Aerospace Technology: An International Journal, vol. 87 no. 4
Type: Research Article
ISSN: 0002-2667

Keywords

Article
Publication date: 3 April 2018

Somayeh Harimi, Azam Marjani and Sadegh Moradi

This paper aims to study the fluid flow and forced convection heat transfer from an isothermal circular cylinder with control rods in the laminar unsteady flow regime.

Abstract

Purpose

This paper aims to study the fluid flow and forced convection heat transfer from an isothermal circular cylinder with control rods in the laminar unsteady flow regime.

Design/methodology/approach

The overset grid method was used for accurate simulation of the unsteady flows around different arrangements of the cylinders. Grid generation for overset grids was performed using a general orthogonal boundary fitted coordinate system. The method of solution was based on a finite volume discretization of the Navier-Stokes equations. Simulations were carried out for the Prandtl numbers of 0.7 and 7.0 with the Reynolds numbers ranging from 60 to 300.

Findings

The results indicate that the performance of multiple control rods depends strongly on the spacing ratio. Furthermore, in a manner similar to the flow patterns, four different thermal regimes were recognized based on the variations of mean Nusselt number versus G/D, as the thermal regimes follow the categories of flow regimes at different diameter ratios. However, for different Prandtl numbers, no single trend of heat transfer variation versus the spacing ratio exists for same regime.

Originality/value

Few studies have been conducted to investigate the heat transfer characteristics from control rods. The results of this study provide a comprehensive knowledge on the dynamical and thermal behavior of the flow around multiple cylinders.

Details

International Journal of Numerical Methods for Heat & Fluid Flow, vol. 28 no. 4
Type: Research Article
ISSN: 0961-5539

Keywords

Article
Publication date: 6 August 2019

Mohammad Haji Mohammadi and Joshua R. Brinkerhoff

Turbomachinery, including pumps, are mainly designed to extract/produce energy from/to the flow. A major challenge in the numerical simulation of turbomachinery is the…

Abstract

Purpose

Turbomachinery, including pumps, are mainly designed to extract/produce energy from/to the flow. A major challenge in the numerical simulation of turbomachinery is the inlet flow rate, which is routinely treated as a known boundary condition for simulation purposes but is properly a dependent output of the solution. As a consequence, the results from numerical simulations may be erroneous due to the incorrect specification of the discharge flow rate. Moreover, the transient behavior of the pumps in their initial states of startup and final states of shutoff phases has not been studied numerically. This paper aims to develop a coupled procedure for calculating the transient inlet flow rate as a part of the solution via application of the control volume method for linear momentum. Large eddy simulation of a four-blade axial hydraulic pump is carried out to calculate the forces at every time step. The sharp interface immersed boundary method is used to resolve the flow around the complex geometry of the propeller, stator and the pipe casing. The effect of the spurious pressure fluctuations, inherent in the sharp interface immersed boundary method, is damped by local time-averaging of the forces. The developed code is validated by comparing the steady-state volumetric flow rate with the experimental data provided by the pump manufacturer. The instantaneous and time-averaged flow fields are also studied to reveal the flow pattern and turbulence characteristics in the pump flow field.

Design/methodology/approach

The authors use control volume analysis for linear momentum to simulate the discharge rate as part of the solution in a large eddy simulation of an axial hydraulic pump. The linear momentum balance equation is used to update the inlet flow rate. The sharp interface immersed boundary method with dynamic Smagorinsky sub-grid stress model and a proper wall model is used.

Findings

The steady-state volumetric flow rate has been computed and validated by comparing to the flow rate specified by the manufacturer at the simulation conditions, which shows a promising result. The instantaneous and time averaged flow fields are also studied to reveal the flow pattern and turbulence characteristics in the pump flow field.

Originality/value

An approach is proposed for computing the volumetric flow rate as a coupled part of the flow solution, enabling the simulation of turbomachinery at all phases, including the startup/shutdown phase. To the best of the authors’ knowledge, this is the first large eddy simulation of a hydraulic pump to calculate the transient inlet flow rate as a part of the solution rather than specifying it as a fixed boundary condition. The method serves as a numerical framework for simulating problems incorporating complex shapes with moving/stationary parts at all regimes including the transient start-up and shut-down phases.

Details

International Journal of Numerical Methods for Heat & Fluid Flow, vol. 29 no. 7
Type: Research Article
ISSN: 0961-5539

Keywords

Article
Publication date: 23 June 2021

Mohammad Mehdi Razzaghi

This study aims to present a moving grid method based on the manipulation of connections.

Abstract

Purpose

This study aims to present a moving grid method based on the manipulation of connections.

Design/methodology/approach

In this study, the grid’s connections were manipulated to simulate a released store’s displacement. The selected model in this research is the EGLIN test case. In the introduced method, connections are modified in specific nodes of the grid. Governing flow equations were solved with the finite volume method. The major characteristic of this technique is using the averaging method for calculating the flux of cells.

Findings

This method maintains the grid’s quality even in large displacements of the released store. The three-dimensional simulation was carried out in transonic and supersonic regimes. Comparison of the results with experimental data were highly satisfactory.

Research limitations/implications

Using this moving grid method is recommended for simulating other models.

Practical implications

Prediction of store trajectory released from air vehicles is one of the most critical issues under study especially in the design of new stores.

Originality/value

The most prominent advantage of this method is maintaining the grid quality simultaneous with large displacements of the released store.

Details

Aircraft Engineering and Aerospace Technology, vol. 93 no. 5
Type: Research Article
ISSN: 1748-8842

Keywords

Article
Publication date: 4 December 2017

Gaohua Li, Xiang Fu and Fuxin Wang

This paper aims to improve the computational efficiency and to achieve high-order accuracy for the computation of helicopter rotor unsteady flows in forward flight during…

179

Abstract

Purpose

This paper aims to improve the computational efficiency and to achieve high-order accuracy for the computation of helicopter rotor unsteady flows in forward flight during the industrial preliminary design stage.

Design/methodology/approach

The integral arbitrary Lagrangian–Eulerian form of unsteady compressible Navier–Stokes equations with low Mach number preconditioned pseudo time terms based on non-inertial frame of reference undergoing rotating and translating was derived and discretized in the framework of multi-block structured finite volume grid using three types of spatial reconstruction schemes, i.e. the third-order accurate monotonic upwind scheme for conservation laws, the fifth-order accurate weighted essentially non-oscillatory and the fifth-order accurate weighted compact nonlinear schemes.

Findings

The results show that the present non-inertial computational method can obtain comparable results with other methods, such as the dynamic overset method, and make sure that the higher-order spatial schemes can significantly improve the tip vortex resolution.

Originality/value

The computational grid used by the present method remained static during the whole unsteady computation process, with only local deformations induced by blade cyclic pitch and other operating motions, which greatly reduced the complexity of grid motion and enhanced the efficiency and robustness.

Details

International Journal of Numerical Methods for Heat & Fluid Flow, vol. 27 no. 12
Type: Research Article
ISSN: 0961-5539

Keywords

Article
Publication date: 3 October 2019

Lisha He, Jianjing Zheng, Yao Zheng, Jianjun Chen, Xuan Zhou and Zhoufang Xiao

The purpose of this paper is to develop parallel algorithms for moving boundary simulations by local remeshing and compose them to a fully parallel simulation cycle for…

Abstract

Purpose

The purpose of this paper is to develop parallel algorithms for moving boundary simulations by local remeshing and compose them to a fully parallel simulation cycle for the solution of problems with engineering interests.

Design/methodology/approach

The moving boundary problems are solved by unsteady flow computations coupled with six-degrees-of-freedom equations of rigid body motion. Parallel algorithms are developed for both computational fluid dynamics (CFD) solution and grid deformation steps. Meanwhile, a novel approach is developed for the parallelization of the local remeshing step. It inputs a distributed mesh after deformation, then marks low-quality elements to be deleted on the respective processors. After that, a parallel domain decomposition approach is used to repartition the hole mesh and then to redistribute the resulting sub-meshes onto all available processors. Then remesh individual sub-holes in parallel. Finally, the element redistribution is rebalanced.

Findings

If the CFD solver is parallelized while the remaining steps are executed in sequential, the performance bottleneck of such a simulation cycle is observed when the simulation of large-scale problem is executed. The developed parallel simulation cycle, in which all of time-consuming steps have been efficiently parallelized, could overcome these bottlenecks, in terms of both memory consumption and computing efficiency.

Originality/value

A fully parallel approach for moving boundary simulations by local remeshing is developed to solve large-scale problems. In the algorithm level, a novel parallel local remeshing algorithm is present. It repartitions distributed hole elements evenly onto all available processors and ensures the generation of a well-shaped inter-hole boundary always. Therefore, the subsequent remeshing step can fix the inter-hole boundary involves no communications.

Details

Engineering Computations, vol. 36 no. 8
Type: Research Article
ISSN: 0264-4401

Keywords

Article
Publication date: 25 August 2021

Lorris Charrier, Mathieu Jubera, Grégoire Pont, Simon Marié, Pierre Brenner and Francesco Grasso

The design of a space launcher requires some considerations about the unsteady loads and heat transfer occurring at the base of the structure. In particular, these…

Abstract

Purpose

The design of a space launcher requires some considerations about the unsteady loads and heat transfer occurring at the base of the structure. In particular, these phenomena are predominant during the early stage of the flight. This paper aims to evaluate the ability of the unstructured, high order finite-volume CFD solver FLUSEPA, developed by Airbus Safran Launchers, to accurately describe these phenomena.

Design/methodology/approach

This paper first performs a steady simulation on a base flow around a four-clustered rocket configuration. Results are compared with NASA experiments and Loci-CHEM simulations. Then, unsteady simulations of supersonic H2/air reacting mixing layer based on the experiment of Miller, Bowman and Mungal are performed. Three meshes with different cells number are used to study the impact of spatial resolution. Instantaneous and time-averaged concentrations are compared with the combined OH/acetone planar laser-induced fluorescence imaging from the experiment.

Findings

FLUSEPA satisfactorily predicts the base heat flux at the base of a four-clustered rocket configuration. NASA Loci-CHEM reactive simulations indicate that afterburning plays an important role and should not be neglected. The unsteady reactive computation of a supersonic mixing layer shows that FLUSEPA is also able to accurately predict flow structures and interactions. However, the complexity of the experiment and the lack of details concerning the facility prevents from obtaining satisfactory converged results.

Originality/value

This study is the first step on the development of a cost-effective method aiming at predicting unsteady loads and heat transfer on space launchers using an unsteady and reactive model for the CDF calculations. It uses original techniques such as conservative CHIMERA-like overset grids, local re-centering of fluxes and local adaptive time-stepping to reduce computational cost while being robust and accurate.

Details

International Journal of Numerical Methods for Heat & Fluid Flow, vol. 31 no. 11
Type: Research Article
ISSN: 0961-5539

Keywords

Article
Publication date: 7 November 2016

Abdelkader Frendi and Michael R. Brown

The purpose of this paper is to carry out an extensive numerical study in order to understand the flow structures and the resulting noise generated by a supersonic…

Abstract

Purpose

The purpose of this paper is to carry out an extensive numerical study in order to understand the flow structures and the resulting noise generated by a supersonic impinging jet on a flat plate. One of the parameters varied in this study is the distance between the jet exit plane and the flat plate.

Design/methodology/approach

Because of the unsteady nature of the problem a time-dependent computation is carried out using the detached eddy simulation turbulence model. The OVERFLOW 2 CFD code was used with a highly resolved grid and small time steps.

Findings

The authors found that as the separation distance increases, the dominant frequencies in the noise spectrum decrease. In addition, the relative strength of the various frequencies to each other changes with changing distance, indicating the changing modes of the jet. The CFD results indicate a strong interaction between the acoustic waves emanating from the impingement plate and the jet plume. This feedback mechanism is responsible for destabilizing the jet shear layer leading to the jet changing modes. The computed near field spectra, convection velocities of the jet vortical structures and mean jet centerline velocity profile are in good agreement with experimental measurements. The results also show very high sound pressure levels all over the impingement plate but especially near the impingement point. These levels, if sustained, are detrimental to both human operators as well as the surrounding structures.

Research limitations/implications

Given the large-scale nature of the computations carried out, it is very costly to run the computations long enough to collect a good, statistically steady time sample to achieve a low frequency bandwidth resolution. Such a long time sample could actually improve the results in terms of frequency resolution and obtained an even better agreement with experiments. Off course there is always the issue of grid resolution as well, but given the good agreement with experiments that the authors obtained, the authors are confident in their results.

Practical implications

The practical implications of the results the authors obtained are significant in that, the authors now know that hybrid RANS-large eddy simulation methods can be used for this complex, unsteady engineering problems. In addition, the results also show the high noise level both on the impingement surface and in the surroundings of the jet. This could have a negative impact on the structural integrity of the flat surface.

Social implications

Noisy environments are never desirable anywhere especially in places where human operations take place. Therefore, given the high noise levels obtained in the simulations and confirmed by experiments, any human presence around the jet will be harmful to hearing and precautions need to be taken.

Originality/value

This is a physics-based study; i.e. understanding the physical phenomena involved in supersonic jet impingement. Of particular interest is the interaction of the jet shear layer with the acoustic waves emanating from the impingement area. This feedback loop is found to be responsible for intensifying the instability of the jet shear layer.

Details

International Journal of Numerical Methods for Heat & Fluid Flow, vol. 26 no. 8
Type: Research Article
ISSN: 0961-5539

Keywords

Article
Publication date: 20 July 2010

C.‐B. Liu, P. Nithiarasu and P.G. Tucker

The purpose of this paper is to numerically solve Eikonal and Hamilton‐Jacobi equations using the finite element method; to use both explicit Taylor Galerkin (TG) and…

Abstract

Purpose

The purpose of this paper is to numerically solve Eikonal and Hamilton‐Jacobi equations using the finite element method; to use both explicit Taylor Galerkin (TG) and implicit methods to obtain shortest wall distances; to demonstrate the implemented methods on some realistic problems; and to use iterative generalized minimal residual method (GMRES) method in the solution of the equations.

Design/methodology/approach

The finite element method along with both the explicit and implicit time discretisations is employed. Two different forms of governing equations are also employed in the solution. The Eikonal equation in its original form is used in the explicit Taylor Galerkin discretisation to save computational time. For implicit method, however, the convection‐diffusion form in its conservation form is used to maintain spatial stability.

Findings

The finite element solution obtained is both accurate and smooth. As expected the implicit method is much faster than the explicit method. Though the proposed finite element solution procedures in serial is slower than the standard search procedure, they are suitable to be used in a parallel environment.

Originality/value

The finite element procedure for Eikonal and Hamilton‐Jacobi equations are attempted for the first time. Though the finite volume and finite difference‐based computational fluid dynamics (CFD) solvers have started employing differential equations for wall distance calculations, it is not common for finite element solvers to use such wall distance calculations. The results presented here clearly show that the proposed methods are suitable for unstructured meshes and finite element solvers.

Details

Engineering Computations, vol. 27 no. 5
Type: Research Article
ISSN: 0264-4401

Keywords

Article
Publication date: 5 October 2015

Kwanchai Chinwicharnam, Edgard David Gomez Ariza, Jean-Marc Moschetta and Chinnapat Thipyopas

The purpose of this paper is to investigate the influence of a propeller slipstream on the aerodynamic characteristics of a fixed-wing micro air vehicle (MAV) by…

Abstract

Purpose

The purpose of this paper is to investigate the influence of a propeller slipstream on the aerodynamic characteristics of a fixed-wing micro air vehicle (MAV) by simplifying a propeller to an actuator disk and an actuator volume.

Design/methodology/approach

A computational fluid dynamic (CFD) approach.

Findings

The simulation flows are found and show that the propeller slipstream changes the flow field around the wing, which improves the aerodynamic performance of the wing. The aerodynamic performance is improved first, when the separation of the boundary flow at the upper surface wing is delayed. Second, the flow region of the boundary layer is boosted close to the wing surface again at a high incidence angle. And finally, the velocity inlet of the wing is increased by the propeller-induced flow.

Research limitations/implications

The incidence angle is in the range of 0-80°with an increment of 20°. The free stream velocity and RPM used are 6 m/s and 5,000 rpm, respectively.

Originality/value

A propeller is simplified to an actuator disk and an actuator volume.

Details

Aircraft Engineering and Aerospace Technology: An International Journal, vol. 87 no. 6
Type: Research Article
ISSN: 0002-2667

Keywords

1 – 10 of 30