Search results

The purpose of this study is to propose a precise and standardized strategy for numerically simulating vehicle aerodynamics.

Error sources in computational fluid dynamics were analyzed. Additionally, controllable experiential and discretization errors, which significantly influence the calculated results, are expounded upon. Considering the airflow mechanism around a vehicle, the computational efficiency and accuracy of each solution strategy were compared and analyzed through numerous computational cases. Finally, the most suitable numerical strategy, including the turbulence model, simplified vehicle model, calculation domain, boundary conditions, grids and discretization scheme, was identified. Two simplified vehicle models were introduced, and relevant wind tunnel tests were performed to validate the selected strategy.

Errors in vehicle computational aerodynamics mainly stem from the unreasonable simplification of the vehicle model, calculation domain, definite solution conditions, grid strategy and discretization schemes. Using the proposed standardized numerical strategy, the simulated steady and transient aerodynamic characteristics agreed well with the experimental results.

Building upon the modified Low-Reynolds Number k-e model and Scale Adaptive Simulation model, to the best of the authors’ knowledge, a precise and standardized numerical simulation strategy for vehicle aerodynamics is proposed for the first time, which can be integrated into vehicle research and design.

The paper aims to expand the scope of application of the lattice Boltzmann method (LBM), especially in the field of aircraft engineering. The traditional LBM is usually applied to incompressible flows at a low Reynolds number, which is not sufficient to satisfy the needs of aircraft engineering. Devoted to tackling the defect, the paper proposes a developed LBM combining the subgrid model and the multiple relaxation time (MRT) approach. A multilayer adaptive Cartesian grid method to improve the computing efficiency of the traditional LBM is also employed.

The subgrid model and the multilayer adaptive Cartesian grid are introduced into MRT-LBM for simulations of incompressible flows at a high Reynolds number. Validated by several typical flow simulations, the numerical methods in this paper can efficiently study the flows under high Reynolds numbers.

Some numerical simulations for the lid-driven flow of cavity, flow around iced GLC305, LB606b and ONERA-M6 are completed. The paper presents the investigation results, indicating that the methods are accurate and effective for the separated flow after icing.

LBM is developed with the addition of the subgrid model and the MRT method. A numerical strategy is proposed using a multilayer adaptive Cartesian grid method and its treatment of boundary conditions. The paper refers to innovative algorithm developments and applications to the aircraft engineering, especially for iced wing simulations with flow separations.

This paper aims to present an unconditionally energy-stable scheme for approximating the convective heat transfer equation.

The scheme stems from the generalized positive auxiliary variable (gPAV) idea and exploits a special treatment for the convection term. The original convection term is replaced by its linear approximation plus a correction term, which is under the control of an auxiliary variable. The scheme entails the computation of two temperature fields within each time step, and the linear algebraic system resulting from the discretization involves a coefficient matrix that is updated periodically. This auxiliary variable is given by a well-defined explicit formula that guarantees the positivity of its computed value.

Compared with the semi-implicit scheme and the gPAV-based scheme without the treatment on the convection term, the current scheme can provide an expanded accuracy range and achieve more accurate simulations at large (or fairly large) time step sizes. Extensive numerical experiments have been presented to demonstrate the accuracy and stability performance of the scheme developed herein.

This study shows the unconditional discrete energy stability property of the current scheme, irrespective of the time step sizes.

This study aims to analyse the non-linear losses of a porous media (stack) composed by parallel plates and inserted in a resonator tube in oscillatory flows by proposing numerical correlations between pressure gradient and velocity.

The numerical correlations origin from computational fluid dynamics simulations, conducted at the microscopic scale, in which three fluid channels representing the porous media are taken into account. More specifically, for a specific frequency and stack porosity, the oscillating pressure input is varied, and the velocity and the pressure-drop are post-processed in the frequency domain (Fast Fourier Transform analysis).

It emerges that the viscous component of pressure drop follows a quadratic trend with respect to velocity inside the stack, while the inertial component is linear also at high-velocity regimes. Furthermore, the non-linear coefficient b of the correlation ax + bx² (related to the Forchheimer coefficient) is discovered to be dependent on frequency. The largest value of the b is found at low frequencies as the fluid particle displacement is comparable to the stack length. Furthermore, the lower the porosity the higher the Forchheimer term because the velocity gradients at the stack geometrical discontinuities are more pronounced.

The main novelty of this work is that, for the first time, non-linear losses of a parallel plate stack are investigated from a macroscopic point of view and summarised into a non-linear correlation, similar to the steady-state and well-known Darcy–Forchheimer law. The main difference is that it considers the frequency dependence of both Darcy and Forchheimer terms. The results can be used to enhance the analysis and design of thermoacoustic devices, which use the kind of stacks studied in the present work.

This new paper aims to extend the authors’ previous contributions about open-loop aircraft flutter test to closed-loop aircraft flutter test by virtue of the proposed direct data–driven strategy. After feeding back the output signal to the input and introducing one feedback controller in the adding feedback loop, two parts, i.e. unknown aircraft flutter model and unknown feedback controller, exist in this closed-loop aircraft flutter system, simultaneously, whose input and output are all corrupted with external noise. Because of the relations between aircraft flutter model parameters and the unknown aircraft model, direct data–driven identification is proposed to identify that aircraft flutter model, then some identification algorithms and their statistical analysis are given through the authors’ own derivations. As the feedback controller can suppress the aircraft flutter or guarantee the flutter response converge to one desired constant value, the direct data–driven control is applied to design that feedback controller only through the observed data sequence directly. Numerical simulation results have demonstrated the efficiency of the proposed direct data–driven strategy. Generally, during our new information age, direct data–driven strategy is widely applied around our living life.

First, consider one more complex closed loop stochastic aircraft flutter model, whose input–output are all corrupted with external noise. Second, for the identification problem of closed-loop aircraft flutter model parameters, new identification algorithm and some considerations are given to the corresponding direct data–driven identification. Third, to design that feedback controller, existing in that closed-loop aircraft flutter model, direct data–driven control is proposed to design the feedback controller, which suppresses the flutter response actively.

A novel direct data–driven strategy is proposed to achieve the dual missions, i.e. identification and control for closed-loop aircraft flutter test. First, direct data–driven identification is applied to identify that unknown aircraft flutter model being related with aircraft flutter model parameters identification. Second, direct data–driven control is proposed to design that feedback controller.

To the best of the authors’ knowledge, this new paper extends the authors’ previous contributions about open-loop aircraft flutter test to closed-loop aircraft flutter test by virtue of the proposed direct data–driven strategy. Consider the identification problem of aircraft flutter model parameters within the presented closed loop environment, direct data–driven identification algorithm is proposed to achieve the identification goal. Direct data–driven control is proposed to design the feedback controller, i.e. only using the observed data to design the feedback controller.

Abnormal vibrations often occur in the liquid oxygen kerosene transmission pipelines of rocket engines, which seriously threaten their safety. Improper handling can result in failed rocket launches and significant economic losses. Therefore, this paper aims to examine vibrations in transmission pipelines.

In this study, a three-dimensional high-pressure pipeline model composed of corrugated pipes, multi-section bent pipes, and other auxiliary structures was established. The fluid–solid coupling method was used to analyse vibration characteristics of the pipeline under various external excitations. The simulation results were visualised using MATLAB, and their validity was verified via a thermal test.

In this study, the vibration mechanism of a complex high-pressure pipeline was examined via a visualisation method. The results showed that the low-frequency vibration of the pipe was caused by fluid self-excited pressure pulsation, whereas the vibration of the engine system caused a high-frequency vibration of the pipeline. The excitation of external pressure pulses did not significantly affect the vibrations of the pipelines. The visualisation results indicated that the severe vibration position of the pipeline thermal test is mainly concentrated between the inlet and outlet and between the two bellows.

The results of this study aid in understanding the causes of abnormal vibrations in rocket engine pipelines.

The causes of different vibration frequencies in the complex pipelines of rocket engines and the propagation characteristics of external vibration excitation were obtained.

The purpose of the paper is to analyse the performances of closures and compressibility corrections classically used in turbulence models when applied to highly-compressible turbulent boundary layers (TBLs) over flat plates.

A direct numerical simulation (DNS) database of TBLs, covering a wide range of thermodynamic conditions, is presented and exploited to perform a priori analyses of classical and recent closures for turbulent models. The results are systematically compared to the “exact” terms computed from DNS.

The few compressibility corrections available in the literature are not found to capture DNS data much better than the uncorrected original models, especially at the highest Mach numbers. Turbulent mass and heat fluxes are shown not to follow the classical gradient diffusion model, which was shown instead to provide acceptable results for modelling the vibrational turbulent heat flux.

The main originality of the present paper resides in the DNS database on which the a priori tests are conducted. The database contains some high-enthalpy simulations at large Mach numbers, allowing to test the performances of the turbulence models in the presence of both chemical dissociation and vibrational relaxation processes.

This paper aims to report the study of a fluid buoy system that includes wave effects, with particular emphasis on validating the numerical results with experimental data.

A fluid–solid coupled algorithm is proposed to describe the motion of a rigid buoy under the effects of waves. The Navier–Stokes equations are solved with the open-source finite volume package Code Saturne, in which a free-surface capture technique and equations of motion for the solid are implemented. An ad hoc experiment on a laboratory scale is built. A buoy is placed into a tank partially filled with water; the tank is mounted into a shake table and subjected to controlled motion that promotes waves. The experiment allows for recording the evolution of the free surface at the control points using the ultrasonic sensors and the movement of the buoy by tracking the markers by postprocessing the recorded videos. The numerical results are validated by comparison with the experimental data.

The implemented free-surface technique, developed within the framework of the finite-volume method, is validated. The best-obtained agreement is for small amplitudes compatible with the waves evolving under deep-water conditions. Second, the algorithm proposed to describe rigid-body motion, including wave analysis, is validated. The numerical body motion and wave pattern satisfactorily matched the experimental data. The complete 3D proposed model can realistically describe buoy motions under the effects of stationary waves.

The novel aspects of this study encompass the implementation of a fluid–structure interaction strategy to describe rigid-body motion, including wave effects in a finite-volume context, and the reported free-surface and buoy position measurements from experiments. To the best of the authors’ knowledge, the numerical strategy, the validation of the computed results and the experimental data are all original contributions of this work.

The purpose of this paper is to present a fast and bare bones implementation of a numerical method for quickly simulating turbulent thermal flows on GPUs. The work also validates earlier research showing that the lattice Boltzmann method (LBM) method is suitable for complex thermal flows.

A dual lattice hydrodynamic (D3Q27) thermal (D3Q7) multiple-relaxation time LBM model capable of thermal DNS calculations is implemented in CUDA.

The model has the same computational performance compared to earlier publications of similar LBM solvers. The solver is validated against three benchmark cases for turbulent thermal flow with available data and is shown to be in excellent agreement.

The combination of a D3Q27 and D3Q7 stencil for a multiple relaxation time -LBM has, to the authors’ knowledge, not been used for simulations of thermal flows. The code is made available in a public repository under a free license.

This paper aims to present a novel approach for computing particle temperatures in simulations coupling computational fluid dynamics (CFD) and discrete element method (DEM) to predict flow and heat transfer in fluidized beds of thermally thick spherical particles.

An improved lumped formulation based on Hermite-type approximations for integrals to relate surface temperature to average temperature and surface heat flux is used to overcome the limitations of classical lumped models. The model is validated through comparisons with analytical solutions for a convectively cooled sphere and experimental data for a fixed particle bed. The coupled CFD-DEM model is then applied to simulate a Geldart D bubbling fluidized bed, comparing the results to those obtained using the classical lumped model.

The validation cases demonstrate that ignoring internal thermal resistance can significantly impact the temperature in cases where the Biot number is greater than 0.1. The results for the fixed bed case clearly demonstrate that the proposed method yields significantly improved outcomes compared to the classical model. The fluidized bed results show that surface temperature can deviate considerably from the average temperature, underscoring the importance of accurately accounting for surface temperature in convective heat transfer predictions and surface processes.

The proposed approach offers a physically more consistent simulation without imposing a significant increase in computational cost. The improved lumped formulation can be easily and inexpensively integrated into a typical DEM solver workflow to predict heat transfer for spherical particles, with important implications for various industrial applications.